Effect of Productive Human Papillomavirus 16 Infection on Global Gene Expression in Cervical Epithelium

Diversity in Proprotein Convertase Reactivity among Human Papillomavirus Types

The cleavage of viral surface proteins by furin is associated with some viruses' high virulence and infectivity. The human papillomavirus (HPV) requires the proteolytic processing of its capsid proteins for activation before entry. Variability in reactivity with furin and other proprotein convertases (PCs) among HPV types was investigated. HPV16, the most prevalent and carcinogenic HPV type, reacted with PCs with the broadest selectivity compared to other types in reactions of pseudoviral particles with the recombinant PCs, furin, PC4, PC5, PACE4, and PC7. Proteolytic preactivation was assessed using a well-established entry assay into PC-inhibited cells based on the green fluorescent protein as a reporter. The inhibition of the target cell PC activity with serpin-based PC-selective inhibitors also showed a diversity of PC selectivity among HPV types. HPV16 reacted with furin at the highest rate compared to the other types in time-dependent preactivation reactions and produced the highest entry values standardized to pseudoviral particle concentration. The predominant expression of furin in keratinocytes and the high reactivity of HPV16 with this enzyme highlight the importance of selectively targeting furin as a potential antiviral therapeutic approach.

Virus-Host Protein Interaction Network of the Hepatitis E Virus ORF2-4 by Mammalian Two-Hybrid Assays

Throughout their life cycle, viruses interact with cellular host factors, thereby influencing propagation, host range, cell tropism and pathogenesis. The hepatitis E virus (HEV) is an underestimated RNA virus in which knowledge of the virus-host interaction network to date is limited. Here, two related high-throughput mammalian two-hybrid approaches (MAPPIT and KISS) were used to screen for HEV-interacting host proteins. Promising hits were examined on protein function, involved pathway(s), and their relation to other viruses. We identified 37 ORF2 hits, 187 for ORF3 and 91 for ORF4. Several hits had functions in the life cycle of distinct viruses. We focused on SHARPIN and RNF5 as candidate hits for ORF3, as they are involved in the RLR-MAVS pathway and interferon (IFN) induction during viral infections. Knocking out (KO) SHARPIN and RNF5 resulted in a different IFN response upon ORF3 transfection, compared to wild-type cells. Moreover, infection was increased in SHARPIN KO cells and decreased in RNF5 KO cells. In conclusion, MAPPIT and KISS are valuable tools to study virus-host interactions, providing insights into the poorly understood HEV life cycle. We further provide evidence for two identified hits as new host factors in the HEV life cycle.

Multidimensional outlook on the pathophysiology of cervical cancer invasion and metastasis

Cervical cancer being one of the primary causes of high mortality rates among women is an area of concern, especially with ineffective treatment strategies. Extensive studies are carried out to understand various aspects of cervical cancer initiation, development and progression; however, invasive cervical squamous cell carcinoma has poor outcomes. Moreover, the advanced stages of cervical cancer may involve lymphatic circulation with a high risk of tumor recurrence at distant metastatic sites. Dysregulation of the cervical microbiome by human papillomavirus (HPV) together with immune response modulation and the occurrence of novel mutations that trigger genomic instability causes malignant transformation at the cervix. In this review, we focus on the major risk factors as well as the functionally altered signaling pathways promoting the transformation of cervical intraepithelial neoplasia into invasive squamous cell carcinoma. We further elucidate genetic and epigenetic variations to highlight the complexity of causal factors of cervical cancer as well as the metastatic potential due to the changes in immune response, epigenetic regulation, DNA repair capacity, and cell cycle progression. Our bioinformatics analysis on metastatic and non-metastatic cervical cancer datasets identified various significantly and differentially expressed genes as well as the downregulation of potential tumor suppressor microRNA miR-28-5p. Thus, a comprehensive understanding of the genomic landscape in invasive and metastatic cervical cancer will help in stratifying the patient groups and designing potential therapeutic strategies.

Contributions of viral oncogenes of HPV-18 and hypoxia to oxidative stress and genetic damage in human keratinocytes

Infection with high-risk human papillomaviruses like HPV-16 and HPV-18 is highly associated with the development of cervical and other cancers. Malignant transformation requires viral oncoproteins E5, E6 and E7, which promote cell proliferation and increase DNA damage. Oxidative stress and hypoxia are also key factors in cervical malignant transformation. Increased levels of reactive species of oxygen (ROS) and nitrogen (RNS) are found in the hypoxic tumor microenvironment, promoting genetic instability and invasiveness. In this work, we studied the combined effect of E5, E6 and E7 and hypoxia in increasing oxidative stress and promoting DNA damage and nuclear architecture alterations. HaCaT cells containing HPV-18 viral oncogenes (HaCaT E5/E6/E7-18) showed higher ROS levels in normoxia and higher levels of RNS in hypoxia compared to HaCaT parental cells, as well as higher genetic damage in hypoxia as measured by γH2AX and comet assays. In hypoxia, HaCaT E5/E6/E7-18 increased its nuclear dry mass and both cell types displayed marked heterogeneity in nuclear dry mass distribution and increased nuclear foci. Our results show contributions of both viral oncogenes and hypoxia to oxidative stress, DNA damage and altered nuclear architecture, exemplifying how an altered microenvironment combines with oncogenic transformation to promote tumor progression.

Human papillomavirus associated cervical lesion: pathogenesis and therapeutic interventions

Human papillomavirus (HPV) is the most prevalent sexually transmitted virus globally. Persistent high-risk HPV infection can result in cervical precancerous lesions and cervical cancer, with 70% of cervical cancer cases associated with high-risk types HPV16 and 18. HPV infection imposes a significant financial and psychological burden. Therefore, studying methods to eradicate HPV infection and halt the progression of precancerous lesions remains crucial. This review comprehensively explores the mechanisms underlying HPV-related cervical lesions, including the viral life cycle, immune factors, epithelial cell malignant transformation, and host and environmental contributing factors. Additionally, we provide a comprehensive overview of treatment methods for HPV-related cervical precancerous lesions and cervical cancer. Our focus is on immunotherapy, encompassing HPV therapeutic vaccines, immune checkpoint inhibitors, and advanced adoptive T cell therapy. Furthermore, we summarize the commonly employed drugs and other nonsurgical treatments currently utilized in clinical practice for managing HPV infection and associated cervical lesions. Gene editing technology is currently undergoing clinical research and, although not yet employed officially in clinical treatment of cervical lesions, numerous preclinical studies have substantiated its efficacy. Therefore, it holds promise as a precise treatment strategy for HPV-related cervical lesions.

Prognostic Value and Immune Infiltration of HPV-Related Genes in the Immune Microenvironment of Cervical Squamous Cell Carcinoma and Endocervical Adenocarcinoma

Mounting evidence has highlighted the immune environment as a critical feature in the development of cervical squamous cell carcinoma and endocervical adenocarcinoma (CESC). However, the relationship between the clinical characteristics of the immune environment and CESC remain unclear. Therefore, the aim of this study was to further characterize the relationship between the tumor and immune microenvironment and the clinical features of CESC using a variety of bioinformatic methods. Expression profiles (303 CESCs and three control samples) and relevant clinical data were obtained from The Cancer Genome Atlas. We divided CESC cases into different subtypes and performed a differential gene expression analysis. In addition, gene ontology (GO) and gene set enrichment analysis (GSEA) were performed to identify potential molecular mechanisms. Furthermore, data from 115 CESC patients from East Hospital were used to help identify the relationship between the protein expressions of key genes and disease-free survival using tissue microarray technology. Cases of CESC (n = 303) were divided into five subtypes (C1-C5) based on their expression profiles. A total of 69 cross-validated differentially expressed immune-related genes were identified. Subtype C4 demonstrated a downregulation of the immune profile, lower tumor immune/stroma scores, and worse prognosis. In contrast, the C1 subtype showed an upregulation of the immune profile, higher tumor immune/stroma scores, and better prognosis. A GO analysis suggested that changes in CESC were primarily enriched nuclear division, chromatin binding, and condensed chromosomes. In addition, GSEA demonstrated that cellular senescence, the p53 signaling pathway, and viral carcinogenesis are critical features of CESC. Moreover, high FOXO3 and low IGF-1 protein expression were closely correlated with decreased clinical prognosis. In summary, our findings provide novel insight into the relationship between the immune microenvironment and CESC. As such, our results may provide guidance for developing potential immunotherapeutic targets and biomarkers for CESC.

Differences in integration frequencies and APOBEC3 profiles of five high-risk HPV types adheres to phylogeny

Persistent infection with Human Papillomavirus (HPV) is responsible for almost all cases of cervical cancers, and HPV16 and HPV18 associated with the majority of these. These types differ in the proportion of viral minor nucleotide variants (MNVs) caused by APOBEC3 mutagenesis as well as integration frequencies. Whether these traits extend to other types remains uncertain. This study aimed to investigate and compare genomic variability and chromosomal integration in the two phylogenetically distinct Alpha-7 and Alpha-9 clades of carcinogenic HPV types. The TaME-seq protocol was employed to sequence cervical cell samples positive for HPV31, HPV33 or HPV45 and combine these with data from a previous study on HPV16 and HPV18. APOBEC3 mutation signatures were found in Alpha-9 (HPV16/31/33) but not in Alpha-7 (HPV18/45). HPV45 had significantly more MNVs compared to the other types. Alpha-7 had higher integration frequency compared to Alpha-9. An increase in integration frequency with increased diagnostic severity was found for Alpha-7. The results highlight important differences and broaden our understanding of the molecular mechanisms behind cervical cancer induced by high-risk HPV types from the Alpha-7 and Alpha-9 clades.

Like Brothers in Arms: How Hormonal Stimuli and Changes in the Metabolism Signaling Cooperate, Leading HPV Infection to Drive the Onset of Cervical Cancer

