Skip to main content

High-risk HPV infection-associated hypermethylated genes in oropharyngeal squamous cell carcinomas



HPV-positive oropharyngeal squamous cell carcinomas (OPSCCs) are sensitive to chemo-radiation therapy and have favorable survival outcomes compared with HPV-negative cancers. These tumors are usually not related to tobacco and alcohol exposure. Therefore, diagnosing HPV-positive OPSCCs for the appropriate disease management is crucial, and no suitable markers are available for detecting early malignancies in HPV-infected tissues. In this study, we attempt to find HPV-specific epigenetic biomarkers for OPSCCs.


A total of 127 surgical samples were analyzed for HPV positivity and promoter methylation of a panel of genes. HPV detection was performed by PCR detection of HPV E6 and E7 viral oncoproteins. In addition, promoter methylation of a total of 8 genes (DAPK, FHIT, RASSF1A, TIMP3, AGTR1, CSGALNACT2, GULP1 and VGF) was analyzed by quantitative-methylation specific PCR (QMSP), and their associations with HPV positivity or RB/p16 expressions were evaluated.


AGTR1 and FHIT were frequently methylated in HPV-positive OPSCC samples with a good area under the curve (AUC over 0.70). In addition, these genes' promoter methylation was significantly associated with p16 positive and RB negative cases, which were the characteristics of OPSCC cases with favorable survival outcomes. Either AGTR1 or FHIT methylated cases were significantly associated with HPV-positive cancers with 92.0% sensitivity (P < 0.001). Also, they had significantly better overall survival (P = 0.047) than both unmethylated cases.


A combination of AGTR1 and FHIT methylation demonstrated a suitable detection marker of OPSCCs derived from the HPV-infected field, familiar with p16-positive and RB-negative phenotypes.

Peer Review reports


Oral Human papillomavirus (HPV) infection is strongly associated with oropharyngeal squamous cell carcinoma (OPSCC) [1, 2]. About 40–80% of oropharyngeal cancers (OCs) are caused by HPV in the USA, whereas the prevalence of HPV infection-associated OCs varies from 20 to 90% in Europe [3]. The incidence of OPSCC typically occurs at tonsil and tongue base [3]. Carcinogenesis of HPV-positive OCs is characterized by TP53-degradation, retinoblastoma (RB) pathway-inactivation and p16-upregulation [3]. Briefly, HPV infection dysregulates the oncoprotein expressions of E6 and E7. E6 protein binds to P53 and promotes its degradation [4, 5]. E7 protein binds and inactivates pRB releasing E2F from the pRB-E2F complex, and E2F activates the cell cycle. Disruption of the pRB-E2F complex also releases p16 from its negative feedback control, resulting in p16 overexpression. Since pRB is degraded and E2F is unbound, p16 no longer has an inhibitory function on pRB and cannot inhibit the cell cycle. By contrast, the carcinogenesis process of HPV-negative OCs is different. Typically, tobacco-related OCs are characterized by TP53 mutation and down-regulation of p16. Thus, HPV-positive and HPV-negative OCs have obviously different kinds of carcinogenesis pathways.

Epigenetic alterations of critical cancer-associated genes are considered as a hallmark of cancer. Van Kempen PM et al. reported a systematic review of differential methylation profiles between HPV-positive and HPV-negative OPSCC [6] in 2014. They summarized differential methylated genes between HPV- positive and HPV-negative cases. Accumulated evidence suggests that promoter methylation more frequently occurs in HPV-positive tumors than HPV-negative tumors. However, no specific panel of methylation markers has yet been identified for potential clinical use.

The primary purpose of this study was to identify differentially methylated candidate genes between HPV-positive OPSCCs and others. They could be used as markers for non-invasive early molecular detection approaches (screening) using bodily fluids (such as saliva and blood) and non-invasive monitoring of recurrent disease. The secondary purpose was to find the correlation of methylation markers with representative HPV-related carcinogenesis pathway factors, p16 and RB.

Materials and methods

Clinical sample collection

A total of 94 patients undergoing surgical resection of primary untreated OPSCCs at The Johns Hopkins Hospital from 1997 to 2008 were included. Written informed consent was obtained from all subjects and their legal guardians recruited under this protocol before participation in the study. All samples were obtained as anonymized materials following the Declaration of Helsinki. They were classified into 50 HPV-positive and 44 HPV-negative cases using the PCR detection method focusing on HPV E6 and E7 regions of HPV type 16 following the previous paper from our lab [7]. In addition, we also used 33 non-neoplastic tonsil tissues from other non-cancer patients cohort as a control. Details of these samples are available in Table 1.

Table 1 Clinicopathological characteristics of study cohort

DNA extraction and bisulfite modification

