Cancer Research CR Helping Patients  Advances in Breast Cancer Research
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
Cancer Research Clinical Cancer Research
Cancer Epidemiology Biomarkers & Prevention Molecular Cancer Therapeutics
Molecular Cancer Research Cancer Prevention Research
Cancer Prevention Journals Portal Cancer Reviews Online
Annual Meeting Education Book Meeting Abstracts Online

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 arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Alazawi, W.
Right arrow Articles by Coleman, N.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Alazawi, W.
Right arrow Articles by Coleman, N.
[Cancer Research 62, 6959-6965, December 1, 2002]
© 2002 American Association for Cancer Research


Molecular Biology and Genetics

Changes in Cervical Keratinocyte Gene Expression Associated with Integration of Human Papillomavirus 161

William Alazawi2, Mark Pett2, Barbara Arch, Laurie Scott, Tom Freeman, Margaret A. Stanley and Nicholas Coleman3

Medical Research Council Cancer Cell Unit, MRC/Hutchison Research Centre, Cambridge CB2 2XZ [W. A., N. C.]; Department of Pathology, University of Cambridge, Cambridge CB2 1QP [M. P., M. A. S., N. C.]; Institute of Public Health, Cambridge CB2 25R [B. A.]; and Medical Research Council Human Genome Mapping Resource Centre, Cambridge CB10 1SB [L. S., T. F.], United Kingdom


    ABSTRACT
 Top
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 REFERENCES
 
Episomal integration is a critical event in human papillomavirus (HPV)-related oncogenesis, although little information is currently available concerning the effect of integration on the host transcriptome. We have used expression microarrays to investigate the effect of integration of HPV16 on gene expression in cervical keratinocytes, using the unique cell line model W12. W12 was generated from a cervical low-grade squamous intraepithelial lesion "naturally" infected with HPV16 and at low passage contains approximately 100 HPV16 episomes/cell. With passage in vitro, integration of viral episomes is associated with the development of phenotypic and genomic abnormalities resembling those seen in cervical neoplastic progression in vivo. We have used the Affymetrix U95A oligonucleotide array that contains probes for 12,600 human transcripts and have identified 85 genes from a range of host cell pathways that show changes in expression levels after integration of HPV16. Whereas some of the genes have previously been implicated in HPV-related oncogenesis in vivo, we have also identified a range of genes not previously described as being involved in cervical neoplastic progression. Interestingly, integration is associated with up-regulation of numerous IFN-responsive genes, in comparison with a baseline of episomally infected cells. These genes include p48, a component of the primary regulator of the IFN response pathway, IFN-stimulated gene factor 3. The physical state of high-risk HPV may substantially influence the response to IFN in infected keratinocytes.


    INTRODUCTION
 Top
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 REFERENCES
 
SCC4 of the cervix evolves through well-defined noninvasive stages, which may be classified using a two-tier system as LG-SIL and HG-SIL. Cervical squamous neoplasia is strongly associated with infection with HPV (1, 2, 3) , and integration of high-risk HPV—such as types 16, 18, 31, and 33—is associated with malignant progression. Whereas in premalignant cervical lesions the HPV genome is typically maintained in its episomal form (4) , integration has been shown in over 80% of HPV16-positive SCCs (4) and in 100% of HPV18-positive SCCs (5 , 6) .

HPV integration is usually accompanied by disruption of the viral E1 and E2 genes, with release of repression of the transforming genes E6 and E7 (7) . It is likely that HPV infection and integration cooperate with changes in the host genome and transcriptome in the multistep process of cervical carcinogenesis (1) . However, little information is currently available concerning changes in host gene expression that are associated with integration of HPV in cervical keratinocytes.

Two recent investigations have used cDNA microarrays to examine changes in gene expression in human foreskin keratinocytes after transfection with high-risk HPV31 (8) and low-risk HPV11 (9) . In both studies, the virus remained in the episomal state. Interestingly, 11 IFN-responsive genes were down-regulated after transfection with HPV31 episomes (8) . In contrast, 2 of these 11 genes [IFN-induced p27 (11.5 kDa) and IFN-induced 56 kDa] were reported to show increased expression levels (and the other 9 no change) after transfection with HPV11 episomes (9) . A third study used cDNA microarrays to demonstrate reduced expression of at least 18 IFN-responsive genes after infection of cervical keratinocytes with retroviruses encoding high-risk HPV16 E6 and E7 and in vitro induction of growth arrest and terminal differentiation (10) .

In the present study, we have investigated changes in host gene expression in cervical keratinocytes that are associated with integration of high-risk HPV16, the HPV type most commonly found in cervical carcinomas (1, 2, 3) . We have used as a model the W12 cell line, which was generated from a cervical LG-SIL and, thus, arose after "natural" cervical infection with HPV16 in vivo (11) . At early passages, W12 is diploid and contains approximately 100 HPV16 episomes/cell (11) . Our group has shown in a parallel study5 that episomal integration (which typically occurs between passages 10 and 22) is accompanied by an increase in the grade of lesion that W12 recapitulates in organotypic culture, enhanced colony-forming efficiency and growth rate, and a stepwise accumulation of host cytogenetic abnormalities that are similar to those seen in cervical neoplasia in vivo (12 , 13) . W12 is, therefore, a valuable and unique reagent, being a naturally infected cervical keratinocyte cell line that represents an accurate model of progression of HPV16-associated cervical neoplasia.

We have used oligonucleotide microarrays to investigate changes in expression of 12,600 human transcripts and expressed sequence tags associated with integration of HPV16 in W12. We have examined three passages of W12 in which we have characterized the physical state of the HPV16 genome by Southern hybridization. In a separate study, all passes of W12 tested were shown by chromosome painting to be polyclonal.5 The passages analyzed were:

(a) W12p10, which contains episomal HPV16 only and is near diploid. We used W12p10 as a baseline for analysis of gene expression changes in higher passages.

