Skip to main content

A comprehensive analysis of penile cancer in the region with the highest worldwide incidence reveals new insights into the disease



Although penile cancer (PC) is uncommon in developed countries, it is widespread in developing countries. The state of Maranhão (Northeast, Brazil) has the highest global incidence recorded for PC, and, despite its socioeconomic vulnerability, it has been attributed to human papillomavirus (HPV) infection. This study aimed to determine the histopathological features, the prevalence of HPV infection, and the immunohistochemical profile of PC in Maranhão.


A retrospective cohort of 200 PC cases were evaluated. HPV detection was performed using nested-PCR followed by direct sequencing for genotyping. Immunohistochemistry (IHC) was performed using monoclonal antibodies anti-p16INK4a, p53, and ki-67.


Our data revealed a delay of 17 months in diagnosis, a high rate of penile amputation (96.5%), and HPV infection (80.5%) in patients from Maranhão (Molecular detection). We demonstrated the high rate of HPV in PC also by histopathological and IHC analysis. Most patients presented koilocytosis (75.5%), which was associated with those reporting more than 10 different sexual partners during their lifetime (p = 0.001). IHC revealed frequent p16INK4a overexpression (26.0%) associated with basaloid (p < 0.001) and high-grade tumors (p = 0.008). Interestingly, p16 appears not to be a better prognostic factor in our disease-free survival analysis, as previously reported. We also demonstrated high ki-67 and p53 expression in a subset of cases, which was related to worse prognostic factors such as high-grade tumors, angiolymphatic and perineural invasion, and lymph node metastasis. We found a significant impact of high ki-67 (p = 0.002, log-rank) and p53 (p = 0.032, log-rank) expression on decreasing patients’ survival, as well as grade, pT, stage, pattern, and depth of invasion (p < 0.05, log-rank).


Our data reaffirmed the high incidence of HPV infection in PC cases from Maranhão and offer new insights into potential factors that may contribute to the high PC incidence in the region. We highlighted the possible association of HPV with worse clinical prognosis factors, differently from what was observed in other regions. Furthermore, our IHC analysis reinforces p16, ki-67, and p53 expression as important diagnosis and/or prognosis biomarkers, potentially used in the clinical setting in emerging countries such as Brazil.

Peer Review reports


Penile cancer (PC) is an uncommon neoplasm in developed countries. However, it has a high incidence in developing countries, mainly in Asia, Africa, and South America [1, 2]. When compared to those of other cancers, the epidemiological, etiological, and histopathological aspects of PC are poorly characterized.

Nevertheless, some risk factors have been proposed for PCs. The disease appears to be more prevalent in older men living in rural zones and those with low education levels [3]. PC has been linked to smoking, alcoholism, phimosis, obesity, poor hygiene of the penis, chronic inflammatory conditions, and human papillomavirus (HPV) infection [4].

According to data from the International Agency for Research on Cancer (IARC), the age-standardized rate (ASR) of PC in the United States is 0.5/100.000 men [5]. It is even lower for Jewish men born in Israel (0.1/100.000 men). The low age-standardized rate reported in Israel is associated with the cultural practice of neonatal circumcision, which appears to be a protective factor against the disease.

In contrast, the highest incidence was reported in Uganda, Africa (5.6 per 100,000 men). According to the IARC, Brazil is ranked thirteenth worldwide for incidence of PC, with a rate of 1.6/100.000 men [5]. However, the IARC only took six cities from Brazil into account to calculate the national incidence: Aracaju, Belo Horizonte, Cuiabá, Fortaleza, Goiânia, and São Paulo. Most of these cities (four out of six) are not in the poorest areas of Brazil (where PC is more common), so these data may not accurately portray the incidence of the disease in Brazil.

Since the prevalence of PC appears to vary considerably between certain regions of Brazil, local studies are essential to better understand the reality of PC in the country. These studies reveal the North and Northeast of Brazil as the regions with the highest incidence of the disease, where according to Korkes et al. 2020, between 2006 and 2018 the incidence of PC in these regions was 5.4 per 100,000 men [6].

In northeastern Brazil, the state of Maranhão has attracted the attention of the scientific community owing to its high prevalence of PC [7]. In 2018, our research group published a retrospective study of 392 patients diagnosed with PC between 2004 and 2014 in Maranhão [8] exceeds the highest national and international rates reported by IARC, revealing an alarming situation of PC in Maranhão.

Our research group and others have aimed to better understand PC features in Maranhão. Over the past five years, these studies have drawn attention to a clinical, epidemiological, and etiological profile distinct from that observed in other studies carried out in Brazil and worldwide, with marked heterogeneity of histological subtypes and high prevalence of HPV infection [7, 9,10,11,12,13].

Therefore, the present study proposes a panoramic analysis of PC in Maranhão, describing the clinical, histopathological, HPV status, and immunohistochemical (IHC) features for diagnostic and/or prognostic importance markers.


Study design

This was a cross-sectional, experimental, and retrospective study. A survey of cases diagnosed with PC between 2006 and 2020 was carried out in two reference hospitals (University Hospital of the Federal University of Maranhão (HU-UFMA) and Hospital do Câncer Aldenora Bello (HCAB)), both located in São Luís, Maranhão, Brazil. Only patients 18 years of age or older with clinical and histopathological PC diagnosis who underwent surgery, and who had tissues available for histological review were included in the study. Patients with clinical and histopathological diagnoses of PC performed under other services and with materials in inadequate conditions for analysis were excluded.

A total of 236 cases of PC (2006–2020) were identified; however, only 200 cases met the inclusion criteria. The study involved a review of clinical and histopathological findings, HPV detection and genotyping, and protein expression profiling by IHC. All clinical data described were collected from medical records. Histopathological and IHC analyzes were performed at the Laboratory of Immunofluorescence and Electronic Microscopy (LIME/HU-UFMA), and molecular assays were performed at the Clinical Research Center (CEPEC) of the HU-UFMA, São Luís, Maranhão, Brazil, and in the Nucleus of Research in Oncology (NPO), at the Federal University of Pará (UFPA), Belém, Pará, Brazil.

Both for histopathological and molecular analysis we used formalin-fixed paraffin-embedded (FFPE) tumor tissue, which was stored under appropriate conditions in the archives of the pathology service of each hospital included in the study.

Histopathological review

Specimens (n = 200) were subjected to hematoxylin and eosin (H&E) staining to confirm the histological diagnosis and classification of tumors according to the criteria of the classification of malignant tumors (TNM) in the American Joint Committee on Cancer (AJCC), 8th edition. All tumors were independently evaluated by three pathologists (SPCM, LOB, and GEBS). In cases of disagreement, the criteria were reviewed a second time.

Morphological changes suggestive of HPV infection

Histological diagnosis of HPV was performed by identifying the koilocytes. For koilocyte analysis, the three main criteria as previously described [14] were considered: (1) presence of perinuclear halo; (2) nuclear atypia, and (3) presence of multinucleation.

Tumor DNA extraction

For tumor DNA extraction from formalin-fixed paraffin-embedded (FFPE) tissue, three 10 µm cuts of each paraffin-embedded tissue were made in a histological microtome and transferred to a 2.0 ml microtube. Then, 1.0 µL of xylol was added to the tube, vortexed, and incubated in a thermoblock for 5 min at 50 °C to remove excess paraffin. The samples were centrifuged at 14,000 rpm for 2 min, and the supernatant was discarded. The xylol-washing process was repeated twice. Then 1.0 mL of absolute ethanol was added to each sample, centrifuged again at 14,000 rpm for 2 min, the supernatant was discarded, and the ethanol wash was repeated two more times. Samples were then incubated at 37ºC for 40 min to dry all alcohol in the sample. DNA isolation was done with the QIAamp DNA FFPE Tissue Kit (Qiagen, Hilden, Germany), according to the manufacturer’s instructions. Each sample was diluted in 50 µL of RNase-free TE, pH 8.0, and quantified using a NanoDropOne® spectrophotometer (Thermo Fisher Scientific, Waltham, USA), with concentrations expressed in ng/µL, and purity analysis with measures 260/280 (between 1.8 and 2.0) and 260/230 (above one). DNA samples were immediately used for molecular detection of HPV or stored at -20 °C.

HPV molecular detection

DNA extraction resulted in only 113 samples within the expected quality standards, so HPV detection was performed only in these cases. Each sample was subjected to PCR (Polymerase Chain Reaction) for the β-globin endogenous gene as second quality control. All samples were β-globin positive. HPV detection was performed by nested-PCR. In the first PCR, a set of generic primers called PGMY09/11 was used, which produced a 450 bp amplicon of the L1-HPV capsid region [15]. In the second PCR, GP5 + /6 + primers were used; which generate a 170 bp amplicon contained in the 450 bp fragment amplified in the first PCR [16]. All primer sequences and conditions for each PCR are described in the tables available in the supplementary material (S1). Amplicons were separated on a 1.5% agarose gel and visualized. Those with amplification for both PGMY09/11 and GP5 + /6 + were considered HPV-positive.

