BACKGROUND
yersinia pestis and bacillus anthracis are two pathogens of significant concern to public health from a biodefense perspective  <cit> .
y. pestis, the causative agent of plague, is a gram-negative, highly communicable coccobacillus that has been responsible for three historic pandemics with high mortality rates  <cit> .
the microorganism possesses a type iii secretion mechanism common to several human, animal and plant pathogens, whereby a series of pathogen-specific structural proteins form a syringe-like structure capable of injecting virulence factors into the mammalian host cell.
these virulence factors then facilitate pathogen use of host nutrients and thwart the host immune response, ultimately causing cell and host death  <cit> .
naturally occurring plague can be transmitted from infected fleas and rodents to humans, and although the pathogen can be phagocytosed, it can also resist destruction by manipulating the host defense mechanism, potentially through antigenic mimicry  <cit> .
y. pestis then multiplies rapidly leading to necrosis of lymph nodes, a condition known as bubonic plague, which can result in death if untreated  <cit> .
in some cases the infection can spread through the blood stream resulting in systemic plague  or to the lungs resulting in the highly contagious and deadly form of the disease known as pneumonic plague.
there are currently no rapid, widely available diagnostic tests for plague, and the most common treatment is streptomycin  <cit> , an antibiotic with adverse effects.
two other species from the genus yersinia are also human pathogens: y. pseudotuberculosis and y. enterocolitica <cit> .
despite their high degree of sequence similarity to y. pestis, these two near neighbors of y. pestis manifest in very different symptoms, ranging from abdominal pain to septicemia in humans, usually caused by infection through contaminated food.
infections caused by y. pseudotuberculosis or y. enterocolitica can be effectively treated with antibiotics and in most cases are self-limiting.
notably, y. pestis is reported to have evolved from y. pseudotuberculosis within the past  <dig>  years  <cit> .
b. anthracis is a gram-positive, rod-shaped spore-forming bacterial pathogen and the causative agent of anthrax  <cit> .
human, livestock, and wildlife mortalities attributable to anthrax occur in numerous regions of the world, although the majority of cases are found in less industrialized nations  <cit> .
three forms of the disease have been described: cutaneous, intestinal and inhalational.
while cutaneous and intestinal forms may be less severe, inhalational anthrax is often fatal without prompt antibiotic treatment  <cit> .
the primary mechanisms of virulence employed by b. anthracis are associated with two virulence plasmids designated pxo <dig> and pxo <dig>  <cit> .
the net effect of these plasmids is virtually unhindered proliferation of b. anthracis within the host, hemorrhaging, cardio-pulmonary collapse, and death.
the regulation of production of host cytokines by both yersinia and b. anthracis has been described previously.
pickering a. k. et.
al.
measured cytokine levels in human dendritic cell supernatant and in mouse peritoneal macrophages exposed to b. anthracis spores  <cit> .
they observed significant increase in tnf-α, il- <dig>  il-1β, il- <dig>  and il- <dig> in human dendritic cell supernatants by  <dig> hours post-exposure.
high levels of il- <dig>  and tnf-α were observed in the supernatant from b. anthracis infected mouse peritoneal macrophages  <cit> .
in a mouse model,  <dig> cytokines, namely il-12p <dig>  tnf, ifn-γ, mcp- <dig>  il- <dig>  and il- <dig>  were increased significantly in mouse lung at  <dig> hours of y. pestis infection  <cit> .
in previous work comparing exposures to different bacterial pathogens, distinct patterns of cytokine expression levels were found that could discriminate the particular host response  <cit> , including while using pathogen-specific lps in whole blood  <cit> .
the hypothesis for the present study is that exposure to diverse bacterial pathogen strains would result in distinct cytokine profiles in the host, with strains from the same species exhibiting more similar profiles than strains from phylogenetically distant species.
a multiplex cytokine protein chip was used, and a multivariate approach was taken that combined expression data on multiple cytokines.
multivariate clustering techniques were used to establish cytokine expression profiles after ex vivo exposure of whole blood to seven pathogens.
methods
bacterial strains and culture conditions
the bacterial strains used in this study include: b. anthracis ames , b. anthracis sterne , y. pestis kim <dig> d <dig> .
y. pestis india/p , and y. pestis nyc , y. pseudotuberculosis serotype  <dig> pb <dig>  and y. enterocolitica wa serovar 0: <dig>  bacteria were grown on tryptose blood agar slants at 26°c for 1- <dig> days and subsequently collected using  <dig> ml of  <dig> m potassium-phosphate, ph  <dig> ;.bacterial densities were measured at od <dig> .
whole blood ex vivo exposure model
human blood was collected from a healthy donor by venipuncture using cpt vacutainer tubes  containing citrate.
informed consent was obtained and our blood collection protocol was approved by the llnl irb committee.
separate cpt tubes were used for the unexposed control and  <dig> different bacterial exposures .
bacteria were added to blood within  <dig> minutes of collection at a multiplicity of infection ratio of 5: <dig>  this ratio was determined against white blood cells in whole blood:∼ <dig> x  <dig> cells/ml.
each whole blood sample was incubated with bacteria for  <dig> hours at 37°c in 5% co <dig> following incubation, plasma was collected by centrifugation at  <dig> x g for  <dig> min at 4°c.
the control plasma was obtained in the same way and treated with  <dig>  m potassium-phosphate as a mock exposure.
these plasma samples were used for cytokine measurements.
cytokine immunoassays with protein arrays
the measurements of cytokines were performed using zyomyx protein profiling biochips .
these protein arrays allow the simultaneous quantification of  <dig> biologically relevant cytokines, as determined by zyomyx, inc: il-1α, il-1β, il- <dig>  il- <dig>  il- <dig>  il- <dig>  il- <dig>  il- <dig>  il- <dig>  il- <dig>  il- <dig>  il- <dig>  il- <dig>  il- <dig>  il- <dig>  tnfα, tnfβ, eotaxin, mcp- <dig>  mcp- <dig>  trail, cd <dig>  mig, sicam- <dig>  ip- <dig>  cd <dig>  tgf-β, gm-csf, gcsf, ifn-γ.
each cytokine assay was optimized for the zyomyx protein profiling biochip based on many factors including the availability of antibodies and the sensitivity and specificity of antibody-cytokine interactions.
each protein array chip is designed with  <dig> independent microfluidic channels that allow up to  <dig> samples to be loaded into isolated regions of an array.
antibodies specific for  <dig> analytes were arrayed in each channel, and each antibody was arrayed in redundancy on  <dig> pillars within the channel.
accordingly, a cytokine measurement represents the average of  <dig> measurements.
all immunoassay steps, including sample loading, washing, and detection, were performed with a fully automated biochip processing station .
eight protein array chips were used in these experiments.
two chips were used for generating calibration curves with a calibration standard kit containing  <dig> analytes .
sample  was injected into each channel of the protein array chips.
standard solutions were applied to two channels of each chip for chip-to-chip normalization.
triplicates of control and pathogen-exposed plasmas were applied randomly to four channels of  <dig> protein array chips.
protein arrays were scanned at  <dig> nm with zyomyx scanner  <dig> after immunoassays.
zyomyx data reduction software was used for normalization, calculation of calibration curves.
dixon’s test was used to remove outliers, and the median feature intensity was background subtracted.
concentrations of cytokines in plasma samples were determined by a four parameter logistic model.
cluster analysis of cytokine data
multiple hierarchical clustering methods were used to group the pathogen exposures based on the multivariate cytokine expression profiles induced in a host infection model system.
first, hierarchical agglomerative clustering  <cit>  was applied to group the control and the seven pathogen-exposed samples based on their cytokine concentration profiles.
each of the eight samples was characterized by the multivariate vector of its average log <dig> concentration over the cytokines.
the samples were then clustered based on the following distance measures between the samples and between the clusters.
distance between two samples was defined using two distance metrics:
· euclidean distance
· correlation distance:
distance between two clusters was defined using three methods:
· complete linkage : the largest distance between members of the clusters
· single linkage : the smallest distance between members of the clusters
· average linkage : the average distance between members of the clusters
given a pair of distance metrics between samples and clusters, the algorithm was initialized with the eight samples forming eight different clusters and then processed iteratively by joining the two most similar clusters.
the tree was built starting from the individual samples, using an agglomerative  approach.
the resulting hierarchy of clusters was displayed as a dendrogram.
these traditional clustering methods provide a quick, exploratory overview of the data.
however, these methods do not estimate the optimal number of clusters in the data; rather, the clustering is performed exhaustively from the lowest possible level of the hierarchy where each sample forms its own cluster, to the highest level where all samples are grouped into one cluster.
in addition to the traditional hierarchical agglomerative clustering method, the hierarchical ordered partitioning and collapsing hybrid  algorithm was also applied to the cytokine measurements  <cit> .
in contrast with the previous approaches where the tree was built starting from the individual samples as the leaf nodes, hopach used a hybrid divisive-agglomerative approach: it started from the root cluster containing all the samples , then divided the root down to leaf nodes, with an extra collapsing  step after each iteration that combined similar clusters.
based on the correlation distance between samples, hopach determined the split that minimized a measure of cluster homogeneity called the median split silhouette.
while computationally more expensive than the previous methods, hopach was expected to perform better because of its dynamic approach to update and potentially revise the clusters at every step of the iteration.
furthermore, hopach also estimated the optimal number of clusters from the data, and thus offered another advantage over the previous methods.
computations were performed in the r computing environment  and the hopach package  <cit> .
RESULTS
cytokine levels were examined using an ex vivo model, termed weem for whole blood x vivoexposure model.
individual samples of anti-coagulated human blood were incubated with b. anthracis ames, b. anthracis sterne, y. pestis kim <dig> d <dig>  y. pestis nyc, y. pestis india/p, y. enterocolitica, and y. pseudotuberculosis.
after  <dig> hours exposure, blood cells were removed by low-speed centrifugation and concentrations of  <dig> cytokines in the plasma were measured with protein arrays.
concentrations of fourteen cytokines, gcsf, ifnγ, gm-csf, il- <dig>  il- <dig>  il- <dig>  il- <dig>  il- <dig>  il- <dig>  il- <dig>  il- <dig>  mcp- <dig>  tgfβ, and tnfβ were below the limit of detection in this study.
the following  <dig> cytokines were detected: eotaxin, il- <dig>  il- <dig>  il- <dig>  il-1α, il-1β, il- <dig>  il- <dig>  ip- <dig>  mcp- <dig>  mig, tnfα, trail, scd <dig>  scd <dig>  and sicam- <dig> .
to determine if there were significant differences among the levels of cytokines in the control and pathogen exposed plasma samples, f-tests were performed.
for thirteen of these  <dig> cytokines, all three replicates were detected and these cytokines were subjected to f-tests.
statistical analysis indicated that  <dig> cytokines  had differentially elevated expression profiles following different bacterial exposures.
figure  <dig> shows the concentrations  of these cytokines in the control and bacteria exposed plasma samples.
the f-tests revealed that the other five cytokines containing complete datasets, trail, scd <dig>  scd <dig>  mig, and sicam- <dig>  had no significant difference between bacterial exposures and the mock-exposed control.
moreover, there was a great variation in absolute concentrations between cytokines.
for example, the concentrations of tnfα, scd <dig>  and sicam- <dig> were as high as  <dig> x  <dig> - <dig> pg/ml, whereas il- <dig> was much lower, about  <dig> pg/ml.
marked differences in induced cytokine patterns between b. anthracis and yersinia exposures were found.
also, the levels of induction of these cytokines differed among the different bacteria.
for example, yersinia species induced much higher cytokine response than b. anthracis for il-1α, il-1β, il- <dig>  and tnf-α .
the two strains of b. anthracis bacteria induced different levels of il-1β and tnf-α , including  <dig>  times higher concentration of il-1β and  <dig>  times higher concentration of tnf-α for the sterne strain than the ames strain of b. anthracis.
these differences were statistically significant .
to discriminate y. pestis exposure from near neighbors, il- <dig> levels can be used, showing cytokine concentrations following y. enterocolitica exposure and y. pseudotuberculosis exposure that are on average 5-fold higher and 2-fold higher, respectively, than after y. pestis exposure .
il- <dig> differential expression was specific to the yersinia spp.
because exposure to b. anthracis strains showed comparable il- <dig> levels to that in unexposed control.
the hopach algorithm estimated the number of clusters as five, and grouped the samples based on their host cytokine expression profiles as follows: 1) y. pestis , 2) y. pseudotuberculosis, 3) y. enterocolitica, 4) b. anthracis , and 5) control .
the closer the pathogen-exposed samples are within the tree on the left, the more similar they are.
height of the branches indicates the distance between the successive nodes in the clustering.
the method separated the b. anthracis and yersinia infected blood samples.
in addition, the cytokine profile of the mock-exposed control was more similar to the pattern produced by b. anthracis exposure than to the profile elicited by yersinia.
results of the hierarchical clustering when using the euclidean distance between samples depended on the distance metric between clusters.
the three methods for determining the distance between clusters  all established three major clusters: 1) y. pestis and near neighbors, 2) b. anthracis, and c) control.
differences between the results occurred when the yersinia cluster was further divided.
the average linkage method, consistent with figure  <dig>  formed a subgroup of the three y. pestis strains, then grouped them first with y. pseudotuberculosis followed by y. enterocolitica.
complete and single linkage methods, however, first grouped the attenuated virulent strain of y. pestis  with the more virulent strain , both clinical isolates from human plague cases, and then clustered them with y. pseudotuberculosis, followed by the attenuated y. pestis , and lastly with y. enterocolitica.
this is interesting from an evolutionary perspective because it has been proposed that y. pestis evolved from y. pseudotuberculosis within the last  <dig>  years, and thus these two pathogens are more closely related  <cit> .
when using hierarchical clustering with the correlation distance between the samples, the final clusters were independent of the distance metric between clusters, and agreed with the tree structure in figure  <dig>  the complete, single, and average linkage methods all resulted in the following major clusters: 1) yersinia, 2) b. anthracis, and 3).
control.
within the yersinia cluster, y. pestis  was closest to y. pestis , followed by y. pestis , y. pseudotuberculosis, and y. enterocolitica.
discussion
the hopach clustering method  produced five distinctly separated clusters: 1) y. pestis , 2) y. pseudotuberculosis, 3) y. enterocolitica, 4) b. anthracis , and 5) control.
this result is consistent with the findings using the correlation distance and the euclidean distance with average linkage.
in addition, hopach estimated the optimal number of clusters as five.
that is, the yersinia subcluster is best if it is divided into the three clusters specified by 1) through 5) above.
y. enterocolitica forms its own cluster, and so does y. pseudotuberculosis.
y. pestis , y. pestis , and y. pestis  are grouped into one cluster.
further subdivisions lead to an overall clustering with inferior quality.
in addition to clustering the cytokine expression profiles across bacterial treatments, figure  <dig> also groups the cytokines themselves and clusters the proteins based on their similarities across the pathogen exposures and reorders them accordingly.
interestingly, the three pro-inflammatory cytokines il-1β, tnfα, and il- <dig> clustered closely, and so did the three chemokines mcp- <dig>  ip- <dig>  and il- <dig>  although these  <dig> cytokines do not cluster as a single group, they do cluster at a branch further away from the leaf node, which includes il- <dig> and scd <dig>  to make a larger group of  <dig> proteins.
several of these proteins are involved in inflammatory conditions, such as il-1beta, tnfα, il- <dig>   <cit>  and have been shown to be upregulated in cell culture and animal model specifically exposed to biothreat agents  <cit> .
increased expression of il- <dig> and tnfα clustered together in a study involving mouse splenic cd11b + cells following sub-lethal y. enterocolitica infection  <cit> .
in addition, several of the cytokines in this cluster, namely tnf-alpha, il1-beta, il- <dig>  and mcp- <dig> are expressed higher in exposed whole blood as compared to control in this study and in whole blood exposure to lps from several other gram negative bacterial pathogens  <cit> .
in addition to expression differences, the absence of detected cytokine expression can also be helpful in discriminating pathogen exposure.
the multiplex detection of  <dig> cytokines in this study revealed the early phase cytokine expression profiles in human plasma following exposures to b. anthracis , y. pestis , y. pseudotuberculosis, and y. enterocolitica.
the expression levels of  <dig> cytokines, il-1α, il-1β, il- <dig>  il- <dig>  il- <dig>  ip- <dig>  mcp- <dig>  and tnfα were significantly different from that of unexposed control .
although the focus of our work was to show that cytokine expression profiling can discriminate between different pathogen exposures in a human whole blood ex vivo model, these results also represent an initial attempt to characterize the full cytokine response to each individual pathogen.
our preliminary study using a single exposure protocol at a single time post-exposure will need to be supplemented with more thorough investigation in order to determine the usefulness of using cytokine levels for diagnosing pathogen exposure.
however, the single time point chosen,  <dig> hours, is sufficient to detect proteomic changes and has been used in previous studies examining cytokine levels  <cit> .
this time point represents a start towards a more complete temporal study, as has been done with gene expression patterns for two of the pathogens studied here  <cit> .
in addition, studies that provide expression patterns for a single cytokine using multiple time points will also be needed to make the results of this paper clinically useful, such as has been done by, cooper and coworkers, who examined il-12p <dig> and il-12p <dig> levels following different growth conditions and exposure levels for a time course of y. pestis exposed dendritic cells  <cit> .
the results of the current work shows a similar expression pattern trend to this previous work, in which, y. pestis induces il-12p <dig> and at a substantially higher level than il-12p <dig>
our results showed that the expression levels of  <dig> chemokines, il- <dig>  mcp- <dig> and ip- <dig>  were induced by both yersinia and b. anthracis exposures.
no significant differences were found for these cytokines between yersinia and b. anthracis exposures.
il- <dig>  mcp- <dig> and ip- <dig> are chemokines that enable the migration of leukocytes from the blood to the site of inflammation.
il- <dig> is a key chemokine regulating neutrophil recruitment  <cit> .
the essential involvement of il- <dig> in acute inflammation was demonstrated by neutralizing il- <dig> with its antibody.
when highly specific antibody against il- <dig> was administered in acute inflammatory reactions induced by several stimuli including lipopolysaccharide, neutrophil infiltration was blocked  <cit> .
mcp- <dig> is known for its ability to act as potent chemoattractant and activator of monocytes/macrophages as well as nk cells but not neutrophils  <cit>  .
ip- <dig> has no chemotactic activity for neutrophils but attracts monocytes, nk, and t cells to the site of infection and regulates t cell maturation  <cit> .
it was reported previously that elevated il- <dig> and mcp- <dig> were secreted by human epithelial cells after y. enterocolitica infection, but not ip- <dig>  <cit> .
human dendritic cells, infected with b. anthracis spores, secreted high level of il- <dig> at  <dig>  hours  <cit> .
in our study, the fold increase of il- <dig> was much greater than mcp- <dig> and ip- <dig> .
for example, the induction of il- <dig> by ames strain of b. anthracis was  <dig> fold, while mcp- <dig> was  <dig> fold and ip- <dig> was  <dig>  fold .
this result may indicate that il- <dig> is a dominant chemokine in early response  and neutrophils are the major player in early inflammatory response.
here we compared cytokines induced by b. anthracis and yersinia exposures.
overall, yersinia exposure induced higher levels of il-1α, il-1β, il- <dig>  il- <dig> and tnfα than b. anthracis exposure, suggesting these cytokines could be used to develop an assay for discriminating yersinia spp.
from b. anthracis exposures.
the vaccine strain  of b. anthracis induced higher levels of il-1β and tnfα than the virulent strain  , suggesting these cytokines can contribute to a biomarker panel to discriminate if a particular isolate of b. anthracis is virulent.
there was also a difference in induction of il- <dig> between y. pestis and near neighbors , suggesting this cytokine is a candidate biomarker for discriminating the virulence of yersinia species.
these data regarding il- <dig> expression following yersinia spp.
exposure are in agreement with published literature that shows y. enterocolitica and y. pestis can elicit statistically different levels of il- <dig> expression  <cit> .
differences in il- <dig> induction may be due to differences in the lcrv protein among yersinia spp.
<cit> .
the different cytokine profiles induced by b. anthracis and yersinia here may be partially due to different surface antigens on the outermost part of these pathogens and the manner in which these bacteria were grown.
lipopolysaccharide , the main constituent of the outer membrane of gram-negative bacteria, and peptidoglycan , the major cell wall component of gram-positive bacteria, have been reported to elicit markedly different immune responses  <cit> .
however, virulence factors, such as b. anthracis lethal toxin and yersinia virulence antigen, lcrv, may also play important roles in differential cytokine induction.
this view is supported by numerous reports that b. anthracis toxin and virulence factors of yersinia bacteria  modulate host cytokine responses  <cit> .
while the various clustering methods resulted in slightly different final hierarchies, all were consistent in separating the unexposed control from the samples exposed to b. anthracis or to the y. pestis and near neighbors.
agreement on this level among the various clustering procedures lends more confidence to the overall results.
on a more detailed level, the methods grouped slightly differently the samples exposed to the y. pestis and near neighbors, which indicates that these samples cannot be unequivocally separated based on the current data and additional biomarkers or a larger sample set would be needed.
the most advanced hopach method estimated the optimal number of clusters in the data as five, corresponding to the unexposed control, and the four species: b. anthracis, y. pseudotuberculosis, y. enterocolitica, and y. pestis  .
information gained from the targeted protein array data for host response complements genomic  <cit> , and other proteomic studies  <cit>  of host-pathogen interactions.
the success of the weem and computational method to distinguish pathogen exposure, based on host response in this initial study, is encouraging and suggests a number of possibilities for future studies to refine the findings.
comparative analysis, such as the current work, can potentially reveal the critical pathogenic mechanism and host innate immune responses during infection as was previously shown for y. pestis and y. pseudotuberculosis <cit> .
opportunities include using statistical hypothesis tests based on analysis of variance to assess the significance of the observed differences among the host-pathogen cytokine concentration profiles, as well as performing follow-up studies to focus more on the y. pestis and near neighbor cluster.
in addition, the methods can be extended to investigate host responses to diverse pathogens in multiple host model systems to cross validate the significance of the biomarkers to distinguish pathogen exposures.
CONCLUSIONS
results from this study suggest that cytokine arrays coupled with statistical clustering methods can distinguish exposures to pathogens, including multiple strains of y. pestis, y. pseudotuberculosis, y. enterocolitica, and b. anthracis.
these methods differentiate both near neighbors and distant evolutionary microbes based on host response data.
the distinct cytokine profiles also provide insight into both the host response and virulence mechanisms of diverse pathogens.
in summary, characterization of host responses based on cytokine profiles has translational application, potentially providing the identification of infectious diseases and leading toward the ultimate goal of presymptomatic detection via sentinel surveillance of pathogen exposure and appropriate treatment.
competing interests
the authors declare that they have no competing interests.
authors’ contributions
bc participated in the design of the study, conducted experiments, and drafted and finalized the manuscript.
km conducted experiments.
if participated in the design of the study, performed the statistical analysis, and helped draft the manuscript.
pl helped draft the manuscript.
sm conceived of the study, obtained funding, and participated in its design and coordination and helped to draft the manuscript.
all authors read and approved the final manuscript.
acknowledgments
we thank david callender, jonathan e. forman, and renee tobias from zyomyx for their assistance with the biochip analyses.
we thank todd corzett for useful comments on the manuscript and brent ricks for technical assistance.
this work was performed under the auspices of the us department of energy by the university of california, lawrence livermore national laboratory under contract no.
w-7405-eng- <dig>  with support from the department of homeland security .
the authors would also like to thank psw rce animal resources and laboratory services core u54-ai <dig>  ucrl-jrnl- <dig>

 yersinia pestis and bacillus anthracis are two pathogens of significant concern to public health from a biodefense perspective. y. pestis, the causative agent of plague, is a gram-negative, highly communicable coccobacillus that has been responsible for three historic pandemics with high mortality rates. the microorganism possesses a type iii secretion mechanism common to several human, animal and plant pathogens, whereby a series of pathogen-specific structural proteins form a syringe-like structure capable of injecting virulence factors into the mammalian host cell. these virulence factors then facilitate pathogen use of host nutrients and thwart the host immune response, ultimately causing cell and host death. naturally occurring plague can be transmitted from infected fleas and rodents to humans, and although the pathogen can be phagocytosed, it can also resist destruction by manipulating the host defense mechanism, potentially through antigenic mimicry. y. pestis then multiplies rapidly leading to necrosis of lymph nodes, a condition known as bubonic plague, which can result in death if untreated. in some cases the infection can spread through the blood stream resulting in systemic plague or to the lungs resulting in the highly contagious and deadly form of the disease known as pneumonic plague. there are currently no rapid, widely available diagnostic tests for plague, and the most common treatment is streptomycin, an antibiotic with adverse effects. two other species from the genus yersinia are also human pathogens: y. pseudotuberculosis and y. enterocolitica. despite their high degree of sequence similarity to y. pestis, these two near neighbors of y. pestis manifest in very different symptoms, ranging from abdominal pain to septicemia in humans, usually caused by infection through contaminated food. infections caused by y. pseudotuberculosis or y. enterocolitica can be effectively treated with antibiotics and in most cases are self-limiting. notably, y. pestis is reported to have evolved from y. pseudotuberculosis within the past years. b. anthracis is a gram-positive, rod-shaped spore-forming bacterial pathogen and the causative agent of anthrax. human, livestock, and wildlife mortalities attributable to anthrax occur in numerous regions of the world, although the majority of cases are found in less industrialized nations. three forms of the disease have been described: cutaneous, intestinal and inhalational. while cutaneous and intestinal forms may be less severe, inhalational anthrax is often fatal without prompt antibiotic treatment. the primary mechanisms of virulence employed by b. anthracis are associated with two virulence plasmids designated pxo and pxo. the net effect of these plasmids is virtually unhindered proliferation of b. anthracis within the host, hemorrhaging, cardio-pulmonary collapse, and death. the regulation of production of host cytokines by both yersinia and b. anthracis has been described previously. pickering a. k. et. al. measured cytokine levels in human dendritic cell supernatant and in mouse peritoneal macrophages exposed to b. anthracis spores. they observed significant increase in tnf-α, il- il-1β, il- and il- in human dendritic cell supernatants by hours post-exposure. high levels of il- and tnf-α were observed in the supernatant from b. anthracis infected mouse peritoneal macrophages. in a mouse model, cytokines, namely il-12p tnf, ifn-γ, mcp- il- and il- were increased significantly in mouse lung at hours of y. pestis infection. in previous work comparing exposures to different bacterial pathogens, distinct patterns of cytokine expression levels were found that could discriminate the particular host response, including while using pathogen-specific lps in whole blood. the hypothesis for the present study is that exposure to diverse bacterial pathogen strains would result in distinct cytokine profiles in the host, with strains from the same species exhibiting more similar profiles than strains from phylogenetically distant species. a multiplex cytokine protein chip was used, and a multivariate approach was taken that combined expression data on multiple cytokines. multivariate clustering techniques were used to establish cytokine expression profiles after ex vivo exposure of whole blood to seven pathogens. methods bacterial strains and culture conditions the bacterial strains used in this study include: b. anthracis ames, b. anthracis sterne, y. pestis kim d. y. pestis india/p, and y. pestis nyc, y. pseudotuberculosis serotype pb and y. enterocolitica wa serovar 0: bacteria were grown on tryptose blood agar slants at 26°c for 1- days and subsequently collected using ml of m potassium-phosphate, ph;.bacterial densities were measured at od. whole blood ex vivo exposure model human blood was collected from a healthy donor by venipuncture using cpt vacutainer tubes containing citrate. informed consent was obtained and our blood collection protocol was approved by the llnl irb committee. separate cpt tubes were used for the unexposed control and different bacterial exposures. bacteria were added to blood within minutes of collection at a multiplicity of infection ratio of 5: this ratio was determined against white blood cells in whole blood:∼ x cells/ml. each whole blood sample was incubated with bacteria for hours at 37°c in 5% co following incubation, plasma was collected by centrifugation at x g for min at 4°c. the control plasma was obtained in the same way and treated with m potassium-phosphate as a mock exposure. these plasma samples were used for cytokine measurements. cytokine immunoassays with protein arrays the measurements of cytokines were performed using zyomyx protein profiling biochips. these protein arrays allow the simultaneous quantification of biologically relevant cytokines, as determined by zyomyx, inc: il-1α, il-1β, il- il- il- il- il- il- il- il- il- il- il- il- il- tnfα, tnfβ, eotaxin, mcp- mcp- trail, cd mig, sicam- ip- cd tgf-β, gm-csf, gcsf, ifn-γ. each cytokine assay was optimized for the zyomyx protein profiling biochip based on many factors including the availability of antibodies and the sensitivity and specificity of antibody-cytokine interactions. each protein array chip is designed with independent microfluidic channels that allow up to samples to be loaded into isolated regions of an array. antibodies specific for analytes were arrayed in each channel, and each antibody was arrayed in redundancy on pillars within the channel. accordingly, a cytokine measurement represents the average of measurements. all immunoassay steps, including sample loading, washing, and detection, were performed with a fully automated biochip processing station. eight protein array chips were used in these experiments. two chips were used for generating calibration curves with a calibration standard kit containing analytes. sample was injected into each channel of the protein array chips. standard solutions were applied to two channels of each chip for chip-to-chip normalization. triplicates of control and pathogen-exposed plasmas were applied randomly to four channels of protein array chips. protein arrays were scanned at nm with zyomyx scanner after immunoassays. zyomyx data reduction software was used for normalization, calculation of calibration curves. dixon’s test was used to remove outliers, and the median feature intensity was background subtracted. concentrations of cytokines in plasma samples were determined by a four parameter logistic model. cluster analysis of cytokine data multiple hierarchical clustering methods were used to group the pathogen exposures based on the multivariate cytokine expression profiles induced in a host infection model system. first, hierarchical agglomerative clustering was applied to group the control and the seven pathogen-exposed samples based on their cytokine concentration profiles. each of the eight samples was characterized by the multivariate vector of its average log concentration over the cytokines. the samples were then clustered based on the following distance measures between the samples and between the clusters. distance between two samples was defined using two distance metrics: · euclidean distance · correlation distance: distance between two clusters was defined using three methods: · complete linkage: the largest distance between members of the clusters · single linkage: the smallest distance between members of the clusters · average linkage: the average distance between members of the clusters given a pair of distance metrics between samples and clusters, the algorithm was initialized with the eight samples forming eight different clusters and then processed iteratively by joining the two most similar clusters. the tree was built starting from the individual samples, using an agglomerative approach. the resulting hierarchy of clusters was displayed as a dendrogram. these traditional clustering methods provide a quick, exploratory overview of the data. however, these methods do not estimate the optimal number of clusters in the data; rather, the clustering is performed exhaustively from the lowest possible level of the hierarchy where each sample forms its own cluster, to the highest level where all samples are grouped into one cluster. in addition to the traditional hierarchical agglomerative clustering method, the hierarchical ordered partitioning and collapsing hybrid algorithm was also applied to the cytokine measurements. in contrast with the previous approaches where the tree was built starting from the individual samples as the leaf nodes, hopach used a hybrid divisive-agglomerative approach: it started from the root cluster containing all the samples, then divided the root down to leaf nodes, with an extra collapsing step after each iteration that combined similar clusters. based on the correlation distance between samples, hopach determined the split that minimized a measure of cluster homogeneity called the median split silhouette. while computationally more expensive than the previous methods, hopach was expected to perform better because of its dynamic approach to update and potentially revise the clusters at every step of the iteration. furthermore, hopach also estimated the optimal number of clusters from the data, and thus offered another advantage over the previous methods. computations were performed in the r computing environment and the hopach package.  cytokine levels were examined using an ex vivo model, termed weem for whole blood x vivoexposure model. individual samples of anti-coagulated human blood were incubated with b. anthracis ames, b. anthracis sterne, y. pestis kim d y. pestis nyc, y. pestis india/p, y. enterocolitica, and y. pseudotuberculosis. after hours exposure, blood cells were removed by low-speed centrifugation and concentrations of cytokines in the plasma were measured with protein arrays. concentrations of fourteen cytokines, gcsf, ifnγ, gm-csf, il- il- il- il- il- il- il- il- mcp- tgfβ, and tnfβ were below the limit of detection in this study. the following cytokines were detected: eotaxin, il- il- il- il-1α, il-1β, il- il- ip- mcp- mig, tnfα, trail, scd scd and sicam-. to determine if there were significant differences among the levels of cytokines in the control and pathogen exposed plasma samples, f-tests were performed. for thirteen of these cytokines, all three replicates were detected and these cytokines were subjected to f-tests. statistical analysis indicated that cytokines had differentially elevated expression profiles following different bacterial exposures. figure shows the concentrations of these cytokines in the control and bacteria exposed plasma samples. the f-tests revealed that the other five cytokines containing complete datasets, trail, scd scd mig, and sicam- had no significant difference between bacterial exposures and the mock-exposed control. moreover, there was a great variation in absolute concentrations between cytokines. for example, the concentrations of tnfα, scd and sicam- were as high as x - pg/ml, whereas il- was much lower, about pg/ml. marked differences in induced cytokine patterns between b. anthracis and yersinia exposures were found. also, the levels of induction of these cytokines differed among the different bacteria. for example, yersinia species induced much higher cytokine response than b. anthracis for il-1α, il-1β, il- and tnf-α. the two strains of b. anthracis bacteria induced different levels of il-1β and tnf-α, including times higher concentration of il-1β and times higher concentration of tnf-α for the sterne strain than the ames strain of b. anthracis. these differences were statistically significant. to discriminate y. pestis exposure from near neighbors, il- levels can be used, showing cytokine concentrations following y. enterocolitica exposure and y. pseudotuberculosis exposure that are on average 5-fold higher and 2-fold higher, respectively, than after y. pestis exposure. il- differential expression was specific to the yersinia spp. because exposure to b. anthracis strains showed comparable il- levels to that in unexposed control. the hopach algorithm estimated the number of clusters as five, and grouped the samples based on their host cytokine expression profiles as follows: 1) y. pestis, 2) y. pseudotuberculosis, 3) y. enterocolitica, 4) b. anthracis, and 5) control. the closer the pathogen-exposed samples are within the tree on the left, the more similar they are. height of the branches indicates the distance between the successive nodes in the clustering. the method separated the b. anthracis and yersinia infected blood samples. in addition, the cytokine profile of the mock-exposed control was more similar to the pattern produced by b. anthracis exposure than to the profile elicited by yersinia. results of the hierarchical clustering when using the euclidean distance between samples depended on the distance metric between clusters. the three methods for determining the distance between clusters all established three major clusters: 1) y. pestis and near neighbors, 2) b. anthracis, and c) control. differences between the results occurred when the yersinia cluster was further divided. the average linkage method, consistent with figure formed a subgroup of the three y. pestis strains, then grouped them first with y. pseudotuberculosis followed by y. enterocolitica. complete and single linkage methods, however, first grouped the attenuated virulent strain of y. pestis with the more virulent strain, both clinical isolates from human plague cases, and then clustered them with y. pseudotuberculosis, followed by the attenuated y. pestis, and lastly with y. enterocolitica. this is interesting from an evolutionary perspective because it has been proposed that y. pestis evolved from y. pseudotuberculosis within the last years, and thus these two pathogens are more closely related. when using hierarchical clustering with the correlation distance between the samples, the final clusters were independent of the distance metric between clusters, and agreed with the tree structure in figure the complete, single, and average linkage methods all resulted in the following major clusters: 1) yersinia, 2) b. anthracis, and 3). control. within the yersinia cluster, y. pestis was closest to y. pestis, followed by y. pestis, y. pseudotuberculosis, and y. enterocolitica. discussion the hopach clustering method produced five distinctly separated clusters: 1) y. pestis, 2) y. pseudotuberculosis, 3) y. enterocolitica, 4) b. anthracis, and 5) control. this result is consistent with the findings using the correlation distance and the euclidean distance with average linkage. in addition, hopach estimated the optimal number of clusters as five. that is, the yersinia subcluster is best if it is divided into the three clusters specified by 1) through 5) above. y. enterocolitica forms its own cluster, and so does y. pseudotuberculosis. y. pestis, y. pestis, and y. pestis are grouped into one cluster. further subdivisions lead to an overall clustering with inferior quality. in addition to clustering the cytokine expression profiles across bacterial treatments, figure also groups the cytokines themselves and clusters the proteins based on their similarities across the pathogen exposures and reorders them accordingly. interestingly, the three pro-inflammatory cytokines il-1β, tnfα, and il- clustered closely, and so did the three chemokines mcp- ip- and il- although these cytokines do not cluster as a single group, they do cluster at a branch further away from the leaf node, which includes il- and scd to make a larger group of proteins. several of these proteins are involved in inflammatory conditions, such as il-1beta, tnfα, il- and have been shown to be upregulated in cell culture and animal model specifically exposed to biothreat agents. increased expression of il- and tnfα clustered together in a study involving mouse splenic cd11b + cells following sub-lethal y. enterocolitica infection. in addition, several of the cytokines in this cluster, namely tnf-alpha, il1-beta, il- and mcp- are expressed higher in exposed whole blood as compared to control in this study and in whole blood exposure to lps from several other gram negative bacterial pathogens. in addition to expression differences, the absence of detected cytokine expression can also be helpful in discriminating pathogen exposure. the multiplex detection of cytokines in this study revealed the early phase cytokine expression profiles in human plasma following exposures to b. anthracis, y. pestis, y. pseudotuberculosis, and y. enterocolitica. the expression levels of cytokines, il-1α, il-1β, il- il- il- ip- mcp- and tnfα were significantly different from that of unexposed control. although the focus of our work was to show that cytokine expression profiling can discriminate between different pathogen exposures in a human whole blood ex vivo model, these results also represent an initial attempt to characterize the full cytokine response to each individual pathogen. our preliminary study using a single exposure protocol at a single time post-exposure will need to be supplemented with more thorough investigation in order to determine the usefulness of using cytokine levels for diagnosing pathogen exposure. however, the single time point chosen, hours, is sufficient to detect proteomic changes and has been used in previous studies examining cytokine levels. this time point represents a start towards a more complete temporal study, as has been done with gene expression patterns for two of the pathogens studied here. in addition, studies that provide expression patterns for a single cytokine using multiple time points will also be needed to make the results of this paper clinically useful, such as has been done by, cooper and coworkers, who examined il-12p and il-12p levels following different growth conditions and exposure levels for a time course of y. pestis exposed dendritic cells. the results of the current work shows a similar expression pattern trend to this previous work, in which, y. pestis induces il-12p and at a substantially higher level than il-12p our results showed that the expression levels of chemokines, il- mcp- and ip- were induced by both yersinia and b. anthracis exposures. no significant differences were found for these cytokines between yersinia and b. anthracis exposures. il- mcp- and ip- are chemokines that enable the migration of leukocytes from the blood to the site of inflammation. il- is a key chemokine regulating neutrophil recruitment. the essential involvement of il- in acute inflammation was demonstrated by neutralizing il- with its antibody. when highly specific antibody against il- was administered in acute inflammatory reactions induced by several stimuli including lipopolysaccharide, neutrophil infiltration was blocked. mcp- is known for its ability to act as potent chemoattractant and activator of monocytes/macrophages as well as nk cells but not neutrophils. ip- has no chemotactic activity for neutrophils but attracts monocytes, nk, and t cells to the site of infection and regulates t cell maturation. it was reported previously that elevated il- and mcp- were secreted by human epithelial cells after y. enterocolitica infection, but not ip-. human dendritic cells, infected with b. anthracis spores, secreted high level of il- at hours. in our study, the fold increase of il- was much greater than mcp- and ip-. for example, the induction of il- by ames strain of b. anthracis was fold, while mcp- was fold and ip- was fold. this result may indicate that il- is a dominant chemokine in early response and neutrophils are the major player in early inflammatory response. here we compared cytokines induced by b. anthracis and yersinia exposures. overall, yersinia exposure induced higher levels of il-1α, il-1β, il- il- and tnfα than b. anthracis exposure, suggesting these cytokines could be used to develop an assay for discriminating yersinia spp. from b. anthracis exposures. the vaccine strain of b. anthracis induced higher levels of il-1β and tnfα than the virulent strain, suggesting these cytokines can contribute to a biomarker panel to discriminate if a particular isolate of b. anthracis is virulent. there was also a difference in induction of il- between y. pestis and near neighbors, suggesting this cytokine is a candidate biomarker for discriminating the virulence of yersinia species. these data regarding il- expression following yersinia spp. exposure are in agreement with published literature that shows y. enterocolitica and y. pestis can elicit statistically different levels of il- expression. differences in il- induction may be due to differences in the lcrv protein among yersinia spp.. the different cytokine profiles induced by b. anthracis and yersinia here may be partially due to different surface antigens on the outermost part of these pathogens and the manner in which these bacteria were grown. lipopolysaccharide, the main constituent of the outer membrane of gram-negative bacteria, and peptidoglycan, the major cell wall component of gram-positive bacteria, have been reported to elicit markedly different immune responses. however, virulence factors, such as b. anthracis lethal toxin and yersinia virulence antigen, lcrv, may also play important roles in differential cytokine induction. this view is supported by numerous reports that b. anthracis toxin and virulence factors of yersinia bacteria modulate host cytokine responses. while the various clustering methods resulted in slightly different final hierarchies, all were consistent in separating the unexposed control from the samples exposed to b. anthracis or to the y. pestis and near neighbors. agreement on this level among the various clustering procedures lends more confidence to the overall results. on a more detailed level, the methods grouped slightly differently the samples exposed to the y. pestis and near neighbors, which indicates that these samples cannot be unequivocally separated based on the current data and additional biomarkers or a larger sample set would be needed. the most advanced hopach method estimated the optimal number of clusters in the data as five, corresponding to the unexposed control, and the four species: b. anthracis, y. pseudotuberculosis, y. enterocolitica, and y. pestis. information gained from the targeted protein array data for host response complements genomic, and other proteomic studies of host-pathogen interactions. the success of the weem and computational method to distinguish pathogen exposure, based on host response in this initial study, is encouraging and suggests a number of possibilities for future studies to refine the findings. comparative analysis, such as the current work, can potentially reveal the critical pathogenic mechanism and host innate immune responses during infection as was previously shown for y. pestis and y. pseudotuberculosis. opportunities include using statistical hypothesis tests based on analysis of variance to assess the significance of the observed differences among the host-pathogen cytokine concentration profiles, as well as performing follow-up studies to focus more on the y. pestis and near neighbor cluster. in addition, the methods can be extended to investigate host responses to diverse pathogens in multiple host model systems to cross validate the significance of the biomarkers to distinguish pathogen exposures.  results from this study suggest that cytokine arrays coupled with statistical clustering methods can distinguish exposures to pathogens, including multiple strains of y. pestis, y. pseudotuberculosis, y. enterocolitica, and b. anthracis. these methods differentiate both near neighbors and distant evolutionary microbes based on host response data. the distinct cytokine profiles also provide insight into both the host response and virulence mechanisms of diverse pathogens. in summary, characterization of host responses based on cytokine profiles has translational application, potentially providing the identification of infectious diseases and leading toward the ultimate goal of presymptomatic detection via sentinel surveillance of pathogen exposure and appropriate treatment. competing interests the authors declare that they have no competing interests. authors’ contributions bc participated in the design of the study, conducted experiments, and drafted and finalized the manuscript. km conducted experiments. if participated in the design of the study, performed the statistical analysis, and helped draft the manuscript. pl helped draft the manuscript. sm conceived of the study, obtained funding, and participated in its design and coordination and helped to draft the manuscript. all authors read and approved the final manuscript. acknowledgments we thank david callender, jonathan e. forman, and renee tobias from zyomyx for their assistance with the biochip analyses. we thank todd corzett for useful comments on the manuscript and brent ricks for technical assistance. this work was performed under the auspices of the us department of energy by the university of california, lawrence livermore national laboratory under contract no. w-7405-eng- with support from the department of homeland security. the authors would also like to thank psw rce animal resources and laboratory services core u54-ai ucrl-jrnl- 
