JVI Figure table search 04
Home Help [Feedback] [For Subscribers] [Archive] [Search] [Contents]
This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrowReprints and Permissions
Right arrow Copyright Information
Right arrow Books from ASM Press
Right arrow MicrobeWorld
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Chang, Y. E.
Right arrow Articles by Laimins, L. A.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Chang, Y. E.
Right arrow Articles by Laimins, L. A.

 Previous Article  |  Next Article 

Journal of Virology, May 2000, p. 4174-4182, Vol. 74, No. 9
0022-538X/00/$04.00+0
Copyright © 2000, American Society for Microbiology. All rights reserved.

Microarray Analysis Identifies Interferon-Inducible Genes and Stat-1 as Major Transcriptional Targets of Human Papillomavirus Type 31

Yijan E. Chang and Laimonis A. Laimins*

Department of Microbiology-Immunology, Northwestern University Medical School, Chicago, Illinois 60611

Received 9 November 1999/Accepted 2 February 2000


    ABSTRACT
Top
Abstract
Introduction
Materials and Methods
Results
Discussion
References

Human papillomaviruses (HPVs) infect keratinocytes and induce proliferative lesions. In infected cells, viral gene products alter the activities of cellular proteins, such as Rb and p53, resulting in altered cell cycle response. It is likely that HPV gene products also alter expression of cellular genes. In this study we used microarray analysis to examine the global changes in gene expression induced by high-risk HPV type 31 (HPV31). Among 7,075 known genes and ESTs (expressed sequence tags) tested, we found that 178 were upregulated and 150 were downregulated twofold or more in HPV31 cells compared to normal human keratinocytes. While no specific pattern could be deduced from the list of genes that were upregulated, downregulated genes could be classified to three groups: genes that are involved in the regulation of cell growth, genes that are specifically expressed in keratinocytes, and genes whose expression is increased in response to interferon stimulation. The basal level of expression of several interferon-responsive genes was found to be downregulated in HPV31 cells by both microarray analysis and Northern blot analysis in different HPV31 cell lines. When cells were treated with alpha or gamma interferon, expression of interferon-inducible genes was impaired. At high doses of interferon, the effects were less pronounced. Among the genes repressed by HPV31 was the signal transducer and activator of transcription (Stat-1), which plays a major role in mediating the interferon response. Suppression of Stat-1 expression may contribute to a suppressed response to interferon as well as immune evasion.


    INTRODUCTION
Top
Abstract
Introduction
Materials and Methods
Results
Discussion
References

Human papillomaviruses (HPVs) are small DNA tumor viruses which infect epithelial tissue and induce hyperproliferative lesions. Infection by high-risk genital HPVs (HPV type 16 [HPV16], HPV18, HPV31, HPV33, and HPV54) is associated with the development of anogenital cancers, while the low-risk HPVs (HPV6 and HPV11) more commonly induce benign genital warts (21, 29, 34, 66). Papillomavirus infection can be divided into several stages. Primary infection occurs at the basal layer of epidermis; following entry, HPV genomes are replicated in S phase as extrachromosomal elements. Genome copy number is maintained at a constant level in these cells, and a low level of transcripts is expressed from the early promoter. As infected keratinocytes divide, they stratify and differentiate. Concurrent with differentiation, HPV genomes are amplified and the late promoter initiates capsid gene expression (8, 13, 23). This results in encapsidation of viral DNA and production of infectious viruses.

Since HPVs encode only 8 to 10 proteins, they must rely extensively on cellular factors to regulate viral transcription and replication. The interaction of viral proteins with the host cell factors is therefore essential for the productive life cycle. Several papillomavirus proteins can directly activate transcription, while others act indirectly by altering the activity of cellular factors. Members of the former group include the E2 protein, which can activate transcription of reporters with multimerized upstream E2 binding sites in transient assays (4, 36). The latter group includes the two oncoproteins, E6 and E7. E6 binds a cellular ubiquitin ligase, which then targets the transactivator p53 for degradation (22, 49, 50, 62). E7 binds the Rb family of proteins, resulting in altered regulation of E2F-inducible genes which control S-phase entry (5, 9, 35, 41). It seems likely that during an HPV infection, alterations in cellular gene expression occur.

Many studies investigating the functions of viral proteins have not been done within the context of the entire viral genome and may not reflect what occurs at physiological concentration as in viral infection. In addition, cervical cancer cell lines that harbor integrated HPV genomes have been used in many studies, and these may not accurately reflect virus-cell interactions during the normal course of infection. Recently, a tissue culture model that simulates the latent stage of HPV infection has been developed. Cloned HPV31 and HPV18 genomes have been used to transfect normal human keratinocytes (NHKs), resulting in cell lines that exhibit characteristics of basal keratinocytes infected with HPV (12, 13, 40). Once induced to differentiate, these cell lines are able to activate late viral functions characteristic of a productive infection (13). This system has been used extensively for studying virus-cell interactions under physiologically relevant conditions.

In recent years, the development of microarrays, or gene chips, has provided a powerful tool to study complicated biological process which results in altered global gene expression (16, 25, 51, 65). Fluorescence-labeled cDNAs derived from two samples are hybridized to a microarray which contains thousands of oligonucleotides or cDNAs. Quantitative measurement of binding allows one to determine the changes in the transcriptional profile of many genes simultaneously. In this study, we used this methodology to examine the effect of high-risk HPV31 on cellular gene expression in keratinocytes. Transcriptional profiles in HPV31-positive cells were compared to those of matched NHKs, using microarray analysis. We found that HPV31 causes major reductions in the expression of genes involved in cell proliferation and interferon response, as well as genes specifically expressed in keratinocytes.


    MATERIALS AND METHODS
Top
Abstract
Introduction
Materials and Methods
Results
Discussion
References