Despite all precautionary actions and the possibility of using vaccinations to counteract infections caused by human papillomaviruses (HPVs), HPV-related cancers still account for approximately 5% of all carcinomas. Worldwide, many women are still excluded from adequate health care due to their social position and origin. Therefore, immense efforts in research and therapy are still required to counteract the challenges that this disease entails. The special thing about an HPV infection is that it is not only able to trick the immune system in a sophisticated way, but also, through genetic integration into the host genome, to use all the resources available to the host cells to complete the replication cycle of the virus without activating the alarm mechanisms of immune recognition and elimination. The mechanisms utilized by the virus are the metabolic, immune, and hormonal signaling pathways that it manipulates. Since the virus is dependent on replication enzymes of the host cells, it also intervenes in the cell cycle of the differentiating keratinocytes and shifts their terminal differentiation to the uppermost layers of the squamocolumnar transformation zone (TZ) of the cervix. The individual signaling pathways are closely related and equally important not only for the successful replication of the virus but also for the onset of cervical cancer. We will therefore analyze the effects of HPV infection on metabolic signaling, as well as changes in hormonal and immune signaling in the tumor and its microenvironment to understand how each level of signaling interacts to promote tumorigenesis of cervical cancer.

Subtyping of Human Papillomavirus-Positive Cervical Cancers Based on the Expression Profiles of 50 Genes

Background

Human papillomavirus-positive (HPV+) cervical cancers are highly heterogeneous in molecular and clinical features. However, the molecular classification of HPV+ cervical cancers remains insufficiently unexplored.

Methods

Based on the expression profiles of 50 genes having the largest expression variations across the HPV+ cervical cancers in the TCGA-CESC dataset, we hierarchically clustered HPV+ cervical cancers to identify new subtypes. We further characterized molecular, phenotypic, and clinical features of these subtypes.

Results

We identified two subtypes of HPV+ cervical cancers, namely HPV+G1 and HPV+G2. We demonstrated that this classification method was reproducible in two validation sets. Compared to HPV+G2, HPV+G1 displayed significantly higher immune infiltration level and stromal content, lower tumor purity, lower stemness scores and intratumor heterogeneity (ITH) scores, higher level of genomic instability, lower DNA methylation level, as well as better disease-free survival prognosis. The multivariate survival analysis suggests that the disease-free survival difference between both subtypes is independent of confounding variables, such as immune signature, stemness, and ITH. Pathway and gene ontology analysis confirmed the more active tumor immune microenvironment in HPV+G1 versus HPV+G2.

Conclusions

HPV+ cervical cancers can be classified into two subtypes based on the expression profiles of the 50 genes with the largest expression variations across the HPV+ cervical cancers. Both subtypes have significantly different molecular, phenotypic, and clinical features. This new subtyping method captures the comprehensive heterogeneity in molecular and clinical characteristics of HPV+ cervical cancers and provides potential clinical implications for the diagnosis and treatment of this disease.

Metabolic reprogramming in cervical cancer and metabolomics perspectives

Cumulative studies have shown that metabolic reprogramming is a hallmark of malignant tumors. The emergence of technological advances, such as omics studies, has strongly contributed to the knowledge of cancer metabolism. Cervical cancer is among the most common cancers in women worldwide. Because cervical cancer is a virus-associated cancer and can exist in a precancerous state for years, investigations targeting the metabolic phenotypes of cervical cancer will enhance our understanding of the interference of viruses on host cells and the progression of cervical carcinogenesis. The purpose of this review was to illustrate metabolic perturbations in cervical cancer, the role that human papillomavirus (HPV) plays in remodeling cervical cell metabolism and recent approaches toward application of metabolomics in cervical disease research. Cervical cancer displays typical cancer metabolic profiles, including glycolytic switching, high lactate levels, lipid accumulation and abnormal kynurenine/tryptophan levels. HPV, at least in part, contributes to these alterations. Furthermore, emerging metabolomics data provide global information on the metabolic traits of cervical diseases and may aid in the discovery of biomarkers for diagnosis and therapy.

Comprehensive analysis of ceRNA networks in HPV16- and HPV18-mediated cervical cancers reveals XIST as a pivotal competing endogenous RNA