toward this end, cytokine profiles in human blood were investigated using a human whole blood ex vivo exposure model, called weem. samples of whole blood from healthy volunteers were incubated with seven pathogens including yersinia pseudotuberculosis, yersinia enterocolitica, bacillus anthracis, and multiple strains of yersinia pestis, and multiplexed protein expression profiling was conducted on supernatants of these cultures with an antibody array to detect cytokines simultaneously. cluster analysis, based on host response cytokine profiles, indicates that distinct patterns of immunomodulatory proteins are induced by the different pathogen exposures and these patterns may enable further development into biomarkers for diagnosing pathogen exposure. levels of cytokines, il-1α, il-1β, il- il- il- ip- mcp- and tnfα, were significantly up-regulated in plasma after bacterial exposures of hours. in addition, the y. pestis and yersinia near neighbors were differentiated from the b. anthracis strains. rapid detection and therapeutic intervention for infectious and emerging diseases is a major scientific goal in biodefense and public health.
 rapid detection and therapeutic intervention for infectious and emerging diseases is a major scientific goal in biodefense and public health. toward this end, cytokine profiles in human blood were investigated using a human whole blood ex vivo exposure model, called weem.  samples of whole blood from healthy volunteers were incubated with seven pathogens including yersinia pseudotuberculosis, yersinia enterocolitica, bacillus anthracis, and multiple strains of yersinia pestis, and multiplexed protein expression profiling was conducted on supernatants of these cultures with an antibody array to detect cytokines simultaneously. levels of cytokines, il-1α, il-1β, il- il- il- ip- mcp- and tnfα, were significantly up-regulated in plasma after bacterial exposures of hours. statistical clustering was applied to group the pathogens based on the host response protein expression profiles. the nearest phylogenetic neighbors clustered more closely than the more distant pathogens, and all seven pathogens were clearly differentiated from the unexposed control. in addition, the y. pestis and yersinia near neighbors were differentiated from the b. anthracis strains.  cluster analysis, based on host response cytokine profiles, indicates that distinct patterns of immunomodulatory proteins are induced by the different pathogen exposures and these patterns may enable further development into biomarkers for diagnosing pathogen exposure. 