(b) W12p22, which contains integrated HPV16 only and is near tetraploid. W12p22 shows cytogenetic abnormalities that are not seen in W12p10 or W12p14, most notably a high-level gain of chromosome 5p, which is present in approximately 50% of metaphases of W12p22 and becomes more frequent in later passages.5 Comparison of the transcriptome between W12p22 and W12p10 would enable us to identify changes in cervical keratinocyte gene expression resulting from the complete process of integration.

(c) W12p14, which contains mixed episomal and integrated HPV16 and is near tetraploid. In the present study, we show that there is an approximately 70% reduction in HPV16 episomal signal by Southern blot in W12p14 compared with W12p10. Comparison of gene expression between W12p14 and W12p10 would allow us to examine effects of early disruption of the HPV16 episome (which we also show in the present study to result in increased levels of HPV16 E7 protein and loss of HPV16 E2 mRNA) and the critical switch in the host genome from diploidy to tetraploidy. Comparison of gene expression between W12p22 and W12p14 would enable us to examine effects of complete loss of HPV16 episomes and the associated cytogenetic abnormalities (e.g., high-level gain of chromosome 5p).

All passages of W12 analyzed were grown in monolayer tissue culture, which selects for cells with a proliferative advantage and limits differentiation (14) . These conditions enable W12 to serve as a model of HPV16 infection of basal cervical squamous cells. This is of particular importance, because derepression of viral E7 and E6 genes in the basal layers of the cervix, as occurs after episomal integration, is an important early event in cervical carcinogenesis (15) (16) .

We have identified genes from several host cellular pathways that show changes in expression levels during the process of integration of HPV16 episomes in cervical keratinocytes. Several of the changes have been demonstrated previously to occur in cervical neoplasia in vivo, further confirming the validity of W12 as a model system for cervical oncogenesis. We have also identified numerous other, previously unreported changes in gene expression that may also contribute to cervical oncogenesis in vivo.


    MATERIALS AND METHODS
 Top
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 REFERENCES
 
Monolayer Cell Culture and Extraction of Nucleic Acids.
W12 keratinocytes were grown in Glasgow’s modified Eagle’s medium supplemented with FCS and EGF, in the presence of irradiated murine 3T3 J2 fibroblast feeder cells, as described previously (11 , 17) . Culture conditions were the same for all passages of W12 examined. Normal ectocervical cells were obtained after hysterectomy for disease unrelated to the cervix and were cultured as described previously (17) . At approximately 80% confluence, total RNA and DNA were extracted from the keratinocytes using the TRIzol extraction kit (Life Technologies, Inc., Grand Island, NY). All samples of W12 investigated were polyclonal.5 In all cases, RNA and DNA were obtained from at least three independently grown tissue culture flasks.

Southern Analysis of the Physical State of HPV16.
DNA (10 µg) from W12 (at passes 10, 14, 15, 22, 34, 44, and 55) and 1 µg of DNA from the cervical SCC cell line CaSki (18) were digested with BamHI or HindIII and electrophoresed through 0.8% agarose. For lanes loaded with 10 µg of DNA, we used copy number controls representing 1, 10, and 100 copies of linearized HPV16/diploid cell. After partial depurination, denaturation, and neutralization, DNA was transferred to Hybond-N+ nylon membrane (Amersham, Little Chalfont, United Kingdom). Prehybridization and hybridization of 32P-labeled HPV16 probe were performed at 65°C using Rapid-Hyb buffer (Amersham) for 1 h and 2.5 h, respectively. The probe for hybridization was prepared by excision and purification of full-length HPV16 from the pspHPV16 plasmid, followed by labeling with {alpha}32P-dCTP using random priming. After hybridization membranes were washed and exposed to autoradiographic film (Amersham) at -70°C using intensifying screens.

Western Blotting for HPV16 E7 Protein.
Total protein was extracted from at least three separate flasks of W12 (at passes 10, 14, 23, 35, 43, and 56), the cervical SCC cell line CaSki (18) , and normal ectocervical keratinocytes at approximately 80% confluence. The protein extraction buffer contained 250 mM NaCl, 0.1% NP40, 50 mM HEPES, and one protease inhibitor tablet (Roche, Basel, Switzerland) per 25 ml. Extracts were stored in 10–20% glycerol at -70°C. Twenty micrograms of total protein was run on a 12.5% polyacrylamide gel, blotted, and analyzed by hybridization of a sheep polyclonal antibody against HPV16 E7 (Ref. 19 ; a kind gift from Cantab Pharmaceuticals, Cambridge, United Kingdom).

RT-PCR for HPV16 E2 mRNA.
Total RNA was extracted from at least three separate flasks of W12 (at passes 10, 15, 24, and 55), the cervical SCC cell lines CaSki (18) and SiHa (20) , and normal ectocervical keratinocytes. As a positive control, we used the plasmid pspHPV16 (11) , which contains the full-length HPV16 genome. The RNA was reverse transcribed using a Superscript reverse transcription kit (Life Technologies, Inc.), including a reverse transcriptase-negative control for each sample. The resultant samples were amplified using primers for E2 (forward, 5'-TGCGATGGATCCATGGAGACTCTTTGCCAACG-3'; reverse, 5'-TGCGATGGATTCTCATATAGACATAAATCCAG-3'), and the products were run on a 1% agarose gel. The expected product size was 1139 bp.

Hybridization to Affymetrix Arrays.
Biotin-labeled target cRNA was generated from total RNA from W12p10, W12p14, and W12p22. We generated two assay replicates from total RNA from W12p10 and W12p14 and three assay replicates from total RNA from W12p22. Generation of cDNA, biotin labeling, and fragmentation of cRNA were performed according to Affymetrix (Santa Clara, CA) protocols. In brief, double-stranded cDNA was generated using the T7 (-dT)24 primer 5'-GGCCAGTGAATTGTAATACGACTCACTATAGGGAGGCGG-3' and Superscript II reverse transcriptase (Life Technologies, Inc., Gaithersburg, MD). Biotin-labeled cRNA was generated from this using the Enzo BioArray HighYield RNA Transcript labeling kit (Enzo, New York, NY).