Hematoxylin and eosin-stained sections were histologically examined to evaluate the presence of tumor cells. Sample sections showing more than 70% of tumor cells were used for DNA extraction. Microdissected tissues were digested with 1% SDS and 50 μg/ml proteinase K (Boehringer Mannheim, Germany) at 48 °C overnight. Phenol and chloroform extraction and ethanol precipitation of DNA were performed as previously described [8]. Then, sodium bisulfite-mediated conversion of unmethylated cytosines in DNA was performed by EpiTect Bisulfite Kit (Qiagen, Venlo, Netherlands).

Candidate genes

A total of 8 genes were selected for promoter methylation analysis. Among these eight genes, four genes are highly methylated in head and neck squamous cell carcinomas (HNSCCs) as others, and we reported previously [DAPK [9], FHIT [10], RASSF1A [11, 12] and TIMP3 [13, 14]]. The remaining four novel methylation markers (AGTR1 [15], CSGALNACT2, GULP1 [16, 17] and VGF [18, 19]) were chosen from our recent high throughput “methylome” study based on their cancer relevance from microarray data analyses [16, 18, 20]. Chromosomal locus, proposed function and known association with cancer of these genes were summarized in Supplementary Table S1.

Quantitative methylation-specific PCR (QMSP)

Bisulfite-modified genomic DNA samples served as templates for QMSP. Primers and probes were listed in Supplementary Table S2. Serial dilutions (90–0.009 ng) of CpG methyltransferase (New England BioLabs, Ipswich, MA) were used to construct a calibration curve for each plate. Also, negative controls (human leukocyte genomic DNA from a healthy donor) and multiple water blanks were placed [21]. Amplification reactions were carried out in triplicate in a final volume of 20 μl containing 3 μl bisulfite-modified DNA, 600 nmol/l forward and reverse primers, 200 nmol/l probes, 0.6 unit Platinum Taq D (Invitrogen), dATP, dCTP, dGTP, and dTTP in a concentration of 200 μmol/l each and 6.7 mmol/l MgCl2. Amplification reactions were carried out in 384-well plates in a 7900HT Fast Real-Time PCR System (Life Technologies, Carlsbad, CA) and were analyzed by the Sequence Detector System software (SDS 2.3; Life technologies). The relative level of methylated DNA for each gene in each sample was determined as a ratio of the QMSP value of the amplified gene to ACTB, multiplied by 1000 for easy tabulation, following the protocol we previously published [22].

Immunohistochemical staining (IHC)

IHC for RB was performed using the G3-245 mouse monoclonal antibody (BD Biosciences, San Jose, CA) at a dilution of 1:2000 after 20 min of steaming in 10 mM citrate buffer. Labeling was visualized using the DAKO LSAB kit (DAKO, Carpinteria, CA). The percentage of neoplastic cell nuclei labeling in each case was assessed. Slides were considered RB negative when no neoplastic cell nucleus was stained. Detailed procedures and standardized IHC photos for RB staining were shown in a previous paper [23]

IHC for p16 was performed using the monoclonal antibody 16PO4 at a dilution of 1:100 (Cell Marque Inc, Hot Springs, AR). Positive p16 labeling was defined as the presence of nuclear and cytoplasmic reactivity. The percentage of cells showing nuclear and cytoplasmic labeling was recorded. The detailed procedure and standardized IHC photos for p16 staining were shown in the previous paper [24].

For both antibodies, negative controls were performed by omitting the primary antibody incubation step. The staining intensity was evaluated by a senior pathologist and scored as 1 (Strong over 50%), 2 (Stained in 20–50%), 3 (Weak), or 4 (Negative). In addition, we defined p16 or RB expression level as follows: Positive (scores 1 and 2) and Negative (scores 3 and 4) because sample distributions were dramatically changed between HPV positive and negative samples by this threshold (Supplementary Table S4).


Continuous variables were analyzed by the Mann–Whitney U test, and categorical variables were analyzed by Fisher’s exact test, two-tailed. Disease-specific survival was defined as the time from surgery to the date of cancer death. Those who remained alive were censored at the last date the subject was known to be alive. The Association of gene methylation with OS was evaluated using the Cox proportional hazards model with hazard ratios and 95% confidence intervals estimated for multivariable analysis. The individual gene’s methylation level cut-off value was determined by maximizing the sensitivity and specificity using Receiver Operating Characteristics (ROC) curve analysis. All statistical analyses were performed using JMP 16 software (SAS Institute, Cary, U.S.A.). The level of statistical significance was set at P < 0.05.


Characteristics of the study cohort

A total of 127 surgical samples were included in this study. The incidence of HPV-positive and HPV-negative cases is significantly different considering gender, smoking history, and alcohol consumption (Table 1). These findings seem typical of HPV-positive OPSCCs, which are known to be seen in young men without tobacco or alcohol use compared to HPV-negative OPSCCs.

