- Research article
- Open Access
- Open Peer Review
Genomic profiling identifies common HPV-associated chromosomal alterations in squamous cell carcinomas of cervix and head and neck
BMC Medical Genomics volume 2, Article number: 32 (2009)
It is well known that a persistent infection with high-risk human papillomavirus (hrHPV) is causally involved in the development of squamous cell carcinomas of the uterine cervix (CxSCCs) and a subset of SCCs of the head and neck (HNSCCs). The latter differ from hrHPV-negative HNSCCs at the clinical and molecular level.
To determine whether hrHPV-associated SCCs arising from different organs have specific chromosomal alterations in common, we compared genome-wide chromosomal profiles of 10 CxSCCs (all hrHPV-positive) with 12 hrHPV-positive HNSCCs and 30 hrHPV-negative HNSCCs. Potential organ-specific alterations and alterations shared by SCCs in general were investigated as well.
Unsupervised hierarchical clustering resulted in one mainly hrHPV-positive and one mainly hrHPV-negative cluster. Interestingly, loss at 13q and gain at 20q were frequent in HPV-positive carcinomas of both origins, but uncommon in hrHPV-negative HNSCCs, indicating that these alterations are associated with hrHPV-mediated carcinogenesis. Within the group of hrHPV-positive carcinomas, HNSCCs more frequently showed gains of multiple regions at 8q whereas CxSCCs more often showed loss at 17p. Finally, gains at 3q24-29 and losses at 11q22.3-25 were frequent (>50%) in all sample groups.
In this study hrHPV-specific, organ-specific, and pan-SCC chromosomal alterations were identified. The existence of hrHPV-specific alterations in SCCs of different anatomical origin, suggests that these alterations are crucial for hrHPV-mediated carcinogenesis.
In the pathogenesis of uterine cervical cancer the necessary and causal involvement of high-risk types of the human papillomavirus (hrHPV) is widely accepted and supported by strong epidemiological and molecular evidence . HrHPV is present in virtually all cervical carcinomas and the viral oncogenes E6 and E7 are consistently expressed in cervical cancers and precancers. Deregulated expression of these oncogenes in the basal, dividing cells of the epithelium interferes with cell cycle control due to their ability to induce degradation of the tumour suppressor proteins p53 and pRb, respectively. This results in uncontrolled cell proliferation and accumulation of specific (epi)genetic changes in the host cell genome, driving progression to a malignant phenotype [2, 3]. In a previous study we have used array-based comparative genomic hybridisation (array CGH) to determine frequent chromosomal alterations in cervical cancer, which included gains at 1q, 3q and 20q and losses at 8q, 10q, 11q, and 13q . The necessity of these and other additional (epi)genetic alterations in the carcinogenic process is illustrated by the fact that their frequency increases with increasing severity of cervical disease.
On the other hand, head and neck squamous cell carcinoma (HNSCC) is known to be mainly caused by well-established life-style related habits, such as tobacco and excessive alcohol consumption. However, besides the influence of these life-style carcinogens, high risk human papillomavirus (hrHPV) is present in 15–35% of HNSCCs and has been suggested to be a separate aetiological factor in head-and-neck carcinogenesis [5–8]. Several studies have shown that hrHPV-positive HNSCCs are associated with a better clinical outcome [7–9]. Moreover, molecular differences were found between hrHPV-positive and hrHPV-negative HNSCCs, supporting the idea of two separate carcinogenic pathways to HNSCC, one determined by life-style carcinogens and the other by hrHPV [5, 6, 10–13].
In a previous study, using array CGH, we identified a number of chromosomal alterations specific for hrHPV-negative HNSCCs that were absent in hrHPV-positive HNSCCs, including loss at 3p, 5q, and 9p, and amplifications at 11q . The hrHPV-positive HNSCCs were characterised by a lower level of chromosomal alterations, none of which were hrHPV-specific. To investigate potential organ-independent, hrHPV-associated chromosomal alterations, genomic profiles of cervical SCCs (CxSCCs), hrHPV-positive and hrHPV-negative HNSCCs were compared using sophisticated clustering and statistical approaches [4, 13]. In addition, we also investigated the presence of organ-specific alterations and alterations shared by all SCCs included in this study.
We used chromosomal profiles of 30 hrHPV-negative HNSCCs, 12 hrHPV-positive HNSCCs and 10 CxSCCs all of which were previously described (Table 1) [4, 13, 14]. To avoid possible confounding of our results, we excluded a subset of hrHPV-negative HNSCCs described in Smeets et al, which showed little to no alterations and should therefore be considered a separate group .
CGH BAC microarrays produced at the Microarray facility of the VU Medical Center were used. These arrays included the 1 Mb resolution Sanger BAC clone set and a subset of clones from the Children's Hospital Oakland Research Institute (CHORI). Spots were quantified using ImaGene 5.6.1 software (BioDiscovery Ltd, Marina del Rey, CA, USA) with default settings for the flagging of bad quality spots.
The entire dataset described here is available from the Gene Expression Omnibus (GEO, http://www.ncbi.nlm.nih.gov/projects/geo/) through series accession numbers GSE6473 (CxSCC) and GSE12020 (HNSCC).
This study followed the ethical guidelines of the Institutional Review Board of the VU University Medical Center and informed consent was obtained from all patients included.
Array CGH analysis
Calling of gains and losses
BAC clones were positioned along the genome according to the May 2004 freeze. After exclusion of clones with one or more flagged spots, the average of the triplicate spots was calculated for each BAC clone. Log2 ratios were normalised per spotted sub-array by subtraction of the median value of all BAC clones spotted within that sub-array. Segmentation and subsequent calling of gained, amplified and lost regions was done using CGHCall, an automated calling algorithm. Segments with a probability score of ≥ 0.5 were considered gained, amplified or lost .
Reduction of dataset into chromosomal regions
We used the CGHregions algorithm to reduce our dataset to chromosomal regions, accepting maximally 0.1% information loss (Threshold = 0.001). It was shown by Van de Wiel et al that the use of regions instead of single BAC clones improved the effectiveness of subsequent statistical analyses and facilitated interpretation of the results .
The samples were clustered by means of a modified version of WECCA . WECCA is a hierarchical clustering method tailor-made for called aCGH data. The modified version accommodates the use of call probabilities instead of calls. The use of the call probabilities in the unsupervised analysis will give a more subtle picture of the similarities and differences between the samples. The modified version of WECCA defines the distance between two features as the symmetric Kullback-Leibler divergence. The distance between the call probability profiles of two samples is then defined as the average of these divergences over all features. In the construction of the dendrogram we used Ward's linkage as it yields compact and well-separated clusters.
The association between clustering results and HPV status was determined by chi-square testing. The average total number of altered regions was compared between HPV-positive and HPV-negative carcinomas using the non-parametric Mann Whitney test. Two-sided p-values below 0.05 were considered statistically significant. Alteration patterns between HPV-positive (HNSCCs and CxSCCs) and HPV-negative tumours as well as between HNSCCs and CxSCCs were compared using a binomial differential proportion test. The test procedure includes a permutation-based false discovery rate (FDR) correction for multiple testing, needed to discriminate real differences from chance effects . An FDR below 0.10 was considered statistically significant.
Gene ontology analysis
To interpret the biological significance of the genes that are located at altered chromosomal regions of interest, a gene ontology analysis was performed using Ingenuity Pathways Analysis (Ingenuity Systems®, Redwood City, USA). Biological/molecular functions were considered to be significantly overrepresented if they contained more than 1 gene and the Benjamini-Hochberg corrected p-value was p < 0.10 .
hrHPV-positive carcinomas cluster together
To obtain an overview of the similarities between samples, unsupervised hierarchical clustering was performed. This method enabled us to determine in an unbiased manner whether chromosomal profiles of hrHPV-positive HNSCCs were more closely related to hrHPV-negative HNSCCs or hrHPV-positive CxSCCs.
As is shown in Figure 1, two clusters emerged. Cluster 1 contained 24 samples, 18 of which were hrHPV-positive (75%) and 6 were hrHPV-negative. Cluster 2 contained 28 samples of which 4 samples were hrHPV-positive and 24 were hrHPV-negative (86%). This association between cluster assignment and hrHPV status was statistically significant (p < 0.0001). The hrHPV-positive cluster 1 included both samples of cervical (n = 8) and head and neck origin (n = 10), indicating similarities between the chromosomal profiles of hrHPV-induced carcinomas of different anatomical origins. Within cluster 1, however, 7 out of 8 cervical samples formed a separate sub-cluster. This suggests that organ-specific alterations exist as well in hrHPV-positive HNSCCs and CxSCCs.
hrHPV-associated loss at chromosome 13q and gain at 20q
To assess the differences between all three sample groups, the frequency of gains (including amplifications) and losses was analysed for all chromosomal regions (Figure 2). In general, hrHPV-negative carcinomas showed significantly more altered regions than hrHPV-positive carcinomas (p = 0.022).
To determine hrHPV-specific chromosomal alterations, the frequency of alterations was compared between hrHPV-positive (12 HNSCCs and 10 CxSCCs) and hrHPV-negative tumours (30 HNSCCs) for all chromosomal regions. Regions showing a significant difference (FDR < 0.10) in this comparison are shown in Table 2. Interestingly, loss of 13q21.1-21.33 and gain of 20p12.1-q13.33 were significantly more frequent in HPV-positive SCCs compared to hrHPV-negative HNSCCs (66.7% and 58.5% in hrHPV-positive SCCs compared to 33% and 24% in hrHPV-negative SCCs, respectively). As was also shown in our previous study, loss of regions at 3p and 5q, and gains/amplifications of a small region at 11q (CCND1 locus) were specific for hrHPV-negative HNSCCs . A region at 8p showed loss in hrHPV-negative HNSCCs, gain in hrHPV-positive HNSCCs and no alteration in CxSCCs. Gain at chromosome 8q was more frequent in hrHPV-negative HNSCCs, which was due to absence of this alteration in CxSCCs.
Our unsupervised classification results indicated that, within our hrHPV-positive cluster, CxSCCs formed a separate group. To identify potential organ-specific alterations, we therefore also compared the frequency of alterations for all regions between hrHPV-positive HNSCCs (n = 12) and CxSCCs (n = 10) (Table 3). Gains at 3q and losses at 17p were significantly more frequent in CxSCCs than HNSCCs (FDR < 0.10). On the other hand, HNSCCs showed frequent gains at chromosome 8q and losses at 11q. It is important to note that even though gain at 3q and loss at 11q were significantly different between hrHPV-positive HNSCCs and CxSCCs, these alterations were frequent in all sample groups (>50%). Significant differences for these regions can mainly be explained by the fact that the size of the exact altered region differed between groups. Therefore, the smallest regions of overlap between all samples, namely 3q24-29 and 11q22.3-25, may represent general alterations in carcinomas derived from squamous epithelium.
All results described above are summarised in a Venn diagram, showing a general overview of the frequently altered chromosomal arms specific to or common between the (sub) groups (Figure 3).
Biological functions of genes located within hrHPV-specific chromosomal alterations
As described above gain at 20p12.1-q13.33 and loss at 13q21.1-21.33 were significantly more frequent in hrHPV-positive SCCs compared to hrHPV-negative ones. In fact, a 4.5 megabase (Mb) region on chromosome 20 (20q11.21-q11.23), and a 2 Mb region on chromosome 13 (13q21.1) formed the smallest regions of overlap (SRO) at these respective loci when all hrHPV-positive carcinomas were considered (Figure 4A and 4B). The SRO at chromosome 20q contains 78 genes and the one at chromosome 13q contains 6 genes [see Additional file 1]. Within the SRO at chromosome 20, two genes reside, i.e. NCOA6 and RBM39, which showed elevated expression in hrHPV16 E7 expressing cells in vitro . Other cancer-related genes located within this SRO include PIGU, E2F1, and DNMT3B. The SRO on chromosome 13 encompasses the PCDH17 gene and a cluster of five identical loci all of which are predicted to encode proline-rich proteins that contain several dopamine D4 receptor signatures.
Subsequent gene ontology analysis of all genes located in these SROs identified a number of significantly overrepresented GO biological/molecular functions, including cell cycle, cell-to-cell signalling and interaction, cellular growth and proliferation, and a number of cellular maintenance functions (i.e. DNA replication, recombination and repair, gene expression, cellular function and maintenance, cellular assembly and organisation, cellular compromise, cell death, cell morphology, cellular development, nucleic acid metabolism, and lipid metabolism).
In the present study we compared genome-wide chromosomal profiles of hrHPV-negative HNSCCs, hrHPV-positive HNSCCs and hrHPV-positive CxSCCs, to determine whether hrHPV-associated carcinomas of different origins have similar chromosomal signatures. In addition, potential organ-specific alterations were determined within the total group of hrHPV-positive SCCs.
Unsupervised hierarchical clustering resulted in a separate hrHPV-positive cluster, indicating similarities in the chromosomal profiles of hrHPV-induced carcinomas. Subsequent supervised statistical analysis identified a number of hrHPV-associated chromosomal alterations, including gains at 20p13-q13.33 and losses at 13q21.1-21.33, which were frequent (>50%) in hrHPV-positive carcinomas, but were only rarely observed in hrHPV-negative carcinomas. In contrast, lost regions at 3p and 5q as well as amplifications at 11q13.3 (CCND1 locus) were frequent in hrHPV-negative carcinomas, but not in hrHPV-positive carcinomas. In our previous study comparing only hrHPV-positive and hrHPV-negative HNSCCs, the same regions were identified as well as being specific for hrHPV-negative HNSCCs .
In the present study we found that gains at 20q and losses at 13q are specific for hrHPV-positive carcinomas of different anatomical origins. Our previous study, including only hrHPV-positive and hrHPV-negative HNSCCs, did not identify these alterations as specific for hrHPV-positive HNSCCs. However, in the present study a larger number of hrHPV-negative hrHNSCCs and HPV-positive SCCs was included. Furthermore, a sophisticated, objective calling method was presently used to determine gained and lost regions . Interestingly, in a number of studies describing hrHPV E6 and/or E7 mediated immortalisation of keratinocytes of multiple anatomical origins, gains of chromosome 20q and losses of chromosome 13q were consistently observed [20–23]. Low-level gains of chromosome 20q are suggested to be caused by E7 expression and consequent inactivation of the pRb pathway in epithelial cells [20, 21, 24]. Deletion of part of the long arm of chromosome 13 is frequently found in a number of hrHPV-associated tumours, including cervical, anal and head and neck carcinomas [4, 13, 25]. Sabbir et al showed that loss of 13q in HNSCCs was associated with the presence of hrHPV, which is in agreement with our findings .
In a recent study Pyeon et al showed that the gene expression patterns of hrHPV-positive HNSCCs and CxSCCs differed yet shared many changes compared to hrHPV-negative HNSCCs . Our study shows that the same holds true on a chromosomal level. Interestingly, 28% of the genes Pyeon et al found to be differentially expressed between hrHPV-positive and hrHPV-negative carcinomas is located within the chromosomal regions identified in this study and showed expression changes concordant with the chromosomal alterations. Of these genes, 39% was located at chromosome 1p, 25% at 5q, and 14% at 3p. The other genes were located at 11q, 18q and 21q. Only one gene, SYCP2, was located within the hrHPV-specific chromosomal alterations found in this study (20q), but was not located within our SRO. To the best of our knowledge none of the genes are known to directly interact with hrHPVE6 and/or E7 [27, 28]. Pathway analysis of all genes overlapping with our findings identified cell cycle/proliferation as most overrepresented biological function, which is in concordance with the observations made by Pyeon et al .
Pathway analysis of all genes located within the hrHPV-associated SROs at chromosome 20q and 13q found in this study, again underlined the importance of cell cycle (replication and proliferation) related genes in hrHPV-mediated carcinogenesis. This may be related to the continuous E7-regulated E2F1 activation and is accompanied by changes in overall cellular maintenance systems, such as nucleic acid metabolism, as was also found in HPV16 E7 expressing epithelial cells in vitro . Two genes, NCOA6 and RBM39, overlapped between this in vitro study and our results, warranting further investigation of their role in hrHPV-mediated transformation. NCOA6 encodes a transcriptional coactivator interacting with basal transcription factors, histone acetyltransferases, and methyltransferases. RBM39 encodes an RNA binding protein and possible splicing factor and may act as a transcriptional coactivator for the AP-1 transcription activator complex and estrogen receptors. Other cancer-related genes at 20q include E2F1, which is specifically targeted by hrHPV-mediated degradation of pRb, PIGU, which may play a role in cell cycle control and was identified as an oncogene in bladder cancer , and DNMT3B, a de novo DNA methyl transferase. We previously showed that DNMT3B is amplified in the cervical cancer cell line SiHa and found a correlation between increased DNMT3B gene copy numbers and elevated mRNA expression in 78% of CxSCCs . The increased DNMT3B levels are most likely related to the high frequency of tumour suppressor gene promoter hypermethylation events during hrHPV-mediated carcinogenesis [30–32]. The SRO at chromosome 13q includes one known gene, PCDH17, a member of the protocadherin family, which is a subfamily of the cadherin superfamily. At present no reports are available describing interactions between hrHPV and protocadherins. However hrHPV presence has been related to decreased E-cadherin and subsequent impaired immune response [33–35].
Next to the hrHPV-related common events, organ-specific alterations for CxSCCs and HNSCCs were identified as well. CxSCCs more frequently showed loss at 17p, while HNSCCs were characterised by frequent gains at chromosome 8q, harbouring the oncogene c-Myc .
Gains at 3q and losses at 11q were found frequently in all SCCs included in this study, suggesting their involvement in carcinogenesis of squamous epithelial cells in general. Interestingly, we previously showed that gains of 3q were highly frequent in CxSCCs but not in adenocarcinomas of the same organ, further emphasising that this alteration may be specific for squamous epithelium .
Together with results from previous studies, our findings support a causal role for hrHPV in the development of a subset of HNSCCs. Consequently, hrHPV-positive and hrHPV-negative HNSCCs should be regarded as different disease entities requiring different diagnostic and therapeutic approaches. The fact that hrHPV-associated SCCs of different organs have chromosomal alterations in common, suggests that these alterations are crucial for hrHPV-induced carcinogenesis. Diagnostic and/or therapeutic targets based on these alterations may therefore be relevant to hrHPV-associated SCCs of all anatomical origins.
high-risk human papillomavirus
squamous cell carcinoma of the uterine cervix
squamous cell carcinoma of the head and neck
squamous cell carcinoma
comparative genomic hybridisation
bacterial artificial chromosome
false discovery rate
smallest region of overlap
zur Hausen H: Papillomaviruses and cancer: from basic studies to clinical application. Nat Rev Cancer. 2002, 2: 342-350. 10.1038/nrc798.
Steenbergen RD, de Wilde J, Wilting SM, Brink AA, Snijders PJ, Meijer CJ: HPV-mediated transformation of the anogenital tract. J Clin Virol. 2005, 32 (Suppl 1): S25-S33. 10.1016/j.jcv.2004.11.019.
zur Hausen H: Papillomaviruses causing cancer: evasion from host-cell control in early events in carcinogenesis. J Natl Cancer Inst. 2000, 92: 690-698. 10.1093/jnci/92.9.690.
Wilting SM, Snijders PJ, Meijer GA, Ylstra B, Ijssel van den PR, Snijders AM, Albertson DG, Coffa J, Schouten JP, Wiel van de MA, et al: Increased gene copy numbers at chromosome 20q are frequent in both squamous cell carcinomas and adenocarcinomas of the cervix. J Pathol. 2006, 209: 220-230. 10.1002/path.1966.
Braakhuis BJ, Snijders PJ, Keune WJ, Meijer CJ, Ruijter-Schippers HJ, Leemans CR, Brakenhoff RH: Genetic patterns in head and neck cancers that contain or lack transcriptionally active human papillomavirus. J Natl Cancer Inst. 2004, 96: 998-1006.
Dahlgren L, Mellin H, Wangsa D, Heselmeyer-Haddad K, Bjornestal L, Lindholm J, Munck-Wikland E, Auer G, Ried T, Dalianis T: Comparative genomic hybridization analysis of tonsillar cancer reveals a different pattern of genomic imbalances in human papillomavirus-positive and -negative tumors. Int J Cancer. 2003, 107: 244-249. 10.1002/ijc.11371.
Gillison ML, Koch WM, Capone RB, Spafford M, Westra WH, Wu L, Zahurak ML, Daniel RW, Viglione M, Symer DE, et al: Evidence for a causal association between human papillomavirus and a subset of head and neck cancers. J Natl Cancer Inst. 2000, 92: 709-720. 10.1093/jnci/92.9.709.
McKaig RG, Baric RS, Olshan AF: Human papillomavirus and head and neck cancer: epidemiology and molecular biology. Head Neck. 1998, 20: 250-265. 10.1002/(SICI)1097-0347(199805)20:3<250::AID-HED11>3.0.CO;2-O.
Hoffmann M, Gorogh T, Gottschlich S, Lohrey C, Rittgen W, Ambrosch P, Schwarz E, Kahn T: Human papillomaviruses in head and neck cancer: 8 year-survival-analysis of 73 patients. Cancer Lett. 2005, 218: 199-206. 10.1016/j.canlet.2004.09.027.
Pyeon D, Newton MA, Lambert PF, den Boon JA, Sengupta S, Marsit CJ, Woodworth CD, Connor JP, Haugen TH, Smith EM, et al: Fundamental differences in cell cycle deregulation in human papillomavirus-positive and human papillomavirus-negative head/neck and cervical cancers. Cancer Res. 2007, 67: 4605-4619. 10.1158/0008-5472.CAN-06-3619.
Schlecht NF, Burk RD, Adrien L, Dunne A, Kawachi N, Sarta C, Chen Q, Brandwein-Gensler M, Prystowsky MB, Childs G, et al: Gene expression profiles in HPV-infected head and neck cancer. J Pathol. 2007, 213: 283-293. 10.1002/path.2227.
Slebos RJ, Yi Y, Ely K, Carter J, Evjen A, Zhang X, Shyr Y, Murphy BM, Cmelak AJ, Burkey BB, et al: Gene expression differences associated with human papillomavirus status in head and neck squamous cell carcinoma. Clin Cancer Res. 2006, 12: 701-709. 10.1158/1078-0432.CCR-05-2017.
Smeets SJ, Braakhuis BJ, Abbas S, Snijders PJ, Ylstra B, Wiel van de MA, Meijer GA, Leemans CR, Brakenhoff RH: Genome-wide DNA copy number alterations in head and neck squamous cell carcinomas with or without oncogene-expressing human papillomavirus. Oncogene. 2006, 25: 2558-2564. 10.1038/sj.onc.1209275.
Smeets SJ, Brakenhoff RH, Ylstra B, van Wieringen WN, Wiel van de MA, Leemans CR, Braakhuis BJ: Genetic classification of oral and oropharyngeal carcinomas identifies subgroups with a different prognosis. Cell Oncol. 2009.
Wiel van de MA, Kim KI, Vosse SJ, van Wieringen WN, Wilting SM, Ylstra B: CGHcall: calling aberrations for array CGH tumor profiles. Bioinformatics. 2007, 23: 892-894. 10.1093/bioinformatics/btm030.
Wiel van de MA, van Wieringen WN: CGHRegions: Dimension Reduction for Array CGH Data with Minimal Information Loss. Cancer Informatics. 2007, 2: 55-63.
van Wieringen WN, Wiel van de MA, Ylstra B: Weighted clustering of called array CGH data. Biostatistics. 2008, 9: 484-500. 10.1093/biostatistics/kxm048.
Manduchi E, Grant GR, McKenzie SE, Overton GC, Surrey S, Stoeckert CJ: Generation of patterns from gene expression data by assigning confidence to differentially expressed genes. Bioinformatics. 2000, 16: 685-698. 10.1093/bioinformatics/16.8.685.
Benjamini Y, Drai D, Elmer G, Kafkafi N, Golani I: Controlling the false discovery rate in behavior genetics research. Behav Brain Res. 2001, 125: 279-284. 10.1016/S0166-4328(01)00297-2.
Klingelhutz AJ, Qian Q, Phillips SL, Gourronc FA, Darbro BW, Patil SR: Amplification of the chromosome 20q region is associated with expression of HPV-16 E7 in human airway and anogenital epithelial cells. Virology. 2005, 340: 237-244. 10.1016/j.virol.2005.06.027.
Jin Y, Zhang H, Tsao SW, Jin C, Lv M, Strombeck B, Wiegant J, Wan TS, Yuen PW, Kwong YL: Cytogenetic and molecular genetic characterization of immortalized human ovarian surface epithelial cell lines: consistent loss of chromosome 13 and amplification of chromosome 20. Gynecol Oncol. 2004, 92: 183-191. 10.1016/j.ygyno.2003.09.007.
Savelieva E, Belair CD, Newton MA, DeVries S, Gray JW, Waldman F, Reznikoff CA: 20q gain associates with immortalization: 20q13.2 amplification correlates with genome instability in human papillomavirus 16 E7 transformed human uroepithelial cells. Oncogene. 1997, 14: 551-560. 10.1038/sj.onc.1200868.
Steenbergen RD, Hermsen MA, Walboomers JM, Meijer GA, Baak JP, Meijer CJ, Snijders PJ: Non-random allelic losses at 3p, 11p and 13q during HPV-mediated immortalization and concomitant loss of terminal differentiation of human keratinocytes. Int J Cancer. 1998, 76: 412-417. 10.1002/(SICI)1097-0215(19980504)76:3<412::AID-IJC20>3.0.CO;2-B.
Cuthill S, Agarwal P, Sarkar S, Savelieva E, Reznikoff CA: Dominant genetic alterations in immortalization: role for 20q gain. Genes Chromosomes Cancer. 1999, 26: 304-311. 10.1002/(SICI)1098-2264(199912)26:4<304::AID-GCC4>3.0.CO;2-1.
Heselmeyer K, du MS, Blegen H, Friberg B, Svensson C, Schrock E, Veldman T, Shah K, Auer G, Ried T: A recurrent pattern of chromosomal aberrations and immunophenotypic appearance defines anal squamous cell carcinomas. Br J Cancer. 1997, 76: 1271-1278.
Sabbir MG, Roy A, Mandal S, Dam A, Roychoudhury S, Panda CK: Deletion mapping of chromosome 13q in head and neck squamous cell carcinoma in Indian patients: correlation with prognosis of the tumour. Int J Exp Pathol. 2006, 87: 151-161. 10.1111/j.0959-9673.2006.00467.x.
Laughlin-Drubin ME, Munger K: The human papillomavirus E7 oncoprotein. Virology. 2009, 384: 335-344. 10.1016/j.virol.2008.10.006.
Wise-Draper TM, Wells SI: Papillomavirus E6 and E7 proteins and their cellular targets. Front Biosci. 2008, 13: 1003-1017. 10.2741/2739.
Guo Z, Linn JF, Wu G, Anzick SL, Eisenberger CF, Halachmi S, Cohen Y, Fomenkov A, Hoque MO, Okami K, et al: CDC91L1 (PIG-U) is a newly discovered oncogene in human bladder cancer. Nat Med. 2004, 10: 374-381. 10.1038/nm1010.
Henken FE, Wilting SM, Overmeer RM, van Rietschoten JG, Nygren AO, Errami A, Schouten JP, Meijer CJ, Snijders PJ, Steenbergen RD: Sequential gene promoter methylation during HPV-induced cervical carcinogenesis. Br J Cancer. 2007, 97: 1457-1464. 10.1038/sj.bjc.6604055.
Wentzensen N, Sherman ME, Schiffman M, Wang SS: Utility of methylation markers in cervical cancer early detection: appraisal of the state-of-the-science. Gynecol Oncol. 2009, 112: 293-299. 10.1016/j.ygyno.2008.10.012.
Whiteside MA, Siegel EM, Unger ER: Human papillomavirus and molecular considerations for cancer risk. Cancer. 2008, 113: 2981-2994.
Al Moustafa AE, Kassab A, Darnel A, Yasmeen A: High-risk HPV/ErbB-2 interaction on E-cadherin/catenin regulation in human carcinogenesis. Curr Pharm Des. 2008, 14: 2159-2172. 10.2174/138161208785740216.
Caberg JH, Hubert PM, Begon DY, Herfs MF, Roncarati PJ, Boniver JJ, Delvenne PO: Silencing of E7 oncogene restores functional E-cadherin expression in human papillomavirus 16-transformed keratinocytes. Carcinogenesis. 2008, 29: 1441-1447. 10.1093/carcin/bgn145.
Hubert P, Caberg JH, Gilles C, Bousarghin L, Franzen-Detrooz E, Boniver J, Delvenne P: E-cadherin-dependent adhesion of dendritic and Langerhans cells to keratinocytes is defective in cervical human papillomavirus-associated (pre)neoplastic lesions. J Pathol. 2005, 206: 346-355. 10.1002/path.1771.
Issing WJ, Wustrow TP, Heppt WJ: Oncogenes related to head and neck cancer. Anticancer Res. 1993, 13: 2541-2551.
The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1755-8794/2/32/prepub
We would like to thank the Mapping Core and Map Finishing groups of the Wellcome Trust Sanger Institute for initial clone supply and verification. This work was partly supported by the Centre for Medical Systems Biology (CMSB) in the framework of the Netherlands Genomic Initiative (NGI).
The authors declare that they have no competing interests.
SMW and SJS performed array experiments and data analysis and drafted the manuscript. RDMS and BJMB participated in the design of the study and the drafting of the manuscript. WNvW and MAvdW contributed to the statistical analyses. BY and GAM were involved in the development of the micorarray platform and provided the facilities for the microarray experiments. PJFS, CRL, CJLMM and RHB contributed to the conception of the study and critically revised the manuscript. All authors read and approved the final manuscript.
Saskia M Wilting, Serge J Smeets, Boudewijn JM Braakhuis and Renske DM Steenbergen contributed equally to this work.