Previous Article | Next Article ![]()
Journal of Virology, March 2008, p. 2324-2329, Vol. 82, No. 5
0022-538X/08/$08.00+0 doi:10.1128/JVI.01930-07
Copyright © 2008, American Society for Microbiology. All Rights Reserved.

Tanja M. Ruokoranta,
,
A. Marika Alapuranen,
and
Dennis H. Bamford*
Department of Biological and Environmental Sciences and Institute of Biotechnology, Viikki Biocenter, P.O. Box 56, 00014 University of Helsinki, Finland
Received 4 September 2007/ Accepted 3 December 2007
|
|
|---|
|
|
|---|
PRR1 is an icosahedral, single-stranded-RNA (ssRNA) phage (24) that is related to well-studied members (MS2, Qbeta, etc.) of the Leviviridae family (13). PRR1 is unique among members of the Leviviridae family in that it infects a wide range of gram-negative bacteria (8, 24), but only if they harbor an IncP-type conjugative plasmid. The PRR1 genome contains 3,573 nucleotides with only four protein-encoding genes, namely, those for the polymerase, the maturation and coat proteins, and the lytic factor (28). The PRR1 replication cycle and its effects on host cell physiology are described elsewhere (G. Daujotaitë, R. Daugelavièius, and D. H. Bamford, submitted for publication). The optimal bacterial host for PRR1 is Pseudomonas aeruginosa (24), a gram-negative, opportunistic human pathogen (15, 16, 18). At 6.3 million base pairs (containing 5,570 open reading frames [ORFs]), the P. aeruginosa genome is one of the largest bacterial genomes that has been sequenced (31). More than 9% of the known ORFs are classified as transcriptional regulators or two-component systems, reflecting the importance of transcriptional regulation in the ability of P. aeruginosa to respond and adapt to a myriad of environments (31). Interestingly, over 300 quorum-sensing control genes (
6% of the P. aeruginosa genome) have also been identified (29, 34). In this study, P. aeruginosa O1 (PAO1) DNA microarrays were utilized to identify genes whose expression levels were affected by infection with PRR1. When we monitored gene expression changes at various time points after infection, we found that viral infection did not grossly affect cell metabolism until the time of lysis. The genes most affected by PRR1 infection were grouped into three major functional classes, namely, transport, energy production, and protein synthesis. To determine whether these changes in gene expression are a general effect of phage infection or are specific to infection by PRR1, the data set was compared with a similar study of PRD1, a well-characterized icosahedral dsDNA model virus that utilizes the same cell surface receptor complex as PRR1 (1, 4, 9, 22, 23, 25).
|
|
|---|
Samples for RNA isolation were taken from both the infected and uninfected cultures at time points corresponding to 5, 10, 15, 30, and 50 min postinfection (p.i.). Additional samples were collected from the control culture at the time of infection (i.e., the 0-min sample). At each time point, two different samples were collected, one for RNA isolation and one for viable cell counts and A600 measurements. All samples were processed immediately. Three independent replicates of the experiment were performed.
Isolation of total RNA. Total RNA was isolated from 1.9 x 1010 cells (10 ml of the cell culture) with two phenol and three ether extractions. In the first extraction, the phenol (pH 4.3; Sigma) was preheated to 95 to 100°C and the sample was supplemented with 1 ml lysis buffer (0.5 M Tris-HCl, pH 7.5, 0.2 M EDTA, and 10% sodium dodecyl sulfate) to ensure rapid cell lysis. The nucleic acids were precipitated from the aqueous phase with sodium acetate and ethanol. Further purification was performed using an RNeasy Midi kit (Qiagen) and on-column DNase (Qiagen) digestion according to the manufacturer's instructions. The purified RNA samples were stored at –80°C. The concentration and purity of the RNA samples were determined by measuring the A260/A280 ratio. These readings were in the range of 2.0 to 2.2 for all samples. The integrity of the total RNA was confirmed by denaturing electrophoresis on formaldehyde-agarose gels, followed by staining with ethidium bromide.
Microarray data analysis. cDNA synthesis, labeling, hybridizations, and data collection were performed by CTL Bio Services (Rockville, MD) according to standard Affymetrix GeneChip procedures. Labeled cDNA was used to probe an Affymetrix GeneChip comprised of 5,900 spots that represented 5,769 probes for the PAO1 genome, including 5,570 predicted ORFs and rRNA and tRNA genes, as well as 199 probe sets corresponding to 100 intergenic regions. In addition, 14 probes served as controls, including genes from Bacillus subtilis, Arabidopsis thaliana, and Saccharomyces cerevisiae. The raw data and the annotated PAO1 genome were imported into GeneSpring GX software, version 7.3 (Agilent Technologies). GeneSpring GX was used to normalize the data prior to further processing. Several custom-made scripts were used to calculate the changes in gene expression levels within experimental replicates and between the infected samples and the control (uninfected) samples. Expression level changes were considered significant if they differed at least twofold and had a Student's t test P value of <0.05. We chose gene selection criteria that were lenient compared with those of our previous array experiments (27) because PAO1 displayed low variability in gene expression levels compared with cells infected with PRD1.
Microarray validation and qRT-PCR. Differential expression of selected genes was examined by quantitative real-time reverse transcription-PCR (qRT-PCR). In each case, qRT-PCR was performed on the same three biological replicates that were used for the microarray experiments. In addition, one to three technical replicates were carried out for each sample. RNA samples were treated with DNase (RQ1 DNase; Promega) and purified with an RNeasy MinElute cleanup kit (Qiagen). cDNA was generated using random hexamers and Moloney murine leukemia virus reverse transcriptase (Finnzymes). qRT-PCR was performed using a Dynamo SYBR green 2-step qRT-PCR kit (Finnzymes) according to the manufacturer's instructions. Primers were designed with the assistance of PrimerExpress software (Applied Biosystems). The following primer sequences were used: for PA4407 (ftsZ), 5'-CGAGCAGCAGTCGGTGAACTAC-3' (forward) and 5'-GTGAGACTGGTTGCGCATCA-3' (reverse); for PA3126 (ibpA), 5'-TCGGCCCTGCGTAATGAG-3' (forward) and 5'-GATACTCGTCGTCACCGTGCT-3' (reverse); for PA3183 (zwf), 5'-GAGCCGCCGCACTACATC-3' (forward) and 5'-TGGTCCTTGGTCATCACTTGC-3' (reverse); for PA0447 (gcdH), 5'-CGGCATCTCCGACGAATTC-3' (forward) and 5'-GTGCCCTCGTAGGTGTTCACC-3' (reverse); and for PA4268 (rpsL), 5'-CCAACGGTTTCGAGGTTTCC-3' (forward) and 5'-CGCTGTGCTCTTGCAGGTT-3' (reverse). All samples were normalized to the cell division gene ftsZ to control for the amount of cDNA in each reaction mix. Measurements were obtained using an ABI 7000 sequence detection system (Applied Biosystems) and a program of 15 min at 95°C and 40 cycles of 10 s at 94°C, 30 s at 60°C, and 30 s at 72°C. qPCR product dissociation curves were used to verify the specificity of the amplified product. Melting curves for PCR specificity analysis were obtained following the completion of the amplification reactions, using the following conditions: 15 s at 95°C, 20 s at 60°C, and 15 s at 95°C. Raw data were analyzed using ABI Prism 7000 software (Applied Biosystems), and relative quantification was performed using Microsoft Excel. The 0-min samples were used as calibration samples for the cycle threshold method (relative quantifications).
|
|
|---|
![]() View larger version (10K): [in a new window] |
FIG. 1. Experimental protocol. P. aeruginosa O1 cells were cultured in Luria-Bertani broth, with aeration, to exponential growth phase. The cultures were divided into two parts, and one of the subcultures was infected with phage PRR1 (multiplicity of infection of 10). Turbidity values (OD600) from three independent experiments, plotted as a function of time (dashed lines, infected cells; solid lines, uninfected cells), are depicted. The time point for virus addition is indicated by an arrow. Samples were collected for total RNA isolation at six time points (0, 5, 10, 15, 30, and 50 min p.i. [black circles]) and at equivalent time points for uninfected cultures.
|
![]() View larger version (22K): [in a new window] |
FIG. 2. Verification of microarray results by qPCR analysis. Comparisons of mRNA levels were performed with noninfected and PRR1-infected cells by microarray analysis and qPCR for selected P. aeruginosa genes. Genes are indicated at the left side of each panel. Solid lines, uninfected x-min samples versus uninfected 0-min sample (Ex/E0); dashed lines, infected x-min samples versus uninfected 0-min sample (Px/E0).
|
|
View this table: [in a new window] |
TABLE 1. Global changes in PAO1 gene expression during PRR1 infectiona
|
PRR1-infected cultures have been shown to produce up to 1,000 infectious particles per cell. They also produce semicrystalline, intracellular aggregates, estimated to contain as many as 4,000 viral particles (Daujotaitë et al., submitted). It is also possible that a fraction of the viral particles produced are noninfectious. Thus, the total number of viral particles produced by each cell may be considerably larger than the estimated number of infectious particles produced. However, owing to the small mass (
4 MDa) of members of the Leviviridae family, the fraction of cellular synthesis of macromolecules directed to virus production is likely still very small. For example, in the case of PRD1 (
66 MDa), approximately 200 particles are produced per cell. Since the mass ratio between PRD1 and PRR1 is approximately 18, the approximately 4,000 PRR1 particles per cell would be roughly equal in terms of viral mass to the mass produced by a PRD1-infected cell. This total viral mass corresponds to only about 10% of the total cellular macromolecule pool (27). Consequently, PRR1, like PRD1, appears to replicate in the cell with little interference to the host cell life cycle until late in infection (30 min p.i.). This view differs from a previous model that proposed that the entire synthesis capacity and cellular macromolecule pool are harnessed for virus production (12). However, our observation that both cell growth and host gene expression are virtually unchanged during the first 30 min p.i. is consistent with the notion that PRR1, as well as PRD1 (27), has minimal effects on the host cell. In fact, even as cells succumb to PRR1-directed lysis, the number of genes that show expression level changes is fewer than the number of genes that change as cells prepare to enter stationary phase. Thus, for most of the PRR1 infectious cycle, its effects on host cell physiology appear to be very small.
Gene-specific changes due to cellular growth and viral infection. We next compiled a list of genes that showed the greatest expression level changes and compared their annotated functions. This analysis revealed that the genes most dramatically affected by PRR1 infection tend to be involved in one of three major cellular functions, i.e., transport, energy production, and protein synthesis (Table 2). Genes involved in these processes represent 69 of the 190 genes whose expression levels changed in response to PRR1 infection. Further analysis of the microarray data set revealed that in almost all cases, several genes in a single operon were coregulated. Notably, PRR1 infection caused the up-regulation of the heat shock protein gene ibpA beginning at 30 min p.i. Late in the infection, the clpC (encoding the ATP-binding subunit of an ATPase with chaperone activity), rpo, and yhb (both involved in transcription) genes were all down-regulated. In contrast, at equivalent time points in uninfected control cells, these three genes were all up-regulated (Table 2, "Others"). When the uninfected cells reached late logarithmic growth phase (at 50 min p.i.), the expression of genes encoding basic metabolic pathways (e.g., glycolysis, tricarboxylic acid cycle, Enter-Doudoroff pathway) was down-regulated. However, only minor changes to this pattern were observed in infected cells at an equivalent time point. The general trend was that such pathways were further down-regulated.
|
View this table: [in a new window] |
TABLE 2. Regulation of selected gene groupsa
|
(ii) Energy production/conversion. Expression of the ccoO (cytochrome oxidase) gene was down-regulated late in infection but was up-regulated in uninfected control cells. The nir and pse genes were regulated similarly in both infected and uninfected cells. The most prominent effect of PRR1 infection was found among the nar genes, which were up-regulated.
(iii) Ribosomal proteins. As a whole, PRR1 infection had the most significant effect on the genes that encode ribosomal proteins. Up to 38 genes in this category were affected. In uninfected control cells, several ribosomal protein genes were up-regulated at the 15-min time point, followed by a marked down-regulation at 50 min. Surprisingly, in PRR1-infected cells, this strong down-regulation of ribosomal protein genes occurred at the 15-min time point. A similar pattern of expression was observed for the tsf and tuf genes, which encode translation elongation factors that are used as replication cofactors by the viral polymerase in leviviral infections (33). At this time point, viral replication is complete and coat protein synthesis has just begun (as shown for phage f2 [5, 32]). The increase in culture turbidity terminated at 30 min p.i. (Fig. 1). These observations are consistent with the notion that only a small fraction of the cell's protein synthesis capacity is needed to produce virions late in infection and, therefore, that no boost in ribosome biogenesis is required.
Comparison of PRR1- and PRD1-induced host gene expression changes.
PRR1 is a small icosahedral ssRNA virus, whereas PRD1 is a complex, membrane-containing icosahedral dsDNA virus. However, both depend on the same membrane-associated conjugative tra complex for cellular entry and, consequently, share the same broad host range (10). The durations of the viral life cycle in liquid medium are approximately the same for these viruses (cell lysis starts at 50 to 60 min p.i.), and in both cases, similar total virion masses are produced by infected cells. A surprising observation with the PRD1 system was that the virus did not cause significant changes to host cell metabolism (only
8% of genes were found to be down- or up-regulated due to viral infection, compared with
9% of genes which were changed due to physiological changes) (27). In the case of PRR1, an even smaller number of host genes responded to viral infection (
4% due to virus infection and
5% due to physiological changes). To determine whether the host responds universally to viral infections, we compared the subset of genes that responded to PRR1 infection to those that responded to PRD1 infection. Overall, we found very little overlap between specific genes that showed expression level changes between these two systems. However, both viruses caused down-regulation of protein synthesis, albeit at different time points (15 min p.i. for PRR1 and 30 min p.i. for PRD1). A large number of heat shock genes were strongly up-regulated during PRD1 infection, whereas in the case of PRR1, only ipbA expression was significantly increased. Thus, this analysis suggests that the down-regulation of protein synthesis may be the only universal response to viral infection and may be required to maintain the energy needed for virus production, since protein synthesis is an energy-expensive process.
|
|
|---|
However, we have provided evidence that the PRR1 phage, like PRD1, replicates within bacterial cells with minimal effects on the major biosynthetic pathways of the host, other than protein synthesis. Since no drastic changes were observed, it seems that host cells act merely as passive vessels for virus multiplication. In this respect, PRR1 is even more invisible to the host than PRD1. This analysis also demonstrates that each virus affects the host in unique ways, since the set of genes altered by PRR1 shows little overlap with those altered by PRD1.
This work was supported by the Finnish Center of Excellence Program (2006-2011) of the Academy of Finland (grants 1213467 and 1213992 to D.H.B.).
Published ahead of print on 12 December 2007. ![]()
J.J.R. and T.M.R. contributed equally to this work. ![]()
Present address: Vactech Oy, Biokatu 8, 33520 Tampere, Finland. ![]()
Present address: Roal Oy, P.O. Box 57, 05201 Rajamäki, Finland. ![]()
|
|
|---|
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»