Julie T. Collins,
Timothy C. Barnett, and
June R. Scott*
Department of Microbiology and Immunology, Emory University School of Medicine, Atlanta, Georgia 30322
Received 30 June 2005/ Accepted 6 October 2005
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
For GAS to cause disease in such diverse host tissues, this organism must be able to withstand a variety of environmental factors and grow in different microenvironments. GAS can face many stress conditions within the human host, such as the increased salt concentration on the skin, the low pH of the vagina, and the pH fluctuations within the nasopharynx. Also, the bacteria must be able to overcome the stress conditions that arise as a consequence of the GAS infection itself, such as high temperature resulting from fever and decreased pH as a result of abscess formation.
Most bacteria use two-component regulatory systems to sense the various environments they encounter and to alter their gene expression in response to these signals (25, 42). In general, these systems consist of a membrane-bound histidine kinase and a cytoplasmic response regulator. In response to a specific environmental signal, the histidine kinase autophosphorylates and transfers the newly acquired phosphoryl group to an aspartate residue on the response regulator protein. The phosphorylated response regulator undergoes a conformational rearrangement that allows it to bind to the promoter region of target genes, resulting in a change in the level of transcription. Many histidine kinases also have phosphatase activity that catalyzes the reverse reaction; i.e., the phosphate is removed from the response regulator, causing it to release the DNA. The CovRS two-component system has been shown by microarray analysis to regulate
15% of the M1 GAS genome, including genes encoding several virulence factors (21). Unlike most response regulators that activate gene expression, CovR functions predominately as a repressor of transcription (16, 21). We have shown previously that this system plays a role in the GAS stress response (10). Inactivation of CovR by its cognate sensor kinase, CovS, is required for growth of GAS under heat, acid, and salt stress. A CovS mutant, in either M6 or M3, does not grow under these stress conditions unless a second site mutation occurs that inactivates CovR, suggesting that CovR represses genes whose products are needed for growth under stress (10). The goal of this work was to identify the CovR-repressed gene or genes responsible for this phenotype.
The best-studied stress response genes in GAS and other gram-positive organisms are those required for growth at elevated temperatures. Heat shock proteins primarily function either as chaperones that promote the refolding of denatured proteins or as proteases that degrade misfolded proteins. Heat shock genes are commonly divided into five classes. The class I heat shock genes, which include the chaperones GroEL and DnaK, are under the control of the repressor HrcA. The class II heat shock genes are transcriptionally regulated by an alternative sigma factor,
B. There are no
B homologs encoded in the GAS genome sequences, so this class of heat shock gene either is not present in GAS or is under different regulatory control (4, 18, 45). Class III heat shock genes, which include the clp ATPases, are repressed by CtsR, which is found in most low-G+C gram-positive organisms (8, 50). Class IV genes are those that are not under the control of HrcA, CtsR, or
B (52). Class V genes are heat stress genes that are under the control of the two-component regulatory system CssRS (11).
Many heat shock genes are also induced in response to acid in several gram-positive bacteria. For example, in addition to being induced by high temperature, acid shock also induces groEL and dnaK in Streptococcus mutans, while in Lactococcus lactis, the levels of both class I- and class III-regulated gene products are increased at low pH (19, 26, 31). In addition to genes regulated by both heat and acid, some genes are specifically upregulated only in acidic conditions. For example, in many species of streptococci, the activities of the F1Fo ATPase and arginine deiminase are induced in response to acid (1, 46, 50). The F1Fo ATPase seems to be critical for the protection of bacteria in acidic environments by exporting excess protons from the cell (40). Arginine deiminase activity leads to the generation of ammonia under acidic conditions to help maintain pH homeostasis (22).
In this work, we have used microarray analysis of transcription to gain insight into the heat and acid stress regulons of GAS. Because we used a different strain of GAS and different growth conditions from those previously used by Graham et al. (21), we also performed microarray analysis to identify the CovR regulon of our M type 6 strain of GAS. The global gene expression data we obtained were used to identify CovR-repressed genes that are activated in response to heat and acid. We demonstrate here that one such gene, which we name rscA (regulated by stress and Cov), is required for optimal growth of GAS at 40°C.
| MATERIALS AND METHODS |
|---|
|
|
|---|
Microarray construction. Our microarrays consist of 2,328 unique 70-mer oligonucleotide probes designed and synthesized by Operon Technologies (Valencia, CA) using a proprietary algorithm to target unique nonrepetitive open reading frames (TIGR annotation as of 2 December 2002) present in the sequenced genomes of a strain of serotype M1 (SF370; 1,932 probes), M3 (MGAS315; 165 probes), and M18 (MGAS8232; 231 probes). Oligonucleotides were spotted in triplicate onto UltraGAPS2 microarray slides (catalog number 40017; Corning) by Microarrays Inc. (Nashville, TN). Slides were rehydrated, cross-linked, and blocked in a prehybridization solution (50% formamide, 5x SSC [1x SSC is 0.15 M NaCl plus 0.015 M sodium citrate], 0.1% sodium dodecyl sulfate [SDS], 0.1 mg/ml bovine serum albumin) according to the manufacturer's instructions (Corning).
cDNA labeling and hybridization conditions. Total RNA was isolated as described previously (35). Labeled cDNA was generated from RNA samples by incorporation of aminoalloyl-dUTP (Fairplay microarray labeling kit, catalog number 252003; Stratagene) using 3 µg random hexamers (catalog number 48190-011; Invitrogen) in each reaction mixture. The cDNA was then directly labeled with Cy3 or Cy5 monoreactive dyes (catalog numbers PA23001 and PA25001; Amersham) according to the manufacturer's instructions. Labeled cDNA was purified using QiaQuick PCR purification columns (catalog number 28104; QIAGEN) and concentrated using Microcon YM30 columns (catalog number 42410; Amicon). The amount of incorporated label was determined by measuring the absorbance at 550 nm (Cy3) or 650 nm (Cy5).
For each set of microarray experiments, a reference cDNA sample consisted of pooled Cy5-labeled cDNA generated from all RNA samples to be examined. The test samples consisted of cDNA generated from RNA from each condition that was labeled with Cy3 in triplicate and pooled before hybridization. Therefore, the 37°C arrays for the heat shock experiments consisted of Cy5-labeled reference cDNA (pooled cDNA samples from 37°C, 40°C/5 min, 40°C/10 min, and 40°C/20 min) and Cy3-labeled test cDNA (37°C cDNA labeled in triplicate and pooled). For each array, an aliquot of reference cDNA and test cDNA containing 60 pmol each Cy dye was combined in 60 µl hybridization buffer 1 (catalog number 8861; Ambion), heated at 95°C for 5 min, and immediately centrifuged at 12,000 x g for 2 min at room temperature before being applied to the array. The probe was allowed to hybridize to the array at 42°C for at least 16 h. Following hybridization, the arrays were washed once with 2x SSC-0.1% SDS, once with 0.1x SSC-0.1% SDS, and twice with 0.1x SSC.
Microarray data analysis. The washed arrays were scanned on a GenePix 4000A scanner (Axon instruments). The resulting images were inspected to remove values from bad spots (for example, those affected by scratches) and analyzed using GenePix Pro 3.0 software. For each of the hybridizations, signals below a threshold level (the mean background level of the whole experiment plus 1 standard deviation) were set to the threshold level to avoid spurious expression ratios. The data for each gene from each of three individual spots on the microarray were averaged and treated as a single value for the experiment. The data were then filtered to remove values for genes that were not present in strain JRS4 (the strain examined) as identified independently by hybridization of Cy3-labeled genomic DNA to the array.
The expression value for each gene was then expressed as a fraction of the total amount of expression on each array (global normalization) in order to make comparisons between arrays valid. Normalized data were then log2 transformed, and statistical analysis was performed using Significance Analysis of Microarrays software (49).
Microarray data have been deposited in the Gene Expression Omnibus database (http://www.ncbi.nlm.nih.gov/geo) (accession no. GSE 3367).
Probe construction and RPAs.
Probes (
250 bp in length) for regions within the coding sequence of indicated genes were constructed as described previously (10). The sequences of primer pairs used to amplify specific probes from the JRS4 chromosome will be provided upon request. Probes were labeled with [
-32P]UTP as previously described (10). T7 polymerase was used in the in vitro transcription reactions for the ska, groEL, ctsR, and atpB probes. SP6 polymerase was used in the in vitro transcription reactions for the clpC, purC, and rscA probes. RNase protection assays (RPAs) were performed as described previously (10). Transcript levels were measured by densitometry using ImageQuant software.
Insertional inactivation of rscA. An internal region of rscA was PCR amplified from the JRS4 chromosome using primers YdaG-S1 (gccgcgcgaattcTCATGGGTTTGGTGATGCGTCAA) and YdaG-A1 (gccgcgcgaattcTTGGCTCATCATCATTAGGATAT) (lowercase letters indicate bases added; underlined letters are introduced restriction sites). The PCR product was then digested with EcoRI and ligated into EcoRI-digested pSK-Erm, a vector that cannot replicate in GAS (6). The resulting plasmid, pEU7248, was then transformed into electrocompetent JRS4 cells, and erythromycin-resistant transformants were selected. The resulting strain is JRS1819. The insertional inactivation of rscA was confirmed by PCR amplification of the junctions between the chromosome and the plasmid using the following primer pairs: 5' junction, ggccgggatccATATGCTAACTTATAAAGG (ydaG-F) and GTAAAACGACGGCCAG (M13F); 3' junction, ggccgctgcagTTTTCATCACTACTCCTC (ydaG-R) and CAGGAAACAGCTATGA (M13R). The same primer sets were used to confirm the rscA insertion mutation following growth of JRS1819 at 40°C.
Construction of the rscA-complementing plasmid pJRS334. To clone rscA and rscB under the constitutively expressed lactococcal promoter P23 (39), a PCR fragment containing the coding regions of rscA and rscB and their ribosome-binding sites was amplified from JRS4 using primers ydaG-F-EcoRI (ggccgggaattcATATGCTAACTTATAAAGG) and ydaG-R-BamHI (ggccggggatccTTATTCAAAGACAAATTG). To construct pJRS334, the 3.5-kb PCR product was digested with EcoRI and BamHI and ligated into EcoRI/BamHI-digested pJRS7759. To construct pJRS7759, P23 was digested out of pORI23 (39) with PstI/EcoRI and ligated into PstI/EcoRI-digested pLZ12 (A. Benitez and S. Chancey, unpublished data). pLZ12 is a chloramphenicol-resistant shuttle vector that is able to replicate in GAS (13).
Temperature shift growth curves. Strains were grown at 37°C until they reached the early exponential phase, at which time the cultures were diluted 1:20 in THY that had been prewarmed to 40°C. Growth was monitored with a colorimeter (Klett-Summerson).
| RESULTS |
|---|
|
|
|---|
(ii) The CovR regulon of strain JRS4. To identify stress response genes that are repressed by CovR, we first used microarrays to determine the global transcriptional profiles of a wild-type M6 strain (JRS4) and an isogenic strain with a deletion of covR (JRS948) (16). These strains were grown at 37°C in THY broth until late in the exponential growth phase, when covR expression is maximal (16), and total RNA was harvested. To minimize variation that may occur due to differences in culture conditions and RNA isolation procedures, we harvested RNA individually from four independent cultures of both JRS4 and JRS948. All RNA samples were used to make cDNA that was then labeled with the appropriate monoreactive dye. Independent cDNA preparations (test) were hybridized individually to at least two array slides in combination with the reference RNA standard (total pooled RNA). Data were then normalized as described in Materials and Methods and analyzed using Significance Analysis of Microarrays software (49). The genes regulated by CovR in JRS4 are listed in the supplemental material (see Tables S1 and S2 in the supplemental material).
Previously identified, well-studied CovR-repressed genes and operons were shown to be upregulated in the absence of CovR in our arrays, thus validating our results. These genes include ska, sda, and grab as well as the has, sag, and covR operons. However, there are significant differences in the genes identified as being CovR regulated between our results and those previously reported by Graham et al. for an M1 strain (21). Possible explanations for the differences between these two sets of data include differences in growth conditions and the time at which the RNA was harvested (Graham et al. grew their strains in 5% CO2, while we grew ours in ambient air), variations in stringency and/or standardization used for analysis, and strain-specific differences. For example, in our arrays, we found that many phage-associated genes, which are not present in the M1 strain, are activated (either directly or indirectly) by CovR.
(iii) Global effects of high temperature and low pH on gene expression in GAS. The above-described microarray results indicated that, although expressed at a detectable level, the well-studied stress response genes such as groEL, clpC, and dnaK showed little difference in transcript level in the CovR strain compared to that in the CovR+ strain (1.1-, 1.2-, and 0.7-fold differences, respectively). Consequently, we felt that it was unlikely that these genes are responsible for the growth defect of covS mutant strains under stress conditions. Therefore, we looked for new genes that were differentially expressed in a wild-type strain following either an increase in temperature or a decrease in pH. For the heat response experiment, wild-type strain JRS4 was grown to early logarithmic phase at 37°C in THY broth, and RNA was isolated from a portion of the culture. The remaining culture was centrifuged at 37°C, and the cell pellet was resuspended in prewarmed THY broth and incubated at 40°C. RNA was isolated from the cells 5, 10, and 20 min after the temperature shift.
Both heat-activated and heat-repressed genes were identified, and they included those genes predicted to encode proteins involved in such functions as fatty acid biosynthesis, metabolism, regulation, transport, and protein synthesis (data not shown). Because RNA was harvested from only one culture at each time point, these data should not be considered conclusive. However, the array analysis indicated that the well-characterized heat shock genes groEL, ctsR, clpC, and purC were activated at 40°C, and we confirmed this activation in response to heat using RPAs (Fig. 1A to D). The transcript level of groEL increased approximately threefold at 5, 10, and 20 min following the shift to 40°C. The RPAs also showed that ctsR was activated
5-fold after 5 min,
3-fold after 10 min, and
2-fold after 20 min at 40°C compared to that at 37°C. Likewise, transcript levels of clpC increased
3-fold at 5 and 10 min, and these levels fell back to
2-fold above the 37°C value 20 min following the temperature shift. Under these same conditions, there was a
2-fold increase in purC transcript at all 40°C time points. As a control, we also assayed the ska transcript levels by RPA. As expected, we found by both microarray analysis and RPA that ska was not regulated by heat (Fig. 1E).
|
Acid-activated genes included those that encode the F1Fo ATPase, which is a major component of the acid shock response in many organisms (50). The microarray results for one member of this ATPase operon, atpB, were confirmed by RPA using newly isolated RNA (Fig. 2A). Quantitation of the transcript levels indicated in the RPA showed a
3-fold increase in atpB transcription levels at 5, 10, and 20 min following the shift to pH 6.0. ska was used again as a control, and expression was unaffected by shifting the culture to low pH as determined both by microarray analysis and RPA (Fig. 2B).
|
|
Since a CovS mutant is unable to grow under either acid or heat stress in the presence of CovR, we reasoned that genes that are activated under both stress conditions might be involved in a response that is repressed by CovR in the absence of stress (10). Thus, we further investigated the members of class 1, the CovR-repressed genes that are activated by both heat and acid stress (Table 1). We identified four genes in this class: SPy0229, SPy0230, SPy0570, and SPy1404. SPy1404 encodes a protein of unknown function with no known homologs. SPy0570 is predicted to encode an exporter because of its homology to members of the major facilitator superfamily, including LmrB of Bacillus subtilis (28% identity and 45% similarity). Transcriptional regulation of SPy0570 by heat, acid, and CovR was confirmed by RPAs (data not shown).
SPy0229 and SPy0230 form an operon and encode products homologous to the MDR1 family of ATP-binding cassette (ABC) transporters. The closest relatives of SPy0229 and SPy0230 appear to be YdaG and YdaB of L. lactis (53% identity and 69% similarity and 55% identity and 72% similarity, respectively), which are also encoded by adjacent genes. In L. lactis, YdaG and YdbA function to export toxic substances such as ethidium bromide and daunomycin (33).
CovR is a repressor of SPy0229, a gene activated in response to temperature and pH.
To confirm our microarray results, we performed RPAs with newly isolated RNA from strain JRS4 to look at SPy0229 transcript levels following exposure to increased temperature and decreased pH (Fig. 3). There was a
4-fold increase in the amount of SPy0229 transcript following 5 min at pH 6.0, and this level decreased to
2-fold after 10 and 20 min at a low pH (Fig. 3A). Similarly, when cells were subjected to 40°C for 5 or 10 min, SPy0229 transcript levels increased (
1.5-fold) and then began to decrease by 20 min (Fig. 3B). These results confirm that transcription of SPy0229 increases transiently under conditions of heat and acid stress. Transient activation of this type is characteristic of many heat shock genes in other bacteria (43).
|
|
|
|
|
To be sure that the phenotype attributed to rscA resulted from the insertion mutation, a complementation experiment was performed. The gene directly downstream of rscA, SPy0230 (which we named rscB), has homology with ydbA of L. lactis and is thus likely to form a protein that complexes with RscA to form a functional transporter. The rscA insertion mutation in JRS1819 is likely to be polar. Therefore, we complemented the rscA insertion mutation in JRS1819 with pJRS334, a plasmid that contains both rscA and rscB under the lactococcal P23 promoter (39). We also constructed a control strain by transforming JRS1819 with the same vector without the rscA and rscB genes, pJRS7759. The EOP at 40°C of JRS1819 carrying the vector alone was similar to that of the rscA mutant (Fig. 7). However, the EOP of the complemented strain at 40°C was restored to that of the wild-type strain, confirming that the mutation in rscAB is responsible for the inability of JRS1819 to grow at high temperatures (Fig. 7).
| DISCUSSION |
|---|
|
|
|---|
Stress specificity of CovR-regulated GAS stress response genes. The purpose of this work was to identify stress genes that are CovR regulated, so we focused our analyses on those genes. Although we anticipated that many CovR-regulated stress response genes would be activated under both acid and heat stress conditions, we found that transcription of most of these genes was affected by only one of the two stress conditions studied. This suggests that the ability of CovR-regulated genes to respond differentially to different stresses may play an important role in the pathogenic potential of the GAS, and the identification of these differentially regulated genes may help identify factors needed at different stages of GAS disease, when the organism is exposed to a variety of conditions in the human host.
The specificity of the stress response of CovR-regulated genes may be mediated by several mechanisms. For instance, an additional stress-activated "pathway" may function to regulate CovR-specific promoters independently of CovR. This was observed to be the case for rscAB, since even in the absence of CovR, transcription is partially activated by lowering the pH (Fig. 5). This second pathway may be mediated by a regulatory protein whose synthesis is stress specific or by the degree of local DNA supercoiling at the promoter. In some bacteria, the degree of local superhelicity can be altered in response to low pH, high temperature, and increased osmolarity, which in turn may alter the level of transcription from certain promoters (2, 3, 5, 15, 24, 27, 30, 48).
In addition, the degree of CovR phosphorylation may be stress dependent. Stress conditions appear to result in the conversion of CovS to a phosphatase that dephosphorylates CovR, leading to the derepression of some CovR-regulated promoters (10). It seems likely that different stresses may affect the phosphatase/kinase ratio of CovS. Thus, the amount of phosphorylated CovR (CovR-P) may depend on the environmental condition in which the GAS is grown. Some promoters appear to have a greater differential affinity for CovR and CovR-P than others (20, 23); thus, the CovR/CovR-P ratio may lead to differential responses of some CovR-regulated promoters to specific stress conditions.
CovR-repressed proteins activated by heat and acid stress have homology to exporters. To identify genes responsible for the inability of a CovS mutant to grow at 40°C and at pH 6.0, we investigated the four genes that we found to be repressed by CovR and activated in response to both heat and acid. The product of one of these genes, SPy1404, has no homology to known proteins. However, the other two operons identified, rscAB and SPy0570, encode proteins homologous to members of each of the two major exporter protein families: the MDR1 class of the ABC transporter family and the major facilitator superfamily, respectively. Proteins in these two families are usually responsible for multidrug resistance in many different organisms, including gram-positive and gram-negative bacteria and humans (reviewed in reference 38). Furthermore, members of these protein families export unrelated small molecules besides antibiotics and other drugs (38). We found that SPy0570 is not required for growth at high temperature or low pH (data not shown), although it is possible that this may be due to the presence of a functionally redundant protein in GAS. In contrast, the MDR1 homolog RscA, while not required for growth at pH 6.0, is essential for GAS growth at 40°C.
MDR1 exporters. The members of the MDR1 class of the ABC transporter family mediate export using the energy of ATP hydrolysis. Proteins in this group, including RscA and RscB, contain at least two ATP-binding domains and two transmembrane domains (38). These proteins localize to the cytoplasmic membrane and function as either homodimers or heterodimers to move amphipathic substances including drugs and lipids across this membrane (reviewed in reference 12). The protein encoded by MDR1 itself, P glycoprotein, which has 28% identity and 47% similarity to both RscA and RscB, functions as a homodimer in humans. It exports a diverse assortment of amphipathic drugs, including anthracyclines, actinomycin D, valinomycin, and gramicidin, as well as lipids and steroids (12, 53). This ABC transporter also limits the success of cancer treatment due to its ability to export many drugs used for chemotherapy (32).
Many MDR1 family members are found in bacterial species where they are responsible for inherent resistance to multiple drugs that are not chemically related (37). In gram-negative bacteria, the MDR1 homolog msbA is an essential gene because, in addition to conferring drug resistance, it is responsible for the transport of lipid A, a major outer membrane constituent, across the cytoplasmic membrane (14, 41, 54). It seems likely that all MDR1 family members of different organisms are able to perform the same functions because of their very close structural relationship. Furthermore, the similarity of the MDR1 proteins of bacteria and humans is so great that LmrA from L. lactis can functionally be substituted not only for MsbA of E. coli but also for the human P glycoprotein (41, 51).
Like RscA, some other MDR1 class members have a demonstrated role in the bacterial response to heat stress. For example, MsbA of E. coli was first identified as a multicopy suppressor of a mutant that lacks the HtrB heat shock protein (28, 29). A deletion in htrB leads to the accumulation of lethal amounts of such precursor lipids as tetra-acylated lipid A and glycerophospholipids in the inner membrane at 42°C. The rescue of growth of an htrB mutant by the addition of msbA on a multicopy plasmid at high temperatures suggests that the MsbA exporter facilitates the removal of these precursor lipids from the inner membrane (36, 54). Since the MDR1 family members can export lipids, these transporters may protect the cell from heat stress by transporting lipids with specific fatty acid compositions to yield the appropriate membrane composition and fluidity for growth at high temperatures (7, 34, 47, 54).
In addition to their role in increasing thermotolerance, some MDR1 ABC transporters are involved in resistance to other stresses. For example, in the bacterium Oenococcus oeni, which is responsible for producing the flavors of wine, OmrA, a RscA homolog, is required for resistance of this organism to the low pH, ethanol, tannins, and free fatty acids found in wine. Production of OmrA is induced at high temperatures and by osmotic stress, consistent with its role in the stress response (7).
Potential roles for RscAB in GAS. Based on their homology to the MDR1 class of ABC transporters, along with the functional similarities between the members of these classes of proteins, we can propose several possibilities for the requirement of RscAB for GAS growth only under stress conditions. For example, like MsbA, RscAB may affect the movement of specific membrane lipids and therefore alter the fatty acid composition of the cytoplasmic membrane. Membrane structure and fluidity appear to be important for cell growth at high temperatures (7, 34, 47, 54).
Alternatively, the need for RscAB only when GAS is subjected to stress conditions might result from a functional requirement for these proteins to export toxic products. Such toxins might either accumulate under specific stress conditions or be produced constitutively during GAS growth but exert toxicity only under stress conditions. If RscAB is needed to export toxic products that are always produced in GAS, a homolog with a redundant function is expected to remove this toxin under nonstress conditions but should be nonfunctional under stress. A candidate operon for this function is SPy1148/SPy1149, which encodes predicted proteins that are 29% identical and 49% (SPy1148) or 50% (SPy1149) similar to RscAB. It is possible that these proteins are inactive at 40°C or are not synthesized at that temperature, making the export of a specific toxin, and thus growth of GAS, dependent on RscAB at 40°C. Likewise, proteins encoded by SPy1148/SPy1149 may be functional at a low pH, thus alleviating the need for RscAB under this condition. Further investigation of the survivors that we isolated when the rscA strain was plated at 40°C should help us understand the role of this exporter in the growth of GAS at high temperatures.
| ACKNOWLEDGMENTS |
|---|
This work was supported by Public Health Service grant R01-AI20723 from the National Institutes of Health. T.L.D. and J.T.C. were supported in part by NIH training grant T32-AI07470.
| FOOTNOTES |
|---|
Supplemental material for this article may be found at http://jb.asm.org/. ![]()
Present address: Division of TB Elimination, National Center for HIV, STD, and TB Prevention, Centers for Disease Control and Prevention, Atlanta, Ga. ![]()
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| Appl. Environ. Microbiol. | Infect. Immun. | Eukaryot. Cell |
|---|---|---|
| Mol. Cell. Biol. | J. Virol. | Microbiol. Mol. Biol. Rev. |
| ALL ASM JOURNALS |