cRNA from each assay replicate of each passage of W12 was hybridized to an Affymetrix human U95A oligonucleotide array slide for 16 h at 45°C in a hybridization mixture (Affymetrix). Each array contained probes for 12,600 transcripts. After hybridization the slides were washed and stained using fluidics protocols described in the Affymetrix manual.

Analysis of Microarray Data.
Gene expression was assessed using algorithms designed by Affymetrix in the Microarray Suite software package (21) . For each sample, genes were classified as either "present," "marginal," or "absent," with a numerical indication of average intensity difference between the perfect match and central-base mismatch features. A global normalization algorithm and linear regression scaling factor were applied to normalize data from different hybridizations. Multiple parameters (21) were used to generate numerical fold changes reflecting the relative abundance of each transcript in each pair of samples. Changes in expression were also assigned into one of five discrete categories: "decrease," "marginal decrease," "no change," "marginal increase," and "increase."

Affymetrix data were analyzed further using Microarray Suite Empirical and Statistical expression algorithms (Affymetrix) and ChipMining Tool software (Gareth Denyer, Sydney, Australia). When comparing gene expression between passages of W12, we excluded genes from further analysis if they were in the "no change," "marginal decrease," or "marginal increase" categories and if their numerical fold change was between 2 and -2. Each replicate sample was compared with each sample from the other passage, and only gene expression changes that were identified consistently were analyzed further. These exclusion criteria were similar to those reported by others (22) .

Analysis of Gene Distribution.
We investigated whether there was overall clustering of the chromosomal locations of genes showing differences in expression levels between passages of W12. We compared the observed distribution of these genes with that expected under the null hypothesis that the genes were evenly distributed throughout the genome. Under this hypothesis, the number of genes located on each chromosome arm was determined from mapping information provided by Affymetrix. A multinomial distribution was used and significance was assessed using the {chi}2 goodness-of-fit test with a 5% significance level.

When evidence of clustering was detected, a post hoc analysis was carried out to identify whether particular arms were observed as the location of a significant number of genes. In this analysis, only arms representing the location of three or more genes were investigated. Significance was assessed by comparing observed numbers with those expected under the null hypothesis that the number of genes on an arm were binomially distributed, and that genes were evenly distributed throughout the genome. A 2.5% significance level was used to help correct for problems of multiple testing that may have been incurred.


    RESULTS
 Top
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 REFERENCES
 
Southern Analysis of the Physical State of HPV16.
For W12p10 BamHI (a HPV16 single cutter) produced a band of 7.9 kb (Fig. 1A)Citation , representing full-length HPV16 episomes at approximately 100 copies/cell, consistent with previous observations (11) . Digestion with HindIII (a HPV16 noncutter) revealed bands consistent with monomeric episomal HPV16, including low molecular weight supercoiled monomers (Fig. 1B)Citation . Also present were high molecular weight (>30 kb) bands suggesting the presence of head-to-tail multimeric episomal DNA, which was subsequently confirmed by Southern analysis after two dimensional gel electrophoresis (data not shown). Thus, W12p10 harbors full-length HPV16 episomes in monomeric and head-to-tail multimeric forms.



View larger version (82K):
[in this window]
[in a new window]
[Download PPT slide]
 
Fig. 1. Southern analysis of the physical state of HPV16 in W12 (10 µg of DNA from each passage) and CaSki (1 µg) restricted with either BamHI (A) or HindIII (B). For lanes loaded with 10 µg of DNA, we used copy number controls representing 1, 10, and 100 copies of linearized HPV16/diploid cell (right lanes). Also shown are size markers in kilobase. W12p10 contains approximately 100 copies/cell of full-length HPV16 episomes in monomeric and head-to-tail multimeric forms, the latter being represented by the high molecular weight (>30 kb) bands in the HindIII digest. The episomal bands show a rapid reduction in intensity in the higher passages studied, with concomitant emergence of bands representing integrated HPV16. There is a single host integration site containing approximately 10 copies of a single type of viral integrant, which is maintained at the same copy number to pass 55. The minor bands in the BamHI digest of W12p10 were due to incomplete digestion of episomes, rather than representing virus-host junctional fragments. The minor bands were of the same size as those seen in the HindIII digest of W12p10 and different to the two BamHI bands produced by integration (passes 14–55). They were abrogated following longer digestion with BamHI (data not shown).

 
Analysis of W12 at the higher passages used in this study revealed a rapid reduction in the intensity of bands representing episomal DNA, with approximately 70% loss in W12p14, 90% loss in W12p15, and absence in W12p22 (Fig. 1)Citation . There was concurrent emergence of bands representing integrated HPV16 DNA. We observed one band after HindIII digestion and two bands after BamHI digestion, indicating a single host integration site containing a single type of viral integrant with one BamHI site. The integrated genomes were maintained at the same copy number in subsequent passages. The cells containing integrated HPV16 were tetraploid,5 and copy number controls suggested approximately 10 copies of integrated HPV16/cell, consistent with amplification of viral DNA after integration of the HPV16 genome (23) .

Expression of HPV16 E7 Protein.
HPV16 E7 was detectable by Western blotting in W12 at all passages between W12p10 and W12p56 (Fig. 2)Citation . Episomal integration was associated with an increase in levels of detectable HPV16 E7 protein, with greater levels of expression in W12p14 compared with W12p10 and in W12p23 compared with W12p14 (Fig. 2)Citation . HPV16 E7 protein was also detectable in the cervical SCC cell line CaSki but was absent from normal ectocervical keratinocytes.



View larger version (31K):
[in this window]
[in a new window]
[Download PPT slide]
 
Fig. 2. Western blot showing that HPV16 E7 is detectable in W12 at all passages from W12p10 to W12p56. Episomal integration (which occurs between W12p10 and W12p23) is associated with an increase in levels of detectable HPV16 E7 protein. HPV16 E7 is also detectable in the cervical SCC cell line CaSki, but not in normal ectocervical keratinocytes (NCx).

 
Expression of HPV16 E2 mRNA.
HPV16 E2 mRNA was detectable by RT-PCR in W12p10 (episomes only), but was lost during the process of integration in W12p15 and W12p24 (Fig. 3)Citation . HPV16 E2 was detectable in the cervical SCC cell line CaSki, in which the full-length early region is retained in the integrated HPV16 DNA (24) , but was absent from the cervical SCC cell line SiHa, in which the E2 gene is deleted, and from normal ectocervical keratinocytes.