Cell culture. NHKs were isolated from human foreskin tissue as described previously (48). Wild-type HPV31 genome was used to immortalize NHKs as described elsewhere (13). The established HPV lines were maintained in serum containing medium (E medium) supplemented with mouse epidermal growth factor (EGF; 5 ng/ml; Collaborative Biomedical Products, Bedford, Mass.) in the presence of murine 3T3 J2 fibroblast feeders (38, 39). Two independently derived HPV lines were used in these studies; the T31 line was used in the microarray analysis, and the LKP31 line (13) was used in Northern and Western blot analysis. NHKs were maintained in E medium in the presence of feeders and EGF for microarray, Northern blot, or Western blot analyses. The growth conditions for NHKs and HPV31-containing cells were identical for all assays. Three NHK lines were used in these studies; TP was used in microarray analysis, and 407 and JP were used in other analyses. All cell lines tested had been frozen and thawed. 3T3 J2 fibroblast feeder cells were maintained in Dulbecco modified Eagle medium supplemented with 10% calf serum (Gibco BRL, Grand Island, N.Y.). Alpha and gamma interferon (IFN-alpha and -gamma ) were purchased from BioSource (Camarillo, Calif.).

Microarray analysis. NHKs or HPV cells were grown in E medium supplemented with mouse EGF in the presence of fibroblast feeders and harvested at 80% confluence. Fibroblast feeders were first removed by versene (phosphate-buffered saline [Gibco BRL], with 0.5 mM EDTA) prior to RNA isolation. Total RNA was isolated by lysing cells in Trizol reagent (Gibco BRL), and poly(A) RNA was affinity purified on an Oligotex column (Qiagen, Valencia, Calif.). Generation of cDNA, fluorescent labeling, and hybridization to the gene chip were performed by Genome Systems (St. Louis, Mo.). Briefly, mRNA was isolated and reverse transcribed with 5' Cy3- or Cy5-labeled random 9-mer (Operon Technologies, Inc., Alameda, Calif.). The paired reactions were combined and purified with a TE-30 column (Clontech, Palo Alto, Calif.). Fluorescently labeled probe was then applied to the array for hybridization at 60°C for 6.5 h. After hybridization, the glass slides were washed with decreasing ionic strength and scanned with 10-µm resolution to detect Cy3 and Cy5 fluorescence. We examined 7,075 human genes and ESTs (expressed sequence tags) on the array Human UniGEM V (Incyte Pharmaceuticals, Inc., Palo Alto, Calif.) and analyzed the results with the software GEMtool 2.4. A gridding and region detection algorithm was used to determine each element. The area surrounding each element image was used to calculate a local background and was subtracted from the total element signal. Background-subtracted element signals were used to calculate Cy3/Cy5 ratios. The average of the resulting total Cy3 and Cy5 signal gives a ratio that was used to balance or normalize the signals.

Northern blot analysis. NHKs and LKP31 cells were cultured to subconfluence in E medium supplemented with EGF. Fibroblast feeders were removed, and total RNA was isolated as described above. Equal amounts of total RNA as determined by UV absorbance were separated by 1% agarose formaldehyde gel electrophoresis and transferred to a nylon membrane (Zeta Probe; Bio-Rad, Hercules, Calif.) in 10× SSC (1× SSC is 0.15 M NaCl plus 0.015 M sodium citrate). The membrane was cut into strips and hybridized with cDNA probes labeled with [alpha 32-P]dCTP (Amersham, Buckinghamshire, England) by High Prime (Boehringer Mannheim, Indianapolis, Ind.). Quantitative analysis was done with a PhosphorImager (Molecular Dynamics, Sunnyvale, Calif.).

Western blot analysis. NHKs or LKP31 cells were grown to 80% confluence in E medium supplemented with EGF, and the cells were lysed in lysis buffer (10 mM Tris HCl [pH 7.4], 150 mM NaCl, 1% deoxycholate, 1% Nonidet P-40, 0.1% sodium dodecyl sulfate [SDS], 5 mM EDTA, 1% Triton X-100, 1 µg of leupeptin per ml, 1 mM phenylmethylsulfonyl fluoride, 0.5 mM sodium orthovanadate, 20 mM p-nitrophenyl phosphate) on ice for 10 min. Fibroblast feeders were removed as described above before harvesting of keratinocyte lysates. The clear lysate was obtained by centrifugation to remove insoluble cell debris at 12,000 × g for 5 min at 4°C. Equal amounts of protein lysates were separated by SDS-10% polyacrylamide gel electrophoresis (PAGE) and transferred to a polyvinylidene difluoride membrane (Immobilon-P; Millipore, Bedford, Mass.). The membrane was then blocked with 5% nonfat dry milk in wash buffer (phosphate-buffered saline, 0.1% Tween 20) for 1 h and incubated with primary antibody (mouse anti-Stat-1alpha /beta ; 1 µg/ml; Santa Cruz Biotechnology, Inc., Santa Cruz, Calif.) for another hour. The membrane was washed for 1 h with four changes of wash buffer and then incubated with a sheep anti-mouse secondary antibody conjugated with horseradish peroxidase (Amersham). After a second wash as described above, the signal was visualized by enhanced chemiluminescence (Amersham). The levels of Stat-1 protein were determined by densitometric analysis.


    RESULTS
Top
Abstract
Introduction
Materials and Methods
Results
Discussion
References

Microarray analysis of gene expression in normal and HPV31-transfected keratinocytes. To identify genes whose expression is altered by HPV31 during conditions similar to infection of basal cells, we examined a matched set of normal and HPV31-transfected human keratinocytes derived from a common donor. Efficient infection of tissue culture cells by HPV is not feasible, and so we examined transfected keratinocytes that stably maintain viral DNA as episomes (13). Human foreskin keratinocytes, derived from a single anonymous donor (designated TP), were transfected with recircularized cloned HPV31 along with a neomycin selectable marker. Following G418 selection, resistant colonies were pooled and expanded as combined culture (designated T31) to ensure that any observed effects were not specific to a single clonal cell line. The episomal state of viral DNA in the T31 pooled culture was confirmed by Southern analysis, and the expression of HPV31 early and late genes was found to be similar to that previously described for HPV31-positive cell lines (data not shown) (13, 48).

