Skip to main content

Genomic profiling identifies common HPV-associated chromosomal alterations in squamous cell carcinomas of cervix and head and neck



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.

Peer Review reports


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 [1]. 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 [4]. 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 [58]. Several studies have shown that hrHPV-positive HNSCCs are associated with a better clinical outcome [79]. 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, 1013].

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 [13]. 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.


CGH microarrays

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 [14].

Table 1 Summary of clinical data of carcinomas included in this study.

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, 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 [15].

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 [16].

Clustering analysis

The samples were clustered by means of a modified version of WECCA [17]. 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.

Statistical analysis

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 [18]. 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 [19].


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.

Figure 1
figure 1

Result from the unsupervised hierarchical clustering analysis. Cluster 1 contains in majority hrHPV-positive carcinomas, as is indicated by the black boxes in the legend underneath the heatmap (p < 0.0001).

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).

Figure 2
figure 2

Frequency plots for all 3 sample groups. The frequency of gains (positive axis) and losses (negative axis) are shown for A. hrHPV-negative HNSCCs, B. hrHPV-positive HNSCCs and C. CxSCCs for chromosome 1–22.

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 [13]. 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.

Table 2 Significantly different chromosomal alterations between hrHPV+ and hrHPV- carcinomas.

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.

Table 3 Significantly different chromosomal alterations between hrHPV-positive HNSCCs and CxSCCs.

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).

Figure 3
figure 3

Summary of common and specific chromosomal alterations in the different sample groups. Chromosome arms showing frequent alterations (>50%) in one or more sample groups (hrHPV-positive HNSCCs; hrHPV-negative HNSCCs; CxSCCs) are shown in a Venn diagram. ↑ indicates gain; ↓ indicates loss

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 [20]. 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.

Figure 4
figure 4

Genomic coordinates of A. losses at chromosome 13 and B. gains at chromosome 20 are shown for all hrHPV-positive carcinomas. Chromosomal alterations in CxSCCs are shown by dashed lines and alterations in hrHPV-positive HNSCCs by solid lines. In C. the smallest regions of overlap (SROs) between hrHPV-positive carcinomas at chromosome 13 and 20 are summarised.

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 [13].

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 [15]. 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 [2023]. 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 [26].

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 [10]. 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 [10].

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 [20]. 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 [29], 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 [4]. The increased DNMT3B levels are most likely related to the high frequency of tumour suppressor gene promoter hypermethylation events during hrHPV-mediated carcinogenesis [3032]. 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 [3335].

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 [36].

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 [4].


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


gene ontology.


  1. zur Hausen H: Papillomaviruses and cancer: from basic studies to clinical application. Nat Rev Cancer. 2002, 2: 342-350. 10.1038/nrc798.

    Article  CAS  PubMed  Google Scholar 

  2. 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.

    Article  CAS  PubMed  Google Scholar 

  3. 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.

    Article  CAS  PubMed  Google Scholar 

  4. 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.

    Article  CAS  PubMed  Google Scholar 

  5. 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.

    Article  CAS  PubMed  Google Scholar 

  6. 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.

    Article  CAS  PubMed  Google Scholar 

  7. 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.

    Article  CAS  PubMed  Google Scholar 

  8. 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.

    Article  CAS  PubMed  Google Scholar 

  9. 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.

    Article  CAS  PubMed  Google Scholar 

  10. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. 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.

    Article  CAS  PubMed  Google Scholar 

  12. 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.

    Article  CAS  PubMed  Google Scholar 

  13. 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.

    Article  CAS  PubMed  Google Scholar 

  14. 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.

    Google Scholar 

  15. 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.

    Article  PubMed  Google Scholar 

  16. Wiel van de MA, van Wieringen WN: CGHRegions: Dimension Reduction for Array CGH Data with Minimal Information Loss. Cancer Informatics. 2007, 2: 55-63.

    Google Scholar 

  17. 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.

    Article  PubMed  Google Scholar 

  18. 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.

    Article  CAS  PubMed  Google Scholar 

  19. 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.

    Article  CAS  PubMed  Google Scholar 

  20. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. 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.

    Article  CAS  PubMed  Google Scholar 

  22. 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.

    Article  CAS  PubMed  Google Scholar 

  23. 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.

    Article  CAS  PubMed  Google Scholar 

  24. 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.

    Article  CAS  PubMed  Google Scholar 

  25. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Laughlin-Drubin ME, Munger K: The human papillomavirus E7 oncoprotein. Virology. 2009, 384: 335-344. 10.1016/j.virol.2008.10.006.

    Article  Google Scholar 

  28. Wise-Draper TM, Wells SI: Papillomavirus E6 and E7 proteins and their cellular targets. Front Biosci. 2008, 13: 1003-1017. 10.2741/2739.

    Article  CAS  PubMed  Google Scholar 

  29. 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.

    Article  CAS  PubMed  Google Scholar 

  30. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. 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.

    Article  CAS  PubMed  Google Scholar 

  32. Whiteside MA, Siegel EM, Unger ER: Human papillomavirus and molecular considerations for cancer risk. Cancer. 2008, 113: 2981-2994.

    Article  CAS  PubMed  Google Scholar 

  33. 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.

    Article  CAS  PubMed  Google Scholar 

  34. 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.

    Article  CAS  PubMed  Google Scholar 

  35. 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.

    Article  CAS  PubMed  Google Scholar 

  36. Issing WJ, Wustrow TP, Heppt WJ: Oncogenes related to head and neck cancer. Anticancer Res. 1993, 13: 2541-2551.

    CAS  PubMed  Google Scholar 

Pre-publication history

Download references


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).

Author information

Authors and Affiliations


Corresponding author

Correspondence to Renske DM Steenbergen.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

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.

Electronic supplementary material


Additional file 1: BAC clones and genes included in the SROs of the hrHPV-specific chromosomal alterations at chromosomes 13q and 20q. All BAC clones and genes located within the identified SROs at chromosome 20q and 13q are listed here. (DOC 26 KB)

Authors’ original submitted files for images

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Wilting, S.M., Smeets, S.J., Snijders, P.J. et al. Genomic profiling identifies common HPV-associated chromosomal alterations in squamous cell carcinomas of cervix and head and neck. BMC Med Genomics 2, 32 (2009).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: