Alternative splicing of human papillomavirus type-16 E6/E6* early mRNA is coupled to EGF signaling via Erk1/2 activation

Regulation of epithelial growth factor receptors by the oncoprotein E5 during the HPV16 differentiation-dependent life cycle

Human papillomavirus (HPV) 16 infection initiates upon viral entry into the basal cells of the epithelium. The virus manipulates signaling pathways to complete its life cycle, which depends on cellular differentiation. The virus expresses the oncoproteins E5, E6, and E7 to promote immune evasion, cell cycle progression, apoptosis inhibition, and viral replication. The least studied viral oncoprotein is E5 (16E5), which can regulate epithelial growth factor receptor (GFR) signaling pathways. GFRs such as transforming growth factor-beta receptor (TGFBR), epidermal growth factor receptor (EGFR), and keratinocyte growth factor receptor (KGFR) have essential roles in cell growth, differentiation, and proliferation. These receptors obtain their ligands from the microenvironment, and once activated, regulate cellular behavior in the epithelium. GFRs therefore represent valuable targets for the virus to establish and maintain a cellular environment supportive of infection. The ability of 16E5 to regulate proliferation and differentiation varies through the differentiating epithelium, making it necessary to adequately describe the association between 16E5 and GFRs. Here we summarize the regulation of GFR signaling pathways by 16E5, discuss the roles of stromal growth factors, and outline unresolved questions over cellular differentiation and proliferation during the HPV life cycle.

Alternative splicing in the genome of HPV and its regulation

Persistent infection with high-risk human papillomavirus (HR-HPV) is the main cause of cervical cancer. These chronic infections are characterized by high expression of the HPV E6 and E7 oncogenes and the absence of the L1 and L2 capsid proteins. The regulation of HPV gene expression plays a crucial role in both the viral life cycle and rare oncogenic events. Alternative splicing of HPV mRNA is a key mechanism in post-transcriptional regulation. Through alternative splicing, HPV mRNA is diversified into various splice isoforms with distinct coding potentials, encoding multiple proteins and influencing the expression of HPV genes. The spliced mRNAs derived from a donor splicing site within the E6 ORF and one of the different acceptor sites located in the early mRNA contain E6 truncated mRNAs, named E6*. E6* is one of the extensively studied splicing isoforms. However, the role of E6* proteins in cancer progression remains controversial. Here, we reviewed and compared the alternative splicing events occurring in the genomes of HR-HPV and LR-HPV. Recently, new HPV alternative splicing regulatory proteins have been continuously discovered, and we have updated the regulation of HPV alternative splicing. In addition, we summarized the functions of known splice isoforms from three aspects: anti-tumorigenic, tumorigenic, and other cancer-related functions, including not only E6*, but also E6^E7, E8^E2, and so on. Comprehending their contributions to cancer development enhances insights into the carcinogenic mechanisms of HPV and explores the potential utility of alternative splicing in the diagnosis and treatment of cervical cancer.

Small molecule NMD and MDM2 inhibitors synergistically trigger apoptosis in HeLa cells

The nonsense-mediated mRNA decay (NMD) pathway and the p53 pathway, linked to tumorgenesis, are also promising targets for cancer treatment. NMD plays an important role in RNA quality control, while the p53 pathway is involved in cancer suppression. However, their individual and combined effects on cervical cancer are poorly understood. In this study, we evaluated the impacts of NMD inhibitor, Mouse double minute 2 homolog (MDM2) inhibitor, and their combination on cell apoptosis, cell cycle, and p53 target genes in human papillomavirus-18-positive HeLa cells. Our findings revealed that XR-2 failed to activate p53 or induce apoptosis in HeLa cells, whereas SMG1 (serine/threonine-protein kinase 1) inhibitor repressed cell proliferation at high concentrations. Notably, the combination of these 2 agents significantly inhibited cell proliferation, arrested the cell cycle, and triggered cell apoptosis. Mechanistically, MDM2 inhibitor and NMD inhibitor likely exert a synergistically through the truncated E6 protein. These results underscore the potential of employing a combination of MDM2 inhibitor and NMD inhibitor as a promising candidate for the clinical treatment of human papillomavirus-infected tumors.

HPV and RNA Binding Proteins: What We Know and What Remains to Be Discovered

Papillomavirus gene regulation is largely post-transcriptional due to overlapping open reading frames and the use of alternative polyadenylation and alternative splicing to produce the full suite of viral mRNAs. These processes are controlled by a wide range of cellular RNA binding proteins (RPBs), including constitutive splicing factors and cleavage and polyadenylation machinery, but also factors that regulate these processes, for example, SR and hnRNP proteins. Like cellular RNAs, papillomavirus RNAs have been shown to bind many such proteins. The life cycle of papillomaviruses is intimately linked to differentiation of the epithelial tissues the virus infects. For example, viral late mRNAs and proteins are expressed only in the most differentiated epithelial layers to avoid recognition by the host immune response. Papillomavirus genome replication is linked to the DNA damage response and viral chromatin conformation, processes which also link to RNA processing. Challenges with respect to elucidating how RBPs regulate the viral life cycle include consideration of the orchestrated spatial aspect of viral gene expression in an infected epithelium and the epigenetic nature of the viral episomal genome. This review discusses RBPs that control viral gene expression, and how the connectivity of various nuclear processes might contribute to viral mRNA production.

Viruses exploit growth factor mechanisms to achieve augmented pathogenicity and promote tumorigenesis

Cellular homeostasis is regulated by growth factors (GFs) which orchestrate various cellular processes including proliferation, survival, differentiation, motility, inflammation and angiogenesis. Dysregulation of GFs in microbial infections and malignancies have been reported previously. Viral pathogens exemplify the exploitation of host cell GFs and their signalling pathways contributing to viral entry, virulence, and evasion of anti-viral immune responses. Viruses can also perturb cellular metabolism and the cell cycle by manipulation of GF signaling. In some cases, this disturbance may promote oncogenesis. Viral pathogens can encode viral GF homologues and induce the endogenous biosynthesis of GFs and their corresponding receptors or manipulate their activity to infect the host cells. Close investigation of how viral strategies exploit and regulate GFs, a will shed light on how to improve anti-viral therapy and cancer treatment. In this review, we discuss and provide insights on how various viral pathogens exploit different GFs to promote viral survival and oncogenic transformation, and how this knowledge can be leveraged toward the design of more efficient therapeutics or novel drug delivery systems in the treatment of both viral infections and malignancies.

RNA elements that control human papillomavirus mRNA splicing-targets for therapy?

Human papillomaviruses (HPVs) cause more than 4.5% of all cancer in the world and more than half of these cases are attributed to human papillomavirus type 16 (HPV16). Prophylactic vaccines are available but antiviral drugs are not. Novel targets for therapy are urgently needed. Alternative RNA splicing is extensively used by HPVs to express all their genes and HPV16 is no exception. This process must function to perfection since mis-splicing could perturb the HPV gene expression program by altering mRNA levels or by generating dysfunctional mRNAs. Cis-acting RNA elements on the viral mRNAs and their cognate cellular trans-acting factors control papillomavirus RNA splicing. The precise but delicate nature of the splicing process renders splicing sensitive to interference. As such, papillomavirus RNA splicing is a potential target for therapy. Here we summarize our current understanding of cis-acting HPV16 RNA elements that control HPV16 mRNA splicing via cellular proteins and discuss how they may be exploited as targets for therapy to papillomavirus infections and cancer.

Differentiated Oral Epithelial Cells Support the HPV Life Cycle

Human Papillomavirus (HPV) associated oral disease continues to increase, both in the context of immune competence and of immune suppression. There are few models of oral HPV infection and current models are laborious. We hypothesized that differentiated oral epithelial cells could support the HPV life cycle. Clinical HPV16 cloned episomes were introduced into differentiated oral epithelial cells (OKF6tert1). Viral and cellular gene expression was assessed in the presence or absence of sodium butyrate, a differentiating agent that moved the cells to full terminal differentiation. Detection of keratin 10, cross-linked involucrin, and loricrin in the presence and absence of sodium butyrate confirmed terminal differentiation. Increasing sodium butyrate concentrations in the absence of HPV, were associated with decreased suprabasal markers and increased terminal differentiation markers. However, in the presence of HPV and of increasing sodium butyrate concentrations, both mitotic and suprabasal markers were increased and the terminal differentiation marker, loricrin, decreased. In this unique differentiated state, early and late viral gene products were detected including spliced mRNAs for E6*, E1^E4, and L1. E7 and L1 proteins were detected. The ratio of late (E1^E4) to early (E6/E7) transcripts in HPV16+ OKF6tert1 cells was distinct compared to HPV16+ C33a cells. Consistent with permissive HPV replication, DNA damage responses (phospho-chk2, gamma-H2AX), HPV E2-dependent LCR transactivation, and DNase-resistant particles were detected and visualized by transmission electron microscopy. In sum, monolayers of differentiated immortalized oral epithelial cells supported the full HPV life cycle. HPV may optimize the differentiation state of oral epithelial cells to facilitate its replication.

hnRNP K induces HPV16 oncogene expression and promotes cervical cancerization

PURPOSE

This study aims to explore the expression of hnRNP K in cervical carcinogenesis and to investigate the regulatory role of hnRNP K on HPV16 oncogene expression as well as biological changes in cervical cancer cells.

METHODS

In total 1042 subjects, including 573 with the normal cervix and 469 with different grades of cervical lesions were enrolled in this study to explore the association between hnRNP K and HPV16 oncogene expression in cervical carcinogenesis. Additionally, the Gene Omnibus (GEO) database was used to analyze hnRNP K mRNA expression in cervical cancerization. Meanwhile, the effects of hnRNP K on cell biological functions and HPV16 oncogene expression were investigated in Siha cells. Moreover, Function analyses were conducted using Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) databases after ChIP-seq.

RESULTS

hnRNP K was highly expressed in cervical cancer and precancerous lesions, and positively correlated with HPV16 E6, but negatively correlated with HPV16 E2 and HPV16 E2/E6 ratio. hnRNP K induced cell proliferation, inhibited apoptosis and caused cell cycle arrest in the S phase, and particularly increased HPV16 E6 protein expression.

CONCLUSION

This study revealed that hnRNP K overexpression has important warning significance for the malignant transformation of cervical lesions, and could be used as a potential therapeutic target for inhibiting the carcinogenicity of HPV16 and prevention of cervical carcinogenesis.

HPV-CCDC106 integration promotes cervical cancer progression by facilitating the high expression of CCDC106 after HPV E6 splicing

Human papillomavirus (HPV) integration and high expression of HPV oncogenes (E6 and E7) are important mechanisms for HPV carcinogenesis in cervical cancer. However, the relationship between HPV integration and HPV E6 spliced transcripts, as well as the underlying mechanisms of HPV integration in carcinogenesis after HPV E6 splicing remains unclear. We analyzed HPV-coiled-coil domain containing 106 (CCDC106) integration samples to characterize the roles of HPV integration, E6 spliceosome I (E6*I), and high CCDC106 expression in cervical carcinogenesis. We found that E6 was alternatively spliced into the E6*I transcript in HPV-CCDC016 integration samples with low p53 expression, in contrast to the role of E6*I in preventing p53 degradation in cervical cancer cells. In addition, CCDC106 was highly expressed after HPV-CCDC106 integration, and interacted with p53, resulting in p53 degradation and cervical cancer cell progression in vitro and in vivo. Importantly, when E6*I was highly expressed in cervical cancer cells, overexpression of CCDC106 independently degraded p53 and promoted cervical cancer cell progression. In this study, we explored the underlying mechanisms of HPV-CCDC106 integration in HPV carcinogenesis after HPV E6 splicing, which should provide insight into host genome dysregulation in cervical carcinogenesis.

Diverse roles of heterogeneous nuclear ribonucleoproteins in viral life cycle

Understanding the host-virus interactions helps to decipher the viral replication strategies and pathogenesis. Viruses have limited genetic content and rely significantly on their host cell to establish a successful infection. Viruses depend on the host for a broad spectrum of cellular RNA-binding proteins (RBPs) throughout their life cycle. One of the major RBP families is the heterogeneous nuclear ribonucleoproteins (hnRNPs) family. hnRNPs are typically localized in the nucleus, where they are forming complexes with pre-mRNAs and contribute to many aspects of nucleic acid metabolism. hnRNPs contain RNA binding motifs and frequently function as RNA chaperones involved in pre-mRNA processing, RNA splicing, and export. Many hnRNPs shuttle between the nucleus and the cytoplasm and influence cytoplasmic processes such as mRNA stability, localization, and translation. The interactions between the hnRNPs and viral components are well-known. They are critical for processing viral nucleic acids and proteins and, therefore, impact the success of the viral infection. This review discusses the molecular mechanisms by which hnRNPs interact with and regulate each stage of the viral life cycle, such as replication, splicing, translation, and assembly of virus progeny. In addition, we expand on the role of hnRNPs in the antiviral response and as potential targets for antiviral drug research and development.

Construction of an immune-related ceRNA network in cervical cancer based on HPV E6 splicing

Background

Cervical cancer is one of the leading causes of cancer-related deaths worldwide. The unspliced human papillomavirus (HPV) E6 plays an important role in tumor progression and immune regulation. Improved immunotherapy implementation might benefit from a better knowledge of HPV E6 splicing-related immune gene expressions and immunocyte infiltration in cervical cancer. This study aimed to identify the potential therapeutic and prognostic roles of unspliced/spliced E6 ratio (E6 ratio) in cervical cancer.

Methods

Data from the TCGA were used to analyze the E6 condition and clinical information. Nomogram and K-M analysis were used to analyze assess the prognostic significance, IOBR was used to investigate immunological infiltrates. Functions and pathway enrichment analysis of DEGs were investigated through GO analysis and KEGG pathway analysis, respectively. A core module was taken from the competitive endogenous RNA (ceRNA) network and used to build a lncRNA-miRNA-mRNA network. QT-qPCR was used to detect the expression of genes. CCK-8, colony formation, wound healing and migration assays were used to detect cell functions.

Results

Our study found that HPV E6 ratio had significantly correlation with overall survival. In cervical cancer, a high E6 ratio was adversely linked with infiltrating levels of aDC, M1 macrophages, monocytes, NKT, and Tgd. High E6 ratio phenotypes were shown to be implicated in immune response regulation, cell adhesion, and Wnt signaling pathways, according to functional enrichment analysis. Subsequently, we constructed an immune-related ceRNA network based on E6 splicing in cervical cancer, including three lncRNA (LINC00943, LIFR-AS1, DANT2, and RASSF8-AS1), four miRNA (miR-205-5p, miR-181d-5p, miR-222-3p, and miR-221-3p), and seven mRNA (FGFR1, PRLR, CXCL2, ISG20, ISG15, SDC1, and NR2F2). Among them, CXCL2, SDC1, and miR-221-3p were associated with survival and immune cell infiltration.

Conclusions

These data imply that a high E6 ratio in cervical cancer contributes to the immune-related ceRNA network, resulting in a low amount of infiltrating effector immune cells and tumor growth. As a result, the E6 ratio might be employed as a biomarker in cervical cancer to determine prognosis and treatment success.

The role of RNA-binding proteins in the processing of mRNAs produced by carcinogenic papillomaviruses

Human papillomaviruses (HPV) are epitheliotropic DNA tumor viruses that are prevalent in the human population. A subset of the HPVs termed high-risk HPVs (HR-HPVs) are causative agents of anogenital cancers and head-and-neck cancers. Cancer is the result of persistent high-risk HPV infections that have not been cleared by the immune system of the host. These infections are characterized by dysregulated HPV gene expression, in particular constitutive high expression of the HPV E6 and E7 oncogenes and absence of the highly immunogenic viral L1 and L2 capsid proteins. HPVs make extensive use of alternative mRNA splicing to express its genes and are therefore highly dependent on cellular RNA-binding proteins for proper gene expression. Levels of RNA-binding proteins are altered in HPV-containing premalignant cervical lesions and in cervical cancer. Here we review our current knowledge of RNA-binding proteins that control HPV gene expression. We focus on RNA-binding proteins that control expression of the E6 and E7 oncogenes since they initiate and drive development of cancer and on the immunogenic L1 and L2 proteins as there silencing may contribute to immune evasion during carcinogenesis. Furthermore, cellular RNA-binding proteins are essential for HPV gene expression and as such may be targets for therapy to HPV infections and HPV-driven cancers.

Emerging roles of hnRNP A2B1 in cancer and inflammation

Heterogeneous nuclear ribonucleoproteins (hnRNPs) are a group of RNA-binding proteins with important roles in multiple aspects of nucleic acid metabolism, including the packaging of nascent transcripts, alternative splicing, transactivation of gene expression, and regulation of protein translation. As a core component of the hnRNP complex in mammalian cells, heterogeneous nuclear ribonucleoprotein A2B1 (hnRNP A2B1) participates in and coordinates various molecular events. Given its regulatory role in inflammation and cancer progression, hnRNP A2B1 has become a novel player in immune response, inflammation, and cancer development. Concomitant with these new roles, a surprising number of mechanisms deemed to regulate hnRNP A2B1 functions have been identified, including post-translational modifications, changes in subcellular localization, direct interactions with multiple DNAs, RNAs, and proteins or the formation of complexes with them, which have gradually made hnRNP A2B1 a molecular target for multiple drugs. In light of the rising interest in the intersection between cancer and inflammation, this review will focus on recent knowledge of the biological roles of hnRNP A2B1 in cancer, immune response, and inflammation.

