Previous Article | Next Article ![]()
Journal of Virology, May 2009, p. 4316-4325, Vol. 83, No. 9
0022-538X/09/$08.00+0 doi:10.1128/JVI.02075-08
Copyright © 2009, American Society for Microbiology. All Rights Reserved.

Department of Microbiology, Molecular Biology, and Biochemistry and Center for Reproductive Biology, University of Idaho, Moscow, Idaho
Received 2 October 2008/ Accepted 6 February 2009
|
|
|---|
|
|
|---|
As a member of the betaherpesvirus family, HCMV has a linear, double-stranded DNA genome of approximately 230 kb, containing 198 known open reading frames (ORFs) (gene image map for human herpesvirus 5 [http://www.oralgen.lanl.gov/oralgen/bacteria/hhv5/]). These ORFs are temporally expressed in three main stages (19, 47, 48). Immediate-early (IE) gene expression begins soon after the viral genome is deposited into the nucleus of the host cell. There, IE gene products drive the expression of early (E) genes, which primarily code for proteins responsible for replicating the viral genome. Once viral genome replication has progressed sufficiently, late (L) gene expression commences. L genes encode components of the capsid, the tegument, and the envelope of the virions.
To ensure precise timing of gene expression and a favorable environment in which to replicate, transcriptional cross talk between the virus and its host cell has evolved. In this cross talk, the virus utilizes cellular transcription factors such as AP-1, USF, E2F, TFIID, and Elk-1 to enhance viral gene expression (10, 24, 25, 46, 54). The viral IE transcription factors IE1 and IE2 indirectly alter host gene expression through protein-protein interactions, leading to multiple outcomes including the transition of cells from G0 to S phase, inhibition of cellular DNA replication, and prevention of apoptosis (1, 37, 55, 58).
As a key cell cycle regulator, the p53 protein is known to regulate the transcription of a variety of genes, including those involved in cell cycling, DNA repair, and programmed cell death (reviewed in reference 40). Besides transcriptional activation, p53 is an important player in the DNA repair machinery of the cell, working primarily via protein-protein interactions. It is stabilized in a normal cell following stress conditions such as UV irradiation or gamma-irradiation, hypoxia, starvation, extreme heat, or viral infection (reviewed in reference 34). This stabilization leads to the activation of the p53 protein and to downstream effects on the expression of cellular proteins, leading to cell cycle arrest or apoptosis. p53 can activate genes encoding p21/Waf1/Cip1, 14-3-3
, GADD45, and PCNA to arrest the cell cycle and allow DNA damage repair, or the apoptosis-promoting Bax, Noxa, PERP, PIG3, and Fas genes can be upregulated by p53 stabilization (26, 44). p53 is also known to function as a transcriptional inhibitor, and other groups have documented the repression of genes after p53 induction (31, 32). In addition, the C-terminal nonspecific DNA binding domain can modulate the binding of p53 to its specific binding site (BS) by protein modifications such as phosphorylation, methylation, or sumoylation, thus altering gene expression by conferring DNA structure specificity (11).
p53 contains four main functional domains (from the N to the C terminus): a protein-protein interaction domain, a specific DNA binding domain, a homotetramerization domain, and a nonspecific DNA binding domain. The protein-protein interaction domain is located on the N terminus of the protein. This domain interacts with proteins such as MDM2, TFIID, and the adenovirus protein E1B. It also serves as a transcriptional transactivation domain (9, 29, 57). The binding of MDM2 and E1B to this domain suppresses the transactivating function of p53, whereas the binding of the TATA box binding protein TFIID facilitates its function.
The specific DNA interaction domain is crucial to the transactivating ability of p53. Transcriptional activation or repression occurs when p53 binds to its responsive element (RE). The RE consists of the sequence RRRCWWGYYY (where R is a purine, Y is a pyrimidine, and W stands for A or T), followed by another RE located 0 to 13 nucleotides from the first (16). General features of strong REs include the conservation of the internal C and G residues and the presence of no more than three nonconsensus bases within the sequence. The RE configuration containing an AT sequence between the C and G residues at the middle of each half-site exhibits the strongest p53 binding (39, 53). Usually, transcriptional activation/inhibition of genes takes place within a few thousand nucleotides of the RE (reviewed in reference 28).
The third domain consists of a C-terminal helix-plus-basic-region motif used for the homotetramerization of the p53 protein, since p53 is most active in transcription as a tetramer (51). Besides the homotetramerization domain, the C terminus also contains the fourth functional domain, the nonspecific DNA binding domain, which is considered responsible for the DNA damage-sensing abilities of p53 (33). The nuclear import and export signals are also located in the C terminus.
Several studies involving p53 and HCMV have shown that p53 levels are rapidly stabilized following infection with HCMV in several cell types, including fibroblasts, vascular smooth muscle cells, human umbilical vein endothelial cells, and astrocytes (23, 27, 30, 36, 49). Although p53 levels are elevated upon viral infection, the cellular downstream targets of p53 are not activated (4, 23). However, p53 appears to play a functional role in the viral life cycle of HCMV. Huang and colleagues have shown that p53 is involved in the regulation of the viral protein UL94 (56). In addition, we have found that cells infected with HCMV in the absence of p53 produce fewer infectious viral particles (7). Examination of viral protein production (by Western blotting) and trafficking (by immunofluorescence) revealed delays in these processes in cells lacking p53. A lag in the formation of replication centers, potentially delaying viral encapsidation, was also observed (7).
Our earlier studies also showed that p53 was sequestered into the viral replication compartments and that its specific DNA binding domain was required for the sequestration of the protein into these centers (17, 41). When the HCMV genome was examined for p53 REs, 21 potential sites were found (41). Chromatin immunoprecipitation (ChIP) analyses revealed that p53 was temporally bound to some of these REs within the viral genome. Several of these REs were located in or near genes known to be essential for viral reproduction in tissue culture (15).
In this study, microarray analyses revealed that the expression levels of 22 viral genes were affected by the absence of p53. We have confirmed this observation for the 11 most highly affected genes by real time reverse transcriptase, or quantitative, PCR (q-PCR). In addition, q-PCR analyses at a multiplicity of infection (MOI) of 5 uncovered p53-dependent regulation of four of six genes shown to have p53 bound to their corresponding p53 BSs located in the vicinity of their start codons (41). q-PCR time courses for four viral genes regulated at a low MOI in microarray analyses indicated that, at both low and high MOIs, p53 exerted its greatest influence on gene expression at 24 h postinfection (hpi). To determine whether direct interaction of p53 with the defined BSs was affecting gene expression during infection, luciferase reporters were constructed using UL65, UL75, and UL122 upstream sequences. Varying effects were seen depending on the time postinfection (p.i.) and the particular site tested. Both UL75 and UL122 upstream sequences showed that the presence of intact p53 BSs, as well as cellular effects due to the p53 knockout environment, influenced promoter activity. These influences are most probably due to a combination of direct and indirect interactions of p53 with the virus and its genome. We propose that the presence of p53 is important for both IE and E viral gene expression and that it enhances the fitness, survival, and reproduction of the virus within the host. We also provide an explanation for the delay in viral replication and reproduction in 53–/– cells.
|
|
|---|
Virus infection. Cells were synchronized in G0 by serum starvation for 60 h (7), trypsinized, and replated onto 100-mm-diameter dishes at a density of 1 x 106 cells per plate for microarray and q-PCR analyses. Cells were then infected with the HCMV strain Towne 3 h after reseeding at an MOI of 1 for the microarray studies and at MOIs of 1 and 5 for q-PCR experiments. Virus-containing media were replaced by fresh media at 3 hpi; the cells were harvested by trypsinization; and the pellets were stored at –80°C until use.
Microarray analysis. For microarray analysis, HFF, THF, and p53–/– fibroblasts were infected as described above and harvested at 12 and 24 hpi. RNA was extracted from the harvested cells with the SV Total RNA isolation kit (Promega) and was concentrated in a SpeedVac concentrator. RNA was amplified with the Superscript RNA amplification kit (Invitrogen), and the amplified RNA was labeled with Cy5 using the ULS labeling kit for CombiMatrix arrays (Kreatech Biotechnology).
The labeled RNAs obtained at 12 and 24 hpi from the infected HFF, THF, and p53–/– cells were hybridized to six custom-designed CombiMatrix microarray slides using protocols provided by CombiMatrix. Each of the four identical arrays on a slide consisted of 2,240 features (305 control features, 1,265 cellular gene features, and 670 viral gene features). Every gene on the array was represented by at least three features.
Three arrays on a slide were hybridized with three different Cy5-labeled RNA samples, gathered as biological replicates from three independent infection experiments. The fourth array served as a control and was hybridized with an equimolar mixture of the three biological replicate RNA samples used for the other three arrays on each slide. The data from these control arrays were consistent with the measurements of the three biological replicates but were not used for further calculations.
The fluorescence of each feature on the hybridized slides was read on an Axon GenePix 4200A microarray scanner (Axon, Sunnyvale, CA), and raw fluorescence values were obtained using the CombiMatrix microarray imager. The slides were then stripped using the protocol provided by CombiMatrix and were reprobed with the same RNAs used in the first round of hybridization. This stripping procedure was repeated twice to provide technical replicates, reducing protocol handling and device variability.
Data analysis. The raw data were analyzed as follows. First, a mean value was calculated for each feature on the array from the three technical replicates. The variance stabilization (22) was calculated separately for the 12-hpi and 24-hpi values. The variance-stabilized values were used for six crosswise comparisons of the three cell types at the two time points (HFF versus THF, HFF versus p53–/– cells, and THF versus p53–/– cells, at 12 hpi and 24 hpi). The differences between the logarithmic values of corresponding features on different slides were calculated, and the t value of a two-tailed Student t test with a 95% confidence interval was used to determine the P value of every feature. The t value was generated from the three biological replicate values. For each individual gene, the average regulation of all features of the gene was calculated. The expression of a gene was considered significantly different in a particular comparison based on two criteria: at least 50% of the gene's features had to be significantly different in the t test (P < 0.05), and the gene had to have a >2-fold difference between the cell types, averaged over all the features of the gene.
q-PCR. RNA was extracted from cells with the SV Total RNA isolation kit (Promega) and was reverse transcribed with Superscript II (Invitrogen) or ImProm-II reverse transcriptase (Promega). RNA samples were used for reverse transcription only if no contaminating genomic DNA could be detected by using oligonucleotide primers for IE1 in a 40-cycle PCR. q-PCR was performed on an ABI Prism 7900 system (Applied Biosystems) by using the Power SYBR green PCR master mix (Applied Biosystems) in 15-µl reaction mixtures for 40 PCR cycles. Calculations were based on absolute starting quantities, using reactions specific for glucose-6-phosphate dehydrogenase (G6PD) as normalization controls (45). Each measurement consisted of two biological replicates, gathered from independent infection experiments. In addition, each biological replicate consisted of three technical replicates, and their averaged values were used for downstream calculations.
Molecular cloning. All plasmids used for q-PCR were generated with the pCR8 TOPO cloning kit (Invitrogen). Taq polymerase was used to amplify the insert fragments used as standards in the q-PCR experiments from the cDNA of HCMV-infected HFF. (cDNA was generated as described in the preceding section.) The sequences of the oligonucleotides used to amplify the insert fragments are available upon request.
For luciferase assays, plasmids carrying firefly luciferase under the control of either 2 kb of UL75 upstream sequence, 2 kb of UL122 upstream sequence, or 1 kb of UL65 upstream sequence were generated by amplifying the corresponding sequences from viral genomic DNA using Phusion polymerase (New England Biolabs) (primer information is available upon request). The constructs were designed with 600 bp 5' to the p53 BS located closest to the start codon of these genes. For the UL65 upstream sequence, the distance between the p53 BS and the start codon was 400 bp; for the UL75 and UL122 upstream sequences, this distance was 1,400 bp. Amplicons were cloned into ptubFLuc (a kind gift from Gustavo Arrizabalaga) via the NsiI and KpnI sites present in both the oligonucleotide primers and 5' of the luciferase ORF in ptubFLuc. The p53 BS in each of these plasmids was mutated by site-specific PCR mutagenesis using QuikChange (Stratagene) (primer information is available upon request). Site-specific mutations in these plasmids were checked by sequencing. As a normalization control in luciferase assays, a plasmid carrying Renilla luciferase under the control of the simian virus 40 (SV40) enhancer element was constructed. This plasmid was created by ligating the BglII/blunt-ended-HindIII 400-bp fragment of pGL4.13 (Promega) into the BglII/blunt-ended-PstI sites of pRL-CMV (Promega). This procedure exchanged the CMV enhancer for the SV40 enhancer element in this vector. Blunt ends were generated using T4 polymerase (Fermentas).
Luciferase assays.
A total of 2 x 106 THF and p53–/– cells were transiently transfected with 10 µg of one of the firefly luciferase plasmids mentioned above and 2 µg of the Renilla luciferase vector by electroporation at 330 V, 2,500 µF, and 75
in a BTX ECM630 electrocell manipulator. Cells were seeded in 60-mm-diameter plates after electroporation and were infected with HCMV at an MOI of 5 at 48 h posttransfection. Cells were harvested at 24 or 72 hpi by trypsinization. Luciferase measurements were performed using the Dual Luciferase reporter assay system (Promega). Graphs show the mean ratio between the firefly luciferase and Renilla luciferase readings from two independent experiments.
Microarray data accession number. The microarray data, including the list of genes on the arrays and the sequences of the spotted oligonucleotides, have been deposited at the GEO database under accession number GPL7474.
|
|
|---|
![]() View larger version (31K): [in a new window] |
FIG. 1. Distributions of microarray features of cross-comparisons of HHF, THF, and p53–/– cells by the log2 of the n-fold regulation (on the abscissa) and the negative log10 of the P value generated from Student's t test (on the ordinate), depicted as volcano plots. Each feature is represented by a dot, with black dots representing cellular genes, red dots representing controls, and green dots representing viral genes. Significant features are represented by dots above the horizontal line at 1.3 (indicating a P value of <0.05) and left or right, respectively, of the vertical line at –1 or the vertical line at 1 (indicating a >2-fold change in one or the other cell type). The relative numbers of features in quadrants are given as percentages.
|
A gene was determined to be significantly different in a particular comparison if the averaged difference over all of this gene's features was >2-fold and at least 50% of its features (with every gene being represented by at least three features on the array) had P values of <0.05. (For example, gene X had five features on the array. The changes determined for these features were 2.5-fold, 2.6-fold, 2.3-fold, 1.9-fold, and 2.8-fold, which averaged to 2.42-fold, meeting the first significance criterion. However, the t tests for each of gene X's five features resulted in two features with P values of <0.05 and three with P values of >0.05, i.e., only 40% of gene X's features were considered to be significantly different. Therefore, gene X would not be regarded as significantly changed. If three rather than two of gene X's features had had P values of <0.05, gene X would have been considered significantly regulated.) These conditions were chosen to minimize the number of false-positive results. Because our experiments concentrated on viral genes and their expression, cellular genes were not further analyzed in this study.
At 12 hpi, by application of the significance criteria defined above, no genes showed a significant difference in any of the comparisons. At 24 hpi, we found 19 genes to be significantly different in the THF-versus-p53–/– cell comparison and 14 genes to be significantly different in the HFF-versus-p53–/– cell comparison. Eleven genes were shared between HFF and THF compared to p53–/– cells (UL30, UL54, UL63, UL65, UL69, UL75, UL110, UL122, IRL4, IRS1, and TRS1) (Table 1). Eight genes were unique in the THF-versus-p53–/– cell comparison (UL48, UL68, UL84, UL100, UL106, UL109, UL111, and RL5A), and three genes were unique in the HFF-versus-p53–/– cell comparison (UL14, UL16, and US2). The 22 separate genes found in these comparisons cover
10% of all 198 known ORFs in the HCMV genome. When HFF and THF were compared, only US2 showed a significant difference. The expression of US2, a glycoprotein responsible for proteasomal degradation of major histocompatibility complex proteins, might therefore be dependent on the cell line (52). The average ratio of the expression level in the HFF-versus-p53–/– cell comparison to the expression level in the THF-versus-p53–/– cell comparison for the 11 shared genes was 1.2, indicating that, with the exception of US2, these two cell lines behaved similarly with respect to p53-influenced viral gene expression.
|
View this table: [in a new window] |
TABLE 1. Significantly different genes in cross-comparisons in the microarray analyses
|
In q-PCR experiments (Fig. 2) performed at an MOI of 1, the largest changes were found in UL64, UL65, UL63, U75, UL110, and UL122: 13.9-fold, 11.0-fold, 7.0-fold, 6.5-fold, 4.0-fold, and 3.6-fold, respectively (Fig. 2). With the exception of UL64, these genes all belong to group 1. Taking into account the large standard deviations of some measurements shown in Fig. 2, differences of <3-fold in q-PCR experiments were regarded as no change. All group 3 genes analyzed as controls (UL18, UL41, and US15) were considered unchanged in p53–/– cells compared to THF at an MOI of 1. The transcription of the five group 2 genes (UL45, UL57, UL97, UL105, and US31) was also unchanged at an MOI of 1 in q-PCR experiments. This finding was consistent with our microarray results. The changes in the expression of the other six group 1 genes (UL48, UL54, UL69, UL84, TRS1, and US2), although considered significantly different in the microarray studies, fell below the threshold of threefold, and thus their expression was considered not significantly different in p53–/– cells versus THF (the average change in the expression of these genes was 1.5-fold in the q-PCR experiments and 2.6-fold in the microarray studies). However, the group 1 genes with the highest changes in the microarray studies showed similar results in q-PCR analysis.
![]() View larger version (22K): [in a new window] |
FIG. 2. q-PCR experiments conducted at 24 hpi on 20 viral genes at MOIs of 1 (light shaded bars) and 5 (dark shaded bars) in THF and p53–/– cells. Each bar represents the average from two independent experiments, and error bars indicate the minimum and maximum values of the single measurements. The expression of the first 11 genes (UL48 to US2) was shown to be significantly different in the microarray experiments, while the 6 genes in the middle group (UL45 to US31) were chosen to verify the results of the ChIP experiments performed by our group previously (41). The three genes in the last group (UL18 to US15) served as control genes, with no p53 BS within 3 kb of the start codons. The ordinate is scaled in a log10 manner, with negative values indicating downregulation in p53–/– cells.
|
|
View this table: [in a new window] |
TABLE 2. Transcript levels of genesa relative to G6PD transcript levels in THF at MOIs of 1 and 5 at 24 hpi
|
5-fold beginning at 48 hpi. UL122 showed no difference between cell types beginning at 72 hpi. This confirmed that the presence of p53 influenced viral gene expression more significantly at the IE and E stages of infection. These results also suggest that as an infection progresses, viral gene expression becomes less dependent on the presence of p53.
![]() View larger version (24K): [in a new window] |
FIG. 3. Expression kinetics of UL63, UL65, UL75, and UL122 as determined by q-PCR over a time course of 96 hpi at MOIs of 1 (light shaded bars) and 5 (dark shaded bars). Each bar represents the average from two independent experiments, and error bars indicate the minimum and maximum values of the single measurements. The ordinate is scaled in a log10 manner, with negative values indicating downregulation in p53–/– cells.
|
The first sequence analyzed was the UL65 upstream sequence, which showed differences of <1.5-fold between cell types or sequences with wild-type versus mutated p53 BSs. UL65 was therefore regarded as unchanged at both 24 and 72 hpi (Fig. 4). Although we had not seen binding of p53 to this site in the ChIP assays (41), this finding was rather surprising. UL63 and UL65 were among the genes with the highest changes in the microarray and q-PCR studies. The luciferase activity of the UL75 upstream sequence was, to a small degree, dependent on the integrity of the p53 BS, with 3.7-fold and 2-fold differences in expression between the wild-type and mutant constructs in THF at 24 hpi and 72 hpi, respectively. However, the differences observed were not due to p53 binding, since the largest differences (4.6-fold) were seen in p53–/– cells at 24 hpi. The changes seen here appear to be due to indirect effects of p53. The third sequence we analyzed was the 2-kb UL122 upstream sequence. The intact p53 BS negatively influenced the promoter activity of the UL122 gene in THF at 72 hpi, as evidenced by the 1.9-fold increase in activity after mutation of the site. This was not seen in p53–/– cells. Luciferase activity was 3-fold to 4.6-fold lower in p53–/– cells than in THF with both the wild-type and the mutated construct. This could point to binding of activating transcription factors other than p53 to this site. Overall, although the three p53 BSs we analyzed behaved differently, the observed direct influence of p53 on viral gene expression was only minor. Therefore, the major changes in transcriptional activation observed in these experiments were caused by secondary effects due to the p53 knockout environment. Our data showed that while regulation of viral gene expression was influenced by the presence of p53, intact p53 BSs in promoters of viral genes are of only minor significance to the regulation of genes observed in microarray and q-PCR analyses.
![]() View larger version (22K): [in a new window] |
FIG. 4. Luciferase assays using upstream sequences of UL65, UL75, or UL122 to drive the expression of firefly luciferase. Asterisks in construct designations indicate the presence of a mutated p53 BS in these sequences. Renilla luciferase expression was used as a normalization control, yielding relative light units (RLU) as the dimension on the ordinate. Each bar represents the average from two independent experiments, and error bars represent the minimum and maximum values of the single measurements.
|
|
|
|---|
The genes found to be significantly changed were classified based on their gene expression kinetics as defined by Chambers and colleagues (8). These genes included 4 IE genes (UL122, IRS1, TRS1, and UL110, corresponding to 50% of known IE genes), 10 E and E-L genes (UL16, UL54, UL63, UL69, UL75, UL84, UL100, UL106, IRL4, and US2, corresponding to 10% of known genes with E and E-L kinetics), and 6 genes with L kinetics (UL14, UL48, UL65, UL68, UL109, and UL111, corresponding to 13% of known L genes). Chambers and colleagues did not provide a timetable for the expression of the genes analyzed, but ganciclovir was used to distinguish L from IE and E genes at 72 hpi (8). Considering that replication centers are clearly visible as early as 24 hpi in HFF and THF, this does not preclude the (perhaps unstable) expression of L genes at 24 hpi (7). No information is available on the gene expression kinetics of UL30 and RL5A. Interestingly, all IE genes identified were described as transcriptional activators or enhancers. No IE apoptosis inhibitors, such as UL36 and UL37, were affected by the presence of p53 (1). This may explain the delay in viral gene expression at the IE and E stages of infection seen previously in p53–/– cells (7).
One of the most striking differences we had observed in earlier studies was the delay in viral DNA accumulation in p53–/– cells (7). q-PCR experiments revealed that the transcript levels of UL54, the catalytic subunit of the viral DNA polymerase, were downregulated twofold at an MOI of 1 (which was considered insignificant) and fourfold at an MOI of 5 in p53–/– cells. This was a moderate downregulation compared to that of other genes, such as UL63 (46-fold) or UL65 (55-fold), at an MOI of 5. However, it could help to account for the delay in viral replication when coupled with our previous data showing that other replication proteins, such as UL44 and UL57, were downregulated at the protein level in p53–/– cells (7; also unpublished data). Starisky and Hayward showed that the gene products of UL84, UL122, and UL69 (a regulator at the pre- or posttranscriptional level) are crucial for viral replication (43). In high-MOI infections, these three genes were also downregulated in p53–/– cells, as determined by q-PCR. The decreases in the levels of UL54, UL69, UL84, and UL122 transcripts could likely impair and delay viral replication in p53–/– cells (7). In addition, the transcriptional expression of five envelope glycoproteins was found to be decreased in p53–/– cells in our microarray study: UL16 (an immunoevasin), US2, UL75 (envelope glycoprotein H), UL100, and RL5A (which belongs to the RL11 glycoprotein family) (14). The damaged appearance of virions shed from p53–/– cells, observed in our unpublished electron microscopy studies, and the poor infectivity of these shed virions may be explained by the lack of glycoprotein transcripts. Indeed, previous protein studies have shown delays in tegument, capsid, and glycoprotein expression in p53–/– cells (7 and our unpublished results). Since the functions of only 12 of the 22 genes revealed in our microarray studies are known, many more processes vital to viral reproduction, such as protein particle trafficking within the cell (42), remodeling of the nuclear membrane (6), processing of replicated DNA (3), or immune response modifications (reviewed in reference 35), could, in an in vivo situation, be impaired by the observed changes in gene expression.
To relate the low-MOI microarray experiments more directly to our earlier experiments performed at a high MOI, we performed q-PCR analyses at both high and low MOIs. We found that transcription was significantly changed in a larger number of genes when infections were performed at an MOI of 5. In fact, almost all group 1 and group 2 genes were affected by the absence of p53 protein at a high MOI. It is possible that a decrease in the levels of IE transcripts in p53–/– cells results in a disproportionate reduction in the level of E transcription from the more numerous input viral genomes present at a higher MOI due to transcriptional feedback mechanisms. In contrast to these group 1 and 2 genes, the group 3 control genes showed similar expression patterns at MOIs of 1 and 5. At both MOIs, the expression of these three genes was considered unregulated.
A plausible mechanism for the transcriptional activity of p53 in the viral genome is direct binding to the viral DNA. Fourteen of the 21 p53 BSs in the viral genome were loaded with p53 protein at various times p.i. in our earlier ChIP experiments (41). To test several of these sites, luciferase reporter constructs were engineered from upstream sequences of three viral genes, containing both wild-type and mutant p53 BSs. Reporter assays were then conducted in the presence and absence of p53 to investigate if the transcript changes seen in q-PCR and microarray analyses were due to intact p53 BSs and/or to secondary effects of the absence of p53. We tested the upstream sequences of UL65, UL75, and UL122. A unique cluster of three overlapping sites exists in the UL64 ORF that may affect the transcription of UL63, UL64 itself, and UL65. Although we had not observed binding of p53 to this site in our ChIP experiments, at a high MOI these three genes showed the largest downregulation in our q-PCR experiments. Further evaluation of p53 BSs in the viral genome revealed 10 additional sites, with one mismatch in each RE. These were in addition to our original group, each of which contained one perfect RE and one with as many as three mismatches. The BS in the UL75 upstream sequence represents 1 of the 10 newly identified potential BSs. The BS in the UL122 upstream sequence represents 1 of the 21 sites previously identified. This site was shown previously to bind p53 during at least one time point during infection (41).
Transcription from the UL65 upstream sequence with its unique BS was not affected at all by the presence of p53 or its BS in these assays. In spite of the significant regulation of UL63 and UL65 in both microarray and q-PCR analyses, this result was consistent with our ChIP data, which did not show binding of this site at any time point analyzed. It should be noted that the UL65 upstream sequence consisted of only 1 kb, whereas the UL75 and UL122 sequences were 2 kb. Thus, the binding of other viral or cellular transcription factors further upstream may be necessary to explain the effects on transcription observed in q-PCR and microarray analyses. It is also possible that a larger genomic context of the p53 BS may be required in order to observe direct regulation of this site by p53. Regardless, in this context, the transcriptional effects on UL63 and UL65 appear to be due primarily to indirect effects of the absence of p53.
Luciferase assays using the UL75 upstream sequence revealed that the integrity of the p53 BS was of minor importance for transcriptional activity. However, p53 binding did not account for these changes. In addition, these assays indicated a potential time-dependent effect of the presence of p53 on the UL75 late promoter.
In luciferase assays using the UL122 upstream sequence, the downregulation observed in microarray and q-PCR experiments was reflected in decreased transcriptional activity in the p53 knockout environment in both wild-type and mutant BS configurations. It was surprising that mutation of the p53 BS upregulated the transcription of the UL122 upstream sequence. This could indicate the presence of an activating transcription factor site overlapping the p53 BS. Although the effects observed in these reporter assays were small, it should be noted that the two- to fourfold changes measured were consistent with the findings of previous studies of the effects of p53 on UL94 gene expression (56).
In future experiments, we will focus on generating recombinant viruses with mutated p53 BSs, since we are fully aware that the genomic context of these upstream sequences could be extremely important. We will also analyze the expression of the UL63 and UL65 proteins, among others, since mRNA expression and protein expression do not necessarily correlate.
In summary, this study demonstrates that p53 is a crucial player in driving viral gene expression at IE and E times p.i., and it presents an explanation for the delayed and reduced production of virions in p53–/– cells. It also contributes to further defining the role that p53 plays in the viral life cycle and expands on our earlier findings that an intact p53 DNA binding domain is essential for its sequestration and that p53 is directly bound to the viral genome at multiple times p.i. Contrary to our expectation that p53 would be implicated directly in viral transcription control, the current results indicate that the binding of p53 to the viral genome may be a mechanism the virus uses to sequester p53 in the viral replication centers, thus preventing p53 from activating its normal cellular targets.
We thank Matt Settles (Washington State University) for friendly assistance with programming and data processing and Dan Streblow (Oregon Health and Science University) for help with the initial design of experiments.
Published ahead of print on 18 February 2009. ![]()
|
|
|---|
by tumor suppressor p53: involvement of the ligand-binding domain and histone deacetylase activity. Mol. Endocrinol. 16:402-410.This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2010 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»