IL-17A produced by neutrophils protects against pneumonic plague through orchestrating IFN-γ-activated macrophage programming.

Innate immune cells, including neutrophils and macrophages, are critically involved in host antimicrobial defense responses. Intrinsic regulatory mechanisms controlling neutrophil and macrophage activities are poorly defined. In this study, we found that IL-17A, a natural signal factor, could provide protection against early pneumonic plague inflammation by coordinating the functions of neutrophils and programming of macrophages. The IL-17A level is promptly increased during the initial infection. Importantly, abrogation of IL-17A or IL-17AR significantly aggravated the infection, but mIL-17A treatment could significantly alleviate inflammatory injury, revealing that IL-17A is a critical requirement for early protection of infection. We also demonstrated that IL-17A was predominantly produced by CD11b(+)Ly6G(+) neutrophils. Although IL-17A could not significantly affect the antimicrobial responses of neutrophils, it could target the proinflammatory macrophage (M1) programming and potentiate the M1's defense against pneumonic plague. Mechanistically, IFN-γ treatment or IFN-γ-activated M1 macrophage transfer could significantly mitigate the aggravated infection of IL-17A(-/-) mice. Finally, we showed that IL-17A and IFN-γ could synergistically promote macrophage anti-infection immunity. Thus, our findings identify a previously unrecognized function of IL-17A as an intrinsic regulator in coordinating neutrophil and macrophage antimicrobial activity to provide protection against acute pneumonic plague.

Historical variations in mutation rate in an epidemic pathogen, Yersinia pestis.

The genetic diversity of Yersinia pestis, the etiologic agent of plague, is extremely limited because of its recent origin coupled with a slow clock rate. Here we identified 2,326 SNPs from 133 genomes of Y. pestis strains that were isolated in China and elsewhere. These SNPs define the genealogy of Y. pestis since its most recent common ancestor. All but 28 of these SNPs represented mutations that happened only once within the genealogy, and they were distributed essentially at random among individual genes. Only seven genes contained a significant excess of nonsynonymous SNP, suggesting that the fixation of SNPs mainly arises via neutral processes, such as genetic drift, rather than Darwinian selection. However, the rate of fixation varies dramatically over the genealogy: the number of SNPs accumulated by different lineages was highly variable and the genealogy contains multiple polytomies, one of which resulted in four branches near the time of the Black Death. We suggest that demographic changes can affect the speed of evolution in epidemic pathogens even in the absence of natural selection, and hypothesize that neutral SNPs are fixed rapidly during intermittent epidemics and outbreaks.

Insight into bacterial virulence mechanisms against host immune response via the Yersinia pestis-human protein-protein interaction network.

A Yersinia pestis-human protein interaction network is reported here to improve our understanding of its pathogenesis. Up to 204 interactions between 66 Y. pestis bait proteins and 109 human proteins were identified by yeast two-hybrid assay and then combined with 23 previously published interactions to construct a protein-protein interaction network. Topological analysis of the interaction network revealed that human proteins targeted by Y. pestis were significantly enriched in the proteins that are central in the human protein-protein interaction network. Analysis of this network showed that signaling pathways important for host immune responses were preferentially targeted by Y. pestis, including the pathways involved in focal adhesion, regulation of cytoskeleton, leukocyte transendoepithelial migration, and Toll-like receptor (TLR) and mitogen-activated protein kinase (MAPK) signaling. Cellular pathways targeted by Y. pestis are highly relevant to its pathogenesis. Interactions with host proteins involved in focal adhesion and cytoskeketon regulation pathways could account for resistance of Y. pestis to phagocytosis. Interference with TLR and MAPK signaling pathways by Y. pestis reflects common characteristics of pathogen-host interaction that bacterial pathogens have evolved to evade host innate immune response by interacting with proteins in those signaling pathways. Interestingly, a large portion of human proteins interacting with Y. pestis (16/109) also interacted with viral proteins (Epstein-Barr virus [EBV] and hepatitis C virus [HCV]), suggesting that viral and bacterial pathogens attack common cellular functions to facilitate infections. In addition, we identified vasodilator-stimulated phosphoprotein (VASP) as a novel interaction partner of YpkA and showed that YpkA could inhibit in vitro actin assembly mediated by VASP.