View larger version (35K):
[in this window]
[in a new window]
[Download PPT slide]
 
Fig. 3. RT-PCR showing that HPV16 E2 is present in W12p10 (episomes only) but is lost during the process of integration in W12p15 and W12p24. HPV16 E2 is detectable in the cervical SCC cell line CaSki, in which the full-length early region is retained in the integrated HPV16 but is absent from the cervical SCC cell line SiHa, in which the E2 gene is deleted, and from normal ectocervical keratinocytes (NCx). Each reaction was performed in duplicate, with and without reverse transcriptase (left and right lanes of each pair, respectively). The plasmid pspHPV16 (11) , which contains the full-length HPV16 genome, was used as a positive control and we included a no-template negative control.

 
Microarray Analysis of Host Gene Expression.
We analyzed two replicates of labeled RNA from W12p10, two replicates of labeled RNA from W12p14, and three replicates of labeled RNA from W12p22. The mean number of genes called as "marginal" or "present" was 3573 (range, 3480–3667) in W12p10, 4166 (range, 4026–4306) in W12p14, and 3981 (range, 3727–4454) in W12p22. The reproducibility of the Affymetrix arrays was confirmed by comparing expression levels of all genes in all available assay replicates from each passage of W12. In this analysis, fewer than 0.04% of genes were described as having changed expression levels by more than 2-fold.

Each replicate sample was compared with each sample from the other two passages. The fold changes for expression levels in W12p14 versus W12p10 represent the mean of four comparisons, and the fold changes for expression levels in W12p22 versus W12p10 and W12p22 versus W12p14 represent the mean of six comparisons.

The stylized Venn diagram in Fig. 4Citation depicts the patterns of changes in gene expression levels (identified using the exclusion criteria set by us) during the process of HPV16 integration in W12. Sets A-D are those of primary interest in analyzing these data:



View larger version (21K):
[in this window]
[in a new window]
[Download PPT slide]
 
Fig. 4. Stylized Venn diagram depicting the patterns of changes in gene expression levels during the process of HPV16 integration in W12. The principal sets contain the genes showing differences between different passages of W12. Regions of overlap contain genes common to two or more comparisons. The Table gives the numbers of genes that show increased or decreased expression levels at higher passage, subdivided by set.

 
(a) Set A contains those genes that show a consistent and continued change throughout both steps of the integration event examined (i.e., W12p14 versus W12p10 and W12p22 versus W12p14, as well as W12p22 versus W12p10).

(b) Set B contains genes for which expression changes in the early stages of the integration process (W12p14 versus W12p10), with persistence of the change at completion of integration (W12p22 versus W12p10) but without additional change between W12p14 and W12p22.

(c) Set C contains genes showing changes in expression during the later stages of the integration process (W12p22 versus W12p14), without change between W12p10 and W12p14, producing an overall change in expression levels between W12p22 and W12p10.

(d) Set D contains genes showing changes in expression during integration overall (W12p22 versus W12p10), but without demonstrable changes in the early or late stages of the integration process.

(e) Sets E and G, respectively, represent genes showing changes in expression in W12p14 versus W12p10 and W12p22 versus W12p14. However, these genes were not detected as showing significant changes overall during the process of integration (W12p22 versus W12p10) and were, therefore, not considered further by us. We predicted that set F would be empty (because any gene falling into set E and set G should fall into set A), and this was indeed the case.

We observed more changes in gene expression in the early stages of the integration process (between W12p14 and W12p10, when 70% of HPV16 episomes are lost) than in the later stages of the integration process (between W12p22 and W12p14, when the remaining 30% of episomes are lost; Fig. 4Citation ). Sets A-D contained 85 genes, of which 53 showed an increase in expression levels and 32 showed a decrease in expression levels. These 85 genes are listed in Table 1Citation , where they are divided according to the sets of the Venn diagram in Fig. 4Citation and then ranked by mean fold change between W12p22 and W12p10. We observed increased levels of gene expression of up to 12.3-fold and reduced levels of gene expression of up to 11.6-fold. Expression changes for 5 of the 85 genes were confirmed by RT-PCR and densitometry, using RNA from cells grown separately from those used in the microarray experiments (data not shown). These genes were CIG49 (increased, set C, GenBank accession no. AF026939); BST-2 (increased, set C, GenBank accession no. D28137); IFN-induced 56 kDa (increased, set C, GenBank accession no. M24594); expressed in osteoblasts (increased, set C, GenBank accession no. AB000115) and elafin (decreased, set A, GenBank accession no. L10343).


View this table:
[in this window]
[in a new window]

 
Table 1 Summary of gene expression changes

The genes are from sets A-D in Fig. 4Citation , sorted by set, and ranked by mean fold change between W12p22 and W12p10. The GenBank accession numbers are those supplied by Affymetrix.

 
Many IFN-responsive genes showed increased levels of expression in W12p22 compared with W12p10 (Table 1)Citation . Eight of these, showing fold changes of up to 10.1, were among 11 IFN-responsive genes reported to show reduced levels of expression in foreskin keratinocytes transfected with HPV31 episomes, compared with a normal keratinocyte baseline (8) . These genes included p48, a component of the primary regulator of the IFN response pathway, IFN-stimulated gene factor 3, which showed a 3.9-fold increase in expression levels in W12p22 compared with W12p10.

