
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
Imaging, Diagnosis, Prognosis |
Authors' Affiliations: Departments of 1 Pathology and 2 Gynecology, Leiden University Medical Center, Leiden, the Netherlands
Requests for reprints: Marjon A. de Boer, Department of Gynecology, Leiden University Medical Center, Albinusdreef 2, Postbus 9600, 2300 RC Leiden, the Netherlands. Phone: 31-71-5266596; Fax: 31-71-5248158; E-mail: M.A.de_boer{at}LUMC.nl.
| Abstract |
|---|
|
|
|---|
Experimental Design: We studied the number of viral DNA copies and the level of HPV E6/E7 mRNA expression in 75 HPV 16positive or HPV 18positive International Federation of Gynecology and Obstetrics stage Ib and IIa cervical cancer patients. Measurements were done with quantitative PCR. DNA copy number analysis was done on pure tumor cell samples enriched with flow sorting. mRNA expression data were compensated for the percentage of tumor cells included.
Results: The number of viral DNA copies was not predictive of survival in cervical cancer patients. In contrast, high HPV E6/E7 mRNA expression was strongly related to an unfavorable prognosis (P = 0.006). In a multivariate Cox model for overall survival, including all known prognostic variables and stratified for HPV type, the level of E6/E7 mRNA expression was an independent prognostic indicator, second only to lymph node status. No correlation was observed between DNA copy number and the level of HPV E6/E7 mRNA expression, which reflects that not all DNA copies are equally transcriptionally active.
Conclusions: Cervical cancer patients with high HPV E6/E7 oncogene mRNA expression have a worse survival independently from established prognostic factors.
HPV expresses two viral oncoproteins: E6 and E7. These proteins bind to and inactivate the tumor suppressor proteins p53 and pRb, respectively, causing deregulation of the cell cycle (3). The E6 and E7 proteins are transcribed from one shared transcript. In addition to the full-length E6/E7 mRNA, high-risk HPV types generate spliced transcripts referred to as E6*. For HPV 16, a second spliced transcript called E6*II is also generated. The function of the spliced transcripts E6*I and E6*II is not well understood. Previously, it was thought that the E7 protein could only be transcribed from the spliced mRNA; however, more recent studies provided evidence that E7 is synthesized from both the spliced and the full-length transcripts (4, 5).
The number of HPV DNA copies per cell, viral load, was suggested to correlate with disease stage, showing increasing copy numbers from mild dysplasia to cervical cancer (610). However, results are conflicting because of the variation in sampling techniques and different methods used to calculate viral load (2). Similarly, the presence of transcripts of the E6 and E7 genes was claimed to correlate with severity of cervical intraepithelial neoplasia (4, 1113). Most of these studies were done nonquantitatively, and most importantly, they did not adjust for the percentage of HPV-infected cells in the total cell mass examined, including normal epithelial cells and stromal cells as well. To overcome this final problem, we designed a method to study tumor cells exclusively. Samples were analyzed by flow cytometry to determine the percentage of tumor cells used in mRNA expression analysis. Subsequently, HPV DNA copy number measurements were done on flow-sorted tumor cells as described previously (14). Both DNA and mRNA measurements were done with quantitative PCR, and results were normalized to the expression of housekeeping genes. We assess the prognostic value of DNA copy number and the level of full-length and spliced E6/E7 mRNA expression in HPV 16infected or HPV 18infected cervical cancer patients.
| Materials and Methods |
|---|
|
|
|---|
Sample preparation. H&E-stained slides were reviewed, and formalin-fixed, paraffin-embedded tissue blocks were macrodissected to remove normal and precancerous epithelial tissue. Subsequently, 60-µm sections were cut for DNA copy number and mRNA expression analysis. Flow cytometry was done using triple staining with keratin (to identify epithelial tumor cells), vimentin (to identify stromal cells), and propidium iodide (for DNA staining; ref. 14). The percentage of tumor cells was calculated followed by separate flow sorting of keratin- and vimentin-positive cell populations.
HPV DNA copy number. Keratin-positive, exclusively tumor cells were used for HPV DNA measurement. Vimentin-positive cells were used as a negative control. Cells were centrifuged and resuspended at a concentration of 1,000/µL. DNA was extracted by incubation with proteinase K (3 mg/mL) overnight at 56°C.
A HPV type-specific quantitative Taqman PCR was done to assess the number of HPV DNA copies using primers for the E7 gene (16). In addition, we amplified HBG2 (hemoglobin ß chain, Genbank accession no. U01317) and DCC (deleted in colon cancer, Genbank accession no. NM005215) genes to allow an accurate compensation for the total DNA input and quality. Amplicons were detected with fluorescent Taqman probes (Eurogentec, Seraing, Belgium). Primer and probe sequences are listed in Table 1
. We added DNA from
200 cells to the PCR mix (qPCR Core kit, Eurogentec) and did reactions in a final volume of 25 µL. The following thermocycler profile was used: 95°C for 2 min; 40 cycles of 15 s at 95°C, 60 s at 60°C. To calculate the viral DNA copy number, we included cell lines established from cervical cancer with known viral copy numbers (SiHa, 2 copies of HPV 16; HeLa, 40 copies of HPV 18; ref. 17) that were fixed in formalin, embedded in paraffin, and prepared identically to patient samples. Standard curves were prepared from cell line DNA. The absolute viral load was calculated as follows: viral load = F x (SQHPV / SQcontrol genes), where SQHPV is the starting quantity of HPV DNA, SQcontrol genes is the starting quantity of the two control genes (geometric mean of HBG2 and DCC starting quantity), and F is the replication factor. The replication factor was calculated from the known viral copy number of the cell lines.
|
Primers were designed using Beacon Designer 3 software (Premier Biosoft International, Palo Alto, CA) with amplicons as small as possible (80-120 bp). Intron-spanning primers were used for the spliced transcripts, and known polymorphic positions were avoided. PCR products were sequenced to confirm that the target gene was amplified. HPV primer sequences are described in Table 2 . The primer sequences of the normalization genes were described previously (18). We used 0.2 µL cDNA template in a total volume of 25 µL containing 0.3 µmol/L of each forward and reverse primer and the components of the qPCR Core kit for SYBR Green I (Eurogentec). The same cycle profile as in the DNA analysis was used. SYBR Green fluorescence was measured at the end of the elongation phase, and a melting curve was generated after the last amplification cycle.
|
600 DNA copies of HPV 16, was also done. To ensure contaminating DNA had been removed, a control PCR was done before reverse transcription (no-reverse transcriptase control) using a nonintron-spanning primer set. We measured the expression of four normalization genes (UBC, GAPD, CFL1, and ACTB; ref. 18). The most stably expressed genes were identified by using "NormFinder." UBC and ACTB proved most stable and were used for the normalization of expression data (18). The following equation was used to calculate the expression of the target genes: ExpHPV = SQHPV / mean (SQACTB, SQUBC) x % tumor cells, where ExpHPV is the relative expression of the HPV target transcript and SQ is the starting quantities of a transcript as calculated from the standard curves. The percentage of tumor cells were calculated from flow cytometry analysis. The equation was based on the assumptions that HPV is expressed only in tumor cells and that the expression of normalization genes is similar in tumor and nontumor cells. Quantitative real-time PCRs were done with an iCycler iQ and analyzed using Optical System Software version 3.0a (Bio-Rad, Richmond, CA). Each measurement was done in duplicate. Statistical analysis. Statistical analysis was done using the Statistical Package for the Social Sciences 10.0.7 software package. Real-time data were log transformed before applying statistics. Differences between groups were evaluated with the Student's t test. Correlation analysis was done using Pearson's correlation. Univariate analysis of overall 3-year survival was done using the Kaplan-Meier method and the log-rank test. The Cox proportional hazard model was used for multivariate survival analysis. All tests were two sided and the significance level was set to 5%, corresponding to 95% confidence intervals.
| Results |
|---|
|
|
|---|
|
mRNA expression. HPV E6/E7 transcripts were identified in 72 of 75 HPV 16positive tumors. In one case where transcripts were not detected, the percentage of tumor cells (5%) were low, and in the other two, a low amount of mRNA was isolated. These three samples were excluded from further analysis of mRNA levels. Importantly, all control reactions lacking reverse transcriptase were negative. In HPV 16positive tumors, the full-length transcript was identified in 50 of 52 cases, whereas the spliced transcripts E6*I and E6*II were identified in all cases. In HPV 18positive tumors, we detected the full-length transcripts in all samples and the HPV 18 E6* transcript in 18 of 20 samples.
The expression levels of all three HPV 16 E6/E7 transcripts were significantly correlated: full-length E6 versus E6*I, R = 0.418, P = 0.002 (Fig. 1A ); full-length E6 versus E6*II, R = 0.438, P = 0.001 (Fig. 1B); and E6*I versus E6*II, R = 0.652, P < 0.0001 (Fig. 1C). Similarly, the expression levels of HPV 18 full-length transcripts and E6* transcripts were correlated: R = 0.623, P = 0.003 (Fig. 1D). mRNA expression was not associated with the age of the samples, FIGO stage, histopathology, lymph node status, tumor size, or depth of infiltration.
|
|
|
|
| Discussion |
|---|
|
|
|---|
In contrast to previous studies, we restricted both DNA and mRNA analyses to tumor cells. DNA copy number analysis was limited to tumor cells by flow sorting, and the quantitative PCR was normalized for total DNA input. In mRNA quantitative PCR analysis, we compensated for the percentage of tumor cells estimated by flow cytometry. Several studies described an association between viral DNA copy number and the level of mRNA expression, but no compensation for the percentage of tumor cells was done (2123). The amount of HPV-infected cells included in the analysis might have largely influenced the detected levels of both DNA and mRNA. In addition, many studies on precancerous lesions, which were done on cervical exfoliated cells, described a correlation between the quantity of viral DNA (69) or viral mRNA (16, 24) and the severity of cervical dysplasia. Those studies did not compensate for the amount of HPV-infected dysplastic cells, and results could well be the effect of lesion size as suggested previously (25). We showed in the present study that the contamination of stromal and normal epithelial cells can be substantial, with a percentage of tumor cells ranging from 2% to 80%. The same effect might be present in precancerous lesions. Therefore, viral DNA load and viral mRNA measurements should be compensated for the percentage of HPV-infected cells included in the analyzed sample, or results should be interpreted with caution.
Apart from restricting the analysis to tumor cells, several other measures were taken to ensure validity. We treated the mRNA samples with DNase to make sure genomic DNA was not being amplified with the transcripts of interest, and a control reaction lacking reverse transcriptase ensured that contaminating DNA was not present. We tested internal control genes, and the two most stably expressed genes were used for normalization of the data. Outcome measures were related to cell lines that were treated the same as patient samples. We used formalin-fixed, paraffin-embedded tumor tissue in which mRNA is partly fragmented, yet the relative gene expression is equivalent to fresh tissue (26). In previous studies, HPV oncogene transcripts could not be detected in all tumors (2729). We detected all three E6/E7 transcripts in almost every sample, indicating that previous methods might have been insufficiently sensitive. Some specific patterns of expression of the different transcripts have been described to relate to clinical variables of cervical cancer (27, 28). In our study, the levels of expression of the three different HPV oncogene transcripts were well correlated, and no specific expression patterns were noticed.
The finding that patients with a high level of HPV E6/E7 mRNA expression have a shorter survival can be explained by the cell cycle deregulating actions of E6 and E7, which bind to and inactivate the tumor suppressor proteins p53 and pRb, respectively, and cause proliferative cell growth (3). Other HPV-induced effects could contribute to the worse survival as well. The process of tumor invasion, which requires changes in epithelial cell-stroma interaction, is influenced by HPV through the induction of vascular endothelial growth factor, permitting the tumor to acquire blood supply (30). In addition, the inhibition of transforming growth factor-ß is influenced by HPV, resulting in increased metastatic potential (31).
The fact that not all DNA copies are transcriptionally equally active might also be a way of immune escape by producing as little antigen as possible. Yet, the relation between the quantity of HPV mRNA expression and immune response has never been studied. Still, from the present study, it can be concluded that, even if antitumor immune response is HPV antigen dose dependent, the effect on the immune response is weaker than the proliferative effect of HPV oncogene expression on tumor growth. The association that we found between HPV oncogene expression and prognosis of cervical cancer patients suggests that a therapy directed against E6 and E7 mRNA expression might be effective. Although none of these molecular therapies have been clinically tested, several successful in vitro studies have been done (3).
In conclusion, we found that a high level of HPV 16 and HPV 18 E6/E7 mRNA expression was an independent predictor of unfavorable prognosis in cervical cancer, whereas the number of HPV DNA copies per cell was not. Importantly, when the analysis was limited to lymph nodenegative patients, HPV E6/E7 mRNA expression was still of significant prognostic value. Our study underscores the importance of HPV E6/E7 mRNA expression and indicates that, apart from an essential role in carcinogenesis, HPV plays a critical role in the progression of cervical cancer.
| Acknowledgments |
|---|
| 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.
Received 6/28/06; revised 9/ 8/06; accepted 10/17/06.
| References |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
S. Bellone, K. El-Sahwi, E. Cocco, F. Casagrande, M. Cargnelutti, M. Palmieri, E. Bignotti, C. Romani, D.-A. Silasi, M. Azodi, et al. Human Papillomavirus Type 16 (HPV-16) Virus-Like Particle L1-Specific CD8+ Cytotoxic T Lymphocytes (CTLs) Are Equally Effective as E7-Specific CD8+ CTLs in Killing Autologous HPV-16-Positive Tumor Cells in Cervical Cancer Patients: Implications for L1 Dendritic Cell-Based Therapeutic Vaccines J. Virol., July 1, 2009; 83(13): 6779 - 6789. [Abstract] [Full Text] [PDF] |
||||
![]() |
W. E. Corver, A. Middeldorp, N. T. ter Haar, E. S. Jordanova, M. van Puijenbroek, R. van Eijk, C. J. Cornelisse, G. J. Fleuren, H. Morreau, J. Oosting, et al. Genome-wide Allelic State Analysis on Flow-Sorted Tumor Fractions Provides an Accurate Measure of Chromosomal Aberrations Cancer Res., December 15, 2008; 68(24): 10333 - 10340. [Abstract] [Full Text] [PDF] |
||||
![]() |
P. E. Castle, C. Meyers, S. Alam, and M. J. Conway How Does Tobacco Smoke Contribute to Cervical Carcinogenesis? J. Virol., June 15, 2008; 82(12): 6084 - 6086. [Full Text] [PDF] |
||||
![]() |
I. Kraus, C. Driesch, S. Vinokurova, E. Hovig, A. Schneider, M. von Knebel Doeberitz, and M. Durst The Majority of Viral-Cellular Fusion Transcripts in Cervical Carcinomas Cotranscribe Cellular Sequences of Known or Predicted Genes Cancer Res., April 1, 2008; 68(7): 2514 - 2522. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| 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 |