Human Papillomavirus Infection in Penile Cancer: Multidimensional Mechanisms and Vaccine Strategies

Penile cancer (PC) is a rare male malignant tumor, with early lymph node metastasis and poor prognosis. Human papillomavirus (HPV) plays a key role in the carcinogenesis of PC. This review aims to summarize the association between HPV infection and PC in terms of virus–host genome integration patterns (the disrupted regions in the HPV and PC genome), genetic alterations, and epigenetic regulation (methylation and microRNA modification) occurring in HPV and PC DNA, as well as tumor immune microenvironment reprogramming. In addition, the potential of HPV vaccination strategies for PC prevention and treatment is discussed. Understanding of the HPV-related multidimensional mechanisms and the application of HPV vaccines will promote rational and novel management of PC.


Introduction
Penile cancer (PC) is a rare malignant tumor in the male genitourinary system, with an estimated 36,068 new cases and 13,211 deaths worldwide in 2020 [1].PC patients often present with inguinal/pelvic lymph nodes or distant metastatic lesions, leading to extremely poor prognosis with a 5-year overall survival rate of 20-35% [2][3][4].Due to its low incidence and lack of systematic studies, the carcinogenic mechanisms of PC remain poorly understood, which immensely hampers the development of new therapeutic regimens for this intractable malignancy.Therefore, it is necessary to summarize the current research on PC pathogenesis in detail to probe novel treatment strategies.
So far, two major carcinogenic pathways have been identified as underlying causes of PC: one related to inflammation, phimosis, and lichen sclerosis, and the other to human papillomavirus (HPV) infection [5].The latest World Health Organization tumor classification of PC had suggested categorizing penile squamous cell carcinomas based on those not associated with HPV and those associated with HPV, based on their HPV infection status [6].Approximately 50.8% of PCs could be classified as HPV-related tumors [7].Moreover, many studies have found that HPV is the essential pathogenic factor of PC, and HPV status is correlated with clinical characteristics and prognosis of PC patients after receiving chemotherapy, radiotherapy, or immunotherapy [8][9][10], which is similar to other HPV-related carcinomas (vulva cancer, vagina, etc.) [11][12][13].Thus, understanding the underlying carcinogenic mechanisms of HPV infection in PC is critical for the prevention and treatment of this rare virus-induced tumor.
This review emphasizes the virus-host genome integration patterns, DNA alterations, epigenetic regulation, and tumor microenvironment reprogramming caused by HPV infection in PC, as well as prophylactic/therapeutic vaccines targeting HPV.secretion, thus providing the possibility of immunotherapy for HPV16-positive PC [30].Another study comparing the HPV detection performance between assays targeting L1 and those targeting E6/E7 in 50 PC samples showed satisfactory results with high agreement (over 90% for HPV16), implying infrequent disruption events occurring in the L1, E6, or E7 regions [31].
In HPV-induced cancers, the inactivation of p53 and pRb is usually attributed to the presence of HPV E2 disruption.However, the results of some other PC studies are inconsistent with this assumption [32][33][34].Therefore, other disruptive sites in the HPV/PC genome and corresponding gene regulation should also be understood in detail.Currently, two studies have analyzed the integration patterns among the Chinese population.Huang et al. conducted a high-throughput viral integration analysis in 41 HPV-positive PC tissues and found 35 cases harboring virus-host DNA integration.They revealed another pattern in which the E2 region was not involved in the integration event, as it is in cervical cancer, with 31.4% of PC cases retaining an intact E2 gene, suggesting that the inhibitory effect of E6/E7 on tumor suppressors derived from E2 disruption might not be as pronounced in PC.Interestingly, E1, the region close to E2, was found to be the most enriched with breaking events.In addition, this work indicates that several viral integration-involved genes (KLF5, CADM2, etc.) and pathways (MAPK, JAK/STAT, etc.) in the host possess carcinogenic potential, and that integration-induced down-regulation of KLF5 and CADM2 may influence the proliferative and invasive capacity of PC cells in vitro [35].Another similar study performed HPV hybridization capture sequencing on 103 PC samples and found that the most frequent subtype was HPV16, with the highest frequency of virushost DNA recombination events in this subtype.They also detected massive E1 break points instead of E2 or L1/E6/E7 disruption, which caused 19 and 42 integration sites in intergenic and intragenic regions among 17 host chromosomes [36].As the E1 is an ATP-dependent helicase important in HPV DNA replication [37], the frequent integration in the E1 region observed in both of the above studies appears to inhibit the HPV replication cycle.Meanwhile, several other articles reported that E1 disruption was associated with viral immortalization capacity and poor prognosis in HPV-related cancers [38,39].Thus, more research concerning the function of E1 should be conducted to explain how disruption of this region influences the development and progression of PC.
Table 1 presents a list of HPV-PC genome integration patterns [23,[28][29][30]35,36].Increasing histopathologic grading (average 2.7 integration sites per patient in grading G1 patients, and average 3.83 integration sites in grading G2 patients) [36] 4. DNA Alterations 4.1.HPVDNA Variations Although HPVs have been classified into hundreds of "subtypes", more specific classifications within a certain subtype are established based on differences in genome sequence, referred to as "HPV variants" [40].These variants often have distinct distributions in different geographic regions and are named according to their specific geographical origin [41].Currently, the DNA sequence heterogeneity of the HPV16 E6, LCR, and L1 genes has been extensively investigated.As a result, five branches of phylogenetic variants have been defined: European (E), Asian (As), Asian American (AA), African-1 (Af-1), and African-2 (Af-2), based on sequence variation using 408 cervical cancer specimens from 22 countries across five continents [42].Subsequently, these general categories have been validated in specific countries, such as As variants in China [43,44] and E variants in Portugal [45].
The alterations in the HPV16 DNA sequence have been implicated in the malignant transformation of squamous epithelium and in the progression of cancer.In cervical cancer, the frequency of AA and E variants was found to be significantly higher compared to normal cervical samples (50.3% vs. 1.1%) [46].Specifically, the L83V, a specific E variant of HPV16 E6, was found to be more prevalent in high-grade cervical malignancies [47].Furthermore, patients with vulvar squamous cell carcinoma who harbored the E-G131 variant had poorer survival than those with the E-p or E-G350 variants [48].
Several studies have focused on HPV variants in PC.Tornesello et al. identified European and non-European variants in 44.4% and 55.6% of 18 HPV16 PC samples from the Italian population, with the AA variant achieving the highest frequency [49].This result suggests that non-European variants, specifically AA, may play a more important role in PC tumorigenesis than European variants, similar to what has been observed in cervical carcinoma [50][51][52].Kalantari et al. also found a higher proportion of AA variants (53%) compared to E variants (47%) in 17 PC tissues from Brazilian patients [29].Interestingly, another case report from Japan indicated the presence of an As variant (T178G) in the HPV16 E6 gene in a married couple (wife with vulvar carcinoma and husband with PC) [53].However, the dominance of non-European variants was not observed in Mexican patients with PC, as they had a higher prevalence of E variants (92%) than AA variants (8%) [54].This contradiction could be attributed to racial differences in these studies, as well as other factors such as sample size and variations in testing methods.Therefore, larger studies enrolling PC patients from different regions worldwide are needed to assess the association between HPV16 variants and PC development.
It is worth noting that specific sequence mutations might occur in a particular variant of the virus, resulting in an increased risk of oncogenesis.Tornesello et al. revealed LCR, E6, and E7 mutations in a certain HPV16 variant, Af-1.These mutations were identified at nucleotide positions 7714 (T→A), 7461 (C→G), 7489 (G→A), 7521 (G→A), 7764 (C→T), 7786 (C→T), and 7833 (G→T) in LCR, positions 132 (G→C), 143 (C→G), 145 (G→T), 286 (T→A), 289 (A→G), and 335 (C→T) in E6, as well as positions 789 (T→C) and 795 (T→G) in E7.Moreover, the LCR in mutated Af-1 had more promoter activity and mutated LCRdriven E6/E7 showed higher transforming activity in vitro than those in classical Af-1, highlighting the critical role of sequence alterations in PC pathogenesis [55].Similarly, pointmutation-induced LCR enhancer activity has also been detected in oral cancer with HPV infection [56].Hence, we consider that the enhanced expression of E6/E7 caused by specific mutations in the LCR may stimulate PC tumorigenesis.However, more studies are needed to identify the detailed and precise DNA mutational patterns in specific HPV16 variants.