A dog-associated primary pneumonic plague in Qinghai Province, China.

BACKGROUND: Primary pneumonic plague (PPP) caused by Yersinia pestis is the most threatening clinical form of plague. An outbreak was reported in July 2009 in Qinghai Province, China. METHODS: This outbreak was investigated by clinical, epidemiological, bacteriological, and immunological methods. Multilocus variable number tandem repeat analysis (MLVA) was used to track the source of the outbreak. RESULTS: The index case, a patient with PPP, contaminated 11 close contacts. All the 12 cases, including the index patient, experienced sudden onset of fever, headache, and productive coughing with bloody sputum. Three of them died. Nevertheless, another 61 direct and 256 indirect contacts were not infected during the 2-week quarantine. Antibodies to F1 antigen were detected in 9 survival cases, with a 4-fold increase in titers in serum samples collected at different periods. Seven strains of Y. pestis were isolated from dogs and patients. Field investigation and MLVA of the isolated strains revealed that this outbreak was started by a deceased dog. CONCLUSION: Dogs are believed to be an indicator animal for plague surveillance, but their association with PPP is rare. Our results provide evidence for this possibility, which suggests the public health significance of dogs as a source of plague.

The complete genome sequence and proteomics of Yersinia pestis phage Yep-phi.

Yep-phi, a lytic phage of Yersinia pestis, was isolated in China and is routinely used as a diagnostic phage for the identification of the plague pathogen. Yep-phi has an isometric hexagonal head containing dsDNA and a short non-contractile conical tail. In this study, we sequenced the Yep-phi genome (GenBank accession no. HQ333270) and performed proteomics analysis. The genome consists of 38 ,616 bp of DNA, including direct terminal repeats of 222 bp, and is predicted to contain 45 ORFs. Most structural proteins were identified by proteomics analysis. Compared with the three available genome sequences of lytic phages for Y. pestis, the phages could be divided into two subgroups. Yep-phi displays marked homology to the bacteriophages Berlin (GenBank accession no. AM183667) and Yepe2 (GenBank accession no. EU734170), and these comprise one subgroup. The other subgroup is represented by bacteriophage ΦA1122 (GenBank accession no. AY247822). Potential recombination was detected among the Yep-phi subgroup.

Phenotypic and transcriptional analysis of the osmotic regulator OmpR in Yersinia pestis.

BACKGROUND: The osmotic regulator OmpR in Escherichia coli regulates differentially the expression of major porin proteins OmpF and OmpC. In Yersinia enterocolitica and Y. pseudotuberculosis, OmpR is required for both virulence and survival within macrophages. However, the phenotypic and regulatory roles of OmpR in Y. pestis are not yet fully understood. RESULTS: Y. pestis OmpR is involved in building resistance against phagocytosis and controls the adaptation to various stressful conditions met in macrophages. The ompR mutation likely did not affect the virulence of Y. pestis strain 201 that was a human-avirulent enzootic strain. The microarray-based comparative transcriptome analysis disclosed a set of 224 genes whose expressions were affected by the ompR mutation, indicating the global regulatory role of OmpR in Y. pestis. Real-time RT-PCR or lacZ fusion reporter assay further validated 16 OmpR-dependent genes, for which OmpR consensus-like sequences were found within their upstream DNA regions. ompC, F, X, and R were up-regulated dramatically with the increase of medium osmolarity, which was mediated by OmpR occupying the target promoter regions in a tandem manner. CONCLUSION: OmpR contributes to the resistance against phagocytosis or survival within macrophages, which is conserved in the pathogenic yersiniae. Y. pestis OmpR regulates ompC, F, X, and R directly through OmpR-promoter DNA association. There is an inducible expressions of the pore-forming proteins OmpF, C, and × at high osmolarity in Y. pestis, in contrast to the reciprocal regulation of them in E. coli. The main difference is that ompF expression is not repressed at high osmolarity in Y. pestis, which is likely due to the absence of a promoter-distal OmpR-binding site for ompF.