HPV-positive OPSCC-specific methylated genes

We established an optimal cut-off value for each of the eight genes using 50 HPV-positive OPSCC cases and 77 others, including 44 HPV-negative OPSCCs and 33 controls by ROC curve analysis (Fig. 1, Supplementary Figure S1). All HPV-positive OPSCC cases and others were dichotomized into high and low methylated cases based on the cut-off value. The cut-off value and AUC for each of the eight genes are listed in Table 2. All genes, except for CSGALNACT2, GULP1 and RASSF1A, showed significantly higher methylation frequencies in HPV-positive OPSCC cases than others. Mainly, AGTR1 (92% and 48%, P < 0.001) and FHIT (68% and 25%, P < 0.001) were frequently methylated in HPV-positive OPSCCs compared to others with fair AUC values (> 0.700). Interestingly, GULP1 and RASSF1A promoter methylation frequency was significantly higher in HPV-negative cases than in HPV-positive cases.

Fig. 1
figure 1

Scattered plots of QMSP values of tested genes in control samples (n = 33), HPV negative OPSCCs (n = 44) and HPV positive OPSCCs (n = 50)

Table 2 HPV positive OPSCC specific methylation markers

The association of promoter methylation of candidate genes with other clinicopathological parameters is shown in Supplementary Table S3. AGTR1 and RASSF1A methylation were associated with advanced age, while CSGALNACT2, FHIT and GULP1 were associated with race. No drinking history affects DAPK methylation, and VGF methylation influences tumor grades 1–2.

Univariate analysis of disease-specific survival of 94 OPSCCs

The Cox proportional hazards model was used for univariate analysis to find prognostic factors of disease-specific survival of 94 OPSCCs. Age, race, smoking, alcohol consumption, HPV negative, p16 negative and RB positive were significant factors of disease-specific survival (Table 3). Some of the Kaplan–Meier curves are shown in Supplementary Figure S2. These factors were significantly correlated with HPV-negative OPSCCs but not with HPV-positive OPSCCs. None of the methylation markers was a significant prognostic factor, neither for HPV-positive nor HPV-negative OPSCCs. In the multivariable analysis, Age over 60 years old and positive RB were the significant prognostic factor of disease-specific survival.

Table 3 Univariate analysis of disease-specific survival of 94 OPSCCs (Cox proportional hazards model)

Association of p16 and RB expression with HPV infection status

Available immunohistochemical staining (IHC) data of p16 and RB were 89/94 (94.7%) and 54/94 cases (57.4%) respectively. The results of IHC were categorized into four groups. 1: Staining over 50% cells, 2: Staining in 20–50% cells, 3: Staining in 1–20% cells, and 4: Negative. Details of staining scores and the number of samples in each category are available in Supplementary Table S4a. We categorized ‘groups 1–2’ as IHC positive and ‘groups 3–4’ as IHC negative for easy tabulation. Forty-seven out of 89 cases (52.8%) were positive for p16 staining, and 25/54 (46.3%) cases were RB-positive staining (Supplementary Table S4b). Notably, 41/45 (91.1%) HPV-positive cases showed p16-positive staining, while 26/28 (92.9%) HPV-positive cases were RB-negative. As we defined HPV-positive cases in this study by PCR detection of HPV E6 and E7 mRNA detection method, a small number of HPV-positive cases might not indicate HPV-positive by IHC. Conversely, HPV-negative cases showed p16 negative staining (38/44, 86.4%) and RB positive (23/26, 88.5%) for most cases. As expected, HPV-positive cases were positively correlated with p16 expression, while HPV-negative cases were positively correlated with RB expression.

Association of candidate gene promoter methylation with p16 or RB expression

To determine whether there is any correlation of p16 or RB expression with candidate gene methylation, the distribution of QMSP values among IHC positive and IHC groups was compared by the Mann–Whitney U test. Promoter methylation of AGTR1, DAPK, FHIT and TIMP3 genes significantly associated with p16 positive cases (P = 0.045, P < 0.001, P < 0.001, P = 0.030). In contrast, CSGALNACT2, DAPK and FHIT promoter methylations were significantly associated with RB negative cases (P = 0.031, P < 0.001, P = 0.002, Supplementary Table S5). Table 4 summarizes the features of 8 candidate methylation markers. Previously reported 4 cancer-specific methylated genes(DAPK, FHIT and TIMP3) were HPV-positive associated OPSCC markers. Among the novel 4 genes tested in this study, AGTR1 was strongly correlated with HPV-positive OPSCC characteristics.

Table 4 Characteristic summary of 8 methylated genes in OPSCCs

Candidate gene combinations specific to HPV-positive cancers