Cervical cancer (CC) is one of the most common cancers in women worldwide, being closely related to high-risk human papillomavirus (HR-HPVs). After a particular HR-HPV infects a cervical cell, transcriptional changes in the host cell are expected, including the regulation of lncRNAs, miRNAs, and mRNAs. Such transcripts may work independently or integrated in complex molecular networks - as in competing endogenous RNA (ceRNA) networks. In our research, we gathered transcriptome data from samples of HPV16/HPV18 cervical squamous cell carcinoma and endocervical adenocarcinoma (CESC), from The Cancer Genome Atlas (TCGA) project. Using GDCRNATools, we identified ceRNA networks that differentiate HPV16- from HPV18-mediated CESC. For HPV16-CESC, three lncRNA-mRNA co-expressed pairs were reported, all led by the X-inactive specific transcript (XIST): XIST | DLG5, XIST | LGR4, and XIST | ZNF81. The XIST | LGR4 and XIST | ZNF81 pairs shared 11 miRNAs, suggesting an increased impact on their final biological effect. XIST also stood out as an important lncRNA in HPV18-CESC, leading 35 of the 42 co-expressed pairs. Some mRNAs, such as ADAM9 and SLC38A2, emerged as important players in the ceRNA regulatory networks due to sharing a considerable amount of miRNAs with XIST. Furthermore, some XIST-associated axes, namely XIST | miR-23a-3p | LGR4 and XIST | miR-30b-5p or miR-30c-5p or miR-30e-5p I ADAM9, had a significant impact on the overall survival of HPV16- and HPV18-CESC patients, respectively. Together, these data suggest that XIST has an important role in HPV-mediated tumorigenesis, which may implicate different molecular signatures between HPV16 and HPV18-associated tumors.

Expression of Viral microRNAs in Serum and White Blood Cells of Cows Exposed to Bovine Leukemia Virus

Bovine leukemia virus (BLV) affects the health and productivity of cattle. The virus causes abnormal immune function and immunosuppression. MicroRNAs (miRNAs) are involved in gene expression, having been associated with stress and immune response, tumor growth, and viral infection. The objective of this study was to determine the expression of circulating miRNAs produced by BLV in animals exposed to the virus. Sera from 14 animals were collected to establish IgG reactivity to BLV by ELISA, where seven animals were seropositive and seven were seronegative for BLV exposure. White blood cells (WBC) from each animal were also collected and miRNAs were identified by sequencing from sera and WBC. The seropositive group had higher counts of BLV miRNAs when compared to seronegative group in sera and WBC. Blv-miR-1-3p, blv-miR-B2-5p, blv-miR-B4-3p, and blv-miR-B5-5p were statistically significant (P < 0.00001) in serum with an average of 7 log2 fold difference between seropositive and seronegative groups. Blv-miR-B1-3p, blv-miR-B1-5p, blv-miR-B3, blv-miR-B4-3p, blv-miR-B4-5p, blv-miR-B5-5p were statistically significant (P < 1.08e⁻⁹) in WBC with an average of 7 log2 fold difference between the seropositive and the seronegative groups. Blv-miR-B2-3p and blv-miR-B2-5p were also statistically significant in WBC (P < 2.79e-17), with an average of 27 log2 fold difference between the seropositive and the seronegative groups. There were 18 genes identified as being potential targets for blv-miR-B1-5p, and 3 genes for blv-miR-B4-5p. Gene ontology analysis indicated that the target genes are mainly involved in the response to stress and in the immune system process. Several of the identified genes have been associated with leukemia development in humans and cattle. Differential expression of genes targeted by BLV miRNAs should be evaluated to determine their effect in BLV replication.

Markers of Angiogenesis, Lymphangiogenesis, and Epithelial-Mesenchymal Transition (Plasticity) in CIN and Early Invasive Carcinoma of the Cervix: Exploring Putative Molecular Mechanisms Involved in Early Tumor Invasion

The establishment of a proangiogenic phenotype and epithelial-to-mesenchymal transition (EMT) are considered as critical events that promote the induction of invasive growth in epithelial tumors, and stimulation of lymphangiogenesis is believed to confer the capacity for early dissemination to cancer cells. Recent research has revealed substantial interdependence between these processes at the molecular level as they rely on common signaling networks. Of great interest are the molecular mechanisms of (lymph-)angiogenesis and EMT associated with the earliest stages of transition from intraepithelial development to invasive growth, as they could provide the source of potentially valuable tools for targeting tumor metastasis. However, in the case of early-stage cervical cancer, the players of (lymph-)angiogenesis and EMT processes still remain substantially uncharacterized. In this study, we used RNA sequencing to compare transcriptomes of HPV(+) preinvasive neoplastic lesions and early-stage invasive carcinoma of the cervix and to identify (lymph-)angiogenesis- and EMT-related genes and pathways that may underlie early acquisition of invasive phenotype and metastatic properties by cervical cancer cells. Second, we applied flow cytometric analysis to evaluate the expression of three key lymphangiogenesis/EMT markers (VEGFR3, MET, and SLUG) in epithelial cells derived from enzymatically treated tissue specimens. Overall, among 201 differentially expressed genes, a considerable number of (lymph-)angiogenesis and EMT regulatory factors were identified, including genes encoding cytokines, growth factor receptors, transcription factors, and adhesion molecules. Pathway analysis confirmed enrichment for angiogenesis, epithelial differentiation, and cell guidance pathways at transition from intraepithelial neoplasia to invasive carcinoma and suggested immune-regulatory/inflammatory pathways to be implicated in initiation of invasive growth of cervical cancer. Flow cytometry showed cell phenotype-specific expression pattern for VEGFR3, MET, and SLUG and revealed correlation with the amount of tumor-infiltrating lymphocytes at the early stages of cervical cancer progression. Taken together, these results extend our understanding of driving forces of angiogenesis and metastasis in HPV-associated cervical cancer and may be useful for developing new treatments.