HPV + PC Genome Alterations
The accumulation of spontaneous somatic mutations in normal cells can lead to the development of corresponding cancers [57].Specifically, persistent high-risk HPV infection has been identified as a potential trigger for genomic instability and the accumulation of somatic mutations [58,59].A study investigating the association between HPV infection and host gene mutations in head and neck squamous cell carcinoma revealed two distinct subtypes, HPV-KRT and HPV-IMU, which appear to be associated with different mutation patterns.HPV-IMU shows a tendency towards chr16q loss, while HPV-KRT is more prone to chr3q gain and PIK3CA mutation [60].Another study also demonstrated the potential impact of HPV status on genetic alterations in cervical and head and neck cancers, where the LRP1B mutation was found to be more frequently observed in samples with HPV A9 or 16 integrations [61].
For PC, some works have also been conducted to analyze the role of viral infection in host genomic mutations.A case report first reported a c-ras Ha missense mutation at codon 61 in an advanced lymph node metastasis of HPV18 + PC, occurring 7 years after the initial lesion.This finding suggests that this specific genetic alteration may occur during long-term HPV infection [62].Furthermore, point mutations were detected in p53 exons 6, 7, and 8 in HPV + PC, showing that high-risk type HPVs and mutated p53 may coexist in these tumors [63].Silva et al. investigated genome-wide profiling of copy number alterations (CNAs) using array comparative genomic hybridization analysis in 28 HPV + PC cases, and revealed 38 altered cytobands with 2314 CNAs, 89.5% of which were located at the HPV integrated positions in the host genome and found that the gain on 14q12 was correlated with HPV16 and multiple infections [64].Another study involving 30 Latin American PC samples with HPV infection suggested that mutations in classical cancer-associated (TP53, NOTCH1, CDKN2A, etc.) and other novel genes (CIC, KMT2C, CR1, etc.) were common, and the majority of them localized at HPV integration sites.Additionally, five cancer-associated genes (NOTCH1, MYC, NUMA1, PLAG1, and RAD21) were identified as the most frequently with gain, and SMARCA4 with loss [65].Nazha et al. found KMT2C mutations and FGF3 amplification in 33% and 30.8% of HPV16/18 + samples, respectively, but not in HPV16/18 -tumors.Moreover, 30.8% of HPV16/18 + tissues were classified as having a high tumor mutational burden (≥10 mutations/Mb) [66].Other studies have identified additional genomic alterations in HPV + PC patients, including heterozygous deletions and single nucleotide change in STK11 [67], as well as mutations in the TERT promoter [68].These aforementioned studies indicate that the host DNA alterations, possibly caused by HPV infection, may contribute to the development of PC.
Table 2 presents a list of the DNA alterations in HPV DNA and the host genome in HPV + PC [29,49,[53][54][55][62][63][64][65][66][67].DNA methylation, which has been recognized as a hallmark of cancer, is a nonmutational epigenetic reprogramming mechanism [69].This modification occurs primarily in CpG islands near or within gene promoters and causes the addition of a methyl group at a specific position on the cytosine ring-carbon 5 site, thereby suppressing transcriptional activation of the corresponding gene [70].Several tumor suppressor genes commonly exhibit methylation, such as APC in esophageal adenocarcinoma, BRCA1 in breast cancer, and ER in prostate cancer.Methylation usually inhibits the expression of these genes, thus inducing carcinogenesis, tumor progression and treatment resistance [71][72][73].
Methylation events also occur in HPV DNA.In cervical cancer, a higher incidence of HPV16 DNA methylation (47%) was detected at E2-binding sites (E2BS-1, E2BS-2) and the Sp1-binding site (Sp1BS) of the LCR region (regions involved in regulating E6/E7), compared to in normal tissues and cervical intraepithelial neoplasia (CIN) [74].Another study suggested a high methylation rate (56.1%) of CpGs in E2BS-1 near the P97 promoter in the presence of intact E2 proteins, leading to the loss of function of E2 proteins that repress E6/E7 expression [75].Clinical data further support the value of HPV DNA methylation.A prospective cohort study showed a significant increase in the 10-year cumulative incidence rate of CIN grade 3 (CIN3) or severe (CIN3+) in women with a methylation panel (CpG 5602, 6650, 7034, 7461, 31, and 37) of HPV 16, which exhibited 85.7% of sensitivity and 78.4% of specificity for cumulative CIN3+ [76].Kalantari et al. analyzed the methylation of 95 HPV16 molecules derived from 19 PC samples and found abundant CpG islands methylated in PC samples, with 58% in L1 (position 7091, 7136, and 7145) and 22% in the 5 part of the LCR (position 7270, 7428, 7434, 7455, and 7461).In addition, they identified several DNA methylation events occurring in five CpGs (position 7535, 7554, 7677, 7683, and 7695) that overlapped with the transcriptional enhancer and six CpGs (position 7862, 31, 37, 43, 52, and 58) that were part of the E2-and SP1-binding sites of the E6 promotor, suggesting similar epigenetic mechanisms associated with HPV16 infection between PC and cervical cancer [29].Interestingly, the study also suggests that five PC samples with a single copy of the integrated HPV16 genomes exhibit high methylation in the L1 region.This finding indicates a potential correlation between HPV DNA methylation and virushost genome integration.Consistently, methylation events of HPV16 DNA have been recognized as biomarkers of HPV genome integration into the host chromosome in cervical cancer and vulval intraepithelial neoplasia [77][78][79].The detailed mechanism underlying this intriguing association will be further investigated in the future.

HPV + PC Genome Methylation
Methylation events occur in HPV + PC cellular genome in addition to in HPV DNA.Generally, p16 ink4a (a cyclin-dependent kinase inhibitor) is inhibited by Rb protein [80,81].Thus, HPV E7-induced blocking of Rb function can increase p16 ink4a expression in PC via a feed-back control [82][83][84].Even though up-regulation of p16 ink4a expression has been observed during HPV infection, several studies have demonstrated that the p16 ink4a gene can be methylated.The proportion of p16 ink4a methylation events in HPV + PC samples ranged from 10% to 46.8% [85][86][87], which indicated that methylated p16 ink4a may be involved in tumorigenesis and progression in PC.Since the p16 ink4a protein induces cell cycle arrest [88], presumedly, methylation of the p16 ink4a gene counteracts the proliferation inhibition caused by elevated p16 ink4a expression due to the HPV E7/Rb axis, thereby promoting PC progression.
Other studies have provided additional evidence regarding the methylation of host DNA in PC patients with HPV infection.One of these studies reported that 88.0% of PC cases infected with high-risk HPV had elevated global 5-methylcytosine levels (>60%), while 30.4% had increased global 5-hydroxymethylcytosine levels (>60%), suggesting that 5-methylcytosine, rather than 5-hydroxymethylcytosine, may serve as a stable epigenetic marker [89].The findings in other cancers suggest the prognostic and drug resistance value of increased global 5-methylcytosine [90,91].Kuasne et al. analyzed the genome-wide methylation profiles in 25 PC samples and found distinct levels of methylation between the HPV + (7 samples) and HPV − (18 samples) groups.Specifically, 65 genes were found to be hypermethylated in HPV + cases, including CD70, HN1, FZD5, and FSCN1, which are involved in tumor migration, invasion, immunomodulation and other processes [92][93][94][95][96]. Another study detected methylation of thrombospondin-1 (TSP-1) and RAS association domain family 1A (RASSF1A) in 54.5% and 27.3% of HPV + PC tissues, respectively, and established a correlation between these two hypermethylation events and pathological variables [97].It has been reported that inhibition of TSP-1 and RASSF1 can induce growth in gastric cancer cells and enhance colorectal cancer stemness [98,99].Therefore, methylation of these two genes may promote HPV + PC progression through similar pathways.Moreover, epigenetic silencing of FHIT via methylation has been observed in 100% of PC samples with viral infection [100] and the suppression of FHIT appears to activate tumorigenesis in hepatic lesions [101].These aforementioned studies highlight the potential epigenetic effects on HPV + PC.   5.2.miRNA/mRNA Regulation in HPV + PC MicroRNAs (miRNAs), a class of small endogenous non-coding RNAs, promote cancer development by negatively regulating gene expression [102].Studies in the field of HPVassociated cancers, including cervical, head and neck, and anal carcinomas, have confirmed that HPV infection mediates miRNA regulation and influences cancer cell proliferation and migration [103][104][105][106].
Several studies have revealed changes in miRNA levels in PC patients infected with HPV.MiR-146a, known for its tumor suppressor function through binding to the 3' UTRs of EGFR and NOTCH1 mRNA in breast cancer and glioma [107,108], was found to present lower levels in HPV + PC cases compared to those in HPV − PC cases, and to correlate with an expected elevation of EGFR expression.Further in vitro experiments verified a negative relationship between the levels of HPV16 E6 and miR-146a, which affects the proliferation of keratinocytes [109].TP53 and RB1 are acknowledged downstream regulators of HPV E6 and E7 in PC.Silva et al. identified 13 up-regulated miRNAs in 22 HPV + PC tissues, 11 of which could target the TP53 and RB1 genes, suggesting an epigenetic mechanism of HPV infection for TP53/RB1 suppression [110].In addition, this group analyzed miRNA regulation from the perspective of CNVs and identified 269 miRNAs in 30 altered cytobands with CNAs and 898 genes that could be potentially targeted by these miRNAs.Furthermore, five pathways are affected by the miRNA/mRNA interactions: the Hippo signaling pathway, lysine degradation, mucin-type O-Glycan biosynthesis, prion disease, and proteoglycans in cancer.They also highlighted potential PC-targeting drugs, such as cetuximab and celecoxib [64].Canto et al. identified 4942 DNA variants occurring in 3277 genes in 30 HPV + PC samples, with 5.5% of the alterations found within the UTRs.These alterations included miRNA targets in the UTRs of four cancer-associated genes (CARD11, CSMD3, KDR, and TLX3).It was speculated that miRNA-directed down-regulated expression of these genes might be impaired due to UTR variants during HPV infection [65].Thus, epigenetic regulation by HPV-associated miRNAs contributes to disease progression in PC.
Table 4 presents a list of the miRNA/mRNA regulations in HPV + PC [64,65,109,110].

Immune Microenvironment Reprogramming of PC with HPV Infection
It has been reported that HPV infection affects the immune microenvironment in cervical cancer through HPV-stimulated immune evasion via the reduction in immunostimulatory IFN-α of natural killer cells and the elevation of M1 to M2 transition in macrophages [111,112].Additionally, it was found that HPV-16 DNA load was associated with the infiltration of CD8 + and PD1 + tumor infiltrating lymphocytes, which conferred a favorable prognosis for anal and oropharyngeal cancer [113,114].Recently, several researchers have focused on understanding the immune regulation associated with HPV to develop immunotherapies for PC.
Tumor infiltrating immune cells are the key component of the immune microenvironment.An early study analyzed the influence of HPV on the amount of Langerhans cells, a kind of immature dendritic cell that functions as immune surveillance, in penile epithelial cells.The results showed a significant reduction in Langerhans cells occurring in HPVinfected genital tracts, while local treatment with 5-fluorouracil restored a normal density of Langerhans cells, suggesting that HPV may promote immune evasion [115].Another study found a significantly lower CD4 + /CD8 + ratio in HPV-positive lymph nodes in PC patients, exhibiting an obvious expansion of CD8 + cytotoxic T lymphocytes [116].Chu et al. explored the immune phenotype characteristics including T cells (CD4 + , CD8 + , or FOXP3 + ) and macrophages (CD68 + or CD206 + ) and found that PC samples with high-risk subtypes of HPV harbored a higher density of stromal Granzyme B + CD8 + T cells compared to HPV − tumors, indicating a highly immunogenic property of HPV [117].Similarly, Lohneis et al. found a higher number of tumor-infiltrating T cells (T helper 1 and cytotoxic T cells) in HPV + PC tissues rather than in HPV − PC samples, further confirming that HPV can induce immune activation.Intriguingly, a type of immune-suppressing T cell, CD4 + Foxp3 + Treg cells, was detected more frequently in virus-infected tissues, indicating the complexity of HPV-induced immune regulation [118].However, some studies did not observe significant differences in tumor-infiltrating immune cells between HPV + and HPV − cases, which could be attributed to patient heterogeneity, sample size, testing methods, and other factors [119,120].
Immunotherapy, including programmed cell death protein 1 (PD-1)/programmed cell death-ligand 1 (PD-L1) inhibitors, has become a preferred option in combination with chemotherapy or targeted treatments [121].PD-L1 expression has been reported as positively associated with tumor growth, progression, and poor prognosis [122].Therefore, it is important to investigate whether HPV infection affects PD-1/PD-L1 expression.Three studies received consistent results: PD-L1 expression appears to be lower in HPV + PC cells than in HPV − PC cells.Specifically, the rate of PD-L1 positivity was 32.7-53.8% in HPV + PC samples, whereas it was 49.4-64.4% in those without HPV infection [119,120,123].These findings partially explain why HPV + PC patients generally have a better prognosis than HPV − PC patients [8][9][10].Additionally, Chu et al. found that the HPV + PC exhibited a higher density of intratumor PD-1-positive T cells than the HPV − tumors, indicating immune evasion and higher immunotherapeutic potential for HPV + PC [117].In fact, a retrospective study using a PD-1 inhibitor in advanced PC patients showed potentially superior anti-tumor efficacy in cases with HPV infection [124].The aforementioned works suggest that immune checkpoint inhibitors may confer benefits to HPV + PC patients due to activation of PD1 + cytotoxic T lymphocytes.Certainly, this assumption needs further validation in large-scale clinical trials.
Figure 1 presents the multidimensional mechanisms regarding HPV infection in PC.
perior anti-tumor efficacy in cases with HPV infection [124].The aforementioned works suggest that immune checkpoint inhibitors may confer benefits to HPV + PC patients due to activation of PD1 + cytotoxic T lymphocytes.Certainly, this assumption needs further validation in large-scale clinical trials.
Figure 1 presents the multidimensional mechanisms regarding HPV infection in PC.

HPV Vaccine in PC
Considering the nature of virus-induced carcinogenesis in HPV-associated malignancies, vaccines targeting HPV may be an effective strategy for the prophylaxis of these cancers.HPV vaccination has been recommended worldwide for females and in certain countries for males to prevent HPV infection and the development of HPV-related neoplasms.Numerous studies have confirmed a significant reduction in the incidence of HPV infection and HPV-related carcinomas in the cervix, anus, and oropharynx in vaccinated individuals compared to those in the unvaccinated population [125][126][127].Therefore, it is crucial to assess the preventive effects of HPV vaccination on viral infection and carcinogenesis in the penis.
Routine HPV vaccination is currently recommended for boys and men by the U.S. Advisory Committee on Immunization Practices [128].Recently, Winer et al. reported the results of their investigation in a subset of the male population-men having anal or oral sex with men and transgender women-using either the quadrivalent or nonavalent HPV vaccine.Inspiringly, the application of the HPV vaccine effectively prevented penile HPV infection in participants vaccinated at age ≤ 18 years, with 85% of those vaccinated harboring the observed efficacy against HPV 6, 11, 16, and 18 subtypes.Nonetheless, no significant difference in the incidence of HPV infection in the penis was found between all vaccinated persons (12.1%) and those with no or unknown vaccination (15.6%) [129].These data reflect a promising application of the HPV vaccine at a lower age for the prevention of penile virus infection, but a long-term follow-up will be important to further determine whether HPV vaccination is associated with decreased PC incidence.Another large clinical trial enrolled 4065 healthy males to assess the impact of a quadrivalent HPV vaccine on protection against viral infection and corresponding diseases, including penile intraepithelial neoplasia and PC.The results showed 90.4% efficacy against the HPV 6, 11, 16, or 18 subtypes and 100% efficacy against the HPV 16 subtype (the most dominant HPV genotype in PC) in the per-protocol population.Additionally, there were no cases of penile/perianal/perineal intraepithelial neoplasia in the vaccine group, but three cases in the control group, hinting to some extent that HPV vaccination can substantially reduce the incidence of PC [130].
Other studies also suggested similar results-potential effectiveness of HPV vaccine against penile/perianal/perineal HPV infections and corresponding lesions [131,132].Surely, it is noted that HPV infection rate in PC is not as high as cervical cancer and the carcinogenic mechanisms are divided into HPV-related and non-HPV-related mechanisms, so we guess the usage of HPV vaccine may be suitable for certain population with the high risk of HPV infection, such as men having sex with men and transgender women.
In addition to its prophylactic effects, HPV vaccination may confer therapeutic action against HPV-related cancers.Therapeutic vaccinations targeting HPV 16/18 E6 and E7 proteins have showed clinical responses in anal and cervical intraepithelial neoplasia [133,134].For PC, a preclinical study demonstrated HPV-L1-protein-induced T cell reactivity in lymphocytes from the lymph nodes of two HPV + PC patients.Strong activation of CD4 + T h and CD8 + T lymphocytes against the HPV vaccine was detected, and dose-dependent responses with HLA-DR, a marker of T cell activation, were observed against gardasil at 0.1, 0.3, and 1.0 µg/mL, supporting T cell stimulation depending on the HPV vaccine [116].In addition, another brief report presented an immunologic response and satisfactory efficacy against penile intraepithelial neoplasia using the HPV vaccine in combination with imiquimod, exhibiting a synergistic immunomodulatory effect on lesion clearance in the penis [135].However, therapeutic HPV vaccines for PC are still in the preclinical stage and more research is needed on vaccine strategies for PC treatment.

Conclusions
HPV, particularly genotype 16, is involved in the malignant transformation of normal penile epithelium and the progression of PC in multiple ways.Regarding virus-host genome integration, the E1 region is disrupted more frequently in PC, and the integration positions in the PC genome tend to indicate the corresponding regulation of various tumor suppressors or oncogenes.In terms of DNA alterations, HPV DNA and the PC genome exhibit distinct characteristics.Non-European variants are more common in PC, sometimes with specific point mutations in the Af-1 region.DNA alterations in the HPV + PC genome, such as CNAs and point mutations (c-ras Ha , P53, NOTCH1, etc.), are potentially affected by viral infection.Epigenetic changes are also significant in PC development.CpG methylation has been observed within the regions of HPV L1 and LCR, including the E2BS-and Sp1binding sites.Specific methylation events and up-or down-regulation of miRNAs in HPV + PC often induce tumor migration and therapy resistance via regulating the expression of corresponding genes.From the perspective of the PC immune microenvironment, anti-tumor immune cell types (especially CD8 + T cells) increase with the up-regulation of PD1 expression on T cells.This enhancement potentially facilitates the efficacy of immunotherapy in PC patients.Finally, HPV vaccination may serve as an effective strategy for preventing HPV infection and PC development, and as a treatment option for HPV + PC patients.
Surely, some limitations exist in this review.First, the number of eligible studies in some sections (HPV DNA methylations and HPV DNA mutations) is relatively inadequate due to the low incidence of PC, leading to potentially insufficient conclusions.Second, PC populations and the detecting technologies of HPV-PC genome integration vary from study to study, which may bias the results in disagreement.Third, the epigenetic effects from HPV infection have only been detected in some studies, which is not enough to clarify and validate the impact of methylation and miRNA/mRNA regulation on PC progression.
In future studies, a more detailed investigation of the integration patterns of the HPV-PC genome, PC gene methylation, and miRNA/mRNA regulation, and consequently the expression regulation of host genes and downstream molecular changes during HPV infection, is warranted to gain a deeper insight into the pathogenesis of HPV-related PC.Moreover, the research focusing on the development of therapeutic vaccines targeting specific HPV DNA regions is indispensable for novel treatment options of PC patients.

Figure 1 .
Figure 1.Presents a summary of the multidimensional mechanisms regarding HPV infection in PC from the aspects of virus-host genome integration patterns, gene alterations, epigenetic regulation, and tumor-immune microenvironment reprogramming.

Figure 1 .
Figure 1.Presents a summary of the multidimensional mechanisms regarding HPV infection in PC from the aspects of virus-host genome integration patterns, gene alterations, epigenetic regulation, and tumor-immune microenvironment reprogramming.

Table 2 .
DNA alterations in HPV DNA and HPV + PC DNA.

Table 3 .
Methylation regulations in HPV DNA and HPV + PC DNA.