To identify genes whose expression was altered by HPV31 gene products, we used microarray analysis of 7,075 sequence-verified ESTs which included about 4,000 known genes (Genome Systems). Poly(A)-selected RNA was isolated from T31 and TP cells, and fluorescent dye-tagged cDNAs (Cy3 and Cy5, respectively) were generated. These fluorescent cDNAs were simultaneously hybridized to probe sequences on the microarray, and the amount of fluorescence seen with the individual dyes was determined by confocal microscopy. The differential expression of each EST was then calculated from the relative intensity of the Cy3-versus-Cy5 fluorescent signal. The plot of the differential expression of the 7,075 genes is shown in Fig. 1. Overall, most ESTs were not altered in expression by HPV31. Approximately 2.5% (i.e., 178) of ESTs examined were increased in expression between 2- and 3-fold, with no gene activated more than 3.2-fold. In the same analysis, 150 genes were found to be repressed at least twofold by HPV31 gene products. Twenty genes were repressed between 11.9- and 3.2-fold, while 130 genes were repressed 3.2- to 2.0-fold. Genes whose expression was altered more than 2.3-fold are listed in Table 1. The complete list of genes examined and their level of induction or suppression will be available at the web site http://bugs.mimnet.nwu.edu/laimins_lab/chang.


View larger version (57K):
[in this window]
[in a new window]
  		FIG. 1.   Global comparison of gene expression in HPV31 and NHKs. Each dot corresponds to the Cy3 (x axis) and Cy5 (y axis) fluorescent intensity of one single element on the microarray. Twofold, 5-fold, and 10-fold changes in expression are indicated as parallel lines. Dots that represent more than 11.5-fold changes are internal controls for hybridization.

                              
View this table:
[in this window]
[in a new window]
  		TABLE 1.   Genes whose expression was altered more than 2.3-fold

The majority of genes whose expression was increased by more than twofold by HPV31 are ESTs of unknown identity (Table 1). Several known genes such as those encoding ubiquitin carboxyl-terminal esterase, cathepsin H, adenomatosis polyposis coli, and corticotropin-releasing hormone were increased approximately threefold in expression. Overall, no obvious patterns of activation of specific gene families by HPV31 could be identified. Interestingly, no E2F-regulated gene examined was found to be activated more than twofold. For example, cyclin E was found to be altered slightly (1.1- to 1.3-fold) in HPV 31 cells.

Expression of interferon-inducible genes is repressed by HPV31 gene products. In contrast to the upregulated genes, those genes whose expression was repressed by HPV proteins could be grouped into three categories: regulators of cell growth, keratinocyte-specific genes, and factors mediating the interferon response. The p21 gene, whose expression is regulated by p53, was found to be reduced 2.6-fold, consistent with the increased turnover of p53 induced by high-risk E6 (10, 47, 50). Two other genes, encoding Mad and transgelin, were also downregulated in HPV cells. Mad is a negative regulator of Myc and Max activities (2, 24), while transgelin is highly expressed in senescent cells (56). Several genes encoding keratinocyte-specific proteins such as SPRK, human small proline-rich protein II (SprII), defensin, desmocollin 2, desmoplakin, and stratifin, were also downregulated by HPV gene products (Table 1) (26, 30, 37, 42, 58, 64).

The most pronounced reductions in gene expression were observed for interferon-responsive genes. At least 10 of the 67 genes downregulated more than 2.3-fold were interferon inducible and included those with the most significant decreases. The interferon-responsive genes whose basal levels of expression were repressed by HPV31 included those encoding myxovirus (influenza virus) resistance 1, IFN-alpha -inducible 11.5-kDa protein, interferon-inducible protein 56, 2',5'-oligoadenylate synthetase 2, Stat-1, interferon-inducible protein 1-8U, Staf50, myxovirus resistance 2, interferon-inducible leucine zipper protein IFP35, interferon-inducible guanylate binding protein 1, and interferon-inducible protein 54 (1, 18, 19, 33, 44, 45, 52, 57, 60, 61).

We next wanted to confirm that genes whose expression was reduced by HPV31 by microarray analysis were also reduced when examined by Northern analysis. In addition, it was important to ensure that the changes observed were not specific to the matched set of keratinocytes examined by microarray. Equal amounts of total RNAs from another HPV31-positive cell line (LKP31) as well as RNA from a second NHK donor (designated 406) were isolated and examined by Northern analysis. The results for Northern analysis of seven of the repressed genes are shown in Fig. 2, and the levels of reduction were consistent with those observed in the microarray analysis. Subsequently, we examined the expression levels of MxA and Stat-1 in three other transfected keratinocyte lines containing episomal HPV31 DNA, as well as two more NHK lines, and obtained similar results (data not shown). In our studies, the repression of interferon-inducible genes was observed in five different HPV31 lines and five different NHK lines with passage number ranging from 4 to 13 (HPV31 cells) and from 2 to 5 (NHKs).


View larger version (56K):
[in this window]
[in a new window]
  		FIG. 2.   Level of expression in genes suppressed in HPV31 cells demonstrated by microarray analysis is confirmed by Northern blot analysis. Total RNAs (20 µg) from NHKs or LKP31 cells were separated by formaldehyde agarose electrophoresis, transferred to a membrane, and hybridized with various cDNA probes. The slightly faster mobility of p21 and Stat-1 signals in LKP31 cells is likely due to electrophoretic artifact since it was not observed in separate experiments.

Repression of interferon-inducible genes by HPV31 results in an altered interferon response. One of the major changes in gene expression identified by our microarray analysis was the repression by HPV31 of a large number of interferon-inducible genes compared to NHKs. This repression was observed in the absence of any interferon treatment. We wanted to determine if this basal repression would result in an impaired response of HPV31-positive cells to interferons. This was first examined by treating LKP31 cells and NHKs with IFN-alpha and determining the expression of a representative interferon-inducible gene, MxA, which encodes the interferon-responsive myxovirus resistance gene 1. MxA is primarily induced by IFN-alpha (1), and its basal level of expression was lower in LKP31 cells than in NHKs (Fig. 3A, lanes 1 and 6). Subconfluent cultures of LKP31 cells and NHKs were treated with IFN-alpha at one of two different concentrations (100 or 500 U/ml), and MxA expression was analyzed by Northern analysis. As seen in Fig. 3A, MxA RNA was induced upon interferon treatment in NHKs as well as in LKP31 cells, but the level of induction was reduced in LKP31 cells, especially with the lower dose of IFN-alpha tested (100 U/ml [Fig. 3B]). Upon treatment with higher levels (1,000 U/ml) of IFN-alpha , the induction of MxA RNA in LKP31 cells was also reduced, especially at earlier time points (Fig. 3C, lanes 1, 2, 6, and 7). However, after the initial delay, the amounts of MxA RNA in LKP31 cells and NHK reached comparable levels at 7 h posttreatment (Fig. 3D).