The Relationship between Estrogen-Related Signaling and Human Papillomavirus Positive Cancers

High risk-human papillomaviruses (HPVs) are known carcinogens. Numerous reports have linked the steroid hormone estrogen, and the expression of estrogen receptors (ERs), to HPV-related cancers, although the exact nature of the interactions remains to be fully elucidated. Here we will focus on estrogen signaling and describe both pro and potentially anti-cancer effects of this hormone in HPV-positive cancers. This review will summarize: (1) cell culture-related evidence, (2) animal model evidence, and (3) clinical evidence demonstrating an interaction between estrogen and HPV-positive cancers. This comprehensive review provides insights into the potential relationship between estrogen and HPV. We suggest that estrogen may provide a potential therapeutic for HPV-related cancers, however additional studies are necessary.

Modelling human papillomavirus biology in oropharyngeal keratinocytes

Most human papillomavirus (HPV) positive head and neck cancers arise in the tonsil crypts; deep invaginations at the tonsil surface that are lined with reticulated epithelium infiltrated by immune cells. As in cervical HPV infections, HPV16 is the most prevalent high-risk type in the oropharyngeal cancers, and a genital-oral route of infection is most likely. However, the natural history of HPV-driven oropharyngeal pathogenesis is an enigma, although there is evidence that it is different to that of cervical disease. It is not known if the virus establishes a productive or abortive infection in keratinocytes of the tonsil crypt, or if viral infections progress to cancer via a neoplastic phase, as in cervical HPV infection. The HPV DNA is more frequently found unintegrated in the cancers of the oropharynx compared to those that arise in the cervix, and may include novel HPV-human DNA hybrids episomes. Here, we review current understanding of HPV biology in the oropharynx and discuss the cell-based systems being used to model the HPV life cycle in tonsil keratinocytes and how they can be used to inform on HPV-driven neoplastic progression in the oropharynx.

This article is part of the theme issue ‘Silent cancer agents: multi-disciplinary modelling of human DNA oncoviruses’.

Estrogen Attenuates the Growth of Human Papillomavirus-Positive Epithelial Cells

Human papillomaviruses cause around 5% of all human cancers, yet there are no specific antiviral therapeutic approaches available for combatting these cancers. These cancers are currently treated with standard chemoradiation therapy (CRT). Specific antiviral reagents are desperately required, particularly for HPV+HNSCC whose incidence is increasing and for which there are no diagnostic tools available for combatting this disease. Using data from The Cancer Genome Atlas (TCGA), we and others determined that the estrogen receptor alpha (ERα) is overexpressed in HPV+HNSCC and that elevated levels are associated with an improved disease outcome. This has led to the proposal that estrogen treatment could be a novel therapeutic approach for combatting HPV+cancers. Here, we demonstrate that estrogen attenuates the growth of HPV+epithelial cells using multiple mechanisms, supporting the idea that estrogen has potential as a therapeutic agent for the treatment of HPV+HNSCC.

Human papillomaviruses (HPVs) are small, double-stranded DNA viruses that are significant risk factors in the development of cancer, and HPV accounts for approximately 5% of all worldwide cancers. Recent studies using data from The Cancer Genome Atlas (TCGA) have demonstrated that elevated levels of estrogen receptor alpha (ERα) are associated with improved survival in oropharyngeal cancers, and these elevated receptor levels were linked with human papillomavirus-positive cancers (HPV+cancers). There has been a dramatic increase in HPV-related head and neck squamous cell carcinomas (HPV+HNSCCs) over the last 2 decades, and therapeutic options for this ongoing health crisis are a priority; currently, there are no antiviral therapeutics available for combatting HPV+cancers. During our TGCA studies on head and neck cancer, we had also discovered the overexpression of ERα in HPV+cancers. Here, we demonstrate that 17β-estradiol (estrogen) attenuates the growth/cell viability of HPV+cancers in vitro, but not HPV-negative cancer cells. In addition, N/Tert-1 cells (foreskin keratinocytes immortalized with human telomerase reverse transcriptase [hTERT]) containing human papillomavirus 16 (HPV16) have elevated levels of ERα and growth sensitivity after estrogen treatment compared with parental N/Tert-1 cells. Finally, we demonstrate that there are potentially two mechanisms contributing to the attenuation of HPV+ cell growth following estrogen treatment. First, estrogen represses the viral transcriptional long control region (LCR) downregulating early gene expression, including E6/E7. Second, expression of E6 and E7 by themselves sensitizes cells to estrogen. Overall, our results support the recent proposal that estrogen could be exploited therapeutically for the treatment of HPV-positive oral cancers.

