Rosario Castro-Oropeza
Patricia Piña-Sánchez- Molecular Oncology Laboratory, Oncology Research Unit, Oncology Hospital, IMSS National Medical Center, Mexico City, Mexico
Human Papillomavirus (HPV) is an oncogenic virus that causes the highest number of viral-associated cancer cases and deaths worldwide, with more than 690,000 new cases per year and 342,000 deaths only for cervical cancer (CC). Although the incidence and mortality rates for CC are declining in countries where screening and vaccination programs have been implemented, other types of cancer in which HPV is involved, such as oropharyngeal cancer, are increasing, particularly in men. Mutational and transcriptional profiles of various HPV-associated neoplasms have been described, and accumulated evidence has shown the oncogenic capacity of E6, E7, and E5 genes of high-risk HPV. Interestingly, transcriptomic analysis has revealed that although a vast majority of the human genome is transcribed into RNAs, only 2% of transcripts are translated into proteins. The remaining transcripts lacking protein-coding potential are called non-coding RNAs. In addition to the transfer and ribosomal RNAs, there are regulatory non-coding RNAs classified according to size and structure in long non-coding RNAs (lncRNAs), circular RNAs (circRNAs), and small RNAs; such as microRNAs (miRNAs), piwi-associated RNAs (piRNAs), small nucleolar RNAs (snoRNAs) and endogenous short-interfering RNAs. Recent evidence has shown that lncRNAs, miRNAs, and circRNAs are aberrantly expressed under pathological conditions such as cancer. In addition, those transcripts are dysregulated in HPV-related neoplasms, and their expression correlates with tumor progression, metastasis, poor prognosis, and recurrence. Nuclear lncRNAs are epigenetic regulators involved in controlling gene expression at the transcriptional level through chromatin modification and remodeling. Moreover, disruption of the expression profiles of those lncRNAs affects multiple biological processes such as cell proliferation, apoptosis, and migration. This review highlights the epigenetic alterations induced by HPV, from infection to neoplastic transformation. We condense the epigenetic role of non-coding RNA alterations and their potential as biomarkers in transformation’s early stages and clinical applications. We also summarize the molecular mechanisms of action of nuclear lncRNAs to understand better their role in the epigenetic control of gene expression and how they can drive the malignant phenotype of HPV-related neoplasia. Finally, we review several chemical and epigenetic therapy options to prevent and treat HPV-associated neoplasms.
1 Epidemiology of Human Papillomavirus Neoplasms
In 2020, 19.3 million new cancer cases and 10 million deaths were estimated worldwide (Sung et al., 2021). About 13% of cancer cases are attributed to infections (2.5 million), being Human Papillomavirus (HPV) the second most frequent infectious agent with 31% of cases (de Martel et al., 2020). The HPV-associated neoplasms are the cancer of the cervix, anus, vulva, penis, vagina, and head and neck squamous cell carcinomas (HNSCC). The health problem of HPV-associated neoplasms is concentrated mainly in women since 80% is due to cervical cancer (CC) (de Martel et al., 2020). The main contributions to the knowledge and control of HPV-associated cancer have been developed from CC studies. Persistent high-risk HPV (HR-HPV) infection is responsible for nearly 100% of cervical carcinomas, anal carcinomas, and to a lesser extent HNSCCs (Figure 1).
FIGURE 1. HPV-attributable cancer. (A) The graph shows the proportion of cancer attributable to HPV worldwide (B) Number of new cancer cases attributable to HPV. Data from de Martel et al. (2020).
Recent data estimated that CC alone accounted for more than 600,000 new cases and more than 340,000 deaths (Sung et al., 2021). HPV infection is the most common sexually transmitted infection; it is estimated that at least 80% of people have been infected with HPV before 45 years (Chesson et al., 2014).
HPV belongs to the Papillomaviridae family, composed of 52 genera (ICTV 2020). According to the International Human Papillomavirus Reference Center, 222 distinct genotypes are registered (de Villiers et al., 2012). The International Agency for Research on Cancer (IARC) classifies 12 genotypes (16, 18, 31, 33, 35, 35, 39, 45, 45, 51, 52, 56, 58, and 59) as type 1 carcinogens, genotype 68 as a possible carcinogen and several more as probable carcinogens (26, 53, 66, 67, 70, 73, 82, 30, 34, 69, 85, and 9) (IARC 2012). Due to epidemiological and carcinogenic evidence, high-risk genotypes are recognized (the 12 genotypes classified as carcinogens type 1, and the genotype 68), constituting 96% of CC cases (Arbyn et al., 2014).
According to data collected worldwide, the 10 most common genotypes in CC are 16, 18, 45, 33, 58, 52, 31, 35, 39, and 59 (Bruni et al., 2021). While in HNSCC the most frequent HPVs in decreasing order are: 16, 18, 33, 35, 52, 45, 39, 58, 31, 53, and 56 (Tumban 2019). Differences in the frequency of genotypes have been described, which may be related to ethnicity and detection methods. However, because of their high frequency, HPV16 and 18 genotypes have been mostly studied from a functional point of view. A consistent finding is that HPV16 has been the most frequent genotype in cancer at different anatomic sites. The oncogenic potential of HPV16 is associated with persistent infection and a high viral load. Also, different variants of the same genotype possess distinct carcinogenic capabilities.
Recent metagenomics and next-generation sequencing (NGS) studies describe a wide variety of HPV genotypes at various anatomical sites in both healthy tissues and malignant neoplasms. In 69% of healthy subjects, the presence of some HPV genotype has been identified, either in the intestine, mouth, skin, or vagina, being more frequent in the skin (61%) and vagina (41%) (Ma et al., 2014). In anogenital warts, HPV 6, 7, and 11 have been identified as the most frequent genotypes (Xu H et al., 2019). HPV genotypes show greater diversity in low-grade squamous intraepithelial lesions (LSIL) than HSIL, with HPV39 dominating in LSIL and HPV16 in HSIL (Shen-Gunther et al., 2021).
2 Structure and Function of Human Papillomavirus Genome
Human papillomaviruses are epitheliotrophic species-specific viruses. HPV viral cycle is closely related to the differentiation of stratified squamous epithelium. About its genomic traits, HPV has a circular double-stranded genome of approximately 8,000 bp. In general, its genomic structure is composed of three regions: the early region (E), which encodes proteins responsible for the regulation of transcription and replication; the late region (L), which encodes structural proteins; and the long control region (LCR) or upstream regulatory region (URR). Transcription is polycistronic as in other viral species. In HPV16, the P97 and P670 promoters are the major transcription regulators. The P97 promoter controls the expression of the E6, E7, E1, E2, E4, and E5 polycistron in undifferentiated cells. The second promoter regulates the transcription of E1^E4, E5, L1, and L2 and is activated in differentiated cells (Wooldridge and Laimins 2008). The general functions and interactions of HPV-encoded proteins are briefly described below.
E1: This viral protein has the largest and most conserved ORF of the papillomavirus. The main function of E1 is to bind to the origin of replication to maintain viral replication and is responsible for maintaining the episomal state. E1 interacts with proteins of the replication machinery (DNA polymerase, RPA, Topoisomerase I and II, PCNA). Interestingly, it also interacts with Histone H1 to displace it from DNA. In addition, E1 has been shown to bind with epigenetic regulators such as HATs, HDACs, and SWI/SNF for chromatin remodeling (Bergvall et al., 2013).
E2: E2 is a key regulatory protein encoded by all papillomaviruses. Interestingly, the short forms of splicing can act as repressors of E2. E2 can interact with specific viral DNA sequences, mainly methylation-dependent sequences, such as E2BS sites located within the LCR region. In addition, E2 is very important for the recruitment of E1 to the replication origin. E2 protein has a repressor role in E6 and E7 transcription, which has been associated with increased apoptosis and decreased cell proliferation (Blachon et al., 2005), in part by the consequent increase of p53, so that the disruption of E2 allows the overexpression of E6 and E7 oncoproteins, which promotes the development of the neoplasm.
On the other hand, E1 and E2 are known to be involved in the process of viral genome insertion into the host genome. However, only the E1 gene is involved in low-grade cervical intraepithelial lesions, and in high-grade lesions, the E1 and E2 genes are also disrupted. The hinge region of the E2 gene has been reported to be unstable making it possibly the most common break when HPV DNA integrates into the host genome (Cricca et al., 2009). Therefore, E2 and E6 have been used as surrogate markers to discriminate CIN2+ in CC (Chen et al., 2020).
E2 functions as a histone code reader by interacting with the bromodomain containing 4 (Brd4) and interacting with acetylated histones in chromatin. E2 has been shown to bind to other chromatin regulators such as HAT (p300) HDAC 1, 2, and 3, KDM5C, SWI/SNF, TIP60, EP400, and SMCX (McBride 2013; Fontan et al., 2020).
E8^E2: E2 splicing isoforms are transcribed from the promoter E8 located in E1. E8^E2 is a DNA-binding protein, which has been characterized as a transcriptional and viral replication repressor through the recruitment of the NCoR platform (Dreer et al., 2016). This isoform participates in the maintenance of the episomal state of the virus and limits viral replication.
E4: E4 localizes to the central region of E2 and is synthesized as a fusion protein of the E1^E4 transcript from promoter 8. E4 is the most abundant protein in the infected epithelium. This protein participates in amplification, release, and transmission of the viral genome. It also promotes the viral productive phase associated with epithelial differentiation-dependent post-translational modifications and participates in cytokeratins collapse and in G2 cell cycle arrest. It has been shown that E4 expression decreases as CIN grades increase until it disappears. E4 also correlates with low levels of methylation in the CADM1, MAL, or miR124-2 genes (Zummeren et al., 2018). These genes are implicated in cervical carcinogenesis and are discussed later in this review. E4 also participates in chromatin regulation by preventing the interaction of MCM2 and MCM7 with chromatin (Wilson et al., 2005; Doorbar 2013).
E5: E5 is mainly localized in the intracellular membranes of the endoplasmic reticulum and Golgi apparatus. In these compartments, it binds with vacuolar ATPase to decrease endosomal acidification, which leads to the degradation of cell surface receptors. Studies have shown that E5 plays an essential role in carcinogenesis, as it is critical in cell transformation and immune modulation by interacting with the major histocompatibility complex class I (MHC-class I). In addition, this protein is involved in several biological processes such as differentiation-dependent viral cycle, DNA synthesis in supra-basal layers, viral amplification, and maintenance of the episomal state, and in the expression of late genes (DiMaio and Petti 2013). E5 activates EGFR in a ligand-independent manner, thus affecting signaling pathways such as MAPK, PI3K-AKT. E5 expression is associated with increased cyclooxygenase 2 (COX2) (Ilahi and Bhatti 2020).
E6 y E7: One of the best described activities by HPV-HR proteins in the transformation process is the inactivation of the tumor suppressor genes p53 and pRB by E6 and E7, respectively (Vats et al., 2021). However, many cellular proteins are affected by direct or indirect interaction with these viral oncoproteins.
E6 is epigenetically regulated by CTCF and can generate distinct isoforms by splicing. Splicing products (E6*) antagonize the effect of full-length E6. E6 participates in immortalization, evasion of apoptosis, evasion of the immune response, the inhibition of the interferon response, participates in epigenetic reprogramming, and regulates chromatin structure through interaction with CBP and P300 (Vande Pol and Klingelhutz 2013; Estevao et al., 2019; Mac and Moody 2020; Scarth et al., 2021; Vats et al., 2021).
E7 is epigenetically regulated by CTCF. It is considered the main protein with transforming activity. It participates in the reprogramming of the cellular environment to facilitate viral replication and reduces the expression of MHC I in conjunction with E5. Like E6, E7 interacts with multiple proteins to regulate epigenetic mechanisms. E7 has been shown to interact with and induce the activity of DNMT1, DNMT3a, EZH2, CBP, and P300. In addition, E7 modulates viral replication by interacting with HDAC1 and HDAC2 via Mi2β. E7 induces KDM6A and KDM6B expression, leading to host cell reprogramming by re-activation of HOX and p16INK4A genes. E7 also post-transcriptionally increases SETD2 levels and regulates H3K36me3 marks on viral chromatin (Burgers et al., 2007; Roman and Munger 2013; Saha et al., 2017; Gautam et al., 2018; Pal and Kundu 2019; Mac and Moody 2020; Vats et al., 2021).
Studies show that enhancer of zeste homolog 2 (EZH2) is specifically activated by E7 from HPV16 through the release of E2F. EZH2 is an enzymatic subunit that mediates H3K27me3 and contributes to cell proliferation (G1-S) involving several mechanisms such as the stimulation of Cyclin D1 and E, as well as retinoic acid signaling or the WNT pathway. In addition, EZH2 can suppress apoptosis through the silencing of the miRNA-31 promoter, which has among its regulatory targets the antiapoptotic gene E2F6 (Holland et al., 2008; Zhang et al., 2014).
Despite EZH2 increase by E6 and E7, the H3K27me3 mark does not increase; this is associated with the Ser21 phosphorylation, which is mediated by AKT. This modification prevents the methylation enzymatic activity while it increases the expression of histone demethylase KDM6A, associated with activation of HOX genes (Hyland et al., 2011).
Late ORF: L1 is the major capsid protein, which, can self-assemble into virus-like particles (VLP) and is mainly involved in the cell infection process. The L1 region of HPV is susceptible to DNA methylation which, interestingly, has been associated with lesion grade and viral integration. For example, methylation of CpG 5698 and 5617 is associated with high-grade lesions and CC (Turan et al., 2006; Buck et al., 2013; Wang and Roden 2013; Torres-Rojas et al., 2018).
LCR: Regarding the LCR region, it is a viral sequence that does not encode proteins, but it is one of the most important regions of the virus since it contains the origin of replication and interacts with multiple epigenetic regulators and transcription factors. Comprising approximately 10% of the genome (Bernard 2013).
3 Changes in the Human Papillomavirus Epigenome
Human papillomavirus DNA associates with proteins similar to H2A, H2B, H3, and H4 histones. Electron microscopy analyses revealed that the complete HPV genome contains 32 nucleosomes, each measuring about 12 nm in diameter (Stunkel and Bernard 1999; Burley et al., 2020). Interestingly, HPV chromatin is also epigenetically regulated through various histone post-translational modifications (acetylation, ubiquitination, methylation, and phosphorylation among others) and viral DNA methylation (Mac and Moody 2020).
Changes in chromatin are closely linked to the epithelial differentiation process. Viral promoters contact each other by three-dimensional folding of nucleosomes at different stages of the viral cycle. Some authors have proposed that chromatin modifications occur mainly in the proximal upstream regulatory region of E6 in undifferentiated cells, whereas in differentiated cells, it occurs in the late promoter region (del Mar Pena and Laimins 2001). However, some studies show that late promoters are not fully repressed in basal cells (Wooldridge and Laimins 2008).
Initially, when HPV infects undifferentiated cells of the basal layer of the epithelium, the viruses remain in an epigenetically repressed episomal state, with enrichment of H3K27Me3 and low levels of H3K4Me3. Consequently, there is a low viral load (between 50 and 100 viral copies per cell), which favors a deficient activation of the immune response (Ferguson et al., 2021).
Chromatin immunoprecipitation assays on HPV31 revealed transcriptionally active chromatin marks such as acetylated H3K4me2 in HPV promoter regions p97 (early in LCR) and p742 (late in E7). Interestingly, these epigenetic marks increase after differentiation triggering E6/E7 transcription. In contrast, undifferentiated cells carry repressive marks such as H3K9me2 and H3K27me2 (Wooldridge and Laimins 2008).
Notably, the p97 promoter is regulated from the 5′ end in the LCR region; therein lies the origin of replication, the keratinocyte early enhancer (KE) and an auxiliary enhancer (AE) that enhances KE activation (Kanaya et al., 1997).
Transcriptional activation marks recruit factors such as C/EBP-β and C/EBP-α to the LCR and late promoter sequences. Notably, the KE region is also an enhancer of the late promoter, so it has been proposed that differentiation alone might be enough to mildly induce late transcription without active viral replication (Wooldridge and Laimins 2008) (Figure 2).
FIGURE 2. HPV Chromatin. The figure illustrates the HPV chromatin regulation. The inner-circle indicates the distribution of the viral genes. The tags indicate the position of the corresponding chromatin mark. In addition, the characteristic loop between CTCF and the LCR region is illustrated. See text for more detail.
Some studies have shown that negative regulation of the HPV-16/18 E6-E7 promoter requires the binding of YY1 to the promoter KE sequence (Lace et al., 2009). YY1-directed epigenetic silencing involves the recruitment of Polycomb repressor complexes 1 and 2 (PRC1 and PRC2) to viral chromatin leading to the enrichment of H3K27me3 and H2AK199Ub marks. Interestingly, epithelial differentiation inhibits YY1 expression, resulting in positive regulation of E6 and E7 (Mac and Moody 2020). Notably, YY1 regulates HPV 16, 18, 11, and 8 promoter sequences; however, the mechanism of action has not been fully elucidated (Lace et al., 2009).
HPV infection is initially established in undifferentiated basal cells, so there is no expression of E6 and E7. This transcriptional state is also dependent on CCCTC binding factor (CTCF). Studies in HPV 18 have described that CTCF can interact with the E2 ORF and, in turn, with YY1, whereby CTCF-YY1 mediates a three-dimensional interaction between the E2 ORF and the LCR region to establish a repressive chromatin loop (Pentland et al., 2018). Epithelial differentiation disrupts the CTCF-dependent chromatin loop (Figure 2). Studies in HPV18 showed that abrogation of CTCF binding to the E2 ORF increases the transcriptional activity of the p102 early promoter, increasing the expression of E6/E7 proteins (Ferguson et al., 2021).
Recent findings have shown that viral integration events frequently trigger host chromatin changes, which enhance the carcinogenic process in HPV+ tumors. Interestingly, CTCF is a key regulator in these processes, since mutations that disrupt the binding sites as well as the introduction of new sites for CTCF can reorganize chromatin and modify transcriptional regulation, allowing tumorigenesis beyond the transcriptional regulation of E6/E7 (Karimzadeh et al., 2022). In addition, sequencing studies elucidated that CTCF may play an important role in the regulation of alternative splicing of viral genes (Ferguson et al., 2021).
Studies using C33a cells infected with HPV16 have revealed that the histone deacetylase sirtuin 1 (SIRT1) is recruited to the HPV replication origin to become part of the functional E1-E2 complex. Loss-of-function assays of SIRT1 by CRISPR-Cas9 showed increased replication of E1-E2 due to a growing stability of E2 caused by an increase in acetylation (Das et al., 2017). Interestingly, Dipon Das et al. showed that SIRT1 post-transcriptionally regulates the DNA repair protein Werner helicase (WRN) recruited to the E1-E2 complex to regulate replication fidelity. Like SIRT1, deletion of WRN also results in uncontrolled viral replication (Das et al., 2019) (Figure 2). By chromatin immunoprecipitation (ChIP) assays, SIRT1 was found to bind to multiple regions of the HPV31 genome except for the L1 region. Particularly in undifferentiated cells, SIRT1 binds to the LCR, erasing the H1K26Ac and H4K16Ac acetylation marks, thus repressing HPV31 late gene expression (Langsfeld et al., 2015) (Figure 2).
Interestingly, HPV31 E7 has been reported to increase SETD2 protein levels. SETD2 is a histone methyltransferase that writes the H3K36me3 mark (active transcription) on viral chromatin. SETD2 can collaborate in the viral cycle regulation by recruiting of H3K36me3 readers to HPV. Furthermore, experiments in CIN612 cells showed that SETD2 writes the H3K36me3 mark abundantly in E2, E4, and E5 ORFs and not in E6 and E7 (Figure 2) (Gautam et al., 2018).
Likewise, the nuclear receptors have an essential role in viral chromatin regulation since they allow the recruitment of chromatin remodelers such as BRG1 and INI1 (Figure 2). These proteins relate to SWI/SNF complex that can recruit the MLL1 and SETD1A enzymes to write the H3K4me1/3 marks, these complexes can also recruit HAT proteins such as p300 (Groves et al., 2016).
3.1 Viral DNA Methylation
The first study to analyze viral DNA methylation, reported changes in the LCR with potential gene regulatory implications (Burnett and Sleeman 1984). Subsequent studies determined that these methylation changes are associated with the viral life cycle and the cell differentiation. The upstream regions of LCR (3′of L1) and E6 have been of particular interest in some methylation studies (Johannsen and Lambert 2013). Vinokurova and von Knebel Doeberitz (2011) determined CpG methylation changes on the LCR of HPV16 isolated from epithelia by microdissection. In epithelium lacking cytopathic alterations, CpGs remain methylated, whereas, in low-grade lesions or permissive infections, changes in methylation are related to the epithelial differentiation layer, keeping the E2-binding site (E2BS) of the promoter region demethylated. In contrast, E2BS1 is methylated in high-grade lesions and is associated with early promoter activation (P97) and viral integration (Vinokurova and von Knebel Doeberitz 2011). Regarding invasive carcinoma, Chaiwongkot et al. (2013) identified differences in methylation as a function of viral integration status, specifically at the E2SB1 (7453 and 7459), E2SB3 (37, 43) and E2SB4 (52, 58) sites of HPV16 LCR; methylation being lower in carcinomas with the integrated virus than in episomal carcinomas (Figure 3). In addition, viral genome copy number appears to be related to methylation status, as higher methylation levels were identified at LCR in cases with higher integrated copy number (Chaiwongkot et al., 2013).
FIGURE 3. Changes of methylation in LCR of HPV 16. (A) LCR region binding sites. (B) Schematic illustrating LCR methylation changes of HPV 16 isolated from different epithelial layers. The letters “S,” “I,” and “B” refer to superficial, intermediate, and basal layers respectively, and at different stages of viral infection and transformation. Dark spots show methylated CpG clear spots show unmethylated CpG. Data from Vinokurova and von Knebel Doeberitz 2011, Chaiwongkot et al. (2013).
Piyathilake et al. (2014) reported an association between low risk CIN2+, and high methylation levels in the HPV 16 promoter and enhancer regions. Additionally, women with high folate concentrations and high methylation have a lower risk of being diagnosed with CIN2+. A similar effect was observed with vitamin B12, suggesting that folate and vitamin B12 intake may reduce HPV hypomethylation associated with the risk of developing CIN2. That is, two of the most important methyl donor micronutrients are associated with site-specific methylation of the promoter (position 31, 37, 43, 43, 52, and 58) and enhancer (7862) regions of the E6 gene, which is associated with a low probability of CIN2+ diagnosis (Piyathilake et al., 2014) (Figure 3).
Several authors have reported increased methylation in the L1 region associated with CIN3 or CC. Recently Bee et al. identified that methylation of the CpG 5615 site of L1 of HPV52 is most frequently associated with NIC3, and SNPs associated with the C clade (Bee et al., 2021). Additionally, increased DNA methylation has been observed in the L1, L2, E2, and E5 regions in high-grade lesions. As the degree of lesion progresses grade advances, increased methylation at L1 is more evident in HPV16 than in HPV18 (Torres-Rojas et al., 2018). Concerning viral integration status, it is more frequent in HPV18 cases than in HPV16 (Torres-Rojas et al., 2018). These data indicate that the viral methylation/integration profile may be related to genotype and viral species. In this regard, in cervical cancer the 3′ LCR methylation and E1/E2 region integrity are higher in HPV16 than in HPV18 and HPV45 (Amaro-Filho et al., 2018). It has been shown that viral integration is dispensable for malignant transformation and the alternative mechanism for E6 and E7 deregulation could be DNA methylation at the E2SB sites in HPV16 (Cheung et al., 2013).
4 Epigenetic Alterations in Human Papillomavirus-Associated Neoplasms
Diverse evidence has shown that patients with HNSCC may have different responses to treatment depending on HPV status, mainly due to transcriptomic changes involving genes related direct or indirectly by viral presence (Pyeon et al., 2007; Lohavanichbutr et al., 2009; Mendez-Matias et al., 2021).
The characterization of epigenetic alterations in neoplastic transformation is a significant event. Unlike genetic alterations, epigenetic changes are potentially reversible. Although the mechanisms that lead to epigenetic alterations are not fully understood, it is relevant that they could be used as diagnostic and prognostic biomarkers as well as therapeutic targets. This section will focus on some of the epigenetic alterations in precursor lesions and invasive cancer and their potential use as biomarkers.
Sartor et al. (2011), identified by global methylation analysis that HPV-positive cell lines have higher levels of methylation in LINE-1 genic regions compared to negative ones, 77% and 45%, respectively. It was reported that the 19q13 region has high methylation rates; as well as some target genes of the polycomb complex (CCNA1, PREX1, RUNX2, and SPON2). Several differentially methylated regions (DMR) have been found in tumors, as well as in genes of great biological importance such as ESR1, and DCC (important in HNSCC) (Sartor et al., 2011).
Lechner et al. (2013), identified an HPV + HNSCC-associated methylation signature, regardless of anatomic site (Lechner et al., 2013). Additionally, two HPV+ subgroups were identified, associated with high levels of methylation and poor prognosis. Other studies have also identified two subgroups within HPV+ based on expression and viral integration status (Ren et al., 2020).
In contrast, Brennan et al. characterized five subtypes of HNSCC based on methylation profiles, one of them associated with HPV that corresponds to the previously described atypical expression profile and is associated with better survival (Brennan et al., 2017). Nakawaga et al. (2020), identified four subtypes of HNSCC; two of them are associated with HPV, another presents high levels of methylation in promoter regions and correlates with a better prognosis, and the last one has intermediate levels of methylation (Nakagawa et al., 2020).
Ren et al. (2018), identified a set of 17 DMRs in OPSCC-HPV+, however, only seven regions are capable of discriminating HPV+ OPSCC from healthy tissue (KCNA3, EMBP1, CCDC181, ELMO1, C1QL3, MIR129–2, and ZNF137P) (Ren et al., 2018).
Misawa et al. (2020), identified that methylation of 10 genes is associated with recurrence in OPSCC HPV+, this signature includes ATP2A1, MRGPRF, GNMT, GPT, LYNX1, MAL, MGC16275, CALML5, DNAJC5G, and LY6D, being the last three from circulating tumor DNA (Misawa et al., 2020).
Recently, Berglund et al. (2022), identified DMRs in HNSCC HPV16; the most significant were: SYCP2, MSX2, HLTF, PITX2 and GRAMD4. Interestingly, the methylation profiles with genotypes of the A7 species such as HPV 31, 33 and 35 are compatible with those presented in HPV16, and a different methylation profile was found in cases with genotypes belonging to other viral species (Berglund et al., 2022).
Comparing with normal tissue, 86 DMRs were identified in anal carcinoma and 36 in CC, of which 17 DMRs are shared in both malignancies (ASCL1, ATP10A, CCDC81, DPP10, FMN2, MARCH11, MIR129-2, PAX1, PEX5L, RYR2, SORCS3, T, WDR17, ZIK1, ZNF154, ZNF177, ZNF529/ZNF382) (Siegel et al., 2021).
Alternatively, different DMRs and differentially hydroxymethylated regions (DhMR) in genes involved in the Hippo pathway were identified in NIC 3 and CC. Interestingly, the number of DMR increases as the degree of injury progresses; however, the DhMRs decrease from NIC 3 to CC.
Some genes from the DMRs were associated with survival (DES, MAL, MTIF2, PIP5K1A, RPS6KA6, ANGEL2, MPP, and PAPSS2) (Han et al., 2021). Gagliardi et al. (2020) identified a greater number of DMRs in cervical tumors with genotypes of the A9 vs. A7 species. In addition, tumors with A7 genotypes have higher H3K4me1 (associated with enhancers) and H3K27Ac (associated with active promoters), while tumors infected with genotypes of the A9 species have an increase in H3K4me3 (associated with active promoters). This could have implications in the expression profiles and tumor aggressiveness (Gagliardi et al., 2020).
Globally, a greater number of methylated regions have been reported in HPV-positive cell lines and malignancies than in negative ones (Sartor et al., 2011; Ren et al., 2018). However, in HPV-positive squamous cell carcinoma of the penis (PeSCC), it was reported that 77% of the variable methylation positions (MVPs) are hypomethylated (Feber et al., 2015). Among the 960 MVPs associated with HPV, there are WNT pathway genes such as GRAMD4 and GPX5.
It is noteworthy that 30 MVPs were validated in independent cohorts in HNSCC and CC, managing to stratify the HPV+ cases from the negative ones according to epigenetic signature; additionally, significant differences were shown in survival at 5 years, in HNSCC (81% vs. 38%) and in CCSCC (77% and 50%), which was a better predictor of survival than the sole presence of HPV (Feber et al., 2015).
In other PeSCC studies, 65 differentially methylated genes were identified, with inverse correlation with expression in HPV+ cases. Among which are CD70, HN1, FZD5, FSCN1, and PRR16 (Kuasne et al., 2015).
The 5-hydroxymethylation of the cytokine (5hmC) associated with loss of transcriptional repression has been identified in lower levels in HPV + HNSCC compared to HPV- cases. Among the genes associated with hyper-5hmC in HPV+ tumors are CDKN2A and enrichment of pathway genes such as desmosomes, NF-kB, oxidoreductase activity, and mesenchymal cell differentiation where genes such as SNAI2, BMP2, SMAD2 and TGFB2 are found (Liu et al., 2020).
Many cellular genes have been proposed as diagnostic and prognostic biomarkers in CC and its precursor lesions (Hesselink et al., 2014; Kalantari et al., 2014; Verhoef et al., 2014; Boers et al., 2016). Global methylation profiles analysis in 12 cancer types identified 4 CC-specific markers cg07211381 (RAB3C), cg12205729 (GABRA2), cg20708961 (ZNF257), and cg26490054 (SLC5A8) with 96.2% sensitivity and 95.2% specificity. Additionally, this methylation profile can distinguish between tumors and normal tissues (Xu W et al., 2019).
Nowadays, the methylome offers several candidates with potential as clinical biomarkers and therapeutic targets in HPV+ neoplasms. In addition, non-coding RNAs have been highlighted in the scientific and therapeutic landscape. Interestingly, some ncRNAs are transcribed by the own HPV.
5 The Pathogenetic Role of ncRNAs in HPV+ Cancers
The invention of new sequencing technologies, bioinformatics tools, as well as new experimental strategies, have made possible the discovery of new functional molecules of non-coding RNA (ncRNA) (Panda 2018; Huang et al., 2019). It is well known that HPV can transcribe non-coding RNAs such as circular RNAs (circRNAs) and microRNAs (miRNAs). Furthermore, it has been suggested that the altered expression of some long-noncoding RNAs (lncRNAs) in some types of cancer is a direct consequence of the expression of the viral oncoproteins E6 or E7. Therefore, we performed a bibliographic search for ncRNAs encoded directly by HPV.
In this section, we begin describing the HPV non-coding transcriptome, and finally, we highlight the function of lncRNAs induced by the viral oncoproteins in the host cell.
5.1 Human Papillomavirus Non-coding Transcriptome
Recent studies have shown that within the HPV genome, there are regions that contain circRNA and miRNA genes (Figure 4).
FIGURE 4. HPV-encoded miRNAs. (A) Illustration shows the miRNAs coding site inside the genome of different HPV subtypes. The predicted biological processes for some HPV16 miRNAs are depicted in the figure. The circE7 synthesis site is also shown. (B) The scheme represents the potential interaction sites of miRNAs in the viral genome of HPV-16 and HPV-41. The predicted target viral sequences location for each miRNA are shown in the same color as the corresponding miRNA sequence (blue, red, or purple).
CircRNAs are single stranded, covalently closed, circular RNA molecules. CircRNAs can be generated from several genomic regions (intergenic, intronic, and coding). Regarding the biosynthesis of circRNAs, two biosynthesis models have been proposed, and both involve back-splicing (Miao et al., 2021).
Although they have been little studied, there is evidence showing that they perform essential cellular functions. These molecules can act as miRNA sponges, or in some cases, they might indeed have coding potential, as is the case of circE7 (Panda 2018). HPV16 cervical carcinoma cells have high expression levels of circE7 (Zhao et al., 2019) (Figure 4A), and this circRNA is primarily localized in the cytoplasm and requires post-transcriptional modification of N6-methyladenosine (m6A) in the UTR to associate with polysomes and enhance E7 protein levels thus promoting cancer cell growth. These new findings show how these types of nucleic acids transcribed from HPV aid in malignant transformation and increase oncogenic potential in cells.
There are few published studies regarding the miRNAs transcribed from the HPV genome. However, it has been shown that viral miRNAs modulate cellular and viral gene expression, promoting carcinogenesis and evading the host’s immune system. Shun-Long Weng et al. (2018) performed bioinformatic analyzes to identify HPV-encoded miRNAs. They analyzed genomes from different high-risk HPV genotypes to identify viral pre-miRNAs based on phylogenetic information and structural evolution. Notably, HPV16 showed coding potential for HPV16-miR-1, HPV16-miR-2, and HPV16-miR-3, which are located inside the ORFs of E6, E1, and L2, respectively (Weng et al., 2018) (Figure 4A). The authors also identified possible targets of these miRNAs; however, these findings require experimental validations.
The predicted targets of HPV16-miR-1 (ARID5B, ZEB2, THBS1, and STAT5B) are mainly involved with biological processes such as cell migration, cell motility, response to hypoxia, and regulation of cell adhesion. For HPV16-miR-3, the predicted targets (SYNE1, GATA6, GULP1, PDE1B, IGFBP4, PRELP, and MYH1) are associated with cell death, development, and cell differentiation (Figure 4A). Regarding HPV16-miR-2, the predicted targets are AFF3, FRMD7, IGDCC4, MYRIP, NRN1, PMP22, and RBPMS; however, these genes do not associate with cervical cancer progression (Weng et al., 2018). Strikingly, some of these miRNAs could be targets for new antiviral drugs, and miRNAs specific to certain HPV subtypes could serve as possible biomarkers. In addition, the authors found putative miRNAs in other HPV genotypes such as HPV-18 (1 miRNA), HPV-31 (3 miRNAs), HPV-45 (2 miRNAs), and HPV-52 (1 miRNA) (Weng et al., 2018) (Figure 4A).
Through sequencing analysis, Qian et al. (2013) discovered nine HPV-transcribed miRNAs. These miRNAs were successfully validated by qPCR and in situ hybridization in cervical tissue samples and cervical cell lines (Qian et al., 2013). Two of the validated miRNAs were encoded by HPV 16 (HPV16-mir-H1 and HPV16-mir-H2), one by HPV 38 (HPV38-mir-H1), and one by HPV 68 (HPV68-mir-H1) (Figure 4A).
The authors also found both human and viral target genes of miRNAs encoded by HPV16. HPV16-miR H1 and miR H2 have unique targets in the human genome, and 15 shared target genes (CDC2L6, EIF2C1, IMPAD1, BNC2, SNX27, TNRC6B, BACH2, CYP26B1, DDX19B, FGF7, PBRM1, PHACTR2, RBM3, RGS7BP, TEAD1). Regarding the HPV16 viral genome, authors found four target sequences for HPV16-miR-H1-1, two in the E5 region and two in the L1 gene (Qian et al., 2013) (Figure 4B).
HPV16-miR-H2-1 presented two targets in the viral genome, one in the LCR region and the other in the L1 gene (Figure 4B). A gene ontology analysis suggests that HPV16-miR-H1-1 is involved in processes related to the M phase of the cell cycle, immune system development, adhesion, and migration. These processes are crucial for viral infection, activation of immune cells, and other mechanisms crucial for neoplastic development (Qian et al., 2013; Gallo et al., 2020).
Subsequent studies conducted by Elina Virtanen et al. (2016) evaluated the expression of some miRNAs encoded by HPV 16 reported by Qian et al. (2013). The authors extracted RNA from HPV 16 samples obtained from paraffin-embedded cervical tissue, cell lines, liquid cytology, and cervical cells obtained from colposcopy. Most HPV-encoded microRNAs were detected in all samples; however, HPV16-miR-H2 was not detected in these assays (Virtanen et al., 2016).
Bioinformatic analyzes and multiple experimental strategies performed by Chirayil et al. (2018) showed the potential of various HPV genotypes to transcribe miRNAs. The most significant miRNAs are transcribed from HPV41 (hpv41-miR-H1), HPV17 (hpv17-miR-H1), HPV37 (pv37-miR-H1) and even an HPV from the Fringilla coelebs FcPV1 (fcpv1-miR-F1 and fcpv1-miR-F2) (Figures 4A,B). Interestingly, the authors suggested that FcPV and HPV41 miRNAs may regulate viral gene expression as they found potential target sequences for these miRNAs in HPV (Chirayil et al., 2018).
The research field on miRNAs and cirRNAs encoded from HPV encourages new therapeutic strategies based on RNA for various medical applications. However, it would be essential to advance in greater depth in knowing the biological functions of these new molecules.
5.2 Long Non-Coding RNAs Regulated by E6/E7 Oncoproteins
LncRNAs are transcripts of more than 200 nucleotides long, which lack coding potential. LncRNAs have been classified in several ways, for example, according to their position concerning neighboring encoding genes (divergent, convergent, intergenic, overlapping—sense/antisense, intronic, among others) (Schmitz et al., 2016) or according to their cis or trans transcription site (Tong and Yin 2021). The diversity of interactions that lncRNAs have with cellular components (RNA-DNA, RNA- Protein, RNA-RNA), have function to classified lncRNAs according to their function as decoys, scaffolds, enhancers, guides, sponges, or competing endogenous (Guttman and Rinn 2012; Castro-Oropeza et al., 2018).
LncRNAs acting as Decoys, induce or prevent proper regulation of gene expression through association with DNA-binding proteins. LncRNA scaffolds assemble various structural complexes to finally adapt them to specific transcriptional regulatory regions (St Laurent et al., 2015). LncRNAs with guide function associates with proteins to direct them to specific regions of the genome and thus cause an effect at the transcriptional or chromatin level. LncRNA enhancers can organize the nuclear architecture to form enhancer contacts thus promoting gene expression. Sponge refers to lncRNAs containing responsive elements to miRNAs, to sequester them and thus prevent their function on their target mRNAs. On the other hand, lncRNAs can act as competing endogenous RNAs, and bind to an mRNA to give it stability by recruiting proteins, thus preventing degradation (Rinn and Chang 2012). Generally, the function of lncRNAs is associated with their subcellular localization, some may regulate specific functions in the nucleus, while others are exported to the cytoplasm, where they exert a different function. The lncRNAs located in the nucleus can modify chromatin by binding epigenetic factors to DNA to regulate transcription and RNA processing. Gene expression regulation is an intricated process that involves epigenetic modifications such as DNA methylation, histone modifications, chromatin remodeling, and ncRNA-based mechanisms. Although the role of nuclear lncRNAs in epigenetic regulation has been extensively investigated, there is little information about the nuclear lncRNAs involved in the epigenetic regulation of HPV-associated neoplasms.
Studies have shown that the presence of HPV in cancerous tissues stimulates the differential expression of lncRNAs. For example, Ma et al. found 107 differentially expressed lncRNAs when comparing normal tissue vs. HPV-negative HNSCC tumors. On the other hand, in the comparision of normal tissue vs. HPV-positive HNSCC, the number of lncRNAs with differential expression increased to 801 transcripts (Ma et al., 2017). Furthermore, studies by Yang et al. compared the transcriptome of HPV-positive cervix cancer cell lines (SiHa and HeLa) vs. an HPV-negative cell line (C33). Interestingly, the authors found changes in the expression of lncRNAs and mRNAs, suggesting that the altered transcripts could be related to HPV-induced proliferation (Yang et al., 2016). These data suggest that HPV involvement is significant in altering lncRNAs in cells.
It should be noted that some studies have proposed lncRNAs expression signatures in cancer associated with HPV infection. However, the molecular mechanisms of action of these molecules have not yet been elucidated. Also, it is still unknown if these expression changes are an indirect consequence of viral infection or if ncRNAs and HPV oncoproteins directly regulate them.
Viral proteins can alter the biological activities of lncRNAs either by binding to them directly or indirectly, thus changing their ability to interact with proteins or other nucleic acids (Sharma and Munger 2020c). HPV, like other viruses, is capable of dramatically altering the expression profile of host cell lncRNAs. Several investigations have identified numerous lncRNAs whose expression is modulated by HPV infection primarily through E6 and E7 proteins (Sharma and Munger 2020c).
Interestingly, lncRNAs have been found to enhance virus invasion and host cell reprogramming to allow sustained infection. Conversely, there are also antiviral lncRNAs that regulate the innate and adaptive immune response to eliminate viral particles. In the following section, we will detail some oncogenic and tumor suppressor lncRNAs that are regulated by E6 or E7 and that influence some features leading to malignancy (Figure 5) (Ouyang et al., 2016).
FIGURE 5. lncRNAs regulated by E6/E7 oncoproteins. The figure shows the molecular mechanisms of lncRNAs regulated by the E6/E7 oncoproteins. Blue circles refer to oncogenic lncRNAs and red circles to tumor suppressors.
5.2.1 Oncogenic Long Non-Coding RNAs
lncRNA SNHG8 -A recent study revealed that lncRNA SNHG8 is highly expressed in HPV-induced cervical cancer cells compared to HPV-negative cells (Qu et al., 2020). Notably, this lncRNA has molecular functions in other types of cancer. For example, it acts as a miR-634 sponge in breast cancer, enabling the expression of the ZBTB20 oncogene, thus triggering breast cancer progression (Xu et al., 2021). Also, in hepatocellular carcinoma, the lncRNA SNHG8 acts as a sponge for the tumor suppressor miR-1493 (Dong et al., 2018). Interestingly, in cervical cancer (HPV+), the lncRNA SNHG8 has an epigenetic mechanism. It is mainly located in the nucleus and interacts with EZH2 to promote H3K27me3 enrichment on the RECK (Reversion inducing cysteine-rich protein with Kazal motifs) promoter, thus suppressing RECK expression. Therefore, the functional implications of lncRNA SNHG8 expression include proliferation, migration, and apoptosis inhibition in cervical cancer cells (Figure 5) (Qu et al., 2020).
lncRNA CCEPR -The Cervical Carcinoma Expressed PCNA Regulatory LncRNA (CCEPR) is highly expressed in cervical cancer (Sharma and Munger 2018). For this reason, CCEPR is also known as cervical carcinoma high-expressed long non-coding RNA 1 (CCHE1) (Chen et al., 2017). Also, its high expression levels correlate with tumor size and poor survival prognosis in patients with this type of cancer. In a keratinocyte model that CCEPR is regulated by HPV16 E6 through a p53-independent mechanism (Sharma and Munger 2018). Interestingly, the viral E6 protein expression increased CCEPR, similar to HPV16-positive cervical cancer cells. In addition, the expression of both E6 and E7 further potentiates the expression of this lncRNA.
CCEPR localizes in the nucleus, where it probably interacts with proteins and nucleic acids rather than modulating PCNA levels in the cytoplasm (Sharma and Munger 2018). Notably, in oral squamous cell carcinoma (OSCC), CCEPR localizes to the nucleus and cytoplasm and acts as an oncogene to promote OSCC progression by increasing PAK2 expression through competitively binding miR-922 in OSCC cells (Figure 5) (Wang et al., 2020). Therefore, we consider it interesting to evaluate the presence of HPV in patients with OSCC to find out if this expression could also be related to the presence of the viral oncoproteins.
LncRNA-FAM83H-AS1 -lncRNA FAM83H-AS1 is overexpressed in cervical cancer cells and HPV16-positive head and neck cancer cell lines compared to HPV-16-negative cells (Barr et al., 2019). Notably, the presence of HPV16 correlates with elevated levels of FAM38H-AS1 from the early stages of carcinogenesis. FAM38H-AS1 is regulated by HPV through the presence of p300 in a p53-independent mechanism (Figure 5). Downregulation of this lncRNA impacts on proliferation, migration, and apoptosis (Barr et al., 2019).
LncRNA-FANCI-2 -It has been described that patients with aberrant expression of Fanconi anemia-linked genes have a high incidence of HPV-high risk-induced neoplasms. In a study by Haibin Liu et al. (2021) the authors showed that the expression of lnc-FANCI-2 and FANCI are positively co-expressed in HPV-18-infected cervical cancer cells (Liu et al., 2021). A cis-regulation of these genes has been proposed; however, the precise mechanism remains unknown. E6 and E7 oncoproteins inhibit miR-29-a expression, culminating in an increment in the expression of the transcription factor YY1. Interestingly, in CaSki cells, YY1 interacts with E7 to bind to the lnc-FANCI-2 promoter and induces transcriptional activation (Figure 3) (Liu et al., 2021).
LncRNA-EBIC -In cervical cancer cells, the expression of E6 and E7 stimulates the expression of EZH2-binding lncRNA (lncRNA-EBIC), also known as TMPOP2 (Thymopoietin pseudogene 2) (Sun et al., 2014; He et al., 2019). Mechanistically, EBIC interacts in the cell nucleus with EZH2, a central component of the PCR2 complex that catalyzes the H3K27me3 histone modification (Sun et al., 2014). Through this mechanism, lnc-EBIC exerts oncogenic properties as it participates in the silencing of E-cadherin to promote cell invasion and metastasis in cervical cancer.
The function of lnc-EBIC in the presence of HPV viral oncoproteins has been recently investigated. Hongpeng He et al. (2019) found that EBIC acts in the cytoplasm as an oncogene by inhibiting the function of some tumor suppressor miRNAs, i.e., it acts as an endogenous competitor of miR-375 and miR-139, whose function is to degrade E6 and E7 mRNAs (He et al., 2019). In this fashion, lnc-EBIC contributes to the tumoral phenotype of cervical cancer cells (Figure 5).
Regarding the transcriptional regulation of lnc-EBIC, p53-responsive elements were found in its promoter region. ChIP studies have shown a decreased enrichment of p53 in the promoter region of this lncRNA in some types of cancer (He et al., 2019). In addition, research by Jun Wang et al. has shown that exogenous expression of E7 from HPV16/18 significantly promotes lnc-EBIC expression, causing increased cell proliferation, migration, invasion, and inhibition of apoptosis through the TAL1/lnc-EBIC/KLHDC7B signaling axis (Wang et al., 2021).
LINC01101 y LINC00277 -Iancu et al. (2017) performed microarray assays based on E7 inhibition in HeLa (Iancu et al., 2017). The authors found many differentially expressed lncRNAs, including LINC01101 and LINC00277 (Figure 5), which were positively upregulated in the cells due to the lack of E7. In addition, they showed that high-risk HPV-positive samples have lower expression of LINC01101 and LINC00277 (Iancu et al., 2017). Moreover, low expression of LINC01101 correlated with FIGO stage and lymph node metastasis (Iancu et al., 2017).
The function and mechanism of action of LINC00277 (known as EWSAT1) appear to be tissue specific. For example, in Ewing’s sarcoma, EWSAT1 is positively regulated by EWS-FLI1 in primary human pediatric mesenchymal progenitor cells, and its high expression correlates with increased cell proliferation (Marques Howarth et al., 2014). Furthermore, Peng Song and Shu-Cheng Yin have shown that EWSAT1 is overexpressed in human nasopharyngeal carcinoma cells and correlates with poor survival (Song and Yin 2016). Mechanistically, it has been revealed that EWSAT1 indirectly regulates cyclin D1 expression, as it functions as a miRNA sponge for miR-326 and miR330-5p, which target cyclin D1 for degradation (Song and Yin 2016). Studies in osteosarcoma have shown that EWSAT1 is related to cellular phenomena such as proliferation, migration, invasion, and metastasis. As for its mechanism of action, EWSAT1 regulates the expression of the lncRNA MEG3 (Sun et al., 2016).
CRNDE -Studies by Harden et al. (2017), have suggested that the lncRNAs CRNDE (Colorectal Neoplasia Differentially Expressed), DANCR and TINCR, are transcribed due to HPV16 E6/E7 expression and could contribute to cancer development. For example, high expression of DANCR and low expression of TINCR could overlap with processes involved in epithelial differentiation (Harden et al., 2017; Sharma and Munger 2020c).
We believe that it would be essential to delve deeper into the regulation of the expression of these lncRNAs and their association with HPV, as they exert a very prominent role in cancer regulation. For example, lnc-CRNDE has multiple functions in the cell. For instance, it could function as ceRNA, RNA-decoy, scaffold and miRNA sponge in various types of cancer (Lu et al., 2020). Even lnc-CRNDE has been highlighted for its role in epigenetic regulation because it interacts with EZH2 to repress DUSP5 and CDKN1A. Specifically, in cervical cancer, which is a cancer type that develops due to HPV persistent infection, Ding et al. (2017) have shown that lnc-CRNDE is overexpressed in cervical cancer tissues and functions as an oncogene since its expression correlates with increased cell proliferation, advanced FIGO stage and poor patient prognosis (Ding et al., 2017).
As for its mechanism of action, in cervical cancer CRNDE binds to PUMA, a modulator of apoptosis, to inhibit its expression. As a result, there is an increase in the expression of Bcl-xL, an anti-apoptotic protein that plays a crucial role in cell death (Figure 5). Interestingly, Bcl-xL is overexpressed in patients with chemoresistance and shorter disease-free intervals (Zhang and Fan 2019).
Further studies in cervical cancer by Lu Ren et al. (2021) determined that CRNDE is a sponge for miR-4262, and interestingly, this miRNA targets ZEB1, a regulator of epithelial to mesenchymal transition (Figure 5). Thus, the CRNDE/miR-4262/ZEB1 signaling axis accelerates cervical cancer progression by increasing cell survival, EMT, migration, and invasion (Ren et al., 2021).
lncRNA-SNHG12 -The lncRNA-SNHG12 (Homo sapiens small nucleolar RNA host gene 12) has been studied in several types of cancer; however, studies by Lai et al. (2020), show that SNHG12 is highly expressed in CC cells infected with HPV. Functional assays showed that SNHG12 contributes to the malignant behavior of CC since SNHG12 participates in the activation of the ERK/Slug signaling pathway, which leads to the promotion of epithelial-mesenchymal transition, an increase in proliferation, invasion, and inhibition of apoptosis (Lai et al., 2020). HPV16-E6 and E7 proteins can drive SNHG12 expression through a c-Myc activation-dependent mechanism; however, the specific mechanism has not been fully described (Figure 5).
5.2.2 Suppressing Long Non-Coding RNAs
LncRNA-DINO -Recent studies show that the Damage Induced lncRNA (DINO) is a 951 nucleotides long lncRNA located on chromosome six and divergently transcribed to the CDKN1A gene (Schmitt et al., 2016). Schmitt and coworkers have shown that DINO is transcribed under p53 regulation upon DNA damage stimuli and other stress factors using human cancer cell lines. DINO, once transcribed, forms a complex with p53 to increase its stability and enhance the signaling of this protein (Schmitt et al., 2016).
Interestingly, it has been demonstrated that DINO expression is increased in cells expressing the HPV-16 E7 protein (Sharma and Munger 2020b). It is also a critical mediator of p53 stabilization and activation. It has become evident that HPV16-E7 positive cells are more sensitive to cell death under various metabolic stress conditions; this association has been attributed mainly to p53 activation (Figure 5) (Sharma and Munger 2020b).
Additionally, DINO expression presents a regulation mediated by epigenetic derepression (H3K27 demethylation) through KDM6A. In addition, in cervical cancer cells, it has been shown that DINO activation occurs through a pathway that is independent of ATM/CHK2 (Sharma and Munger 2020a). Recent studies show that DINO is silenced in human cancers (Liu et al., 2019; Marney et al., 2021), and some evidence shows that a region of the DINO/CDKN1A locus is hypermethylated, leading to the silencing of DINO but not CDKN1A (Marney et al., 2021). Furthermore, silencing DINO in a model of primary HFKs E7+ HPV16+ using shRNAs decreases sensitivity to doxorubicin (Sharma and Munger 2020b).
Studies by Marney et al. (2021) have shown that loss of one or two DINO alleles alters the p53-mediated signaling and apoptosis pathway, resulting in an inefficient tumor suppressor phenotype. These insights are of great interest as new therapeutic strategies could aim to activate p53 through the reversal of epigenomic silencing of DINO.
LncRNA-PRINS -To understand HPV-linked HNSCC progression, Magda Kopczynska et al. (2020) performed analyses using data from The Cancer Genome Atlas (TCGA) (Kopczynska et al., 2020). In this analysis, the authors included 81 HPV (−) and 40 HPV (+) samples from patients with HNSCC. The results showed differentially expressed lncRNAs in HPV (+) such as PRINS, CDKN2B-AS1, TTTY14, TTTY15, MEG3, and H19 vs. HPV (−) patients (Kopczynska et al., 2020).
Interestingly, the Psoriasis susceptibility-related RNA Gene Induced by Stress lncRNA (lnc-PRINS) was overexpressed in HVP (+) patients and was associated with better overall survival and disease-free survival. Lnc-PRINS expression levels correlate with antiviral response genes, e.g., CRTSS, TLR8, IRF5, CCL5, CD40, NOD, CARD9, PYCARD, PSTPIP, and IFNAR1. In addition, lnc-PRINS modulates gene expression of transcripts involved in the inflammatory response, growth factors, proinflammatory genes, chemokines, and immunostimulators among others, so lnc-PRINS may stimulate the immune response of patients with HPV + HNSCC (Kopczynska et al., 2020). Some studies have indicated that tumors from HPV (+) patients have an increased number of immune cells, which may explain why HPV + HNSCC patients have a better response to therapy than HPV-negative patients (Figure 5).
6 Implications for Triage and Prognosis
DNA methylation status has been proposed in CC screening programs as a triage strategy for women with positive HPV results and cases with atypical squamous cells of undetermined significance (ASCUS). If not diagnosed and treated on time, the HPV persistent infection can generate precursor lesions that become cervical cancer. For this reason, once HPV+ cases have been identified in the screening programs, a second test (triage) should be carried out to differentiate between cancer precursor lesions that could potentially progress (persistent infection) from those that do not (transient infection). At this second triage, determining DNA methylation of a group of cellular or viral genes has enormous diagnostic utility for HPV (+) cases.
Initial studies postulated the methylation of CADM1 and MAL to detect cervical intraepithelial neoplasia (CIN) 2 and 3, being more sensitive as a triage strategy than cytology (Verhoef et al., 2015). Methylation of the transcription factor ZNF582 is associated with the presence of CIN3 in ASCUS patients (Liou et al., 2015), and ZNF582, together with PAX1, showed high sensitivity and specificity for HSIL detection (Tian et al., 2017). Several genes have been proposed whose methylation is a potential biomarker of precursor lesions; among them, we found POUF4 (Pun et al., 2015; Kocsis et al., 2017), ASCL1, LHX8, and ST6GALNAC5 (Verlaat et al., 2018).
Among the clinical studies that validate methylation profiles as biomarkers of CIN2, CIN3, and CC, is the S5 trial. S5 evaluates the methylation of EPB41L3 and the late region of HPV16, −18, −31, and −33 with a sensitivity of 74% and a specificity of 90% (Lorincz et al., 2016; Cook et al., 2019). Additionally, in cases with ASCUS, S5 can discriminate between low-grade (>CIN1 LSIL) and high-grade (CIN2+, HSIL) lesions (Hernandez-Lopez et al., 2019; Gu et al., 2020). Recently, a multicenter analysis demonstrated a sensitivity of 99.8% for detecting CIN3 and invasive carcinoma (Banila et al., 2022). The EPB41L3, L1, L2, and E2 (HPV16) panel has also been evaluated for the early diagnosis of oropharyngeal carcinoma. The use of his panel has a sensitivity of 70% and a specificity of 91%, and it only requires a gargle sample (Giuliano et al., 2020).
Subsequent studies confirmed the usefulness of EPB41L3, SOX1, and HS3ST2 as a triage method for HSIL detection (Clarke et al., 2017). Studies in the Latin American population have modified the cut-off point of this assay (3.1, 3.7) to improve sensitivity and specificity (Hernandez-Lopez et al., 2019; Ramirez et al., 2021).
Two panels of methylation markers to detect CIN2+ have been validated in the Dutch and Chinese populations (C13ORF18/EBP41L3/JAM3 y C13ORF18/ANKRD18CP/JAM3). These assays showed a sensitivity of 79% and 76% and a specificity of 57% and 65%, respectively (Li N et al., 2021). Other studies evaluated the ASCL1, LHX8, ST6GALNAC5, GHSR, ZIC1, and SST genes. Interestingly, the ASCL1/LHX8 genes showed greater specificity in identifying HPV-HR+ CIN3 with a similar sensitivity to HPV16/18 genotyping (Dick et al., 2020; Verhoef et al., 2022). This panel has been evaluated in urine samples to discriminate between LSIL and HSIL, managing to distinguish non-tumor tissue and low-grade lesions versus high-grade lesions, showing promising results, so it is suggested as a strategy to increase screening and triage (Van Keer et al., 2021).
Furthermore, long-term studies (14 years) showed that the cumulative incidence of CC in HPV+ samples is higher using FAM19A4 and miRNA-124 methylation profiles compared with a positive cytological result as a triage strategy (Luttmer et al., 2016; Clarke et al., 2017; Del Mistro et al., 2017; De Strooper et al., 2018; Dick et al., 2019). This assessment indicates that a negative result in FAM19A4/miR124-2 methylation could rule out the existence of CC (Vink et al., 2021).
Recent studies show that the FAM19A4/miR124-2 assay performed better than HPV16, 18, 31, 33, and 45 genotyping in HPV-positive women with borderline or moderate dyskaryosis, ASC, or LSIL, since the absolute risk of CIN3+ was higher (33%) than with other methods. Interestingly, combining the FAM19A4/miR124-2 assay with HPV16/18 increased detection by 40%. This strategy would allow more significant differentiation of patients with CIN3+ (Dick et al., 2022). Additionally, the methylation profile of FAM19A4/miR124-2 is constant in lesions caused by HPV 16 and 18 or by other high-risk genotypes (Leeman et al., 2018; Leeman et al., 2019). These genes were evaluated in European multicenter studies and have been assessed as a good triage option and an alternative to cytology (Bonde et al., 2021).
Other miRNAs such as miR-15b-5p and miR-375 have been evaluated together with FAM19A4, and this approach manages to distinguish CIN3 cases (Babion et al., 2018). In addition, the panel: C13orf18, EPB41L3 and JAM3, was validated for the CIN3 and CC detection (van Leeuwen et al., 2019). Assays for detecting malignant neoplasms associated with HPV-HR (+) based on epigenetic signatures are commercially available, for instance, GynTect® and QIAsure®. GynTect® is a molecular test for cervical carcinoma that analyzes the methylation profiles of 6 markers (ASTN1, ZNF671, DLX1, ITGA4, RXFP3, SOX17) (Schmitz et al., 2017; Schmitz et al., 2018). QIAsure® is a multiplex real-time methylation-specific PCR-based (qMSP) assay that detects hypermethylation of FAM19A4 and mir124-2 (Floore et al., 2019). It should be noted that the panel is consistent in CC, even with rare histotypes and negative HPV-HR carcinomas (Vink et al., 2019). A comparative study of both kits showed that GynTec has higher specificity than QIASure to detect CIN2+ (87% vs. 67%) and CIN3 (84% vs. 68%) (Dippmann et al., 2020).
Other commercially available panels for methylation evaluation are Precursor M Gold (CADM1, MAL, and miR124-2), CONFIDENCE Marker (POU4F3), and CERVI-M (PAX1) (Kremer et al., 2021). PAX1, in particular, has been evaluated in the Taiwanese population, finding comparable performance with cytology and better specificity than HPV16/18 genotyping as a triage strategy for the detection of CIN3+ in HPV-HR+ women (Chang et al., 2021). Interestingly, changes in viral DNA methylation have been associated with the transition from transient to persistent infection associated with cancer precursor lesions. Also, several studies have proposed the methylation of L1/L2 regions of carcinogenic HPV type 1 for the CIN3 and adenocarcinoma in situ detection (Clarke et al., 2018).
A meta-analysis identified the DNA methylation of CADM1, MAL, MIR-124-2, FAM19A4, POU4F3, EPB41L3, PAX1, SOX1, and L1/L2 from HPV16 as a triage strategy for HPV-AR+ cases (Kelly et al., 2019).
These studies indicate that in CC, the HPV test and triage based on DNA methylation, identify with higher specificity and sensitivity, healthy vs. precursor lesions. All those methylation biomarkers validated, could be extended to other HPV (+) neoplasms to early detection.
7 Epigenetic Therapy and Chemoprevention
The dynamic and reversible nature of epigenetic modifications makes them ideal therapeutic targets for treating HPV-associated cancers. Several studies have analyzed the effect of epigenetic therapies and chemopreventive agents in both cell lines and HPV+ neoplasms, including CC and HNSCC.
The potential of different natural compounds derived from plant sources such as phytochemicals has been evaluated as chemopreventive or adjuvants in anticancer therapies for several neoplasms. Exploring the potential therapeutic use of these compounds is particularly important since there are no specific treatment strategies for HPV + HNSCC.
Accumulating evidence suggests that dietary intake of vitamin B and folates have preventive effects against CC (Piyathilake et al., 2014; Mahmoud and Ali 2019). Another therapeutic resource for HPV neoplasms is the systemic administration of retinoids such as vitamin A (Oren-Shabtai et al., 2021), due to the retinoid deficiency is associated with the methylation of the retinoic acid receptor β 2 (RARβ2), a tumor suppressor. In a transgenic mouse model, it was shown that the expression of HPV E7 and the lack of dietary retinoids induce proliferation, hyperkeratinization, development of precursor lesions, and cancer in situ (Ocadiz-Delgado et al., 2021).
Analyses using molecular modeling have predicted that genistein, a naturally occurring compound found in soybean, interacts with several DNMT and HDAC family members. This compound reduces the expression and enzymatic activity of DNMTs and HDACs in Hela cells. Genistein is also known to activate the expression of tumor suppressor genes such as MGMT, RARb, p21, E-cadherin, DAPK1, FTHI, RUNX3, CDH1, PTEN, SOCS51, and other genes related to cell cycle regulation, migration, and inflammation such as CCNB1, TWIST1, MMP14, TERT, AKT1, PTPRR, FOS, and ILIA (Sundaram et al., 2018; Sundaram et al., 2019). These data suggest that genistein could have chemopreventive and chemotherapeutic properties since it exhibits anti-inflammatory, antiangiogenic, antiproliferative, and pro-apoptotic activities.
Curcumin (diferuloylmethane), a polyphenol found in the rhizome of curcuma plants, has been reported to inhibit the transcription of E6 and E7, thus restoring the expression of p53, pRB, and PTEN (Maher et al., 2011). In CC cell lines, it was shown that curcumin improves paclitaxel-induced apoptosis through the NF-κB-p53-caspase-3 (Dang et al., 2015). Curcumin also suppresses proliferation and exhibits anticancer effects in oral squamous cell carcinoma (OSCC). Interestingly, this compound can chemo-sensitize CSCs by abrogating stemness and inhibiting growth more effectively in HVP+ cells. These pieces of research propose that curcumin is a strong candidate for therapeutic approaches, especially in HPV+ tumors (Bano et al., 2018). Mechanistically, curcumin also acts by inhibiting HDAC 1, 2, 3, 4, 5, 6, 8, and 11, resulting in increased histone acetylation levels (Soflaei et al., 2018).
The FDA has approved several HDAC inhibitors as therapeutics with anticancer potential. Accumulating evidence has demonstrated that HDAC inhibitors are helpful for the treatment of hematological malignancies such as multiple myeloma and leukemia. Meanwhile, the efficiency of HDAC inhibitors in solid tumors is currently being assessed.
In vitro assays show that HDAC inhibitors such as sodium butyrate and trichostatin A can reverse the HPV-induced proliferation. These inhibitors neutralize the viral oncoproteins, increasing p21 and p27 expression, resulting in the Cyclin-dependent kinases suppression activity, arrest in the G1 phase of the cell cycle, and therefore, p53-independent and p73-dependent apoptosis is induced (Finzer et al., 2001; Finzer et al., 2004). Also, in CC cell lines showed that hydralazine and valproic acid exhibited an antiproliferative effect dependent on upregulation and stability of p53 by acetylation. However, they also increased the expression of E6/E7 (de la Cruz-Hernandez et al., 2007). Interestingly, when these drugs were administered to patients, the expression of viral oncoproteins was not increase (NCT00404326). Later, it was identified that combining both drugs increased radiation cytotoxicity, enhanced by cisplatin (Mani et al., 2014). Subsequent phase III clinical trials in patients with advanced CC showed a more prolonged median survival (10 months) when treated with chemotherapy + valproic acid + hydralazine than chemotherapy + placebo (6 months). However, toxicities such as thrombocytopenia were higher in the experimental group (Coronel et al., 2011). Other HDACs inhibitors such as vorinostat, belinostat, and panobinostat reduce E6 and E7 activity resulting in apoptosis and abrogating viral replication (Banerjee et al., 2018).
It has been postulated that the failure of HDAC inhibitors is due to the activation of the LIFRα oncogenic signaling pathways since EC359 (inhibitor of LIFRα) attenuates side effects provoked by HDAC inhibitors (Li M et al., 2021). Another aspect that should be further evaluated is the epigenetic classification of tumors and their relationship with the response to HDAC inhibitors.
Regarding the combination of HDAC inhibitors with antivirals to eradicate HPV infections, hitherto, there are no HPV-specific antiviral agents approved for clinical use, however, machine learning and therapeutic repositioning are emerging strategies to identify antiviral molecules. Lin et al. (2021) identified 57 antiviral drugs with potential interaction with E1, E2, E4, E5, E6, E7 L1, and E8^E2 (Lin et al., 2021). However, the effectiveness of these candidates and their impact on epigenetic alterations have not yet evaluated.
In addition, HDAC inhibitors have been used in combination with other drugs that block the activity of immune checkpoint proteins, for example, Chidamide is a benzamide that selectively suppresses the activity of class I and class IIb HDACs, increasing the acetylation of histones H3/H4, thus promoting transcription. Chidamide is currently in phase Ib/II clinical trials (NCT04651127). This inhibitor sensitizes advanced and metastatic, persistent, or recurrent cervical cancer tumors for further treatment with PD-1 inhibitors to boost the immune response against cancer cells.
In transformed follicular lymphoma, chidamide impacts on the PI3K-AKT pathway (Zhong et al., 2021) and in lymphoma it causes G0/G1 arrest, declining the cyclin D and c-myc expression, and increasing p53 and p21. Additionally, chidamide induces caspase-dependent apoptosis (Yuan et al., 2019).
PEVO trial evaluates the response to the combination of pembrolizumab and vorinostat, HDAC inhibitor, in several tumors recurrent and metastatic squamous cell carcinomas, HPV+ vs. HPV−. This study might provide clues to determine how HPV presence influences tumors’ therapeutic response and could also identify predictive and pharmacodynamic biomarkers (de Guillebon et al., 2021).
Other therapies include the epigenetic factor DNMT. 5-azacytidine and decitabine are drugs approved by the FDA to inhibit DNMTs for myelodysplastic syndrome and acute myeloid leukemia treatment. In HPV + solid tumors, has been evaluated the potential of these drugs as monotherapy or in combination with chemotherapy or immunotherapy.
In murine models and clinical trials in HNSCC HPV+ patients (NCT02178072), the 5-azacitidine treatment decreases proliferation, expression of viral genes, the activity of metalloproteases and induces p53-dependent apoptosis (Biktasova et al., 2017). In another clinical trial, intravenous administration of decitabine was evaluated for patients with HPV+ HNSCC anogenital carcinomas at high risk of recurrence, after radiotherapy (NCT04252248). However, there are still no conclusive results from this trial. Interestingly, it has been proposed to evaluate the effect of orally administered decitabine, in combination with durvalumab for patients with oral cavity, oropharynx, hypopharynx, or larynx recurrent or metastatic carcinomas. In this clinical trial, the response in HPV+ vs. HPV− patients will be compared, evaluating the APOBEC expression, the interferon pathway activation, proliferation, and apoptosis (NCT03019003).
An additional epigenetic factor of great relevance is EZH2. Since EZH2 is more expressed in HPV+ OPSCC tumors than in HPV− cases (Idris et al., 2016), there are therapies employed to inhibit EZH2, such as siRNAs, 3-deazaneplanocin A, GSK-343, and EPZ005687. Immunotherapy has undoubtedly improved the life expectancy of patients with malign neoplasms such as melanoma. Nevertheless, there are still several types of tumors where the impact of immunotherapy remains to be proven. The combination of immunotherapy and epitherapy has recently been proposed in different neoplasms.
The epigenetic effects of chemopreventives and therapeutic agents in HPV-associated neoplasms are still a subject of study from a functional and clinical perspective.
8 Final Remarks
Unlike genetic alterations, epigenetic alterations are potentially reversible. Epigenetic, systemic or targeted therapy is an area with great potential for the treatment of several neoplasms. These therapies must be directed according to the biological characteristics of tumor, in order to provide the greatest benefit and reduce adverse effects. Epigenetic alterations are a tool for monitoring therapeutic responses. The control of neoplasms associated with HPV is a feasible challenge to achieve since there are prophylactic vaccines, screening, and triage methods, which are focused on cervical cancer, however, in other neoplasms these strategies can also have an impact favorable. The discovery of epigenetic alterations induced by HPV, and new viral non-coding RNAs and methylation signatures, opens the possibilities for their application as triage biomarkers, and for designing drugs to modify the epigenetic landscape for improved patient survival and quality of life based on comprehensive treatments.
Author Contributions
The authors, PP-S and RC-O contributed to bibliography review and writing of this article. We declare that the images and graphics are original and were created from the articles cited.
Funding
This article is part of the CONACyT-funded project Salud 2017- C01-290427. Project R-2017-3602-31.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
Amaro-Filho, S. M., Pereira Chaves, C. B., Felix, S. P., Basto, D. L., de Almeida, L. M., and Moreira, M. A. M. (2018). HPV DNA Methylation at the Early Promoter and E1/E2 Integrity: A Comparison between HPV16, HPV18 and HPV45 in Cervical Cancer. Papillomavirus Res. 5, 172–179. doi:10.1016/j.pvr.2018.04.002
Arbyn, M., Tommasino, M., Depuydt, C., and Dillner, J. (2014). Are 20 Human Papillomavirus Types Causing Cervical Cancer? J. Pathol. 234, 431–435. doi:10.1002/path.4424
Babion, I., Snoek, B. C., Novianti, P. W., Jaspers, A., van Trommel, N., Heideman, D. A. M., et al. (2018). Triage of High-Risk HPV-Positive Women in Population-Based Screening by miRNA Expression Analysis in Cervical Scrapes; a Feasibility Study. Clin. Epigenet 10, 76. doi:10.1186/s13148-018-0509-9
Banerjee, N. S., Moore, D. W., Broker, T. R., and Chow, L. T. (2018). Vorinostat, a Pan-HDAC Inhibitor, Abrogates Productive HPV-18 DNA Amplification. Proc. Natl. Acad. Sci. U. S. A. 115, E11138–E47. doi:10.1073/pnas.1801156115
Banila, C., Lorincz, A. T., Scibior‐Bentkowska, D., Clifford, G. M., Kumbi, B., Beyene, D., et al. (2022). Clinical Performance of Methylation as a Biomarker for Cervical Carcinoma In Situ and Cancer Diagnosis: A Worldwide Study. Intl J. Cancer 150, 290–302. doi:10.1002/ijc.33815
Bano, N., Yadav, M., and Das, B. C. (2018). Differential Inhibitory Effects of Curcumin between HPV+ve and HPV-Ve Oral Cancer Stem Cells. Front. Oncol. 8, 412. doi:10.3389/fonc.2018.00412
Barr, J. A., Hayes, K. E., Brownmiller, T., Harold, A. D., Jagannathan, R., Lockman, P. R., et al. (2019). Long Non-coding RNA FAM83H-AS1 Is Regulated by Human Papillomavirus 16 E6 Independently of P53 in Cervical Cancer Cells. Sci. Rep. 9, 3662. doi:10.1038/s41598-019-40094-8
Bee, K. J., Gradissimo, A., Chen, Z., Harari, A., Schiffman, M., Raine-Bennett, T., et al. (2021). Genetic and Epigenetic Variations of HPV52 in Cervical Precancer. Int. J. Mol. Sci. 22, 6463. doi:10.3390/ijms22126463
Berglund, A., Muenyi, C., Siegel, E. M., Ajidahun, A., Eschrich, S. A., Wong, D., et al. (2022). Characterization of Epigenomic Alterations in HPV16+ Head and Neck Squamous Cell Carcinomas. Cancer Epidemiol. Biomarkers Prev. 31 (4), 858–869. doi:10.1158/1055-9965
Bergvall, M., Melendy, T., and Archambault, J. (2013). The E1 Proteins. Virology 445, 35–56. doi:10.1016/j.virol.2013.07.020
Bernard, H.-U. (2013). Regulatory Elements in the Viral Genome. Virology 445, 197–204. doi:10.1016/j.virol.2013.04.035
Biktasova, A., Hajek, M., Sewell, A., Gary, C., Bellinger, G., Deshpande, H. A., et al. (2017). Demethylation Therapy as a Targeted Treatment for Human Papillomavirus-Associated Head and Neck Cancer. Clin. Cancer Res. 23, 7276–7287. doi:10.1158/1078-0432.ccr-17-1438
Blachon, S., Bellanger, S., Demeret, C., and Thierry, F. (2005). Nucleo-cytoplasmic Shuttling of High Risk Human Papillomavirus E2 Proteins Induces Apoptosis. J. Biol. Chem. 280, 36088–36098. doi:10.1074/jbc.m505138200
Boers, A., Wang, R., van Leeuwen, R. W., Klip, H. G., de Bock, G. H., Hollema, H., et al. (2016). Discovery of New Methylation Markers to Improve Screening for Cervical Intraepithelial Neoplasia Grade 2/3. Clin. Epigenet 8, 29. doi:10.1186/s13148-016-0196-3
Bonde, J., Floore, A., Ejegod, D., Vink, F. J., Hesselink, A., Ven, P. M., et al. (2021). Methylation Markers FAM19A4 and miR124 ‐2 as Triage Strategy for Primary Human Papillomavirus Screen Positive Women: A Large European Multicenter Study. Int. J. Cancer 148, 396–405. doi:10.1002/ijc.33320
Brennan, K., Koenig, J. L., Gentles, A. J., Sunwoo, J. B., and Gevaert, O. (2017). Identification of an Atypical Etiological Head and Neck Squamous Carcinoma Subtype Featuring the CpG Island Methylator Phenotype. EBioMedicine 17, 223–236. doi:10.1016/j.ebiom.2017.02.025
Bruni, L., Albero, G., Serrano, B., Mena, M., Gómez, D., Muñoz, J., et al. (2021). “Human Papillomavirus and Related Diseases in the World,” in ICO/IARC Information Centre on HPV and Cancer (HPV Information Centre).
Buck, C. B., Day, P. M., and Trus, B. L. (2013). The Papillomavirus Major Capsid Protein L1. Virology 445, 169–174. doi:10.1016/j.virol.2013.05.038
Burgers, W. A., Blanchon, L., Pradhan, S., Launoit, Y. d., Kouzarides, T., and Fuks, F. (2007). Viral Oncoproteins Target the DNA Methyltransferases. Oncogene 26, 1650–1655. doi:10.1038/sj.onc.1209950
Burley, M., Roberts, S., and Parish, J. L. (2020). Epigenetic Regulation of Human Papillomavirus Transcription in the Productive Virus Life Cycle. Semin. Immunopathol. 42, 159–171. doi:10.1007/s00281-019-00773-0
Burnett, T. S., and Sleeman, J. P. (1984). Uneven distribution of methylation sites within the human papillomavirus la genome: possible relevance to viral gene expression. Nucl. Acids Res. 12, 8847–8860. doi:10.1093/nar/12.23.8847
Castro-Oropeza, R., Melendez-Zajgla, J., Maldonado, V., and Vazquez-Santillan, K. (2018). The Emerging Role of lncRNAs in the Regulation of Cancer Stem Cells. Cell Oncol. 41, 585–603. doi:10.1007/s13402-018-0406-4
Chaiwongkot, A. S., Vinokurova, C., Pientong, T., Ekalaksananan, B., Kongyingyoes, P., Kleebkaow, B., et al. (2013). Differential Methylation of E2 Binding Sites in Episomal and Integrated HPV 16 Genomes in Preinvasive and Invasive Cervical Lesions. Int. J. Cancer 132, 2087–2094. doi:10.1002/ijc.27906
Chang, C.-L., Ho, S.-C., Su, Y.-F., Juan, Y.-C., Huang, C.-Y., Chao, A.-S., et al. (2021). DNA Methylation Marker for the Triage of hrHPV Positive Women in Cervical Cancer Screening: Real-World Evidence in Taiwan. Gynecol. Oncol. 161, 429–435. doi:10.1016/j.ygyno.2021.02.011
Chen, L., Dong, B., Gao, H., Xue, H., Pan, D., and Sun, P. (2020). HPV-16 E2/E6 and POU5F1B as Biomarkers to Determine Cervical High-Grade Squamous Lesions and More. Jir 13, 813–821. doi:10.2147/jir.s278911
Chen, Y., Wang, C. X., Sun, X. X., Wang, C., Liu, T. F., and Wang, D. J. (2017). Long Non-coding RNA CCHE1 Overexpression Predicts a Poor Prognosis for Cervical Cancer. Eur. Rev. Med. Pharmacol. Sci. 21, 479–483.
Chesson, H. W., Dunne, E. F., Hariri, S., and Markowitz, L. E. (2014). The Estimated Lifetime Probability of Acquiring Human Papillomavirus in the United States. Sex. Transm. Dis. 41, 660–664. doi:10.1097/olq.0000000000000193
Cheung, J. L. K., Cheung, T.-H., Yu, M. Y., and Chan, P. K. S. (2013). Virological Characteristics of Cervical Cancers Carrying Pure Episomal Form of HPV16 Genome. Gynecol. Oncol. 131, 374–379. doi:10.1016/j.ygyno.2013.08.026
Chirayil, R., Kincaid, R. P., Dahlke, C., Kuny, C. V., Dälken, N., Spohn, M., et al. (2018). Identification of Virus-Encoded microRNAs in Divergent Papillomaviruses. PLoS Pathog. 14, e1007156. doi:10.1371/journal.ppat.1007156
Clarke, M. A., Gradissimo, A., Schiffman, M., Lam, J., Sollecito, C. C., Fetterman, B., et al. (2018). Human Papillomavirus DNA Methylation as a Biomarker for Cervical Precancer: Consistency across 12 Genotypes and Potential Impact on Management of HPV-Positive Women. Clin. Cancer Res. 24, 2194–2202. doi:10.1158/1078-0432.ccr-17-3251
Clarke, M. A., Luhn, P., Gage, J. C., Bodelon, C., Dunn, S. T., Walker, J., et al. (2017). Discovery and Validation of Candidate Host DNA Methylation Markers for Detection of Cervical Precancer and Cancer. Int. J. Cancer 141, 701–710. doi:10.1002/ijc.30781
Cook, D. A., Krajden, M., Brentnall, A. R., Gondara, L., Chan, T., Law, J. H., et al. (2019). Evaluation of a Validated Methylation Triage Signature for Human Papillomavirus Positive Women in the HPV FOCAL Cervical Cancer Screening Trial. Int. J. Cancer 144, 2587–2595. doi:10.1002/ijc.31976
Coronel, J., Cetina, L., Pacheco, I., Trejo-Becerril, C., González-Fierro, A., de la Cruz-Hernandez, E., et al. (2011). A Double-Blind, Placebo-Controlled, Randomized Phase III Trial of Chemotherapy Plus Epigenetic Therapy with Hydralazine Valproate for Advanced Cervical Cancer. Preliminary Results. Med. Oncol. 28 (Suppl. 1), S540–S546. doi:10.1007/s12032-010-9700-3
Cricca, M., Venturoli, S., Leo, E., Costa, S., Musiani, M., and Zerbini, M. (2009). Disruption of HPV 16 E1 and E2 Genes in Precancerous Cervical Lesions. J. Virological Methods 158, 180–183. doi:10.1016/j.jviromet.2009.01.005
Dang, Y.-P., Yuan, X.-Y., Tian, R., Li, D.-G., and Liu, W. (2015). Curcumin Improves the Paclitaxel-Induced Apoptosis of HPV-Positive Human Cervical Cancer Cells via the NF-Κb-P53-Caspase-3 Pathway. Exp. Ther. Med. 9, 1470–1476. doi:10.3892/etm.2015.2240
Das, D., Smith, N., Wang, X., and Morgan, I. M. (2017). The Deacetylase SIRT1 Regulates the Replication Properties of Human Papillomavirus 16 E1 and E2. J. Virol. 91, e00102–17. doi:10.1128/JVI.00102-17
Das, D., Bristol, M. L., Smith, N. W., James, C. D., Wang, X., Pichierri, P., et al. (2019). Correction for Morgan et al., "Werner Helicase Control of Human Papillomavirus 16 E1-E2 DNA Replication Is Regulated by SIRT1 Deacetylation. mBio 10, e01635–19. doi:10.1128/mbio.01635-19
de Guillebon, E., Jimenez, M., Mazzarella, L., Betsou, F., Stadler, P., Peták, I., et al. (2021). Combining Immunotherapy with an Epidrug in Squamous Cell Carcinomas of Different Locations: Rationale and Design of the PEVO Basket Trial. ESMO Open 6, 100106. doi:10.1016/j.esmoop.2021.100106
de la Cruz-Hernández, E., Pérez-Cárdenas, E., Contreras-Paredes, A., Cantú, D., Mohar, A., Lizano, M., et al. (2007). The Effects of DNA Methylation and Histone Deacetylase Inhibitors on Human Papillomavirus Early Gene Expression in Cervical Cancer, an In Vitro and Clinical Study. Virol. J. 4, 18. doi:10.1186/1743-422x-4-18
de Martel, C., Georges, D., Bray, F., Ferlay, J., and Clifford, G. M. (2020). Global Burden of Cancer Attributable to Infections in 2018: a Worldwide Incidence Analysis. Lancet Glob. Health 8, e180–e190. doi:10.1016/s2214-109x(19)30488-7
De Strooper, L. M. A., Berkhof, J., Steenbergen, R. D. M., Lissenberg‐Witte, B. I., Snijders, P. J. F., Meijer, C. J. L. M., et al. (2018). Cervical Cancer Risk in HPV‐positive Women after a Negative FAM19A4/mir124‐2 Methylation Test: A Post Hoc Analysis in the POBASCAM Trial with 14 Year Follow‐up. Int. J. Cancer 143, 1541–1548. doi:10.1002/ijc.31539
del Mar Peña, L. M., and Laimins, L. A. (2001). Differentiation-dependent Chromatin Rearrangement Coincides with Activation of Human Papillomavirus Type 31 Late Gene Expression. J. Virol. 75, 10005–10013. doi:10.1128/JVI.75.20.10005-10013.2001
Del Mistro, A., Frayle, H., Rizzi, M., Fantin, G., Ferro, A., Angeletti, P. M., et al. (2017). Methylation Analysis and HPV Genotyping of Self-Collected Cervical Samples from Women Not Responding to Screening Invitation and Review of the Literature. PLoS One 12, e0172226. doi:10.1371/journal.pone.0172226
de Villiers, M., Dillner, J., and Eklund, C. (2012). International HPV Reference Center. Stockholm, Sweden: Karolinska Institutet. Available at: https://ki.se/en/labmed/international-hpv-reference-center.
Dick, S., Kremer, W. W., De Strooper, L. M. A., Lissenberg-Witte, B. I., Steenbergen, R. D. M., Meijer, C. J. L. M., et al. (2019). Long-term CIN3+ Risk of HPV Positive Women after Triage with FAM19A4/miR124-2 Methylation Analysis. Gynecol. Oncol. 154, 368–373. doi:10.1016/j.ygyno.2019.06.002
Dick, S., Verhoef, L., De Strooper, L. M., Ciocănea-Teodorescu, I., Wisman, G. B. A., Meijer, C. J., et al. (2020). Evaluation of Six Methylation Markers Derived from Genome-wide Screens for Detection of Cervical Precancer and Cancer. Epigenomics 12, 1569–1578. doi:10.2217/epi-2019-0331
Dick, S., Vink, F. J., Heideman, D. A. M., Lissenberg-Witte, B. I., Meijer, C. J. L. M., and Berkhof, J. (2022). Risk-stratification of HPV-Positive Women with Low-Grade Cytology by FAM19A4/miR124-2 Methylation and HPV Genotyping. Br. J. Cancer 126, 259–264. doi:10.1038/s41416-021-01614-4
DiMaio, D., and Petti, L. M. (2013). The E5 Proteins. Virology 445, 99–114. doi:10.1016/j.virol.2013.05.006
Ding, J., Li, J., Wang, H., Tian, Y., Xie, M., He, X., et al. (2017). Long Noncoding RNA CRNDE Promotes Colorectal Cancer Cell Proliferation via Epigenetically Silencing DUSP5/CDKN1A Expression. Cell Death Dis. 8, e2997. doi:10.1038/cddis.2017.328
Dippmann, C., Schmitz, M., Wunsch, K., Schütze, S., Beer, K., Greinke, C., et al. (2020). Triage of hrHPV-Positive Women: Comparison of Two Commercial Methylation-specific PCR Assays. Clin. Epigenet 12, 171. doi:10.1186/s13148-020-00963-w
Dong, J., Teng, F., Guo, W., Yang, J., Ding, G., and Fu, Z. (2018). lncRNA SNHG8 Promotes the Tumorigenesis and Metastasis by Sponging miR-149-5p and Predicts Tumor Recurrence in Hepatocellular Carcinoma. Cell Physiol. Biochem. 51, 2262–2274. doi:10.1159/000495871
Doorbar, J. (2013). The E4 Protein; Structure, Function and Patterns of Expression. Virology 445, 80–98. doi:10.1016/j.virol.2013.07.008
Dreer, M., Fertey, J., van de Poel, S., Straub, E., Madlung, J., Macek, B., et al. (2016). Interaction of NCOR/SMRT Repressor Complexes with Papillomavirus E8^E2C Proteins Inhibits Viral Replication. PLoS Pathog. 12, e1005556. doi:10.1371/journal.ppat.1005556
Estêvão, D., Costa, N. R., Gil da Costa, R. M., and Medeiros, R. (2019). Hallmarks of HPV Carcinogenesis: The Role of E6, E7 and E5 Oncoproteins in Cellular Malignancy. Biochimica Biophysica Acta (BBA) - Gene Regul. Mech. 1862, 153–162. doi:10.1016/j.bbagrm.2019.01.001
Feber, A., Arya, M., de Winter, P., Saqib, M., Nigam, R., Malone, P. R., et al. (2015). Epigenetics Markers of Metastasis and HPV-Induced Tumorigenesis in Penile Cancer. Clin. Cancer Res. 21, 1196–1206. doi:10.1158/1078-0432.ccr-14-1656
Ferguson, J., Campos-León, K., Pentland, I., Stockton, J. D., Günther, T., Beggs, A. D., et al. (2021). The Chromatin Insulator CTCF Regulates HPV18 Transcript Splicing and Differentiation-dependent Late Gene Expression. PLoS Pathog. 17, e1010032. doi:10.1371/journal.ppat.1010032
Finzer, P., Krueger, A., Stöhr, M., Brenner, D., Soto, U., Kuntzen, C., et al. (2004). HDAC Inhibitors Trigger Apoptosis in HPV-Positive Cells by Inducing the E2F-P73 Pathway. Oncogene 23, 4807–4817. doi:10.1038/sj.onc.1207620
Finzer, P., Kuntzen, C., Soto, U., zur Hausen, H., and Rösl, F. (2001). Inhibitors of Histone Deacetylase Arrest Cell Cycle and Induce Apoptosis in Cervical Carcinoma Cells Circumventing Human Papillomavirus Oncogene Expression. Oncogene 20, 4768–4776. doi:10.1038/sj.onc.1204652
Floore, A., Hesselink, A., Oštrbenk, A., Alcaniz, E., Rothe, B., Pedersen, H., et al. (2019). Intra‐ and Inter‐laboratory Agreement of the FAM19A4/mir124‐2 Methylation Test: Results from an International Study. J. Clin. Lab. Anal. 33, e22854. doi:10.1002/jcla.22854
Fontan, C. T., Das, D., Bristol, M. L., James, C. D., Wang, X., Lohner, H., et al. (2020). Human Papillomavirus 16 (HPV16) E2 Repression of TWIST1 Transcription Is a Potential Mediator of HPV16 Cancer Outcomes. mSphere 5, e00981–20. doi:10.1128/msphere.00981-20
Gagliardi, A., Porter, V. L., Zong, Z., Bowlby, R., Titmuss, E., Namirembe, C., et al. (2020). Analysis of Ugandan Cervical Carcinomas Identifies Human Papillomavirus Clade-specific Epigenome and Transcriptome Landscapes. Nat. Genet. 52, 800–810. doi:10.1038/s41588-020-0673-7
Gallo, A., Miceli, V., Bulati, M., Iannolo, G., Contino, F., and Conaldi, P. G. (2020). 'Viral miRNAs as Active Players and Participants in Tumorigenesis. Cancers (Basel) 12, 358. doi:10.3390/cancers12020358
Gautam, D., Johnson, B. A., Mac, M., and Moody, C. A. (2018). SETD2-dependent H3K36me3 Plays a Critical Role in Epigenetic Regulation of the HPV31 Life Cycle. PLoS Pathog. 14, e1007367. doi:10.1371/journal.ppat.1007367
Giuliano, A. R., Nedjai, B., Lorincz, A. T., Schell, M. J., Rahman, S., Banwait, R., et al. (2020). Methylation of HPV 16 and EPB41L3 in Oral Gargles: Associations with Oropharyngeal Cancer Detection and Tumor Characteristics. Int. J. Cancer 146, 1018–1030. doi:10.1002/ijc.32570
Groves, I. J., Knight, E. L. A., Ang, Q. Y., Scarpini, C. G., and Coleman, N. (2016). HPV16 Oncogene Expression Levels during Early Cervical Carcinogenesis Are Determined by the Balance of Epigenetic Chromatin Modifications at the Integrated Virus Genome. Oncogene 35, 4773–4786. doi:10.1038/onc.2016.8
Gu, Y.-Y., Zhou, G.-N., Wang, Q., Ding, J.-X., and Hua, K.-Q. (2020). Evaluation of a Methylation Classifier for Predicting Pre-cancer Lesion Among Women with Abnormal Results between HPV16/18 and Cytology. Clin. Epigenet 12, 57. doi:10.1186/s13148-020-00849-x
Guttman, M., and Rinn, J. L. (2012). Modular Regulatory Principles of Large Non-coding RNAs. Nature 482, 339–346. doi:10.1038/nature10887
Han, Y., Ji, L., Guan, Y., Ma, M., Li, P., Xue, Y., et al. (2021). An Epigenomic Landscape of Cervical Intraepithelial Neoplasia and Cervical Cancer Using Single-Base Resolution Methylome and Hydroxymethylome. Clin. Transl. Med. 11, e498. doi:10.1002/ctm2.498
Harden, M. E., Prasad, N., Griffiths, A., and Munger, K. (2017). Modulation of microRNA-mRNA Target Pairs by Human Papillomavirus 16 Oncoproteins. mBio 8, e02170–1. doi:10.1128/mbio.02170-16
He, H., Liu, X., Liu, Y., Zhang, M., Lai, Y., Hao, Y., et al. (2019). Human Papillomavirus E6/E7 and Long Noncoding RNA TMPOP2 Mutually Upregulated Gene Expression in Cervical Cancer Cells. J. Virol. 93, e01808–18. doi:10.1128/JVI.01808-18
Hernández-López, R., Lorincz, A. T., Lorincz, A. T., Torres-Ibarra, L., Reuter, C., Scibior-Bentkowska, D., et al. Study Group (2019). 'Methylation Estimates the Risk of Precancer in HPV-Infected Women with Discrepant Results between Cytology and HPV16/18 Genotyping. Clin. Epigenetics 11, 140. doi:10.1186/s13148-019-0743-9
Hesselink, A. T., Heideman, D. A., Steenbergen, R. D., Gok, M., van Kemenade, F. J., Wilting, S. M., et al. (2014). 'Methylation Marker Analysis of Self-Sampled Cervico-Vaginal Lavage Specimens to Triage High-Risk HPV-Positive Women for Colposcopy. Int. J. Cancer 135, 880–886. doi:10.1002/ijc.28723
Holland, D., Hoppe-Seyler, K., Schuller, B., Lohrey, C., Maroldt, J., Durst, M., et al. (2008). 'Activation of the Enhancer of Zeste Homologue 2 Gene by the Human Papillomavirus E7 Oncoprotein. Cancer Res. 68, 9964–9972. doi:10.1158/0008-5472.can-08-1134
Huang, J. T., Chen, J. N., Gong, L. P., Bi, Y. H., Liang, J., Zhou, L., et al. (2019). 'Identification of Virus-Encoded Circular RNA. Virology 529, 144–151. doi:10.1016/j.virol.2019.01.014
Hyland, P. L., McDade, S. S., McCloskey, R., Dickson, G. J., Arthur, K., McCance, D. J., et al. (2011). 'Evidence for Alteration of EZH2, BMI1, and KDM6A and Epigenetic Reprogramming in Human Papillomavirus Type 16 E6/E7-Expressing Keratinocytes. J. Virol. 85, 10999–11006. doi:10.1128/jvi.00160-11
Iancu, I. V., Anton, G., Botezatu, A., Huica, I., Nastase, A., Socolov, D. G., et al. (2017). 'LINC01101 and LINC00277 Expression Levels as Novel Factors in HPV-Induced Cervical Neoplasia. J. Cell Mol. Med. 21, 3787–3794. doi:10.1111/jcmm.13288
IARC, International Agency Research Cancer. (2012). "Biological Agents: A Review of Human Carcinogens." In, edited by IARC Monographs on the Evaluation of Carcinogenic Risk to Humans. Lyon, France: IARC.
ICTV, International Committee on Taxonomy of Viruses (2020). Virus Taxonomy 2020 Release. Available at: https://talk.ictvonline.org/taxonomy/.
Idris, S., Lindsay, C., Kostiuk, M., Andrews, C., Cote, D. W., O'Connell, D. A., et al. (2016). 'Investigation of EZH2 Pathways for Novel Epigenetic Treatment Strategies in Oropharyngeal Cancer. J. Otolaryngol. Head. Neck Surg. 45, 54. doi:10.1186/s40463-016-0168-9
Ilahi, N. E., and Bhatti, A. (2020). 'Impact of HPV E5 on Viral Life Cycle via EGFR Signaling. Microb. Pathog. 139, 103923. doi:10.1016/j.micpath.2019.103923
Johannsen, E., and Lambert, P. F. (2013). 'Epigenetics of Human Papillomaviruses. Virology 445, 205–212. doi:10.1016/j.virol.2013.07.016
Kalantari, M., Osann, K., Calleja-Macias, I. E., Kim, S., Yan, B., Jordan, S., et al. (2014). 'Methylation of Human Papillomavirus 16, 18, 31, and 45 L2 and L1 Genes and the Cellular DAPK Gene: Considerations for Use as Biomarkers of the Progression of Cervical Neoplasia. Virology 448, 314–321. doi:10.1016/j.virol.2013.10.032
Kanaya, T., Kyo, S., and Laimins, L. A. (1997). 'The 5' Region of the Human Papillomavirus Type 31 Upstream Regulatory Region Acts as an Enhancer Which Augments Viral Early Expression through the Action of YY1. Virology 237, 159–169. doi:10.1006/viro.1997.8771
Karimzadeh, M., Arlidge, C., Rostami, A., Lupien, M., Bratman, S. V., and Hoffman, M. M. (2022). Human Papillomavirus Integration Transforms Chromatin to Drive Oncogenesis. New York: bioRxiv.
Kelly, H., Benavente, Y., Pavon, M. A., De Sanjose, S., Mayaud, P., and Lorincz, A. T. (2019). 'Performance of DNA Methylation Assays for Detection of High-Grade Cervical Intraepithelial Neoplasia (CIN2+): a Systematic Review and Meta-Analysis. Br. J. Cancer 121, 954–965. doi:10.1038/s41416-019-0593-4
Kocsis, A., Takacs, T., Jeney, C., Schaff, Z., Koiss, R., Jaray, B., et al. (2017). 'Performance of a New HPV and Biomarker Assay in the Management of hrHPV Positive Women: Subanalysis of the Ongoing Multicenter TRACE Clinical Trial (N > 6,000) to Evaluate POU4F3 Methylation as a Potential Biomarker of Cervical Precancer and Cancer. Int. J. Cancer 140, 1119–1133. doi:10.1002/ijc.30534
Kopczynska, M., Kolenda, T., Guglas, K., Sobocinska, J., Teresiak, A., Blizniak, R., et al. (2020). PRINS lncRNA Is a New Biomarker Candidate for HPV Infection and Prognosis of Head and Neck Squamous Cell Carcinomas. Diagn. (Basel) 10, 762. doi:10.3390/diagnostics10100762
Kremer, W. W., Steenbergen, R., Heideman, D., Kenter, G. G., and Meijer, C. (2021). The Use of Host Cell DNA Methylation Analysis in the Detection and Management of Women with Advanced Cervical Intraepithelial Neoplasia: a Review. BJOG 128, 504–514. doi:10.1111/1471-0528.16395
Kuasne, H., Colus, I. M., Busso, A. F., Hernandez-Vargas, H., Barros-Filho, M. C., Marchi, F. A., et al. (2015). 'Genome-wide Methylation and Transcriptome Analysis in Penile Carcinoma: Uncovering New Molecular Markers. Clin. Epigenetics 7, 46. doi:10.1186/s13148-015-0082-4
Lace, M. J., Yamakawa, Y., Ushikai, M., Anson, J. R., Haugen, T. H., and Turek, L. P. (2009). 'Cellular Factor YY1 Downregulates the Human Papillomavirus 16 E6/E7 Promoter, P97, In Vivo and In Vitro from a Negative Element Overlapping the Transcription-Initiation Site. J. Gen. Virol. 90, 2402–2412. doi:10.1099/vir.0.012708-0
Lai, S. Y., Guan, H. M., Liu, J., Huang, L. J., Hu, X. L., Chen, Y. H., et al. (2020). 'Long Noncoding RNA SNHG12 Modulated by Human Papillomavirus 16 E6/E7 Promotes Cervical Cancer Progression via ERK/Slug Pathway. J. Cell Physiol. 235, 7911–7922. doi:10.1002/jcp.29446
Langsfeld, E. S., Bodily, J. M., and Laimins, L. A. (2015). 'The Deacetylase Sirtuin 1 Regulates Human Papillomavirus Replication by Modulating Histone Acetylation and Recruitment of DNA Damage Factors NBS1 and Rad51 to Viral Genomes. PLoS Pathog. 11, e1005181. doi:10.1371/journal.ppat.1005181
Lechner, M., Fenton, T., West, J., Wilson, G., Feber, A., Henderson, S., et al. (2013). 'Identification and Functional Validation of HPV-Mediated Hypermethylation in Head and Neck Squamous Cell Carcinoma. Genome Med. 5, 15. doi:10.1186/gm419
Leeman, A., Del Pino, M., Marimon, L., Torne, A., Ordi, J., Ter Harmsel, B., et al. (2019). 'Reliable Identification of Women with CIN3+ Using hrHPV Genotyping and Methylation Markers in a Cytology-Screened Referral Population. Int. J. Cancer 144, 160–168. doi:10.1002/ijc.31787
Leeman, A., Ebisch, R. M. F., Kasius, A., Bosgraaf, R. P., Jenkins, D., van de Sandt, M. M., et al. (2018). Defining hrHPV Genotypes in Cervical Intraepithelial Neoplasia by Laser Capture Microdissection Supports Reflex Triage of Self-Samples Using HPV16/18 and FAM19A4/miR124-2 Methylation. Gynecol. Oncol. 151, 311–318. doi:10.1016/j.ygyno.2018.09.006
Li, M., Viswanadhapalli, S., Santhamma, B., Pratap, U. P., Luo, Y., Liu, J., et al. (2021). 'LIFR Inhibition Enhances the Therapeutic Efficacy of HDAC Inhibitors in Triple Negative Breast Cancer. Commun. Biol. 4, 1235. doi:10.1038/s42003-021-02741-7
Li, N., Hu, Y., Zhang, X., Liu, Y., He, Y., van der Zee, A. G. J., et al. (2021). 'DNA Methylation Markers as Triage Test for the Early Identification of Cervical Lesions in a Chinese Population. Int. J. Cancer 148, 1768–1777. doi:10.1002/ijc.33430
Lin, H. H., Zhang, Q. R., Kong, X., Zhang, L., Zhang, Y., Tang, Y., et al. (2021). 'Machine Learning Prediction of Antiviral-HPV Protein Interactions for Anti-HPV Pharmacotherapy. Sci. Rep. 11, 24367. doi:10.1038/s41598-021-03000-9
Liou, Y. L., Zhang, Y., Liu, Y., Cao, L., Qin, C. Z., Zhang, T. L., et al. (2015). 'Comparison of HPV Genotyping and Methylated ZNF582 as Triage for Women with Equivocal Liquid-Based Cytology Results. Clin. Epigenetics 7, 50. doi:10.1186/s13148-015-0084-2
Liu, H., Xu, J., Yang, Y., Wang, X., Wu, E., Majerciak, V., et al. (2021). Oncogenic HPV Promotes the Expression of the Long Noncoding RNA Lnc-FANCI-2 through E7 and YY1. Proc. Natl. Acad. Sci. U. S. A. 118, e2014195118. doi:10.1073/pnas.2014195118
Liu, Q., Xiao, Y., Cai, P., Li, J., and Li, D. (2019). Long Noncoding RNA DINO (Damage Induced Noncoding) Represses the Development of Gastric Cancer by Modulating P21 and Bcl-2 Associated X Protein (Bax) Expression. J. Cell Biochem. 1, 1. doi:10.1002/jcb.28394
Liu, S., de Medeiros, M. C., Fernandez, E. M., Zarins, K. R., Cavalcante, R. G., Qin, T., et al. (2020). '5-Hydroxymethylation Highlights the Heterogeneity in Keratinization and Cell Junctions in Head and Neck Cancers. Clin. Epigenetics 12, 175. doi:10.1186/s13148-020-00965-8
Lohavanichbutr, P., Houck, J., Fan, W., Yueh, B., Mendez, E., Futran, N., et al. (2009). 'Genomewide Gene Expression Profiles of HPV-Positive and HPV-Negative Oropharyngeal Cancer: Potential Implications for Treatment Choices. Arch. Otolaryngol. Head. Neck Surg. 135, 180–188. doi:10.1001/archoto.2008.540
Lorincz, A. T., Brentnall, A. R., Scibior-Bentkowska, D., Reuter, C., Banwait, R., Cadman, L., et al. (2016). 'Validation of a DNA Methylation HPV Triage Classifier in a Screening Sample. Int. J. Cancer 138, 2745–2751. doi:10.1002/ijc.30008
Lu, Y., Sha, H., Sun, X., Zhang, Y., Wu, Y., Zhang, J., et al. (2020). 'CRNDE: an Oncogenic Long Non-coding RNA in Cancers. Cancer Cell Int. 20, 162. doi:10.1186/s12935-020-01246-3
Luttmer, R., De Strooper, L. M., Dijkstra, M. G., Berkhof, J., Snijders, P. J., Steenbergen, R. D., et al. (2016). 'FAM19A4 Methylation Analysis in Self-Samples Compared with Cervical Scrapes for Detecting Cervical (Pre)cancer in HPV-Positive Women. Br. J. Cancer 115, 579–587. doi:10.1038/bjc.2016.200
Ma, X. S, Sheng, J., Wu, Y., Jiang, X., Gao, X., Cen, J., et al. (2017). LncRNAs as an Intermediate in HPV16 Promoting Myeloid-Derived Suppressor Cell Recruitment of Head and Neck Squamous Cell Carcinoma. Oncotarget 8, 42061–42075.
Ma, Y., Madupu, R., Karaoz, U., Nossa, C. W., Yang, L., Yooseph, S., et al. (2014). 'Human Papillomavirus Community in Healthy Persons, Defined by Metagenomics Analysis of Human Microbiome Project Shotgun Sequencing Data Sets. J. Virol. 88, 4786–4797. doi:10.1128/jvi.00093-14
Mac, M., and Moody, C. A. (2020). Epigenetic Regulation of the Human Papillomavirus Life Cycle. Pathogens 9, 483. doi:10.3390/pathogens9060483
Maher, D. M., Bell, M. C., O'Donnell, E. A., Gupta, B. K., Jaggi, M., and Chauhan, S. C. (2011). 'Curcumin Suppresses Human Papillomavirus Oncoproteins, Restores P53, Rb, and PTPN13 Proteins and Inhibits Benzo[a]pyrene-Induced Upregulation of HPV E7. Mol. Carcinog. 50, 47–57. doi:10.1002/mc.20695
Mahmoud, A. M., and Ali, M. M. (2019). Methyl Donor Micronutrients that Modify DNA Methylation and Cancer Outcome. Nutrients 11, 608. doi:10.3390/nu11030608
Mani, E., Medina, L. A., Isaac-Olive, K., and Duenas-Gonzalez, A. (2014). Radiosensitization of Cervical Cancer Cells with Epigenetic Drugs Hydralazine and Valproate. Eur. J. Gynaecol. Oncol. 35, 140–142.
Marney, C. B., Anderson, E. S., Adnan, M., Peng, K. L., Hu, Y., Weinhold, N., et al. (2021). 'p53-intact Cancers Escape Tumor Suppression through Loss of Long Noncoding RNA Dino. Cell Rep. 35, 109329. doi:10.1016/j.celrep.2021.109329
Marques Howarth, M., Simpson, D., Ngok, S. P., Nieves, B., Chen, R., Siprashvili, Z., et al. (2014). 'Long Noncoding RNA EWSAT1-Mediated Gene Repression Facilitates Ewing Sarcoma Oncogenesis. J. Clin. Invest 124, 5275–5290. doi:10.1172/jci72124
McBride, A. A. (2013). The Papillomavirus E2 Proteins. Virology 445, 57–79. doi:10.1016/j.virol.2013.06.006
Mendez-Matias, G., Velazquez-Velazquez, C., Castro-Oropeza, R., Mantilla-Morales, A., Ocampo-Sandoval, D., Burgos-Gonzalez, A., et al. (2021). Prevalence of HPV in Mexican Patients with Head and Neck Squamous Carcinoma and Identification of Potential Prognostic Biomarkers. Cancers (Basel) 13, 5602. doi:10.3390/cancers13225602
Miao, Q., Ni, B., and Tang, J. (2021). Coding Potential of circRNAs: New Discoveries and Challenges. PeerJ 9, e10718. doi:10.7717/peerj.10718
Misawa, K., Imai, A., Matsui, H., Kanai, A., Misawa, Y., Mochizuki, D., et al. (2020). Identification of Novel Methylation Markers in HPV-Associated Oropharyngeal Cancer: Genome-wide Discovery, Tissue Verification and Validation Testing in ctDNA. Oncogene 39, 4741–4755. doi:10.1038/s41388-020-1327-z
Nakagawa, T., Matsusaka, K., Misawa, K., Ota, S., Fukuyo, M., Rahmutulla, B., et al. (2020). Stratification of HPV-Associated and HPV-Negative Oropharyngeal Squamous Cell Carcinomas Based on DNA Methylation Epigenotypes. Int. J. Cancer 146, 2460–2474. doi:10.1002/ijc.32890
Ocadiz-Delgado, R., Serafin-Higuera, N., Alvarez-Rios, E., Garcia-Villa, E., Tinajero-Rodriguez, M., Rodriguez-Uribe, G., et al. (2021). Vitamin A Deficiency in K14E7HPV Expressing Transgenic Mice Facilitates the Formation of Malignant Cervical Lesions. APMIS 129, 512–523. doi:10.1111/apm.13159
Oren-Shabtai, M., Snast, I., Lapidoth, M., Sherman, S., Noyman, Y., Mimouni, D., et al. (2021). Topical and Systemic Retinoids for the Treatment of Genital Warts: A Systematic Review and Meta-Analysis. Dermatology 237, 389–395. doi:10.1159/000511398
Ouyang, J., Hu, J., and Chen, J. L. (2016). lncRNAs Regulate the Innate Immune Response to Viral Infection. Wiley Interdiscip. Rev. RNA 7, 129–143. doi:10.1002/wrna.1321
Pal, A., and Kundu, R. (2019). 'Human Papillomavirus E6 and E7: The Cervical Cancer Hallmarks and Targets for Therapy. Front. Microbiol. 10, 3116. doi:10.3389/fmicb.2019.03116
Panda, A. C. (2018). 'Circular RNAs Act as miRNA Sponges. Adv. Exp. Med. Biol. 1087, 67–79. doi:10.1007/978-981-13-1426-1_6
Pentland, I., Campos-Leon, K., Cotic, M., Davies, K. J., Wood, C. D., Groves, I. J., et al. (2018). Disruption of CTCF-YY1-dependent Looping of the Human Papillomavirus Genome Activates Differentiation-Induced Viral Oncogene Transcription. PLoS Biol. 16, e2005752. doi:10.1371/journal.pbio.2005752
Piyathilake, C. J., Macaluso, M., Chambers, M. M., Badiga, S., Siddiqui, N. R., Bell, W. C., et al. (2014). 'Folate and Vitamin B12 May Play a Critical Role in Lowering the HPV 16 Methylation-Associated Risk of Developing Higher Grades of CIN. Cancer Prev. Res. (Phila) 7, 1128–1137. doi:10.1158/1940-6207.capr-14-0143
Pun, P. B., Liao, Y. P., Su, P. H., Wang, H. C., Chen, Y. C., Hsu, Y. W., et al. (2015). Triage of High-Risk Human Papillomavirus-Positive Women by Methylated POU4F3. Clin. Epigenetics 7, 85. doi:10.1186/s13148-015-0122-0
Pyeon, D., Newton, M. A., Lambert, P. F., den Boon, J. A., Sengupta, S., Marsit, C. J., et al. (2007). 'Fundamental Differences in Cell Cycle Deregulation in Human Papillomavirus-Positive and Human Papillomavirus-Negative Head/neck and Cervical Cancers. Cancer Res. 67, 4605–4619. doi:10.1158/0008-5472.can-06-3619
Qian, K., Pietila, T., Ronty, M., Michon, F., Frilander, M. J., Ritari, J., et al. (2013). Identification and Validation of Human Papillomavirus Encoded microRNAs. PLoS One 8, e70202. doi:10.1371/journal.pone.0070202
Qu, X., Li, Y., Wang, L., Yuan, N., Ma, M., and Chen, Y. (2020). 'LncRNA SNHG8 Accelerates Proliferation and Inhibits Apoptosis in HPV-Induced Cervical Cancer through Recruiting EZH2 to Epigenetically Silence RECK Expression. J. Cell Biochem. 121, 4120–4129. doi:10.1002/jcb.29646
Ramirez, A. T., Sanchez, G. I., Nedjai, B., Agudelo, M. C., Brentnall, A. R., Cuschieri, K., et al. Asc-Us-Col Trial Group (2021). 'Effective Methylation Triage of HPV Positive Women with Abnormal Cytology in a Middle-Income Country. Int. J. Cancer 148, 1383–1393. doi:10.1002/ijc.33314
Ren, L., Yang, S., Cao, Q., and Tian, J. (2021). 'CRNDE Contributes Cervical Cancer Progression by Regulating miR-4262/ZEB1 Axis. Onco Targets Ther. 14, 355–366. doi:10.2147/ott.s263505
Ren, S., Gaykalova, D. A., Guo, T., Favorov, A. V., Fertig, E. J., Tamayo, P., et al. (2020). 'HPV E2, E4, E5 Drive Alternative Carcinogenic Pathways in HPV Positive Cancers. Oncogene 39, 6327–6339. doi:10.1038/s41388-020-01431-8
Ren, S., Gaykalova, D., Wang, J., Guo, T., Danilova, L., Favorov, A., et al. (2018). 'Discovery and Development of Differentially Methylated Regions in Human Papillomavirus-Related Oropharyngeal Squamous Cell Carcinoma. Int. J. Cancer 143, 2425–2436. doi:10.1002/ijc.31778
Rinn, J. L., and Chang, H. Y. (2012). 'Genome Regulation by Long Noncoding RNAs. Annu. Rev. Biochem. 81, 145–166. doi:10.1146/annurev-biochem-051410-092902
Roman, A., and Munger, K. (2013). The Papillomavirus E7 Proteins. Virology 445, 138–168. doi:10.1016/j.virol.2013.04.013
Saha, S. S., Chowdhury, R. R., Mondal, N. R., Roy, S., and Sengupta, S. (2017). 'Expression Signatures of HOX Cluster Genes in Cervical Cancer Pathogenesis: Impact of Human Papillomavirus Type 16 Oncoprotein E7. Oncotarget 8, 36591–36602. doi:10.18632/oncotarget.16619
Sartor, M. A., Dolinoy, D. C., Jones, T. R., Colacino, J. A., Prince, M. E., Carey, T. E., et al. (2011). Genome-wide Methylation and Expression Differences in HPV(+) and HPV(-) Squamous Cell Carcinoma Cell Lines Are Consistent with Divergent Mechanisms of Carcinogenesis. Epigenetics 6, 777–787. doi:10.4161/epi.6.6.16216
Scarth, J. A., Patterson, M. R., Morgan, E. L., and Macdonald, A. (2021). The Human Papillomavirus Oncoproteins: a Review of the Host Pathways Targeted on the Road to Transformation. J. Gen. Virol. 102, 001540. doi:10.1099/jgv.0.001540
Schmitt, A. M., Garcia, J. T., Hung, T., Flynn, R. A., Shen, Y., Qu, K., et al. (2016). An Inducible Long Noncoding RNA Amplifies DNA Damage Signaling. Nat. Genet. 48, 1370–1376. doi:10.1038/ng.3673
Schmitz, M., Eichelkraut, K., Schmidt, D., Zeiser, I., Hilal, Z., Tettenborn, Z., et al. (2018). Performance of a DNA Methylation Marker Panel Using Liquid-Based Cervical Scrapes to Detect Cervical Cancer and its Precancerous Stages. BMC Cancer 18, 1197. doi:10.1186/s12885-018-5125-8
Schmitz, M., Wunsch, K., Hoyer, H., Scheungraber, C., Runnebaum, I. B., Hansel, A., et al. (2017). Performance of a Methylation Specific Real-Time PCR Assay as a Triage Test for HPV-Positive Women. Clin. Epigenetics 9, 118. doi:10.1186/s13148-017-0419-2
Schmitz, S. U., Grote, P., and Herrmann, B. G. (2016). Mechanisms of Long Noncoding RNA Function in Development and Disease. Cell Mol. Life Sci. 73, 2491–2509. doi:10.1007/s00018-016-2174-5
Sharma, S., and Munger, K. (2018). Expression of the Cervical Carcinoma Expressed PCNA Regulatory (CCEPR) Long Noncoding RNA Is Driven by the Human Papillomavirus E6 Protein and Modulates Cell Proliferation Independent of PCNA. Virology 518, 8–13. doi:10.1016/j.virol.2018.01.031
Sharma, S., and Munger, K. (2020a). Expression of the Long Noncoding RNA DINO in Human Papillomavirus-Positive Cervical Cancer Cells Reactivates the Dormant TP53 Tumor Suppressor through ATM/CHK2 Signaling. mBio 11, e01190–20. doi:10.1128/mBio.01190-20
Sharma, S., and Munger, K. (2020b). KDM6A-Mediated Expression of the Long Noncoding RNA DINO Causes TP53 Tumor Suppressor Stabilization in Human Papillomavirus 16 E7-Expressing Cells. J. Virol. 94, e02178–19. doi:10.1128/JVI.02178-19
Sharma, S., and Munger, K. (2020c). The Role of Long Noncoding RNAs in Human Papillomavirus-Associated Pathogenesis. Pathogens 9, 289. doi:10.3390/pathogens9040289
Shen-Gunther, J., Xia, Q., Cai, H., and Wang, Y. (2021). HPV DeepSeq: An Ultra-fast Method of NGS Data Analysis and Visualization Using Automated Workflows and a Customized Papillomavirus Database in CLC Genomics Workbench. Pathogens 10, 1026. doi:10.3390/pathogens10081026
Siegel, E. M., Ajidahun, A., Berglund, A., Guerrero, W., Eschrich, S., Putney, R. M., et al. (2021). 'Genome-wide Host Methylation Profiling of Anal and Cervical Carcinoma. PLoS One 16, e0260857. doi:10.1371/journal.pone.0260857
Soflaei, S. S., Momtazi-Borojeni, A. A., Majeed, M., Derosa, G., Maffioli, P., and Sahebkar, A. (2018). Curcumin: A Natural Pan-HDAC Inhibitor in Cancer. Curr. Pharm. Des. 24, 123–129. doi:10.2174/1381612823666171114165051
Song, P., and Yin, S. C. (2016). 'Long Non-coding RNA EWSAT1 Promotes Human Nasopharyngeal Carcinoma Cell Growth In Vitro by Targeting miR-326/-330-5p. Aging (Albany NY) 8, 2948–2960. doi:10.18632/aging.101103
St Laurent, G., Wahlestedt, C., and Kapranov, P. (2015). The Landscape of Long Noncoding RNA Classification. Trends Genet. 31, 239–251. doi:10.1016/j.tig.2015.03.007
Stunkel, W., and Bernard, H. U. (1999). 'The Chromatin Structure of the Long Control Region of Human Papillomavirus Type 16 Represses Viral Oncoprotein Expression. J. Virol. 73, 1918–1930. doi:10.1128/JVI.73.3.1918-1930.1999
Sun, L., Yang, C., Xu, J., Feng, Y., Wang, L., and Cui, T. (2016). 'Long Noncoding RNA EWSAT1 Promotes Osteosarcoma Cell Growth and Metastasis through Suppression of MEG3 Expression. DNA Cell Biol. 35, 812–818. doi:10.1089/dna.2016.3467
Sun, N. X., Ye, C., Zhao, Q., Zhang, Q., Xu, C., Wang, S. B., et al. (2014). 'Long Noncoding RNA-EBIC Promotes Tumor Cell Invasion by Binding to EZH2 and Repressing E-Cadherin in Cervical Cancer. PLoS One 9, e100340. doi:10.1371/journal.pone.0100340
Sundaram, M. K., Ansari, M. Z., Al Mutery, A., Ashraf, M., Nasab, R., Rai, S., et al. (2018). Genistein Induces Alterations of Epigenetic Modulatory Signatures in Human Cervical Cancer Cells. Anticancer Agents Med. Chem. 18, 412–421. doi:10.2174/1871520617666170918142114
Sundaram, M. K., Unni, S., Somvanshi, P., Bhardwaj, T., Mandal, R. K., Hussain, A., et al. (2019). Genistein Modulates Signaling Pathways and Targets Several Epigenetic Markers in HeLa Cells. Genes (Basel) 10, 955. doi:10.3390/genes10120955
Sung, H., Ferlay, J., Siegel, R. L., Laversanne, M., Soerjomataram, I., Jemal, A., et al. (2021). Global Cancer Statistics 2020: GLOBOCAN Estimates of Incidence and Mortality Worldwide for 36 Cancers in 185 Countries. CA Cancer J. Clin. 71 (3), 209–249. doi:10.3322/caac.21660
Tian, Y., Yuan Wu, N. Y., Liou, Y. L., Yeh, C. T., Cao, L., Kang, Y. N., et al. (2017). 'Utility of Gene Methylation Analysis, Cytological Examination, and HPV-16/18 Genotyping in Triage of High-Risk Human Papilloma Virus-Positive Women. Oncotarget 8, 62274–62285. doi:10.18632/oncotarget.19459
Tong, C., and Yin, Y. (2021). 'Localization of RNAs in the Nucleus: Cis- and Trans- Regulation. RNA Biol. 18, 2073–2086. doi:10.1080/15476286.2021.1894025
Torres-Rojas, F. I., Alarcon-Romero, L. D. C., Leyva-Vazquez, M. A., Ortiz-Ortiz, J., Mendoza-Catalan, M. A., Hernandez-Sotelo, D., et al. (2018). 'Methylation of the L1 Gene and Integration of Human Papillomavirus 16 and 18 in Cervical Carcinoma and Premalignant Lesions. Oncol. Lett. 15, 2278–2286. doi:10.3892/ol.2017.7596
Tumban, E. (2019). A Current Update on Human Papillomavirus-Associated Head and Neck Cancers. Viruses 11, 922. doi:10.3390/v11100922
Turan, T., Kalantari, M., Calleja-Macias, I. E., Cubie, H. A., Cuschieri, K., Villa, L. L., et al. (2006). Methylation of the Human Papillomavirus-18 L1 Gene: a Biomarker of Neoplastic Progression? Virology 349, 175–183. doi:10.1016/j.virol.2005.12.033
Van Keer, S., van Splunter, A. P., Pattyn, J., De Smet, A., Herzog, S. A., Van Ostade, X., et al. (2021). 'Triage of Human Papillomavirus Infected Women by Methylation Analysis in First-Void Urine. Sci. Rep. 11, 7862. doi:10.1038/s41598-021-87329-1
van Leeuwen, R. W., Ostrbenk, A., Poljak, M., van der Zee, A. G. J., Schuuring, E., and Wisman, G. B. A. (2019). 'DNA Methylation Markers as a Triage Test for Identification of Cervical Lesions in a High Risk Human Papillomavirus Positive Screening Cohort. Int. J. Cancer 144, 746–754. doi:10.1002/ijc.31897
Vande Pol, S. B., and Klingelhutz, A. J. (2013). 'Papillomavirus E6 Oncoproteins. Virology 445, 115–137. doi:10.1016/j.virol.2013.04.026
Vats, A., Trejo-Cerro, O., Thomas, M., and Banks, L. (2021). Human Papillomavirus E6 and E7: What Remains? Tumour Virus Res. 11, 200213. doi:10.1016/j.tvr.2021.200213
Verhoef, L., Bleeker, M. C. G., Polman, N., Steenbergen, R. D. M., Meijer, Cjlm., Melchers, W. J. G., et al. (2022). 'Performance of DNA Methylation Analysis of ASCL1, LHX8, ST6GALNAC5, GHSR, ZIC1 and SST for the Triage of HPV-Positive Women: Results from a Dutch Primary HPV-Based Screening Cohort. Int. J. Cancer 150, 440–449. doi:10.1002/ijc.33820
Verhoef, V. M., Heideman, D. A., van Kemenade, F. J., Rozendaal, L., Bosgraaf, R. P., Hesselink, A. T., et al. (2014). 'Methylation Marker Analysis and HPV16/18 Genotyping in High-Risk HPV Positive Self-Sampled Specimens to Identify Women with High Grade CIN or Cervical Cancer. Gynecol. Oncol. 135, 58–63. doi:10.1016/j.ygyno.2014.08.003
Verhoef, V. M., van Kemenade, F. J., Rozendaal, L., Heideman, D. A., Bosgraaf, R. P., Hesselink, A. T., et al. (2015). 'Follow-up of High-Risk HPV Positive Women by Combined Cytology and Bi-marker CADM1/MAL Methylation Analysis on Cervical Scrapes. Gynecol. Oncol. 137, 55–59. doi:10.1016/j.ygyno.2015.01.550
Verlaat, W., Snoek, B. C., Heideman, D. A. M., Wilting, S. M., Snijders, P. J. F., Novianti, P. W., et al. (2018). 'Identification and Validation of a 3-Gene Methylation Classifier for HPV-Based Cervical Screening on Self-Samples. Clin. Cancer Res. 24, 3456–3464. doi:10.1158/1078-0432.ccr-17-3615
Vink, F. J., Lissenberg-Witte, B. I., Meijer, Cjlm., Berkhof, J., van Kemenade, F. J., Siebers, A. G., et al. (2021). 'FAM19A4/miR124-2 Methylation Analysis as a Triage Test for HPV-Positive Women: Cross-Sectional and Longitudinal Data from a Dutch Screening Cohort. Clin. Microbiol. Infect. 27, 125–25. doi:10.1016/j.cmi.2020.03.018
Vink, F. J., Meijer, C. J. L. M., Clifford, G. M., Poljak, M., Oštrbenk, A., and Petry, K. U. (2019). FAM19A4/miR124‐2 Methylation in Invasive Cervical Cancer: A Retrospective Cross‐sectional Worldwide Study. Int. J. Cancer 147, 1215–1221. doi:10.1002/ijc.32614
Vinokurova, S., and von Knebel Doeberitz, M. (2011). 'Differential Methylation of the HPV 16 Upstream Regulatory Region during Epithelial Differentiation and Neoplastic Transformation. PLoS One 6, e24451. doi:10.1371/journal.pone.0024451
Virtanen, E., Pietila, T., Nieminen, P., Qian, K., and Auvinen, E. (2016). 'Low Expression Levels of Putative HPV Encoded microRNAs in Cervical Samples. Springerplus 5, 1856. doi:10.1186/s40064-016-3524-3
Wang, J. W., and Roden, R. B. (2013). 'L2, the Minor Capsid Protein of Papillomavirus. Virology 445, 175–186. doi:10.1016/j.virol.2013.04.017
Wang, J., Xiang, F., Liu, X., Ma, X., Cai, X., Yang, Y., et al. (2021). HPV E7 Affects the Function of Cervical Cancer Cells via the TAL1/lncEBIC/KLHDC7B axis. Oncol. Rep. 45, 51. doi:10.3892/or.2021.8002
Wang, Y., Tong, J., Lin, H., Ma, L., and Xu, Y. (2020). CCHE1 Accelerated the Initiation of Oral Squamous Cell Carcinoma through Enhancing PAK2 Expression by Sponging miR-922. J. Oral Pathol. Med. 49, 636–644. doi:10.1111/jop.12995
Weng, S. L., Huang, K. Y., Weng, J. T., Hung, F. Y., Chang, T. H., and Lee, T. Y. (2018). Genome-wide Discovery of Viral microRNAs Based on Phylogenetic Analysis and Structural Evolution of Various Human Papillomavirus Subtypes. Brief. Bioinform 19, 1102–1114. doi:10.1093/bib/bbx046
Wilson, R., Fehrmann, F., and Laimins, L. A. (2005). 'Role of the E1--E4 Protein in the Differentiation-dependent Life Cycle of Human Papillomavirus Type 31. J. Virol. 79, 6732–6740. doi:10.1128/jvi.79.11.6732-6740.2005
Wooldridge, T. R., and Laimins, L. A. (2008). 'Regulation of Human Papillomavirus Type 31 Gene Expression during the Differentiation-dependent Life Cycle through Histone Modifications and Transcription Factor Binding. Virology 374, 371–380. doi:10.1016/j.virol.2007.12.011
Xu, H., Ling, Y., Xi, Y., Ma, H., Wang, H., Hu, H. M., et al. (2019). Viral Metagenomics Updated the Prevalence of Human Papillomavirus Types in Anogenital Warts. Emerg. Microbes Infect. 8, 1291–1299. doi:10.1080/22221751.2019.1661757
Xu, W., Xu, M., Wang, L., Zhou, W., Xiang, R., Shi, Y., et al. (2019). 'Integrative Analysis of DNA Methylation and Gene Expression Identified Cervical Cancer-specific Diagnostic Biomarkers. Signal Transduct. Target Ther. 4, 55. doi:10.1038/s41392-019-0081-6
Xu, X., Xie, Q., Xie, M., Zeng, Y., and Liu, Q. (2021). 'LncRNA SNHG8 Serves as an Oncogene in Breast Cancer through miR-634/ZBTB20 Axis. Cancer Manag. Res. 13, 3017–3028. doi:10.2147/cmar.s270128
Yang, L., Yi, K., Wang, H., Zhao, Y., and Xi, M. (2016). 'Comprehensive Analysis of lncRNAs Microarray Profile and mRNA-lncRNA Co-expression in Oncogenic HPV-Positive Cervical Cancer Cell Lines. Oncotarget 7, 49917–49929. doi:10.18632/oncotarget.10232
Yuan, X. G., Huang, Y. R., Yu, T., Jiang, H. W., Xu, Y., and Zhao, X. Y. (2019). 'Chidamide, a Histone Deacetylase Inhibitor, Induces Growth Arrest and Apoptosis in Multiple Myeloma Cells in a Caspase-dependent Manner. Oncol. Lett. 18, 411–419. doi:10.3892/ol.2019.10301
Zhang, J. J., and Fan, L. P. (2019). 'Long Non-coding RNA CRNDE Enhances Cervical Cancer Progression by Suppressing PUMA Expression. Biomed. Pharmacother. 117, 108726. doi:10.1016/j.biopha.2019.108726
Zhang, Q., Padi, S. K., Tindall, D. J., and Guo, B. (2014). 'Polycomb Protein EZH2 Suppresses Apoptosis by Silencing the Proapoptotic miR-31. Cell Death Dis. 5, e1486. doi:10.1038/cddis.2014.454
Zhao, J., Lee, E. E., Kim, J., Yang, R., Chamseddin, B., Ni, C., et al. (2019). 'Transforming Activity of an Oncoprotein-Encoding Circular RNA from Human Papillomavirus. Nat. Commun. 10, 2300. doi:10.1038/s41467-019-10246-5
Zhong, M., Tan, J., Pan, G., Jiang, Y., Zhou, H., Lai, Q., et al. (2021). 'Preclinical Evaluation of the HDAC Inhibitor Chidamide in Transformed Follicular Lymphoma. Front. Oncol. 11, 780118. doi:10.3389/fonc.2021.780118
Keywords: HPV, epigenetic, therapy, viral oncoprotein, regulation, biomarker, miRNA, LncRNA
Citation: Castro-Oropeza R and Piña-Sánchez P (2022) Epigenetic and Transcriptomic Regulation Landscape in HPV+ Cancers: Biological and Clinical Implications. Front. Genet. 13:886613. doi: 10.3389/fgene.2022.886613
Received: 28 February 2022; Accepted: 24 May 2022;
Published: 14 June 2022.
Edited by:
Ángeles Carlos Reyes, Instituto Nacional de Enfermedades Respiratorias(INER), MexicoReviewed by:
Pietro Taverna, Constellation Pharmaceuticals, United StatesStephanie I. Nuñez-Olvera, National Autonomous University of Mexico, Mexico
Copyright © 2022 Castro-Oropeza and Piña-Sánchez. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Patricia Piña-Sánchez, patricia_1307@yahoo.com.mx