Regulatory effects of cAMP receptor protein (CRP) on porin genes and its own gene in Yersinia pestis.

BACKGROUND: The cAMP receptor protein (CRP) is a global bacterial regulator that controls many target genes. The CRP-cAMP complex regulates the ompR-envZ operon in E. coli directly, involving both positive and negative regulations of multiple target promoters; further, it controls the production of porins indirectly through its direct action on ompR-envZ. Auto-regulation of CRP has also been established in E. coli. However, the regulation of porin genes and its own gene by CRP remains unclear in Y. pestis. RESULTS: Y. pestis employs a distinct mechanism indicating that CRP has no regulatory effect on the ompR-envZ operon; however, it stimulates ompC and ompF directly, while repressing ompX. No transcriptional regulatory association between CRP and its own gene can be detected in Y. pestis, which is also in contrast to the fact that CRP acts as both repressor and activator for its own gene in E. coli. It is likely that Y. pestis OmpR and CRP respectively sense different signals (medium osmolarity, and cellular cAMP levels) to regulate porin genes independently. CONCLUSION: Although the CRP of Y. pestis shows a very high homology to that of E. coli, and the consensus DNA sequence recognized by CRP is shared by the two bacteria, the Y. pestis CRP can recognize the promoters of ompC, F, and X directly rather than that of its own gene, which is different from the relevant regulatory circuit of E. coli. Data presented here indicate a remarkable remodeling of the CRP-mediated regulation of porin genes and of its own one between these two bacteria.

Use of protein microarray to identify gene expression changes of Yersinia pestis at different temperatures.

Yersinia pestis is a bacterium that is transmitted between fleas, which have a body temperature of 26 °C, and mammalian hosts, which have a body temperature of 37 °C. To adapt to the temperature shift, phenotype variations, including virulence, occur. In this study, an antigen microarray including 218 proteins of Y. pestis was used to evaluate antibody responses in a pooled plague serum that was unadsorbed, adsorbed by Y. pestis cultivated at 26 °C, or adsorbed by Y. pestis cultivated at 26 and 37 °C to identify protein expression changes during the temperature shift. We identified 12 proteins as being expressed at 37 °C but not at 26 °C, or expressed at significantly higher levels at 37 °C than at 26 °C. The antibodies against 7 proteins in the serum adsorbed by Y. pestis cultivated at 26 and 37 °C remained positive, suggesting that they were not expressed on the surface of Y. pestis in LB broth in vitro or specifically expressed in vivo. This study proved that protein microarray and antibody profiling comprise a promising technique for monitoring gene expression at the protein level and for better understanding pathogenicity, to find new vaccine targets against plague.

Transcriptional regulation of ompF2, an ompF paralogue, in Yersinia pestis.

A regulatory circuit composed of three porins (OmpF, OmpC, and OmpX) and two transcriptional regulators (OmpR and CRP) has previously been characterized in Yersinia pestis . In this follow-up study, OmpF2, an OmpF paralogue, was integrated into this regulatory circuit. Only basal expression was detected for ompF2 in the wild-type strain under different osmotic conditions. The ompF2 transcription was dramatically enhanced with increasing medium osmolarity in the ompR null mutant background. The CRP regulator had no regulatory effect on ompF2 under the growth conditions tested.

[Purification of recombinant H-NS protein of Yersinia pestis and characterization of its DNA-binding activity].