IMPORTANCE Human papillomaviruses cause around 5% of all human cancers, yet there are no specific antiviral therapeutic approaches available for combatting these cancers. These cancers are currently treated with standard chemoradiation therapy (CRT). Specific antiviral reagents are desperately required, particularly for HPV+HNSCC whose incidence is increasing and for which there are no diagnostic tools available for combatting this disease. Using data from The Cancer Genome Atlas (TCGA), we and others determined that the estrogen receptor alpha (ERα) is overexpressed in HPV+HNSCC and that elevated levels are associated with an improved disease outcome. This has led to the proposal that estrogen treatment could be a novel therapeutic approach for combatting HPV+cancers. Here, we demonstrate that estrogen attenuates the growth of HPV+epithelial cells using multiple mechanisms, supporting the idea that estrogen has potential as a therapeutic agent for the treatment of HPV+HNSCC.

Comparative transcriptome analysis of the human endocervix and ectocervix during the proliferative and secretory phases of the menstrual cycle

Despite extensive studies suggesting increased susceptibility to HIV during the secretory phase of the menstrual cycle, the molecular mechanisms involved remain unclear. Our goal was to analyze transcriptomes of the endocervix and ectocervix during the proliferative and secretory phases using RNA sequencing to explore potential molecular signatures of susceptibility to HIV. We identified 202 differentially expressed genes (DEGs) between the proliferative and secretory phases of the cycle in the endocervix (adjusted p < 0.05). The biofunctions and pathways analysis of DEGs revealed that cellular assembly and epithelial barrier function in the proliferative phase and inflammatory response/cellular movement in the secretory phase were among the top biofunctions and pathways. The gene set enrichment analysis of ranked DEGs (score = log fold change/p value) in the endocervix and ectocervix revealed that (i) unstimulated/not activated immune cells gene sets positively correlated with the proliferative phase and negatively correlated with the secretory phase in both tissues, (ii) IFNγ and IFNα response gene sets positively correlated with the proliferative phase in the ectocervix, (iii) HIV restrictive Wnt/β-catenin signaling pathway negatively correlated with the secretory phase in the endocervix. Our data show menstrual cycle phase-associated changes in both endocervix and ectocervix, which may modulate susceptibility to HIV.

Suppression of Stromal Interferon Signaling by Human Papillomavirus 16

Human papillomaviruses (HPVs) infect squamous epithelia and cause several important cancers. Immune evasion is critical for viral persistence. Fibroblasts in the stromal microenvironment provide growth signals and cytokines that are required for proper epithelial differentiation, maintenance, and immune responses and are critical in the development of many cancers. In this study, we examined the role of epithelial-stromal interactions in the HPV16 life cycle using organotypic (raft) cultures as a model. Rafts were created using uninfected human foreskin keratinocytes (HFKs) and HFKs containing either wild-type HPV16 or HPV16 with a stop mutation to prevent the expression of the viral oncogene E5. Microarray analysis revealed significant changes in gene expression patterns in the stroma in response to HPV16, some of which were E5 dependent. Interferon (IFN)-stimulated genes (ISGs) and extracellular matrix remodeling genes were suppressed, the most prominent pathways affected. STAT1, IFNAR1, IRF3, and IRF7 were knocked down in stromal fibroblasts using lentiviral short hairpin RNA (shRNA) transduction. HPV late gene expression and viral copy number in the epithelium were increased when the stromal IFN pathway was disrupted, indicating that the stroma helps control the late phase of the HPV life cycle in the epithelium. Increased late gene expression correlated with increased late keratinocyte differentiation but not decreased IFN signaling in the epithelium. These studies show HPV16 has a paracrine effect on stromal innate immunity, reveal a new role for E5 as a stromal innate immune suppressor, and suggest that stromal IFN signaling may influence keratinocyte differentiation.IMPORTANCE The persistence of high-risk human papillomavirus (HPV) infections is the key risk factor for developing HPV-associated cancers. The ability of HPV to evade host immunity is a critical component of its ability to persist. The environment surrounding a tumor is increasingly understood to be critical in cancer development, including immune evasion. Our studies show that HPV can suppress the expression of immune-related genes in neighboring fibroblasts in a three-dimensional (3D) model of human epithelium. This finding is significant, because it indicates that HPV can control innate immunity not only in the infected cell but also in the microenvironment. In addition, the ability of HPV to regulate stromal gene expression depends in part on the viral oncogene E5, revealing a new function for this protein as an immune evasion factor.

Transregulation of microRNA miR-21 promoter by AP-1 transcription factor in cervical cancer cells

Background

Gene expression profiles have demonstrated that miR-21 expression is altered in almost all types of cancers and it has been classified as an oncogenic microRNA. Persistent HPV infection is the main etiologic agent in cervical cancer and induces genetic instability, including disruption of microRNA gene expression. In the present study, we analyzed the underlying mechanism of how AP-1 transcription factor can active miR-21 gene expression in cervical cancer cells.

Methods