For HPV genotyping, the second PCR amplicons of HPV-positive cases were sequenced in a single direction using the GP5 + primer. Automatic sequencing was performed on the ABI PRISM 3500XL Genetic Analyzer, using the BigDye™ Terminator v3.1 Cycle Sequencing kit (Applied Biosystems™). The sequences obtained were uploaded to the MEGA 6.0 and compared with those from the GeneBank/NCBI database using the Basic Local Alignment Search Tool (BLAST) for genotyping (


IHC analysis was performed on 3 µm-thick sections from FFPE tumor tissue from 173 patients, using the EnVision FLEX Mini Kit High pH System (Agilent Technologies, Santa Clara, CA, USA). Antigen retrieval was performed using PT Link equipment (Agilent Technologies) following the manufacturer's instructions. Slides were next washed in Dako Wash Buffer, blocking of endogenous peroxidase was done for 5 min, primary monoclonal antibody for each evaluated marker was applied for 60 min, polymer application (20 min), and staining in chromogen 3,3’ diaminobenzidine (DAB) for 10 min, counterstaining with Harris hematoxylin (5 min), and analysis under an optical microscope. Table 1 describes the specifications of the monoclonal antibodies used in this study. For IHC, first, we performed immunostaining for ki-67, and 173 samples had positive internal control (lymphocytes). Cases that were negative for ki-67 internal control (n = 27) were excluded.

Table 1 Monoclonal antibodies used in immunohistochemical assays

Immunostaining was measured semi-quantitatively by two pathologists independently (GEBS and SPCM), and when in disagreement, the slides were reviewed together for consensus. For p16INK4a, the Cubilla et al. (2011) criteria were used, which proposed the following categories: (0) complete absence of p16INK4a staining in all tumor cells; (1) spotty, patchy, and discontinuous staining in some of the tumor cells; (2) more extensive but discontinuous staining pattern with small clusters of positive tumor cells; (3) cytoplasmic or nuclear staining in all tumor cells, except in hyperkeratotic or parakeratotic areas when present [17]. Only cases with pattern 3 were considered positive for p16INK4a.

For p53, cases with < 20% of tumor cells stained were considered negative, 20-50% were classified as positive ( +), and > 50% of tumor cells were positive (+ +) [18]. For ki-67, cases were classified as low (< 40.5%) and high (> 40.5%) expression, being performed globally (analyzing the entire tumoral slide extent) and in hotspots of proliferating cells [19].

Disease-free survival

For DFS analysis, only cases with clear information in the medical record (beginning and end of follow-up, and clinical outcome) were included. Only 134 cases met the inclusion criteria. We defined the time of diagnosis as the beginning of follow-up and noted the events of recurrence and/or lymph node metastasis during the analysis period. We considered recurrence when the tumor returned after the first surgical treatment and a period during which it could not be detected.

Statistical analysis

Data were tabulated using Excel (Microsoft Office 2019) and statistically analyzed using the SPSS® software v. 22 (IBM Corp., Armonk, NY, USA). The association analysis between clinical, histopathological, immunohistochemical, and HPV infection variables was performed using bivariate tests (chi-square test or Fisher's exact test). The Kaplan–Meier method was used to estimate the DFS, and the log-rank test was applied to compare the survival curves. For all analyses, results with p < 0.05 were considered significant.


Social, demographic, and clinical data

In this study, we evaluated 200 patients with clinical and pathological diagnoses of PC treated in the state of Maranhão, Brazil, between 2006 and 2020. Of these, 40 (20.0%) were admitted to the HU-UFMA and 160 (80.0%) were admitted to the HCAB. Figure 1 describes the primary socio-behavioral and clinical characteristics of these patients and their geographical distribution in Maranhão. Their ages ranged from 23 to 95 years (mean 61 + 16 years old), and there was a prevalence of men with low education (75.0% never attended school or had incomplete primary education).

Fig. 1
figure 1

Panoramic view of PC in Maranhão, Northeast Brazil. (A) Map of Brazil highlighting the state of Maranhão in green; (B) Distribution of PC cases in Maranhão according to the five main regions of the state; (C) Primary socio-behavioral and clinical characteristics of the 200 PC cases from Maranhão analyzed in this study

There was a high frequency of patients with risky sexual behavior, with 54.1% reporting more than 10 different sexual partners during their lifetime (mostly having sex without condoms), and 58.7% practicing zoophilia (sex with animals). Among those who reported a history of sexually transmitted infections (STIs), gonorrhea was the most frequent (67.6%). Only one patient tested positive for the human immunodeficiency virus (HIV). Among those who reported phimosis, 43.0% underwent circumcision in adulthood, usually preceded by surgical treatment of the tumor.

Surgery and histopathology

Surgical treatment was performed in all cases, mainly partial penis amputation (72.5%). The main affected topography was the glans (91.9%), either individually (44.4%), or simultaneously with the foreskin and/or coronal sulcus (27.3%) or penile body (20.2%). The tumor size ranged from 0.6 cm to > 10.0 cm, the majority being between 2.0 and 5.0 cm (61.9%). Tumors were unifocal in 89.9% of cases. An ulcerated aspect occurred in 30.8%, or ulcer-vegetating in 25.8%, and with an expansive pattern of growth in 75.3%.

In H&E analysis, all tumors were squamous cell carcinomas (SCC), poorly differentiated in 49.0% of cases, and mainly of the usual (39.5%) and warty (29.0%) subtypes. The presence of carcinoma in situ associated with invasive lesion was observed in 82.0% of cases, and lichen sclerosus in 27.3%. Angiolymphatic and perineural involvement was present in 35.5% and 36.5% of cases, respectively, and sarcomatoid dedifferentiation in 17.1% of cases.

We also evaluated the thickness and depth of tumor invasion. We considered “thickness” as the measure between the largest point on the tumor surface and the deepest point of the infiltration front and as “depth” the distance between the basement membrane and the deepest point of infiltration. Slight majorities had a thickness > 5.0 cm (61.0%) and depth < 5.0 cm (52.4%). According to the TNM classification, these tumors were mainly in stage pT3 (48.0%), pNx (92.5%), and pMx (100.0%). Pathological stage IIb at diagnosis was seen in 44.0% of cases. The complete characterization of the histopathological parameters is presented in Table 2.

Table 2 Histopathological aspects evaluated according to each subgroup

The main complaints reported were the appearance of a nodule/tumor in the penis (46.5%); ulcerated wound (21.0%), secretion (12.3%), nodule with presence of secretion (7.9%), dysuria and/or oliguria (7.0%), and other symptoms, such as a burning sensation, discomfort in the penile region, or itching (5.3%). These patients took an average of 17 months to seek medical help from the onset of symptoms (ranging from 1 to 96 months; two cases were excluded as they reported a delay of 180 and 240 months) and usually with locally advanced tumors (pT3-T4, 51.0%).

Additionally, we identified histological alterations suggestive of HPV infection (koilocytosis) in 78.5% of cases (Fig. 2). The presence of koilocytosis was associated with individuals who reported more than 10 lifetime sexual partners (p = 0.001), warty tumors (p < 0.0001), HPV-associated histological subtypes (p < 0.0001), and p53 negativity (p = 0.007).

Fig. 2
figure 2

Histological section of penile cancer in HE representing koilocytes (magnification 100x)

Physical examination of palpable lymph nodes was performed in 102 of the 200 cases evaluated (obtained from medical records). Of these, 61 (59.8%) had palpable lymph nodes unilaterally or bilaterally, 47 (46.1%) at admission, and 14 (13.7%) after tumor amputation. 74 cases underwent inguinal lymphadenectomy, which was performed unilaterally in 32 cases (43.2%) and bilaterally in 42 cases (56.8%). The presence of lymph node metastasis was confirmed in the histopathological analysis in 47 cases (63.5%), and 13 (27.6%) were simultaneous with primary tumor amputation/diagnosis. Extracapsular involvement was observed in 37 (78.7%) of the 74 patients with lymph node metastasis. The lymph node involvement data are presented in Table 3.

Table 3 Clinical and histopathological features of lymph node metastasis in penile cancer patients

HPV Molecular Detection

In Fig. 3 we describe the main findings of HPV detection and genotyping, showing that 80.5% (91/113) of the tumors were positive for HPV. Determination of the viral type (genotyping) was possible in only 57.1% of the positive cases (52/91). The HR-HPV genotypes detected were 16, 35, 51, 52, 53, 56, 58, 59, 66, and 73; and LR-HPV (low-risk human papillomavirus) were 6, 11, 30, 44, 70, and 74. Type 16 was the most prevalent (75% of positive cases).

Fig. 3
figure 3

Human papillomavirus (HPV) detection and genotyping data. A Global frequency of HPV in penile cancer; B Viral type of infection present in positive cases; C Separation of positive cases according to the oncogenic risk genotype detected; D Frequency of HR-HPV16 infection in positive cases; E Type of infection in HR-HPV16-positive cases

HPV positivity was associated with tumors located in the glans (p = 0.042) and the presence of carcinoma in situ associated (p = 0.029), while the absence of HPV was associated with tumors in advanced stages (pT3-pT4) (p = 0.042).

p16INK4a, p53, and ki-67 Protein Expression

Protein expression analysis of p16INK4a, p53, and ki-67 by IHC was performed in 173 cases. About 26.0% of cases (45/173) overexpressed p16INK4a and 74.0% (128/173) were negative. The p16INK4a overexpression (positive, score 3) were associated with basaloid (p < 0.0001), poorly differentiated (p = 0.008), and high proliferative tumors (high ki-67, p < 0.0001). Furthermore, the absence of p16INK4a overexpression was associated with tumors without angiolymphatic invasion (p = 0.031). The full description of the statistical analyses for koilocytosis, molecular detection of HPV, and p16INK4a expression is described in the supplementary material (S2).

For p53, 55.0% of the cases were negative (95/173), and 45.0% were positive. Among the positive cases, 56.4% were ( +) (44/78) and 43.6% were (+ +) (34/78). The absence of p53 expression was associated with warty tumors (p = 0.001), presence of carcinoma in situ (p = 0.0.41), HPV-associated histology (p < 0.0001), absence of perineural invasion (p = 0.036), koilocytosis (p = 0.002), and tumors with an expansive growth pattern (p = 0.002). In contrast, p53 positivity was associated with poorly differentiated tumors (p < 0.0001), pT3-pT4 tumors (p = 0.031), stage II disease (p = 0.020), and lymph node metastasis (p = 0.046).

For ki-67, 26.6% of the cases had high expression (46/173), and 73.4% had low expression (127/173). When analyzing only the hotspot zones, 57.8% of tumors showed high expression (100/173), and 42.2% exhibited low ki-67 expression (73/173). The high expression of ki-67, either global or in hotspot, was associated with basaloid tumors (p < 0.0001), poorly differentiated tumors (p < 0.0001), presence of angiolymphatic (global p = 0.001; hotspot p = 0.011), and perineural invasion (global p = 0.049; hotspot p < 0.0001), pT3-pT4 tumors (global p = 0.031; hotspot p = 0.014), stage III-IV disease (global p = 0.038; hotspot p = 0.002), and lymph node metastasis (global p = 0.004; hotspot p < 0.0001).

In contrast, low ki-67 expression was associated with warty tumors (p = 0.017), with an absence of sarcomatoid transformation (global p = 0.006; hotspot p = 0.001), an expansive invasion pattern (global p = 0.001; hotspot p < 0.0001), small tumors (global p = 0.034), negativity for p53 (global p = 0.013; hotspot p = 0.002), and with no p16INK4a overexpression (global p < 0.0001; hotspot p < 0.0001).

Figure 4 shows the immunostaining profiles for each marker analyzed (p16INK4a, p53, and ki-67, and a full description of the statistical analyses for p53 and ki-67 is described in the supplementary material (S3).

Fig. 4
figure 4

Immunostaining patterns of p16INK4a, p53, and ki-67. A p16INK4a positive control; B p16INK4a negative control; C penile cancer showing score 1 for p16INK4a; D Score 2 for p16INK4a; E Score 3 for p16INK4a; F p53 positive control; G p53 negative control; H PC showing no p53 expression; I 40% of tumor cells immunostained for p53; J 90% of tumor cells immunostained for p53; (K) tonsil used as ki-67 positive control; L tonsil with no ki-67 antibody applied and used as the negative control; M 20% of cells immunostained for ki-67; N 40% of cells immunostained for ki-67; O 90% of cells immunostained for ki-67. Magnification 100x (A-O)

Disease-free survival analysis

For the DFS analysis, we obtained survival data for 134 (67.0%) of the 200 cases evaluated. Approximately 66 patients (33.0%) were excluded because they did not have enough information in their medical records for follow-up analysis or were lost to follow-up. The follow-up time ranged from 1.0 to 114 months (mean 20.96 + 21.32). By the end, 53 cases (53/134) had local recurrence and/or lymph node metastasis (39.5%), and 81 cases (81/134) had no evidence of the disease after surgical treatment (60.5%).

We performed the Kaplan–Meier test to estimate the DFS and the log-rank test to compare the survival curves according to clinical, histopathological, HPV infection parameters (koilocytosis, molecular detection, and p16INK4a), and p53 and ki-67 protein expression. In Fig. 5, we describe the impact of a delayed diagnosis (log-rank p = 0.056) and age on the DFS of these patients (log-rank p = 0.106), however, we did not observe a statistically significant difference. In Fig. 6 (A-F), we show the DFS curves according to (A) degree of differentiation (grade), (B) primary tumor (pT), (C) stage, (D) sarcomatoid transformation, (E) pattern of invasion, and (F) depth of invasion. In Fig. 7, we describe PCs with angiolymphatic and perineural invasion and their respective DFS. We observed a statistically significant difference between the survival curves for all parameters shown in Figs. 6 and 7. In Fig. 8 (A-F), we describe DFS according to the presence of koilocytes, molecular detection of HPV, and protein expression of p16INK4a, p53, and ki67 (global and hotspot), where we observed a statistically significant difference only with p53 (log-rank p = 0.032) and ki-67 expression (log-rank global p = 0.002; hotspot p < 0.0001).

Fig. 5
figure 5

Association of disease-free survival with the delay between the onset of symptoms and diagnosis by patients with PC in Maranhão, Brazil

Fig. 6
figure 6

Disease-free survival according to histological features. A grade; B primary tumor (pT); C stage; D presence of sarcomatoid transformation; E pattern of invasion and F depth of invasion

Fig. 7
figure 7

Association of disease-free survival (DFS) with the presence of angiolymphatic and perineural invasion. A left, tumor vascular involvement; right, its impact on the DFS of patients with PC; B left, tumor perineural involvement; right, its impact on the DFS of patients with PC. Magnification 400x (A-B)

Fig. 8
figure 8

Association of disease-free survival with A the presence of koilocytes; B HPV molecular detection; C expression of p16INK4a; D p53 expression; E global expression of ki-67; and F expression of ki-67 in hotspot


Recurrent Profile of PC in Maranhão, Brazil

PC is a neglected and aggressive disease with a sizable impact on quality of life [20]. Its incidence is high in developing countries such as Brazil, particularly in the North and Northeast regions, representing a serious local public health problem [1, 8]. In 2018, our research group found that the state of Maranhão in Northeast Brazil has the highest incidence of PC ever recorded in the world [8]. In 2020, we characterized the sociodemographic features of patients with PC in the region [7].

In this study, we evaluated 200 PC cases diagnosed in Maranhão, integrating clinical-histopathological, immunohistochemical, and HPV infection for a better understanding of disease etiology. The cases were mostly composed of men in their 60 s, the age at which the risk of developing PC is higher; however, it still frequently occurs in younger men [1, 9, 21, 22].

Men under 60 years old accounted for 42.3% of all cases evaluated, with the youngest being only 23 years old. In two studies previously published by our group, 19.7% and 22% of the men were less than 40 years old, respectively [7, 8]. In the present study, we observed a decrease in the survival time of younger individuals (18–40 years old) in relation to those over 40 years old. According to the literature, younger men seem to develop mostly HPV-associated PC, being more prone to relapse, infiltrative patterns, and perineural invasion [23, 24].

We note the prevalence of individuals from the North of Maranhão, where the reference centers are located (in the capital, São Luís). The southern region seems to account for the lowest percentage of cases and this may be associated with the evasion of these patients to cities in neighboring states, such as Piauí (closer than to the capital), or to the southern region of Brazil (looking for better therapeutics approaches). The large migration of PC patients from poor regions to underdeveloped areas is commonly reported in Brazil [1], suggesting that the incidence of PC in Maranhão is possibly even higher than that described by Coelho et al. (2018) [8].

Socioeconomic vulnerability seems to be one of the main risk factors for PC [25]. In the present study, the high prevalence of the disease in men with low education, those with a family income less than a Brazilian minimum monthly wage, and residents/workers from rural areas, reinforce this profile. However, it is important to remember that knowledge about PC and its risk factors seems to be scarce, even among highly educated men [26].

Exploring the profile of these patients, we observed that most of them reported poor genital hygiene habits, an important and known risk factor for PC [27]. Proper penile hygiene has been one of the main prophylactic measures for PC, and when not performed, it favors the emergence of a favorable microenvironment for the proliferation of uropathogenic microorganisms, such as HPV [28]. In addition, it increases the chances of the onset of episodes of chronic inflammation (balanitis and lichen sclerosus), which may also contribute to PC development [29,30,31].

In addition to the level of education, phimosis also seems to contribute to inadequate cleaning of the penis due to the difficulty or inability in exposing the glans. This condition was described in 70.5% of the cases evaluated in our study and is frequently reported in PC patients. In general, the prevalence of phimosis in PC is approximately 25% to 75%, [32] while it occurs only in 3.4% of healthy adult men [33]. Some studies suggest that phimosis increases the risk of developing PC by 3.5 times, with no protective effect if dealt with after the neonatal stage [3, 34, 35].

Furthermore, we observed a high frequency of risky sexual behavior (54.1%), which can significantly contribute to the acquisition of sexually transmitted infections, such as HPV [36, 37]. A high prevalence (58.7%) of patients who had sex with animals (zoophilia) was shocking data but had been reported before. Vieira et al. (2020) described the practice of zoophilia in 60% of the cases in Maranhão [7]. Both our data and those published by Vieira et al. (2020) corroborate the findings of Zequi et al. (2012), who demonstrated via univariate and multivariate analyses that zoophilia is an important risk factor for PC, but does not seem to impact pathological and clinical prognoses [38].

Upon evaluating the primary symptoms of these patients, we found that the medical records mainly described the appearance of a nodule/tumor in the penis (46.5%) or a foul-smelling ulcerative wound (21.0%). These symptoms are consistent with the main manifestations of the disease, which include changes in skin thickness or color, a wound not having healed in four weeks, rash, bleeding from the penis or under the foreskin, and a smelly odor [39]. Vieira et al. (2020) described the presence of pruritus as the primary initial symptom in patients with PC in Maranhão (26.6%) and discussed its association with lichen sclerosus, despite only 6.0% of cases having this condition in their study [7]. Interestingly, in our study, 27.3% of the cases had associated lichen sclerosus, and the presence of pruritus was observed in 20.2% of the cases, corroborating the findings of Nasca et al. (1999) [40] and Vieira et al. (2020) [7].

As the penis is an external organ, the PC symptoms are classic and easily recognizable. However, we noted an average delay of 17 months between the onset of symptoms and seeking medical care, similar to what was previously described in Maranhão by Vieira et al. (2020), who described a delay of 18.9 months [7]. Late diagnosis is generally attributed to a lack of knowledge about the symptoms of the disease, associated stigmas, and fear of seeking medical help [41]. These factors are intensified in Maranhão, where a large part of the population is socioeconomically vulnerable [42]. The region has the second lowest Human Development Index (HDI) in Brazil (0.639), has rural socioeconomic characteristics, and has historically been marked by significant social inequality and extreme poverty [42].

In this context, access to basic health, especially in rural areas, is scarce, making it difficult to carry out follow-up and early diagnosis of several diseases, including PC. All these factors contribute to an advanced stage at diagnosis, a high rate of metastasis, limitations in the therapeutic approach, and a low overall survival rate. Population-based analyses from Europe and the United States have shown that, unlike other types of tumors, tumor-specific survival rates for PC have not shown any improvement since 1990 [43, 44].

Our findings reinforce a recurrent social-economic profile of patients with PC [7, 8, 45, 46] and call attention to an alarming situation of this disease in poor regions, especially in Maranhão [7, 8, 47]. We suggest that the sum of factors that simultaneously favor disease onsets, such as late diagnosis and advanced stages, may contribute significantly to the high incidence of PC in Maranhão, Brazil, and public policies aimed at remedying these factors in the region are urgent.

High frequency of HPV-associated PC in Maranhão: differences in disease behavior

The high prevalence of HPV-positive tumors reported in the last five years has been one of the most notable findings in Maranhão. The first study on the detection and genotyping of HPV in local cases was conducted by De Sousa et al. (2015) [48] and Martins et al. (2018) [9], where HPV was reported in 63.15% and 89.1% of tumors analyzed, respectively. Recently, another study described the presence of HPV in 96.4% of cases, the highest frequency ever recorded for PC [13].

All these studies carried out in the region showed a frequency of HPV-positive tumors above the global average of approximately 46.9—50.8% [49,50,51], and describe HR-HPV 16 as the most common viral type associated with PC. To better understand this situation, in the present study, we used morphological, immunohistochemical, and molecular criteria to investigate the involvement of HPV in the etiology of PC. Based on this, when we performed the molecular detection of HPV, we found a global frequency of HPV infection in 80.5% of PCs, mainly by HR-HPV (84.5%), especially type 16 (75%). We found an association between global HPV positivity and tumors located in the glans, those with carcinoma in situ, and pT1. In fact, glans has been reported as a frequent site of penile HPV infection (35.8%) and is more likely to be infected with multiple HPV types than other sites [52].

In DFS analysis, we did not observe any difference in survival curves according to the HPV status. In contrast, in a cohort from Maranhão, Martins et al. (2018) found better DFS in HPV-positive PC patients than in HPV-negative PC patients [9]; this result has already been reported for PC and other HPV-associated tumors, such as the head and neck [53,54,55]. Glegoire et al. (1995) demonstrated an association between HPV positivity and worse prognostic factors, such as high staging and aggressive growth patterns [56]; Wiener et al. (1992) found that the overall survival and DFS did not differ between groups in Kaplan–Meier analysis [57].

Histopathological analysis revealed morphological alterations suggestive of HPV infection (koilocytosis) in 78.5% of the evaluated cases. This corroborates previous findings in Maranhão, where koilocytosis was observed in 74.5% of the PC analyzed [9]. The presence of koilocytosis was associated with men reporting more than 10 lifetime sexual partners, tumors with HPV-associated histology, and warty subtype. These findings illuminate the pathophysiology of HPV-associated PC and that in men with risky sexual behavior [58].

Koifman et al. (2011) observed koilocytosis in 88.4% of PC cases evaluated, with high specificity (90%), but low sensitivity in detecting the virus [59]. De-Paula et al. (2007) observed the presence of koilocytosis in 63.1% of the PC samples evaluated, being associated with tumors with low or moderate staging and better DFS [60]. In Maranhão, Nascimento et al. (2020) described angiolymphatic invasion and absence of koilocytosis as predictive factors for lymph node metastasis [10], corroborating previous findings by Ornellas et al. (2008) [61]. However, in the present study, we found no difference in DFS between the groups with and without koilocytosis, as well as no association between koilocytosis and lymph node metastasis.

The presence of koilocytosis is well accepted as pathognomonic evidence for possible the presence of HPV, and has been widely used in screening for cases of viral infection [62]. However, koilocytes are not always present; therefore, it is necessary to use more sensitive and specific methods that also provide information about the viral type and oncogenic risk [63].

To better understand aspects related to HPV infection, we searched the literature for data on the global prevalence of the virus in the general population. In Brazil, the prevalence of HPV in individuals without clinical symptoms appears to vary from 30 to 75%, depending on the region studied [64]. Interestingly, Maranhão represents one of the highest rates recorded in the country (59.1%), according to sampling carried out in the capital (São Luís) [65]. It is believed that in the interior regions of the state, where most PC patients come from, these rates may be even higher, but there are no studies with a representative sample of these locations.

The high prevalence of healthy individuals infected with HPV in Maranhão has been attributed to factors such as early initiation of sexual activity (average age, 15.5) and the high number of sexual partners [65]. In the present study, we did not assess the onset of sexual activity, but we observed a statistically significant association between risky sexual behavior and the presence of koilocytosis in PC patients, which corroborates the epidemiological HPV aspects of this disease.

Knowledge of the regional prevalence of HPV infection is important for understanding the local etiology of PC, as the incidence of HPV-associated tumors seems to be proportional to the prevalence of viral infection in healthy individuals [66]. In our series, when we searched for tumors that were negative for koilocytosis, HPV infection, and p16INK4a overexpression, only one case was not positive for any of these three parameters. Interestingly, this was the first case report of pseudoangiosarcomatous squamous cell carcinoma on the penis, an uncommon and highly aggressive variant of squamous cell carcinoma [67].

In addition to koilocyte and HPV molecular detection, we performed p16INK4a protein expression analysis, which is widely used as an indirect marker of HR-HPV infection and has been associated with basaloid and high-grade tumors [9, 68, 69]. However, the absence of p16INK4a overexpression does not exclude infection by LR-HPV or even the absence of HPV infection [17, 70].

We found p16INK4a overexpression in 26.0% of PC cases. In the present study, all basaloid tumors showed p16INK4a overexpression, reinforcing its secondary use as a tool in the differential diagnosis of this histological subtype [9, 17, 70]. Furthermore, the lack of p16INK4a overexpression was associated with the absence of angiolymphatic invasion.

According to Cubilla et al. (2011), p16INK4a overexpression is closely correlated with the presence of the HR-HPV genotype in PC, as observed in cervical carcinoma and other HPV-related tumors [17, 70]. However, in the present study, no association between p16INK4a overexpression and HR-HPV infection was observed. We suggest that limitations in HPV detection methods may explain these findings, such as: (1) we used DNA samples obtained from FFPE tissue, unlike Martins et al. (2018), who used fresh tumor samples; and (2) viral genotyping was not done for all positive cases.

Despite the high frequency of cases with koilocytes and HPV-positivity, the frequency of p16INK4a overexpression was low compared to that described by Martins et al. (2018) in Maranhão, where 40% of the cases had overexpression of p16INK4a [9]. We speculate that the use of different antibody clones for IHC of p16INK4a might underlie this difference. Martins et al. (2018) described four basaloid tumors in a total of 55 cases, while here we had eight basaloid tumors in approximately 173 cases.

We also did not observe a statistically significant difference in DFS with the level of p16INK4a expression. However, patients with p16INK4aoverexpression lived on average less (30.1 months) when compared to patients without p16INK4a overexpression (46.8 months). These findings do not corroborate with previous studies where PC men with HPV or p16INK4a overexpression generally have a significantly more favorable DFS compared to those without HPV or p16INK4a overexpression [71]. Furthermore, a recent meta-analysis showed that p16INK4a overexpression was significantly associated with and was an independent factor for better DFS in PC patients (HR = 0.45, 95% CI: 0.30–0.67, P < 0.001) [72].

Given these controversial data, we performed DFS analysis for p16INK4a overexpression according to the following criteria: (1) basaloid tumors were excluded; and only (2) HR-HPV positive cases; (3) usual subtype tumors; (4) warty subtype tumors; and (5) stage I-II and III-IV tumors. However, in all analyses, there was no significant difference in DFS, and the mean survival time of the p16INK4a overexpressing group remained lower than that of no overexpressed p16INK4a.

Unfortunately, Martins et al. (2018) did not describe survival data based on p16INK4a expression, making it difficult to compare their findings in Maranhão with those observed in the present study [9]. However, similar to what was observed in our study, in a study conducted in Northern Brazil (state of Amazonas), no differences were observed between the survival curves of PC patients with p16INK4a overexpression (log-rank p = 0.753) or molecular HPV positivity (log-rank p = 0.979) [73]. According to Ornellas & Ornellas (2018), it is uncertain whether PCs involving HPV infection have better survival profiles than those without HPV infection [74]. Our data and those published in Amazonas suggest a different impact of p16INK4a overexpression in PC patients from these Brazilian regions.

For histopathological classification, all the evaluated tumors were classified as squamous cell carcinoma (SCC), the most common type of PC [23]. According to the WHO [75], these tumors are grouped into HPV-associated subtypes (basaloid, papillary-basaloid, warty, warty-basaloid, and clear cell carcinoma) and non-associated HPV (usual, pseudohyperplastic, pseudoglandular, verrucous pure or cuniculatum, papillary, adenosquamous, and sarcomatoid).

In the present study, most cases were the usual subtype (39.5%), followed by the warty (29.8%). These frequencies differ from other studies in the area, such as that carried out by Cubilla et al. (2001), where out of 61 tumors evaluated, 59% corresponded to the usual subtype and only 10% corresponded to warty tumors [76]. In 2010, Chaux et al. compared the distribution of different histological subtypes of PC from regions of high (Paraguay, 144 cases) and low (USA, 157 cases) incidence and observed that the prevalence of histological subtypes in the Paraguayan and American series was as follow: usual, 49.3% and 46.5%; verrucous, 8.3% and 7.6%; papillary NOS, 7.6% and 5.7%; warty, 6.9% and 8.3%; basaloid, 4.2% and 7.0%; sarcomatoid, 0.7% and 0.6%; adenosquamous, 3.5% and 0.6%; and mixed, 19.4% and 23.6%, respectively [77].

In 2018, Martins et al. observed that 47.3% of the PC cases evaluated in Maranhão were usual and 29.1% were warty subtypes [9]. Subsequently, Vieira et al. observed 40.0% of usual tumors and 33.0% of warty tumors [7]. Another important aspect is the high frequency of mixed tumors (with characteristics of more than one histological subtype), which in the present study was approximately 32.0%, higher than that described by Martins et al. 2018 (16.3%) and Vieira et al. 2020 (25.0%). Considering mixed tumors according to the presence or absence of an HPV-associated component, similar to the classification described by Canete-Portillo et al. (2020) [78], all mixed tumors in the present study had a basaloid, warty-basaloid, or warty component associated with another component (usual). Thus, the overall prevalence of HPV-associated histology in the cases evaluated in our study go up to 77.5%.

The high frequency of warty tumors in the present study and in others carried out in Maranhão, as well as the prevalence of mixed tumors with HPV subtypes associated, differs drastically from that observed in other studies, and supports the high prevalence of HPV-associated PC in the region [76, 77].

Expression of p53 is Related to Worse Prognosis and Lymph node metastasis

When we evaluated the p53 expression profile, negative cases were more prevalent (55.0%) than positive cases (45.0%). The statistical analysis revealed that p53 positivity was associated with worse prognostic factors such as high-grade tumors (G3), pT3-pT4, those in stage II, and those with lymph node metastasis, corroborating previous studies with PC and other tumor types [79,80,81].

Previously, all poorly differentiated PCs analyzed by Marinescu et al. (2016) were p53 positive [82], and Lopes et al. (2002) found that p53 expression was an independent factor for lymph node metastasis (relative risk 4.8, 95% CI 1.6 to 14.9); the 5 and 10-year overall survival rates were 65.2% and 54.6% in PCs with negative p53, and 30.2% and 26.4%, respectively, in those with positive p53 [18]. In this study, p53 expression was associated with lower DFS when compared to negative cases.

According to the DGIdb database, the TP53 gene is the target of 193 therapeutic drugs [83], many of which are already used for head and neck tumors, which appear to have a mutational signature very similar to PC [84]. Overexpression of p53 is frequently associated with disruptive mutations and oncogenic activation in HPV-negative PC, possibly in cooperation with other mutations in genes of the HER-PTEN-AKT pathway [85].

In a study conducted by Jacob et al. (2019), most patients with metastatic PC had mutations in TP53 and were HPV-negative [86] This suggests that TP53 mutations appear to be the main factor associated with metastasis in patients with advanced PC, whereas HPV-related oncogenic mechanisms appear to be more effective in the early stages of the disease and are gradually lost [85, 86]. These findings may explain the pathophysiological differences observed in this study and in other studies carried out in Maranhão, where most cases are diagnosed in the advanced stages of the disease.

In contrast to what we observed for p53 positivity, the absence of p53 expression was associated with warty tumors, presence of carcinoma in situ, tumors with HPV-associated histology, absence of perineural invasion, koilocytosis, and expansive growth pattern. These findings corroborate the main pathway for HPV action, where the viral oncoprotein E6 binds to the p53 protein and inhibits its tumor suppressor function [87]. In addition, a previous study in Maranhão found that 81.0% of PC cases were negative for p53 expression by IHC and 83.3% presented TP53 gene downregulation [13].

Ki-67 Expression as a Marker of Tumor Aggressiveness in PC, Including Perineural and Angiolymphatic Invasion

Most tumors analyzed in the present study had a low global expression of ki-67 (73.4%), but high expression in hotspot zones (57.8%). Here, we observed an association between high ki-67 expression and basaloid and high-grade tumors, presence of angiolymphatic and perineural invasion, pT3-pT4 tumors, stage III-IV, and lymph node metastasis. The association between high expression of ki-67 and worse prognostic factors has been previously described in PC [19, 88, 89].

Ki-67 analysis has been used as a good parameter to assess the proliferative index and aggressiveness of numerous tumors [90, 91]. In PC, Stankiewicz et al. (2012) found an association between high expression of ki-67, basaloid tumors, and HPV infection, corroborating the results of the present study, in which all basaloid tumors had high ki-67 expression [89]. In contrast, we observed that warty tumors were associated with low ki-67 expression (global and hotspot). This difference corroborates the clinical behavior of each histological subtype, where basaloid tumors generally have fast growth and generally poor prognosis, while warty tumors have slow growth and are more often a good prognosis [88].

A strong association between high ki-67 expression and poorly differentiated tumors has been described previously in PC [89]. Berdjis et al. (2005) observed an inverse relationship between ki-67 expression and tumor differentiation (p < 0.0005), but they did not observe a relationship with advanced staging or lymph node metastasis, unlike our findings [19].

The association between high ki-67 expression and angiolymphatic and perineural invasion is another important finding of our study and was previously described for PC. Angiolymphatic and perineural invasions are known to be associated with a worse prognosis, and its detection is sometimes influenced by the observer (pathologist) bias, suggesting that IHC for ki-67 can refine the diagnosis and prompt one to search for neural or vascular invasion in the histopathological analysis [92, 93].

Protzel et al. (2007) and Zhu et al. (2007) described a strong association between high ki-67 expression and lymph node metastasis in PC [79, 88], and Warli et al. (2020) argued that this association appears to be independent of grade and stage [94]. Our findings also demonstrated an association between high ki-67 expression and the occurrence of lymph node metastasis in PC patients, indicating the role of ki-67 as a predictive factor for metastasis in these patients.

Despite the association with worse prognostic factors and lymph node metastasis, we did not observe a statistically significant difference between the DFS rate according to the ki-67 expression profile, similar to the study published by Stankiewicz et al. (2012) [89]. We hypothesized that the cumulative effect of other factors such as staging, angiolymphatic invasion, and perineural invasion in our study may have influenced DFS patients more than ki-67 high expression. For our studies, future multivariate analyses are required to elucidate these aspects. However, high ki-67 expression has been reported to be associated with worse DFS in PC (log-rank p = 0.0098) [88].

Considerations and perspectives

PC is a neglected, heterogeneous, aggressive, and deadly disease that is typically not diagnosed early. There is a high incidence of PC in the Maranhão state, and little is known about the local etiological profile of these tumors. Most patients are socioeconomically vulnerable and come from the state's countryside. There are few centers for treating the disease in the region, most of which are located in the capital, São Luís, making early diagnosis and monitoring of these patients challenging.

Here, we discussed the impact of the geographic isolation of the state capital on the displacement of these patients, suggesting that the actual incidence of the disease is more significant than that described previously [8]. Our findings reinforce a profile of the disease in the state consistent with other studies and bring new insights into possible factors that, together, contribute to Maranhão having the highest incidence of PC in the world, such as the long delay in diagnosis and the high prevalence of HPV infection.

According to Ornellas & Ornellas (2018), when PC is considered globally, only a tiny proportion is described as HPV-associated. Therefore, increased men's education and prevention strategies such as condom use, hygienic measures, smoking cessation, and avoidance of chronic inflammatory states can considerably impact the pathogenesis of precancerous lesions of the penis [74]. In Maranhão, where a significant number of PCs appear to be HPV-associated and most patients are socioeconomically vulnerable, both HPV immunization and patient education could impact disease rates.

Our data discuss HPV infection in a multi-methodological way, integrating histological, molecular, and immunohistochemical aspects and bringing new insights into the role of HPV on histological presentation and clinical behavior of these tumors. An interesting question involves the histological and prognostic variation observed in HPV-associated tumors. For example, why do warty and basaloid tumors, both HPV-associated, have widely different prognoses? Why does p16INK4a overexpression appear to be related to decreased disease-free survival and what is its impact on overall and cancer-specific survival in PC patients? To answer these questions we need to explore the correlation between CDKN2a (p16INK4a) copy number, mRNA expression, methylation profile, and protein expression (IHC).

Another important approach that we are currently investigating is the biological profile of viral oncoproteins (E6 and E7), both by gene (mRNA) and protein expression. These data linked to the characterization of p16INK4a can reveal the real status of HPV infection in these patients, whether it is active and impacting tumor progression or part of the HPV-positive cases are just a background due to the high rate of HPV infection in the general population, especially in Maranhão (as discussed).

The joint analysis of this information can help to better characterize each of the cited markers and their real role in the patient's prognosis. The perspective of our group has been, given the vast majority of HPV-positive cases, to investigate the differences within this group. We are currently performing all these analyzes on a prospective cohort of 60 cases from Maranhão, involving histological (H&E, IHC), genomic (exome, RNA-seq, miRNA-seq), and metagenomic tools to unravel these questions. All these approaches ongoing rely on the expertise of geneticists, physicians, pathologists, and oncologists with experience in PC to perform an integrative analysis between genomic (multi-omics) and the clinical aspects of PC. We hope with this approach to also elucidate the role of changes in p53 in PC carcinogenesis, evaluating somatic mutations and gene expression patterns.

At the time, the associations found between p16INK4a, p53, and ki-67 expression have clinical relevance and seem to impact the DFS of these patients. The test used to assess these protein expression profiles (IHC) is offered free of charge by the Brazilian Unified Health System (SUS) and can be used in histopathological routines to help diagnostic and prognostic classification of these patients.

Despite the difficulties inherent in the PC study, the present study strongly points to a clinical-histopathological and etiological profile different from that observed in other studies in Brazil and worldwide. It highlights the importance of local and collaborative research aimed at understanding these tumors' etiological aspects, behavior, and progression, especially in places with high incidences, such as Brazil. We emphasize that, due to limitations in current protocols for the management of PC patients, early diagnosis becomes even more decisive for better therapeutic conduct. This reality is expected to change since recently, the Brazilian consensus statement for managing PC patients in low- and middle-income countries has been proposed [95].


In the present study, we presented important clinical-histopathological and pathophysiological features of PC. However, some limitations may reduce the clinical impact of the data obtained and must be considered. First, the retrospective aspect of the study makes it difficult to have reliable control over the variables evaluated, creating gaps in the data. Second, the use of DNA from FFPE tumor samples may increase HPV-PCR false negatives and the absence of genotyping data for all HPV-positive tumors limits the stratification of cases according to the HPV genotype. Third, this study had many censored patients in DFS analysis due to loss of follow-up and there is no data regarding cancer-specific and/or overall survival. Despite this, the data obtained from a large Brazilian cohort of PC present a wealth of information that significantly contributes to a better understanding of this neglected disease. All limitations may be further explored in a future prospective study, as discussed.


Our study obtained important clinical and histopathological data for the characterization of the epidemiological profile of PC in Maranhão, Brazil. Our IHC findings revealed important parameters with prognostic values, including ki-67 and p53 expression associated with lymph node metastasis, grade, pT, stage, pattern, and depth of invasion, as well as ki-67 and p53 expression associated with decreased survival of PC patients. Additionally, p16INK4a overexpression can be used as a marker for basaloid tumors, an HR-HPV-related histological subtype. The IHC is easily accessible in poorer regions such as Brazil, reinforcing its applicability in clinical practice.

Our data reaffirmed the high incidence of human papillomavirus (HPV) infection in PC and highlighted its possible association with worse clinical prognosis factors in the cohort evaluated, differently from what was observed in other regions. The association between HPV histological changes and men with risky sexual behavior draws attention to the urgency of effective public policies to vaccinate men to prevent the disease. This is especially important in Brazil, where there is a high rate of amputations as the first treatment option. Finally, our integrative approach offers new insights into possible factors that together contribute to Maranhão having the highest global incidence of PC ever recorded and highlights the importance of studies that aim to understand the pathophysiological behavior of HPV-related PCs better.

Availability of data and materials

The data used in this study are not publicly available due to privacy and ethical regulations, being made available only when requested to the corresponding author.



Penile cancer


Polymerase chain reaction


Formalin-fixed paraffin embedded


Human papillomavirus


High-risk human papillomavirus


Human immunodeficiency virus


Hematoxilina e eosina


Hospital Universitário da Universidade Federal do Maranhão


Hospital do Câncer Aldenora Bello


Disease-free survival




International Agency for Research on Cancer


Age-standardized rate


Low-risk human papillomavirus


  1. Favorito LA, Nardi AC, Ronalsa M, Zequi SC, Sampio FJB, Glina S. Epidemiologic study on penile cancer in Brazil. Int Braz J Urol. 2008;34:587–91.

    Article  PubMed  Google Scholar 

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

    Article  Google Scholar 

  3. Romero FR, Romero AW, De Almeida RMS, de Oliveira FC, Filho RT. Prevalence and risk factors for penile lesions/anomalies in a cohort of Brazilian men ≥ 40 years of age. Int Braz J Urol. 2013;39:55–62.

    Article  PubMed  Google Scholar 

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

    Article  PubMed  Google Scholar 

  5. Forman, D, Bray, F, Brewster, DH, Gombe Mbalawa, C, Kohler, B, Piñeros, M. Cancer Incidence in Five Continents, Volume X. IARC Scientific Publication No. 164. Lyon, France: International Agency for Research on Cancer (2014).

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

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

    Article  CAS  Google Scholar 

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

  9. Martins VA, Pinho JD, Teixeira Júnior AAL, Nogueira LR, Silva FF, Maulen VE, et al. P16INK4a expression in patients with penile cancer. PLoS ONE. 2018;13(10): e0205350.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Nascimento ADMTD, Pinho JD, Teixeira Júnior AAL AAL, Larges JS, Soares FM, Calixto JRR, et al. Angiolymphatic invasion and absence of koilocytosis predict lymph node metastasis in penile cancer patients and might justify prophylactic lymphadenectomy. Medicine (Baltimore). 2020;99(9):e19128.

    Article  Google Scholar 

  11. Pinho JD, Silva GEB, Júnior AALT, de Castro Belfort MR, Macedo JM, da Cunha IW, et al. MIR-107, MIR-223-3P and MIR-21-5P Reveals Potential Biomarkers in Penile Cancer. Asian Pacific J Cancer Prev. 2020;21:391–7.

    Article  CAS  Google Scholar 

  12. Vieira CB, Teixeira-Júnior AAL, Feitoza L, Pinho JD, Calixto JRR, Do Nascimento FSMS, et al. A cohort study among 402 patients with penile cancer in Maranhão, Northeast Brazil with the highest worldwide incidence. BMC Res Notes. 2020;13:10–2.

    Article  Google Scholar 

  13. Macedo J, Silva E, Nogueira L, Coelho R, da Silva J, dos Santos A, et al. Genomic profiling reveals the pivotal role of hrHPV driving copy number and gene expression alterations, including mRNA downregulation of TP53 and RB1 in penile cancer. Mol Carcinog. 2020;59:604–17.

    Article  CAS  PubMed  Google Scholar 

  14. Abadi MA, Ho GY, Burk RD, Romney SL, Kadish AS. Stringent criteria for histological diagnosis of koilocytosis fail to eliminate overdiagnosis of human papillomavirus infection and cervical intraepithelial neoplasia grade 1. Hum Pathol. 1998;29(1):54–9.

    Article  CAS  PubMed  Google Scholar 

  15. Coutlée F, Gravitt P, Kornegay J, Hankins C, Richardson H, Lapointe N, et al. Use of PGMY primers in L1 consensus PCR improves detection of human papillomavirus DNA in genital samples. J Clin Microbiol. 2002;40(3):902–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. Jacobs MV, Snijders PJF, Van den Brule AJC, Helmerhorst TJM, Meijer CJLM, Walboomers JMM. A general primer GP5+/GP6+-mediated PCR-enzyme immunoassay method for rapid detection of 14 high-risk and 6 low-risk human papillomavirus genotypes in cervical scrapings. J Clin Microbiol. 1997;35:791–5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

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

  18. Lopes A, Bezerra ALR, Antonio C, Pinto L, Serrano SV, De Mello CA, et al. p53 as a new prognostic factor for lymph node metastasis in penile carcinoma: Analysis of 82 patients treated with amputation and bilateral lymphadenectomy. J Urol. 2002;168:81–6.

    Article  PubMed  Google Scholar 

  19. Berdjis N, Meye A, Nippgen J, Dittert D, Hakenberg O, Baretton GB, et al. Expression of Ki-67 in squamous cell carcinoma of the penis. BJU Int. 2005;96:146–8.

    Article  CAS  PubMed  Google Scholar 

  20. Torbrand C, Wigertz A, Drevin L, Folkvaljon Y, Lambe M, Håkansson U, et al. Socioeconomic factors and penile cancer risk and mortality; a population-based study. BJU Int. 2017;119:254–60.

    Article  PubMed  Google Scholar 

  21. Fonseca, AGD, Pinto, JASDA, Marques, MC, Drosdoski, FS, Fonseca Neto, LORD. Estudo epidemiológico do câncer de pênis no Estado do Pará, Brasil. Revista Pan-Amazônica de Saúde, (2010). 1(2), 85–90.

  22. De Sousa IDB, Vidal FCB, Branco Vidal JPC, De Mello GCF, Do Desterro Soares Brandão Nascimento M, Brito LMO. Prevalence of human papillomavirus in penile malignant tumors: Viral genotyping and clinical aspects. BMC Urol (2015) 15:1–6.

  23. Chaux A, Cubilla AL. Advances in the pathology of penile carcinomas. Hum Pathol. 2012;43:771–89.

    Article  CAS  PubMed  Google Scholar 

  24. Paiva GR, de Oliveira Araújo IB, Athanazio DA, de Freitas LAR. Penile cancer: impact of age at diagnosis on morphology and prognosis. Int Urol Nephrol. 2015;47:295–9.

    Article  PubMed  Google Scholar 

  25. Júnior PFM, Silva EHV, Moura KL, de Aquino YF, Weller M. Increased Risk of Penile Cancer among Men Working in Agriculture. Asian Pac J Cancer Prev. 2018;19(1):237–41.

    Article  PubMed  PubMed Central  Google Scholar 

  26. Abdalla GK, Fajardo EF, Gomes BBM, Bianco TM, Salge AKM, de Carvalho EEV, et al. Analysis of knowledge level in Brasilian students about human papillomavirus infection and development of penile cancer. Asian Pac J Cancer Prev. 2017;18:1371–6.

    Article  PubMed  PubMed Central  Google Scholar 

  27. Morrison BF. Risk factors and prevalence of penile cancer. West Indian Med J. 2014;63:559–60.

    Article  CAS  PubMed  Google Scholar 

  28. Chung JM, Park CS, Lee SD. Microbiology of smegma: Prospective comparative control study. Investig Clin Urol. 2019;60(2):127–32.

    Article  PubMed  PubMed Central  Google Scholar 

  29. Hellberg D, Valentin J, Eklund T, Nilsson S. Penile cancer: Is there an epidemiological role for smoking and sexual behaviour? Br Med J (Clin Res Ed). 1987;295:1306–8.

    Article  CAS  Google Scholar 

  30. Nasca MR, Innocenzi D, Micali G. Penile cancer among patients with genital lichen sclerosus. J Am Acad Dermatol. 1999;41(6):911–4.

    Article  CAS  PubMed  Google Scholar 

  31. Thuret R, Sun M, Budaus L, Abdollah F, Liberman D, Shariat SF, Iborra F, et al. A population-based analysis of the effect of marital status on overall and cancer-specific mortality in patients with squamous cell carcinoma of the penis. Cancer Causes Control. 2013;24:71–9.

    Article  PubMed  Google Scholar 

  32. Heller DS. Lesions and Neoplasms of the Penis: A Review. J Low Genit Tract Dis. 2016;20:107–11.

    Article  PubMed  Google Scholar 

  33. Morris BJ, Matthews JG, Krieger JN. Prevalence of Phimosis in Males of All Ages: Systematic Review. Urology. 2020;135:124–32.

    Article  PubMed  Google Scholar 

  34. De Paula AAP, Netto JCA, da Cruz AD, de Freitas Júnior R. Carcinoma epidermóide do pênis: considerações epidemiológicas, histopatológicas, influência viral e tratamento cirúrgico. Revista Brasileira de Cancerologia. 2005;51(3):243–52.

    Article  Google Scholar 

  35. Morris BJ, Kennedy SE, Wodak AD, Mindel A, Golovsky D, Schrieber L, et al. Early infant male circumcision: Systematic review, risk-benefit analysis, and progress in policy. World J Clin Pediatr. (2017) 8;6(1):89–102. PMID: 28224100

  36. Castellsagué X, Bosch FX, Muñoz N, Meijer CJ, Shah KV, de Sanjose S, et al. Male circumcision, penile human papillomavirus infection, and cervical cancer in female partners. N Engl J Med. 2002;346(15):1105–12.

    Article  PubMed  Google Scholar 

  37. Panatto D, Amicizia D, Trucchi C, Casabona F, Lai PL, Bonanni P, et al. Sexual behaviour and risk factors for the acquisition of human papillomavirus infections in young people in Italy: Suggestions for future vaccination policies. BMC Public Health. 2012;12:1–9.

    Article  Google Scholar 

  38. Zequi S de C, 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 

  39. Skeppner E, Andersson SO, Johansson JE, Windahl T. Initial symptoms and delay in patients with penile carcinoma. Scand J Urol Nephrol. 2012;46:319–25.

    Article  PubMed  Google Scholar 

  40. Nasca MR, Innocenzi D, Micali G. Penile cancer among patients with genital lichen sclerosus. J Am Acad Dermatol. 1999;41:911–4.

    Article  CAS  PubMed  Google Scholar 

  41. Attalla K, Paulucci DJ, Blum K, Anastos H, Moses KA, Badani KK, et al. Demographic and socioeconomic predictors of treatment delays, pathologic stage, and survival among patients with penile cancer: A report from the National Cancer Database. Urol Oncol. 2018;36(1):14.e17-14.e24. Epub 2017 Oct 12 PMID: 29031418.

    Article  Google Scholar 

  42. IBGE - Instituto Brasileiro de Geografia e Estatística. Censo demográfico: resultados preliminares - Maranhão. (2010) 8º Recenseamento Geral do Brasil; v. 1, t. 4.

  43. Verhoeven RHA, Janssen-Heijnen MLG, Saum KU, Zanetti R, Caldarella A, Holleczek B, et al. Population-based survival of penile cancer patients in Europe and the United States of America: No improvement since 1990. Eur J Cancer. 2013;49:1414–21.

    Article  CAS  PubMed  Google Scholar 

  44. Hakenberg OW, Dräger DL, Erbersdobler A, Naumann CM, Jünemann KP, Protzel C. The Diagnosis and Treatment of Penile Cancer. Dtsch Arztebl Int. 2018;115(39):646–52.

    Article  PubMed  Google Scholar 

  45. Pow-Sang MR, Ferreira U, Pow-Sang JM, Nardi AC, Destefano V. Epidemiology and natural history of penile cancer. Urology. 2010;76(2 Suppl 1):S2-6.

    Article  PubMed  Google Scholar 

  46. Chaux A, Amin M, Cubilla AL, Young RH. Metastatic tumors to the penis: a report of 17 cases and review of the literature. Int J Surg Pathol. 2011;19(5):597–606.

    Article  PubMed  Google Scholar 

  47. Sonpavde G, Pagliaro LC, Buonerba C, Dorff TB, Lee RJ, Di Lorenzo G. Penile cancer: current therapy and future directions. Ann Oncol. 2013;24(5):1179–89.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  48. De Sousa ID, Vidal FC, Branco Vidal JP, de Mello GC, do Desterro Soares Brandão Nascimento M, Brito LM. Prevalence of human papillomavirus in penile malignant tumors: viral genotyping and clinical aspects. BMC Urol. (2015). 15(1):13.

  49. Backes DM, Kurman RJ, Pimenta JM, Smith JS. Systematic review of human papillomavirus prevalence in invasive penile cancer. Cancer Causes Control. 2009;20:449–57.

    Article  PubMed  Google Scholar 

  50. Miralles-Guri C, Bruni L, Cubilla AL, Castellsagué X, Bosch FX, De Sanjosé S. Human papillomavirus prevalence and type distribution in penile carcinoma. J Clin Pathol. 2009;62:870–8.

    Article  CAS  PubMed  Google Scholar 

  51. Olesen TB, Sand FL, Rasmussen CL, Albieri V, Toft BG, Norrild B, Munk C, 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  CAS  PubMed  Google Scholar 

  52. Nielson CM, Flores R, Harris RB, Abrahamsen M, Papenfuss MR, Dunne EF, et al. Human papillomavirus prevalence and type distribution in male anogenital sites and semen. Cancer Epidemiol Biomarkers Prev. 2007;6:1107–14.

    Article  Google Scholar 

  53. Ang KK, Harris J, Wheeler R, Weber R, Rosenthal DI, Nguyen-Tân PF, et al. Human papillomavirus and survival of patients with oropharyngeal cancer. New England J Med. 2010;363:24–35.

    Article  CAS  Google Scholar 

  54. Koerber SA, Schoneweg C, Slynko A, Krug D, Haefner MF, Herfarth K, Debus J, et al. Influence of human papillomavirus and p16INK4a on treatment outcome of patients with anal cancer. Radiother Oncol. 2014;113:331–6.

    Article  CAS  PubMed  Google Scholar 

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

    Article  Google Scholar 

  56. Gregoire L, Cubilla AL, Reuter VE, Haas GP, Lancaster WD. Preferential association of human papillomavirus with high-grade histologic variants of penile-invasive squamous cell carcinoma. J Natl Cancer Inst. 1995;87(22):1705–9.

    Article  CAS  PubMed  Google Scholar 

  57. Wiener JS, Effert PJ, Humphrey PA, Yu L, Liu ET, Walther PJ. Prevalence of human papillomavirus types 16 and 18 in squamous-cell carcinoma of the penis: A retrospective analysis of primary and metastatic lesions by differential polymerase chain reaction. Int J Cancer. 1992;50:694–701.

    Article  CAS  PubMed  Google Scholar 

  58. Burchell AN, Winer RL, de Sanjosé S, Franco EL. Chapter 6: Epidemiology and transmission dynamics of genital HPV infection. Vaccine. 2006;24:52–61.

    Article  Google Scholar 

  59. 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(2):231–43.

    Article  PubMed  Google Scholar 

  60. De Paula AA, Netto JC, Freitas R Jr, de Paula LP, Mota ED, Alencar RC. Penile carcinoma: the role of koilocytosis in groin metastasis and the association with disease specific survival. J Urol. (2007) 177(4):1339–43; discussion 1343.

  61. Ornellas AA, Nóbrega BLB, Wei Kin Chin E, Wisnescky A, da Silva PCB, de Santos Schwindt AB. Prognostic Factors in Invasive Squamous Cell Carcinoma of the Penis: Analysis of 196 Patients Treated at the Brazilian National Cancer Institute. J Urol (2008) 180:1354–1359.

  62. Lawson JS, Glenn WK, Heng B, Ye Y, Tran B, Lutze-Mann L, et al. Koilocytes indicate a role for human papilloma virus in breast cancer. Br J Cancer. 2009;101:1351–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  63. Cubie, H. A., K. S. Cuschieri, and C. Y. W. Tong. "Papillomaviruses and polyomaviruses: Warts: warts and cancers; polyomavirus associated nephropathy; progressive multifocal leuco-encephalopathy." Medical Microbiology. (2012). 452–463.

  64. Afonso LA, Carestiato FN, Ornellas AA, Ornellas P, Rocha WM, Cordeiro TI, et al. Human papillomavirus, Epstein-Barr virus, and methylation status of p16ink4a in penile cancer. J Med Virol. 2017;89:1837–43.

    Article  CAS  PubMed  Google Scholar 

  65. Brasil P. Estudo Epidemiológico sobre a prevalência nacional de infecção pelo HPV-. (2017). Available at:

  66. Ginindza TG, Dlamini X, Almonte M, Herrero R, Jolly PE, Tsoka-Gwegweni JM, et al. Prevalence of and associated risk factors for high risk human papillomavirus among sexually active women. Swaziland PLoS One. 2017;12:1–18.

    Article  CAS  Google Scholar 

  67. Barbosa LO, Neto JOB, Teixeira-Júnior AAL, Nogueira LR, de Ribamar Calixto JR, Cunha IW, et al. Pseudoangiosarcomatous squamous cell carcinoma: First case report on penis. Transl Androl Urol. 2021;10:1803–6.

    Article  PubMed  PubMed Central  Google Scholar 

  68. Romagosa C, Simonetti S, López-Vicente L, Mazo A, Lleonart ME, Castellvi J, et al. P16Ink4a overexpression in cancer: A tumor suppressor gene associated with senescence and high-grade tumors. Oncogene. 2011;30:2087–97.

    Article  CAS  PubMed  Google Scholar 

  69. Geibler C, Tahtali A, Diensthuber M, Gassner D, Stover T, Wagenblast J. The role of p16 expression as a predictive marker in HPV positive oral SCCHN - A retrospective single-center study. Anticancer Res. 2013;33:913–6.

    Google Scholar 

  70. Cubilla AL, Lloveras B, Alejo M, Clavero O, Chaux A, Kasamatsu E, et al. The basaloid cell is the best tissue marker for human papillomavirus in invasive penile squamous cell carcinoma: A study of 202 cases from Paraguay. Am J Surg Pathol. 2010;34:104–14.

    Article  PubMed  Google Scholar 

  71. Steinestel J, Al Ghazal A, Arndt A, Schnoeller TJ, Schrader AJ, Moeller P, et al. The role of histologic subtype, p16INK4a expression, and presence of human papillomavirus DNA in penile squamous cell carcinoma. BMC Cancer. 2015;15:1–9.

    Article  CAS  Google Scholar 

  72. Zhang J, Zhang H, Xiu Y, Cheng H, Gu M, Song N. Prognostic significance of P16INK4a expression in penile squamous cell carcinoma: A meta-analysis with trial sequential analysis. Biomed Res Int. 2018;2018:9–12.

    Article  CAS  Google Scholar 

  73. Martins VDCA, Cunha IW, Figliuolo G, Rondon HHMF, de Souza PM, Torres Silva FL, et al. Presence of HPV with overexpression of p16INK4a protein and EBV infection in penile cancer-A series of cases from Brazil Amazon. PLoS ONE. 2020;15(5): e0232474.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  74. Ornellas P, Ornellas AA. HPV vaccination is fundamental for reducing or erradicate penile cancer | Opinion: NO. Int Braz J Urol. 2018;44:862–4.

    Article  PubMed  PubMed Central  Google Scholar 

  75. Moch H, Humphrey PA, Ulbright TM, Reuter VE et al. WHO classification of tumors of the urinary system and male genital organs. International Agency for Research on Cancer (2016). 356 p.

  76. Cubilla AL, Reuter V, Velazquez E, Piris A, Saito S, Young RH. Histologic classification of penile carcinoma and its relation to outcome in 61 patients with primary resection. Int J Surg Pathol. 2001;9(2):111–20.

    Article  CAS  PubMed  Google Scholar 

  77. Chaux A, Lezcano C, Cubilla AL, Tamboli P, Ro J, Ayala A. Comparison of subtypes of penile squamous cell carcinoma from high and low incidence geographical regions. Int J Surg Pathol. 2010;18:268–77.

    Article  PubMed  Google Scholar 

  78. Canete-Portillo S, Velazquez EF, Kristiansen G, Egevad L, Grignon D, Chaux A, et al. Report From the International Society of Urological Pathology (ISUP) Consultation Conference on Molecular Pathology of Urogenital Cancers V: Recommendations on the Use of Immunohistochemical and Molecular Biomarkers in Penile Cancer. Am J Surg Pathol (2020) Jul;44(7):e80-e86.

  79. Zhu Y, Zhou XY, Yao XD, Dai B, Ye DW. The prognostic significance of p53, Ki-67, epithelial cadherin and matrix metalloproteinase-9 in penile squamous cell carcinoma treated with surgery. BJU Int. 2007;100:204–8.

    Article  CAS  PubMed  Google Scholar 

  80. Gunia S, Kakies C, Erbersdobler A, Hakenberg OW, Koch S, May M. Expression of p53, p21 and cyclin D1 in penile cancer: p53 predicts poor prognosis. J Clin Pathol. 2012;65:232–6.

    Article  PubMed  Google Scholar 

  81. Zargar-Shoshtari K, Spiess PE, Berglund AE, Sharma P, Powsang JM, Giuliano A, et al. Clinical Significance of p53 and p16ink4a Status in a Contemporary North American Penile Carcinoma Cohort. Clin Genitourin Cancer. 2016;14:346–51.

    Article  PubMed  Google Scholar 

  82. Marinescu A, Stepan AE, Mărgăritescu C, Marinescu AM, Zăvoi RE, Simionescu CE, et al. P53, p16 and Ki67 immunoexpression in cutaneous squamous cell carcinoma and its precursor lesions. Rom J Morphol Embryol. 2016;57:691–6.

    PubMed  Google Scholar 

  83. The Drug Gene Interection Database (DGIdb). TP53 (2021). [Accessed October 20, 2021].

  84. Chahoud J, Gleber-Netto FO, McCormick BZ, Rao P, Lu X, Guo M, et al. Whole-exome sequencing in penile squamous cell carcinoma uncovers novel prognostic categorization and drug targets similar to head and neck squamous cell carcinoma A C. Clin Cancer Res. 2021;27:2560–70.

    Article  CAS  PubMed  Google Scholar 

  85. Aydin AM, Chahoud J, Adashek JJ, Azizi M, Magliocco A, Ross JS, et al. Understanding genomics and the immune environment of penile cancer to improve therapy. Nat Rev Urol. 2020;17:555–70.

    Article  PubMed  Google Scholar 

  86. Jacob JM, Ferry EK, Gay LM, Elvin JA, Vergilio JA, Ramkissoon S, et al. Comparative genomic profiling of refractory and metastatic penile and nonpenile cutaneous squamous cell carcinoma: Implications for selection of systemic therapy. J Urol. 2019;201:541–8.

    Article  PubMed  Google Scholar 

  87. Akagi K, Li J, Broutian TR, Padilla-Nash H, Xiao W, Jiang B, et a. Genome-wide analysis of HPV integration in human cancers reveals recurrent, focal genomic instability. Genome Res. (2014) 24(2):185–99. Epub 2013 Nov 7.

  88. Protzel C, Knoedel J, Zimmermann U, Woenckhaus C, Poetsch M, Giebel J. Expression of proliferation marker Ki67 correlates to occurrence of metastasis and prognosis, histological subtypes and HPV DNA detection in penile carcinomas. Histol Histopathol. 2007;22:1197–204.

    Article  CAS  PubMed  Google Scholar 

  89. Stankiewicz E, Ng M, Cuzick J, Mesher D, Watkin N, Lam W, Corbishley C, Berney DM. The prognostic value of Ki-67 expression in penile squamous cell carcinoma. J Clin Pathol. 2012;65:534–7.

    Article  PubMed  Google Scholar 

  90. De Azambuja E, Cardoso F, De Castro G, Colozza M, Mano MS, Durbecq V, Sotiriou C, et al. Ki-67 as prognostic marker in early breast cancer: A meta-analysis of published studies involving 12 155 patients. Br J Cancer. 2007;96:1504–13.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  91. Li LT, Jiang G, Chen Q, Zheng JN. Predic Ki67 is a promising molecular target in the diagnosis of cancer (Review). Mol Med Rep. 2015;11:1566–72.

    Article  CAS  PubMed  Google Scholar 

  92. Zekan DS, Dahman A, Hajiran AJ, Luchey AM, Chahoud J, Spiess PE. Prognostic predictors of lymph node metastasis in penile cancer: a systematic review. Int Braz J Urol. 2021;47:943–56.

    Article  PubMed  PubMed Central  Google Scholar 

  93. Dos Santos KM, de Barros FD, de Araujo-Mariz C. Câncer de pênis e fatores associados a metástases em linfonodos inguinais. Research, Society and Development. 2021;10(2):e55810212854–e55810212854.

    Article  Google Scholar 

  94. Warli SM, Siregar GP. Over-expression of Ki-67 as a predictor of lymph node metastasis in penile cancer patients. F1000Research (2020) 9:289.

  95. Soares A, de Carvalho IT, da Fonseca AG, Alencar AM Jr, Leite CHB, Bastos DA, et al. Penile cancer: a Brazilian consensus statement for low- and middle-income countries. J Cancer Res Clin Oncol. 2020;146(12):3281–96.

    Article  PubMed  PubMed Central  Google Scholar 

Download references


The authors are grateful for the collaboration of the Hospital do Câncer Aldenora Bello (HCAB), the University Hospital of the Federal University of Maranhão (HUUFMA), and the Rede Brasileira de Serviços Hospitalares (EBSERH) for the authorization to carry out this research. We are grateful to the team at the Oncology Research Center at the Federal University of Pará (Belém, Brazil), for the training given to J.D.P and A.A.L.T.Jr. in HPV genotyping by sequencing. We also thank the Coordenação de Aperfeiçoamento de Pessoal de Nível Superior—Brasil (CAPES), for supporting the training of A.A.L.T.Jr., S.P.C.M., T.M.S.R., W.E.D. (Master's Program in Adult Health) and M.R.C.B (Master's Program in Health Science). We are grateful to the Conselho Nacional de Desenvolvimento Científico e Tecnológico do Brasil (CNPq), for providing scholarships for A.L.S.N. and D.R.D.D. (Graduation course) and to the Fundação de Amparo à Pesquisa e ao Desenvolvimento Científico e Tecnológico do Maranhão (FAPEMA), for the two projects financed by GEBS (Grant numbers: 01339/17 and 01184/19) and for providing scholarships for A.A.L.T.Jr., T.M.S.R. (Master's Program in Adult Health) and T.B.M.S. (Graduation course).


This research was funded by in part by the Coordenação de Aperfeiçoamento de Pessoal de Nível Superior—Brasil (CAPES), grant number 001, and for Fundação de Amparo à Pesquisa e ao Desenvolvimento Científico e Tecnológico do Maranhão (FAPEMA), grant Numbers 01339/17 and 01184/19.

Author information

Authors and Affiliations



A.A.L.T.Jr. and G.E.B.S. are the principal investigators and wrote the first version of the manuscript. T.B.M.S.; T.M.S.R.; D.R.D.D.; W.E.D.; M. R.C.B.; A.L.S.N.; J.R.R.C.; F.S.M.S.N.; A.M.A.J. and R.G.C.F.C. participated in data collection, care and management of the patients. G.E.B.S.; S.P.C.M. and L.O.B performed pathological analysis and interpretation. J.D.P.; A.S.K. and K.G.D. contributed to statistical analysis. A.A.L.T.Jr.; R.B.R.; W.A.S.Jr. and G.E.B.S. contributed to critical revision of important intellectual content of the manuscript. All authors contributed to the writing and approved the final manuscript.

Corresponding author

Correspondence to Gyl Eanes Barros Silva.

Ethics declarations

Ethics approval and consent to participate

The study was approved by the Research Ethics Committee on Humans from the University Hospital of the Federal University of Maranhão (CEP/HU-UFMA) (CAAE 05918918.9.0000.5086) and was conducted following the Declaration of Helsinki. The Research Ethics Committee on Humans from the University Hospital of the Federal University of Maranhão waived the need for the Patient’s consent due to the retrospective nature of the study (Opinion Number 3.122.045).

Competing interests

The authors declare no conflict of interest.

Additional information

Publisher’s Note

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

Supplementary Information

Additional file 1: Supplementary material (S1) Table 1.

Primers used for HPV detection in penile cancer samples Table 2. PCR mix preparation and thermocycling conditions.

Additional file 2: Supplementary material (S2) Table 1.

Statistical data for koilocytosis, molecular detection of HPV, and p16INK4a expression analysis in penile cancer.

Additional file 3: Supplementary material (S3) Table 1.

Statistical data for p53 and ki-67 protein expression analysis in penile cancer.

Rights and permissions

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

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Teixeira Júnior, A.A.L., da Costa Melo, S.P., Pinho, J.D. et al. A comprehensive analysis of penile cancer in the region with the highest worldwide incidence reveals new insights into the disease. BMC Cancer 22, 1063 (2022).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: