Skip to main content

HPV infection and 5mC/5hmC epigenetic markers in penile squamous cell carcinoma: new insights into prognostics

Abstract

Background

Penile cancer is one of the most aggressive male tumors. Although it is preventable, the main etiologic causes are lifestyle behaviors and viral infection, such as human papillomavirus (HPV). Long-term epigenetic changes due to environmental factors change cell fate and promote carcinogenesis, being an important marker of prognosis. We evaluated epidemiological aspects of penile squamous cell carcinoma (SCC) and the prevalence of HPV infection using high-risk HPV (hrHPV) and p16INK4A expression of 224 participants. Global DNA methylation was evaluated through 5-methylcytosine (5mC) and 5-hydroxymethylcytosine (5hmC).

Results

The incidence of HPV was 53.2% for hrHPV and 22.32% for p16INK4a. hrHPV was not related to systemic or lymph node metastasis and locoregional recurrence, nor influenced the survival rate. P16INK4a seems to be a protective factor for death, which does not affect metastasis or tumor recurrence. Lymph node and systemic metastases and locoregional recurrence increase the risk of death. An increased 5mC mark was observed in penile SCC regardless of HPV infection. However, there is a reduction of the 5hmC mark for p16INK4a + (P = 0.024). Increased 5mC/5hmC ratio (> 1) was observed in 94.2% of penile SCC, irrespective of HPV infection. Despite the increase in 5mC, it seems not to affect the survival rate (HR = 1.06; 95% CI 0.33–3.38).

Conclusions

P16INK4a seems to be a good prognosis marker for penile SCC and the increase in 5mC, an epigenetic mark of genomic stability, may support tumor progression leading to poor prognosis.

Background

Penile cancer is a rare tumor type with an increased worldwide incidence, having 36,068 new cases registered in 2020. The majority of cases occur in regions with a low human development index, where India (10,677), China (4,628) and Brazil (1,698) are the most affected countries [1]. Despite the low incidence compared to other types of male malignancies, penile cancer has a poor prognosis and is associated with high mortality [2] and morbidity [3]. Therapeutic strategies are very limited and, therefore, the main related therapy is partial or total penectomy, which directly affects the emotional and social life of patients. Given narrow options for early diagnosis and non-surgery treatments, restricted funding for medical care, and mutilating treatments resulting in negative effects on well-being, penile cancer can be regarded as a neglected disease [4,5,6].

Different histological types are associated with penile cancer, such as sarcoma, melanoma and basal cell carcinoma [7, 8], nevertheless, penile squamous cell carcinoma (SCC) is reported in 95% of cases worldwide [9, 10]. Multiple risk factors are described, mostly related to lifestyle behaviors, such as promiscuous sexual behavior [11], history of zoophilia [12], poor hygiene [13], psoralen UV-A phototherapy [14], smoking [15] and obesity [16]. Non-circumcision (phimosis) leads to chronic inflammation conditions like posthitis, lichen sclerosus and balanitis xerotic obliterans [9] increasing the risk of developing penile cancer by 22-fold [17]. However, human papillomavirus (HPV) infection in penile cancer is one of the main etiologic causes, especially for squamous cell carcinomas [18].

The prevalence of HPV in penile neoplasia can vary widely depending on the literature and among different regions of the world, ranging from 11 to 87% [19, 20]. Associations of HPV infection and death risk are still unclear, as the results are controversial [21,22,23]. Positive survival prognosis for high-risk HPV (hrHPV) in penile cancer has already been highlighted [24] and Wang and collaborators (2020) [25] have recently suggested a regional lymph node infiltration staging based on the presence of hrHPV. Subtypes of hrHPV are linked to malignant lesions due to the degradation of the cell cycle control proteins P53 and Rb, and the expression of the viral oncoproteins, HPV E6 and HPV E7, which causes evasion to cell death and DNA damage repair, respectively [26]. These viral oncoproteins' expression promotes chromosomal rearrangements, multiple centromeres and aneuploidy [27]. Viral infection, as HPV is often related to an increase in genomic instability by aberrantly reprogramming the epigenome [28]. In different squamous cell carcinomas hrHPV directly modulates enzymes that maintain the conformation of nucleosomes [28, 29] and enzymes responsible for maintaining DNA methylation [30, 31].

Environmentally induced epigenetic changes have been recently added as hallmarks of cancer, which contribute to tumor initiation and progression [32]. Global DNA methylation is characterized by the addition of a methyl at position 5 of cytosine in a dinucleotide CpG (cytosine-phosphate-guanine) resulting in a 5-methylcytosine (5mC) [33], which affects gene expression, chromatin remodeling and genomic stability. DNA methyltransferases (DNMTs) mediate the transfer of the methyl group to DNA, and loss of the 5mC marker can occur either passively by DNA replication or actively by the enzymes known as Ten-Eleven-Translocation (TETs) [34]. These enzymes catalyze the conversion of 5mC into 5-hydroxymethylcytosine (5hmC). Imbalances between 5mC and 5hmC marks cause transcriptional dysregulation of promoters and enhancers, changing the cell fate [35]. The imbalance of 5mc/5hmc dynamics has been described in many types of cancers, but not for penile squamous cell carcinomas.

Since HPV infection may drive epigenetic changes, the characterization of the global DNA methylation and its association with hrHPV DNA and p16INK4a expression will provide better knowledge about the molecular mechanisms related to penile SCC. Therefore, evaluate the prevalence of HPV infection and 5mC and 5hmC epigenetic marks in penile SCC and its association with clinicopathological alterations. As an important mechanism of genomic stability, aberrant epigenetic reprogramming related to penile SCC pathogenesis and viral infection may be used as a biomarker for prognosis and targeted therapies. Determining the prevalence of HPV, one of the most common infections of the reproductive tract, is important for the development of public health strategies for the prevention and treatment of related diseases.

Results

Clinicopathological characteristics and HPV infection