OBJECTIVE: The regulator protein H-NS of Yersinia pestis was expressed using the Escherichia coli BL21lambdaDE3 protein expression system, and its DNA-binding activity was characterized. METHODS: The entire coding region of the hns gene was amplified by PCR from Y. pestis strain 201, and then cloned into the BamHI and Sal sites of the vector pET28a. The recombinant plasmid pET28a-hns was transformed into BL21lambdaDE3. Over-expression of His-H-NS in the LB medium was induced by addition of 1 mM IPTG (isopropyl-b-D-thiogalactoside). The over-expressed protein was purified under native conditions with nickel loaded HiTrap Chelating Sepharose columns (Amersham). The electrophoretic mobility shift assay and DNase I footprinting experiments were carried out to analyze the DNA-binding activity of His-H-NS in vitro. RESULTS: The purified His-H-NS protein could bind to the upstream DNA regions of psaA, psaE and rovA of Y. pestis, and the H-NS-binding sites were determined at the single nucleotide resolution. CONCLUSION: The purified His-H-NS protein could bind to target DNA fragments, suggesting that H-NS would regulate the transcription of relevant genes in Y. pests.

Different strategies for preparation of non-tagged rV270 protein and its efficacy against Yersinia pestis challenge.

OBJECTIVE: LcrV is an important component for the development of a subunit vaccine against plague. To reduce immunosuppressive activity of LcrV, a recombinant LcrV variant lacking amino acids 271 to 326 (rV270) was prepared by different methods in this study. METHODS: A new strategy that produced non-tagged or authentic rV270 protein was designed by insertion of rV270-thrombin-hexahistidine fusion gene into the vector pET24a, or by insertion of hexahistidine-enterokinase-rV270 or hexahistitine-factor Xa-rV270 fusion gene into the vector pET32a. After Co(2+) affinity chromatography, a purification strategy was developed by cleavage of His tag on column, following Sephacryl S-200HR column filtration chromatography. RESULTS: Removal of His tag by thrombin, enterokinase and factor Xa displayed a yield of 99.5%, 32.4% and 15.3%, respectively. Following Sephacryl S-200HR column filtration chromatography, above 97% purity of rV270 protein was obtained. Purified rV270 that was adsorbed to 25% (v/v) Al(OH)3 adjuvant in phosphate-buffered saline (PBS) induced very high titers of antibody to rV270 in BALB/c mice and protected them (100% survival) against subcutaneous challenge with 106 CFU of Y. pestis virulent strain 141. CONCLUSION: The completely authentic rV270 protein can be prepared by using enterokinase or factor Xa, but they exhibited extremely low cleavage activity to the corresponding recognition site. Thrombin cleavage is an efficient strategy to prepare non-tagged rV270 protein and can be easily operated in a large scale due to its relatively low cost and high cleavage efficacy. The recombinant rV270 can be used as a key component to develop a subunit vaccine of plague.

Draft genome sequences of Yersinia pestis isolates from natural foci of endemic plague in China.

To gain insights into the evolutionary origin, emergence, and pathogenicity of the etiologic agent of plague, we have sequenced the genomes of four Yersinia pestis strains isolated from the zoonotic rodent reservoir in foci of endemic plague in China. These resources enable in-depth studies of Y. pestis sequence variations and detailed whole-genome comparisons of very closely related genomes from the supposed site of the origin and the emergence of global pandemics of plague.

High-throughput identification of new protective antigens from a Yersinia pestis live vaccine by enzyme-linked immunospot assay.

Yersinia pestis, the plague pathogen, is a facultative intracellular bacterium. Cellular immunity plays important roles in defense against infections. The identification of T-cell targets is critical for the development of effective vaccines against intracellular bacteria; however, the function of cellular immunity in protection from plague was not clearly understood. In this study, 261 genes from Y. pestis were selected on the basis of bioinformatics analysis and previous research results for expression in Escherichia coli BL21(DE3). After purification, 101 proteins were qualified for examination of their abilities to induce the production of gamma interferon in mice immunized with live vaccine EV76 by enzyme-linked immunospot assay. Thirty-four proteins were found to stimulate strong T-cell responses. The protective efficiencies for 24 of them were preliminarily evaluated using a mouse plague model. In addition to LcrV, nine proteins (YPO0606, YPO1914, YPO0612, YPO3119, YPO3047, YPO1377, YPCD1.05c, YPO0420, and YPO3720) may provide partial protection against challenge with a low dose (20 times the 50% lethal dose [20x LD(50)]) of Y. pestis, but only YPO0606 could partially protect mice from infection with Y. pestis at a higher challenge dosage (200x LD(50)). These proteins would be the potential components for Y. pestis vaccine development.

Yersinia pestis and host macrophages: immunodeficiency of mouse macrophages induced by YscW.

The virulence of the pathogenic Yersinia species depends on a plasmid-encoded type III secretion system (T3SS) that transfers six Yersinia outer protein (Yop) effector proteins into the cytoplasm of eukaryotic cells, leading to disruption of host defence mechanisms. It is shown in this study that Yersinia pestis YscW, a protein of the T3SS injectisome, contributes to the induction of a deficiency in phagocytosis in host macrophages and a reduction in their antigen-presenting capacity. A Y. pestis strain lacking yscW had no effect on uptake by host macrophages. In mice infected with wild-type Y. pestis, the yscW mutant or a complement strain, immunodeficiency was observed in host macrophages compared with those from uninfected mice. However, the phagocytosis and antigen presenting capacities of macrophages infected by yscW mutant strain both in vivo and in vitro were significantly higher than those by wild type strain. Consistent with this finding, when YscW was expressed in the RAW264.7 macrophage cell line, phagocytosis and antigen-presenting capacities were significantly lower than those of the control groups. These results indicate that Y. pestis YscW may directly induce immunodeficiency in murine macrophages by crippling their phagocytosis and antigen-presenting capacities. These data provide evidences to Y. pestis pathogenesis that some proteins in T3SS injectisome, such as YscW protein, might play independent roles in disrupting host defense apart from their known functions.

Gene expression profiling of Yersinia pestis with deletion of lcrG, a known negative regulator for Yop secretion of type III secretion system.

Yersinia pestis injects a set of virulent proteins into the cytosol of eukaryotic cells by a type III secretion system (T3SS). LcrG is a known negative regulator for secretion of Yersinia outer-membrane proteins (Yops) by blocking the secretion apparatus (Ysc) from the inner membrane. To further understand the effect of lcrG deletion on Y. pestis T3SS regulation, transcriptional profiles from the DeltalcrG mutant and wild-type Y. pestis strains were compared. The results showed that although the DeltalcrG mutant was markedly attenuated (600-fold increase of LD(50) in s.c. challenged BALB/c mice), transcriptions of almost all the type III genes were upregulated significantly in the DeltalcrG mutant. The immunoblotting analysis of YopM and LcrV demonstrated that their expressions were also increased in the DeltalcrG mutant in comparison to the wild-type strain. We speculate that, in addition to the negative regulation of the Yop secretion, LcrG could possibly play a negative regulatory role in the transcription of T3SS genes through indirect mechanisms. Furthermore, this report also revealed significant transcriptional changes in the genes encoding cell-envelope-related proteins and a virulence-related transcription factor RovA in the DeltalcrG mutant.

Characterization of Zur-dependent genes and direct Zur targets in Yersinia pestis.

BACKGROUND: The zinc uptake regulator Zur is a Zn2+-sensing metalloregulatory protein involved in the maintenance of bacterial zinc homeostasis. Up to now, regulation of zinc homeostasis by Zur is poorly understood in Y. pestis. RESULTS: We constructed a zur null mutant of Y. pestis biovar microtus strain 201. Microarray expression analysis disclosed a set of 154 Zur-dependent genes of Y. pestis upon exposure to zinc rich condition. Real-time reverse transcription (RT)-PCR was subsequently used to validate the microarray data. Based on the 154 Zur-dependent genes, predicted regulatory Zur motifs were used to screen for potential direct Zur targets including three putative operons znuA, znuCB and ykgM-RpmJ2. The LacZ reporter fusion analysis verified that Zur greatly repressed the promoter activity of the above three operons. The subsequent electrophoretic mobility shift assay (EMSA) demonstrated that a purified Zur protein was able to bind to the promoter regions of the above three operons. The DNase I footprinting was used to identify the Zur binding sites for the above three operons, verifying the Zur box sequence as predicted previously in gamma-Proteobacteria. The primer extension assay was further used to determine the transcription start sites for the above three operons and to localize the -10 and -35 elements. Zur binding sites overlapped the -10 sequence of its target promoters, which was consistent with the previous observation that Zur binding would block the entry of the RNA polymerase to repress the transcription of its target genes. CONCLUSION: Zur as a repressor directly controls the transcription of znuA, znuCB and ykgM-RpmJ2 in Y. pestis by employing a conserved mechanism of Zur-promoter DNA association as observed in gamma-Proteobacteria. Zur contributes to zinc homeostasis in Y. pestis likely through transcriptional repression of the high-affinity zinc uptake system ZnuACB and two alternative ribosomal proteins YkgM and RpmJ2.

Direct and negative regulation of the sycO-ypkA-ypoJ operon by cyclic AMP receptor protein (CRP) in Yersinia pestis.

BACKGROUND: Pathogenic yersiniae, including Y. pestis, share a type III secretion system (T3SS) that is composed of a secretion machinery, a set of translocation proteins, a control system, and six Yop effector proteins including YpkA and YopJ. The cyclic AMP receptor protein (CRP), a global regulator, was recently found to regulate the laterally acquired genes (pla and pst) in Y. pestis. The regulation of T3SS components by CRP is unknown. RESULTS: The sycO, ypkA and yopJ genes constitute a single operon in Y. pestis. CRP specifically binds to the promoter-proximate region of sycO, and represses the expression of the sycO-ypkA-yopJ operon. A single CRP-dependent promoter is employed for the sycO-ypkA-yopJ operon, but two CRP binding sites (site 1 and site 2) are detected within the promoter region. A CRP box homologue is found in site 1 other than site 2. The determination of CRP-binding sites, transcription start site and core promoter element (-10 and -35 regions) promotes us to depict the structural organization of CRP-dependent promoter, giving a map of CRP-promoter DNA interaction for sycO-ypkA-yopJ. CONCLUSION: The sycO-ypkA-yopJ operon is under the direct and negative regulation of CRP in Y. pestis. The sycO-ypkA-yopJ promoter-proximate regions are extremely conserved in Y. pestis, Y. pseudotuberculosis and Y. enterocolitica. Therefore, data presented here can be generally applied to the above three pathogenic yersiniae.

Simultaneous detection of five biothreat agents in powder samples by a multiplexed suspension array.

A suspension array-based multiplexed immunoassay was developed for rapid, sensitive, specific, and simultaneous detection of multiple biothreat-associated agents in powder samples. The 5-plexed immunoassays using sets of 9-plexed coupled fluorescent beads were employed to simultaneously detect five representative biothreat agents, including B. anthracis spore, Y. pestis, SARS-CoV, staphylococcal enterotoxin B (SEB) and ricin from a single powder sample and the feasibility for field samples was demonstrated by both blinded and standard laboratory trials. The detection sensitivity and dynamic range for the five biothreat agents from different powders might be varied depending on the nature of the powder and the feature of the contaminating agent. The limit of detection for Y. pestis, B. anthracis spores, SEB, ricin, SARS-CoV N protein in milk powder was 20 cfu, 111 cfu, 110pg, 5.4 ng and 2 ng per test respectively. Compared to conventional ELISA method, the suspension array has a higher sensitive ability, and can detect five biothreat agents simultaneously with high reproducibility.

Transcriptional profiling of a mice plague model: insights into interaction between Yersinia pestis and its host.

Despite the importance of pneumonic plague caused by Yersinia pestis, a few is known about the interaction between Y. pestis and its host at the molecular level during the pneumonic plague development. In this study, we employed an intranasally challenged plague model in mice for investigating the kinetics of the disease progression by transcriptional profiling of Y. pestis and mice using qRT-PCR and microarray, respectively. The increasing transcription of important virulence genes of Y. pestis and of mice genes involving in immune and inflammatory defensive responses, and responses to stimuli, presents an overview of interaction between Y. pestis and mice during development of pneumonic plague. The early and persisting up-regulation of caf 1, psa A and lcr V in vivo indicated their role in resisting the host innate immune responses. The up-regulation of fur, ybt A and hms H in vivo reflected the ability of Y. pestis for acquiring iron. The transcription regulators, including pho P, oxy R and omp R, were up-regulated during plague development, suggesting their roles in interaction between Y. pestis and mice. Many genes encoding cytokines, such as IL2, IL-1B, CXCL2, CXCL5, CCL20, CD14 and TNFRSF13B, were up-regulated during the infection, confirming the report that they are important mediators to activate host responses to invading pathogens. The up-regulation of some genes encoding important virulent factors of Y. pestis and expression alterations of some genes encoding cytokines in the host reflect the interaction between the pathogen and the host, which will help us better understand plague pathogenesis.