To identify that c-Fos and c-Jun regulate the expression of miR-21 we performed RT-qPCR and western blot assays. We analyzed the interaction of AP-1 with miR-21 promoter by EMSA and ChIP assays and determined the mechanism of its regulation by reporter construct plasmids. We identified the nuclear translocation of c-Fos and c-Jun by immunofluorescence microscopy assays.

Results

We demonstrated that c-Fos and c-Jun proteins are expressed and regulate the expression of miR-21 in cervical cancer cells. DNA sequence analysis revealed the presence of AP-1 DNA-binding sites in the human miR-21 promoter region. EMSA analyses confirmed the interactions of the miR-21 upstream transcription factor AP-1. ChIP assays further showed the binding of c-Fos to AP-1 sequences from the miR-21 core promoter in vivo. Functional analysis of AP-1 sequences of miR-21 in reporter plasmids demonstrated that these sequences increase the miR-21 promoter activation.

Conclusions

Our findings suggest a physical interaction and functional cooperation between AP-1 transcription factor in the miR-21 promoter and may explain the effect of AP-1 on miR-21 gene expression in cervical cancer cells.

Single-cell RNA sequencing reveals cell type-specific HPV expression in hyperplastic skin lesions

Human Papillomavirus infection is highly prevalent worldwide. While most types of HPV cause benign warts, some high-risk types are known to cause cervical cancer, as well as cancer of the oral cavity and head and neck. Persistent cutaneous HPV infection can be particularly problematic in patients with chronic immunosuppression, for example following organ transplantation. Due to unknown mechanisms, these patients may develop numerous warts, as well as present with a dramatically increased skin cancer prevalence. Despite an association between HPV persistence in the epidermis and excessive wart or squamous cancer development, the molecular mechanisms linking immunosuppression, HPV expression and excessive epidermal proliferation have not been determined, largely due to low-sensitivity methodology to capture rare viral transcription events. Here, we use single-cell RNA sequencing to profile HPV-positive skin lesions from an immunosuppressed patient that were found to express the alphapapillomavirus HPV78 in basal keratinocytes, suprabasal keratinocytes and hair follicle stem cells. This method can be applied to detect and investigate HPV transcripts in cutaneous lesions, allowing mechanistic links between immunosuppression-induced HPV life cycle and epidermal hyperproliferation to be uncovered.

Tissue-Specific Gene Expression during Productive Human Papillomavirus 16 Infection of Cervical, Foreskin, and Tonsil Epithelium

Epidemiological data confirm a much higher incidence of high-risk human papillomavirus 16 (HPV16)-mediated carcinogenesis of the cervical epithelium than for other target sites. In order to elucidate tissue-specific responses to virus infection, we compared gene expression changes induced by productive HPV16 infection of cervical, foreskin, and tonsil organotypic rafts. These rafts closely mimic persistent HPV16 infection, long before carcinogenesis sets in. The total number of gene expression changes varied considerably across the tissue types, with only 32 genes being regulated in common. Among them, we confirmed the Kelch-like family protein KLHL35 and the laminin-5 complex to be upregulated and downregulated, respectively, in all the three tissues. HPV16 infection induces upregulation of genes involved in cell cycle control, cell division, mitosis, DNA replication, and DNA damage repair in all the three tissues, indicative of a hyperproliferative environment. In the cervical and tonsil epithelium, we observe significant downregulation of genes involved in epidermis development, keratinocyte differentiation, and extracellular matrix organization. On the other hand, in HPV16-positive foreskin (HPV16 foreskin) tissue, several genes involved in interferon-mediated innate immunity, cytokine signaling, and cellular defenses were downregulated. Furthermore, pathway analysis and experimental validations identified important cellular pathways like STAT1 and transforming growth factor β (TGF-β) to be differentially regulated among the three tissue types. The differential modulation of important cellular pathways like TGF-β1 and STAT1 can explain the sensitivity of tissues to HPV cancer progression.IMPORTANCE Although the high-risk human papillomavirus 16 infects anogenital and oropharyngeal sites, the cervical epithelium has a unique vulnerability to progression of cancer. Host responses during persistent infection and preneoplastic stages can shape the outcome of cancer progression in a tissue-dependent manner. Our study for the first time reports differential regulation of critical cellular functions and signaling pathways during productive HPV16 infection of cervical, foreskin, and tonsil tissues. While the virus induces hyperproliferation in infected cells, it downregulates epithelial differentiation, epidermal development, and innate immune responses, according to the tissue type. Modulation of these biological functions can determine virus fitness and pathogenesis and illuminate key cellular mechanisms that the virus employs to establish persistence and finally initiate disease progression.

HPV16 whole genome minority variants in persistent infections from young Dutch women

BACKGROUND

Chronic infections by one of the oncogenic human papillomaviruses (HPVs) are responsible for near 5% of the global cancer burden and HPV16 is the type most often found in cancers. HPV genomes display unexpected levels of variation when deep-sequenced. Minor nucleotide variations (MNVs) may reveal HPV genomic instability and HPV-related carcinogenic transformation of host cells.

OBJECTIVES

The objective of this study was to investigate HPV16 genome variation at the minor variant level on persisting HPV16 cervical infections from a population of young Dutch women.

STUDY DESIGN

15 HPV16 infections were sequenced using a whole-HPV genome deep sequencing protocol (TaME-seq). One infection was followed over a three-year period, eight were followed over a two-year period, three were followed over a one-year period and three infections had a single sampling point.

RESULTS AND CONCLUSIONS

Using a 1% variant frequency cutoff, we find on average 48 MNVs per HPV16 genome and 1717 MNVs in total when sequencing coverage was >100 × . We find the transition mutation T > C to be the most common, in contrast to other studies detecting APOBEC-related C > T mutation profiles in pre-cancerous and cancer samples. Our results suggest that the relative mutagenic footprint of HPV16 genomes may differ between the infections in this study and transforming lesions. In addition, we identify a number of MNVs that have previously been associated with higher incidence of high-grade lesions (CIN3+) in a population study. These findings may provide a starting point for future studies exploring causality between emerging HPV minor genomic variants and cancer development.

Combined Transcriptome and Proteome Analysis of Immortalized Human Keratinocytes Expressing Human Papillomavirus 16 (HPV16) Oncogenes Reveals Novel Key Factors and Networks in HPV-Induced Carcinogenesis

Human papillomavirus (HPV)-associated cancers still remain a big health problem, especially in developing countries, despite the availability of prophylactic vaccines. Although HPV oncogenes have been intensively investigated for decades, a study applying recent advances in RNA-Seq and quantitative proteomic approaches to a precancerous model system with well-defined HPV oncogene expression alongside HPV-negative parental cells has been missing until now. Here, combined omics analyses reveal global changes caused by the viral oncogenes in a less biased way and allow the identification of novel factors and key cellular networks potentially promoting malignant transformation. In addition, this system also provides a basis for mechanistic research on novel key factors regulated by HPV oncogenes, especially those that are confirmed in vivo in cervical cancer as well as in head and neck cancer patient samples from TCGA data sets.

Although the role of high-risk human papillomaviruses (hrHPVs) as etiological agents in cancer development has been intensively studied during the last decades, there is still the necessity of understanding the impact of the HPV E6 and E7 oncogenes on host cells, ultimately leading to malignant transformation. Here, we used newly established immortalized human keratinocytes with a well-defined HPV16 E6E7 expression cassette to get a more complete and less biased overview of global changes induced by HPV16 by employing transcriptome sequencing (RNA-Seq) and stable isotope labeling by amino acids in cell culture (SILAC). This is the first study combining transcriptome and proteome data to characterize the impact of HPV oncogenes in human keratinocytes in comparison with their virus-negative counterparts. To enhance the informative value and accuracy of the RNA-Seq data, four different bioinformatic workflows were used. We identified potential novel upstream regulators (e.g., CNOT7, SPDEF, MITF, and PAX5) controlling distinct clusters of genes within the HPV-host cell network as well as distinct factors (e.g., CPPED1, LCP1, and TAGLN) with essential functions in cancer. Validated results in this study were compared to data sets from The Cancer Genome Atlas (TCGA), demonstrating that several identified factors were also differentially expressed in cervical squamous cell carcinoma and endocervical adenocarcinoma (CESC) and HPV-positive head and neck squamous cell carcinomas (HNSCs). This highly integrative approach allows the identification of novel HPV-induced cellular changes that are also reflected in cancer patients, providing a promising omics data set for future studies in both basic and translational research.

IMPORTANCE Human papillomavirus (HPV)-associated cancers still remain a big health problem, especially in developing countries, despite the availability of prophylactic vaccines. Although HPV oncogenes have been intensively investigated for decades, a study applying recent advances in RNA-Seq and quantitative proteomic approaches to a precancerous model system with well-defined HPV oncogene expression alongside HPV-negative parental cells has been missing until now. Here, combined omics analyses reveal global changes caused by the viral oncogenes in a less biased way and allow the identification of novel factors and key cellular networks potentially promoting malignant transformation. In addition, this system also provides a basis for mechanistic research on novel key factors regulated by HPV oncogenes, especially those that are confirmed in vivo in cervical cancer as well as in head and neck cancer patient samples from TCGA data sets.

Recent Advances on the Molecular Mechanism of Cervical Carcinogenesis Based on Systems Biology Technologies

Cervical cancer is one of the common malignancies in women worldwide. Exploration of pathogenesis and molecular mechanism of cervical cancer is pivotal for development of effective treatment for this disease. Recently, systems biology approaches based on high-throughput technologies have been carried out to investigate the expression of some genes and proteins in genomics, transcriptomics, proteomics, and metabonomics of cervical cancer. Compared with traditional methods，systems biology technology has been shown to provide large of information regarding prognostic biomarkers and therapeutic targets for cervical cancer. These molecular signatures from system biology technology could be useful to understand the molecular mechanisms of cervical cancer development and progression, and help physicians to design targeted therapeutic strategies for patients with cervical cancer.