Fig. 5Citation shows the known chromosomal locations of the genes contained in sets A-D (Fig. 4Citation and Table 1Citation ). The locations of the genes were not consistent with a multinomial distribution, under the null hypothesis that the genes would be distributed evenly throughout the genome (P = 0.024 for genes increasing in expression, P = 0.03 for genes decreasing in expression, {chi}2 test on 40 degrees of freedom). Hence, it is unlikely that the distribution of the genes was determined by chance alone. We observed clustering of genes showing increased expression after integration, with significant clustering on 8p (four genes, P = 0.006) and 10q (five genes at 10q23–24, P = 0.01). We also observed clustering of genes showing reduced expression after integration, with significant clustering on 10p (three genes at 10p15-p14, P = 0.006) and a trend toward clustering on 20q (three genes, P = 0.045). There was no change in the levels of expression of any gene on chromosome 5p between W12p10 and later passages.



View larger version (13K):
[in this window]
[in a new window]
[Download PPT slide]
 
Fig. 5. Expression karyogram showing the chromosomal locations (where known) of the genes in sets A-D in Fig. 4Citation (80 locations in total). Green dots to the right of each chromosome represent genes showing increased expression at higher passage, and red dots to the left of each chromosome represent genes showing reduced expression at higher passage.

 

    DISCUSSION
 Top
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 REFERENCES
 
We have used W12 to study the changes in expression levels of 12,600 transcripts that occur during integration of HPV16 episomes in cervical keratinocytes. W12 was derived from a naturally infected cervical LG-SIL (11) and represents a valuable and unique model for investigations of the natural history of cervical squamous cell neoplasia. During passage of W12, the viral integration that we have demonstrated is associated with morphological and cytogenetic changes similar to those seen in cervical neoplastic progression in vivo.5 In the present study, we have shown that integration of HPV16 in W12 is associated with loss of expression of E2 mRNA and with increased levels of E7 protein (consistent with loss of E2-mediated repression (25, 26, 27) . It is likely that these events exert direct or indirect effects on the expression of host genes in W12 cells. W12p10 was, therefore, used as a baseline for comparative analysis of changes in gene expression seen at higher passages.

Analysis of W12 by comparative genomic hybridization has shown that progression during passage is associated with the acquisition of nonrandom copy number imbalances.5 In particular, between W12p10 and W12p22 we observed evidence of loss of 10p and gain of 10q as a result of the formation of an isochromosome (10q), which we have confirmed by spectral karyotyping (19) . This isochromosome is then maintained through to high passage. Our present data showing significant clustering of genes up-regulated after HPV integration on 10q and of genes down-regulated after HPV integration on 10p is consistent with our comparative genomic hybridization and spectral karyotyping data.

On the other hand, not all genes showing expression differences in the present study correspond to regions of copy number imbalance detectable in W12 by comparative genomic hybridization. Although such disparities may, in some cases, reflect the limited resolution of comparative genomic hybridization using chromosome targets, epigenetic factors are also likely to underlie the altered expression of some of the genes that we have identified using microarrays. Despite the high level gain of chromosome 5p seen by us in W12p14 and W12p22 but not W12p10,5 there was no change in the levels of expression of genes on chromosome 5p between W12p10 and later passages. This is consistent with recent data showing that chromosomal regions amplified in neoplasms may be transcriptionally silent (28) .

The 85 genes in Table 1Citation are of potential significance in the biology of cervical squamous cell neoplasia in vivo. Several of the genes have been demonstrated in previous studies to change their expression levels during the development of SCC in the cervix. These changes include increased expression at the protein level of thrombomodulin (29) and the antiapoptotic Bcl-2 family member Mcl-1 (Ref. 30 ; both in set B). There are also genes for which changes in expression levels have been shown in malignancies at other sites, but not in the cervix, to date, such as acid ceramidase (set D), which is up-regulated at the mRNA level in prostatic adenocarcinoma (31) ; and Mac2-binding protein (set B), which is induced at the protein level in breast and other adenocarcinomas (32) .

Several genes with immunomodulatory function showed altered expression in W12 after integration. This observation is of interest, in view of the apparent switch from local immune quiescence in cervical LG-SIL to a local cell-mediated immune response in HG-SIL and SCC (17) . We observed up-regulation of HLA class IC (set B), HLA class IF (set D), and the immunostimulatory molecule Mac 2-binding protein (set B; Ref. 33 ). We also observed that numerous IFN-responsive genes were up-regulated following integration. These genes include p48 (set C), a component of the primary regulator of the IFN response pathway, IFN-stimulated gene factor 3. These data suggest that the physical state of high-risk HPV may substantially influence the response to IFN in infected keratinocytes.

Several of the changes in gene expression that we have observed in W12 would contribute to remodeling of the ECM in vivo. Whereas some changes, such as up-regulation of laminin ß1 (set D) and collagen XVII{alpha}1 (set D) and down-regulation of type IV collagenase (set B), would promote deposition of ECM, other changes, such as down-regulation of TIMP1 (set B) and elastase inhibitor (set D), would facilitate ECM degradation and might contribute to stromal invasion. Integration is also associated with reduced expression of elafin (set A), an anti-leukoproteinase with specific anti-elastase activity (34) that is detectable in many types of epithelium. Expression of elafin mRNA and protein in SCC has been shown to correlate inversely with grade, with strong expression in the keratinocytes of well-differentiated SCC but absence from poorly differentiated SCC in the skin (35) and esophagus (36) .

Both IGF-II and EGFR showed high levels of up-regulation after HPV16 integration in W12 (both set B), with fold changes of 12.3 and 4.7, respectively. This is consistent with evidence that autocrine secretion of IGF-II, acting through the IGF-I receptor, participates in EGF-induced mitogenesis in cervical keratinocytes (37) . We also demonstrated up-regulation of an EGF response factor, ERF 2, a zinc finger protein with unknown function (set B), which would be consistent with increased activity in the EGFR signaling pathway after episomal integration. IGF-II and EGFR may, therefore, cooperate in stimulating growth of W12 after integration and may also play a role in cervical neoplastic progression in vivo. Indeed, in one study, cervical HG-SIL and SCC were shown to express higher levels of EGFR and to be associated with elevated serum levels of IGF-II (38) .

The experiments reported here were performed on W12 cells in monolayer culture, to provide a model of HPV16 infection of basal cervical squamous cells. It is likely that additional changes in gene expression would become apparent with greater epithelial differentiation. As much of the HPV life cycle is dependent on the maturation state of host keratinocytes, changes in gene expression would be likely to depend not only on the physical state of HPV but also on the differentiation state of the infected keratinocytes. Indeed, a recent study of the effects of HPV16 E6 and E7 on cervical keratinocytes (10) detected more gene expression changes in differentiating than in proliferating cells.

In conclusion, the W12 cell line has enabled us to study the effect of episomal integration on host gene expression in cervical keratinocytes naturally infected with HPV16, as a component of progression toward malignancy. We have identified functional groups of genes showing changes in expression as a result of integration and have further identified the stages in the integration process at which expression levels change. This data set can be used to direct further studies into mechanisms of progression of premalignant lesions in the cervix, and potentially at other sites.


    ACKNOWLEDGMENTS
 
We are grateful to Richard Stephens (MRC United Kingdom HGMP Resource Centre, Hinxton, Cambridge, United Kingdom) for assistance with data analysis. We also thank Gareth Denyer (University of Sydney, New South Wales, Australia) for the use of his ChipMining Tool software and Cantab Pharmaceuticals (Cambridge, United Kingdom) for providing the anti-HPV16E7 antibody.


    FOOTNOTES
 
The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.

1 Funded by the Medical Research Council. Back

2 Joint first authors. Back

3 To whom requests for reprints should be addressed, at Medical Research Council Cancer Cell Unit, MRC/Hutchison Research Centre, Hills Road, Cambridge CB2 2XZ, United Kingdom. Phone: 44-1223-763285; Fax: 44-1223-763284; E-mail: nc109{at}cam.ac.uk Back

4 The abbreviations used are: SCC, squamous cell carcinoma; ECM, extracellular matrix; EGF, epidermal growth factor; EGFR; EGF receptor, HG-SIL, high-grade squamous intraepithelial lesion; HPV, human papillomavirus; IGF, insulin-like growth factor; LG-SIL, low-grade SIL. Back

5 Personal communications; 19th International Papillomavirus Conference, 2001, Florianopolis, Brazil: Alazawi et al., Changes in the host cell genome and transcriptome associated with integration of HPV16 in cervical keratinocytes; and Dowen et al., Gain of chromosome arm 5p: a recurrent late cytogenetic event in HPV16 and HPV18 expressing cell lines. Back

Received 4/ 4/02. Accepted 10/ 4/02.


    REFERENCES
 Top
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 REFERENCES
 

  1. Schiffman M. H., Bauer H. M., Hoover R. N., Glass A. G., Cadell D. M., Rush B. B., Scott D. R., Sherman M. E., Kurman R. J., Wacholder S., et al Epidemiologic evidence showing that human papillomavirus infection causes most cervical intraepithelial neoplasia [see comments]. J. Natl. Cancer Inst., 85: 958-964, 1993.[Abstract/Free Full Text]
  2. Schiffman M. H., Brinton L. A. The epidemiology of cervical carcinogenesis. Cancer (Phila.), 76: 1888-1901, 1995.[Medline]
  3. Arends M. J., Buckley C. H., Wells M. Aetiology, pathogenesis, and pathology of cervical neoplasia. J. Clin. Pathol., 51: 96-103, 1998.[Abstract]
  4. Kalantari M., Blennow E., Hagmar B., Johansson B. Physical state of HPV16 and chromosomal mapping of the integrated form in cervical carcinomas. Diagn. Mol. Pathol., 10: 46-54, 2001.[Medline]
  5. Berumen J., Casas L., Segura E., Amezcua J. L., Garcia-Carranca A. Genome amplification of human papillomavirus types 16 and 18 in cervical carcinomas is related to the retention of E1/E2 genes. Int. J. Cancer, 56: 640-645, 1994.[Medline]
  6. Park J. S., Hwang E. S., Park S. N., Ahn H. K., Um S. J., Kim C. J., Kim S. J., Namkoong S. E. Physical status and expression of HPV genes in cervical cancers. Gynecol. Oncol., 65: 121-129, 1997.[Medline]
  7. Lazo P. A. The molecular genetics of cervical carcinoma. Br. J. Cancer, 80: 2008-2018, 1999.[Medline]
  8. Chang Y. E., Laimins L. A. Microarray analysis identifies interferon-inducible genes and Stat-1 as major transcriptional targets of human papillomavirus type 31. J. Virol., 74: 4174-4182, 2000.[Abstract/Free Full Text]
  9. Thomas J. T., Oh S. T., Terhune S. S., Laimins L. A. Cellular changes induced by low-risk human papillomavirus type 11 in keratinocytes that stably maintain viral episomes. J. Virol., 75: 7564-7571, 2001.[Abstract/Free Full Text]
  10. Nees M., Geoghegan J. M., Hyman T., Frank S., Miller L., Woodworth C. D. Papillomavirus type 16 oncogenes down-regulate expression of interferon-responsive genes and up-regulate proliferation-associated and NF{kappa}B-responsive genes in cervical keratinocytes. J. Virol., 75: 4283-4296, 2001.[Abstract/Free Full Text]
  11. Stanley M. A., Browne H. M., Appleby M., Minson A. C. Properties of a non-tumorigenic human cervical keratinocyte cell line. Int. J. Cancer, 43: 672-676, 1989.[Medline]
  12. Heselmeyer K., Macville M., Schrock E., Blegen H., Hellstrom A. C., Shah K., Auer G., Ried T. Advanced-stage cervical carcinomas are defined by a recurrent pattern of chromosomal aberrations revealing high genetic instability and a consistent gain of chromosome arm 3q. Genes Chromosomes Cancer, 19: 233-240, 1997.[Medline]
  13. Kirchhoff M., Rose H., Petersen B. L., Maahr J., Gerdes T., Philip J., Lundsteen C. Comparative genomic hybridization reveals non-random chromosomal aberrations in early preinvasive cervical lesions. Cancer Genet. Cytogenet., 129: 47-51, 2001.[Medline]
  14. Stanley M. A., Greenfield I. M. Culture of human cervical epithelial cells Freshney R. I. eds. . Culture of Epithelial Cells, 135-158, Wiley-Liss, New York: New York, NY 1992.
  15. Durst M., Glitz D., Schneider A., zur Hausen H. Human papillomavirus type 16 (HPV 16) gene expression and DNA replication in cervical neoplasia: analysis by in situ hybridization. Virology, 189: 132-140, 1992.[Medline]
  16. Higgins G. D., Uzelin D. M., Phillips G. E., McEvoy P., Marin R., Burrell C. J. Transcription patterns of human papillomavirus type 16 in genital intraepithelial neoplasia: evidence for promoter usage within the E7 open reading frame during epithelial differentiation. J. Gen. Virol., 73: 2047-2057, 1992.[Abstract/Free Full Text]
  17. Coleman N., Greenfield I. M., Hare J., Kruger-Gray H., Chain B. M., Stanley M. A. Characterization and functional analysis of the expression of intercellular adhesion molecule-1 in human papillomavirus-related disease of cervical keratinocytes. Am. J. Pathol., 143: 355-367, 1993.[Abstract]
  18. Pattillo R. A., Hussa R. O., Story M. T., Ruckert A. C., Shalaby M. R., Mattingly R. F. Tumor antigen and human chorionic gonadotropin in CaSki cells: a new epidermoid cervical cancer cell line. Science (Wash DC), 196: 1456-1458, 1977.[Abstract/Free Full Text]
  19. Cottage A., Dowen S., Roberts I., Pett M., Coleman N., Stanley M. Early genetic events in HPV immortalized keratinocytes. Genes Chromosomes Cancer, 30: 72-79, 2001.[Medline]
  20. Friedl F., Kimura I., Osato T., Ito Y. Studies on a new human cell line (SiHa) derived from carcinoma of uterus. I. Its establishment and morphology. Proc. Soc. Exp. Biol. Med., 135: 543-545, 1970.[Medline]
  21. Schadt E. E., Li C., Su C., Wong W. H. Analyzing high-density oligonucleotide gene expression array data. J. Cell. Biochem., 80: 192-202, 2000.[Medline]
  22. Ly D. H., Lockhart D. J., Lerner R. A., Schultz P. G. Mitotic misregulation and human aging. Science (Wash DC), 287: 2486-2492, 2000.[Abstract/Free Full Text]
  23. Wagatsuma M., Hashimoto K., Matsukura T. Analysis of integrated human papillomavirus type 16 DNA in cervical cancers: amplification of viral sequences together with cellular flanking sequences. J. Virol., 64: 813-821, 1990.[Abstract/Free Full Text]
  24. Smits H. L., Cornelissen M. T., Jebbink M. F., van den Tweel J. G., Struyk A. P., Briet M., ter Schegget J. Human papillomavirus type 16 transcripts expressed from viral-cellular junctions and full-length viral copies in CaSki cells and in a cervical carcinoma. Virology, 182: 870-873, 1991.[Medline]
  25. Haugen T. H., Cripe T. P., Ginder G. D., Karin M., Turek L. P. Trans-activation of an upstream early gene promoter of bovine papilloma virus-1 by a product of the viral E2 gene. EMBO J., 6: 145-152, 1987.[Medline]
  26. Dong G., Broker T. R., Chow L. T. Human papillomavirus type 11 E2 proteins repress the homologous E6 promoter by interfering with the binding of host transcription factors to adjacent elements. J. Virol., 68: 1115-1127, 1994.[Abstract/Free Full Text]
  27. Francis D. A., Schmid S. I., Howley P. M. Repression of the integrated papillomavirus E6/E7 promoter is required for growth suppression of cervical cancer cells. J. Virol., 74: 2679-2686, 2000.[Abstract/Free Full Text]
  28. Platzer P., Upender M. B., Wilson K., Willis J., Lutterbaugh J., Nosrati A., Willson J. K., Mack D., Ried T., Markowitz S. Silence of chromosomal amplifications in colon cancer. Cancer Res., 62: 1134-1138, 2002.[Abstract/Free Full Text]
  29. Ordonez N. G. Thrombomodulin expression in transitional cell carcinoma. Am. J. Clin. Pathol., 110: 385-390, 1998.[Medline]
  30. Wei L. H., Kuo M. L., Chen C. A., Chou C. H., Cheng W. F., Chang M. C., Su J. L., Hsieh C. Y. The anti-apoptotic role of interleukin-6 in human cervical cancer is mediated by up-regulation of Mcl-1 through a PI 3-K/Akt pathway. Oncogene, 20: 5799-5809, 2001.[Medline]
  31. Seelan R. S., Qian C., Yokomizo A., Bostwick D. G., Smith D. I., Liu W. Human acid ceramidase is overexpressed but not mutated in prostate cancer. Genes Chromosomes Cancer, 29: 137-146, 2000.[Medline]
  32. Iacobelli S., Arno E., Sismondi P., Natoli C., Gentiloni N., Scambia G., Giai M., Cortese P., Panici P. B., Mancuso S. Measurement of a breast cancer associated antigen detected by monoclonal antibody SP-2 in sera of cancer patients. Breast Cancer Res. Treat., 11: 19-30, 1988.[Medline]
  33. Ullrich A., Sures I., D’Egidio M., Jallal B., Powell T. J., Herbst R., Dreps A., Azam M., Rubinstein M., Natoli C., et al The secreted tumor-associated antigen 90K is a potent immune stimulator. J. Biol. Chem., 269: 18401-18407, 1994.[Abstract/Free Full Text]
  34. Ameshima S., Ishizaki T., Demura Y., Imamura Y., Miyamori I., Mitsuhashi H. Increased secretory leukoprotease inhibitor in patients with nonsmall cell lung carcinoma. Cancer (Phila.), 89: 1448-1456, 2000.[Medline]
  35. Alkemade H. A., van Vlijmen-Willems I. M., van Haelst U. J., van de Kerkhof P. C., Schalkwijk J. Demonstration of skin-derived antileukoproteinase (SKALP) and its target enzyme human leukocyte elastase in squamous cell carcinoma. J. Pathol., 174: 121-129, 1994.[Medline]
  36. Yamamoto S., Egami H., Kurizaki T., Ohmachi H., Hayashi N., Okino T., Shibata Y., Schalkwijk J., Ogawa M. Immunohistochemical expression of SKALP/elafin in squamous cell carcinoma of the oesophagus. Br. J. Cancer, 76: 1081-1086, 1997.[Medline]
  37. Steller M. A., Delgado C. H., Zou Z. Insulin-like growth factor II mediates epidermal growth factor-induced mitogenesis in cervical cancer cells. Proc. Natl. Acad. Sci. USA, 92: 11970-11974, 1995.[Abstract/Free Full Text]
  38. Mathur S. P., Mathur R. S., Young R. C. Cervical epidermal growth factor-receptor (EGF-R) and serum insulin- like growth factor II (IGF-II) levels are potential markers for cervical cancer. Am. J. Reprod. Immunol., 44: 222-230, 2000.



This article has been cited by other articles:


Home page
J. Clin. Microbiol.Home page
M. Jiang, J. G. Baseman, L. A. Koutsky, Q. Feng, C. Mao, N. B. Kiviat, and L. F. Xi
Sequence Variation of Human Papillomavirus Type 16 and Measurement of Viral Integration by Quantitative PCR
J. Clin. Microbiol., March 1, 2009; 47(3): 521 - 526.
[Abstract] [Full Text] [PDF]


Home page
Cancer Res.Home page
K. L. Dall, C. G. Scarpini, I. Roberts, D. M. Winder, M. A. Stanley, B. Muralidhar, M. T. Herdman, M. R. Pett, and N. Coleman
Characterization of Naturally Occurring HPV16 Integration Sites Isolated from Cervical Keratinocytes under Noncompetitive Conditions
Cancer Res., October 15, 2008; 68(20): 8249 - 8259.
[Abstract] [Full Text] [PDF]


Home page
Cancer Res.Home page
S. Vinokurova, N. Wentzensen, I. Kraus, R. Klaes, C. Driesch, P. Melsheimer, F. Kisseljov, M. Durst, A. Schneider, and M. von Knebel Doeberitz
Type-Dependent Integration Frequency of Human Papillomavirus Genomes in Cervical Lesions
Cancer Res., January 1, 2008; 68(1): 307 - 313.
[Abstract] [Full Text] [PDF]


Home page
CarcinogenesisHome page
M.T. Herdman, M. R. Pett, I. Roberts, W. O.F. Alazawi, A. E. Teschendorff, X.-Y. Zhang, M. A. Stanley, and N. Coleman
Interferon-{beta} treatment of cervical keratinocytes naturally infected with human papillomavirus 16 episomes promotes rapid reduction in episome numbers and emergence of latent integrants
Carcinogenesis, November 1, 2006; 27(11): 2341 - 2353.
[Abstract] [Full Text] [PDF]


Home page
J. Clin. Pathol.Home page
S-M A Kulmala, S M Syrjanen, U B Gyllensten, I P Shabalova, N Petrovichev, P Tosi, K J Syrjanen, and B C Johansson
Early integration of high copy HPV16 detectable in women with normal and low grade cervical cytology and histology
J. Clin. Pathol., May 1, 2006; 59(5): 513 - 517.
[Abstract] [Full Text] [PDF]


Home page
Proc. Natl. Acad. Sci. USAHome page
M. R. Pett, M. T. Herdman, R. D. Palmer, G. S. H. Yeo, M. K. Shivji, M. A. Stanley, and N. Coleman
Selection of cervical keratinocytes containing integrated HPV16 associates with episome loss and an endogenous antiviral response
PNAS, March 7, 2006; 103(10): 3822 - 3827.
[Abstract] [Full Text] [PDF]


Home page
BioinformaticsHome page
K. Willbrand, F. Radvanyi, J.-P. Nadal, J.-P. Thiery, and T. M. A. Fink
Identifying genes from up-down properties of microarray expression series
Bioinformatics, October 15, 2005; 21(20): 3859 - 3864.
[Abstract] [Full Text] [PDF]


Home page
J. Virol.Home page
W. G. Hubert
Variant Upstream Regulatory Region Sequences Differentially Regulate Human Papillomavirus Type 16 DNA Replication throughout the Viral Life Cycle
J. Virol., May 15, 2005; 79(10): 5914 - 5922.
[Abstract] [Full Text] [PDF]


Home page
J. Virol.Home page
K. Munger, A. Baldwin, K. M. Edwards, H. Hayakawa, C. L. Nguyen, M. Owens, M. Grace, and K. Huh
Mechanisms of Human Papillomavirus-Induced Oncogenesis
J. Virol., November 1, 2004; 78(21): 11451 - 11460.
[Full Text] [PDF]


Home page
Cancer Res.Home page
M. R. Pett, W. O. F. Alazawi, I. Roberts, S. Dowen, D. I. Smith, M. A. Stanley, and N. Coleman
Acquisition of High-Level Chromosomal Instability Is Associated with Integration of Human Papillomavirus Type 16 in Cervical Keratinocytes
Cancer Res., February 15, 2004; 64(4): 1359 - 1368.
[Abstract] [Full Text] [PDF]


Home page
J. Biol. Chem.Home page
S. Baars, A. Bachmann, A. Levitzki, and F. Rosl
Tyrphostin AG 555 Inhibits Bovine Papillomavirus Transcription by Changing the Ratio between E2 Transactivator/Repressor Function
J. Biol. Chem., September 26, 2003; 278(39): 37306 - 37313.
[Abstract] [Full Text] [PDF]


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 arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Alazawi, W.
Right arrow Articles by Coleman, N.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Alazawi, W.
Right arrow Articles by Coleman, N.


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
Cancer Research Clinical Cancer Research
Cancer Epidemiology Biomarkers & Prevention Molecular Cancer Therapeutics
Molecular Cancer Research Cancer Prevention Research
Cancer Prevention Journals Portal Cancer Reviews Online
Annual Meeting Education Book Meeting Abstracts Online