[Evaluation of immunization protection efficacy of plague subunit vaccine].

OBJECTIVE: To evaluate the protective efficacy of plague subunit vaccine, BALB/c mice, guinea pigs and rabbits were used in this study. METHODS: Groups of mice (10 per group), guinea pigs (14 per group) and rabbits (6 per group) were immunized with F1 + rV270 vaccine, EV76 vaccine and alum adjuvant by intramuscular route, respectively. Serum antibody titres of mice, guinea pigs and rabbits were determined by ELISA and the immunized animals were challenged with 10(6) CFU of Y. pestis strain 141 at the 8th week after the primary immunization. RESULTS: The immunized mice, guinea pigs or rabbits with subunit vaccine developed anti-F1 IgG titre of 41 587.3 +/- 2.1, 11 543.7 +/- 2.1 or 522.4 +/- 22.4 and elicited statistical anti-F1 IgG titre difference among them (F = 17.58, P < 0.01). The immunized mice, guinea pigs or rabbits with subunit vaccine had anti-rV270 IgG titre of 15 748.7 +/- 1.6, 12.6 +/- 1.4 or 1648.0 +/- 5.0 and induced statistical anti-rV270 IgG titre difference among them (F value was 16.34, P < 0.01). There was significant anti-F1 IgG titre difference among mice, guinea pigs and rabbits immunized with EV76 vaccine that developed anti-F1 IgG titre of 913.4 +/- 4.5, 937.0 +/- 2.0 or 342.0 +/- 12.0 (F = 23.67, P < 0.01), whereas the immunized mice, guinea pigs and rabbits with EV76 vaccine developed anti-rV270 IgG titre of 12.0 +/- 1.0, 447.0 +/- 10.0, 40.0 +/- 11.0 and there was no anti-rV270 IgG titre difference between them (F = 2.20, P = 0.1314). The immunized mice with subunit vaccine developed significantly higher anti-F1 IgG titres than immunized guinea pigs and rabbits (q value was 30.57 and 19.04, respectively, P < 0.01), and there were no anti-F1 IgG titre differences between the immunized guinea pigs and rabbits (q = 0.04, P = 0.8485). The immunized mice with subunit vaccine developed significantly higher anti-rV270 IgG titres than immunized guinea pigs and rabbits (q value was 27.10 and 19.49, respectively, P < 0.01), and there were no anti-rV270 IgG titre differences between the immunized guinea pigs and rabbits with the subunit vaccine (q = 0.25, P = 0.6187). The immunized mice with EV76 elicited higher anti-F1 IgG titres than immunized guinea pigs and rabbits (q value was 40.67 and 29.10, respectively, P < 0.01), whereas there was no difference of F1 IgG titer between immunized guinea pigs and rabbits (q = 0.06, P = 0.8098). The immunized mice, guinea pigs and rabbits with subunit vaccine provided 100% (10/10), 86% (12/14) and 100% (5/5) protection against 10(6) CFU Y. pestis of challenge, respectively. The immunized mice, guinea pigs and rabbits with EV76 vaccine gave 100% (6/6), 93% (13/14) and 100% (6/6) protection against 10(6) CFU Y. pestis of challenge respectively. CONCLUSION: BALB/c mice is the best small animal model for valuation of protective efficacy of plague subunit vaccine. The guinea pigs showed a high individual variation for this purpose. The rabbits can be used as an alternative model for evaluating plague subunit vaccine.