Regulation of human papillomavirus gene expression by splicing and polyadenylation

RNA Binding Proteins that Control Human Papillomavirus Gene Expression

The human papillomavirus (HPV) life cycle is strictly linked to the differentiation program of the infected mucosal epithelial cell. In the basal and lower levels of the epithelium, early genes coding for pro-mitotic proteins and viral replication factors are expressed, while terminal cell differentiation is required for activation of late gene expression and production of viral particles at the very top of the epithelium. Such productive infections are normally cleared within 18–24 months. In rare cases, the HPV infection is stuck in the early stage of the infection. Such infections may give rise to cervical lesions that can progress to cancer, primarily cancer of the uterine cervix. Since cancer progression is strictly linked to HPV gene expression, it is of interest to understand how HPV gene expression is regulated. Cis-acting HPV RNA elements and cellular RNA-binding proteins control HPV mRNA splicing and polyadenylation. These interactions are believed to play a particularly important role in the switch from early to late gene expression, thereby contributing to the pathogenesis of HPV. Indeed, it has been shown that the levels of various RNA binding proteins change in response to differentiation and in response to HPV induced cervical lesions and cancer. Here we have compiled published data on RNA binding proteins involved in the regulation of HPV gene expression.

Human Papillomavirus 16 Oncoprotein Expression Is Controlled by the Cellular Splicing Factor SRSF2 (SC35)

High-risk human papillomaviruses (HR-HPV) cause anogenital cancers, including cervical cancer, and head and neck cancers. Human papillomavirus 16 (HPV16) is the most prevalent HR-HPV. HPV oncogenesis is driven by two viral oncoproteins, E6 and E7, which are expressed through alternative splicing of a polycistronic RNA to yield four major splice isoforms (E6 full length, E6*I, E6*II, E6*X). The production of multiple mRNA isoforms from a single gene is controlled by serine/arginine-rich splicing factors (SRSFs), and HPV16 infection induces overexpression of a subset of these, SRSFs 1, 2, and 3. In this study, we examined whether these proteins could control HPV16 oncoprotein expression. Small interfering RNA (siRNA) depletion experiments revealed that SRSF1 did not affect oncoprotein RNA levels. While SRSF3 knockdown caused some reduction in E6E7 expression, depletion of SRSF2 resulted in a significant loss of E6E7 RNAs, resulting in reduced levels of the E6-regulated p53 proteins and E7 oncoprotein itself. SRSF2 contributed to the tumor phenotype of HPV16-positive cervical cancer cells, as its depletion resulted in decreased cell proliferation, reduced colony formation, and increased apoptosis. SRSF2 did not affect transcription from the P97 promoter that controls viral oncoprotein expression. Rather, RNA decay experiments showed that SRSF2 is required to maintain stability of E6E7 mRNAs. These data show that SRSF2 is a key regulator of HPV16 oncoprotein expression and cervical tumor maintenance.

IMPORTANCE

Expression of the HPV16 oncoproteins E7 and E6 drives HPV-associated tumor formation. Although increased transcription may yield increased levels of E6E7 mRNAs, it is known that the RNAs can have increased stability upon integration into the host genome. SR splicing factors (SRSFs) control splicing but can also control other events in the RNA life cycle, including RNA stability. Previously, we demonstrated increased levels of SRSFs 1, 2, and 3 during cervical tumor progression. Now we show that SRSF2 is required for expression of E6E7 mRNAs in cervical tumor but not nontumor cells and may act by inhibiting their decay. SRSF2 depletion in W12 tumor cells resulted in increased apoptosis, decreased proliferation, and decreased colony formation, suggesting that SRSF2 has oncogenic functions in cervical tumor progression. SRSF function can be targeted by known drugs that inhibit SRSF phosphorylation, suggesting a possible new avenue in abrogating HPV oncoprotein activity.

Acetylation of intragenic histones on HPV16 correlates with enhanced HPV16 gene expression

We report that many histone modifications are unevenly distributed over the HPV16 genome in cervical cancer cells as well as in HPV16-immortalized keratinocytes. For example, H3K36me3 and H3K9Ac that are common in highly expressed cellular genes and over exons, were more common in the early than in the late region of the HPV16 genome. In contrast, H3K9me3, H4K20me3, H2BK5me1 and H4K16Ac were more frequent in the HPV16 late region. Furthermore, a region encompassing the HPV16 early polyadenylation signal pAE displayed high levels of histone H3 acetylation. Histone deacetylase (HDAC) inhibitors caused a 2- to 8-fold induction of HPV16 early and late mRNAs in cervical cancer cells and in immortalized keratinocytes, while at the same time increasing the levels of acetylated histones in the cells and on the HPV16 genome specifically. We concluded that increased histone acetylation on the HPV16 genome correlates with increased HPV16 gene expression.

Heterogeneous Nuclear Ribonucleoprotein C Proteins Interact with the Human Papillomavirus Type 16 (HPV16) Early 3'-Untranslated Region and Alleviate Suppression of HPV16 Late L1 mRNA Splicing

In order to identify cellular factors that regulate human papillomavirus type 16 (HPV16) gene expression, cervical cancer cells permissive for HPV16 late gene expression were identified and characterized. These cells either contained a novel spliced variant of the L1 mRNAs that bypassed the suppressed HPV16 late, 5'-splice site SD3632; produced elevated levels of RNA-binding proteins SRSF1 (ASF/SF2), SRSF9 (SRp30c), and HuR that are known to regulate HPV16 late gene expression; or were shown by a gene expression array analysis to overexpress the RALYL RNA-binding protein of the heterogeneous nuclear ribonucleoprotein C (hnRNP C) family. Overexpression of RALYL or hnRNP C1 induced HPV16 late gene expression from HPV16 subgenomic plasmids and from episomal forms of the full-length HPV16 genome. This induction was dependent on the HPV16 early untranslated region. Binding of hnRNP C1 to the HPV16 early, untranslated region activated HPV16 late 5'-splice site SD3632 and resulted in production of HPV16 L1 mRNAs. Our results suggested that hnRNP C1 controls HPV16 late gene expression.

Association of human papillomavirus with Fanconi anemia promotes carcinogenesis in Fanconi anemia patients

Fanconi anemia (FA) is a rare recessive disorder associated with chromosomal fragility. FA patients are at very high risk of cancers, especially head and neck squamous cell carcinomas and squamous cell carcinomas caused by infection of human papillomaviruses (HPVs). By integrating into the host genome, HPV oncogenes E6 and E7 drive the genomic instability to promote DNA damage and gene mutations necessary for carcinogenesis in FA patients. Furthermore, E6 and E7 oncoproteins not only inhibit p53 and retinoblastoma but also impair the FANC/BRCA signaling pathway to prevent DNA damage repair and alter multiple signals including cell-cycle checkpoints, telomere function, cell proliferation, and interference of the host immune system leading to cancer development in FA patients. In this review, we summarize recent advances in unraveling the molecular mechanisms of FA susceptibility to HPV-induced cancers, which facilitate rational preventive and therapeutic strategies.

E6^E7, a novel splice isoform protein of human papillomavirus 16, stabilizes viral E6 and E7 oncoproteins via HSP90 and GRP78

Transcripts of human papillomavirus 16 (HPV16) E6 and E7 oncogenes undergo alternative RNA splicing to produce multiple splice isoforms. However, the importance of these splice isoforms is poorly understood. Here we report a critical role of E6^E7, a novel isoform containing the 41 N-terminal amino acid (aa) residues of E6 and the 38 C-terminal aa residues of E7, in the regulation of E6 and E7 stability. Through mass spectrometric analysis, we identified that HSP90 and GRP78, which are frequently upregulated in cervical cancer tissues, are two E6^E7-interacting proteins responsible for the stability and function of E6^E7, E6, and E7. Although GRP78 and HSP90 do not bind each other, GRP78, but not HSP90, interacts with E6 and E7. E6^E7 protein, in addition to self-binding, interacts with E6 and E7 in the presence of GRP78 and HSP90, leading to the stabilization of E6 and E7 by prolonging the half-life of each protein. Knocking down E6^E7 expression in HPV16-positive CaSki cells by a splice junction-specific small interfering RNA (siRNA) destabilizes E6 and E7 and prevents cell growth. The same is true for the cells with a GRP78 knockdown or in the presence of an HSP90 inhibitor. Moreover, mapping and alignment analyses for splicing elements in 36 alpha-HPVs (α-HPVs) suggest the possible expression of E6^E7 mostly by other oncogenic or possibly oncogenic α-HPVs (HPV18, -30, -31, -39, -42, -45, -56, -59, -70, and -73). HPV18 E6^E7 is detectable in HPV18-positive HeLa cells and HPV18-infected raft tissues. All together, our data indicate that viral E6^E7 and cellular GRP78 or HSP90 might be novel targets for cervical cancer therapy.

HPV16 is the most prevalent HPV genotype, being responsible for 60% of invasive cervical cancer cases worldwide. What makes HPV16 so potent in the development of cervical cancer remains a mystery. We discovered in this study that, besides producing two well-known oncoproteins, E6 and E7, seen in other high-risk HPVs, HPV16 produces E6^E7, a novel splice isoform of E6 and E7. E6^E7, in addition to self-interacting, binds cellular chaperone proteins, HSP90 and GRP78, and viral E6 and E7 to increase the steady-state levels and half-lives of viral oncoproteins, leading to cell proliferation. The splicing cis elements in the regulation of HPV16 E6^E7 production are highly conserved in 11 oncogenic or possibly oncogenic HPVs, and we confirmed the production of HPV18 E6^E7 in HPV18-infected cells. This study provides new insight into the mechanism of splicing, the interplay between different products of the polycistronic viral message, and the role of the host chaperones as they function.

HeLa nucleic acid contamination in the cancer genome atlas leads to the misidentification of human papillomavirus 18

We searched The Cancer Genome Atlas (TCGA) database for viruses by comparing non-human reads present in transcriptome sequencing (RNA-Seq) and whole-exome sequencing (WXS) data to viral sequence databases. Human papillomavirus 18 (HPV18) is an etiologic agent of cervical cancer, and as expected, we found robust expression of HPV18 genes in cervical cancer samples. In agreement with previous studies, we also found HPV18 transcripts in non-cervical cancer samples, including those from the colon, rectum, and normal kidney. However, in each of these cases, HPV18 gene expression was low, and single-nucleotide variants and positions of genomic alignments matched the integrated portion of HPV18 present in HeLa cells. Chimeric reads that match a known virus-cell junction of HPV18 integrated in HeLa cells were also present in some samples. We hypothesize that HPV18 sequences in these non-cervical samples are due to nucleic acid contamination from HeLa cells. This finding highlights the problems that contamination presents in computational virus detection pipelines.

IMPORTANCE

Viruses associated with cancer can be detected by searching tumor sequence databases. Several studies involving searches of the TCGA database have reported the presence of HPV18, a known cause of cervical cancer, in a small number of additional cancers, including those of the rectum, kidney, and colon. We have determined that the sequences related to HPV18 in non-cervical samples are due to nucleic acid contamination from HeLa cells. To our knowledge, this is the first report of the misidentification of viruses in next-generation sequencing data of tumors due to contamination with a cancer cell line. These results raise awareness of the difficulty of accurately identifying viruses in human sequence databases.

Pathogenesis of human papillomavirus-associated mucosal disease

Human papillomaviruses (HPVs) are a necessary cause of carcinoma of the cervix and other mucosal epithelia. Key events in high-risk HPV (HRHPV)-associated neoplastic progression include persistent infection, deregulated expression of virus early genes in basal epithelial cells and genomic instability causing secondary host genomic imbalances. There are multiple mechanisms by which deregulated virus early gene expression may be achieved. Integration of virus DNA into host chromosomes is observed in the majority of cervical squamous cell carcinomas (SCCs), although in ∼15% of cases the virus remains extrachromosomal (episomal). Interestingly, not all integration events provide a growth advantage to basal cervical epithelial cells or lead to increased levels of the virus oncogenes E6 and E7, when compared with episome-containing basal cells. The factors that provide a competitive advantage to some integrants, but not others, are complex and include virus and host contributions. Gene expression from integrated and episomal HRHPV is regulated through host epigenetic mechanisms affecting the virus long control region (LCR), which appear to be of functional importance. New approaches to treating HRHPV-associated mucosal neoplasia include knockout of integrated HRHPV DNA, depletion of virus transcripts and inhibition of virus early gene transcription through targeting or use of epigenetic modifiers. Copyright © 2014 Pathological Society of Great Britain and Ireland. Published by John Wiley & Sons, Ltd.

The cellular bromodomain protein Brd4 has multiple functions in E2-mediated papillomavirus transcription activation

The cellular bromodomain protein Brd4 functions in multiple processes of the papillomavirus life cycle, including viral replication, genome maintenance, and gene transcription through its interaction with the viral protein, E2. However, the mechanisms by which E2 and Brd4 activate viral transcription are still not completely understood. In this study, we show that recruitment of positive transcription elongation factor b (P-TEFb), a functional interaction partner of Brd4 in transcription activation, is important for E2’s transcription activation activity. Furthermore, chromatin immunoprecipitation (ChIP) analyses demonstrate that P-TEFb is recruited to the actual papillomavirus episomes. We also show that E2’s interaction with cellular chromatin through Brd4 correlates with its papillomavirus transcription activation function since JQ1(+), a bromodomain inhibitor that efficiently dissociates E2-Brd4 complexes from chromatin, potently reduces papillomavirus transcription. Our study identifies a specific function of Brd4 in papillomavirus gene transcription and highlights the potential use of bromodomain inhibitors as a method to disrupt the human papillomavirus (HPV) life cycle.

Biochemical and proteomic analysis of spliceosome factors interacting with intron-1 of human papillomavirus type-16

The human papillomavirus type 16 (HPV-16) E6/E7 spliced transcripts are heterogeneously expressed in cervical carcinoma. The heterogeneity of the E6/E7 splicing profile might be in part due to the intrinsic variation of splicing factors in tumor cells. However, the splicing factors that bind the E6/E7 intron 1 (In-1) have not been defined. Therefore, we aimed to identify these factors; we used HeLa nuclear extracts (NE) for in vitro spliceosome assembly. The proteins were allowed to bind to an RNA/DNA hybrid formed by the In-1 transcript and a 5'-biotinylated DNA oligonucleotide complementary to the upstream exon sequence, which prevented interference in protein binding to the intron. The hybrid probes bound with the nuclear proteins were coupled to streptavidin magnetic beads for chromatography affinity purification. Proteins were eluted and identified by mass spectrometry (MS). Approximately 170 proteins were identified by MS, 80% of which were RNA binding proteins, including canonical spliceosome core components, helicases and regulatory splicing factors. The canonical factors were identified as components of the spliceosomal B-complex. Although 35-40 of the identified factors were cognate splicing factors or helicases, they have not been previously detected in spliceosome complexes that were assembled using in vivo or in vitro models.

Impact of inhibitors and L2 antibodies upon the infectivity of diverse alpha and beta human papillomavirus types

The licensed human papillomavirus (HPV) vaccines elicit type-restricted immunity but do not target cutaneous HPV types of the beta genus that are associated with non-melanoma skin cancer in immune-compromised patients, and it is unclear if these diverse types share a common mechanism of infection. Residues 11-88 of minor capsid protein L2 contain cross-protective epitopes, and vaccination with concatamers of this region derived from as many as eight alpha HPV (L2 α11-88x8) is being developed as an alternative prophylactic vaccine with potentially broader efficacy. There is also interest in developing broadly protective topical microbicides, such as carrageenan or heparin that block HPV receptor interactions, or small molecule inhibitors of infection. Here we have examined several inhibitors of HPV infection and antisera to L2 α11-88x8 for their breadth of activity against infection by 34 HPV types from within both the alpha and beta families using pseudovirions (PsV) carrying a luciferase reporter as surrogates for native virus. We observed that both heparin and carrageenan prevented infection by mucosatropic HPV types, but surprisingly PsV of several epidermotropic alpha4 and beta HPV types exhibited increased infectivity especially at low inhibitor concentrations. Furin and γ-secretase inhibitors and L2 α11-88x8 antiserum blocked infection by all HPV PsV types tested. These findings suggest that the distinct tropism of mucosal and cutaneous HPV may reflect distinct cell surface receptor interactions, but a common uptake mechanism dependent upon furin and γ-secretase proteolytic activities. Carrageenan, which is being tested as a vaginal microbicide, broadly inhibited infection by the high-risk mucosatropic HPV PsV, but not most skin tropic alpha and beta HPV. Vaccination with an L2 multimer derived exclusively from alpha papillomavirus sequences induced antibodies that broadly neutralized PsV of all 34 HPVs from within both the alpha and beta families, suggesting each displays conserved L2 neutralizing epitopes.

Ménage à trois: an evolutionary interplay between human papillomavirus, a tumor, and a woman

Cervical cancer is the third most common cancer in women with human papillomavirus (HPV) being a key etiologic factor of this devastating disease. In this article, we describe modern advances in the genomics and transcriptomics of cervical cancer that led to uncovering the key gene drivers. We also introduce, herein, a model of cervical carcinogenesis that explains how the interplay between virus, tumor, and woman results in the selection of clones that simultaneously harbor genomic amplifications for genes that drive cell cycle, antiviral response, and inhibit cell differentiation. The new model may help researchers understand the controversies in antiviral therapy and immunogenetics of this cancer and may provide a basis for future research directions in early diagnostics and personalization of therapy.

The papillomavirus E2 proteins

The papillomavirus E2 proteins are pivotal to the viral life cycle and have well characterized functions in transcriptional regulation, initiation of DNA replication and partitioning the viral genome. The E2 proteins also function in vegetative DNA replication, post-transcriptional processes and possibly packaging. This review describes structural and functional aspects of the E2 proteins and their binding sites on the viral genome. It is intended to be a reference guide to this viral protein.

Suppression of HPV-16 late L1 5'-splice site SD3632 by binding of hnRNP D proteins and hnRNP A2/B1 to upstream AUAGUA RNA motifs

Human papillomavirus type 16 (HPV-16) 5′-splice site SD3632 is used exclusively to produce late L1 mRNAs. We identified a 34-nt splicing inhibitory element located immediately upstream of HPV-16 late 5′-splice site SD3632. Two AUAGUA motifs located in these 34 nt inhibited SD3632. Two nucleotide substitutions in each of the HPV-16 specific AUAGUA motifs alleviated splicing inhibition and induced late L1 mRNA production from episomal forms of the HPV-16 genome in primary human keratinocytes. The AUAGUA motifs bind specifically not only to the heterogeneous nuclear RNP (hnRNP) D family of RNA-binding proteins including hnRNP D/AUF, hnRNP DL and hnRNP AB but also to hnRNP A2/B1. Knock-down of these proteins induced HPV-16 late L1 mRNA expression, and overexpression of hnRNP A2/B1, hnRNP AB, hnRNP DL and the two hnRNP D isoforms hnRNP D37 and hnRNP D40 further suppressed L1 mRNA expression. This inhibition may allow HPV-16 to hide from the immune system and establish long-term persistent infections with enhanced risk at progressing to cancer. There is an inverse correlation between expression of hnRNP D proteins and hnRNP A2/B1 and HPV-16 L1 production in the cervical epithelium, as well as in cervical cancer, supporting the conclusion that hnRNP D proteins and A2/B1 inhibit HPV-16 L1 mRNA production.

Current understanding of the role of the Brd4 protein in the papillomavirus lifecycle

The Brd4 protein is an epigenetic reader that is central to regulation of cellular transcription and mitotic bookmarking. The transcription and replication proteins of many viruses interact with Brd4. We describe the multiple roles of Brd4 in the papillomavirus lifecycle.