Previous Article | Next Article ![]()
Journal of Clinical Microbiology, June 2006, p. 2007-2018, Vol. 44, No. 6
0095-1137/06/$08.00+0 doi:10.1128/JCM.02630-05
Copyright © 2006, American Society for Microbiology. All Rights Reserved.
National Microbiology Laboratory, Public Health Agency of Canada, Winnipeg, Manitoba, Canada
Received 19 December 2005/ Returned for modification 3 March 2006/ Accepted 7 April 2006
|
|
|---|
|
|
|---|
Both serovars Typhi and Paratyphi A are highly adapted and can cause only systemic disease in humans. Other subspecies I S. enterica serovars can result in nonsystemic disease symptoms in humans and may be commensal or opportunistic pathogens in other warm-blooded animals. Sequencing of the serovar Typhi (strains CT18 and Ty2) and serovar Paratyphi A (ATCC 9150) genomes has revealed significant similarities between the two pathogens (8, 25, 32). Comparative genomic hybridization (CGH) experiments and phylogenetic analyses with Salmonella serovars have also demonstrated the genetic relatedness of serovars Typhi and Paratyphi A, even though they are members of different serogroups (serogroups D1 and A, respectively) (4, 34). Genome degradation has led both pathogens to independently accumulate a high proportion of pseudogenes, most of which are required by other Salmonella serovars to colonize and invade the gastrointestinal epithelium (25). Gene silencing from pseudogene formation, along with other loss-of-function mutations, has resulted in the adaptation of serovars Typhi and Paratyphi A to a human-specific niche in the last few thousand years (3, 19, 25, 32). The genetic determinants responsible for their systemic nature during human infection have not been readily identified. Serovars Paratyphi B and C are antigenically and genetically distinct from serovars Typhi and Paratyphi A (34, 40), and while all Paratyphi serotypes can result in typhoid-like illness in humans, serovars B and C can also result in zoonotic infections.
Diagnosis of enteric fever has traditionally been based on blood culture or the Widal test for serum antibodies against O-somatic and H-flagellar antigens, although the latter suffers from substantial variations in interlaboratory specificity and sensitivity (33). Molecular immunology kits that have improved sensitivity and specificity over those of the Widal test are available (31). Once pure cultures have been isolated from clinical samples, serovars Paratyphi A and Typhi are readily identified by Kaufmann-White serotyping (antigenic formulas, 1,2,12:a:[1, 5] and 9,12,[Vi]:d:, respectively) and by biochemical tests based on the differential production of H2S and decarboxylation of lysine, with serovar Paratyphi A strains typically being negative for both of these characteristics on lysine-iron agar. PCR-based techniques have been developed to discriminate serovar Typhi from other Salmonella serovars; these techniques target the Vi antigen-encoding gene and the flagellin antigen fliC-d gene (9, 13, 14, 23, 42). Individually, these targets were not entirely specific for serovar Typhi, including some Vi-negative strains of serovar Typhi that are endemic (1), and neither of these loci are encoded by serovar Paratyphi A. The use of these loci along with O- and H-antigen encoding genes in a five-locus multiplex PCR assay discriminated both serovars Typhi and Paratyphi A from a large panel of S. enterica serovars (14). Real-time PCR strategies, however, are a significant improvement over standard PCR methods, leading to more rapid, sensitive, and potentially quantitative results. Recently, real-time PCR methods have been successful in identifying the Vi-antigen gene (9) and estimating serovar Typhi bacterial loads in blood samples (24). Our goal was to identify genetic traits that are characteristic for both serovar Typhi and serovar Paratyphi A by using comparative sequence analyses and to develop real-time PCR and liquid microsphere suspension assays for the molecular identification of these human-adapted typhoidal Salmonella strains.
|
|
|---|
|
View this table: [in a new window] |
TABLE 1. Bacterial strains used in this study
|
PCR and sequencing. Template DNA was prepared by centrifuging 1 ml of log-phase cultures grown in brain heart infusion broth, resuspending the pellet in 1 ml of TE buffer (10 mM Tris-HCl, 1 mM EDTA, pH 8.0) (Sigma, St. Louis, MO), and boiling for 10 min. The boiled cell debris was pelleted, and the supernatant was removed and used as the DNA template in real-time and standard PCRs.
Standard PCR was performed with Platinum High Fidelity Taq (Invitrogen, Burlington, Ontario, Canada) following the manufacturer's directions. The oligonucleotides are described in Table 2. PCR conditions were as follows: initial denaturation at 94°C for 5 min and 30 cycles of denaturation at 94°C for 30 s, annealing at 50°C for 30 s, and extension at 68°C for 30 s, with a final extension at 68°C for 5 min. The PCR products were purified by using a QIAquick PCR purification kit (QIAGEN, Mississauga, Ontario, Canada) and were sequenced by using the same primers used to generate this template. Sequencing was performed on an ABI 3730 instrument (Applied Biosystems, Foster City, CA).
|
View this table: [in a new window] |
TABLE 2. Oligonucleotides used in this study
|
For 5'-nuclease real-time PCR, primers and probes (Table 2) were designed with Applied Biosystems Primer Express software, version 2.0. TaqMan Universal PCR Master Mix and No AmpErase UNG (Applied Biosystems) were used as the amplification mixture, with the final concentrations of the 5'-nuclease probes and primers (Operon Biotechnologies Inc., Huntsville, AL) being 125 nM and 800 nM, respectively. The total reaction volume was 25 µl, including 2.5 µl of template. The DNA was amplified in a SmartCycler 2.0 instrument by an initial denaturation at 95°C for 10 min and 40 cycles of denaturation at 95°C for 15 s and a single annealing-extension step at 60°C for 60 s. Fluorescence was detected during the annealing-extension step, and the threshold level was set at 30 fluorescence units.
Liquid microsphere suspension arrays. Allelic discrimination of trpS and sseC was achieved after PCR amplification of biotin-labeled target DNA from representative strains with a GeneAmp 9700 thermocycler (Applied Biosystems), primers GIL259 and GIL260-L (trpS) or primers GIL286 and GIL287-L (sseC) (Table 2), and the thermocycling parameters detailed above, except that 35 cycles were performed. The PCR mixtures were purified with QIAquick DNA purification kits (QIAGEN) and eluted with 50 µl of EB buffer (QIAGEN). Oligonucleotides GIL260-L and GIL287-L contain four bases with phosphorothioate linkages, as well as a biotin molecule, all at the 5' end. Between the two strands of the target DNA, the strands produced from GIL259 and GIL286 (trpS and sseC "sense" strands) are sensitive to T7 exonuclease digestion, whereas the antisense strands are protected due to the phosphorothioate linkages (12, 29). DNA digestion was performed by mixing 43 µl of purified PCR product with 5 µl of buffer 4 and 2 µl of T7 exonuclease (both from New England Biolabs, Ipswich, MA) (20 U total) and incubating at 37°C for 1 h. T7 exonuclease was inactivated by adding 2 µl of 0.5 M EDTA (Ambion, Austin, TX). Selective degradation ensures elimination of the unlabeled target DNA strand, thereby preventing reannealing between the two target DNA strands during hybridization, which, if left intact, would limit the intended hybridization to that between the biotin-labeled strand and the trpS or sseC allele-specific probes coupled to microspheres in subsequent steps.
Oligonucleotide probes were designed by matching the sense strand in regions characteristic for individual allele subtypes. Oligonucleotides were screened for potential secondary structures or cross-hybridization between probes by using SBEprimer software (16). The oligonucleotide probes were synthesized with a 5' C-12 amine and coupled to xMAP-carboxylated fluorescently coded microspheres (Luminex Corporation, Austin, TX). Microsphere sets 103 and 108 were coupled to oligonucleotides DOB75 and DOB78, respectively (Table 2); and hybridization of biotin-labeled trpS and sseC target DNA strands to the capture probe-coupled microspheres and flow cytometry were performed in triplicate, as described previously (12). The positive cutoff value for the sseC assay was chosen to be a value 10 times greater than the value for the negative control; for the trpS assay the positive cutoff value was chosen to be a value 2.5 times greater than the value for the negative control, accounting for the lower overall signal strength of the true-positive signals.
Bioinformatics. Initial screening of polymorphic loci was performed with BLASTn, and position-specific iterated basic local alignments of the STY4217 to STY4222 gene products were performed with PSI-BLAST (www.ncbi.nlm.nih.gov/BLAST/). Pairwise global DNA sequence alignments were performed with Align (http://www.ebi.ac.uk/emboss/align/), multiple-sequence alignments were completed with ClustalW (www.ebi.ac.uk/clustalw/) and Boxshade (www.ch.embnet.org), neighbor-joining trees were constructed with the Hasegawa-Kishono-Yano (HKY85) distance correction with SplitsTree4 (15), and genetic diversity statistics were calculated with DnaSP 4.10.3 (37). Split decomposition analysis was performed with SplitsTree4 by using the alignment inputs created by ClustalW, and the calculations used only parsimony-informative sites. Both neighbor-joining and split decomposition trees were calculated by using only those segments of each locus for which data for all strains examined were available (i.e., regions from the complete genome data outside of the amplified segments of each locus were not included). Artemis (38) was used for calculation of the G+C content and visualization of the annotated features of the Salmonella chromosomal segments.
Nucleotide sequence accession numbers. The sequence data from this study were deposited in GenBank under accession numbers DQ320510 to DQ320558 and DQ451529 to DQ451530.
|
|
|---|
![]() View larger version (54K): [in a new window] |
FIG. 1. Comparative genomic analyses of the S. enterica core locus trpS, which encodes tryptophanyl-tRNA synthetase. (A) Phylogeny of clinically significant S. enterica serovars based upon a neighbor-joining tree; the scale of the distance score is represented by a horizontal bar. (B) Split decomposition of trpS; recombination between loci is indicated when the topology resembles a network, rather than the single branching points seen in normal tree topologies; the scale of the distance score is represented by a horizontal bar. (C) Multiple-sequence alignment of the trpS segment containing potential serovar-specific polymorphisms; strain identification is indicated in parentheses, and where sequence data from multiple strains from a single serovar were identical, no strain identifier was included; the 5'-nuclease probe is identified with a horizontal line under the corresponding sequence, and the microsphere-coupled probe is identified with a horizontal line labeled with a circle. The GenBank accession numbers for previously sequenced loci are presented in the Materials and Methods.
|
|
View this table: [in a new window] |
TABLE 3. Genetic diversity of trpS, sseC, and sseF
|
![]() View larger version (88K): [in a new window] |
FIG. 2. Comparative genomic analyses of the S. enterica pathogenicity island 2 locus sseC. (A) Phylogeny and (B) split decomposition of sseC; see the Fig. 1 legend for further details. (C) Multiple-sequence alignment of the sseC segment containing potential serovar-specific polymorphisms; strain identifications are indicated in parentheses, and where sequence data from multiple strains from a single serovar were identical, no strain identifier was included; LUX real-time PCR primers are identified by a horizontal lines with half arrowheads under the corresponding sequence, and the microsphere-coupled probe is identified with a horizontal line labeled with a circle. The GenBank accession numbers for previously sequenced loci are presented in the Materials and Methods.
|
![]() View larger version (92K): [in a new window] |
FIG. 3. Comparative genomic analyses of the S. enterica pathogenicity island 2 locus sseF. (A) Phylogeny and (B) split decomposition of sseF; see the Fig. 1 legend for further details. (C) Multiple-sequence alignment of the sseF segment containing potential serovar-specific polymorphisms; strain identifications are indicated in parentheses, and where sequence data from multiple strains from a single serovar were identical, no strain identifier was included; PCR primers are identified by horizontal lines with half arrowheads under the corresponding sequence. The GenBank accession numbers for previously sequenced loci are presented in the Materials and Methods.
|
PCR-based detection of serovars Typhi and Paratyphi A. To confirm the presence of the STY4217 to STY4222 region in additional serovar Typhi and Paratyphi A strains and to examine additional serovars, we developed LUX real-time PCR primers for STY4221 (Table 2; Fig. 4). The LUX system requires only one fluorogenic primer, which is self-quenching through the formation of a hairpin loop, and one unlabeled primer (28). As a more economical alternative, standard PCR primer pairs were also developed for the detection of STY4220 (186-bp product) and STY4221 (264-bp product), and our panel of strains were examined for the presence of these genes (Table 4). Both STY4220 and STY4221 were detected only in serovars Typhi, Paratyphi A, and Sendai; and all strains of these particular serovars were positive. Notably, S. enterica serovar Sendai is also a human-adapted serovar that can result in systemic, typhoid-like disease symptoms (40), although it is rarely observed in clinical laboratories. Within the previously published CGH data, these genes were observed to be absent in serovar Sendai reference strain SARB58 (34); however, the accuracy of a PCR assay directed against a single locus is possibly higher than that of an individual probe in a CGH experiment comprising a whole genome. A total of seven serovar Sendai strains (originating from Canada, the United States, and China) were screened, and all strains were positive by the STY4220 and STY4221 PCR assays. Sequence analysis of a STY4221 amplicon overlapping the LUX primer-binding sites identified a single conserved polymorphism in the serovar Sendai and Paratyphi A strains compared to the sequences of the serovar Typhi strains, but this position was outside of the primer-binding sites (data not shown).
![]() View larger version (23K): [in a new window] |
FIG. 4. Allelic discrimination of serovar Typhi and typhoidal serovar-specific polymorphisms and detection of loci unique to human-adapted typhoidal serovars by real-time PCR. (A) 5'-Nuclease real-time PCR characterization of trpS alleles; (B) LUX real-time PCR detection of STY4221; (C) LUX real-time PCR characterization of sseC alleles. The strains examined are indicated in the legend insets. NTC, no-template control. Positive experimental reactions and the no-template control reaction in each data set have solid lines connecting the datum points, whereas the negative experimental reactions are indicated with dashed lines. In the sseC LUX assay, positive experimental reactions where the product was detected more than seven cycles later than serovar Typhi was detected are indicated with shaded markers at each datum point. The primers are described in Table 2.
|
|
View this table: [in a new window] |
TABLE 4. Summary of allelic discrimination for trpS, sseC, and sseF by PCR, real-time PCR, and liquid suspension microsphere arrays
|
To determine the success of these primer and probe designs, real-time or standard PCRs were performed with template DNA from serovars Typhi and Paratyphi A. An expanded panel of serovars with known pathogenic or genetic similarities to serovars Typhi and Paratyphi A, as well as strains representing the most commonly observed S. enterica serovars in Canadian clinical laboratories, was also examined (Table 4; Fig. 4). The trpS probe was successful in differentiating serovar Typhi from all other serovars examined except for a single Paratyphi B isolate, S1583. The LUX-sseC reaction was positive with all serovar Typhi, Paratyphi A, and Sendai strains. A LUX-sseC product was also generated with all serovar Agona, Oranienburg, Reading, Javiana, and Paratyphi C strains examined and a single serovar Paratyphi B isolate, S1583 (Fig. 4; Table 4); however, in all instances the product formed more than seven cycles after the serovar Typhi product formed. The sseF PCR was positive with serovars Typhi, Paratyphi A, Sendai, Infantis, and Paratyphi B isolate S1583 (Table 4).
Liquid microsphere suspension arrays for trpS and sseC. To investigate supplemental or improved allelic discrimination methods, microsphere-oligonucleotide probe conjugates targeting serovar Typhi-specific alleles (trpS) and human-adapted typhoidal serovars (serovars Typhi, Paratyphi A, and Sendai; sseC alleles) were developed (Fig. 1, 2, and 5). The Luminex microsphere suspension array technology uses allele-specific oligonucleotide probes conjugated to fluorescently coded microspheres to capture soluble DNA in a liquid phase and characterize hybridization partners with flow cytometry. Target DNA for the trpS and sseC loci was generated by using biotinylated PCR primers with internal phosphorothioate linkages (Table 2) (see Materials and Methods), and the region of each locus amplified for liquid microsphere suspension assays corresponded to the same region subjected to sequence analysis.
![]() View larger version (14K): [in a new window] |
FIG. 5. Allelic discrimination of serovar Typhi and human-adapted typhoidal serovar-specific polymorphisms by liquid microsphere suspension arrays. Biotin-labeled template (trpS [A] and sseC [B]) was amplified from strains of the indicated serotypes and incubated with fluorescently coded microspheres coupled with an oligonucleotide probe targeting the serovar-specific alleles. The background fluorescence contributed by the microsphere-probe mixture was determined by using a no-template control (TE buffer), and the standard errors are indicated by vertical lines on each bar.
|
Comparative analyses of the trpS, sseC, sseF, and STY4217 to STY4222 sequences. To confirm the genotypes of all serovars producing positive reactions (serovars Typhi, Paratyphi A, Sendai, Agona, Javiana, Paratyphi C, Infantis, Oranienburg, Reading, and Paratyphi B strain S1583), the corresponding regions of trpS, sseC, and sseF were sequenced (Fig. 1 to 3). We also sequenced these loci of additional strains of serovars Typhimurium, Muenster, Heidelberg, and Dublin to provide insight into the relationship between the genotype and the PCR results, as well the phylogenetic associations among S. enterica serovars. Primers to amplify and sequence each locus were designed for conserved sites identified within multiple-sequence alignments of the available S. enterica sequence data (Table 2). The selected region of each allele comprised the putative serovar-specific polymorphisms and other variable sites; and the amplicon lengths for trpS, sseC, and sseF were 334 bp, 295 bp, and 278 bp, respectively. These data indicated that serovars producing positive PCR or microsphere array reactions encoded sites identical or similar to the targeted alleles of trpS, sseC, and sseF (Fig. 1 to 3).
The genetic diversity and number of synonymous and nonsynonymous mutations were calculated for the trpS, sseC, and sseF loci by using the sequence data from the amplicons produced in this study and the publicly available complete coding sequence (CDS) data from reference strains (Table 3). At the trpS, sseC, and sseF loci, the sequenced amplicon represented 17%, 19%, and 31% of the complete CDSs, respectively; but these segments accounted for 31%, 49%, and 60% of the total polymorphic sites encoded in the complete CDSs, respectively (Table 3). Notably, the cellular function of each gene product was reflected in the ratio of the number of nonsynonymous mutations to the number of synonymous mutations (dN/dS). At the trpS locus, which encodes a core bacterial determinant essential for the translation module, a single nonsynonymous site was observed, whereas the sseC and sseF loci, both of which encode secreted virulence proteins, had more nonsynonymous sites than synonymous sites (Table 3). The higher dN/dS likely reflects the adaptive nature of these gene products to host structures during pathogenesis. Notably, none of the targeted alleles of trpS (serovar Typhi) or sseC and sseF (serovars Typhi, Paratyphi A, and Sendai) have been subjected to recombination, as suggested by split decomposition analysis (Fig. 1 to 3).
The island from STY4217 to STY4222 was inserted between yhiI and yhiN genes (present as adjacent CDSs in S. enterica serovar Typhimurium strain LT2, STM3587 to STM3588), and this genetic layout is similar to that observed in serovar Typhi strain Ty2 (CDSs t3930 to t3934) and serovar Paratyphi A strain ATCC 9150 (CDSs SPA3439 to SPA3443). Bioinformatic analyses of this region did not identify any obvious virulence factors or indicate the possible origin. Each of the individual coding sequences from STY4217 to STY4222 had other possible orthologues encoded by enteric proteobacteria, as detected by PSI-BLAST analysis (data not shown); but no currently sequenced genome encoded contiguous coding sequences that had a gene order similar to that of this operon or a protein sequence identity with this operon (other than the serovar Typhi strain Ty2 and the seorvar Paratyphi A genomes). Because of the the conservation and exclusivity of the island from STY4217 to STY4222 to the human-adapted typhoidal Salmonella strains, the island provides tremendous promise for use in the identification of systemic serovars by molecular tests that target this region.
|
|
|---|
Serovar Typhi is one the most genetically distinct but homogeneous serovars of S. enterica, as observed by multilocus enzyme electrophoresis (36, 40), multilocus sequence typing (19), protein profiling (11), and plasmid analysis (22). Comparative genomic analyses with DNA microarrays have also identified that serovars Typhi, Paratyphi A, and Sendai have similar genetic contents but are diverse from the other S. enterica serovars (4, 34). Serovars Sendai and Paratyphi A also share many serological characteristics, with both serovars having similar O antigens and identical phase 1 and 2 flagellar antigens (40). The genetic relatedness of serovars Typhi, Paratyphi A, and Sendai was represented in the characteristic polymorphisms of the sseC and sseF alleles, as well as the presence of the entire island from STY4217 to STY4222. The sseC and sseF loci are among the most polymorphic in pathogenicity island 2; but the alleles of each locus encoded by serovars Typhi, Paratyphi A, and Sendai were 99.6 to 100% identical. These alleles were each distinct from those encoded by other serovars, paralleling the overall genetic relatedness of S. enterica serovars; and the prevalence of nonsynonymous mutations in sseC and sseF might be reflective of host adaptation. The cellular function contributed by the island from STY4217 to STY4222 was not obvious after comparison to other annotated bacterial genomes; but coinheritance of this region suggests that serovars Typhi, Paratyphi A, and Sendai have a common ancestor. The two polymorphic sites observed in trpS and considered unique to serovar Typhi are unlikely to contribute to the pathogenesis of this serovar, but these sites signified divergence between serovars Typhi and Paratyphi A outside of the pseudogenes. Notably, recombination between serovars was indicated at each locus examined by split decomposition analysis but was not indicated at the serovar Typhi-encoded alleles of trpS nor at the serovar Typhi-, Paratyphi A-, and Sendai-encoded alleles of sseC and sseF. A lack of recombination might indicate functional constraint on further mutation of these alleles. Comparative genomic analyses of a group of serovars with related pathogenicity traits (such as the human-adapted typhoidal Salmonella) therefore identified lineage-specific polymorphisms and coding sequences, and these determinants might be responsible for the pathogenicity traits (or are minimally coinherited with the determinants that are responsible).
The PCR-based allelic discrimination assays for trpS, sseC, and sseF detected all intended serovars (serovar Typhi and/or Paratyphi A), but at each locus other strains or serovars had positive reactions. This included a single serovar Paratyphi B strain (S1583) with positive reactions for all three loci; serovars Paratyphi C, Agona, Oranienburg, Reading, and Javiana with positive reactions for sseC; and serovar Infantis with a positive reaction for sseF due to identical or nearly identical sequences at the real-time primer and probe binding sites (Fig. 1 to 3). Notably, the subpopulation of serovars that produced sseC at later cycle thresholds than serovar Typhi (more than seven cycles) each encoded polymorphisms at the 5' ends of the primer-binding sites (in relation to the serovar Typhi sequences). This indicates that LUX primers can tolerate mismatches outside the 3' regions but that product is not formed at the same rate as template without mismatches. The specificity at sseC was improved by targeting an alternate site in the microsphere suspension assay; and by that assay only serovars Typhi, Paratyphi A, Paratyphi B, and Sendai were detected. Serovar Paratyphi B isolate S1583 was d-tartarate positive (indicative of biovar Java). This serovar normally causes nonsystemic disease symptoms (35), so it was surprising that there was considerable similarity to the serovar Typhi-encoded alleles of trpS, sseC, and sseF. Notably, there are extensive genetic differences within strains of serovar Paratyphi B (27, 35, 39), so the observation of a single strain with alleles similar to those of the human-adapted serovars might be possible. Otherwise, use of the island from STY4217 to STY4222 for molecular identification of human-adapted typhoidal serovars resolved the cross-reaction with strain S1583, as this region was not encoded in any serovar Paratyphi B isolate. Additionally, since this region was encoded exclusively in all strains of serovars Typhi, Paratyphi A, and Sendai examined, the STY4220 and STY4221 CDSs were ideal markers for the detection of human-adapted typhoidal Salmonella, whereas serovars Paratyphi A and Sendai would not be detected by use of the current fliC-d- or Vi antigen-specific reagents.
The endemic nature of serovar Typhi throughout a significant proportion of the world, the emergence of antibiotic resistance, and the rising incidence of serovar Paratyphi A in Asia demand the development of sensitive molecular protocols and vaccines suitable for the identification and prevention of typhoidal Salmonella infections other than just those caused by serovar Typhi (2, 30, 47). Although real-time PCR offers a rapid method for the detection of serovars Typhi and Paratyphi A and this platform can readily be deployed in outbreak situations, the expense and availability of the equipment required remain limiting factors in its widespread use in routine surveillance worldwide. Accordingly, standard PCR assays with nonfluorogenic oligonucleotides were also developed in this study. These molecular methods were tested only with pure cultures of S. enterica; and it will be necessary to test these reagents with human clinical samples such as blood, urine, and stool to determine if typhoidal Salmonella can be detected without the initial requirement of microbial culturing or examination of circulating antibodies. If these protocols can successfully identify serovars Typhi, Paratyphi A, and Sendai (especially in blood, which is one of the first sites from which these organisms can be isolated from during illness), then incidences of enteric fever caused by S. enterica could be more accurately diagnosed and the worldwide surveillance of these pathogens would be greatly enhanced.
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»