The participants’ age, HPV infection, pathological classification, and treatments of penile SCC are presented in Table 1. Participants' mean age was 63.8 (± 15.86) years old, and the incidence of SCC was higher in men over 60 years old (59.2%; 133/224), even though the manifestation of the disease was also detected before the age of 60 (33.2%; 74/224) and 40 (7.6%; 17/224). The frequency of HPV infection was 22.3% when evaluated with p16INK4a immunodepression and reached 53.2% when tested by hrHPV hybridization. Low-risk HPV infection was not identified. Tumor staging based on the primary tumor was mainly corpus spongiosum invasion (39.3%, pT2) followed by the first stages including carcinoma in situ (pTis), noninvasive carcinoma (pTa), subepithelial invasion without lymphovascular invasion (pT1a) and with lymphovascular invasion (pT1b) (30.1%).

Table 1 Clinical and pathological aspects, staging and treatment of penile squamous cell carcinoma participants

Systemic metastasis was observed in 16.5% of participants diagnosed with penile SCC, while lymph node metastasis was presented in 40.1% and locoregional recurrence in 26.1%. Intense peritumoral lymphocyte infiltrate (PLI) and intratumoral polymorphonuclear infiltrate (IPI) reached 40.1% and 41.5% of penile SCC, respectively. The main therapeutic strategies were partial amputation 73.2% and total amputation 14.7%, followed by other associated strategies (12.1%) such as prostatectomy, emasculation, and exercise injury. Chemotherapy was the treatment option for 25.8% of men, which consisted of cisplatin associated with 5-Fluorouracil or cisplatin associated with taxol or paclitaxel, as a palliative or neoadjuvant therapy before and after lymphadenectomy. Radiotherapy was used only as a palliative alternative for pain in cases of metastasis (14.1%) (Table 1). Despite the differences in HPV infection diagnosis using p16INK4a or hrHPV, a positive correlation was observed between these markers (P = 0.003, Table 2). Radiotherapy was associated with p16INK4a positive cases (P = 0.0136). All other clinicopathological variables were not associated with HPV markers p16INK4a and hrHPV.

Table 2 Clinical-pathological aspects, staging and therapy of penile SCC patients and p16INK4a and hrHPV correlation

The survival rate was measured by a multivariate-adjusted Cox hazards regression model (Fig. 1A, Additional file 1: Supplementary Table 1). Locoregional recurrence (Hazard Ratio, HR = 5.51, 95% CI 2.20–13.84; P < 0.001), systemic metastasis (HR = 4.40, 95% CI 1.68–11.49, P = 0.003) and lymph node metastasis (HR = 4.65, 95% CI 1.3–15.7, P = 0.013) decreased the survival rate in penile SCC. On the other hand, p16INK4a expression seems to be a protective factor for death (HR = 0.34, 95% CI 0.531–2.05; P = 0.04) and hrHPV (HR = 0.82, 95% CI 0.36–1.86; P = 0.637) did not affect survival even with a HR < 1. The remaining parameters evaluated did not affect the survival rate. However, on the non-adjusted Kaplan–Meier curves metastasis, recurrence, chemotherapy, radiotherapy (Fig. 1B–D, F, M) and staging (Additional file 1: Supplementary Fig. 3E) seems to decrease the time to death (P < 0.05).

Fig. 1
figure 1

Survival analyses. Multivariate-adjusted Cox hazards regression model (A) and Kaplan–Meier survival curves: B lymph node metastasis; C systemic metastasis; D locoregional recurrence; E surgery; F radiotherapy; G p16 INK4a; H hrHPV; I double positive (+), p16INK4a (+), hrHPV (+) and negative (−); J 5-mehtylcitosine, 5mC; L 5-hydroximetilcitosine, 5hmC and M chemotherapy. Surgery 1, partial amputation; Surgery 2, total amputation; Surgery 3, others. *Adjusted with surgery 1

5mC and 5hmC marks in penile SCC

Increased global 5mC (> 60%) mark was observed in 89.2% of penile SCC cases, regardless of HPV infection (Fig. 2A , B). The average of global 5hmC levels varied from lower (< 30%), intermediate (30–60%), and increased (> 60%) hydroxymethylation levels (Fig. 2B). The total cases of penile cancer evaluated for the 5mC mark (Fig. 2C), corresponding to 77.5% ± 19.9, have a higher percentage value than the 5hmC mark (35.9% ± 29.3, Fig. 2D). No association was observed between global DNA methylation marks and HPV infection (Table 2). Also, no differences were observed in the distribution of the 5mC mark for both hrHPV + and p16INK4a + (Fig. 2C), but an increase in the 5hmC mark for p16INK4a negative (P = 0.024) was observed (Fig. 2D). 5mC in double positive HPV was similar to hrHPV and p16INK4a (Fig. 2C) and 5hmC double positive HPV was similar to p16INK4a + (Fig. 2D). The 5mC/5hmC ratio represents the proportion between the marks for each sample (Fig. 3A). The majority of penile SCC showed 5mC/5hmC ratio greater or equal to 1 (94.2%), which was not associated with viral infection (Fig. 3B). Despite increased 5mC mark, methylation level was not associated with survival rate (HR = 1.05, 95% CI 0.33–3.38; P = 0.92), as well as 5hmC (HR = 0.95, 95% CI 0.60–1.50; P = 0.82).

Fig. 2
figure 2

The pattern of 5mC and 5hmC epigenetic marks on penile SCC. A 5mC and 5hmC marks staining in penile squamous cell carcinoma. B Stratification of 5mC and 5hmC mark levels. Distribution of 5mC C and 5hmC D levels in a total number of samples and separated by hrHPV and p16INK4a groups. Data are presented as mean and standard deviation. *P = 0.024

Fig. 3
figure 3

5mC/5hmC ratio distribution. A 5mC/5hmC ratio dispersion by 5mC and 5hmC percentage, y and x-axis, respectively. The higher the color tone, the higher the ratio value. B 5mC/5hmC ratio distribution values by hrHPV and p16INK4a groups

Discussion

Penile cancer is one of the most aggressive male malignancies, although it is one of the easiest to prevent, as the main etiologic causes are lifestyle conditions and viral infection [38]. We observed a high incidence of HPV infection considering hrHPV detection (53.2%) while the incidence reduces using p16INK4a (22.3%). Despite the elevated incidence, HPV infection was not related to systemic or lymph node metastasis and locoregional recurrence and, therefore, does not influence the survival rate. Long-term changes in DNA methylation are characteristic of environmentally induced carcinogenesis as penile cancer [13]. Increased 5mC was observed in penile SCC, which seems to be a stable epigenetic marker, while the 5hmC was lower and more variable. The HPV infection was not related to the overall 5mC/5hmC ratio. Interestingly, 5mC may be a risk factor for poor prognosis and survival rate.