Therapeutic Vaccines for HPV-Associated Oropharyngeal and Cervical Cancer: The Next De-Intensification Strategy?

The rise in human papillomavirus (HPV)-associated oropharyngeal squamous cell carcinoma (OPSCC) has prompted a quest for further understanding of the role of high-risk HPV in tumor initiation and progression. Patients with HPV-positive OPSCC (HPV+ OPSCC) have better prognoses than their HPV-negative counterparts; however, current therapeutic strategies for HPV+ OPSCC are overly aggressive and leave patients with life-long sequalae and poor quality of life. This highlights a need for customized treatment. Several clinical trials of treatment de-intensification to reduce acute and late toxicity without compromising efficacy have been conducted. This article reviews the differences and similarities in the pathogenesis and progression of HPV-related OPSCC compared to cervical cancer, with emphasis on the role of prophylactic and therapeutic vaccines as a potential de-intensification treatment strategy. Overall, the future development of novel and effective therapeutic agents for HPV-associated head and neck tumors promises to meet the challenges posed by this growing epidemic.

Heterogeneous nuclear ribonucleoprotein A/B: an emerging group of cancer biomarkers and therapeutic targets

Heterogeneous nuclear ribonucleoprotein A/B (hnRNPA/B) is one of the core members of the RNA binding protein (RBP) hnRNPs family, including four main subtypes, A0, A1, A2/B1 and A3, which share the similar structure and functions. With the advance in understanding the molecular biology of hnRNPA/B, it has been gradually revealed that hnRNPA/B plays a critical role in almost the entire steps of RNA life cycle and its aberrant expression and mutation have important effects on the occurrence and progression of various cancers. This review focuses on the clinical significance of hnRNPA/B in various cancers and systematically summarizes its biological function and molecular mechanisms.

The splicing factor SF3B4 drives proliferation and invasion in cervical cancer by regulating SPAG5

Regulation of alternative splicing (AS) by the splicing factor 3b (SF3B) family plays an essential role in cancer. However, the biological function of SF3B family members in cervical cancer (CC) needs to be further elucidated. In this study, we found that splicing factor 3b subunit 4 (SF3B4) was highly expressed in CC by bioinformatics analysis using cervical squamous cell carcinoma and endocervical adenocarcinoma (CESC) data from The Cancer Genome Atlas (TCGA). Then, we demonstrated that high expression of SF3B4 promoted proliferation and invasion abilities of CC cells in vitro and in vivo and that reduced expression of SF3B4 performed the opposite effect. Further RNA-seq and AS analysis showed that sperm-associated antigen 5 (SPAG5) was a downstream target gene of SF3B4. Interestingly, SPAG5 expression was decreased after SF3B4 knockdown because of retained introns (RIs) and reduced maturation of SPAG5 pre-mRNA. Importantly, SPAG5 deficiency impaired the oncogenic effects of SF3B4 overexpression on CC cells. In conclusion, SF3B4 promotes CC progression by regulating the effective splicing of SPAG5. SF3B4 could be a promising target for CC.

High-Risk Human Papillomavirus Oncogenic E6/E7 mRNAs Splicing Regulation

High-risk human papillomavirus infection may develop into a persistent infection that is highly related to the progression of various cancers, including cervical cancer and head and neck squamous cell carcinomas. The most common high-risk subtypes are HPV16 and HPV18. The oncogenic viral proteins expressed by high-risk HPVs E6/E7 are tightly involved in cell proliferation, differentiation, and cancerous transformation since E6/E7 mRNAs are derived from the same pre-mRNA. Hence, the alternative splicing in the E6/E7-coding region affects the balance of the E6/E7 expression level. Interrupting the balance of E6 and E7 levels results in cell apoptosis. Therefore, it is crucial to understand the regulation of E6/E7 splice site selection and the interaction of splicing enhancers and silencers with cellular splicing factors. In this review, we concluded the relationship of different E6/E7 transcripts with cancer progression, the known splicing sites, and the identified cis-regulatory elements within high-risk HPV E6/E7-coding region. Finally, we also reviewed the role of various splicing factors in the regulation of high-risk HPV oncogenic E6/E7 mRNA splicing.

It Takes Two to Tango: A Review of Oncogenic Virus and Host Microbiome Associated Inflammation in Head and Neck Cancer

Simple Summary

Certain viruses, specifically, human papillomavirus (HPV) and Epstein–Barr virus (EBV), have been linked with the development of head and neck cancer. In this study, we review the mechanisms by which (these) viruses lead to cellular transformation and a chronic inflammatory state. Given that the head and neck host a rich microbiome (which itself is intrinsically linked to inflammation), we scrutinize the literature to highlight the interplay between viruses, cellular transformation, inflammation, and the local host microbiome in head and neck cancer.

Abstract

While the two primary risk factors for head and neck squamous cell carcinoma (HNSCC) are alcohol and tobacco, viruses account for an important and significant upward trend in HNSCC incidence. Human papillomavirus (HPV) is the causative agent for a subset of oropharyngeal squamous cell carcinoma (OPSCC)—a cancer that is impacting a rapidly growing group of typically middle-aged non-smoking white males. While HPV is a ubiquitously present (with about 1% of the population having high-risk oral HPV infection at any one time), less than 1% of those infected with high-risk strains develop OPSCC—suggesting that additional cofactors or coinfections may be required. Epstein–Barr virus (EBV) is a similarly ubiquitous virus that is strongly linked to nasopharyngeal carcinoma (NPC). Both of these viruses cause cellular transformation and chronic inflammation. While dysbiosis of the human microbiome has been associated with similar chronic inflammation and the pathogenesis of mucosal diseases (including OPSCC and NPC), a significant knowledge gap remains in understanding the role of bacterial-viral interactions in the initiation, development, and progression of head and neck cancers. In this review, we utilize the known associations of HPV with OPSCC and EBV with NPC to investigate these interactions. We thoroughly review the literature and highlight how perturbations of the pharyngeal microbiome may impact host-microbiome-tumor-viral interactions—leading to tumor growth.

HPV16 and HPV18 Genome Structure, Expression, and Post-Transcriptional Regulation

Human papillomaviruses (HPV) are a group of small non-enveloped DNA viruses whose infection causes benign tumors or cancers. HPV16 and HPV18, the two most common high-risk HPVs, are responsible for ~70% of all HPV-related cervical cancers and head and neck cancers. The expression of the HPV genome is highly dependent on cell differentiation and is strictly regulated at the transcriptional and post-transcriptional levels. Both HPV early and late transcripts differentially expressed in the infected cells are intron-containing bicistronic or polycistronic RNAs bearing more than one open reading frame (ORF), because of usage of alternative viral promoters and two alternative viral RNA polyadenylation signals. Papillomaviruses proficiently engage alternative RNA splicing to express individual ORFs from the bicistronic or polycistronic RNA transcripts. In this review, we discuss the genome structures and the updated transcription maps of HPV16 and HPV18, and the latest research advances in understanding RNA cis-elements, intron branch point sequences, and RNA-binding proteins in the regulation of viral RNA processing. Moreover, we briefly discuss the epigenetic modifications, including DNA methylation and possible APOBEC-mediated genome editing in HPV infections and carcinogenesis.

hnRNP G/RBMX enhances HPV16 E2 mRNA splicing through a novel splicing enhancer and inhibits production of spliced E7 oncogene mRNAs

Human papillomavirus type 16 (HPV16) E2 is an essential HPV16 protein. We have investigated how HPV16 E2 expression is regulated and have identifed a splicing enhancer that is required for production of HPV16 E2 mRNAs. This uridine-less splicing enhancer sequence (ACGAGGACGAGGACAAGGA) contains 84% adenosine and guanosine and 16% cytosine and consists of three ‘AC(A/G)AGG’-repeats. Mutational inactivation of the splicing enhancer reduced splicing to E2-mRNA specific splice site SA2709 and resulted in increased levels of unspliced E1-encoding mRNAs. The splicing enhancer sequence interacted with cellular RNA binding protein hnRNP G that promoted splicing to SA2709 and enhanced E2 mRNA production. The splicing-enhancing function of hnRNP G mapped to amino acids 236–286 of hnRNP G that were also shown to interact with splicing factor U2AF65. The interactions between hnRNP G and HPV16 E2 mRNAs and U2AF65 increased in response to keratinocyte differentiation as well as by the induction of the DNA damage response (DDR). The DDR reduced sumoylation of hnRNP G and pharmacological inhibition of sumoylation enhanced HPV16 E2 mRNA splicing and interactions between hnRNP G and E2 mRNAs and U2AF65. Intriguingly, hnRNP G also promoted intron retention of the HPV16 E6 coding region thereby inhibiting production of spliced E7 oncogene mRNAs.

Overexpression of m6A-factors METTL3, ALKBH5, and YTHDC1 alters HPV16 mRNA splicing

We report that overexpression of the m6A-demethylase alkB homolog 5 RNA demethylase (ALKBH5) promoted production of intron retention on the human papillomavirus type 16 (HPV16) E6 mRNAs thereby promoting E6 mRNA production. ALKBH5 also altered alternative splicing of the late L1 mRNA by an exon skipping mechanism. Knock-down of ALKBH5 had the opposite effect on splicing of these HPV16 mRNAs. Overexpression of the m6A-methylase methyltransferase-like protein 3 (METLL3) induced production of intron-containing HPV16 E1 mRNAs over spliced E2 mRNAs and altered HPV16 L1 mRNA splicing in a manner opposite to ALKBH5. Overexpression of the nuclear m6A-“reader” YTH domain-containing protein 1 (YTHDC1), enhanced retention of the E6-encoding intron and promoted E6 mRNA production. We also show that HPV16 mRNAs are bound to YTHDC1 in human cells and that YTHDC1 affected splicing of HPV16 E6/E7 mRNAs produced from the episomal form of the HPV16 genome. Finally, we show that HPV16 mRNAs are m6A-methylated in tonsillar cancer cells. In summary, HPV16 mRNAs are methylated in HPV16-infected tonsillar cancer cells and overexpression of m6A-“writer” METTL3, m6A-“eraser” ALKBH5 and the m6A-“reader” YTHDC1 affected HPV16 mRNA splicing, suggesting that m6A plays an important role in the HPV16 gene expression program, at least in cancer cells.

Identification of heterogenous nuclear ribonucleoproteins (hnRNPs) and serine- and arginine-rich (SR) proteins that induce human papillomavirus type 16 late gene expression and alter L1 mRNA splicing

We have determined the effect of seven serine- and arginine-rich (SR) proteins and 15 heterogenous nuclear ribonucleoproteins (hnRNPs) on human papillomavirus type 16 (HPV16) late gene expression. Of the seven SR proteins analyzed here, SRSF1, SRSF3, and SRSF9 induced HPV16 late gene expression, and five of the SR proteins affected HPV16 L1 mRNA splicing. Of the 15 hnRNP proteins analyzed here, hnRNP A2, hnRNP F, and hnRNP H efficiently induced HPV16 late gene expression, and all of the hnRNPs affected HPV16 L1 mRNA levels or mRNA splicing. Thus, the majority of SR proteins and hnRNPs have the potential to regulate HPV16 L1 mRNA splicing. Strict control of the expression of the immunogenic L1 and L2 capsid proteins may contribute to the ability of HPV16 to cause persistence.

Differential Incidence of Tongue Base Cancer in Male and Female HPV16-Transgenic Mice: Role of Female Sex Hormone Receptors

A growing proportion of oropharyngeal squamous cell carcinomas (OPSCC) are associated with infection by high-risk human papillomavirus (HPV). For reasons that remain largely unknown, HPV+OPSCC is significantly more common in men than in women. This study aims to determine the incidence of OPSCC in male and female HPV16-transgenic mice and to explore the role of female sex hormone receptors in the sexual predisposition for HPV+ OPSCC. The tongues of 30-weeks-old HPV16-transgenic male (n = 80) and female (n = 90) and matched wild-type male (n = 10) and female (n = 10) FVB/n mice were screened histologically for intraepithelial and invasive lesions in 2017 at the Centre for the Research and Technology of Agro-Environmental and Biological Sciences (CITAB), Portugal. Expression of estrogen receptors alpha (ERα) and beta (ERβ), progesterone receptors (PR) and matrix metalloproteinase 2 (MMP2) was studied immunohistochemically. Collagen remodeling was studied using picrosirius red. Female mice showed robust ERα and ERβ expression in intraepithelial and invasive lesions, which was accompanied by strong MMP2 expression and marked collagen remodeling. Male mice showed minimal ERα, ERβ and MMP2 expression and unaltered collagen patterns. These results confirm the association of HPV16 with tongue base cancer in both sexes. The higher cancer incidence in female versus male mice contrasts with data from OPSCC patients and is associated with enhanced ER expression via MMP2 upregulation.

HPV transcript expression affects cervical cancer response to chemoradiation

Persistent HPV infection is causative for the majority of cervical cancer cases; however, current guidelines do not require HPV testing for newly diagnosed cervical cancer. Using an institutional cohort of 88 patients with cervical cancer treated uniformly with standard-of-care chemoradiation treatment (CRT) with prospectively collected clinical outcome data, we observed that patients with cervical tumors containing HPV genotypes other than HPV 16 have worse survival outcomes after CRT compared with patients with HPV 16⁺ tumors, consistent with previously published studies. Using RNA sequencing analysis, we quantified viral transcription efficiency and found higher levels of E6 and the alternative transcript E6*I in cervical tumors with HPV genotypes other than HPV 16. These findings were validated using whole transcriptome data from The Cancer Genome Atlas (n = 304). For the first time to our knowledge, transcript expression level of HPV E6*I was identified as a predictive biomarker of CRT outcome in our complete institutional data set (n = 88) and within the HPV 16⁺ subset (n = 36). In vitro characterization of HPV E6*I and E6 overexpression revealed that both induce CRT resistance through distinct mechanisms dependent upon p53–p21. Our findings suggest that high expression of E6*I and E6 may represent novel biomarkers of CRT efficacy, and these patients may benefit from alternative treatment strategies.

MEK/ERK signaling is a critical regulator of high-risk human papillomavirus oncogene expression revealing therapeutic targets for HPV-induced tumors

Intracellular pathogens have evolved to utilize normal cellular processes to complete their replicative cycles. Pathogens that interface with proliferative cell signaling pathways risk infections that can lead to cancers, but the factors that influence malignant outcomes are incompletely understood. Human papillomaviruses (HPVs) predominantly cause benign hyperplasia in stratifying epithelial tissues. However, a subset of carcinogenic or “high-risk” HPV (hr-HPV) genotypes are etiologically linked to nearly 5% of all human cancers. Progression of hr-HPV-induced lesions to malignancies is characterized by increased expression of the E6 and E7 oncogenes and the oncogenic functions of these viral proteins have been widely studied. Yet, the mechanisms that regulate hr-HPV oncogene transcription and suppress their expression in benign lesions remain poorly understood. Here, we demonstrate that EGFR/MEK/ERK signaling, influenced by epithelial contact inhibition and tissue differentiation cues, regulates hr-HPV oncogene expression. Using monolayer cells, epithelial organotypic tissue models, and neoplastic tissue biopsy materials, we show that cell-extrinsic activation of ERK overrides cellular control to promote HPV oncogene expression and the neoplastic phenotype. Our data suggest that HPVs are adapted to use the EGFR/MEK/ERK signaling pathway to regulate their productive replicative cycles. Mechanistic studies show that EGFR/MEK/ERK signaling influences AP-1 transcription factor activity and AP-1 factor knockdown reduces oncogene transcription. Furthermore, pharmacological inhibitors of EGFR, MEK, and ERK signaling quash HPV oncogene expression and the neoplastic phenotype, revealing a potential clinical strategy to suppress uncontrolled cell proliferation, reduce oncogene expression and treat HPV neoplasia.

Human papillomavirus (HPV) infections occur in differentiating squamous epithelium and induce hyperplasia during the viral replicative cycle. Although HPV oncogene expression is necessary to promote cellular proliferation for viral genome amplification in the middle epithelial layers, oncogene levels are thereafter suppressed to permit differentiation-induced late gene expression in the uppermost epithelial cells. Yet, the mechanisms responsible for controlling HPV oncogene expression are not well understood. Here, we demonstrate that EGFR/MEK/ERK signaling, which is subject to the normal cellular cues of contact inhibition and epithelial tissue differentiation, is a critical regulator of hr-HPV oncogene expression. We found that extrinsic activation of ERK overrides cellular control to promote oncogene expression and the neoplastic phenotype. Many epidemiologically defined risk factors activate the EGFR/MEK/ERK pathway, suggesting a common mechanism whereby they may promote HPV persistence and disease progression. Lastly, we show that HPV oncogene transcription and protein expression remain susceptible to MEK/ERK control in early neoplastic tissues and tumor cells and that targeted inhibition of MEK/ERK signaling might be exploited therapeutically for HPV-induced infections and tumors.

Modeling the epithelial-mesenchymal transition process in a 3D organotypic cervical neoplasia

Here, we proposed an innovative organotypic cervical tumor model able to investigate the bi-directional crosstalk between epithelium and stroma as well as the key disease features of the epithelial-mesenchymal transition (EMT) process in vitro. By using a modular tissue assembling approach, we developed 3D cervical stromal models composed of primary human cervical fibroblasts (HCFs) or cervical cancer-associated fibroblasts (CCAFs) embedded in their own ECM to produce 3D normal cervical-instructed stroma (NCIS) or 3D cervical cancer-instructed stroma (CCIS), respectively. Then, we demonstrate the role of the tumor microenvironment (TME) in potentiating the intrinsic invasive attitude of cervical cancer derived SiHa cells and increasing their early viral gene expression by comparing the SiHa behavior when cultured on NCIS or CCIS (SiHa-NCIS or SiHa-CCIS). We proved the crucial role of the CCAFs and stromal microenvironment in the mesenchymalization of the cancer epithelial cells by analyzing several EMT markers. We further assessed the expression of the epithelial adhesion molecules, matricellular enzymes, non-collagenous proteins as well as ECM remodeling in terms of collagen fibers texture and assembly. This cervical tumor model, closely recapitulating key cervical carcinogenesis features, may provide efficient and relevant support to current approaches characterizing cancer progression and develop new anticancer therapy targeting stroma rather than cancer cells.

Role of Viral Ribonucleoproteins in Human Papillomavirus Type 16 Gene Expression

Human papillomaviruses (HPVs) depend on the cellular RNA-processing machineries including alternative RNA splicing and polyadenylation to coordinate HPV gene expression. HPV RNA processing is controlled by cis-regulatory RNA elements and trans-regulatory factors since the HPV splice sites are suboptimal. The definition of HPV exons and introns may differ between individual HPV mRNA species and is complicated by the fact that many HPV protein-coding sequences overlap. The formation of HPV ribonucleoproteins consisting of HPV pre-mRNAs and multiple cellular RNA-binding proteins may result in the different outcomes of HPV gene expression, which contributes to the HPV life cycle progression and HPV-associated cancer development. In this review, we summarize the regulation of HPV16 gene expression at the level of RNA processing with focus on the interactions between HPV16 pre-mRNAs and cellular RNA-binding factors.

Heterogeneous Nuclear Ribonucleoprotein A1 (hnRNP A1) and hnRNP A2 Inhibit Splicing to Human Papillomavirus 16 Splice Site SA409 through a UAG-Containing Sequence in the E7 Coding Region

Human papillomavirus type 16 (HPV16) belongs to the high-risk-group of HPVs and is causing a variety of anogenital cancers and head and neck cancer. The two HPV16 oncoproteins E6 and E7 prevent apoptosis and promote mitosis and are essential for completion of the HPV16 life cycle and for transformation of the infected cell and maintenance of malignancy. E6 and E7 are produced from two mRNAs that are generated in a mutually exclusive manner by alternative splicing. While E6 protein is made from the unspliced mRNA, E7 is made from the spliced version of the same pre-mRNA. Since sufficient quantities of both E6 and E7 are required for malignant transformation, this intricate arrangement of gene expression renders E6 and E7 expression vulnerable to external interference. Since antiviral drugs to HPV16 are not available, a detailed knowledge of the regulation of HPV16 E6 and E7 mRNA splicing may uncover novel targets for therapy.

Human papillomavirus 16 (HPV16) 5′-splice site SD226 and 3′-splice site SA409 are required for production of the HPV16 E7 mRNAs, whereas unspliced mRNAs produce E6 mRNAs. The E6 and E7 proteins are essential in the HPV16 replication cycle but are also the major HPV16 proteins required for induction and maintenance of malignancy caused by HPV16 infection. Thus, a balanced expression of unspliced and spliced mRNAs is required for production of sufficient quantities of E6 and E7 proteins under physiological and pathophysiological conditions. If splicing becomes too efficient, the levels of unspliced E6 mRNAs will decrease below a threshold level that is no longer able to produce E6 protein quantities high enough to significantly reduce p53 protein levels. Similarly, if splicing becomes too inefficient, the levels of spliced E7 mRNAs will decrease below a threshold level that is no longer able to produce E7 protein quantities high enough to significantly reduce pRb protein levels. To determine how splicing between SD226 and SA409 is regulated, we have investigated how SA409 is controlled by the cellular proteins hnRNP A1 and hnRNP A2, two proteins that have been shown previously to control HPV16 gene expression. We found that hnRNP A1 and A2 interacted directly and specifically with a C-less RNA element located between HPV16 nucleotide positions 594 and 604 downstream of SA409. Overexpression of hnRNP A1 inhibited SA409 and promoted production of unspliced E6 mRNAs at the expense of the E7 mRNAs, whereas overexpression of hnRNP A2 inhibited SA409 to redirect splicing to SA742, a downstream 3′-splice site that is used for generation of HPV16 E6̂E7, E1, and E4 mRNAs. Thus, high levels of either hnRNP A1 or hnRNP A2 inhibited production of the promitotic HPV16 E7 protein. We show that the hnRNP A1 and A2 proteins control the relative levels of the HPV16 unspliced and spliced HPV16 E6 and E7 mRNAs and function as inhibitors of HPV16 E7 expression.

IMPORTANCE Human papillomavirus type 16 (HPV16) belongs to the high-risk-group of HPVs and is causing a variety of anogenital cancers and head and neck cancer. The two HPV16 oncoproteins E6 and E7 prevent apoptosis and promote mitosis and are essential for completion of the HPV16 life cycle and for transformation of the infected cell and maintenance of malignancy. E6 and E7 are produced from two mRNAs that are generated in a mutually exclusive manner by alternative splicing. While E6 protein is made from the unspliced mRNA, E7 is made from the spliced version of the same pre-mRNA. Since sufficient quantities of both E6 and E7 are required for malignant transformation, this intricate arrangement of gene expression renders E6 and E7 expression vulnerable to external interference. Since antiviral drugs to HPV16 are not available, a detailed knowledge of the regulation of HPV16 E6 and E7 mRNA splicing may uncover novel targets for therapy.

An Immunocompetent Mouse Model of HPV16(+) Head and Neck Squamous Cell Carcinoma

The incidence of human papilloma virus (HPV)-associated head and neck squamous cell carcinoma (HNSCC) is increasing and implicated in more than 60% of all oropharyngeal carcinomas (OPSCCs). Although whole-genome, transcriptome, and proteome analyses have identified altered signaling pathways in HPV-induced HNSCCs, additional tools are needed to investigate the unique pathobiology of OPSCC. Herein, bioinformatics analyses of human HPV(+) HNSCCs revealed that all tumors express full-length E6 and identified molecular subtypes based on relative E6 and E7 expression levels. To recapitulate the levels, stoichiometric ratios, and anatomic location of E6/E7 expression, we generated a genetically engineered mouse model whereby balanced expression of E6/E7 is directed to the oropharyngeal epithelium. The addition of a mutant PIK3CA^E545K allele leads to the rapid development of pre-malignant lesions marked by immune cell accumulation, and a subset of these lesions progress to OPSCC. This mouse provides a faithful immunocompetent model for testing treatments and investigating mechanisms of immuno- suppression.

Carper et al. present the ‘‘iKHP’’ mouse, in which HPV16 oncogenes are inducibly activated in vivo in a tissue-specific and temporal manner. Oropharyngeal- specific expression of E6/E7 with PIK3CA^E545K in these mice promotes the development of premalignant lesions marked by immune cell infiltration, but only a subset spontaneously convert to OPSCC.

HPV-EM: an accurate HPV detection and genotyping EM algorithm

Accurate HPV genotyping is crucial in facilitating epidemiology studies, vaccine trials, and HPV-related cancer research. Contemporary HPV genotyping assays only detect < 25% of all known HPV genotypes and are not accurate for low-risk or mixed HPV genotypes. Current genomic HPV genotyping algorithms use a simple read-alignment and filtering strategy that has difficulty handling repeats and homology sequences. Therefore, we have developed an optimized expectation–maximization algorithm, designated HPV-EM, to address the ambiguities caused by repetitive sequencing reads. HPV-EM achieved 97–100% accuracy when benchmarked using cell line data and TCGA cervical cancer data. We also validated HPV-EM using DNA tiling data on an institutional cervical cancer cohort (96.5% accuracy). Using HPV-EM, we demonstrated HPV genotypic differences in recurrence and patient outcomes in cervical and head and neck cancers.

The Role of RNA Splicing Factors in Cancer: Regulation of Viral and Human Gene Expression in Human Papillomavirus-Related Cervical Cancer

The spliceosomal complex components, together with the heterogeneous nuclear ribonucleoproteins (hnRNPs) and serine/arginine-rich (SR) proteins, regulate the process of constitutive and alternative splicing, the latter leading to the production of mRNA isoforms coding multiple proteins from a single pre-mRNA molecule. The expression of splicing factors is frequently deregulated in different cancer types causing the generation of oncogenic proteins involved in cancer hallmarks. Cervical cancer is caused by persistent infection with oncogenic human papillomaviruses (HPVs) and constitutive expression of viral oncogenes. The aberrant activity of hnRNPs and SR proteins in cervical neoplasia has been shown to trigger the production of oncoproteins through the processing of pre-mRNA transcripts either derived from human genes or HPV genomes. Indeed, hnRNP and SR splicing factors have been shown to regulate the production of viral oncoprotein isoforms necessary for the completion of viral life cycle and for cell transformation. Target-therapy strategies against hnRNPs and SR proteins, causing simultaneous reduction of oncogenic factors and inhibition of HPV replication, are under development. In this review, we describe the current knowledge of the functional link between RNA splicing factors and deregulated cellular as well as viral RNA maturation in cervical cancer and the opportunity of new therapeutic strategies.

Significant association between host transcriptome-derived HPV oncogene E6* influence score and carcinogenic pathways, tumor size, and survival in head and neck cancer

BACKGROUND

Human papillomavirus (HPV) oncogenes E6, E7, and shorter isoforms of E6 (E6*) are known carcinogenic factors in head and neck squamous cell carcinoma (HNSCC). Little is known regarding E6* functions.

METHODS

We analyzed RNA-seq data from 68 HNSCC HPV type 16-positive tumors to determine host genes and pathways associated with E6+E7 expression (E6E7) or the percent of full-length E6 (E6%FL). Influence scores of E6E7 and E6%FL were used to test for associations with clinical variables.

RESULTS

For E6E7, we recapitulated all major known affected pathways and revealed additional pathways. E6%FL was found to affect mitochondrial processes, and E6%FL influence score was significantly associated with overall survival and tumor size.

CONCLUSIONS

HPV E6E7 and E6* result in extensive, dose-dependent compensatory effects and dysregulation of key cancer pathways. The switch from E6 to E6* promotes oxidative phosphorylation, larger tumor size, and worse prognosis, potentially serving as a prognostic factor for HPV-positive HNSCC.

The multifarious roles of heterogeneous ribonucleoprotein A1 in viral infections

Viruses are obligate parasites known to interact with a wide variety of host proteins at different stages of infection. Current antiviral treatments target viral proteins and may be compromised due to the emergence of drug resistant viral strains. Targeting viral-host interactions is now gaining recognition as an alternative approach against viral infections. Recent research has revealed that heterogeneous ribonucleoprotein A1, an RNA-binding protein, plays an essential functional and regulatory role in the life cycle of many viruses. In this review, we summarize the interactions between heterogeneous ribonucleoprotein A1 (hnRNPA1) and multiple viral proteins during the life cycle of RNA and DNA viruses. hnRNPA1 protein levels are modulated differently, in different viruses, which further dictates its stability, function, and intracellular localization. Multiple reports have emphasized that in Sindbis virus, enteroviruses, porcine endemic diarrhea virus, and rhinovirus infection, hnRNPA1 enhances viral replication and survival. However, in others like hepatitis C virus and human T-cell lymphotropic virus, it exerts a protective response. The involvement of hnRNPA1 in viral infections highlights its importance as a central regulator of host and viral gene expression. Understanding the nature of these interactions will increase our understanding of specific viral infections and pathogenesis and eventually aid in the development of novel and robust antiviral intervention strategies.

E6/E7 and E6* From HPV16 and HPV18 Upregulate IL-6 Expression Independently of p53 in Keratinocytes

Keratinocyte infection with high-risk human papillomavirus genotypes has been linked to cancer development. In cervix, the alpha HPV16 and HPV18 have been reported as the mayor causative agents of cervical cancer. Oncogenic progression and chronic inflammation are closely related processes, with IL-6 as one of the main pro-inflammatory cytokines involved. However, there are limited studies about the regulation of IL-6 by low and high risk HPVs and the HPV proteins implicated in this modulation. In this work, we report the overexpression of IL-6 in HPV infected cervical cancer derived cell lines (HeLa and SiHa) compared to non-tumorigenic keratinocytes (HaCaT), and in Cervical Intraepithelial Neoplasia grade 1 HPV16 and HPV18 positive cervical samples compared to HPV negative samples without lesions. Moreover, we generated HaCaT keratinocytes that express E5, E6, and E7 from high risk (16 or 18) or low risk (62 and 84) HPVs. E5 proteins do not modify IL-6 expression, while E7 modestly increase it. Interestingly, E6 proteins in HaCaT cells upregulate IL-6 mRNA expression and protein secretion. Indeed, in HaCaT cells that express high risk HPV16E6 or HPV18E6 proteins, only the truncated E6^* isoforms were expressed, showing the stronger IL-6 overexpression, while in HaCaT cells that express low risk HPV62 and HPV84 full length E6 proteins, IL-6 was also upregulated but not so drastically. Since HaCaT cells have a mutated p53 form that is not degraded by the introduction of E6 or E6/E7, it seems that E6/E7 regulate IL-6 by an additional mechanism independent of p53. In addition, basal keratinocytes showed a strong expression of IL-6R using immunohistochemistry, suggesting an autocrine mechanism over proliferative cells. Altogether, IL-6 cytokine expression in keratinocytes is upregulated by E6 and E7 proteins from HPVs 16, 18, 62, and 84, especially by high risk HPV16 and HPV18 E6^*, which may contribute to promote a pro-inflammatory and highly proliferative microenvironment that can persist over time and lead to cervical tumorigenesis.

Human DNA Virus Exploitation of the MAPK-ERK Cascade

The extracellular signal-regulated kinases (ERKs) comprise a particular branch of the mitogen-activated protein kinase cascades (MAPK) that transmits extracellular signals into the intracellular environment to trigger cellular growth responses. Similar to other MAPK cascades, the MAPK-ERK pathway signals through three core kinases—Raf, MAPK/ERK kinase (MEK), and ERK—which drive the signaling mechanisms responsible for the induction of cellular responses from extracellular stimuli including differentiation, proliferation, and cellular survival. However, pathogens like DNA viruses alter MAPK-ERK signaling in order to access DNA replication machineries, induce a proliferative state in the cell, or even prevent cell death mechanisms in response to pathogen recognition. Differential utilization of this pathway by multiple DNA viruses highlights the dynamic nature of the MAPK-ERK pathway within the cell and the importance of its function in regulating a wide variety of cellular fates that ultimately influence viral infection and, in some cases, result in tumorigenesis.

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.

MiR-875 and miR-3144 switch the human papillomavirus 16 E6/E6* mRNA ratio through the EGFR pathway and a direct targeting effect

By employing bioinformatics scanning approaches and luciferase reporter, our previous study showed that two less common human miRNAs, miR-875 and miR-3144, target a conserved site in the genomes of most high-risk human papillomaviruses (HR-HPVs). In this study, we found that the site targeted by miR-875 and miR-3144 overlapped with the 5' alternative splice site of E6E7 transcripts in HPV16. Using HPV16+ SiHa cells, we showed that high levels of miR-875 and miR-3144 reduced the abundance of unspliced E6, while they promoted three E6* spliced transcripts and decreased the expression levels of E6/E7 oncoproteins and epidermal growth factor receptor (EGFR). A potential miR-875 target site was predicted in EGFR. Meanwhile, depletion of EGFR resulted in a failure to promote E6* but maintained the suppression of unspliced E6 driven by miR-875 and miR-3144. The data suggest that these two miRNAs switch the E6/E6* ratio through both the EGFR pathway and direct targeting. Here, we demonstrate for the first time that human miRNAs regulate the HPV splice isoforms. Furthermore, miRNA-875 and miRNA-3144 are only found in vertebrates and Homo sapiens, and the binding site in EGFR is highly conserved in Boreoeutheria. Our findings highlight the tumour-suppressing effect of miRNAs that possibly appeared in the late stage of biological evolution.

ABBREVIATIONS

E5 can be expressed in anal cancer and leads to epidermal growth factor receptor-induced invasion in a human papillomavirus 16-transformed anal epithelial cell line

We created the first human papillomavirus (HPV)-16-immortalized anal epithelial cell line, known as AKC2 cells to establish an in vitro model of HPV-16-induced anal carcinogenesis. Consistent with detection of E6, E7 and E5 expression in anal cancer biopsies, AKC2 cells expressed high levels of all three HPV oncogenes. Also, similar to findings in anal cancer biopsies, epidermal growth factor receptor (EGFR) was overexpressed in AKC2 cells. AKC2 cells exhibited a poorly differentiated and invasive phenotype in three-dimensional raft culture and inhibition of EGFR function abrogated AKC2 invasion. Reducing E5 expression using E5-targeted siRNAs in AKC2 cells led to knockdown of E5 expression, but also HPV-16 E2, E6 and E7 expression. AKC2 cells treated with E5-targeted siRNA had reduced levels of total and phosphorylated EGFR, and reduced invasion. Rescue of E6/E7 expression with simultaneous E5 knockdown confirmed that E5 plays a key role in EGFR overexpression and EGFR-induced invasion.

The Role of E6 Spliced Isoforms (E6*) in Human Papillomavirus-Induced Carcinogenesis

Persistent infections with High Risk Human Papillomaviruses (HR-HPVs) are the main cause of cervical cancer development. The E6 and E7 oncoproteins of HR-HPVs are derived from a polycistronic pre-mRNA transcribed from an HPV early promoter. Through alternative splicing, this pre-mRNA produces a variety of E6 spliced transcripts termed E6*. In pre-malignant lesions and HPV-related cancers, different E6/E6* transcriptional patterns have been found, although they have not been clearly associated to cancer development. Moreover, there is a controversy about the participation of E6* proteins in cancer progression. This review addresses the regulation of E6 splicing and the different functions that have been found for E6* proteins, as well as their possible role in HPV-induced carcinogenesis.

Human papillomavirus oncogenic E6 protein regulates human β-defensin 3 (hBD3) expression via the tumor suppressor protein p53

Human β-defensin-3 (hBD3) is an epithelial cell-derived innate immune regulatory molecule overexpressed in oral dysplastic lesions and fosters a tumor-promoting microenvironment. Expression of hBD3 is induced by the epidermal growth factor receptor signaling pathway. Here we describe a novel pathway through which the high-risk human papillomavirus type-16 (HPV-16) oncoprotein E6 induces hBD3 expression in mucosal keratinocytes. Ablation of E6 by siRNA induces the tumor suppressor p53 and diminishes hBD3 in HPV-16 positive CaSki cervical cancer cells and UM-SCC-104 head and neck cancer cells. Malignant cells in HPV-16-associated oropharyngeal cancer overexpress hBD3. HPV-16 E6 induces hBD3 mRNA expression, peptide production and gene promoter activity in mucosal keratinocytes. Reduction of cellular levels of p53 stimulates hBD3 expression, while activation of p53 by doxorubicin inhibits its expression in primary oral keratinocytes and CaSki cells, suggesting that p53 represses hBD3 expression. A p53 binding site in the hBD3 gene promoter has been identified by using electrophoretic mobility shift assays and chromatin immunoprecipitation (ChIP). In addition, the p63 protein isoform ΔNp63α, but not TAp63, stimulated transactivation of the hBD3 gene and was co-expressed with hBD3 in head and neck cancer specimens. Therefore, high-risk HPV E6 oncoproteins may stimulate hBD3 expression in tumor cells to facilitate tumorigenesis of HPV-associated head and neck cancer.

hnRNP L controls HPV16 RNA polyadenylation and splicing in an Akt kinase-dependent manner

Inhibition of the Akt kinase activates HPV16 late gene expression by reducing HPV16 early polyadenylation and by activating HPV16 late L1 mRNA splicing. We identified ‘hot spots’ for RNA binding proteins at the early polyA signal and at splice sites on HPV16 late mRNAs. We observed that hnRNP L was associated with sequences at all HPV16 late splice sites and at the early polyA signal. Akt kinase inhibition resulted in hnRNP L dephosphorylation and reduced association of hnRNP L with HPV16 mRNAs. This was accompanied by an increased binding of U2AF65 and Sam68 to HPV16 mRNAs. Furthermore, siRNA knock-down of hnRNP L or Akt induced HPV16 gene expression. Treatment of HPV16 immortalized keratinocytes with Akt kinase inhibitor reduced hnRNP L binding to HPV16 mRNAs and induced HPV16 L1 mRNA production. Finally, deletion of the hnRNP L binding sites in HPV16 subgenomic expression plasmids resulted in activation of HPV16 late gene expression. In conclusion, the Akt kinase inhibits HPV16 late gene expression at the level of RNA processing by controlling the RNA-binding protein hnRNP L. We speculate that Akt kinase activity upholds an intracellular milieu that favours HPV16 early gene expression and suppresses HPV16 late gene expression.

Modulation of basal cell fate during productive and transforming HPV-16 infection is mediated by progressive E6-driven depletion of Notch

In stratified epithelia such as the epidermis, homeostasis is maintained by the proliferation of cells in the lower epithelial layers and the concomitant loss of differentiated cells from the epithelial surface. These differentiating keratinocytes progressively stratify and form a self‐regenerating multi‐layered barrier that protects the underlying dermis. In such tissue, the continual loss and replacement of differentiated cells also limits the accumulation of oncogenic mutations within the tissue. Inactivating mutations in key driver genes, such as TP53 and NOTCH1, reduce the proportion of differentiating cells allowing for the long‐term persistence of expanding mutant clones in the tissue. Here we show that through the expression of E6, HPV‐16 prevents the early fate commitment of human keratinocytes towards differentiation and confers a strong growth advantage to human keratinocytes. When E6 is expressed either alone or with E7, it promotes keratinocyte proliferation at high cell densities, through the combined inactivation of p53 and Notch1. In organotypic raft culture, the activity of E6 is restricted to the basal layer of the epithelium and is enhanced during the progression from productive to abortive or transforming HPV‐16 infection. Consistent with this, the expression of p53 and cleaved Notch1 becomes progressively more disrupted, and is associated with increased basal cell density and reduced commitment to differentiation. The expression of cleaved Notch1 is similarly disrupted also in HPV‐16‐positive cervical lesions, depending on neoplastic grade. When taken together, these data depict an important role of high‐risk E6 in promoting the persistence of infected keratinocytes in the basal and parabasal layers through the inactivation of gene products that are commonly mutated in non‐HPV‐associated neoplastic squamous epithelia. © 2017 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of Pathological Society of Great Britain and Ireland.

Host control of human papillomavirus infection and disease

Most human papillomaviruses cause inapparent infections, subtly affecting epithelial homeostasis, to ensure genome persistence in the epithelial basal layer. As with conspicuous papillomas, these self-limiting lesions shed viral particles to ensure population level maintenance and depend on a balance between viral gene expression, immune cell stimulation and immune surveillance for persistence. The complex immune evasion strategies, characteristic of high-risk HPV types, also allow the deregulated viral gene expression that underlies neoplasia. Neoplasia occurs at particular epithelial sites where vulnerable cells such as the reserve or cuboidal cells of the cervical transformation zone are found. Beta papillomavirus infection can also predispose an individual with immune deficiencies to the development of cancers. The host control of HPV infections thus involves local interactions between keratinocytes and the adaptive immune response. Effective immune detection and surveillance limits overt disease, leading to HPV persistence as productive microlesions or in a true latent state.

Protein Sam68 regulates the alternative splicing of survivin DEx3

Messenger RNA alternative splicing (AS) regulates the expression of a variety of genes involved in both physiological and pathological processes. AS of the anti-apoptotic and proliferation-associated survivin (BIRC5) gene generates six isoforms, which regulate key aspects of cancer initiation and progression. One of the isoforms is survivin DEx3, in which the exclusion of exon 3 generates a unique carboxyl terminus with specific anti-apoptotic functions. This isoform is highly expressed in advanced stages of breast and cervical tumors. Therefore, understanding the mechanisms that regulate survivin DEx3 mRNA AS is clearly important. To this end, we designed a minigene (M), and in combination with a series of deletions and site-directed mutations, we determined that the first 22 bp of exon 3 contain cis-acting elements that enhance the exclusion of exon 3 to generate the survivin DEx3 mRNA isoform. Furthermore, using pulldown assays, we discovered that Sam68 is a possible trans-acting factor that binds to this region and regulates exon 3 splicing. This result was corroborated using a cell line in which the Sam68 binding site in the survivin gene was mutated with the CRISPR/Cas system. This work provides the first clues regarding the regulation of survivin DEx3 mRNA splicing.

Alternative splicing: the pledge, the turn, and the prestige : The key role of alternative splicing in human biological systems

Alternative pre-mRNA splicing is a tightly controlled process conducted by the spliceosome, with the assistance of several regulators, resulting in the expression of different transcript isoforms from the same gene and increasing both transcriptome and proteome complexity. The differences between alternative isoforms may be subtle but enough to change the function or localization of the translated proteins. A fine control of the isoform balance is, therefore, needed throughout developmental stages and adult tissues or physiological conditions and it does not come as a surprise that several diseases are caused by its deregulation. In this review, we aim to bring the splicing machinery on stage and raise the curtain on its mechanisms and regulation throughout several systems and tissues of the human body, from neurodevelopment to the interactions with the human microbiome. We discuss, on one hand, the essential role of alternative splicing in assuring tissue function, diversity, and swiftness of response in these systems or tissues, and on the other hand, what goes wrong when its regulatory mechanisms fail. We also focus on the possibilities that splicing modulation therapies open for the future of personalized medicine, along with the leading techniques in this field. The final act of the spliceosome, however, is yet to be fully revealed, as more knowledge is needed regarding the complex regulatory network that coordinates alternative splicing and how its dysfunction leads to disease.

Splicing and Polyadenylation of Human Papillomavirus Type 16 mRNAs

The human papillomavirus type 16 (HPV16) life cycle can be divided into an early stage in which the HPV16 genomic DNA is replicated, and a late stage in which the HPV16 structural proteins are synthesized and virions are produced. A strong coupling between the viral life cycle and the differentiation state of the infected cell is highly characteristic of all HPVs. The switch from the HPV16 early gene expression program to the late requires a promoter switch, a polyadenylation signal switch and a shift in alternative splicing. A number of cis-acting RNA elements on the HPV16 mRNAs and cellular and viral factors interacting with these elements are involved in the control of HPV16 gene expression. This review summarizes our knowledge of HPV16 cis-acting RNA elements and cellular and viral trans-acting factors that regulate HPV16 gene expression at the level of splicing and polyadenylation.

Control of human papillomavirus gene expression by alternative splicing

•

Alternative splicing is a key cellular mechanism controlling HPV gene expression.

•

Many cellular SR proteins and hnRNPs have been identified that bind and control production of viral mRNAs.

•

HPV16 E2 protein controls expression of SR proteins and has splicing-related functions.

•

HPV16 infection through its regulatory effects on splicing factors may significantly alter cellular gene expression and cellular metabolism.

Human papillomaviruses possess circular double stranded DNA genomes of around 8 kb in size from which multiple mRNAs are synthesized during an infectious life cycle. Although at least three viral promoters are used to initiate transcription, viral mRNAs are largely the product of processing of pre-mRNAs by alternative splicing and polyadenylation. The HPV life cycle and viral gene expression are tightly linked to differentiation of the epithelium the virus infects: there is an orchestrated production of viral mRNAs and proteins. In this review we describe viral mRNA expression and the roles of the SR and hnRNP proteins that respectively positively and negatively regulate splicing. We discuss HPV regulation of splicing factors and detail the evidence that the papillomavirus E2 protein has splicing-related activities. We highlight the possibility that HPV-mediated control of splicing in differentiating epithelial cells may be necessary to accomplish the viral replication cycle.

Serine/Arginine-Rich Splicing Factor 3 and Heterogeneous Nuclear Ribonucleoprotein A1 Regulate Alternative RNA Splicing and Gene Expression of Human Papillomavirus 18 through Two Functionally Distinguishable cis Elements

Human papillomavirus 18 (HPV18) is the second most common oncogenic HPV type associated with cervical, anogenital, and oropharyngeal cancers. Like other oncogenic HPVs, HPV18 encodes two major (one early and one late) polycistronic pre-mRNAs that are regulated by alternative RNA splicing to produce a repertoire of viral transcripts for the expression of individual viral genes. However, RNA cis-regulatory elements and trans-acting factors contributing to HPV18 alternative RNA splicing remain unknown. In this study, an exonic splicing enhancer (ESE) in the nucleotide (nt) 3520 to 3550 region in the HPV18 genome was identified and characterized for promotion of HPV18 929^3434 splicing and E1^E4 production through interaction with SRSF3, a host oncogenic splicing factor differentially expressed in epithelial cells and keratinocytes. Introduction of point mutations in the SRSF3-binding site or knockdown of SRSF3 expression in cells reduces 929^3434 splicing and E1^E4 production but activates other, minor 929^3465 and 929^3506 splicing. Knockdown of SRSF3 expression also enhances the expression of E2 and L1 mRNAs. An exonic splicing silencer (ESS) in the HPV18 nt 612 to 639 region was identified as being inhibitory to the 233^416 splicing of HPV18 E6E7 pre-mRNAs via binding to hnRNP A1, a well-characterized, abundantly and ubiquitously expressed RNA-binding protein. Introduction of point mutations into the hnRNP A1-binding site or knockdown of hnRNP A1 expression promoted 233^416 splicing and reduced E6 expression. These data provide the first evidence that the alternative RNA splicing of HPV18 pre-mRNAs is subject to regulation by viral RNA cis elements and host trans-acting splicing factors.

IMPORTANCE

Expression of HPV18 genes is regulated by alternative RNA splicing of viral polycistronic pre-mRNAs to produce a repertoire of viral early and late transcripts. RNA cis elements and trans-acting factors contributing to HPV18 alternative RNA splicing have been discovered in this study for the first time. The identified ESS at the E7 open reading frame (ORF) prevents HPV18 233^416 splicing in the E6 ORF through interaction with a host splicing factor, hnRNP A1, and regulates E6 and E7 expression of the early E6E7 polycistronic pre-mRNA. The identified ESE at the E1^E4 ORF promotes HPV18 929^3434 splicing of both viral early and late pre-mRNAs and E1^E4 production through interaction with SRSF3. This study provides important observations on how alternative RNA splicing of HPV18 pre-mRNAs is subject to regulation by viral RNA cis elements and host splicing factors and offers potential therapeutic targets to overcome HPV-related cancer.

Human papillomavirus infection and induction of neoplasia: a matter of fitness

The aetiologic association between infection with certain human papillomavirus (HPV) types, high-grade squamous neoplasia, and cancer at different epithelial sites is well established. In this review we briefly discuss recent breakthroughs in the regulation of squamous epithelia in homeostasis and disease, and provide a view of how these discoveries modify our understanding of how HPV-induced neoplasia in squamous epithelia is triggered. Taken together, these observations highlight how HPVs have evolved the ability to inactivate the products of genes that are frequently mutated in non-HPV-associated pre-neoplasia and squamous cell carcinoma of sun-exposed skin, and introduce a Darwinian model of clonal evolution of HPV-infected cells. These concepts are considered against our current understanding of transformation zones where HPV-associated cancers occur more frequently, and other sites of non-productive (or abortive) HPV infection.

Nuclear phosphoproteomics analysis reveals that CDK1/2 are involved in EGF-regulated constitutive pre-mRNA splicing in MDA-MB-468 cells

The epidermal growth factor (EGF) receptor (EGFR) pathway is one of the most dysregulated and extensively investigated signaling pathways in human cancers and plays important roles in the regulation of nuclear functions through both cytoplasmic and nuclear EGFR pathways. However, the current understanding of the nuclear phosphorylation responses to activated EGFR pathways remains limited. In the present study, phosphoproteomics analysis revealed the increased phosphorylation of 90 nuclear proteins, primarily involved in RNA processing, pre-mRNA splicing and cell cycle regulation, upon EGF stimulation in MDA-MB-468 cells. Cellular splicing assays of the β-globin (HBB) minigene confirmed that EGF induced constitutive pre-mRNA splicing. Further analysis of phosphoproteomics data identified multiple CDK1/2 substrates in pre-mRNA splicing-related proteins, and both CDK1/2 inhibitors and CDK1/2 knockdowns reduced EGF-regulated pre-mRNA splicing. In conclusion, the results of the present study provide evidence that CDK1/2 participate in the regulation of constitutive pre-mRNA splicing by EGF stimulation in MDA-MB-468 cells.

SIGNIFICANCE

In this study, we successfully carried out a survey of nuclear phosphorylation changes in response to EGF stimulation. The results from the functional category analysis and pre-mRNA splicing assay strongly indicated that EGFR activation increased constitutive pre-mRNA splicing in MDA-MB-468 cells, revealing additional role of EGFR on regulation of mRNA maturation beyond alternative pre-mRNA splicing reported by previous studies. Furthermore, we found that CDK1/2 participated in constitutive pre-mRNA splicing regulation by EGF in MDA-MB-468 cells. Our study provides new knowledge for understanding the regulation of constitutive pre-mRNA splicing by EGF stimulation.