Megha Budhwani
Samuel W. Lukowski
Sandro V. Porceddu
Ian H. Frazer
Janin Chandra- 1Diamantina Institute, Translational Research Institute, The University of Queensland, Woolloongabba, QLD, Australia
- 2Institute for Molecular Bioscience, The University of Queensland, St Lucia, QLD, Australia
- 3Cancer Services, Princess Alexandra Hospital, Woolloongabba, QLD, Australia
- 4Faculty of Medicine, The University of Queensland, Brisbane, QLD, Australia
Human papillomavirus (HPV) infection is associated with a range of malignancies that affect anogenital and oropharyngeal sites. α-HPVs dominantly infect basal epithelial cells of mucosal tissues, where they dysregulate cell division and local immunity. The cervix is one of the mucosal sites most susceptible to HPV infections. It consists of anatomically diverse regions, and the majority of cervical intraepithelial neoplasia and cancers arise within the cervical squamo-columnar junction where undifferentiated basal progenitor cells with stem cell properties are found. The cancer stem cell theory particularly associates tumorigenesis, invasion, dissemination, and metastasis with cancer cells exhibiting stem cell properties. In this perspective, we discuss evidence of a cervical cancer stem cell niche and explore the association of stemness related genes with 5-year survival using a publicly available transcriptomic dataset of a cervical cancer cohort. We report that poor prognosis in this cohort correlates with overexpression of a subset of stemness pathway genes, a majority of which regulate the central Focal Adhesion pathway, and are also found to be enriched in the HPV infection pathway. These observations support therapeutic targeting of stemness genes overexpressed by mucosal cells infected with high-risk HPVs.
Introduction to HPV Virology and Disease
Oncogenic human papillomavirus (HPV) infections are responsible for 5% of the global human cancer burden, resulting in over 634,000 new cancers annually, of which ~90% are in women (Schiffman et al., 2016). Intraepithelial neoplasia and cancers driven by oncogenic HPV infections include those of the cervix, vagina, vulva, penis, and anus, and of the oropharynx where 570,000 new cases of cervical cancer alone and associated 311,000 female deaths were reported in 2018 (Arbyn et al., 2020). HPV infection is also associated with a proportion of bladder cancer (Zapatka et al., 2020). Five genera of HPVs (alpha, beta, gamma, mu, nu) are classified, comprising more than 450 HPV genotypes, but only so called high-risk α-HPVs (HPV16, 18, 31, 33, 35, 39, 45, 51, 52, 56, 58, 59, and 66) are clearly associated with malignancy, and have been classified as group I carcinogens (Burley et al., 2020). α-HPV types 16 and 18 account for 70% of cervical cancers (Tommasino, 2014) while low-risk genotypes are identified in mucosal and cutaneous papillomas (warts) and in apparently normal skin (Devitt et al., 2019) and genital epithelium (Egawa et al., 2015).
HPV16 is an epitheliotropic 60 nm small, non-enveloped, 8-kb circular double-stranded DNA virus and infects basal keratinocytes primarily of mucosal epithelium, where viral DNA is maintained within the nucleus as episomal DNA (Zhou C. et al., 2019). HPV genes include early (E1–E7) non-structural and late (L1-2) structural genes. In proliferating basal cells, the viral DNA is amplified to 50–100 copies, while viral gene expression is limited to low levels of E1 and E2 (Moody and Laimins, 2010). Viral early and late gene expression, however, increases in differentiating keratinocytes. As HPV uses the host's DNA replication machinery to replicate its episome, it utilizes strategies to delay basal keratinocyte differentiation and exit from cell cycle. HPV proteins E6 and E7 prolong keratinocyte survival and proliferation by initiating p53 degradation and continued E2F release (Yim and Park, 2005; Yeo-Teh et al., 2018), thereby enabling HPV episomal replication. L1 and L2 capsid proteins are expressed in terminally differentiated cells, where they assemble into functional viral particles and are shed from the surface (Egawa et al., 2015). Malignant transformation of HPV infected cells usually occurs over many years, and is associated with random integration of viral DNA into the host genome and specific gene mutations (Zapatka et al., 2020).
The Microanatomy of the Cervix Determines Prevalence of Premalignant Transformation
The cervix can be anatomically divided into the columnar endocervix and the squamous ectocervix, which are connected by the squamo-columnar junction (SCJ) within the transformation zone (Herfs et al., 2017; Huang and Rofstad, 2017). The SCJ is a monolayer consisting of two types of specialized epithelial progenitor cells, reserve cells, and cuboidal cells, with stem-like properties, which as undifferentiated basal cells regenerate the endo- and ectocervix (Figure 1A) (Schiffman et al., 2016). Productive α-HPV infection occurs across the genital mucosa, but malignant transformation and cancer is more common in the cervix, and specifically in the SCJ (Herfs et al., 2011, 2012; Doorbar and Griffin, 2019).
Figure 1. (A) While the productive life cycle of HPV occurs in the squamous ectocervix, cervical intraepithelial neoplasia (CIN) and invasive cervical cancer dominantly arise from the squamo-columnar junction (SJL) of the cervix, which contains reserve and cuboidal cells with stem-like properties. Stemness-related pathways such as Focal Adhesion, Notch, Wnt, Hippo, and Hedgehog with cross-talk mechanisms facilitate cell survival, proliferation, epithelial mesenchymal transition, dissemination, and metastasis. (B) Using the Xenabrowser online portal (https://xenabrowser.net/) (Goldman et al., 2019), The Cancer Genome Atlas Cervical Squamous Cell Carcinoma and Endocervical Adenocarcinoma (TCGA-CESC) database was filtered on primary site samples (n = 304) and screened for significant association of high expression of stemness-related genes with poor 5-year (1,825 days) survival using Kaplan-Meier Plots, comparing subjects of the top and bottom 25% of gene expression. The web-based interactive Xenabrowser analysis integrates data of the following datasets: https://tcga.xenahubs.net/download/TCGA.CESC.sampleMap/HiSeqV2.gz (gene expression data), https://tcga.xenahubs.net/download/survival/CESC_survival.txt.gz (survival data) and https://tcga.xenahubs.net/download/TCGA.CESC.sampleMap/CESC_clinicalMatrix (clinical data). Input genes (n = 526) were selected from KEGG pathways (hsa04150, hsa04330, hsa04310, hsa04390, hsa04340). P-values of significantly (p < 0.05) associated genes where overexpression associated with poorer survival were –log10 transformed and plotted. Underlined genes are those not previously associated with cervical cancer prognosis. (C) Significant genes from (B) were analyzed for enriched cancer and viral infection pathways using Enrichr and the KEGG 2019 human database. P-values were –log10 transformed and plotted. (D) STRING-db analysis of genes from (B) was used to estimate protein-protein interaction networks and levels of connectivity. Genes enriched in the KEGG HPV infection pathway (hsa05165) are displayed in red. (E) Heat map of genes ordered based on highest to lowest connectivity estimated by STRING-db. (F) A combined score was established based on the sum of a connectivity score (1: 0–10 connections; 2: 11–20 connections; 3: 21–30 connections) and a p-value score (1: –log10 < 1.5; 2: –log10 = 1.5–2.0; 3: –log10 > 2.0), and displayed as heat map.
In contrast, the SCJ/transformation zones of vulva, vagina, and anus are differentiated and multi-layered, and the cancer ratio of cervical compared to vulvar/vaginal/anal is ~20:1 while the frequency of HPV infection in cervix and anus is similar (Yang et al., 2015), indicating that HPV infection to the cervical SCJ is more likely to result in malignant transformation compared to other anogenital sites. In the oropharynx, the majority of HPV16 associated cancers arise in the tonsillar crypts and base of tongue (Westra, 2012; Morbini et al., 2015; Doorbar and Griffin, 2019). The specialized porous reticulated epithelium of the tonsillar crypts appears to be more susceptible to carcinogenesis due to the presence of a basal cell type with reduced polarization which results in an increased likelihood for persistent infection, a main risk factor for malignancy (Tomaic, 2016).
Cervical SCJ Cells Represent a Privileged Site for HPV Infection
Increased vulnerability of SCJ cells to HPV infection and persistence is likely caused by multiple factors. As the cervical SCJ consists of a monolayer of cells, no physical or chemical trauma is required to acquire access to basal cells for infection. Additionally, due to monthly pH fluctuations associated with menstrual hormonal changes, the cervical mucosa undergoes regular changes in the extracellular environment. Also, the cervical SCJ contains an altered composition of host-defense peptides. Alpha defensins have been shown to inhibit HPV infection (Buck et al., 2006). However, SCJ cells and neoplastic lesions lack expression of the human defensin 5 (HD5) (Hubert et al., 2014), and while HD5 doesn't block HPV entry into cells, it redirects HPV viral particles to the lysosome for degradation, therewith counteracting successful infection (Wiens and Smith, 2017). Tetraspanin has been implicated in post-endocytic trafficking of HPVs by interacting with HPV L1, enabling capsid disassembly (Grassel et al., 2016). Interestingly, tetraspanin is overexpressed in SCJ cells (Herfs et al., 2012), possibly enhancing HPV reproduction. Two cytokeratins, Krt7 and Krt9, are also overexpressed in cervical and anal SCJs, and have been shown to enhance translation of HPV transcripts (Kanduc, 2002; Favia et al., 2004). These mechanisms might collectively help establish infection of HPV within SCJ cells more effectively.
In addition to the above, the SCJ favors HPV persistence as it is an immune-privileged site. The density of professional antigen-presenting cells (APCs) in the cervical SCJ and in HPV associated neoplastic lesions is decreased, and this might be related to the lack of HD5, which acts as a chemoattractant. Furthermore, elevated TGF-β levels inhibit E-cadherin expression, which disrupts APC adhesion and subsequently alters APC maturation and regulator T cell polarization (Jiang et al., 2007; Herfs et al., 2008). Prostaglandin E2 is overexpressed in the cervical SCJ and in neoplastic lesions, and inhibits APC migration. Prostaglandin E2 also promotes immune tolerance by stimulating IL-10 and inhibiting IL-12 production (Herfs et al., 2009). RANKL is also overexpressed in SCJ and promotes development of tolerogenic APCs (Demoulin et al., 2015). Also the density of plasmacytoid DCs (pDCs) and regulatory T cells increases from the ectocervix to metaplasia and (pre)neoplasia, and pDCs exposed to HPV-transformed keratinocyte cell lines were inhibited in maturation and IFN-α secretion, and induced regulatory T cells (Demoulin S. et al., 2015). In a transgenic mouse model expressing HPV16 E7 under the control of the keratinocyte-specific K14 promotor and which displays epithelial hyperplasia and shares a gene signature with human cervical intraepithelial neoplasia (CIN) (Tuong et al., 2018), we observed multiple immune-suppressive mechanisms including alteration of antigen-presenting cells and induction of regulatory T cells (Mittal et al., 2013; Chandra et al., 2016; Bashaw et al., 2017, 2019). Together, these observations might explain why the SCJ epithelium expressing viral oncoproteins is more likely than elsewhere in the genital tract to establish persistent HPV infection and malignant transformation.
Evidence for a Cervical Cancer Stem Cell Niche
The morphology of cervical neoplastic lesions is heterogenous. Well-differentiated squamous lesions tend to arise towards the ectocervix, while lesions with an immature phenotype locate towards the endocervix. However, the unique morphology and genetic signature of cervical SCJ cells is evident in cervical squamous, columnar, and adenosquamous neoplastic lesions, suggesting that most cervical malignancies arise from multipotent cervical SCJ cells (Herfs et al., 2012, 2017; Sato et al., 2017). Further, Krt7+ CIN stage 1 located in the SCJ has a higher risk to progress to cancer than CIN stage 1 located in the ectocervix (Paquette et al., 2016; Mills et al., 2017). SCJ cells have, therefore, been proposed as cervical cancer stem cells (CCSCs). CCSCs can be isolated based on the expression of cytokeratins and CDs and are resistant to cisplatin-chemotherapy and radiotherapy (Mendoza-Almanza et al., 2019). As discussed earlier, cytokeratin expression is increased in CIN and cervical cancer and is associated with metastatic processes (Wang et al., 2006; Ikeda et al., 2008). Biomarkers currently accepted for classification of CCSCs include CD44, CD133, CD49F, ABCG2, ALDH1, OCT4, OPN, and SOX2 (Huang and Rofstad, 2017). Additionally, high level expression of stem cell markers MSI1, ALDH1, and SOX2 in cervical cancer lesions correlates with poor survival (Hou et al., 2015).
Cancer cells inappropriately maintain self-renewal pathways similar to those found in stem cells (SCs), allowing indefinite proliferation, and understanding how these self-renewal pathways are maintained in cancer has therapeutic implications. Cancer stem cells (CSCs) are found in many malignant diseases and are associated with poor survival (Zhu et al., 2020). CSCs tend to remain quiescent, thereby conferring protection from therapies which target mitotic cells, as well as escaping anti-cancer immune responses. Epithelial-mesenchymal transition (EMT) of CSCs is strongly linked with cancer cell metastasis, in which the proportion of CSCs is greatly elevated compared to primary tumors (Aktas et al., 2009).
To enhance our understanding of stemness as a contributor to prognosis of cervical cancer, here we systematically sought association of the expression of genes from key stemness pathways with cervical cancer survival and found that expression of multiple genes contributing to the central Focal Adhesion pathway correlated with poor survival. In the remainder of this perspective, we will discuss these stemness-associated genes and their potential role in HPV infection and cervical cancer progression.
Dysregulated Stemness-Associated Gene Expression in Primary Cervical Tumors
CCSCs share signaling pathways with SCs, including the Hedgehog, Notch, Wnt, and Hippo pathways that together control organ size, cell proliferation, adhesion, and apoptosis. Many of these pathways interact with the Focal Adhesion pathway that regulates tissue homeostasis and is involved in malignant transformation and metastasis through disrupted control of proliferation, differentiation, apoptosis, and migration (Figure 1A). These interactions place this pathway as a significant target for cancer therapy (Zhou J. et al., 2019). We examined the correlation between stemness-associated genes of key stemness KEGG pathways and the 5-year survival of cervical cancer patients using the Xenabrowser online data analysis tool (Goldman et al., 2019) filtered on The Cancer Genome Atlas Cervical Squamous Cell Carcinoma and Endocervical Adenocarcinoma (TCGA-CESC) dataset of primary tumors (n = 304), of which 93% have previously been determined to be HPV positive (Lu et al., 2019). Of 526 analyzed genes, we identified 45 genes where the top quartile of expression level, when compared with the bottom quartile of expression, was associated with poorer 5-year survival (Figure 1B). We also identified 31 genes where the bottom quartile of expression, when compared with the top quartile of expression, correlated with poor survival (data not shown).
To explore potential candidates for targeted therapies, we focused on and henceforth discuss the genes where overexpression correlated with poor survival. Overall, the majority of these genes belonged to the Focal Adhesion pathway, but we also identified genes of other stemness-related pathways including Notch, Hippo, Wnt, and Hedgehog, and many of these gene products signal across pathways (Figure 1B). While overexpression of 18 of these genes has been associated with cervical cancer previously (referenced in Table 1), we also identified 27 new prognosis-associated genes (underscored in Figure 1B). We applied a KEGG pathway enrichment analysis using Enrichr, and observed that the 45 genes were significantly overexpressed in multiple cancers including cancers of lung, stomach, liver, pancreas, breast, colorectum, and bladder (Figure 1C). Enrichr pathway analysis further revealed significant gene expression enrichment in the HPV infection pathway (Figures 1C,D) and other viral infection pathways including HTLV-1, Hepatitis B, and EBV (Figure 1C). We assessed protein-protein interactions for the 45 genes using STRING-db which estimated the number of connections of each gene (Figure 1E). We scored genes according to the number of predicted connections and p-value, where genes with a high level of connectivity and p-value received an overall high score (Figure 1F), as this can potentially guide identification of biomarker or therapeutic targets.
Table 1. Known and novel stemness-associated genes which associate with poor 5-year survival in cervical cancer.
The focal adhesion kinase PTK2, a gene in the core of the Focal Adhesion pathway, and its central activator fibronectin 1 (FN1), as well as upstream integrin subunits ITGB1, ITGB3, and ITGA5, paxillin (PXN), and thrombospondin 1 (THBS1) represented hub nodes next to the common cancer associates VEGFA and EGFR, and all these genes are part of the Focal Adhesion pathway. A range of PTK2 inhibitors are being tested in clinical cancer trials, with 40 studies currently registered (clinicaltrials.gov, 2020), including one study in patients with head and neck cancer. Integrins have been classified as metastasis-promoting genes due to their function in regulating adhesion, proliferation, and migration through cell-cell and cell-matrix interactions (Cooper and Giancotti, 2019). Paxillin contributes to dissemination and metastasis by accumulating in nascent focal complexes in migrating cells (Maziveyi and Alahari, 2017), and thrombospondin 1 regulates cell-matrix interactions by binding to integrins, fibronectins, and laminins, of which we identified LAMA1, LAMB1, LAMC1, and LAMC2 as associated with poor survival. We also identified 3 collagens (COL4A1, COL4A6, and COL9A3) associated with poor survival, and collagens are known to drive tumor budding (Miyake et al., 2017). Because integrins, paxillin, thrombospondin, laminins, and collagens are positioned upstream of multiple signaling cascades which contribute to cancer invasion and metastasis, these represent desired targets to inhibit these pathways. Other genes associated with poor survival and a high connectivity/p-value score included SHC1, ACTN1, PDGFB, ADAM17, DOCK1, LATS1, and NCOR1, and overexpression of many of these genes has been described in association with other cancers.
Interactions Between HPV Viral Proteins and Stemness Associated Proteins
Expression of some of the genes identified above has impact on the HPV viral life cycle. According to our current understanding of α-HPVs, viral gene expression can modulate cellular pathways to either impart basal stem-cell like characteristics to infected cells in the epithelium or infect specialized epithelial stem-cell populations. In any event, the virus interacts with and affects genes and proteins involved in stemness regulation and maintenance. HPV viral gene expression is highly regulated with limited expression in the infected basal epithelial layer where early protein E2 facilitates viral episomal genome partitioning during the division of infected cells. E2 also regulates cell differentiation and has been shown to interact with focal adhesion kinase PTK2B and FN1 (Muller and Demeret, 2012), two of the top scoring stemness pathway genes we identified to be associated with cervical cancer prognosis. Both of these proteins are involved in cellular adhesion, which suggests an early impact of the viral life cycle not only on stemness but also on cohesion of the epithelium perhaps contributing to invasiveness.
During the phase of low HPV viral gene expression, E6 also plays an important role by helping delay the migration of infected epithelial basal cells to the suprabasal layers. Viral genome amplification occurs in the middle epithelial layers, where deregulated expression of E6 and E7 is seen in association with the manifestation of high-grade neoplasia. In addition to suppressing host immune responses, and associating with key cell cycle proteins p53 and Rb, the early HPV viral proteins also bind to other cellular proteins, enabling suprabasal cell cycle entry. As observed in the epithelial basal layers, levels of FN1 and other adhesion proteins are also increased by E6 and E7 expression in suprabasal layers (Hellner et al., 2009). Additionally, E6 expression can interfere with paxillin, affecting the integrity of the epithelium and leading to invasiveness (Estevao et al., 2019). Expression of integrins such as ITGB3 is also upregulated by HPV E6 (Yim et al., 2004), however, the mechanism is yet to be explored. Angiogenic protein VEGF, another top-scoring gene in our analysis, is involved in increasing vascular permeability and is also upregulated by E6 and E7 (Lopez-Ocejo et al., 2000; Li and Cui, 2015). Interestingly, E6 has been further shown to cause degradation of the polarity protein SCRIB (Gheit, 2019), and overexpression of SCRIB is associated with poor cervical cancer prognosis. This observation is not completely surprising as SCRIB overexpression is seen in a number of cancers, with mislocalisation usually favoring a mesenchymal phenotype (Vaira et al., 2011), and perhaps SCRIB overexpression begins in later stages of HPV associated malignant disease.
HPV early proteins also influence cell proliferation. Among the candidate genes identified in Figure 1, YAP1 enhances proliferation and E6 prevents its proteasome-mediated degradation. YAP1 expression in turn promotes EGFR expression and has recently been shown to be essential for HPV-mediated carcinogenesis (He et al., 2015). In addition to indirect regulation of EGFR expression by E6, EGFR expression is also directly controlled by the early viral protein E5 and promotes genome amplification in mid-epithelial layers (Ilahi and Bhatti, 2020).
E5, E6, and E7 have previously been shown to co-operate to drive tumorigenesis. Co-expression of E5 with either E6 or E7 promotes transformation to a greater extent than with either oncoprotein alone (Ilahi and Bhatti, 2020). The association of these viral proteins with stemness-associated genes suggests a potential for a combined targeted therapy approach against cervical cancers.
Conclusion
The biggest challenges in the treatment of HPV-mediated cervical malignancy are associated with acquired local immune-suppression, and disruption of cellular homeostatic processes leading to impaired focal adhesion and break-down of extracellular matrix, promoting cancer cell invasion, dissemination, and metastasis. The cancer stem cell theory suggests that these processes are more likely to occur when pluripotent cells with active stemness signaling pathways undergo transformation. HPV infects mucosal epithelial cells at many different anatomical sites, including distinct sub-anatomic sections of the cervix, but malignant transformation occurs predominantly at the SCJ, which contains unique basal cells with stem-like properties, whereas the productive HPV life cycle occurs mainly in the squamous section of the cervix where no undifferentiated cells are seen. The cervical SCJ/transformation zone microanatomy is unique compared to those of other HPV-infected anogenital regions including vagina, vulva, and anus, and this specialized undifferentiated cervical SCJ is associated with a higher risk for malignant transformation. In our analysis of stemness-pathway genes in association with cervical cancer prognosis, we confirmed 18 genes known to be dysregulated in cervical cancer cells. We also identified 27 novel genes whose overexpression significantly associated with poor survival and a majority of these play a role in the central Focal Adhesion pathway. Furthermore, some of these genes were also associated with the HPV infection pathway, and viral protein interactions with stemness proteins have been demonstrated previously.
Importantly, the data presented here indicates but does not prove a direct contribution of the cervical stem cell niche to cervical cancer survival outcomes, as it was collected from bulk RNA sequencing data of primary tumor biopsies. An alternative interpretation would be that cervical cancers which overexpress genes associated with the Focal Adhesion pathway have further progressed on a trajectory towards invasion and dissemination, leading to poorer survival outcomes. Nevertheless, the pathway overlap with HPV infection confirms importance of stemness associated genes in HPV driven cervical cancer progression and might prove important for cancer therapy target selection. While all drivers of HPV infection and malignant transformation are still to be determined, treatment of cervical, and other HPV-associated cancers may be improved by a holistic approach targeting multiple cancerogenic pathways, as well as by initiating a strong HPV-specific immune response through immunization-based immunotherapy that enhances virus-specific immunity (Frazer and Chandra, 2019).
Data Availability Statement
Publicly available TCGA-CESC datasets were analyzed in this study. This data can be found here: https://xenabrowser.net/.
Ethics Statement
Ethical review and approval was not required for the study on human participants in accordance with the local legislation and institutional requirements. Written informed consent for participation was not required for this study in accordance with the national legislation and the institutional requirements.
Author Contributions
MB and JC conceived the work, analyzed data, and wrote the first draft of the manuscript. SL supervised data analysis and revised the manuscript. SP and IF wrote and revised the manuscript. All authors contributed to the article and approved the submitted version.
Funding
This work was funded in part by grants GNT1071822, 1173927, and 1051092 from the National Health and Medical Research Council of Australia.
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.
References
Aktas, B., Tewes, M., Fehm, T., Hauch, S., Kimmig, R., and Kasimir-Bauer, S. (2009). Stem cell and epithelial-mesenchymal transition markers are frequently overexpressed in circulating tumor cells of metastatic breast cancer patients. Breast Cancer Res. 11:R46. doi: 10.1186/bcr2333
Arbyn, M., Weiderpass, E., Bruni, L., De Sanjose, S., Saraiya, M., Ferlay, J., et al. (2020). Estimates of incidence and mortality of cervical cancer in 2018: a worldwide analysis. Lancet Glob Health 8, e191–e203. doi: 10.1016/S2214-109X(19)30482-6
Bashaw, A. A., Leggatt, G. R., Chandra, J., Tuong, Z. K., and Fraze,R, I. H. (2017). Modulation of antigen presenting cell functions during chronic HPV infection. Papillomavirus Res. 4, 58–65. doi: 10.1016/j.pvr.2017.08.002
Bashaw, A. A., Teoh, S. M., Tuong, Z. K., Leggatt, G. R., Frazer, I. H., and Chandra, J. (2019). HPV16 E7-driven epithelial hyperplasia promotes impaired antigen presentation and regulatory T cell development. J. Invest. Dermatol. 139, 2467–2476.e3. doi: 10.1016/j.jid.2019.03.1162
Buck, C. B., Day, P. M., Thompson, C. D., Lubkowski, J., Lu, W., Lowy, D. R., et al. (2006). Human alpha-defensins block papillomavirus infection. Proc. Natl. Acad. Sci. U.S.A. 103, 1516–1521. doi: 10.1073/pnas.0508033103
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
Campos-Parra, A. D., Padua-Bracho, A., Pedroza-Torres, A., Figueroa-Gonzalez, G., Fernandez-Retana, J., millan-Catalan, O., et al. (2016). Comprehensive transcriptome analysis identifies pathways with therapeutic potential in locally advanced cervical cancer. Gynecol. Oncol. 143, 406–413. doi: 10.1016/j.ygyno.2016.08.327
Cancer Genome Atlas Research Network Albert Einstein College of Medicine, Analytical Biological Services, Barretos Cancer Hospital, Baylor College of Medicine, Beckman Research Institute of City of Hope. (2017). Integrated genomic and molecular characterization of cervical cancer. Nature 543, 378–384. doi: 10.1038/nature21386
Chandra, J., Miao, Y., Romoff, N., and Frazer, I. H. (2016). Epithelium expressing the E7 oncoprotein of HPV16 attracts immune-modulatory dendritic cells to the skin and suppresses their antigen-processing capacity. PLoS ONE 11:e0152886. doi: 10.1371/journal.pone.0152886
Cheng, Y., Ma, D., Zhang, Y., Li, Z., and Geng, L. (2015). Cervical squamous cancer mRNA profiles reveal the key genes of metastasis and invasion. Eur. J. Gynaecol. Oncol. 36, 309–317. doi: 10.12892/ejgo2683.2015
Cooper, J., and Giancotti, F. G. (2019). Integrin signaling in cancer: mechanotransduction, stemness, epithelial plasticity, and therapeutic resistance. Cancer Cell 35, 347–367. doi: 10.1016/j.ccell.2019.01.007
Demoulin, S., Herfs, M., Somja, J., Roncarati, P., Delvenne, P., and Hubert, P. (2015). HMGB1 secretion during cervical carcinogenesis promotes the acquisition of a tolerogenic functionality by plasmacytoid dendritic cells. Int. J. Cancer 137, 345–358. doi: 10.1002/ijc.29389
Demoulin, S. A., Somja, J., Duray, A., Guenin, S., Roncarati, P., Delvenne, P. O., et al. (2015). Cervical (pre)neoplastic microenvironment promotes the emergence of tolerogenic dendritic cells via RANKL secretion. Oncoimmunology 4:e1008334. doi: 10.1080/2162402X.2015.1008334
Devitt, K., Hanson, S. J., Tuong, Z. K., Mcmeniman, E., Soyer, H. P., Frazer, I. H., et al. (2019). Single-cell RNA sequencing reveals cell type-specific HPV expression in hyperplastic skin lesions. Virology 537, 14–19. doi: 10.1016/j.virol.2019.08.007
Doorbar, J., and Griffin, H. (2019). Refining our understanding of cervical neoplasia and its cellular origins. Papillomavirus Res. 7, 176–179. doi: 10.1016/j.pvr.2019.04.005
Egawa, N., Egawa, K., Griffin, H., and Doorbar, J. (2015). Human papillomaviruses; epithelial tropisms, and the development of neoplasia. Viruses 7, 3863–3890. doi: 10.3390/v7072802
Estevao, 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. Biochim. Biophys. Acta Gene. Regul. Mech. 1862, 153–162. doi: 10.1016/j.bbagrm.2019.01.001
Fang, F., Huang, B., Sun, S., Xiao, M., Guo, J., Yi, X., et al. (2018). miR-27a inhibits cervical adenocarcinoma progression by downregulating the TGF-betaRI signaling pathway. Cell. Death Dis. 9:395. doi: 10.1038/s41419-018-0431-2
Favia, G., Kanduc, D., Lo Muzio, L., Lucchese, A., and Serpico, R. (2004). Possible association between HPV16 E7 protein level and cytokeratin 19. Int. J. Cancer 111, 795–797. doi: 10.1002/ijc.20343
Frazer, I. H., and Chandra, J. (2019). Immunotherapy for HPV associated cancer. Papillomavirus Res. 8:100176. doi: 10.1016/j.pvr.2019.100176
Gheit, T. (2019). Mucosal and cutaneous human papillomavirus infections and cancer biology. Front. Oncol. 9:355. doi: 10.3389/fonc.2019.00355
Goldman, M., Craft, B., Hastie, M., and Haussler, D. (2019). The UCSC Xena platform for public and private cancer genomics data visualization and interpretation. bioRvix [preprint]. doi: 10.1101/326470
Grassel, L., Fast, L. A., Scheffer, K. D., Boukhallouk, F., Spoden, G. A., Tenzer, S., et al. (2016). The CD63-syntenin-1 complex controls post-endocytic trafficking of oncogenic human papillomaviruses. Sci. Rep. 6:32337. doi: 10.1038/srep32337
Gruber, G., Hess, J., Stiefel, C., Aebersold, D. M., Zimmer, Y., Greiner, R. H., et al. (2005). Correlation between the tumoral expression of beta3-integrin and outcome in cervical cancer patients who had undergone radiotherapy. Br. J. Cancer 92, 41–46. doi: 10.1038/sj.bjc.6602278
Hazelbag, S., Kenter, G. G., Gorter, A., Dreef, E. J., Koopman, L. A., Violette, S. M., et al. (2007). Overexpression of the alpha v beta 6 integrin in cervical squamous cell carcinoma is a prognostic factor for decreased survival. J. Pathol. 212, 316–324. doi: 10.1002/path.2168
He, C., Mao, D., Hua, G., Lv, X., Chen, X., and Angeletti, P. C. (2015). The hippo/YAP pathway interacts with EGFR signaling and HPV oncoproteins to regulate cervical cancer progression. EMBO Mol. Med. 7, 1426–1449. doi: 10.15252/emmm.201404976
He, M., Cheng, Y., Li, W., Liu, Q., Liu, J., and Huang, J. (2010). Vascular endothelial growth factor C promotes cervical cancer metastasis via up-regulation and activation of RhoA/ROCK-2/moesin cascade. BMC Cancer 10:170. doi: 10.1186/1471-2407-10-170
Hellner, K., Mar, J., Fang, F., Quackenbush, J., and Munger, K. (2009). HPV16 E7 oncogene expression in normal human epithelial cells causes molecular changes indicative of an epithelial to mesenchymal transition. Virology 391, 57–63. doi: 10.1016/j.virol.2009.05.036
Herfs, M., Herman, L., Hubert, P., Minner, F., Arafa, M., Roncarati, P., et al. (2009). High expression of PGE2 enzymatic pathways in cervical (pre)neoplastic lesions and functional consequences for antigen-presenting cells. Cancer Immunol. Immunother. 58, 603–614. doi: 10.1007/s00262-008-0584-4
Herfs, M., Hubert, P., Kholod, N., Caberg, J. H., Gilles, C., Berx, G., et al. (2008). Transforming growth factor-beta1-mediated Slug and Snail transcription factor up-regulation reduces the density of langerhans cells in epithelial metaplasia by affecting E-cadherin expression. Am. J. Pathol. 172, 1391–1402. doi: 10.2353/ajpath.2008.071004
Herfs, M., Hubert, P., Moutschen, M., and Delvenne, P. (2011). Mucosal junctions: open doors to HPV and HIV infections? Trends Microbiol. 19, 114–120. doi: 10.1016/j.tim.2010.12.006
Herfs, M., Soong, T. R., Delvenne, P., and Crum, C. P. (2017). Deciphering the multifactorial susceptibility of mucosal junction cells to HPV infection and related carcinogenesis. Viruses 9:85. doi: 10.3390/v9040085
Herfs, M., Yamamoto, Y., Laury, A., Wang, X., Nucci, M. R., Mclaughlin-Drubin, M. E., et al. (2012). A discrete population of squamocolumnar junction cells implicated in the pathogenesis of cervical cancer. Proc. Natl. Acad. Sci. U.S.A. 109, 10516–10521. doi: 10.1073/pnas.1202684109
Hou, T., Zhang, W., Tong, C., Kazobinka, G., Huang, X., Huang, Y., et al. (2015). Putative stem cell markers in cervical squamous cell carcinoma are correlated with poor clinical outcome. BMC Cancer 15:785. doi: 10.1186/s12885-015-1826-4
Huang, R., and Rofstad, E. K. (2017). Cancer stem cells (CSCs), cervical CSCs and targeted therapies. Oncotarget 8, 35351–35367. doi: 10.18632/oncotarget.10169
Hubert, P., Herman, L., Roncarati, P., Maillard, C., Renoux, V., Demoulin, S., et al. (2014). Altered alpha-defensin 5 expression in cervical squamocolumnar junction: implication in the formation of a viral/tumour-permissive microenvironment. J. Pathol. 234, 464–477. doi: 10.1002/path.4435
Ikeda, K., Tate, G., Suzuki, T., and Mitsuya, T. (2008). Coordinate expression of cytokeratin 8 and cytokeratin 17 immunohistochemical staining in cervical intraepithelial neoplasia and cervical squamous cell carcinoma: an immunohistochemical analysis and review of the literature. Gynecol. Oncol. 108, 598–602. doi: 10.1016/j.ygyno.2007.11.042
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
Jiang, A., Bloom, O., Ono, S., Cui, W., Unternaehrer, J., Jiang, S., et al. (2007). Disruption of E-cadherin-mediated adhesion induces a functionally distinct pathway of dendritic cell maturation. Immunity 27, 610–624. doi: 10.1016/j.immuni.2007.08.015
Kanduc, D. (2002). Translational regulation of human papillomavirus type 16 E7 mRNA by the peptide SEQIKA, shared by rabbit alpha(1)-globin and human cytokeratin 7. J. Virol. 76, 7040–7048. doi: 10.1128/JVI.76.14.7040-7048.2002
Kodama, J., Hashimoto, I., Seki, N., Hongo, A., Yoshinouchi, M., Okuda, H., et al. (2001). Thrombospondin-1 and−2 messenger RNA expression in invasive cervical cancer: correlation with angiogenesis and prognosis. Clin. Cancer Res. 7, 2826–2831.
Li, F., and Cui, J. (2015). Human telomerase reverse transcriptase regulates vascular endothelial growth factor expression via human papillomavirus oncogene E7 in HPV-18-positive cervical cancer cells. Med. Oncol. 32:199. doi: 10.1007/s12032-015-0649-0
Liu, Q., Wang, J., Tang, M., Chen, L., Qi, X., Li, J., et al. (2018). The overexpression of PXN promotes tumor progression and leads to radioresistance in cervical cancer. Future Oncol. 14, 241–253. doi: 10.2217/fon-2017-0474
Lopez-Ocejo, O., Viloria-Petit, A., Bequet-Romero, M., Mukhopadhyay, D., Rak, J., and Kerbel, R. S. (2000). Oncogenes and tumor angiogenesis: the HPV-16 E6 oncoprotein activates the vascular endothelial growth factor (VEGF) gene promoter in a p53 independent manner. Oncogene 19, 4611–4620. doi: 10.1038/sj.onc.1203817
Lu, X., Jiang, L., Zhang, L., Zhu, Y., Hu, W., Wang, J., et al. (2019). Immune signature-based subtypes of cervical squamous cell carcinoma tightly associated with human papillomavirus type 16 expression, molecular features, and clinical outcome. Neoplasia 21, 591–601. doi: 10.1016/j.neo.2019.04.003
Maziveyi, M., and Alahari, S. K. (2017). Cell matrix adhesions in cancer: the proteins that form the glue. Oncotarget 8, 48471–48487. doi: 10.18632/oncotarget.17265
Mendoza-Almanza, G., Ortiz-Sanchez, E., Rocha-Zavaleta, L., Rivas-Santiago, C., Esparza-Ibarra, E., and Olmos, J. (2019). Cervical cancer stem cells and other leading factors associated with cervical cancer development. Oncol. Lett. 18, 3423–3432. doi: 10.3892/ol.2019.10718
Mills, A. M., Paquette, C., Terzic, T., Castle, P. E., and Stoler, M. H. (2017). CK7 Immunohistochemistry as a predictor of CIN1 progression: a retrospective study of patients from the quadrivalent HPV vaccine trials. Am. J. Surg. Pathol. 41, 143–152. doi: 10.1097/PAS.0000000000000747
Mittal, D., Kassianos, A. J., Tran, L. S., Bergot, A. S., Gosmann, C., Hofmann, J., et al. (2013). Indoleamine 2,3-dioxygenase activity contributes to local immune suppression in the skin expressing human papillomavirus oncoprotein e7. J. Invest. Dermatol. 133, 2686–2694. doi: 10.1038/jid.2013.222
Miyake, M., Hori, S., Morizawa, Y., Tatsumi, Y., Toritsuka, M., and Ohnishi, S. (2017). Collagen type IV alpha 1 (COL4A1) and collagen type XIII alpha 1 (COL13A1) produced in cancer cells promote tumor budding at the invasion front in human urothelial carcinoma of the bladder. Oncotarget 8, 36099–36114. doi: 10.18632/oncotarget.16432
Moody, C. A., and Laimins, L. A. (2010). Human papillomavirus oncoproteins: pathways to transformation. Nat. Rev. Cancer 10, 550–560. doi: 10.1038/nrc2886
Morbini, P., Capello, G. L., Alberizzi, P., Benazzo, M., Paglino, C., Comoli, P., et al. (2015). Markers of squamocolumnar junction cells in normal tonsils and oropharyngeal cancer with and without HPV infection. Histol. Histopathol. 30, 833–839. doi: 10.14670/HH-11-590
Muller, M., and Demeret, C. (2012). The HPV E2-host protein-protein interactions: a complex hijacking of the cellular network. Open Virol. J. 6, 173–189. doi: 10.2174/1874357901206010173
Oktay, M. H., Oktay, K., Hamele-Bena, D., Buyuk, A., and Koss, L. G. (2003). Focal adhesion kinase as a marker of malignant phenotype in breast and cervical carcinomas. Hum. Pathol. 34, 240–245. doi: 10.1053/hupa.2003.40
Paquette, C., Mills, A. M., and Stoler, M. H. (2016). Predictive value of cytokeratin 7 immunohistochemistry in cervical low-grade squamous intraepithelial lesion as a marker for risk of progression to a high-grade lesion. Am. J. Surg. Pathol. 40, 236–243. doi: 10.1097/PAS.0000000000000548
Peng, S., Qiu, J., Yang, A., Yang, B., Jeang, J., Wang, J. W., et al. (2016). Optimization of heterologous DNA-prime, protein boost regimens and site of vaccination to enhance therapeutic immunity against human papillomavirus-associated disease. Cell. Biosci. 6:16. doi: 10.1186/s13578-016-0080-z
Sato, A., Ishiwata, T., Matsuda, Y., Yamamoto, T., Asakura, H., Takeshita, T., et al. (2012). Expression and role of nestin in human cervical intraepithelial neoplasia and cervical cancer. Int. J. Oncol. 41, 441–448. doi: 10.3892/ijo.2012.1473
Sato, M., Kawana, K., Adachi, K., Fujimoto, A., Yoshida, M., Nakamura, H., et al. (2017). Regeneration of cervical reserve cell-like cells from human induced pluripotent stem cells (iPSCs): a new approach to finding targets for cervical cancer stem cell treatment. Oncotarget 8, 40935–40945. doi: 10.18632/oncotarget.16783
Schiffman, M., Doorbar, J., Wentzensen, N., De Sanjose, S., Fakhry, C., Monk, B. J., et al. (2016). Carcinogenic human papillomavirus infection. Nat. Rev. Dis. Prim. 2:16086. doi: 10.1038/nrdp.2016.86
Tomaic, V. (2016). Functional roles of E6 and E7 oncoproteins in HPV-induced malignancies at diverse anatomical sites. Cancers. 8:95. doi: 10.3390/cancers8100095
Tommasino, M. (2014). The human papillomavirus family and its role in carcinogenesis. Semin. Cancer Biol. 26, 13–21. doi: 10.1016/j.semcancer.2013.11.002
Trugilo, K. P., Cebinelli, G. C. M., Berti, F. C. B., Okuyama, N. C. M., Cezar-Dos-Santos, F., Sena, M. M., et al. (2019). Polymorphisms in the TGFB1 signal peptide influence human papillomavirus infection and development of cervical lesions. Med. Microbiol. Immunol. 208, 49–58. doi: 10.1007/s00430-018-0557-y
Tuong, Z. K., Noske, K., Kuo, P., Bashaw, A. A., Teoh, S. M., and Frazer, I. H. (2018). Murine HPV16 E7-expressing transgenic skin effectively emulates the cellular and molecular features of human high-grade squamous intraepithelial lesions. Papillomavirus Res. 5, 6–20. doi: 10.1016/j.pvr.2017.10.001
Vaira, V., Faversani, A., Dohi, T., Maggioni, M., Nosotti, M., Tosi, D., et al. (2011). Aberrant overexpression of the cell polarity module scribble in human cancer. Am. J. Pathol. 178, 2478–2483. doi: 10.1016/j.ajpath.2011.02.028
Wang, H. Y., Sun, J. M., Lu, H. F., Shi, D. R., Ou, Z. L., Ren, Y. L., et al. (2006). Micrometastases detected by cytokeratin 19 expression in sentinel lymph nodes of patients with early-stage cervical cancer. Int. J. Gynecol. Cancer. 16, 643–648. doi: 10.1111/j.1525-1438.2006.00381.x
Westra, W. H. (2012). The morphologic profile of HPV-related head and neck squamous carcinoma: implications for diagnosis, prognosis, and clinical management. Head Neck Pathol. 6(Suppl. 1), S48–54. doi: 10.1007/s12105-012-0371-6
Wiens, M. E., and Smith, J. G. (2017). α-defensin HD5 inhibits human papillomavirus 16 infection via capsid stabilization and redirection to the lysosome. mBio 8, e02304–16. doi: 10.1128/mBio.02304-16
Wu, M. P., Tzeng, C. C., Wu, L. W., Huang, K. F., and Chou, C. Y. (2004). Thrombospondin-1 acts as a fence to inhibit angiogenesis that occurs during cervical carcinogenesis. Cancer J. 10, 27–32. doi: 10.1097/00130404-200401000-00007
Yang, C., Xu, X., and Jin, H. (2016). Identification of potential miRNAs and candidate genes of cervical intraepithelial neoplasia by bioinformatic analysis. Eur. J. Gynaecol. Oncol. 37, 469–473. doi: 10.12892/ejgo3131.2016
Yang, E. J., Quick, M. C., Hanamornroongruang, S., Lai, K., Doyle, L. A., Mckeon, F. D., et al. (2015). Microanatomy of the cervical and anorectal squamocolumnar junctions: a proposed model for anatomical differences in HPV-related cancer risk. Mod. Pathol. 28, 994–1000. doi: 10.1038/modpathol.2015.54
Yang, S., Liu, Y., Xia, B., Deng, J., Liu, T., Li, Q., et al. (2016). DLL4 as a predictor of pelvic lymph node metastasis and a novel prognostic biomarker in patients with early-stage cervical cancer. Tumour. Biol. 37, 5063–5074. doi: 10.1007/s13277-015-4312-3
Yeo-Teh, N. S. L., Ito, Y., and Jha, S. (2018). High-risk human papillomaviral oncogenes E6 and E7 target key cellular pathways to achieve oncogenesis. Int. J. Mol. Sci. 19:1706. doi: 10.3390/ijms19061706
Yim, E. K., Meoyng, J., Namakoong, S. E., Um, S. J., and Park, J. S. (2004). Genomic and proteomic expression patterns in HPV-16 E6 gene transfected stable human carcinoma cell lines. DNA Cell. Biol. 23, 826–835. doi: 10.1089/dna.2004.23.826
Yim, E. K., and Park, J. S. (2005). The role of HPV E6 and E7 oncoproteins in HPV-associated cervical carcinogenesis. Cancer Res. Treat. 37, 319–324. doi: 10.4143/crt.2005.37.6.319
Yuan, Y., Min, S. J., Xu, D. Q., Shen, Y., Yan, H. Y., Wang, Y., et al. (2018). Expressions of VEGF and miR-21 in tumor tissues of cervical cancer patients with HPV infection and their relationships with prognosis. Eur. Rev. Med. Pharmacol. Sci. 22, 6274–6279. doi: 10.26355/eurrev_201810_1
Zapatka, M., Borozan, I., Brewer, D. S., Iskar, M., Grundhoff, A., Alawi, M., et al. (2020). The landscape of viral associations in human cancers. Nat. Genet. 52, 320–330. doi: 10.1038/s41588-019-0558-9
Zhou, C., Tuong, Z. K., and Frazer, I. H. (2019). Papillomavirus immune evasion strategies target the infected cell and the local immune system. Front. Oncol. 9:682. doi: 10.3389/fonc.2019.00682
Zhou, J., Yi, Q., and Tang, L. (2019). The roles of nuclear focal adhesion kinase (FAK) on cancer: a focused review. J. Exp. Clin. Cancer Res. 38:250. doi: 10.1186/s13046-019-1265-1
Zhou, Y., Shu, C., and Huang, Y. (2019). Fibronectin promotes cervical cancer tumorigenesis through activating FAK signaling pathway. J. Cell. Biochem. 120, 10988–10997. doi: 10.1002/jcb.28282
Keywords: HPV, cervical cancer, cancer stem cells, Notch, Wnt, Hedgehog, Hippo, Focal Adhesion
Citation: Budhwani M, Lukowski SW, Porceddu SV, Frazer IH and Chandra J (2020) Dysregulation of Stemness Pathways in HPV Mediated Cervical Malignant Transformation Identifies Potential Oncotherapy Targets. Front. Cell. Infect. Microbiol. 10:307. doi: 10.3389/fcimb.2020.00307
Received: 03 March 2020; Accepted: 20 May 2020;
Published: 25 June 2020.
Edited by:
Hasan Zaki, University of Texas Southwestern Medical Center, United StatesReviewed by:
Massimo Costalonga, University of Minnesota Twin Cities, United StatesThorsten Maretzky, The University of Iowa, United States
Copyright © 2020 Budhwani, Lukowski, Porceddu, Frazer and Chandra. 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: Janin Chandra, ai5jaGFuZHJhJiN4MDAwNDA7dXEuZWR1LmF1