The incidence of penile cancer directly impacts public health, as the specialized medical care, psychosocial impact and the risk of cervical cancer in their female partners increase the cost of treatment [66]. Despite significant reductions in hospital admissions in Brazil over the past two decades, the incidence of penile carcinoma is still high and quite unequal across the five regions of the country, whereby the Northwest stands out with the highest incidence and the lowest per capita income in the country [39]. Socioeconomic status and low education levels reflect on penile SCC staging at diagnosis [40, 41] and, a recent study conducted in the states of Maranhão and Rio de Janeiro, Brazil, report that 69.9% of patients were diagnosed with pT2 and pT3, and 87.9% of them underwent surgical removal of the organ. A contrasting scenario is observed in developed countries, such as the USA, where 54.1% of penile SCC staging are pT1 and there were 75.8% of amputated cases [42]. Penile SCC was observed in young men, ranging from 18 to 103 years old, different from other cities in Brazil as the age varied between 23 to 98 years [13, 41], and even more discrepant in developed countries, where the disease is rare before the third decade of life [43, 44].

Despite the correlation between hrHPV and p16INK4a, the incidence of HPV infection was twofold higher using hrHPV than p16INK4a, and only 17% of all cases were double positive for both markers. Considering that the tumor suppressor protein p16INK4a is upregulated by the HPV oncogene [45, 46], HPV DNA detection may be a primary event that is followed by increased expression of this protein. Also, negative testing for p16INK4a expression may be linked to gene expression loss, caused by promoter methylation events and loss of heterozygosity [45, 46], implying a false negative result [49]. A recent systematic review and meta-analysis showed 49% (95% CI 43.1–54.9) of HPV prevalence using qPCR or hybridization analyses (n = 3772 cases in 47 studies) and 42.1% (95% CI 36.4–47.8) using p16INK4a (n = 1296 cases in 23 studies) in SCC [47]. Both techniques are widely used for HPV detection, however, the values may show discrepancies. Therefore, double positive HPV DNA and p16INK4a must be considered for some types of cancer, such as head and neck SCC [74].

Although previously mentioned as a risk factor for penile SCC [10, 23], our present analysis indicates that hrHPV and p16INK4a were not related to disease staging, metastasis, or infiltration of the immune system cells. On the other hand, p16INK4a seems to be a protective factor for death, as well as hrHPV (HR < 1), although it was not significant, an increase in sample size would answer this question. In addition, positive cases for p16INK4a seem to be associated with radiotherapy. Although previous studies associate p16INK4a expression and radiotherapy [75, 76], this might be a random correlation. Increased p16INK4a was correlated with a good prognosis, as lack of nodal metastasis in head and neck SCC [67, 77] and cervical cancer [78]. In penile SCC, Martins and collaborators showed that p16INK4a was not associated with prognosis parameters or survival rate [46].

The presence of lymph node metastasis in the inguinal and iliac regions increased the risk of death, which is the most common metastatic event [50] and the main prognostic marker [51]. Systemic metastases are rarely reported [41, 50] and locoregional recurrences are also less frequent [52], but these two parameters were associated with an increased risk of death by 4.3 and 5.5-fold, respectively. Penectomy is still a widely used treatment for extensive lesions or tumors involving the base and the bulbar urethral part of the penis [13, 41, 53], but our results indicate that this type of surgery was not linked to survival rate. However, morbidity and psychological traumas are detrimental to the quality of life [53]. The absence of minimally invasive therapies and predictive characteristics of penile cancer is an overwhelming burden for medical practice and scientists, who are challenged to search for molecular targets of the disease [54].

Epigenetic mechanisms, like DNA methylation, present great plasticity as the epigenome can be reprogrammed through environmental factors [68, 69]. Aberrant epigenetic reprogramming changes cell landscape given the malignant phenotype, as increased proliferation, resistance to apoptosis, invasion, and metastasis, triggering tumorigenesis [32]. Several studies have been dedicated to the characterization of epimutations in penile SCC [19, 70]. In this study, an increased level of the global 5mC mark was observed regardless of HPV infection. Genomic instability and mutations are important molecular alterations that trigger carcinogenesis. Increased global 5mC (hypermethylation) contributes to genomic stability in tumor progression and, therefore, may be related to tumor invasion and chemotherapy resistance [58, 73]. Unlike our findings, reduced global DNA methylation was reported in other types of the tumors such as colorectal [55] and prostate [56]. As a repressive epigenetic marker, it changes gene expression of specific targets and large chromosomal regions, as repetitive DNA elements [58].

The increased 5mC marker was accompanied by a reduced level of the 5hmc mark in the majority of penile SCC, which was associated with a p16INK4a negative diagnosis. Decreased level of 5hmC mark was observed in head and neck cancers positive for HPV, and virus presence might influence the oxidation process from 5mC to 5hmC in genes of cell junction pathways [65]. 5hmC is an epigenetic marker for active DNA demethylation [57], which depends on TET proteins that mediate this process [71]. Aberrant active DNA demethylation was also observed in non-small cell lung cancer [72]. Furthermore, we observed a complete absence of 5hmC staining in several SCC samples, similar to previous reports about global loss of the 5hmC mark in large subsets of oral squamous cell carcinoma [59] and cervical cancer [60].

Given that 5mC is a protective marker against DNA damage and the hydroxymethylation of 5mC is an active demethylation process, the balance between 5mC and 5hmC marks is directly linked to genomic stability, contributing to cancer development and progression [34]. We observed that a higher 5mC may be a risk factor for death in penile SCC, irrespective of HPV infection, and an increase in sample size may confirm these findings. The 5mC/5hmC mark is an important prognostic and predictive parameter, and we have recently demonstrated the 5mC/5hmC imbalance related to hypercellular bone marrow, dyserythropoiesis and cases of high-risk myelodysplastic syndromes [61]. Moreover, the clinical relevance of 5mC and 5hmC levels was reported as a critical marker for prognosis in colorectal cancer, as increased 5mC was associated with lymph node metastasis [55]. Global epigenetic marks should be deciphered as they indicate characteristics of advanced tumor staging [60] and tumor subtypes biomarker [62], and targeted therapy using DNMT inhibitors [63, 64].

Despite the important findings regarding the HPV incidence and the association with the 5mC/5hmC mark with penile SCC prognosis, other studies should address the epidemiological characteristics of penile cancer and epigenetic reprogramming through lifestyle changes. Our study was limited by the lack of information regarding the participant's background, such as sexual behavior, sociodemographic profile, HPV vaccination, and the identification of other related histological subtypes and tumor topography. Even though we evaluated hrHPV and lrHPV, there is a diversity of viral genotypes that have not been identified individually and should give better knowledge regarding the infection and disease prognosis, as well as the time of exposure to infection that was not tracked. Although partial and total penectomy represents better survival compared to even more invasive ones, these techniques have drastic implications for the individual's personal life, which calls for less invasive therapies. The epidemiologic data, for staging and surgical excision, suggest that medical care is just sought late raising costs in specialized treatment. The elevated incidence of HPV highlights the importance of policies to encourage vaccination to control viral infection in the male population.

Conclusions

We reported an incidence of 53.2% of hrHPV infection in men with penile SCC in the State of Ceará, Northeast Brazil, a region with the lowest per capita income in the country. Despite the increased incidence, HPV infection was not associated with poor prognosis, such as systemic or lymph node metastasis, locoregional recurrence, but it seems to increase the survival rate. The current research also indicates, for the first time, that increased global DNA 5mC and reduced 5hmC marks are characteristics of penile SCC. Despite no statistical difference, increased 5mC may contribute to poor prognosis as the HR was 1.06 in SCC. Although viral infection contributes to the loss of aberrant DNA reprogramming, the 5mC/5hmC ratio was not related to HPV infection, and 5hmC levels were increased in p16INK4a negative samples. This information should be further explored as this data may predict potential clinical relevance for penile cancer prognosis and may suggest new treatment strategies, as hypomethylating agents for epigenetic targeted therapies (epidrugs). The poor prognosis of penile cancer and its relationship with environmentally induced changes reinforces the benefits of primary health care and the elevated incidence of hrHPV highlights the importance of vaccination and continued educational strategies for both prevention and treatment of malignant lesions.

Methods

Participants and study design

This is a retrospective study that included 224 participants with penile SCC who underwent partial or total penectomy, and advantage stage with enlarged prostatectomy and emasculation, without any prior history of chemotherapy or radiotherapy. Penile SCC samples were obtained at Hospital Haroldo Juaçaba, Ceará, Northeast Brazil, and detailed clinicopathological and follow-up data were assessed from 2000 to 2018. The study was approved by our Institutional Review Board (process number 2.427.846).

Available epidemiological data for penile SCC were obtained. The anatomopathological evaluation was performed by two pathologists who were blinded to the clinical data, and the regions corresponding to the neoplasm were marked on the histological slides. Pathological staging was performed according to the eighth edition of the American Joint Committee on Cancer (AJCC) (2017) [36]. Local recurrence in the amputation stump, as well as lymph node and systemic metastasis, was evaluated. The mean follow-up was 29.69 (37.3) months from surgery to the last visit or death. As we had different types of block quality and the amount of material had to be fractionated for pathological analysis and p16INK4a, hrHPV, 5mC, and 5hmC markings, the number of samples available for analysis varied (Table 1 and Additional file 1: Supplementary Fig. 1).

Tissue microarray construction

The Tissue MicroArray (TMA) blocks were composed of representative 2.0 mm cores, in duplicate, from the 224 samples from the epidemiological study cases. Blocks were sectioned to 4 μm thick and mounted on glass slides with an organosilane-based adhesive (3-aminopropyltriethoxy-silane; Sigma Chemical Co®, St Louis, MO, USA). TMA sections were then used for staining in HPV detection assays and analyses of 5mC/5hmC epigenetic markers. All the experiments were analyzed by a blinded observer.

p16INK4a immunoexpression and HPV in situ hybridization

The p16INK4a immunoexpression was carried out using an anti-p16 antibody clone E6H4 (Roche, USA). p16INK4a expression must be ≥ 75% in neoplastic cells to be considered positive in immunohistochemistry (IHC) staining, with continuous and complete cytoplasmic and nuclear staining [37]. Chromogenic in situ hybridization (CISH) was used to distinguish hrHPV and low-risk HPV (lrHPV) forms. For that, Ventana Inform HPV III Family 16 Probe (Ventana Medical Systems, Tucson, AZ) diagnostic kit was used to identify hrHPV (genotypes 16, 18, 31, 33, 35, 39, 45, 51, 52, 56, 58 and 66) and Ventana Inform HPV II Family 6 (Ventana Medical Systems, Tucson, AZ) for lrHPV identification (genotypes 6 and 11). High-grade cervical intraepithelial neoplasia was used as a positive control for hrHPV and condyloma acuminata for lrHPV. Skeletal striated muscle was used as a negative control for both assays. Representative images of the positive and negative markings of the tumor area for each assay used in this work are available in Additional file 1: Supplementary Fig. 2.

Global DNA methylation and hydroxymethylation assays

The TMA sections were stained with anti-5-methylcytosine (anti-5mC) and anti-5-hydroxymethylcytosine (anti-5hmC) antibodies using streptoavidin-biotin IHC. Briefly, after deparaffinization and rehydration, TMA slides were submitted to antigenic recovery in citrate pH 6.0 0.01 M. The slides were, then, blocked by endogenous peroxidase in a 3% hydrogen peroxide solution for 30 min, to block unspecified staining. Subsequently, slides were incubated with primary antibodies (anti-5mC: Abcam–Ab214727ug50; anti-5hmC: Abcam–Ab214728ug50) in a humid and dark chamber at 4 °C overnight. Secondary incubation took place with the DAKO EnVision™ kit (DAKO®, Carpentaria, CA, USA) for 30 min at room temperature and counterstained with Mayer’s hematoxylin. The analysis of 5mC and 5hmC was performed by evaluating the percentage of positivity of 100 cells in the tumor area of each TMA sample, in duplicate, for each target. Counting was made by two independent observers and the percent positivity rate was used to construct the score (% score). As each sample had 5mC and 5hmC scores, the 5mC/5hmC ratio was calculated for all participants. We used brain samples as positive control and the absence of the primary antibody during the incubation as the negative control.

Statistical analyses

Qualitative variables were summarized considering absolute and relative frequencies. The Fisher’s exact test was applied to verify which qualitative variables are associated with p16INK4a and hrHPV and the Wilcoxon nonparametric test for independent variables was applied to verify the relation between quantitative variables and p16INK4a and hrHPV. A multivariate regression analysis of survival data was based on the Cox proportional hazards model using PHREG procedure. The characteristics with reduced sample size as staging (pT4) were not included in the model. The Student's t-test was used to compare 5mC and 5hmC marks and HPV infection. All statistical analyses were performed using SAS® 9.3 software (SAS Institute Inc, University of North Carolina, North Carolina), and P < 0.05 was considered significant. The Kaplan–Meier curves were predicted by log-rank (Mantel-Cox) test performed in GraphPad Prism version 8.4.2. For survival curves, the sample size evaluated to time to death is presented in Table 1.

Availability of data and materials

Not applicable.

References

  1. Sung H, Ferlay J, Siegel RL, Laversanne M, Soerjomataram I, Jemal A, et al. Global cancer statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 Countries. CA Cancer J Clin. 2021;71:209–49.

    Article  PubMed  Google Scholar 

  2. de Souza MAC, Zacchi SR, Viana KCG, de Souza CB, Zandonade E, Amorim MHC. Survival analysis of penile cancer patients treated at a tertiary oncology hospital. Cienc e Saude Coletiva. 2018;23:2479–86.

    Article  Google Scholar 

  3. Eduardo C, Cardona M, García-perdomo HA. Incidence of penile cancer worldwide: systematic review and meta-analysis. Pan Am J Public Heal. 2017;2:1–10.

    Google Scholar 

  4. Bandini M, Ahmed M, Basile G, Watkin N, Master V, Zhu Y, et al. A global approach to improving penile cancer care. Nat Rev Urol. 2022;19:231–9.

    Article  PubMed  Google Scholar 

  5. Brito HO, de Calixto JRR, Medeiros R, da Costa RMG. Comment on DKK1 inhibits canonical Wnt signaling in human papillomavirus-positive penile cancer cells. Transl Oncol. 2022;16:2021–2.

    Article  Google Scholar 

  6. Maddineni SB, Lau MM, Sangar VK. Identifying the needs of penile cancer sufferers: A systematic review of the quality of life, psychosexual and psychosocial literature in penile cancer. BMC Urol. 2009;9:1–6.

    Article  Google Scholar 

  7. Douglawi A, Masterson TA. Penile cancer epidemiology and risk factors: a contemporary review. Curr Opin Urol. 2019;29:145–9.

    Article  PubMed  Google Scholar 

  8. Baraziol R, Schiavon M, Fraccalanza E, De Giorgi G. Melanoma in situ of penis: a very rare entity. Med (United States). 2017;96:3–6.

    Google Scholar 

  9. Bleeker MCG, Heideman DAM, Snijders PJF, Horenblas S, Dillner J, Meijer CJLM. Penile cancer: epidemiology, pathogenesis and prevention. World J Urol. 2009;27:141–50.

    Article  PubMed  CAS  Google Scholar 

  10. Christodoulidou M, Sahdev V, Houssein S, Muneer A. Epidemiology of penile cancer. Curr Probl Cancer [Internet]. 2015;39:126–36. https://doi.org/10.1016/j.currproblcancer.2015.03.010.

    Article  PubMed  Google Scholar 

  11. de Júnior PFM, Silva EHV, Moura KL, de Aquino YF, Weller M. Increased risk of penile cancer among men working in agriculture. Asian Pacific J Cancer Prev. 2018;19:237–41.

    Google Scholar 

  12. de Zequi SC, Guimarães GC, da Fonseca FP, Ferreira U, de Matheus WE, Reis LO, et al. Sex with animals (SWA): behavioral characteristics and possible association with penile cancer. A multicenter study. J Sex Med. 2012;9:1860–7.

    Article  Google Scholar 

  13. Vieira CB, Feitoza L, Pinho J, Teixeira-Júnior A, Lages J, Calixto J, et al. Profile of patients with penile cancer in the region with the highest worldwide incidence. Sci Rep. 2020;10:1–7.

    Article  Google Scholar 

  14. Archier E, Devaux S, Castela E, Gallini A, Aubin F, Le Maître M, et al. Carcinogenic risks of Psoralen UV-A therapy and Narrowband UV-B therapy in chronic plaque psoriasis: a systematic literature review. J Eur Acad Dermatol Venereol. 2012;26:22–31.

    Article  PubMed  CAS  Google Scholar 

  15. Daling JR, Madeleine MM, Johnson LG, Schwartz SM, Shera KA, Wurscher MA, et al. Penile cancer: importance of circumcision, human papillomavirus and smoking in in situ and invasive disease. Int J Cancer. 2005;116:606–16.

    Article  PubMed  CAS  Google Scholar 

  16. Barnes KT, McDowell BD, Button A, Smith BJ, Lynch CF, Gupta A. Obesity is associated with increased risk of invasive penile cancer. BMC Urol [Internet]. 2016;16:7–10. https://doi.org/10.1186/s12894-016-0161-7.

    Article  Google Scholar 

  17. Schoen EJ, Oehrli M, Colby C, Machin G. The highly protective effect of newborn circumcision against invasive penile cancer. Pediatrics. 2000;2:105.

    Google Scholar 

  18. Muñoz N, Bosch FX, de Sanjosé S, Herrero R, Castellsagué X, Shah KV, et al. Epidemiologic classification of human papillomavirus types associated with cervical cancer. N Engl J Med. 2003;348:518–27.

    Article  PubMed  Google Scholar 

  19. Yanagawa N, Osakabe M, Hayashi M, Tamura G, Motoyama T. Detection of HPV-DNA, p53 alterations, and methylation in penile squamous cell carcinoma in Japanese men. Pathol Int. 2008;58:477–82.

    Article  PubMed  CAS  Google Scholar 

  20. Ferreux E, Lont AP, Horenblas S, Gallee MPW, Raaphorst FM, von Knebel DM, et al. Evidence for at least three alternative mechanisms targeting the p16INK4A/cyclin D/Rb pathway in penile carcinoma, one of which is mediated by high-risk human papillomavirus. J Pathol. 2003;201:109–18.

    Article  PubMed  CAS  Google Scholar 

  21. Gu W, Zhang P, Zhang G, Zhou J, Ding X, Wang Q, et al. Importance of HPV in Chinese penile cancer: a contemporary multicenter study. Front Oncol. 2020;10:1–8.

    Article  Google Scholar 

  22. Takamoto D, Kawahara T, Kasuga J, Sasaki T, Yao M, Yumura Y, et al. The analysis of human papillomavirus DNA in penile cancer tissue by in situ hybridization. Oncol Lett. 2018;15:8102–6.

    PubMed  PubMed Central  Google Scholar 

  23. Blomberg M, Friis S, Munk C, Bautz A, Kjaer SK. Genital warts and risk of cancer: a danish study of nearly 50000 patients with genital warts. J Infect Dis. 2012;205:1544–53.

    Article  PubMed  Google Scholar 

  24. Lont AP, Kroon BK, Horenblas S, Gallee MPW, Berkhof J, Meijer CJLM, et al. Presence of high-risk human papillomavirus DNA in penile carcinoma predicts favorable outcome in survival. Int J Cancer. 2006;119:1078–81.

    Article  PubMed  CAS  Google Scholar 

  25. Wang B, Gu W, Wan F, Wei Y, Xiao W, Lu X, et al. Prognosis of the 8th TNM staging system for penile cancer and refinement of prognostication by incorporating high risk human papillomavirus status. J Urol. 2019;203:562–9.

    Article  PubMed  Google Scholar 

  26. Sano D, Oridate N. The molecular mechanism of human papillomavirus-induced carcinogenesis in head and neck squamous cell carcinoma. Int J Clin Oncol. 2016;21:819–26.

    Article  PubMed  CAS  Google Scholar 

  27. Melsheimer P, Vinokurova S, Wentzensen N, Bastert G, Von Knebel DM. DNA aneuploidy and integration of human papillomavirus type 16 E6/E7 oncogenes in intraepithelial neoplasia and invasive squamous cell carcinoma of the cervix uteri. Clin Cancer Res. 2004;10:3059–63.

    Article  PubMed  CAS  Google Scholar 

  28. McLaughlin-Drubin ME, Crum CP, Münger K. Human papillomavirus E7 oncoprotein induces KDM6A and KDM6B histone demethylase expression and causes epigenetic reprogramming. Proc Natl Acad Sci U S A. 2011;108:2130–5.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  29. Bernat A, Avvakumov N, Mymryk JS, Banks L. Interaction between the HPV E7 oncoprotein and the transcriptional coactivator p300. Oncogene. 2003;22:7871–81.

    Article  PubMed  Google Scholar 

  30. Sen P, Ganguly P, Ganguly N. Modulation of DNA methylation by human papillomavirus E6 and E7 oncoproteins in cervical cancer. Oncol Lett. 2018;15:11–22.

    PubMed  Google Scholar 

  31. Lai P, Chun CLAY, Wing PT, Tsun YT, Ngai N, Chiu T, Huang Y. HPV-16 E6 upregulation of DNMT1 through repression of tumor suppressor p53. Oncol Rep. 2010;24:1599–604.

    Google Scholar 

  32. Hanahan D. Hallmarks of cancer: new dimensions. Cancer Discov. 2022;12:31–46.

    Article  PubMed  CAS  Google Scholar 

  33. Steenbergen RDM, Snijders PJF, Heideman DAM, Meijer CJLM. Clinical implications of (epi)genetic changes in HPV-induced cervical precancerous lesions. Nat Rev Cancer. 2014;14:395–405.

    Article  PubMed  CAS  Google Scholar 

  34. Tahiliani M, Koh KP, Shen Y, Pastor WA, Brudno Y, Agarwal S, et al. Conversion of 5-methylcytosine to 5-hydroxymethylcytosine in mammalian DNA by MLL partner TET1. Science. 2009;324:930–5.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  35. Laird A, Thomson JP, Harrison DJ, Meehan RR. 5-hydroxymethylcytosine profiling as an indicator of cellular state. Epigenomics. 2013;5:655–69.

    Article  PubMed  CAS  Google Scholar 

  36. Amin MB, Greene FL, Edge SB, Compton CC, Gershenwald JE, Brookland RK, et al. The eighth edition AJCC cancer staging manual: continuing to build a bridge from a population-based to a more “personalized” approach to cancer staging. CA Cancer J Clin. 2017;67:93–9.

    Article  PubMed  Google Scholar 

  37. Cubilla AL, Lloveras B, Alejo M, Clavero O, Chaux A, Kasamatsu E, et al. Value of p16INK4a in the pathology of invasive penile Squamous cell carcinomas: a report of 202 cases. Am J Surg Pathol. 2011;35:253–61.

    Article  PubMed  Google Scholar 

  38. De Cássio ZS. Penile cancer: the importance of prevention. Int Braz J Urol. 2013;39:611–3.

    Article  Google Scholar 

  39. Korkes F, Rodrigues AFS, Baccaglini W, Cunha FTS, Slongo J, Spiess P, et al. Penile cancer trends and economic burden in the Brazilian public health system. Einstein (Sao Paulo). 2020;18:eAO5577.

    Article  PubMed Central  Google Scholar 

  40. Coelho RWP, Pinho JD, Moreno JS, Garbis DVEO, Do Nascimento AMT, Larges JS, et al. Penile cancer in Maranhão, Northeast Brazil: the highest incidence globally? BMC Urol. 2018;18:1–7.

    Article  Google Scholar 

  41. Koifman L, Vides AJ, Koifman N, Carvalho JP, Ornellas AA. Epidemiological aspects of penile cancer in Rio de Janeiro: evaluation of 230 cases. Int Braz J Urol. 2011;37:231–40.

    Article  PubMed  Google Scholar 

  42. Zhu Y, Gu WJ, Wang HK, Gu CY, Ye DW. Surgical treatment of primary disease for penile squamous cell carcinoma: a surveillance, epidemiology, and end results database analysis. Oncol Lett. 2015;10:84–92.

    Article  Google Scholar 

  43. Sewell J, Ranasinghe W, De Silva D, Ayres B, Ranasinghe T, Hounsome L, et al. Trends in penile cancer: a comparative study between Australia, England and Wales, and the US. Springerplus. Springer International Publishing; 2015;4.

  44. Tanaka K, Kandori S, Nitta S, Chihara I, Kojo K, Nagumo Y, et al. Characteristics of penile cancer in Japan: an analysis of nationwide hospital-based cancer registry data. Int J Urol. 2020;27:538–42.

    Article  PubMed  Google Scholar 

  45. Poetsch M, Hemmerich M, Kakies C, Kleist B, Wolf E, Vom Dorp F, et al. Alterations in the tumor suppressor gene p16 INK4A are associated with aggressive behavior of penile carcinomas. Virchows Arch. 2011;458:221–9.

    Article  PubMed  CAS  Google Scholar 

  46. De Andrade MV, Pinho JD, Júnior AALT, Nogueira LR, Silva FF, Maulen VE, et al. P16INK4a expression in patients with penile cancer. PLoS ONE. 2018;13:1–13.

    Google Scholar 

  47. Olesen TB, Sand FL, Rasmussen CL, Albieri V, Toft BG, Norrild B, et al. Prevalence of human papillomavirus DNA and p16 INK4a in penile cancer and penile intraepithelial neoplasia: a systematic review and meta-analysis. Lancet Oncol. 2019;20:145–58.

    Article  PubMed  CAS  Google Scholar 

  48. Flaherty A, Kim T, Giuliano A, Magliocco A, Hakky TS, Pagliaro LC, et al. Implications for human papillomavirus in penile cancer. Urol Oncol Semin Orig Investig. 2014;32:53.e1-53.e8. https://doi.org/10.1016/j.urolonc.2013.08.010.

    Article  Google Scholar 

  49. Xing B, Guo J, Sheng Y, Wu G, Zhao Y. Human papillomavirus-negative cervical cancer: a comprehensive review. Front Oncol. 2021;10:1–8.

    Article  Google Scholar 

  50. Purkayastha J. Multiple cutaneous metastasis from carcinoma of the penis-report of two cases. Indian J Surg Oncol. 2013;4:73–5.

    Article  PubMed  PubMed Central  Google Scholar 

  51. Leone A, Diorio GJ, Pettaway C, Master V, Spiess PE. Contemporary management of patients with penile cancer and lymph node metastasis. Nat Rev Urol. 2017;14:335–47.

    Article  PubMed  Google Scholar 

  52. da Reis AAS, de Paula LB, de Paula AAP, Saddi VA, da Cruz AD. Aspectos clínico-epidemiológicos associados ao câncer de pênis. Cien Saude Colet. 2010;15:1105–11.

    Article  PubMed  Google Scholar 

  53. Douglawi A, Masterson TA. Updates on the epidemiology and risk factors for penile cancer. Transl Androl Urol. 2017;6:785–90.

    Article  PubMed  PubMed Central  Google Scholar 

  54. Feber A, Arya M, De Winter P, Saqib M, Nigam R, Malone PR, et al. Epigenetics markers of metastasis and HPV-induced tumorigenesis in penile cancer. Clin Cancer Res. 2015;21:1196–206.

    Article  PubMed  CAS  Google Scholar 

  55. Tian Y, Lin A, Gan M, Wang H, Yu D, Lai C, et al. Global changes of 5-hydroxymethylcytosine and 5-methylcytosine from normal to tumor tissues are associated with carcinogenesis and prognosis in colorectal cancer. J Zhejiang Univ Sci B. 2017;18:747–56.

    Article  PubMed Central  CAS  Google Scholar 

  56. Storebjerg TM, Strand SH, Høyer S, Lynnerup A-S, Borre M, Ørntoft TF, et al. Dysregulation and prognostic potential of 5-methylcytosine (5mC), 5-hydroxymethylcytosine (5hmC), 5-formylcytosine (5fC), and 5-carboxylcytosine (5caC) levels in prostate cancer. Clin Epigenet. 2018;10:1–16.

    Article  Google Scholar 

  57. Azzi S, Habib WA, Netchine I. Beckwith-Wiedemann and Russell-Silver Syndromes: from new molecular insights to the comprehension of imprinting regulation. Curr Opin Endocrinol Diabetes Obes. 2014;21:30–8.

    Article  PubMed  CAS  Google Scholar 

  58. Rodríguez-Paredes M, Esteller M. Cancer epigenetics reaches mainstream oncology. Nat Med. 2011;17:330–9.

    Article  PubMed  Google Scholar 

  59. Wang Y, Hu H, Wang Q, Li Z, Zhu Y, Zhang W, et al. The level and clinical significance of 5-hydroxymethylcytosine in oral squamous cell carcinoma: an immunohistochemical study in 95 patients. Pathol Res Pract. 2017;213:969–74. https://doi.org/10.1016/j.prp.2017.04.016.

    Article  PubMed  CAS  Google Scholar 

  60. Wang J, Su Y, Tian Y, Ding Y, Wang X. Characterization of DNA hydroxymethylation profile in cervical cancer. Artif Cells Nanomed Biotechnol. 2019;47:2706–14.

    Article  PubMed  CAS  Google Scholar 

  61. Cavalcante GM, Borges DP, De Oliveira RTG, Furtado CLM, Alves APNN, Sousa AM, et al. Tissue methylation and demethylation influence translesion synthesis DNA polymerases (TLS) contributing to the genesis of chromosomal abnormalities in myelodysplastic syndrome. J Clin Pathol. 2020;2:1–9.

    Google Scholar 

  62. Xu T, Gao H. Hydroxymethylation and tumors: Can 5-hydroxymethylation be used as a marker for tumor diagnosis and treatment? Hum Genom. 2020;14:1–10.

    Article  Google Scholar 

  63. Song CX, Yin S, Ma L, Wheeler A, Chen Y, Zhang Y, et al. 5-Hydroxymethylcytosine signatures in cell-free DNA provide information about tumor types and stages. Cell Res. 2017;27:1231–42.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  64. Li W, Zhang X, Lu X, You L, Song Y, Luo Z, et al. 5-Hydroxymethylcytosine signatures in circulating cell-free DNA as diagnostic biomarkers for human cancers. Cell Res. 2017;27:1243–57. https://doi.org/10.1038/cr.2017.121.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  65. Liu S, de Medeiros MC, Fernandez EM, Zarins KR, Cavalcante RG, Qin T, et al. 5-Hydroxymethylation highlights the heterogeneity in keratinization and cell junctions in head and neck cancers. Clin Epigenet. 2020;12:1–14. https://doi.org/10.1186/s13148-020-00965-8.

    Article  CAS  Google Scholar 

  66. Tulay P, Serakinci N. The role of human papillomaviruses in cancer progression. J Cancer Metastasis Treat. 2016;2:201.

    Article  CAS  Google Scholar 

  67. Fakhry C, Westra WH, Li S, Cmelak A, Ridge JA, Pinto H, et al. Improved survival of patients with human papillomavirus-positive head and neck squamous cell carcinoma in a prospective clinical trial. J Natl Cancer Inst. 2008;100:261–9.

    Article  PubMed  CAS  Google Scholar 

  68. Miranda Furtado CL, Dos Santos Luciano MC, Da Silva Santos R, Furtado GP, Moraes MO, Pessoa C. Epidrugs: targeting epigenetic marks in cancer treatment. Epigenetics. 2019;14:1164–76. https://doi.org/10.1080/15592294.2019.1640546.

    Article  PubMed  PubMed Central  Google Scholar 

  69. Sharma S, Kelly TK, Jones PA. Epigenetics in cancer. Carcinogenesis. 2009;31:27–36.

    Article  PubMed  PubMed Central  Google Scholar 

  70. Kuasne H, de Syllos Cólus IM, Busso AF, Hernandez-Vargas H, Barros-Filho MC, Marchi FA, et al. Genome-wide methylation and transcriptome analysis in penile carcinoma: uncovering new molecular markers. Clin Epigenet. 2015;7:1–10. https://doi.org/10.1186/s13148-015-0082-4.

    Article  CAS  Google Scholar 

  71. Wu X, Zhang Y. TET-mediated active DNA demethylation: mechanism, function and beyond. Nat Rev Genet. 2017;18:517–34. https://doi.org/10.1038/nrg.2017.33.

    Article  PubMed  CAS  Google Scholar 

  72. Liao CG, Liang XH, Ke Y, Yao L, Liu M, Liu ZK, et al. Active demethylation upregulates CD147 expression promoting non-small cell lung cancer invasion and metastasis. Oncogene. 2022;41:1780–94.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  73. Romero-Garcia S, Prado-Garcia H, Carlos-Reyes A. Role of DNA methylation in the resistance to therapy in solid tumors. Front Oncol. 2020;10:1–20.

    Article  Google Scholar 

  74. Braakhuis BJM, Snijders PJF, Keune WJH, Meijer CJLM, Ruijter-Schippers HJ, Leemans CR, et al. Genetic patterns in head and neck cancers that contain or lack transcriptionally active human papillomavirus. J Natl Cancer Inst. 2004;96:998–1006.

    Article  PubMed  CAS  Google Scholar 

  75. Kida K, Terada T, Nobuhiro UWA, Omori Y, Fujii T, Tomita Y, et al. Relationship between p16 expression and prognosis in patients with oropharyngeal cancer undergoing surgery. In Vivo (Brooklyn). 2018;32:927–35.

    Article  CAS  Google Scholar 

  76. Lassen P, Primdahl H, Johansen J, Kristensen CA, Andersen E, Andersen LJ, et al. Impact of HPV-associated p16-expression on radiotherapy outcome in advanced oropharynx and non-oropharynx cancer. Radiother Oncol. 2014;113:310–6. https://doi.org/10.1016/j.radonc.2014.11.032.

    Article  PubMed  Google Scholar 

  77. Hashmi AA, Younus N, Naz S, Irfan M, Hussain Z, Shaikh ST, et al. p16 Immunohistochemical Expression in Head and Neck Squamous Cell Carcinoma: Association With Prognostic Parameters. Cureus. 2020;12.

  78. Shi Q, Xu L, Yang R, Meng Y, Qiu L. Ki-67 and P16 proteins in cervical cancer and precancerous lesions of young women and the diagnostic value for cervical cancer and precancerous lesions. Oncol Lett. 2019;18:1351–5.

    PubMed  PubMed Central  CAS  Google Scholar 

Download references

Acknowledgements

The authors are grateful to the members of the Experimental Oncology Laboratory (LOE), the Drug Research and Development Center (NPDM) at the Federal University of Ceará (UFC), the Cancer Institute of Ceará (ICC), and the University of Fortaleza (UNIFOR) for the support offered. We also would like to thank Suleimy Mazin for performing the statistical analysis and Silvana França dos Santos for laboratorial support.

Funding

Programa Pesquisa para o SUS (PPSUS-process number: 3931545/2017), Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES), Programa de Excelência Acadêmica (CAPES-PROEX) and Programa de Apoio à Pós-Graduação (CAPES-PROAP), Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq, Grant Numbers: 437037/2018-5 (C.L.M.F.) and 434821/2018-7 (C.P.), besides the research sponsorships of C. P. (PQ-1B, Process n°: 303102/2013-6).

Author information

Authors and Affiliations

Authors

Contributions

C.G.H., R.S.S. and D.P.P provided epidemiological data. I.J.L.S.F contributed to organization of the penile cancer sample bank. C.G.H. and R.S.S performed immunohistochemical analysis of p16INK4a, hrHPV and methylation. D.S.P. and A.P.N.N.A contributed to TMA and slides assembly and immunohistochemistry method. C.P., C.L.M.F, A.A.A.M., M.J.B.B., M.V.A.L. and M.O.M.F.C performed project development and fund-raising. C.L.M.F and R.S.S contributed to study design and wrote the manuscript. All authors have approved the final version.

Corresponding author

Correspondence to Cristiana Libardi Miranda Furtado.

Ethics declarations

Ethics approval and consent to participate

This study was approved by our Institutional Review Board (process number 2.427.846). The samples are kept in the custody of the Cancer Institute of Ceará Biobank, as well as the data from the medical records and the histopathological blocks. All patients agreed to the Informed Consent Form.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

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

. Study flowchart, immunohistochemistry and survival rate for penile squamous cell carcinoma.

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 http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) 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

Santos, R.S., Hirth, C.G., Pinheiro, D.P. et al. HPV infection and 5mC/5hmC epigenetic markers in penile squamous cell carcinoma: new insights into prognostics. Clin Epigenet 14, 133 (2022). https://doi.org/10.1186/s13148-022-01360-1

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13148-022-01360-1

Keywords