View larger version (22K):
[in this window]
[in a new window]
  		FIG. 3.   Induction of MxA RNA by IFN-alpha is impaired in HPV31 cells. (A) Northern blot analysis of 8 µg of total RNA from NHK or LKP31 cells treated with IFN-alpha (100 or 500 U/ml) with MxA cDNA as a probe. MxA RNA signals in panel A were quantified by a PhosphorImager, and the results are shown graphically in panel B. The induction of MxA RNA by IFN-alpha reaches comparable levels in LKP31 cells and NHKs at a higher dose (1,000 U/ml) and longer time exposure. (C) Northern blot analysis of 10 µg of total RNA from NHKs and LKP31 cells treated with IFN-alpha with MxA cDNA as a probe. MxA RNA signals in panel C were quantified by a PhosphorImager, and the results are shown graphically in panel D.

Among the interferon-inducible genes whose expression was downregulated by HPV31 was Stat-1. Stat-1 is a major regulator of interferon-stimulated transactivation, and the downregulation of other interferon-inducible genes could be the result of reduced Stat-1 expression (7, 20). In addition, the expression of Stat-1 itself is activated by interferon treatment. We therefore examined Stat-1 protein levels in response to interferon treatment. As seen in Fig. 4A, the levels of Stat-1 protein were significantly reduced in untreated LKP31 cells (lanes 1 and 6). Following treatment with IFN-alpha (100 and 500 U/ml), there was an induction in Stat-1 protein but levels remained reduced in LKP31 cells and did not reach levels seen in NHKs (Fig. 4B). Upon treatment with higher doses of IFN-alpha (1,000 U/ml), the response was still reduced until 20 h posttreatment, when the levels of Stat-1 reached those seen in NHKs (Fig. 4C, lanes 3, 4, 7, and 8; Fig. 4D). These results suggest that the response to IFN-alpha was reduced in HPV31-positive cells. At higher doses and at longer times of exposure to IFN-alpha , the effects were less pronounced.


View larger version (20K):
[in this window]
[in a new window]
  		FIG. 4.   The induction of Stat-1 protein by IFN-alpha is impaired in HPV31 cells. (A) Western blot analysis of equal amounts of protein lysates from NHKs and LKP31 cells treated with IFN-alpha (100 or 500 U/ml) to detect Stat-1 protein level. Stat-1 protein migrates in SDS-PAGE as a doublet with molecular weights of 91,000 and 84,000. Quantitative estimation of the results in panel A was done with a densitometer and shown graphically in panel B. The induction of Stat-1 protein by IFN-alpha reaches comparable levels in LKP31 cells and NHKs with a higher dose (1,000 U/ml) and longer time. (C) Western blot analysis of equal amount of protein lysates from NHKs or LKP31 cells treated with IFN-alpha to detect Stat-1 protein. Quantitative estimation of the results in panel C was done with a densitometer and shown graphically in panel D.

Stat-1 plays a major role in both IFN-alpha /beta and IFN-gamma pathways. Upon exposure to IFN-alpha , Stat-1 is phosphorylated by Tyk2 and Jak1 in the cytoplasm and forms heterodimer with phosphorylated Stat-2. The heterodimer translocates to the nucleus and forms the ISGF-3 (interferon-stimulated gene factor 3) complex with p48, leading to transactivation of IFN-alpha /beta -responsive genes. When cells are treated with IFN-gamma , Stat-1 is phosphorylated by Jak1 and Jak2, forms a homodimer, and translocates to the nucleus to activate gene expression. Since the basal level of Stat-1 is downregulated in HPV31 cells, we next examined if the IFN-gamma response was also altered by HPV31. For this, we examined the induction of Stat-1 protein by IFN-gamma . As shown in Fig. 5, the levels of Stat-1 protein were low in LKP31 cells, and while they increased upon IFN-gamma treatment, they remained significantly below the levels seen in NHKs at earlier time points (Fig. 5A, lanes 1 to 3 and 6 to 8). However, at 18 and 24 h after IFN-gamma treatment, the amount of Stat-1 protein reached a level comparable to that seen in NHKs (Fig. 5A, lanes 4, 5, 9, and 10; Fig. 5B). We conclude that in response to IFN-alpha and -gamma stimulation, HPV31 cells exhibit a lag in activation compared to NHKs. Eventually this lag in response could be restored upon a longer exposure to interferon or at a higher dose of interferon. The interferon response is complex, involving multiple factors. The initial delay observed in HPV31 cells maybe explained, in part, by a reduced basal level of Stat-1.


View larger version (24K):
[in this window]
[in a new window]
  		FIG. 5.   The induction of Stat-1 protein by IFN-gamma is impaired in HPV31 cells. (A) Western blot analysis of equal amounts of protein lysates from NHKs and LKP31 cells treated with IFN-gamma (1,000 U/ml) to detect Stat-1 protein level. The image in panel A was quantified by a densitometer, and the results are shown graphically in panel B.


    DISCUSSION
Top
Abstract
Introduction
Materials and Methods
Results
Discussion
References

In this report, we used microarray analysis to examine the effect of HPV31 genome on the transcription of cellular genes. In cells that stably maintain HPV31 episomes, we observed the expression of 178 ESTs to be upregulated and that of 150 to be downregulated at least twofold among the 7,075 verified sequences screened. These cells are believed to mimic infected basal keratinocytes in vivo, as they can induce late functions upon differentiation. In our analysis, no gene was found to be activated by HPV31 gene products more than 3.2-fold, and among these no easily discernible families could be identified. Since the high-risk E7 proteins bind Rb, it was expected that E2F-responsive genes would be substantially increased in expression in HPV31-positive cells. However, none were found to be activated more than twofold. This could be due to the asynchronous nature of the cells examined or other, more complex reasons.

Among the 150 genes that were repressed by HPV31, three groups were identified. The first group included genes involved in regulation of cell growth. The E6 protein of HPV31, like that of HPV16 and -18, increases the turnover rate of p53, resulting in decreased steady-state levels (55). A major transcriptional target of p53 is the p21waf1 gene, which is a negative regulator of cyclin-dependent kinases (15). Consistent with the reduction in p53 levels, the levels of p21 transcripts were found to be reduced by both microarray and Northern analysis. In addition, Mad, the cellular regulator of Myc activity, was found to be reduced in expression by HPV31. Mad proteins form a heterodimeric complex with Max and antagonize the positive effect of Myc-Max heterodimers (2). Myc has been implicated as an activator of telomerase, and in cells expressing HPV E6 and E7, high levels of telomerase activity have been observed (27, 28, 59, 63). In our studies, no consistent transcriptional activation of Myc was detected, and downregulation of Mad expression may provide an alternative mechanism by which Myc function can be increased in infections by high-risk HPV types. A third gene found to be repressed by HPV31 is the transgelin gene. Transgelin is able to cross-link actin and remodel cytoskeleton (43). The expression of transgelin is often increased in senescing cells and is overexpressed in prematurely aging patients with Werner syndrome (53, 56). Downregulation of transgelin may contribute to loss of contact inhibition and enhanced motility characteristic of HPV-transformed cells (43).

The second class of genes found to be repressed by HPV31 were those which are expressed specifically in keratinocytes. These included genes encoding SprII, a small proline-rich protein found in UV-irradiated keratinocytes (26), and defensin, which is highly expressed in airway epithelia and has been suggested to possess antimicrobial activity (37, 64). Another repressed gene encodes stratafin, a member of family of acidic proteins that are enriched in stratifying epithelia (30). The most interesting of the keratinocyte-specific genes found to be repressed were the desmocollin and desmoplakin genes. It is possible that decreased expression of these factors could lead to less stringent cell-cell junctions that may be advantageous to productive viral life cycle.

The most intriguing set of genes found to be downregulated by HPV31 gene products were the interferon-responsive genes. Among these, Stat-1 is a primary regulator of the interferon response pathway. Stat-1 belongs to a family of proteins which normally remain latent in cytoplasm. Stat-1 is phosphorylated and activated upon interferon stimulation and turns on genes under the control of the ISRE element (7, 20) and the GAS element (7). The Stat-2 gene, also examined on the microarray, showed only marginal reduction in HPV31 cells (-1.1-fold [data not shown]). The low basal level of Stat-1 in HPV31 cells could be responsible for the low basal level of various interferon-inducible genes and contribute to the impaired response to interferon stimulation.

Additional mechanisms for altering the cellular response to interferon by HPV gene products have been previously reported. The HPV16 E6 protein has been shown to bind to and inhibit the function of the IRF3 protein (46). The IRF3 protein is an activator of interferon synthesis, and disruption of this activity interferes with the interferon response. In addition, the HPV16 E7 protein has been shown to directly bind to the p48 subunit of the ISGF-3 complex and prevent its nuclear translocation (3). Many viruses target and debilitate the interferon and immune responses during infection. Examples include adenoviruses, herpes simplex viruses (HSV), and hepatitis C viruses (6, 14, 54), all of which inhibit the activity of the double-stranded RNA-dependent kinase, PKR. In addition, adenovirus E1A binds to p48 and inhibits its nuclear translocation (31, 32). With respect to immune surveillance, HSV can downregulate antigen presentation by interfering with peptide translocation via the TAP proteins (17). Thus, it is not surprising that HPV uses multiple mechanisms to interfere with repression of the interferon response. It is interesting that among the 7,075 ESTs examined in this study, the interferon-inducible genes were identified as major targets of HPV action. This underscores the primary importance of suppressing the interferon response during papillomavirus infection.

The repression of the interferon response by HPV31 may contribute to evasion of HPV-infected cells from immune surveillance by the host. Interferon also stimulates expression of class I and II major histocompatibility antigens as well as LMP2 and -7, which are part of the 20S proteasome degradation machinery used to display of antigens on the cell surface. In our HPV-positive cells, we did not detect increased expression of class I or class II genes, suggesting that viral products may impair the activation in response to viral infection (data not shown). It is therefore possible that repression of Stat-1 expression may contribute to the evasion of the immune response during HPV infections. In our studies, the basal levels of Stat-1 were found to be reduced consistent with a repression of transcription initiation. However, it is equally plausible that this reduction is due to increased turnover of Stat-1 mRNA. It will be important to determine which of the HPV gene products are responsible for the repression of Stat-1, the mechanism of action, and whether this property extends to low-risk HPV types.

Previous studies using microarray analysis of a number of ESTs similar to that used in our study demonstrated that infection by human cytomegalovirus (HCMV) altered expression of 248 genes more than 4-fold, with many activated 6- to 15-fold (65). Interestingly, the expression of most interferon-responsive genes was increased upon HCMV infection. Acute viral infection is known to trigger an interferon response, and it is not surprising that cell lines that stably maintain the HPV genome showed an altered response manifested in this case by suppressed basal levels of interferon-inducible genes. This suppression may contribute to evasion of immune surveillance. In our study, no gene was found to be activated more than 3.2-fold, while 22 genes were repressed between 11- and 3-fold. The less severe effects seen in our studies may reflect differences between acute infection by HCMV and latent infection that is modeled by our stable HPV cell lines. Alternatively, it may be the result of more fundamental differences in the biology of these different viruses. Microarray analysis of genes activated by the platelet-derived growth factor receptor following ligand binding revealed no gene to be activated more than twofold (11). It is possible that the changes in expression of less than twofold by HPV31 can lead to profound changes in the cellular environment; however, we did not investigate these genes in detail in this study. This study has identified possible pathways by which HPV modifies the cellular environment to facilitate viral replication and evade detection by the immune system.


    ACKNOWLEDGMENTS

We thank G. Sen for valuable advice and reagents, and we thank K. Rundell, R. Longnecker, F. Stubenrauch, and W. Hubert for comments on the manuscript. We also thank R. N. Eisenman and C. M. Horvath for Mad and Stat-1 cDNAs and C. Biaggi, W. Hubert, and A. Merchant for the HPV31 immortalized line T31 and for NHK isolates 407 and TP.

This work was supported by grants from the National Cancer Institute (CA59655) to L.A.L. and the Illinois Department of Public Health (96190152) to Y.E.C.


    FOOTNOTES

* Corresponding author. Mailing address: Department of Microbiology-Immunology, Northwestern University Medical School, 303 E. Chicago Ave., Chicago, IL 60611. Phone: (312) 503-0648. Fax: (312) 503-1339. E-mail: lal{at}merle.acns.nwu.edu.


    REFERENCES
Top
Abstract
Introduction
Materials and Methods
Results
Discussion
References

1. Aebi, M., J. Fah, N. Hurt, C. E. Samuel, D. Thomis, L. Bazzigher, J. Pavlovic, O. Haller, and P. Staeheli. 1989. cDNA structures and regulation of two interferon-induced human Mx proteins. Mol. Cell. Biol. 9:5062-5072[Abstract/Free Full Text].
2. Ayer, D. E., L. Kretzner, and R. N. Eisenman. 1993. Mad: a heterodimeric partner for Max that antagonizes Myc transcriptional activity. Cell 72:211-222[CrossRef][Medline].
3. Barnard, P., and N. A. J. McMillan. 1999. The human papillomavirus E7 oncoprotein abrogates signaling mediated by interferon-alpha . Virology 259:305-313[CrossRef][Medline].
4. Bouvard, V., A. Storey, D. Pim, and L. Banks. 1994. Characterization of the human papillomavirus E2 protein: evidence of trans-activation and trans-repression in cervical keratinocytes. EMBO J. 13:5451-5459[Medline].
5. Cheng, S., G. D. Schmidt, T. Murant, T. R. Broker, and L. T. Chow. 1995. Differentiation-dependent up-regulation of the human papillomavirus E7 gene reactivates cellular DNA replication in suprabasal differentiated keratinocytes. Genes Dev. 9:2335-2349[Abstract/Free Full Text].
6. Chou, J., J. J. Chen, M. Gross, and B. Roizman. 1995. Association of a M(r) 90,000 phosphoprotein with protein kinase PKR in cells exhibiting enhanced phosphorylation of translation initiation factor eIF-2 alpha and premature shutoff of protein synthesis after infection with gamma 1 34.5- mutants of herpes simplex virus 1. Proc. Natl. Acad. Sci. USA 92:10516-10520[Abstract/Free Full Text].
7. Darnell, J. E., Jr., I. M. Kerr, and G. R. Stark. 1994. Jak-STAT pathways and transcriptional activation in response to IFNs and other extracellular signaling proteins. Science 264:1415-1421[Abstract/Free Full Text].
8. Dollard, S. C., J. L. Wilson, L. M. Demeter, W. Bonnez, R. C. Reichman, T. R. Broker, and L. T. Chow. 1992. Production of human papillomavirus and modulation of the infectious program in epithelial raft cultures. Genes Dev. 6:1131-1142[Abstract/Free Full Text].
9. Dyson, N., P. M. Howley, K. Munger, and E. Harlow. 1989. The human papillomavirus-16 E7 oncoprotein is able to bind to the retinoblastoma gene product. Science 243:934-937[Abstract/Free Full Text].
10. El-Deiry, W., T. Tokino, V. Velculescu, D. Levy, R. Parson, J. Trent, D. Lin, W. Mercer, K. Kinzler, and B. Vogelstein. 1993. WAF1, a potential mediator of p53 tumor suppression. Cell 75:817-825[CrossRef][Medline].
11. Fambrough, D., K. McClure, A. Kazlauskas, and E. S. Lander. 1999. Diverse signaling pathways activated by growth factor receptors induce broadly overlapping, rather than independent, sets of genes. Cell 97:727-741[CrossRef][Medline].
12. Frattini, M. G., H. B. Lim, J. Doorbar, and L. A. Laimins. 1997. Induction of human papillomavirus type 18 late gene expression and genomic amplification in organotypic cultures from transfected DNA templates. J. Virol. 71:7068-7072[Abstract].
13. Frattini, M. G., H. B. Lim, and L. A. Laimins. 1996. In vitro synthesis of oncogenic human papillomaviruses requires episomal genomes for differentiation-dependent late expression. Proc. Natl. Acad. Sci. USA 93:3062-3067[Abstract/Free Full Text].
14. Ghadge, G. D., P. Malhotra, M. R. Furtado, R. Dhar, and B. Thimmapaya. 1994. In vitro analysis of virus-associated RNA I (VAI RNA): inhibition of the double-stranded RNA-activated protein kinase PKR by VAI RNA mutants correlates with the in vivo phenotype and the structural integrity of the central domain. J. Virol. 68:4137-4151[Abstract/Free Full Text].
15. Harper, J., G. Adami, N. Wei, K. Keyomarsi, and S. Elledge. 1993. The p21 cdk-interacting protein Cip1 is a potent inhibitor of G1 cyclin-dependent kinases. Cell 75:805-816[CrossRef][Medline].
16. Heller, R. A., M. Schena, A. Chai, D. Shalon, T. Bedilion, J. Gilmore, D. E. Woolley, and R. W. Davis. 1997. Discovery and analysis of inflammatory disease-related genes using cDNA microarrays. Proc. Natl. Acad. Sci. USA 94:2150-2155[Abstract/Free Full Text].
17. Hill, A., P. Jugovic, I. York, G. Russ, J. Bennink, J. Yewdell, H. Ploegh, and D. Johnson. 1995. Herpes simplex virus turns off the TAP to evade immunity. Nature 375:411-415[CrossRef][Medline].
18. Horisberger, M. A. 1992. Interferon-induced human protein MxA is a GTPase which binds transiently to cellular proteins. J. Virol. 66:4705-4709[Abstract/Free Full Text].
19. Horisberger, M. A., G. K. McMaster, H. Zeller, M. G. Wathelet, J. Dellis, and J. Content. 1990. Cloning and sequence analyses of cDNAs for interferon- and virus-induced human Mx proteins reveal that they contain putative guanine nucleotide sites: functional study of the corresponding gene promoter. J. Virol. 64:1171-1181[Abstract/Free Full Text].
20. Horvath, C. M., and J. E. Darnell, Jr. 1997. The state of the STATs: recent developments in the study of signal transduction to the nucleus. Curr. Opin. Cell Biol. 9:233-239[CrossRef][Medline].
21. Howley, P. M. 1996. Papillomavirinae: the viruses and their replication., p. 947-978. In B. N. Fields, D. M. Knipe, and P. M. Howley (ed.), Fundamental virology, 3rd ed. Lippincott-Raven Publisher, Philadelphia, Pa.
22. Huibregtse, J. M., M. Scheffner, and P. M. Howley. 1991. A cellular protein mediates association of p53 with the E6 oncoprotein of human papillomavirus types 16 or 18. EMBO J. 10:4129-4135[Medline].
23. Hummel, M., J. B. Hudson, and L. A. Laimins. 1992. Differentiation-induced and constitutive transcription of human papillomavirus type 31b in cell lines containing viral episomes. J. Virol. 66:6070-6080[Abstract/Free Full Text].
24. Hurlin, P. J., K. P. Foley, D. E. Ayer, R. N. Eisenman, D. Hanahan, and J. M. Arbeit. 1995. Regulation of Myc and Mad during epidermal differentiation and HPV-associated tumorigenesis. Oncogene 11:2487-2501[Medline].
25. Iyer, V. R., M. B. Eisen, D. T. Ross, G. Schuler, T. Moore, J. C. F. Lee, J. M. Trent, L. M. Staudt, J. J. Hudson, M. S. Boguski, D. Lashkari, D. Shalon, D. Botstein, and P. O. Brown. 1999. The transcriptional program in the response of human fibroblasts to serum. Science 283:83-87[Abstract/Free Full Text].
26. Kartasova, T., and P. van de Putte. 1988. Isolation, characterization, and UV-stimulated expression of two families of genes encoding polypeptides of related structure in human epidermal keratinocytes. Mol. Cell. Biol. 8:2195-2203[Abstract/Free Full Text].
27. Kiyono, T., S. A. Foster, J. I. Koop, M. J. K., D. A. Galloway, and A. J. Klingelhutz. 1998. Both Rb/p16ink4a inactivation and telomerase activity are required to immortalize human epithelial cells. Nature 396:84-88[CrossRef][Medline].
28. Klingelhutz, A. J., S. A. Foster, and J. K. McDougall. 1996. Telomerase activation by the E6 gene product of human papillomavirus type 16. Nature 380:79-82[CrossRef][Medline].
29. Laimins, L. A. 1993. The biology of human papillomaviruses: from warts to cancer. Infect. Agents Dis. 2:74-86[Medline].
30. Leffers, H., P. Madsen, H. H. Rasmussen, B. Honore, A. H. Andersen, E. Walbum, J. Vandekerckhove, and J. E. Celis. 1993. Molecular cloning and expression of the transformation sensitive epithelial marker stratifin. J. Mol. Biol. 231:982-998[CrossRef][Medline].
31. Leonard, G. T., and G. Sen. 1997. Restoration of interferon responses of adenovirus E1A-expressing HT1080 cell lines by overexpression of p48 protein. J. Virol. 71:5095-5101[Abstract].
32. Leonard, G. T., and G. C. Sen. 1996. Effects of adenovirus E1A protein on interferon-signaling. Virology 224:25-33[CrossRef][Medline].
33. Levy, D., A. Larner, A. Chaudhuri, L. Babiss, and J. E. Darnell, Jr. 1986. Interferon-stimulated transcription: isolation of an inducible gene and identification of its regulatory region. Proc. Natl. Acad. Sci. USA 83:8929-8933[Abstract/Free Full Text].
34. Lowy, D. R., R. Kirnbauer, and J. T. Schiller. 1994. Genital human papillomavirus infection. Proc. Natl. Acad. Sci. USA 91:2436-2440[Abstract/Free Full Text]. (Review.)
35. Martin, L. G., G. W. Demers, and D. A. Galloway. 1998. Disruption of the G1/S transition in human papillomavirus type 16 E7-expressing human cells is associated with altered regulation of cyclin E. J. Virol. 72:975-985[Abstract/Free Full Text].
36. McBride, A. A., H. Romanczuk, and P. M. Howley. 1991. The papillomavirus E2 regulatory proteins. J. Biol. Chem. 266:18411-18414[Free Full Text].
37. McCray, P. B. J., and L. Bentley. 1997. Human airway epithelia express a beta -defensin. Am. J. Respir. Cell Mol. Biol. 16:343-349[Abstract].
38. Meyers, C., M. G. Frattini, J. B. Hudson, and L. A. Laimins. 1992. Biosynthesis of human papillomavirus from a continuous cell line upon epithelial differentiation. Science 257:971-973[Abstract/Free Full Text].
39. Meyers, C., and L. A. Laimins. 1994. In vitro systems for the study and propagation of human papillomaviruses. Curr. Top. Microbiol. Immol. 186:199-215. (Review.)
40. Meyers, C., T. J. Mayer, and M. A. Ozbun. 1997. Synthesis of infectious human papillomavirus type 18 in differentiating epithelium transfected with viral DNA. J. Virol. 71:7381-7386[Abstract].
41. Munger, K., B. A. Werness, N. Dyson, W. C. Phelps, E. Harlow, and P. M. Howley. 1989. Complex formation of human papillomavirus E7 proteins with the retinoblastoma tumor suppressor gene product. EMBO J. 8:4099-4105[Medline].
42. Parker, A. E., G. N. Wheeler, J. Arnemann, S. C. Pidsley, P. Ataliotis, C. L. Thomas, D. A. Rees, A. I. Magee, and R. S. Buxton. 1991. Desmosomal glycoproteins II and III. J. Biol. Chem. 266:10438-10445[Abstract/Free Full Text].
43. Prinjha, R. K., C. E. Shapland, J. J. Hsuan, N. F. Totty, I. J. Mason, and D. Lawson. 1994. Cloning and sequencing of cDNAs encoding the actin cross-linking protein transgelin defines a new family of actin-associated proteins. Cell Motil. Cytoskel. 28:243-255[CrossRef][Medline].
44. Rasmussen, U. B., C. Wolf, M.-G. Mattei, M.-P. Chenard, J.-P. Bellocq, P. Chambon, M.-C. Rio, and P. Basset. 1993. Identification of a new interferon-alpha -inducible gene (p27) on human chromosome 14q32 and its expression in breast carcinoma. Cancer Res. 53:4096-4101[Abstract/Free Full Text].
45. Robinson, P. A., J. J. Marley, A. S. High, and W. J. Hume. 1994. Differential expression of protease inhibitor and small proline-rich protein genes between normal human oral tissue and odontogenic keratinocytes. Arch. Oral Biol. 39:251-259[CrossRef][Medline].
46. Ronco, L. V., A. Y. Karpova, M. Vidal, and P. M. Howley. 1998. Human papillomavirus 16 E6 oncoprotein binds to interferon regulatory factor-3 and inhibits its transcriptional activity. Genes Dev. 12:2061-2072[Abstract/Free Full Text].
47. Ruesch, M. N., and L. A. Laimins. 1998. Human papillomavirus oncoproteins alter differentiation-dependent cell cycle exit on suspension in semisolid medium. Virology 250:19-29[CrossRef][Medline].
48. Ruesch, M. N., F. Stubenrauch, and L. A. Laimins. 1998. Activation of papillomavirus late gene transcription and genome amplification upon differentiation in semisolid medium is coincident with expression of involucrin and transglutaminase but not keratin-10. J. Virol. 72:5016-5024[Abstract/Free Full Text].
49. Scheffner, M., J. M. Huibregtse, R. D. Vierstra, and P. M. Howley. 1993. The HPV-16 E6 and E6-AP complex functions as a ubiquitin-protein ligase in the ubiquitination of p53. Cell 75:495-505[CrossRef][Medline].
50. Scheffner, M., B. A. Werness, J. M. Huibregtse, A. J. Levine, and P. M. Howley. 1990. The E6 oncoprotein encoded by human papillomavirus types 16 and 18 promotes the degradation of p53. Cell 63:1129-1136[CrossRef][Medline].
51. Schena, M., D. Shalon, R. Heller, A. Chai, P. O. Brown, and R. W. Davis. 1996. Parallel human genome analysis: microarray-based expression monitoring of 1000 genes. Proc. Natl. Acad. Sci. USA 93:10614-10619[Abstract/Free Full Text].
52. Schindler, C., X.-Y. Fu, T. Improta, R. Aebersold, and J. E. Darnell, Jr. 1992. Proteins of transcription factor ISGF-3: one gene encodes the 91- and 84-kDa ISGF-3 proteins that are activated by interferon alpha . Proc. Natl. Acad. Sci. USA 89:7836-7839[Abstract/Free Full Text].
53. Shapland, C., J. J. Hsuan, N. F. Totty, and D. Lawson. 1993. Purification and properties of transgelin: a transformation and shape change sensitive actin-gelling protein. J. Cell Biol. 121:1065-1073[Abstract/Free Full Text].
54. Taylor, D. R., S. T. Shi, P. R. Romano, G. N. Barber, and M. M. C. Lai. 1999. Inhibition of the interferon-inducible protein kinase PKR by HCV E2 protein. Science 285:107-110[Abstract/Free Full Text].
55. Thomas, J. T., and L. A. Laimins. 1998. Human papillomavirus oncoproteins E6 and E7 independently abrogate the mitotic spindle checkpoint. J. Virol. 72:1131-1137[Abstract/Free Full Text].
56. Thweatt, R., C. K. J. Lumpkin, and S. Goldstein. 1992. A novel gene encoding a smooth muscle protein is overexpressed in senescent human fibroblasts. Biochem. Biophys. Res. Commun. 187:1-7[CrossRef][Medline].
57. Tissot, C., and N. Mechti. 1995. Molecular cloning of a new interferon-induced factor that represses human immunodeficiency virus type 1 long terminal repeat expression. J. Biol. Chem. 270:14891-14898[Abstract/Free Full Text].
58. Virata, M. L. A., R. M. Wagner, D. A. D. Parry, and K. J. Green. 1992. Molecular structure of the human desmoplakin I and II amino terminus. Proc. Natl. Acad. Sci. USA 89:544-548[Abstract/Free Full Text].
59. Wang, J., L. Y. Xie, S. Allan, D. Beach, and G. J. Hannon. 1998. Myc activates telomerase. Genes Dev. 12:1769-1774[Abstract/Free Full Text].
60. Wang, X., L. M. Johansen, H.-J. Tae, and E. J. Taparowsky. 1996. IFP 35 forms complexes with B-ATF, a member of the AP1 family of transcription factors. Biochem. Biophys. Res. Commun. 229:316-322[CrossRef][Medline].
61. Wathelet, M. G., I. M. Clauss, J. Content, and G. A. Hues. 1988. The IFI-56K and IFI-54K interferon-inducible human genes belong to the same gene family. FEBS Lett. 231:164-171[CrossRef][Medline].
62. Werness, B. A., A. J. Levine, and P. M. Howley. 1990. Association of human papillomavirus types 16 and 18 E6 proteins with p53. Science 248:76-79[Abstract/Free Full Text].
63. Wu, K.-J., C. Grandori, M. Amacker, N. Simon-Vermont, A. Polack, J. Lingner, and R. Dalla-Favera. 1999. Direct activation of TERT transcription by c-MYC. Nat. Genet. 21:220-224[CrossRef][Medline].
64. Yang, D., O. Chertov, S. N. Bykovskaia, Q. Chen, M. J. Buffo, J. Shogan, M. Anderson, J. M. Schroder, J. M. Wang, O. M. Z. Howard, and J. J. Oppenheim. 1999. beta -Defensins: linking innate and adaptive immunity through dendritic and T cell CCR6. Science 286:525-528[Abstract/Free Full Text].
65. Zhu, H., J.-P. Cong, G. Mamtora, T. Gingeras, and T. Shenk. 1998. Cellular gene expression altered by human cytomegalovirus: global monitoring with oligonucleotide arrays. Proc. Natl. Acad. Sci. USA 95:14470-14475[Abstract/Free Full Text].
66. zur Hausen, H., and E. M. de Villiers. 1994. Human papillomaviruses. Annu. Rev. Microbiol. 48:427-447[Medline].

Journal of Virology, May 2000, p. 4174-4182, Vol. 74, No. 9
0022-538X/00/$04.00+0
Copyright © 2000, American Society for Microbiology. All rights reserved.


This article has been cited by other articles: