Figures
Abstract
Background
In Brazil, penile cancer (PC) is not uncommon. The highest incidence of PC is in the North and Northeast of the country. In addition to phimosis, the Human Papillomavirus (HPV) and Epstein-Baar Virus (EBV) infections are also related as risk factors for PC. The overexpression of p16INK4a is a surrogate sensitive marker of HPV infection in PC.
Objectives
To correlate p16INK4a overexpression and HPV infection status with EBV infection in a series of PC patients from the Amazon region.
Methods
Tumor tissues from 47 PC cases were analyzed for the presence of HPV and EBV DNA by PCR. All PC patients were diagnosed between 2013 and 2018 at a public reference cancer center hospital in Manaus, Amazonas—Brazil. HPV was genotyped using E7 HPV16/HPV18 type-specific real-time PCR and the PapilloCheck® HPV-Screening assay. p16INK4a expression was evaluated by immunohistochemistry using the automated Ventana® BenchMark Ultra.
Results
The mean age of patients at the time of diagnosis was 57.4 years ±SD 17.8 ranging from 20 to 90 years old. Most of the patients (64%) came from rural areas of the Amazonas State. Thirty patients had phimosis (64%). Among the patients with phimosis, 43% (13/30) underwent circumcision, three during childhood and 10 in adulthood. 60% of the patients were smokers or ex-smokers. HPV infection was observed in 45% (21/47) of cases. HPV16 was detected in 13 patients (61%). Other HPV types detected were HPV 6, 11, 42, 51, 53, 68 and 44/55. EBV infection was observed in 30% (14/47) of the patients with PC. Co-infection with HPV and EBV was observed in 28% (6/21) cases. p16INK4a was only investigated in 26 samples. The p16INK4a overexpression was observed exclusively in HPV 16 positive cases and four HPV negative cases. In the survival analysis, the follow-up time was 35.4 months/patient. The mortality rate during the follow up time was 38%.
Citation: Martins VdCA, Cunha IW, Figliuolo G, Rondon HHdMF, de Souza PM, Torres Silva FL, et al. (2020) Presence of HPV with overexpression of p16INK4a protein and EBV infection in penile cancer—A series of cases from Brazil Amazon. PLoS ONE 15(5): e0232474. https://doi.org/10.1371/journal.pone.0232474
Editor: Marc O. Siegel, George Washington University, UNITED STATES
Received: December 18, 2019; Accepted: April 15, 2020; Published: May 6, 2020
Copyright: © 2020 Martins et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Data Availability: All relevant data are in the paper and in the Supporting information.
Funding: This work was supported by The Fundação de Amparo à Pesquisa do Estado do Amazonas- FAPEAM for financial support (PROGRAMA UNIVERSAL AMAZONAS N.021/2011-FAPEAM; PROGRAMA PAPAC N.020/2013 – FAPEAM; PROGRAMA PAPAC N.05/2019 – FAPEAM and PROGRAMA PRó ESTADO N.002./2008-FAPEAM). The Rede D’OR- São Luiz provided support in the form of a salary for author Isabela Werneck Cunha IW, but did not have any additional role in the study design, data collection and analysis, decision to publish, or preparation of the manuscript. The specific roles of this author is articulated in the ‘author contributions’ section. The Genomic Health Surveillance Network: Optimization of Assistance and Research in the State of Amazonas – REGESAM is not a funder but a network of researchers in the Genomic field on the State of Amazonas.
Competing interests: The authors have read the journal’s policy and have the following conflicts: author Isabela Werneck Cunha IW is affiliated with The Rede D’OR- São Luiz. This does not alter our adherence to all the PLOS ONE policies on sharing data and materials.
Introduction
The occurrence of Penile Cancer (PC) varies worldwide. In developed countries, PC has a low incidence, corresponding to 0.3–1% of malignant neoplasms in men. In some developing countries, the incidence of PC may be much higher than the global average incidence [1–3]. In Brazil, PC accounts for approximately 2.1% of all tumors in men being the highest incidence reported in the Latin America (2.9–6.8 cases per 100,000 men-years) [3]. We can find regional rate differences along the country. In the North and Northeast of Brazil, PC incidence is five times higher compared to the Midwest, South and Southeast regions [4]. According to the Brazilian National Cancer Institute (INCA), the death rate from PC in the northern region of Brazil has doubled in the last decade, from 0.05% to 0.10% [5].
PC is a multifactorial disease and the risk factors and/or favorable conditions to develop PC are not fully established. Phimosis and smegma accumulation are observed in more than 80% of the patients with PC associated with chronic inflammation process [6,7]. Other factors, such as smoking and sexually transmitted infection, are also related to the onset of neoplasms [7,8]. HPV infection is present in approximately 50% of PC cases and the most prevalent genotype is HPV16 (30%) [7,9–11]. A recent meta-analysis study of 52 studies showed a pooled prevalence of 50.8% (44.8–56.7) of HPV infection in PC with a rate of 68.3% (58.9–77.1) of HPV16 [12].
The role played by HPV in carcinogenesis of the penis appears to be similar to cervical cancer. HPV encodes the E6 and E7 oncogenes which are required for malignant transformation and maintenance of host cells. The viral oncoproteins (E6 and E7) may compromise the regulation of the host cell cycle and lead to an uncontrolled proliferation [13,14]. P16 is a tumor suppressor gene and its protein is physiologically expressed in normal tissues. The inactivation of the retinoblastoma gene (pRb) by HPV E7 results in overexpression of p16INK4a due to the lack of negative feedback loop between pRb and p16 protein [15]. The overexpression of p16INK4a in tumor cells has been shown to correlate with high-risk HPV DNA detection in PC [16].
Epstein-Barr virus (EBV) is another agent associated with PC [17,18]. Inappropriate expression of its latent genes (Latent Membrane Protein) LMP-1, LMP-2A e LMP-2B, involved in cell persistence, may contribute to the development of tumors [19]. EBV is suggested as a viral cofactor rather than a primary carcinogen in Burkitt’s lymphoma, Hodgkin’s disease, nasopharyngeal carcinoma [17,20]. In HPV-associated cancer, the presence of EBV may also act as a viral cofactor [21]. Several studies have shown the presence of EBV and co-infection with HPV in PC but a relationship between PC and EBV is yet to be established [17,18].
In light of the possible roles of HPV and EBV infection in the development of PC, this study investigated p16INK4a expression and HPV and EBV infection in a series of patients with PC from the Brazilian Amazon region.
Materials and methods
Enrollment
A total of 47 patients with PC and no concomitant urological neoplastic diseases participating in the study were attended at the public reference cancer center hospital—Fundação Centro de Controle de Oncologia do Estado do Amazonas/FCECON from 2013 to 2018. The patients were followed at the Urology clinic of the hospital. All patients were surgically treated by total or partial penectomy.
All patients provided written informed consent and were interviewed to fill a questionnaire concerning sociodemographic data and risk factors. Histopathology characteristics of the tumors were obtained from the medical charts.
This study was approved by the internal review board of the Ethics Committee of the FCECON—approval document #2.230.007, August 21, 2017. SISGEN- A5F36C5.
Biological samples collection
At the moment of surgery, three to five mm3 of tissue fragments (mass of 50-150mg) from the tumor were collected and stored in a dry plastic microtube free of DNAse and RNAse. Samples were stored at -30 °C until processed.
DNA extraction
DNA was extracted from the frozen tissue (mass of 20-40mg from different parts of the tumor) using the DNeasy® Blood & Tissue Kit (QIAGEN Inc., USA), according to the manufacturer’s recommendations. The DNA was eluted in a volume of 200 μL UltraPure™ DNase/RNase-Free Distilled Water (Invitrogen Life Technologies, São Paulo, Brazil).
Human β-globin PCR
For DNA extraction quality control, the human β-globin gene was amplified by PCR as previously described in the literature with the following pair of primers: GH20: (5’GAAGAGCCAAGGACAGGTAC’3) and PCO4: (5’CAACTTCATCCACGTTCACC’3) generating a 270 bp DNA fragment [22].
HPV detection
All samples were submitted to generic HPV PCR using the consensus primers (PGMY09/11) which amplifies a 450 bp DNA fragment within the L1 region of mucosal HPVs [23]. Amplification was carried out as previously described using 50–100 ng of DNA in 25 μL of reaction mixture and a thermocycling profile of 1 cycle at 5 min at 95 °C, followed by 40 cycles: 1 min at 95 °C, 1 min at 55°C, and 1 min at 72 °C, with a final extension for 10 min at 72 °C. The PCR products (450 bp DNA) were analyzed on 1.5% agarose gel stained with SYBR™ Safe DNA Gel Stain (Invitrogen Life Technologies, São Paulo, Brazil) for visualization of DNA under UV light and 100 bp DNA ladder was used as molecular weight control pattern. Precautions to avoid contamination were followed. DNA from the HeLa cell line which harbors 10–20 copies of integrated HPV 18 per cell was used as a positive control in all reactions.
E7 HPV16/HPV18 type-specific real-time PCR.
All samples were also submitted to two specific TaqMan based real-time qPCR assays targeting either HPV16/HPV18 E7 gene in an ABI 7300 Real-Time PCR System (Applied Biosystems, Foster City, CA). All samples and controls were run in duplicate.
HPV16-E7.
qPCR assay included the following primers: forward (5’GATGAAATAGATGGTCCAGC3’) and reverse (5’GCTTTGTACGCACAACCGAAGC3’) primers, and the probe (5’FAM-CAAGCAGAACCGGACAG-MGB-NFQ) in a final reaction volume of 25 μL [24]. Each qPCR reaction contained 1X TaqMan master mix (Applied Biosystems, Foster City, CA), 400 nM each of the forward and reverse primers, 200 nM of fluorogenic TaqMan probe, and 50–100 ng of DNA. The amplification conditions consisted of 50 °C for 2 min and 95°C for 10 min, followed by 40 cycles of 95 °C for 15 sec, 55 °C for 1 min, and 60 °C for 1 minute. DNA from a SiHa cell line which contains 1–2 copies of integrated HPV 16 per cell was used as a positive control in all reactions.
HPV18-E7.
qPCR assay included the following primers: forward (5’AAGAAAACGATGAAATAGATGGA3’) and reverse (5’GGCTTCCACCTTACAACACA3’) primers, and a probe (5’VIC-AATCATCAACATTTACCAGCC-MGBNFQ3’) in a final reaction volume of 25 μL, each qPCR contained 1X TaqMan master mix (Applied Biosystems, Foster City, CA), 400 nM each of the forward and reverse primers, 400 nM fluorogenic TaqMan probe, and 50–100 ng of DNA [24]. The amplification conditions consisted of 50°C for 2 min and 95°C for 10 min, followed by 40 cycles of 95°C for 15 sec, 50°C for 1 min, and 60°C for 1 min. DNA from the HeLa cell line which harbors 10–20 copies of integrated HPV 18 per cell was used as a positive control in all reactions.
HPV genotyping—PapilloCheck® HPV-Screening.
All samples that were positive in generic HPV DNA (PGMY09/11) and negative for 16/18 genotypes were submitted to HPV-Screening Test (Greiner Bio-One GmbH, Frickenhausen, Germany) to identify other genotypes. This is a PCR-based DNA microarray system for detection and identification of 24 HPV genotypes, including 16 high-risk HPV genotypes (HPV 16, 18, 31, 33, 35, 39, 45, 51, 52, 56,58, 59, 68, 70, 73, 82), 2 probable high-risk HPV genotypes (HPV 53,66) and 6 low-risk-HPV genotypes (HPV 6, 11, 40, 42, 43, 44/55) [25].
EBV detection
EBV DNA detection was performed as described elsewhere [26]. Briefly, a sensitive multiplex PCR which amplifies 182 bp within the Exons 4/5 from the terminal protein RNA of EBV and a fragment of human β-actin 450 bp as internal control. The following pairs of primers for PCR amplification were used: EP5-AACATTGGCAGCAGGTAAGC and EM3 -ACTTACCAAGTGTCCATAGGAGC for EBV and B-ACT F–TCTACAATGAGCTGCGTGTG and B-ACT R -CATCTCTTGCTCGAAGTC for β-actin. PCR was performed as follows: 1 cycle at 5 min at 95°C, followed by 10 cycles of 30 sec at 95°C, 60 sec at 63°C and subsequently by 30 cycles of 30 sec at 95°C, 30 sec at 60°C, and 30 sec at 72°C, with a final extension for 40 sec at 72°C. PCR products were analyzed on 1.8% agarose gel stained with SYBR™ Safe DNA Gel Stain (Invitrogen Life Technologies, São Paulo, Brazil) for visualization of DNA under UV light and 100 bp DNA ladder was used as molecular weight control pattern. An EBV positive known sample from the laboratory was used as control.
Immunohistochemistry for p16INK4a
Twenty six samples from PC patients, spotted in tissue microarray were submitted to the immunohistochemistry (IHC) assay for qualitative detection of the p16INK4a (21 samples were not available due to pre-analytical factors). IHC assay was performed in an automated system using the Ventana® BenchMark Ultra according to the manufacturer’s instructions. The IHC slides were analyzed by a pathologist. Positivity for p16INK4a was defined as unequivocally nuclear and cytoplasmic staining of at least 70% of the tumor cells [27,28].
Statistical analysis
Data were compiled using Epi Info™ version: 7.2.2.2 and analyzed using Stata application (v.13, StataCorp, 2013, College Station, Texas, USA). Descriptive analyses were performed using frequency tables. In this study, we analyze independently the relationship between the sociodemographic and clinical profile of patients with PC and three different outcomes: (i) p16INK4a overexpression, (ii) HPV infection status and (iii) EBV infection status. Univariate logistic regression models were used to analyze the relationship between these outcomes and sociodemographic and clinical profile of patients with PC and calculate the odds ratios. A p-value <0.05 was considered statistically significant.
Those variables that presented an association at a level of significance of 0.2 in the univariate logistic regression were selected to perform the multiple logistic regression. Multiple logistic regression models were used to identify the independent relationships between clinical characteristics of the patients and the studied outcomes. In this manner, adjusted odds ratios were calculated. A p-value <0.05 was considered statistically significant.
Kaplan-Meier method was used to determine the differences between survival time of patients stratified by HPV and EBV infection, and with p16INK4a overexpression. Groups were compared using the log rank test, taking a p-value less than 0.05 as statistically significant.
Results
HPV and EBV were investigated in tumor tissues from 47 patients with PC, aged 20 to 90 years old (mean 57.4 years ±SD 17.8). Ten patients were below the age of 39 (21%). 64% of the patients came from rural areas of the Amazonas State, 23% from Manaus, the capital city and 13% from other Northern states of Brazil. 79% of the patients with PC were Amerindians descendent (mestizo) and 83% had less than eight years of schooling. 32% of the patients had total penectomy and 68% had partial penectomy. Lymphadenectomy was performed in 13/47 (28%) patients (Table 1).
All cases were diagnosed as squamous cell carcinoma. Of the 47 patients with PC, we were able to review and reclassify only 20 patients according to the 2016 WHO classification [29]. The most common histology subtype was “usual” in 11/20 (55%). We could obtain tumor characteristic from the medical records for only 43 patients. 58% of tumors were located in the glans and foreskin and the predominant pattern of growth was verruciform (48%). Other clinical and pathological parameters can be appreciated in Table 1.
In this study, 36% of the patient had a history of cancer in the family. 60% of the patients were smokers or ex-smokers and 64% had phimosis. 43% were circumcised, mostly in adulthood. 28% of the patients reported a history of sexually transmitted diseases once or several times in their adulthood (Table 2).
The prevalence of HPV is shown in Table 3. HPV DNA was detected in 45% (21/47) of the patients with PC. HPV 16 was the most prevalent genotype 61% (13/21). The other 8 genotypes detected were HPV 6, 11, 42, 51, 53, 68 and 44/55. One patient had multiple co-infections; genotypes 16, 42, 44/55. The distribution of EBV and HPV among the patients with PC is shown in Fig 1. EBV infection was observed in 14 out of 47 patients (30%). Co-infection HPV/EBV was detected in 6 patients, four with HPV 16/EBV and two with HPV 6/EBV and HPV 53/EBV, each. The p16INK4a overexpression was observed in three patients with co-infection HPV/EBV, all HPV 16.
EBV: Epstein-Barr virus; HPV: Human papillomavirus.
Association of p16INK4a overexpression with clinical factors and HPV infection status are shown in Table 4. Overexpression of p16INK4a was found in 12 cases 46% (12/26). Patients with phimosis had 11 times more chance of having overexpression of p16INK4a [OR = 11 (95%CI 1.1–109.7); p = 0.04]. The p16INK4a overexpression was observed in the eight HPV 16 positive cases and in four HPV-. Other clinical factors were not significantly associated with HPV infection and p16INK4a overexpression. HPV infection and p16INK4a overexpression were also related to some histological subtypes, HPV infection being positive in the basaloid subtypes (45%) and negative in usual subtypes (81%) (Table 4). Regarding p16INK4a overexpression, of the 11 nonsmokers patients, eight presented over expression of p16INK4a while only four patients of the 15 smokers were with p16INK4a overexpression [OR = 0.13 (95%CI 0.02–0.780); p = 0.026]. Notably, of the eight nonsmoker’s patients with p16INK4a overexpression, six had HPV 16 genotypes and the remaining two were negative for HPV. Among the four smoker’s patients with p16INK4a overexpression, two were HPV16 and two negative for HPV. HPV+ patients with PC had poorly differentiated carcinomas (Grade III) compared to HPV- patients with PC [OR = 0.07 95%CI 0.01–0.047; p = 0.005].
In the survival analysis, the follow-up time was 35.4 months/patient. The mortality rate during the study period was 38% (18/47). Survival analysis was performed by stratification of the patients into p16INK4a+ vs p16INK4a-, HPV+ vs HPV-, EBV+ vs EBV- and HPV/EBV co-infection vs no infection. Infection status of the deceased and survival patients is shown in S1 Table. There is no evidence of differences in survival of patients according to p16 overexpression (Plog rank = 0.753), HPV infection (Plog rank = 0.979), EBV infection (Plog rank = 0.106) and HPV/EBV co-infection (Plog rank = 0.318) (Fig 2).
- P16 expression and overall mortality
- HPV status and overall mortality
- EBV status and overall mortality
- HPV/EBV co-infection vs no infection and overall mortality
Overall survival depending on p16INK4a expression (A), HPV infection status (B), positive and negative for EBV (C) and HPV/EBV co-infection (D). EBV: Epstein-Barr virus; HPV: Human papillomavirus.
Discussion
In the Northern region of Brazil, the mortality rate from PC has doubled in recent years according to INCA [5]. The state of Amazonas, the largest state in Brazil, covers a geographical area of 1,559,168.12 km2 and has about 4 million inhabitants. Half of the population lives in the capital city, Manaus while the other half is distributed irregularly in the rural regions with low density population with very poor access to health services [31].
The disease affects mainly men between the fifth and seventh decade of life [2,6,32]. In this study, the mean age of the diagnosis of PC was 57.4 years, with a higher prevalence in the age group below 59 years. However, young adults less than 30 years were also diagnosed with PC (9%). Other studies in the country have also reported the early occurrence of PC [4,33–37]. The risk factors and carcinogenesis of PC among young adults is still not established. In this study, young patients were at an advanced stage of the disease. Half of them underwent total penectomy and had metastasis. Of the four young patients, one had co-infection HPV 6/EBV and one HPV 44/55. Two were HPV negative. Interestingly, none of them had high-risk HPV. This may suggest that there are other factors or genetic and molecular changes possibly involved in the development of PC in young adults.
In this study, most of the patients came from the interior of Amazonas and belong to riverside isolated populations with low levels of education and poor access to health services. Besides, they have poor knowledge about the disease and seek diagnosis at an advanced stage. This reality was also observed by Chalya et al. (2015) in Tanzania [38] and in Brazil by other groups [33–35,39,40].
PC usually begins with a superficial or ulcerated lesion on the glans and foreskin, but can also spread through the penile shaft and the scrotum [41]. In the current study, the initial lesion was diagnosed in the glans and/or prepuce in 76% of the cases similar to other studies [4,38,42]. Lesions were predominantly of the verruciform type followed by warts. The accumulation of smegma due to poor hygiene followed by probable irritation of the local is a favorable environment for several infections. Chronic inflammatory processes may progress to the development of lesions and if untreated, can lead to neoplasm [7].
In the present study, 70% (23/33) of the patient showed well and moderately differentiated tumors (Grade I/II) and was not associated with HPV infection (adjusted p-value = 0.006). Patients with HPV presented higher grade tumor. Histopathology grade is an important prognostic factor. Degree of poorly differentiated cell could indicate a worse prognosis of lesions. Our findings reinforced other studies observations [43,44].
Phimosis, a well-known risk factor for the development of PC, was present in 64% of the patients. 43% were submitted to circumcision during adulthood. Adult circumcision is known to have no protective effect against the development of PC [45]. The univariate logistic regression showed that patients with phimosis had 11 times more chance of overexpression of p16INK4a [OR = 11 (95%CI 1.1–109.7); p = 0.04]. Among men circumcised in adulthood, phimosis was strongly associated with development of invasive penile cancer, in concordance with the findings of Daling, et al. (2005) [46].
Of note, the presence of HPV in an individual does not mean that the individual will develop cancer. There are many risk factors that contribute to the development or not such as the environment and the genetic background of the individuals as well as the viral clearance capacity of the individual. PC can be HPV-mediated or not. Individuals HPV- with phimosis and chronic inflammation often develop PC. Genetic and molecular changes associated with HPV- PC leading to disturbance of the p14ARF/MDM2/p53 and/or p16INK4a/cyclin D/Rb pathways have been suggested as plausible mechanisms for the development of PC [47]. One study showed that there is silencing of the p16INK4a gene through promoter hypermethylation in 15% of cases and over-expression of the polycomb group (PcG) gene BMI-1, which targets the INK4A/ARF locus, encoding both p16INK4a and p14ARF, in 10% of cases. Another study have suggested that the inactivation of p14ARF/MDM2/p53 pathway as well as somatic mutation of the p53 gene and over-expression of MDM2 and mutation of p14ARF may lead to the development of PC [48,49].
Emerging interest regarding PC carcinogenesis is the association of oncogenic viruses co-infection. EBV is associated with several malignancies in humans and its involvement in PC is still controversial [50–52]. In this study, the prevalence of EBV was 30% and co-infection with HPV was 29%. Our study differs to the one conducted in Rio de Janeiro where the prevalence of EBV in penile malignancies was 46% and co-infection with HPV was 26% [17]. In cervical cancer, EBV has been suggested as a cofactor that facilitates the integration of the HPV16 genome, contributing to the development of cancer [20,21].
The prevalence of HPV in invasive PC is approximately 45%, ranging from 30% to 75% according to the detection method, the population and type of sample analyzed [11,17,18,40,42,53–55]. High-risk HPV genotype 16 was observed in 61% of HPV+ cases, reinforcing other studies observations [9,56]. HPV 18, the second most common high-risk HPV [11], was not identified in this study. Low incidence of HPV 18 in PC has also been reported in the country [17,40,42,54]. Of note, the prevalence of HPV 18 is also low in the female population of the Amazonas region [57,58]. Interestingly, in one study of the Thailand population, only 1/65 patients with PC had HPV16 but high presence of HPV18 genotype (55%) was detected [59]. This can probably be explained due to geographical distribution of HPV genotypes.
The occurrence of viral co-infection between HPV genotypes 6, 16, 42, 44/55 and EBV was observed. One patient presented EBV/HPV 6 co-infection. Similar findings, EBV/HPV 6 had been described [18] suggesting a probable viral synergism in tumor development due to their similar tropism of epithelial cells [20,21]. The association of EBV and carcinogenesis is still to be demonstrated with EBV genome or virus gene products within the tumor cell population [60]. HPV 6 is classified as low oncogenic risk and is related to condyloma [61]. However, it is a prominent feature in infections in cases of PC and has multiple co-infections with high-risk types [11,62,63].
The overexpression of p16INK4a, a surrogate sensitive marker of HPV in PC [16], still remains to correlate with prognostics [12]. Overexpression of P16INK4a was observed in 12 cases (12/26, 46%). 66% (8/12) were from patients infected with HPV 16, reinforcing the role played by HPV16 in the oncogenic process.
Four HPV negative patients showed overexpression of p16INK4a (4/12, 33%). Bleeker et al. (2009) in a systematic review, assumes that PC would be related to a pathway mediated by HPV infection and another due to different epigenetic changes in the absence of HPV and related to chronic inflammation [47]. Understanding the molecular processes involved in the onset and progression of the disease is fundamental for the prevention and treatment of this mutilating.
In this series of cases, the mortality rate during the study was 38%. Several studies tried to identify prognostic factors to manage the selection of patients at high risk for metastases in PC [64–68]. HPV infections as well as co-infection with EBV and p16INK4a positivity were not predictive of survival of the patients with PC. Lymph node involvement is related to poor prognosis, with a 5-year survival of less than 40% [69,70]. In this study, only 39% of the patients with PC had lymphadenectomy.
Limitations of the study were the difficulties inherent to the non-recording of clinical and histopathological data in medical records. Many biopsy samples fixed in paraffin blocks were missing for the IHC assay for the qualitative detection of the p16INK4a, reducing the sample size. In the state of Amazonas, usually patients seek care at an advanced stage of the disease and as soon as they complete their surgical treatment (partial or total penectomy) they return to their city and there is no follow-up.
HPV vaccination has been shown to be effective for HPV-related cancers and inclusion of young in the immunization programs against HPV is well established [71,72]. In the Amazonas, the vaccination program for girls started in 2013 with the quadrivalent vaccine that protects against genotypes 6, 11, 16 and 18. In Brazil, the Ministry of Health only included young males (between 12 and 13 years old) in the vaccination HPV program in 2017. In the future, we expect that this action may reduce the incidence of PC and other HPV-related.
Conclusions
In summary, our results show that patients with HPV+ PC have in general low grade tumors. Overexpression of p16INK4a was correlated to the detection of HPV 16 DNA, reinforcing that it can be used as a marker to high-risk HPV genotype 16 infection as found in oropharyngeal cancers. EBV infection was observed in one-third of the patients with PC and the co-infection with HPV in a quarter. The knowledge of the etiology of penile cancer is far from definite. The individual role or synergisms of the known oncogenic viruses such as HPV and EBV at the onset of carcinogenic events are still not well defined. However, our data show the profile of these viral oncogenic infections and the reality of this neoplasm in individuals from the Brazilian Amazon.
Supporting information
S1 Protocol. HPV detection and genotyping method and EBV detection.
https://doi.org/10.1371/journal.pone.0232474.s001
(DOCX)
S2 Protocol. Immunohistochemical for p16INK4a protein.
https://doi.org/10.1371/journal.pone.0232474.s002
(DOCX)
S1 Fig. Patterns of p16 expression in penile carcinomas.
Microarray tissue block immunohistochemistry for p16INK4a (from left to right): a. absence; b. strong and diffuse cytoplasmic staining; c. Moderate and focal cytoplasmic staining; d. weak and focal cytoplasmic staining; e. absence of staining; f. strong and diffuse cytoplasmic staining.
https://doi.org/10.1371/journal.pone.0232474.s003
(TIF)
S1 Table. Distribution of HPV, EBV status and p16INK4a overexpression among the patients deceased and survived with penile cancer at Amazon—Brazil.
EBV: Epstein-Barr virus; HPV: Human papillomavirus, n: Absolute frequency; +:positive; -:negative.
https://doi.org/10.1371/journal.pone.0232474.s004
(DOCX)
S1 File. Data collection instrument—Penis cancer.
https://doi.org/10.1371/journal.pone.0232474.s005
(DOCX)
Acknowledgments
We thank the Fundação Centro de Controle de Oncologia do Estado do Amazonas (FCECON), Fundação de Medicina Tropical Doutor Heitor Vieira Dourado (FMT-HVD) and A.C. Camargo Cancer Center for the infrastructural support. The authors gratefully acknowledge Dr. Fernando Augusto Soares and Dra. Stephania Martins Bezerra (A.C. Camargo Cancer Center) who performed the Immunohistochemistry for p16INK4a, Dr. Marcel Heibel for assistance during the inclusion of patients in the study and Dra. Ana Carolina Soares de Oliveira who performed the HPV Genotyping—PapilloCheck® HPV-Screening. The authors would also like to thank the patients that participated in this study.
References
- 1. Parkin DM, Ferlay J, Curado MP, Bray F, Edwards B, et al. (2010) Fifty years of cancer incidence: CI5 I–IX. International journal of cancer 127: 2918–2927. pmid:21351270
- 2. Christodoulidou M, Sahdev V, Houssein S, Muneer A (2015) Epidemiology of penile cancer. Current problems in cancer 39: 126–136. pmid:26076979
- 3. Cardona CEM, García-Perdomo HA (2018) Incidence of penile cancer worldwide: systematic review and meta-analysis. Revista Panamericana de Salud Pública 41: e117.
- 4. Favorito LA, Nardi AC, Ronalsa M, Zequi SC, Sampaio FJ, et al. (2008) Epidemiologic study on penile cancer in Brazil. International braz j urol 34: 587–593. pmid:18986562
- 5.
BRASIL. (2019) Ministério da saúde. Instituto Nacional de Câncer José Alencar Gomes da Silva. Disponível em https://mortalidade.inca.gov.br/MortalidadeWeb/ In: 09/10/2019 A, editor.
- 6. Pow-Sang MR, Benavente V, Pow-Sang JE, Morante C, Meza L, et al. (2002) Cancer of the penis. Cancer control 9: 305–314. pmid:12228756
- 7. Calmon MF, Tasso Mota M, Vassallo J, Rahal P (2011) Penile carcinoma: risk factors and molecular alterations. The Scientific World Journal 11: 269–282. pmid:21298218
- 8. Schabath MB, Villa LL, Lazcano-Ponce E, Salmerón J, Quiterio M, et al. (2011) Smoking and human papillomavirus (HPV) infection in the HPV in Men (HIM) study. Cancer Epidemiology and Prevention Biomarkers.
- 9. Miralles-Guri C, Bruni L, Cubilla A, Castellsague X, Bosch FX, et al. (2009) HPV prevalence and type distribution in penile carcinoma. J Clin Pathol.
- 10. Stratton KL, Culkin DJ (2016) A Contemporary Review of HPV and Penile Cancer. Oncology (Williston Park, NY) 30: 245–249.
- 11. Alemany L, Cubilla A, Halec G, Kasamatsu E, Quirós B, et al. (2016) Role of human papillomavirus in penile carcinomas worldwide. Eur Urol 69: 953–961. pmid:26762611
- 12. Olesen TB, Sand FL, Rasmussen CL, Albieri V, Toft BG, et al. (2019) Prevalence of human papillomavirus DNA and p16INK4a in penile cancer and penile intraepithelial neoplasia: a systematic review and meta-analysis. The lancet oncology 20: 145–158. pmid:30573285
- 13. Schiffman M, Doorbar J, Wentzensen N, De Sanjosé S, Fakhry C, et al. (2016) Carcinogenic human papillomavirus infection. Nature Reviews Disease Primers 2: 16086. pmid:27905473
- 14. Spiess PE, Dhillon J, Baumgarten AS, Johnstone PA, Giuliano AR (2016) Pathophysiological basis of human papillomavirus in penile cancer: key to prevention and delivery of more effective therapies. CA: A Cancer Journal for Clinicians 66: 481–495.
- 15. Chipollini J, Chaing S, Azizi M, Kidd LC, Kim P, et al. (2017) Advances in Understanding of Penile Carcinogenesis: The Search for Actionable Targets. International journal of molecular sciences 18: 1777.
- 16. Tang DH, Clark PE, Giannico G, Hameed O, Chang SS, et al. (2015) Lack of P16 ink4a Over Expression in Penile Squamous Cell Carcinoma is Associated with Recurrence after Lymph Node Dissection. J Urol 193: 519–525. pmid:25261804
- 17. Afonso LA, Moyses N, Alves G, Ornellas AA, Passos MR, et al. (2012) Prevalence of human papillomavirus and Epstein-Barr virus DNA in penile cancer cases from Brazil. Memórias do Instituto Oswaldo Cruz 107: 18–23. pmid:22310531
- 18. Afonso LA, Carestiato FN, Ornellas AA, Ornellas P, Rocha WM, et al. (2017) Human papillomavirus, Epstein-Barr virus, and methylation status of p16ink4a in penile cancer. Journal of medical virology.
- 19. de Lima MAP, Rabenhorst SHB (2006) Associação do vírus Epstein-Barr (EBV) com tumores sólidos. Revista Brasileira de Cancerologia 52: 87–96.
- 20. Szostek S, Zawilinska B, Kopec J, Kosz-Vnenchak M (2009) Herpesviruses as possible cofactors in HPV-16-related oncogenesis. Acta Biochim Pol 56: 337. pmid:19499088
- 21. Guidry J, Scott R (2017) The interaction between human papillomavirus and other viruses. Virus research 231: 139–147. pmid:27826043
- 22. Saiki RK, Scharf S, Faloona F, Mullis KB, Horn GT, et al. (1985) Enzymatic amplification of b-globin genomic sequences and restriction site analysis for diagnosis of sickle cell anemia. Science 230: 1350–1354. pmid:2999980
- 23. Gravitt P, Peyton C, Alessi T, Wheeler C, Coutlee F, et al. (2000) Improved amplification of genital human papillomaviruses. Journal of clinical microbiology 38: 357–361. pmid:10618116
- 24. Veo CA, Saad SS, Fregnani JHT, Scapulatempo-Neto C, Tsunoda AT, et al. (2015) Clinical characteristics of women diagnosed with carcinoma who tested positive for cervical and anal high-risk human papillomavirus DNA and E6 RNA. Tumor Biology 36: 5399–5405. pmid:25677903
- 25. Heard I, Cuschieri K, Geraets D, Quint W, Arbyn M (2016) Clinical and analytical performance of the PapilloCheck HPV-Screening assay using the VALGENT framework. Journal of Clinical Virology 81: 6–11. pmid:27262102
- 26. Markoulatos P, Georgopoulou A, Siafakas N, Plakokefalos E, Tzanakaki G, et al. (2001) Laboratory diagnosis of common herpesvirus infections of the central nervous system by a multiplex PCR assay. Journal of clinical microbiology 39: 4426–4432. pmid:11724856
- 27. Bezerra SM, Chaux A, Ball MW, Faraj SF, Munari E, et al. (2015) Human papillomavirus infection and immunohistochemical p16INK4a expression as predictors of outcome in penile squamous cell carcinomas. Hum Pathol 46: 532–540. pmid:25661481
- 28. Cubilla AL, Lloveras B, Alejo M, Clavero O, Chaux A, et al. (2011) Value of p16INK4a in the pathology of invasive penile squamous cell carcinomas: a report of 202 cases. Am J Surg Pathol 35: 253–261. pmid:21263246
- 29. Moch H, Cubilla AL, Humphrey PA, Reuter VE, Ulbright TM (2016) The 2016 WHO classification of tumours of the urinary system and male genital organs—part A: renal, penile, and testicular tumours. Eur Urol 70: 93–105. pmid:26935559
- 30. Paner GP, Stadler WM, Hansel DE, Montironi R, Lin DW, et al. (2018) Updates in the eighth edition of the tumor-node-metastasis staging classification for urologic cancers. Eur Urol 73: 560–569. pmid:29325693
- 31.
BRASIL. (2019) Ministério da Economia. Instituto Brasileiro de Geografia e Estatística—IBGE. Disponível em:https://www.ibge.gov.br/cidades-e-estados/am.html Acesso em:09/10/209.
- 32.
Salvioni R, Necchi A, Piva L, Colecchia M, Nicolai N. Penile cancer; 2009. Elsevier. pp. 677–685.
- 33.
Fonseca AGd, Pinto JASdA, Marques MC, Drosdoski FS, Fonseca Neto LORd (2010) Estudo epidemiológico do câncer de pênis no Estado do Pará, Brasil. Revista Pan-Amazônica de Saúde 1: 85–90.
- 34.
FIGLIUOLO G, Lima S, Costa S, Silva J, Paiva C, et al. (2015) Perfil clínico-epidemiológico associado a fatores de risco de pacientes com câncer de pênis atendidos em um Hospital de Referência Oncológica em Manaus. Vol. Revista Brasileira de Oncologia Clínica 11.
- 35.
Barros ÉNd Melo MCBd (2009) Câncer de pênis: perfil sócio-demográfico e respostas emocionais à penectomia em pacientes atendidos no Serviço de Psicologia do Hospital de Câncer de Pernambuco. Revista da SBPH 12: 99–111.
- 36. Couto TCd Arruda RMB, Couto MCd Barros FD (2014) Epidemiological study of penile cancer in Pernambuco: experience of two reference centers. International braz j urol 40: 738–744. pmid:25615242
- 37. de Araujo LA, De Paula AAP, de Paula HdSC, Ramos JEP, de Oliveira BR, et al. (2018) Human papillomavirus (HPV) genotype distribution in penile carcinoma: Association with clinic pathological factors. PLoS One 13: e0199557. pmid:29949632
- 38. Chalya PL, Rambau PF, Masalu N, Simbila S (2015) Ten-year surgical experiences with penile cancer at a tertiary care hospital in northwestern Tanzania: a retrospective study of 236 patients. World journal of surgical oncology 13: 71. pmid:25889238
- 39. Barbosa AdA Júnior, Athanázio PRF, Oliveira B (1984) Câncer do pênis: estudo da sua patologia geográfica no Estado da Bahia, Brasil. Revista de saude publica 18: 429–435. pmid:6536113
- 40. de Andrade Martins V, Pinho JD, Junior AALT, Nogueira LR, Silva FF, et al. (2018) P16INK4a expression in patients with penile cancer. PLoS One 13: e0205350. pmid:30312320
- 41. Mosconi AM, Roila F, Gatta G, Theodore C (2005) Cancer of the penis. Crit Rev Oncol Hematol 53: 165–177. pmid:15661566
- 42. de Sousa IDB, Vidal FCB, Vidal JPCB, de Mello GCF, Nascimento MdDSB, et al. (2015) Prevalence of human papillomavirus in penile malignant tumors: viral genotyping and clinical aspects. BMC urology 15: 13. pmid:25887354
- 43. Djajadiningrat RS, Horenblas S, Heideman DAM, Sanders J, de Jong J, et al. (2015) Classic and Nonclassic HLA Class I Expression in Penile Cancer and Relation to HPV Status and Clinical Outcome. Journal of Urology 193: 1245–1251. pmid:25463996
- 44. Martínez-Bailón C, Mantilla-Morales A, Méndez-Matías G, Alvarado-Cabrero I, Maldonado-Rodríguez R, et al. (2019) Human papillomavirus genotypes and P16INK4A expression in squamous penile carcinoma in Mexican patients. BMC Infect Dis 19: 1–8.
- 45. Larke NL, Thomas SL, dos Santos Silva I, Weiss HA (2011) Male circumcision and penile cancer: a systematic review and meta-analysis. Cancer Causes & Control 22: 1097–1110.
- 46. Daling JR, Madeleine MM, Johnson LG, Schwartz SM, Shera KA, et al. (2005) Penile cancer: importance of circumcision, human papillomavirus and smoking in in situ and invasive disease. International journal of cancer 116: 606–616. pmid:15825185
- 47. Bleeker M, Heideman D, Snijders P, Horenblas S, Dillner J, et al. (2009) Penile cancer: epidemiology, pathogenesis and prevention. World J Urol 27: 141–150. pmid:18607597
- 48. Ferreux E, Lont AP, Horenblas S, Gallee MP, Raaphorst FM, et al. (2003) 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. The Journal of Pathology: A Journal of the Pathological Society of Great Britain and Ireland 201: 109–118.
- 49. Soufir N, Queille S, Liboutet M, Thibaudeau O, Bachelier F, et al. (2007) Inactivation of the CDKN2A and the p53 tumour suppressor genes in external genital carcinomas and their precursors. British Journal of Dermatology 156: 448–453. pmid:17300232
- 50. Young LS, Murray PG (2003) Epstein-Barr virus and oncogenesis: from latent genes to tumours. Oncogene 22: 5108. pmid:12910248
- 51. Lam KY, Chan AC, Chan KW, Leung ML, Srivastava G (1995) Absence of Epstein–Barr virus in penile carcinoma. A study of 42 cases using in situ hybridization. Cancer 76: 658–660. pmid:8625161
- 52. Alves G, Macrini CMT, de Souza Nascimento P, Morais JC, Ornellas AA (2004) Detection and expression of Epstein-Barr Virus (EBV) DNA in tissues from penile tumors in Brazil. Cancer Lett 215: 79–82. pmid:15374635
- 53. Levi JE, Rahal P, Sarkis AS, Villa L (1998) Human papillomavirus DNA and p53 status in penile carcinomas. Int J Cancer 76: 779–783. pmid:9626339
- 54. Scheiner MA, Campos MM, Ornellas AA, Chin EW, Ornellas MH, et al. (2008) Human papillomavirus and penile cancers in Rio de Janeiro, Brazil: HPV typing and clinical features. International braz j urol 34: 467–476. pmid:18778498
- 55. Hernandez BY, Goodman MT, Unger ER, Steinau M, Powers A, et al. (2014) Human papillomavirus genotype prevalence in invasive penile cancers from a registry-based United States population. Frontiers in oncology 4.
- 56. Backes DM, Kurman RJ, Pimenta JM, Smith JS (2009) Systematic review of human papillomavirus prevalence in invasive penile cancer. Cancer Causes & Control 20: 449–457.
- 57.
SALES L (2015) Caracterização da infecção por hpv em mulheres que realizam exame de rastreio do câncer do colo do útero do município de Manaus. Amazonas: Universidade Federal do Amazonas.
- 58. Torres KL, Mariño JM, Rocha DAP, de Mello MB, de Melo Farah HH, et al. (2018) Self-sampling coupled to the detection of HPV 16 and 18 E6 protein: A promising option for detection of cervical malignancies in remote areas. PLoS One 13: e0201262. pmid:30036381
- 59. Senba M, Kumatori A, Fujita S, Jutavijittum P, Yousukh A, et al. (2006) The prevalence of human papillomavirus genotypes in penile cancers from northern Thailand. Journal of medical virology 78: 1341–1346. pmid:16927292
- 60. Young LS, Murray PG (2003) Epstein-Barr virus and oncogenesis: from latent genes to tumours. Oncogene 22: 5108–5121. pmid:12910248
- 61. Muñoz N, Bosch FX, de Sanjosé S, Herrero R, Castellsagué X, et al. (2003) Epidemiologic classification of human papillomavirus types associated with cervical cancer. N Engl J Med 2003: 518–527.
- 62. Rubin MA, Kleter B, Zhou M, Ayala G, Cubilla AL, et al. (2001) Detection and typing of human papillomavirus DNA in penile carcinoma: evidence for multiple independent pathways of penile carcinogenesis. Am J Pathol 159: 1211–1218. pmid:11583947
- 63. Dorfman S, Cavazza M, Cardozo J (2006) Penile cancer associated with so-called low-risk human papilloma virus. Report of five cases from rural Venezuela. Tropical doctor 36: 232–233. pmid:17034702
- 64. Bezerra SM, Chaux A, Ball MW, Faraj SF, Munari E, et al. (2015) Human papillomavirus infection and immunohistochemicat p16(INK4a) expression as predictors of outcome in penile squamous cell carcinomas. Human Pathology 46: 532–540. pmid:25661481
- 65. Zargar-Shoshtari K, Spiess PE, Berglund AE, Sharma P, Powsang JM, et al. (2016) Clinical significance of p53 and p16ink4a status in a contemporary North American Penile Carcinoma Cohort. Clinical genitourinary cancer 14: 346–351. pmid:26794389
- 66. Bethune G, Campbell J, Rocker A, Bell D, Rendon R, et al. (2012) Clinical and pathologic factors of prognostic significance in penile squamous cell carcinoma in a North American population. Urology 79: 1092–1097. pmid:22386252
- 67. Ferrandiz-Pulido C, Hernandez-Losa J, Masferrer E, Vivancos A, Somoza R, et al. (2015) Identification of somatic gene mutations in penile squamous cell carcinoma. Genes Chromosomes Cancer 54: 629–637. pmid:26216163
- 68. Bezerra AL, Lopes A, Santiago GH, Ribeiro KC, Latorre MR, et al. (2001) Human papillomavirus as a prognostic factor in carcinoma of the penis: analysis of 82 patients treated with amputation and bilateral lymphadenectomy. Cancer 91: 2315–2321. pmid:11413520
- 69. Micali G, Nasca MR, Innocenzi D, Schwartz RA (2006) Penile cancer. J Am Acad Dermatol 54: 369–391; quiz 391–364. pmid:16488287
- 70. Reddy JP, Pettaway CA, Levy LB, Pagliaro LC, Tamboli P, et al. (2017) Factors associated with regional recurrence after lymph node dissection for penile squamous cell carcinoma. Bju international 119: 591–597. pmid:27753187
- 71. Prue G, Lawler M, Baker P, Warnakulasuriya S (2017) Human papillomavirus (HPV): making the case for ‘Immunisation for All’. Oral diseases 23: 726–730. pmid:27492979
- 72. Michels KB, zur Hausen H (2009) HPV vaccine for all. The Lancet 374: 268–270.