To improve the specificity of HPV-positive OPSCCs, we picked AGTR1 and FHIT methylation markers based on AUC over 0.7 of ROC analysis results (Supplementary Figure S1). Either gene methylated cases were significantly associated with HPV-positive cancers with 92.0% sensitivity and 46.8% specificity, P < 0.001). Also, they had significantly better overall survival (P = 0.047, Supplementary Figure S2).


Viral infections cause human cancers, including liver, pharyngeal and cervical malignancies [25, 26]. Viruses disrupt the host cell biology and epigenetic process to promote replication and induce the deregulation of various gene promoter methylations [27]. As a result, aberrant DNA methylations are accumulated in non-cancerous or pre-cancerous tissues. It is a phenomenon, so-called the ‘epigenetic field for cancerization’ [28]. Hepatitis virus infection causes chronic hepatitis and induces epigenetic alterations associated with hepatoma risk [29]. Furthermore, the infection is proven to activate a natural killer cell-dependent innate immune response in mice and contribute to aberrant epigenetic accumulation [30].

Our findings solidified that promoter methylation of specific genes correlated with HPV viral infection. Here, we reported some novel genes frequently methylated in HPV-positive human OPSCCs. However, we do not know the precise mechanisms of how HPV infection induces these methylations and the biological effect of these methylations in carcinogenesis processes. The previous report suggests that the E6 and E7 oncoproteins of HPV increase DNA methyltransferase 1 (DNMT1) expression through p53 degradation (24), affecting the promoter methylation of relevant genes. Another group found that methylation arises as a host cell's defense mechanism to silence viral DNA [31]. Although we did not check the expression of each methylated gene in our cohort, TCGA data analysis proves that some of the analyzed gene methylations were associated with expression in OPSCCs.

Several DNA methylation markers have been identified in HPV-positive cervical cancers [32, 33]. Therefore, we examined eight candidate genes in OPSCCs in this study and summarized the result in Table 4. Based on optimal cut-off values, we found OPSCC-specific promoter methylation of AGTR1, DAPK, FHIT, TIMP3 and VGF. Among the analyzed genes, promoter methylation of RASSF1A is less frequent in HPV-positive cases than HPV-negative cases, consistent with the previous report that promoter methylation of RASSF1A is inversely correlated with HPV infection in head and neck and cervical cancer [11, 34]. Although the cut-off values we used need to be validated in another similar cohort, our pilot study supports that some of our panel genes have the potential for future clinical use in different clinical contexts.

In HPV-positive HNSCCs, disruptive p53 mutation is not common [35]. Among the available 28 HPV-positive OPSCCs, most cases showed RB down-regulation (26/28, 92.9%). There is an inverse relationship between p16 expression and RB downregulation (out of 26 RB down-regulated cases, 23 expressed p16, and only 3 cases negative for p16 expressions). This result means that p16 + /RB- phenotype is dominant in HPV-positive cases (23/28, 82.1%) following the previous reports [36,37,38]. Also, this subgroup is known as the high Ki-67 index subgroup in Basal-like carcinomas [39]. In our cohort, p16 + /RB- cases (n = 23) were compared with p16-/RB + cases (n = 24). Only negative smoking history (P = 0.002, Fisher’s exact test) was a significant factor of p16 + /RB- phenotype. Disease-specific survival of p16 + /RB- cases was significantly favorable than p16-/RB + cases (HR = 0.31, P = 0.008, 95%CI: 0.12–0.74). The promoter methylations of AGTR1, CSGALNACT2, DAPK, FHIT and TIMP3 were significantly correlated with the p16 + /RB- phenotype (Table 4). However, there is no correlation between promoter methylation of individual genes and patients' survival. It may be interesting in future studies using a larger cohort to see whether promoter methylations of these genes have an additive effect with p16/RB status.

Among four novel methylation markers in OPSCC (AGTR1, CSGALNACT2, GULP1 and VGF), AGTR1 showed a significant correlation with cancer specificity, HPV infection specificity and p16 + /RB- phenotype (Table 4). AGTR1 is a member of the angiotensin group of G-protein-coupled receptors that play a role in vasoconstriction, salt and water retention, cell proliferation and migration [40]. This gene is reportedly overexpressed in a subset of breast cancers [41, 42]. On the other hand, AGTR1 methylation was significantly found in non-small cell lung cancer (NSCLC) [43] and colorectal cancer [15]. They showed that AGTR1 methylation was found in 60% of NSCLC and 65% of colorectal cancers. Mitra et al. revealed that AGTR1 was a target of oncogenic miR-155, and was silenced via miR-155 up-regulation in head and neck cancer [44]. They also speculated that AGTR1 expression was influenced by promoter hypermethylation.

A panel of OPSCC-specific and HPV-specific methylation markers can be available for early oropharyngeal cancer detection. One of the advantages of methylation marker detection is its high sensitivity and rapidity. We previously reported a rapid gene methylation detection procedure for head and neck cancers [45]. It takes less than three hours from sampling to get the result. We also reported that a few markers combination produces a definitive decision of cancer cell existence [46]. Thus, unlike conventional HPV PCR or p16 IHC examination, methylation detection analysis with excellent markers for sampled tissues or salivas has the potential to make a rapid and definitive diagnosis of malignancies originated from the HPV-infected field of cancerization. Furthermore, it may enable recurrent non-invasive examinations for high-risk patients or monitor postoperative advanced cancer patients.

There are some limitations of our study. First, since most cases were symptomatic and advanced and the samples came from surgically resected cases, the results may partially reflect the characteristics of advanced cancers. Second, an optimal cut-off of methylation level may be the cohort-specific value, which needs to be optimized by multiple data set testing. Third, we used a relatively small number of the single institute sample set, and data should be ascertained by more extensive multicenter cohort data. Finally, biological studies must be performed in contrast with HPV infection and promoter methylation.

In summary, we have detected a panel of novel HPV-positive and OPSCC-specific methylation markers. Especially, the combination of ATGR1 methylation and FHIT methylation could be a sensitive biomarker of HPV-positive OPSCCs, which has a characteristic of the p16 positive and RB negative phenotypes.

Availability of data and materials

The datasets generated and analyzed during the current study are not publicly available due to the ethical restrictions on surgery data and histological information. Still, they are available from the corresponding author upon reasonable request.



HPV-positive oropharyngeal squamous cell carcinomas


Quantitative-methylation specific PCR


Area under the curve


Human papillomavirus




Head and neck squamous cell carcinomas


Immunohistochemical staining


  1. D’Souza G, Kreimer AR, Viscidi R, Pawlita M, Fakhry C, Koch WM, Westra WH, Gillison ML. Case-control study of human papillomavirus and oropharyngeal cancer. N Engl J Med. 2007;356(19):1944–56.

    Article  CAS  PubMed  Google Scholar 

  2. Mahal BA, Catalano PJ, Haddad RI, Hanna GJ, Kass JI, Schoenfeld JD, Tishler RB, Margalit DN. Incidence and Demographic Burden of HPV-Associated Oropharyngeal Head and Neck Cancers in the United States. Cancer Epidemiol Biomarkers Prev. 2019;28(10):1660–7.

    Article  PubMed  Google Scholar 

  3. Marur S, D’Souza G, Westra WH, Forastiere AA. HPV-associated head and neck cancer: a virus-related cancer epidemic. Lancet Oncol. 2010;11(8):781–9.

    Article  PubMed  PubMed Central  Google Scholar 

  4. Scheffner M, Werness BA, Huibregtse JM, Levine AJ, Howley PM. The E6 oncoprotein encoded by human papillomavirus types 16 and 18 promotes the degradation of p53. Cell. 1990;63(6):1129–36.

    Article  CAS  PubMed  Google Scholar 

  5. Wiest T, Schwarz E, Enders C, Flechtenmacher C, Bosch FX. Involvement of intact HPV16 E6/E7 gene expression in head and neck cancers with unaltered p53 status and perturbed pRb cell cycle control. Oncogene. 2002;21(10):1510–7.

    Article  CAS  PubMed  Google Scholar 

  6. van Kempen PM, Noorlag R, Braunius WW, Stegeman I, Willems SM, Grolman W. Differences in methylation profiles between HPV-positive and HPV-negative oropharynx squamous cell carcinoma: a systematic review. Epigenetics. 2014;9(2):194–203.

    Article  PubMed  Google Scholar 

  7. Ha PK, Pai SI, Westra WH, Gillison ML, Tong BC, Sidransky D, Califano JA. Real-time quantitative PCR demonstrates low prevalence of human papillomavirus type 16 in premalignant and malignant lesions of the oral cavity. Clin Cancer Res. 2002;8(5):1203–9.

    CAS  PubMed  Google Scholar 

  8. Brait M, Maldonado L, Begum S, Loyo M, Wehle D, Tavora FF, Looijenga LHJ, Kowalski J, Zhang Z, Rosenbaum E, et al. DNA methylation profiles delineate epigenetic heterogeneity in seminoma and non-seminoma. Br J Cancer. 2012;106(2):414–23.

    Article  CAS  PubMed  Google Scholar 

  9. Sanchez-Cespedes M, Esteller M, Wu L, Nawroz-Danish H, Yoo GH, Koch WM, Jen J, Herman JG, Sidransky D. Gene promoter hypermethylation in tumors and serum of head and neck cancer patients. Cancer Res. 2000;60(4):892–5.

    CAS  PubMed  Google Scholar 

  10. Loyo M, Brait M, Kim MS, Ostrow KL, Jie CC, Chuang AY, Califano JA, Liegeois NJ, Begum S, Westra WH, et al. A survey of methylated candidate tumor suppressor genes in nasopharyngeal carcinoma. Int J Cancer. 2011;128(6):1393–403.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Dong SM, Sun DI, Benoit NE, Kuzmin I, Lerman MI, Sidransky D. Epigenetic inactivation of RASSF1A in head and neck cancer. Clin Cancer Res. 2003;9(10 Pt 1):3635–40.

    CAS  PubMed  Google Scholar 

  12. Lo KW, Kwong J, Hui AB, Chan SY, To KF, Chan AS, Chow LS, Teo PM, Johnson PJ, Huang DP. High frequency of promoter hypermethylation of RASSF1A in nasopharyngeal carcinoma. Cancer Res. 2001;61(10):3877–81.

    CAS  PubMed  Google Scholar 

  13. Righini CA, de Fraipont F, Timsit JF, Faure C, Brambilla E, Reyt E, Favrot MC. Tumor-specific methylation in saliva: a promising biomarker for early detection of head and neck cancer recurrence. Clin Cancer Res. 2007;13(4):1179–85.

    Article  CAS  PubMed  Google Scholar 

  14. Sun W, Zaboli D, Wang H, Liu Y, Arnaoutakis D, Khan T, Khan Z, Koch WM, Califano JA. Detection of TIMP3 promoter hypermethylation in salivary rinse as an independent predictor of local recurrence-free survival in head and neck cancer. Clin Cancer Res. 2012;18(4):1082–91.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Carmona FJ, Azuara D, Berenguer-Llergo A, Fernandez AF, Biondo S, de Oca J, Rodriguez-Moranta F, Salazar R, Villanueva A, Fraga MF, et al. DNA methylation biomarkers for noninvasive diagnosis of colorectal cancer. Cancer Prev Res (Phila). 2013;6(7):656–65.

    Article  CAS  Google Scholar 

  16. Maldonado L, Brait M, Izumchenko E, Begum S, Chatterjee A, Sen T, Loyo M, Barbosa A, Poeta ML, Makarev E, et al. Integrated transcriptomic and epigenomic analysis of ovarian cancer reveals epigenetically silenced GULP1. Cancer Lett. 2018;433:242–51.

    Article  CAS  PubMed  Google Scholar 

  17. Hayashi M, Guida E, Inokawa Y, Goldberg R, Reis LO, Ooki A, et al. GULP1 regulates the NRF2-KEAP1 signaling axis in urothelial carcinoma. Sci Signal. 2020;13(645):eaba0443.

  18. Brait M, Maldonado L, Noordhuis MG, Begum S, Loyo M, Poeta ML, Barbosa A, Fazio VM, Angioli R, Rabitti C et al: Association of promoter methylation of VGF and PGP9.5 with ovarian cancer progression. PLoS One 2013, 8(9):e70878.

  19. Hayashi M, Bernert H, Kagohara LT, Maldonado L, Brait M, Schoenberg M, Bivalacqua T, Netto GJ, Koch W, Sidransky D, et al. Epigenetic inactivation of VGF associated with Urothelial Cell Carcinoma and its potential as a non-invasive biomarker using urine. Oncotarget. 2014;5(10):3350–61.

    Article  PubMed  PubMed Central  Google Scholar 

  20. Hoque MO, Kim MS, Ostrow KL, Liu J, Wisman GB, Park HL, Poeta ML, Jeronimo C, Henrique R, Lendvai A, et al. Genome-wide promoter analysis uncovers portions of the cancer methylome. Cancer Res. 2008;68(8):2661–70.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Brait M, Maldonado L, Noordhuis M, Begum S, Loyo M, Poeta ML, Barbosa A, Fazio VM, Angioli R, Rabitti C et al: Association of Promoter Methylation of VGF and PGP9.5 with Ovarian Cancer Progression. PloS one 2013, 8(9):e70878.

  22. Hoque MO, Begum S, Topaloglu O, Chatterjee A, Rosenbaum E, Van Criekinge W, Westra WH, Schoenberg M, Zahurak M, Goodman SN, et al. Quantitation of promoter methylation of multiple genes in urine DNA and bladder cancer detection. J Natl Cancer Inst. 2006;98(14):996–1004.

    Article  CAS  PubMed  Google Scholar 

  23. Yoo GH, Xu HJ, Brennan JA, Westra W, Hruban RH, Koch W, Benedict WF, Sidransky D. Infrequent inactivation of the retinoblastoma gene despite frequent loss of chromosome 13q in head and neck squamous cell carcinoma. Can Res. 1994;54(17):4603–6.

    CAS  Google Scholar 

  24. Chaux A, Pfannl R, Rodríguez IM, Barreto JE, Velazquez EF, Lezcano C, Piris A, Netto GJ, Cubilla AL. Distinctive immunohistochemical profile of penile intraepithelial lesions: a study of 74 cases. Am J Surg Pathol. 2011;35(4):553–62.

    Article  PubMed  Google Scholar 

  25. Hattori N, Ushijima T. Epigenetic impact of infection on carcinogenesis: mechanisms and applications. Genome medicine. 2016;8(1):10.

    Article  PubMed  PubMed Central  Google Scholar 

  26. Sadhukhan P, Ugurlu MT, Hoque MO: Effect of COVID-19 on Lungs: Focusing on Prospective Malignant Phenotypes. Cancers 2020, 12(12).

  27. Milavetz BI, Balakrishnan L. Viral epigenetics. Methods Mol Biol. 2015;1238:569–96.

    Article  PubMed  PubMed Central  Google Scholar 

  28. Ushijima T. Epigenetic field for cancerization. J Biochem Mol Biol. 2007;40(2):142–50.

    CAS  PubMed  Google Scholar 

  29. Hamdane N, Jühling F, Crouchet E, El Saghire H, Thumann C, Oudot MA, Bandiera S, Saviano A, Ponsolles C, Roca Suarez AA, et al. HCV-Induced Epigenetic Changes Associated With Liver Cancer Risk Persist After Sustained Virologic Response. Gastroenterology. 2019;156(8):2313-2329.e2317.

    Article  PubMed  Google Scholar 

  30. Okamoto Y, Shinjo K, Shimizu Y, Sano T, Yamao K, Gao W, Fujii M, Osada H, Sekido Y, Murakami S, et al. Hepatitis virus infection affects DNA methylation in mice with humanized livers. Gastroenterology. 2014;146(2):562–72.

    Article  CAS  PubMed  Google Scholar 

  31. Jimenez-Wences H, Peralta-Zaragoza O, Fernandez-Tilapa G. Human papilloma virus, DNA methylation and microRNA expression in cervical cancer (Review). Oncol Rep. 2014;31(6):2467–76.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Farkas SA, Milutin-Gasperov N, Grce M, Nilsson TK. Genome-wide DNA methylation assay reveals novel candidate biomarker genes in cervical cancer. Epigenetics. 2013;8(11):1213–25.

    Article  CAS  PubMed  Google Scholar 

  33. Leonard SM, Wei W, Collins SI, Pereira M, Diyaf A, Constandinou-Williams C, Young LS, Roberts S, Woodman CB. Oncogenic human papillomavirus imposes an instructive pattern of DNA methylation changes which parallel the natural history of cervical HPV infection in young women. Carcinogenesis. 2012;33(7):1286–93.

    Article  CAS  PubMed  Google Scholar 

  34. Lai HC, Lin YW, Chang CC, Wang HC, Chu TW, Yu MH, Chu TY: Hypermethylation of two consecutive tumor suppressor genes, BLU and RASSF1A, located at 3p21.3 in cervical neoplasias. Gynecologic oncology 2007, 104(3):629–635.

  35. Westra WH, Taube JM, Poeta ML, Begum S, Sidransky D, Koch WM. Inverse relationship between human papillomavirus-16 infection and disruptive p53 gene mutations in squamous cell carcinoma of the head and neck. Clin Cancer Res. 2008;14(2):366–9.

    Article  CAS  PubMed  Google Scholar 

  36. Rezaei M, Mostafaei S, Aghaei A, Hosseini N, Darabi H, Nouri M, Etemadi A, Neill AO, Nahand JS, Mirzaei H, et al. The association between HPV gene expression, inflammatory agents and cellular genes involved in EMT in lung cancer tissue. BMC Cancer. 2020;20(1):916.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. Zhang J, Chen T, Yang X, Cheng H, Späth SS, Clavijo PE, Chen J, Silvin C, Issaeva N, Su X, et al. Attenuated TRAF3 Fosters Activation of Alternative NF-κB and Reduced Expression of Antiviral Interferon, TP53, and RB to Promote HPV-Positive Head and Neck Cancers. Cancer Res. 2018;78(16):4613–26.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. Wang H, Zhang Y, Bai W, Wang B, Wei J, Ji R, Xin Y, Dong L, Jiang X. Feasibility of Immunohistochemical p16 Staining in the Diagnosis of Human Papillomavirus Infection in Patients With Squamous Cell Carcinoma of the Head and Neck: A Systematic Review and Meta-Analysis. Front Oncol. 2020;10: 524928.

    Article  PubMed  PubMed Central  Google Scholar 

  39. Subhawong AP, Subhawong T, Nassar H, Kouprina N, Begum S, Vang R, Westra WH, Argani P. Most basal-like breast carcinomas demonstrate the same Rb-/p16+ immunophenotype as the HPV-related poorly differentiated squamous cell carcinomas which they resemble morphologically. Am J Surg Pathol. 2009;33(2):163–75.

    Article  PubMed  PubMed Central  Google Scholar 

  40. de Ronde JJ, Lips EH, Mulder L, Vincent AD, Wesseling J, Nieuwland M, Kerkhoven R, Vrancken Peeters MJ, Sonke GS, Rodenhuis S, et al. SERPINA6, BEX1, AGTR1, SLC26A3, and LAPTM4B are markers of resistance to neoadjuvant chemotherapy in HER2-negative breast cancer. Breast Cancer Res Treat. 2013;137(1):213–23.

    Article  CAS  PubMed  Google Scholar 

  41. Ateeq B, Tomlins SA, Chinnaiyan AM. AGTR1 as a therapeutic target in ER-positive and ERBB2-negative breast cancer cases. Cell Cycle. 2009;8(23):3794–5.

    Article  CAS  PubMed  Google Scholar 

  42. Rhodes DR, Ateeq B, Cao Q, Tomlins SA, Mehra R, Laxman B, Kalyana-Sundaram S, Lonigro RJ, Helgeson BE, Bhojani MS, et al. AGTR1 overexpression defines a subset of breast cancer and confers sensitivity to losartan, an AGTR1 antagonist. Proc Natl Acad Sci U S A. 2009;106(25):10284–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  43. Guo S, Yan F, Xu J, Bao Y, Zhu J, Wang X, Wu J, Li Y, Pu W, Liu Y, et al. Identification and validation of the methylation biomarkers of non-small cell lung cancer (NSCLC). Clin Epigenetics. 2015;7(1):3.

    Article  PubMed  PubMed Central  Google Scholar 

  44. Mitra S, Mukherjee N, Das S, Das P, Panda CK, Chakrabarti J. Anomalous altered expressions of downstream gene-targets in TP53-miRNA pathways in head and neck cancer. Sci Rep. 2014;4:6280.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  45. Hayashi M, Guerrero-Preston R, Okamura J, Michailidi C, Kahn Z, Li X, Ahn J, Goldsmith M, Koch W. Innovative rapid gene methylation analysis of surgical margin tissues in head and neck cancer. Ann Surg Oncol. 2014;21(9):3124–31.

    Article  PubMed  PubMed Central  Google Scholar 

  46. Hayashi M, Wu G, Roh JL, Chang X, Li X, Ahn J, Goldsmith M, Khan Z, Bishop J, Zhang Z, et al. Correlation of gene methylation in surgical margin imprints with locoregional recurrence in head and neck squamous cell carcinoma. Cancer. 2015;121(12):1957–65.

    Article  CAS  PubMed  Google Scholar 

Download references




This work was supported by a career development award to M. O. Hoque from the Specialized Program of Research Excellence and P50 CA098252 (T-C Wu) National Institutes of Health (NIH)/National Cancer Institute (NCI). The funders had no role in study design, data collection and analysis, publication decisions, or manuscript preparation.

Author information

Authors and Affiliations



Y.I., M.H.and M.O.H. designed the project. Y.I., M.H., S.B. and M.B.performed the experiments. Y.I., M.H., M.N.and M.O.H. analyzed the data. Y.I., M.H. S.B. and M.B. checked and approved all the statistical analyses. D.S., J.C., W.K., W.W. and M.O.H. provide the materials, reagents, and analytic tools. Y.I., M.H., M.B. and M.O.H. wrote the manuscript. The author(s) read and approved the final manuscript.

Corresponding author

Correspondence to Mohammad O. Hoque.

Ethics declarations

Ethics approval and consent to participate

Ethics approval was obtained from the Johns Hopkins Medicine Internal Review Board and performed under research protocol NA_00036235. The study was in accordance with the exemption under the U.S. Department of Health and Human Services policy for protection of human subjects [45 CFR 46.101(b)].

Consent for publication

Not applicable.

Competing interests

The authors declare no potential conflicts of interest.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Supplementary Figure S1.

ROC Curves for determining optimal cut-off values between HPV-positive OPSCCs (n=50) and other samples including HPV-negative OPSCCs (n=44) and control samples (n=33). Supplementary Figure S2. Disease-specific survival curves of 94 OPSCCs depending on various clinicohistological factors, including HPV status.

Additional file 2: Supplementary Table S1.

Candidate promoter methylated genes. Supplementary Table S2. Primers and probes for QMSP. Supplementary Table S3. Association of methylation markers with clinicopathological data (n=94). Supplementary Table S4. Classification of cases based on p16 and RB Staining. Supplementary Table S5. Correlation of QMSP values with p16 and RB expression (Mann-Whitney U test).

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Inokawa, Y., Hayashi, M., Begum, S. et al. High-risk HPV infection-associated hypermethylated genes in oropharyngeal squamous cell carcinomas. BMC Cancer 22, 1146 (2022).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: