The transcription map of human papillomavirus type 18 during genome replication in U2OS cells

PRMT1 Modulates Alternative Splicing to Enhance HPV18 mRNA Stability and Promote the Establishment of Infection

Only persistent HPV infections lead to the development of cancer. Thus, understanding the virus-host interplay that influences the establishment of viral infection has important implications for HPV biology and human cancers. The ability of papillomaviruses to establish in cells requires the strict temporal regulation of viral gene expression in sync with cellular differentiation. This control primarily happens at the level of RNA splicing and polyadenylation. However, the details of how this spatio-temporal regulation is achieved still need to be fully understood. Until recently, it has been challenging to study the early events of the HPV lifecycle following infection. We used a single-cell genomics approach to identify cellular factors involved in viral infection and establishment. We identify protein arginine N-methyltransferase 1 (PRMT1) as an important factor in viral infection of primary human cervical cells. PRMT1 is the main cellular enzyme responsible for asymmetric dimethylation of cellular proteins. PRMT1 is an enzyme responsible for catalyzing the methylation of arginine residues on various proteins, which influences processes such as RNA processing, transcriptional regulation, and signal transduction. In this study, we show that HPV18 infection leads to increased PRMT1 levels across the viral lifecycle. PRMT1 is critical for the establishment of a persistent infection in primary cells. Mechanistically, PRMT1 inhibition leads to a highly dysregulated viral splicing pattern. Specifically, reduced PRMT1 activity leads to intron retention and a change in the E6 and E7 expression ratio. In the absence of PRMT1, viral transcripts are destabilized and subject to degradation via the nonsense-mediated decay (NMD) pathway. These findings highlight PRMT1 as a critical regulator of the HPV18 lifecycle, particularly in RNA processing, and position it as a potential therapeutic target for persistent HPV18 infections.

Chimera RNA transcribed from integrated HPV18 genome with adjacent host genomic region promotes oncogenic gene expression through condensate formation

Most cervical cancers are caused by human papillomavirus (HPV) infection. In HeLa cells, the HPV18 viral genome is integrated at chromosome 8q24.21 and activates transcription of the proto-oncogene c-Myc. However, the mechanism of how the integrated HPV genome and its transcribed RNAs exhibit transcription activation function has not been fully elucidated. In this study, we found that HPV18 transcripts contain an enhancer RNA-like function to activate proximal genes including CCAT1-5L and c-Myc. We showed that the human genome-integrated HPV18 genes are activated by transcription coregulators including BRD4 and Mediator. The transcribed HPV18 RNAs form a liquid-like condensate at chromosome 8q24.21 locus, which in turn accumulates RNA polymerase II. Moreover, we focused on a relatively uncharacterized transcript from the upstream region of CCAT1, named URC. The URC RNA is transcribed as a chimera RNA with HPV18 and is composed of the 3'-untranslated region of the HPV18 transcript. We experimentally showed that the URC contributes to stabilization of HPV18 RNAs by supplying a polyadenylation site for the HPV18 transcript. Our findings suggest that integrated HPV18 at 8q24.21 locus produces HPV18-URC chimera RNA and promotes tumorigenesis through RNA-based condensate formation.

The POU-HD TFs impede the replication efficiency of several human papillomavirus genomes

Human papillomavirus (HPV) is a double-stranded DNA virus that infects cutaneous and mucosal epithelial cells. HPV replication initiates at the origin (ori), located within a noncoding region near the major early promoter. Only two viral proteins, E1 and E2, are essential for replication, with the host cell contributing other necessary factors. However, the role of host cell proteins in regulating HPV replication remains poorly understood. While several binding sites for cellular transcription factors (TFs), such as POU-HD proteins, have been mapped in the regulatory region, their functional importance is unclear. Some POU-HD TFs have been shown to influence replication in a system where E1 and E2 are provided exogenously. In this study, we investigated the impact of several POU-HD TFs on the replication of the HPV5, HPV11, and HPV18 genomes in U2OS cells and human primary keratinocytes. We demonstrated that OCT1, OCT6, BRN5A, and SKN1A are expressed in HPV host cells and that their overexpression inhibits HPV genome replication, whereas knocking down OCT1 had a positive effect. Using the replication-deficient HPV18-E1- genome, we demonstrated that OCT1-mediated inhibition of HPV replication involves modulation of HPV early promoters controlling E1 and E2 expression. Moreover, using Oct6 mutants deficient either in DNA binding or transcriptional regulation, we showed that the inhibition of HPV18 replication is solely dependent on Oct6's DNA binding activity. Our study highlights the complex regulatory roles of POU-HD factors in the HPV replication.

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.

Functions of Papillomavirus E8^E2 Proteins in Tissue Culture and In Vivo

Papillomaviruses (PV) replicate in undifferentiated keratinocytes at low levels and to high levels in differentiated cells. The restricted replication in undifferentiated cells is mainly due to the expression of the conserved viral E8^E2 repressor protein, a fusion protein consisting of E8 and the hinge, DNA-binding, and dimerization domain of E2. E8^E2 binds to viral genomes and represses viral transcription and genome replication by recruiting cellular NCoR/SMRT-HDAC3 corepressor complexes. Tissue culture experiments have revealed that E8^E2 modulates long-term maintenance of extrachromosomal genomes, productive replication, and immortalization properties in a virus type-dependent manner. Furthermore, in vivo experiments have indicated that Mus musculus PV1 E8^E2 is required for tumor formation in immune-deficient mice. In summary, E8^E2 is a crucial inhibitor whose levels might determine the outcome of PV infections.

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.

Transcriptome analysis of cervical cancer exosomes and detection of HPVE6*I transcripts in exosomal RNA

Background

Exosomes play a key role in cell-to-cell communication and are integral component of the tumor microenvironment. Recent observations suggest transfer of RNA through tumor-derived exosomes that can potentially translate into regulatory proteins in the recipient cells. Role of cervical cancer-derived exosomes and their transcript cargo is poorly understood.

Materials and methods

The total RNA of exosomes from HPV-positive (SiHa and HeLa) and HPV-negative (C33a) cervical cancer cell lines were extracted and the transcripts were estimated using Illumina HiSeq X. Further, validation of HPV transcripts were performed using RT-PCR.

Results

3099 transcripts were found to be differentially-exported in HPV-positive vs. HPV-negative exosomes (p value <0.05). Analysis of top 10 GO terms and KEGG pathways showed enrichment of transcripts belonging to axon guidance and tumor innervation in HPV-positive exosomes. Among top 20 overexpressed transcripts, EVC2, LUZP1 and ANKS1B were the most notable due to their involvement in Hh signaling, cellular migration and invasion, respectively. Further, low levels of HPV-specific reads were detected. RT-PCR validation revealed presence of E6*I splice variant of HPV18 in exosomal RNA of HeLa cells. The E6*I transcripts were consistently retained in exosomes obtained from HeLa cells undergoing 5-FU and cisplatin-induced oxidative stress.

Conclusion

Our data suggests the enrichment of poly-A RNA transcripts in the exosomal cargo of cervical cancer cells, which includes pro-tumorigenic cellular RNA and viral transcripts such as HPV E6, which may have clinical utility as potential exosomal biomarkers of cervical cancer.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12885-022-09262-4.

The chromatin insulator CTCF regulates HPV18 transcript splicing and differentiation-dependent late gene expression

The ubiquitous host protein, CCCTC-binding factor (CTCF), is an essential regulator of cellular transcription and functions to maintain epigenetic boundaries, stabilise chromatin loops and regulate splicing of alternative exons. We have previously demonstrated that CTCF binds to the E2 open reading frame (ORF) of human papillomavirus (HPV) 18 and functions to repress viral oncogene expression in undifferentiated keratinocytes by co-ordinating an epigenetically repressed chromatin loop within HPV episomes. Keratinocyte differentiation disrupts CTCF-dependent chromatin looping of HPV18 episomes promoting induction of enhanced viral oncogene expression. To further characterise CTCF function in HPV transcription control we utilised direct, long-read Nanopore RNA-sequencing which provides information on the structure and abundance of full-length transcripts. Nanopore analysis of primary human keratinocytes containing HPV18 episomes before and after synchronous differentiation allowed quantification of viral transcript species, including the identification of low abundance novel transcripts. Comparison of transcripts produced in wild type HPV18 genome-containing cells to those identified in CTCF-binding deficient genome-containing cells identifies CTCF as a key regulator of differentiation-dependent late promoter activation, required for efficient E1^E4 and L1 protein expression. Furthermore, our data show that CTCF binding at the E2 ORF promotes usage of the downstream weak splice donor (SD) sites SD3165 and SD3284, to the dominant E4 splice acceptor site at nucleotide 3434. These findings demonstrate that in the HPV life cycle both early and late virus transcription programmes are facilitated by recruitment of CTCF to the E2 ORF.

Oncogenic human papillomavirus (HPV) infection is the cause of a subset of epithelial cancers of the uterine cervix, other anogenital areas and the oropharynx. HPV infection is established in the basal cells of epithelia where a restricted programme of viral gene expression is required for replication and maintenance of the viral episome. Completion of the HPV life cycle is dependent on the maturation (differentiation) of infected cells which induces enhanced viral gene expression and induction of capsid production. We previously reported that the host cell transcriptional regulator, CTCF, is hijacked by HPV to control viral gene expression. In this study, we use long-read mRNA sequencing to quantitatively map the variety and abundance of HPV transcripts produced in early and late stages of the HPV life cycle and to dissect the function of CTCF in controlling HPV gene expression and transcript processing.

Uncovering the Role of the E1 Protein in Different Stages of Human Papillomavirus 18 Genome Replication

The life cycle of human papillomaviruses (HPVs) comprises three distinct phases of DNA replication: initial amplification, maintenance of the genome copy number at a constant level, and vegetative amplification. The viral helicase E1 is one of the factors required for the initiation of HPV genome replication. However, the functions of the E1 protein during other phases of the viral life cycle are largely uncharacterized. Here, we studied the role of the HPV18 E1 helicase in three phases of viral genome replication by downregulating E1 expression using RNA interference or inducing degradation of the E1 protein via inhibition of casein kinase 2α expression or catalytic activity. We generated a novel modified HPV18 genome expressing Nanoluc and tagged E1 and E2 proteins and created several stable HPV18-positive cell lines. We showed that, in contrast to initial amplification of the HPV18 genome, other phases of viral genome replication involve also an E1-independent mechanism. We characterize two distinct populations of HPV18 replicons existing during the maintenance and vegetative amplification phases. We show that a subset of these replicons, including viral genome monomers, replicate in an E1-dependent manner, while some oligomeric forms of the HPV18 genome replicate independently of E1 function.IMPORTANCE Human papillomavirus (HPV) infections pose serious medical problem. To date, there are no HPV-specific antivirals available due to poor understanding of the molecular mechanisms of virus infection cycle. The infection cycle of HPV involves initial amplification of the viral genomes and maintenance of the viral genomes with a constant copy number, followed by another round of viral genome amplification and new viral particle formation. The viral protein E1 is critical for the initial amplification of the viral genome. However, E1 involvement in other phases of the viral life cycle has remained controversial. In the present study, we show that at least two different replication modes of the HPV18 genome are undertaken simultaneously during the maintenance and vegetative amplification phases, i.e., replication of the majority of the HPV18 genome proceeds under the control of the host cell replication machinery without E1 function, whereas a minority of the genome replicates in an E1-dependent manner.

Differential phosphorylation determines the repressor and activator potencies of GLI1 proteins and their efficiency in modulating the HPV life cycle

The Sonic Hedgehog (Shh) signalling pathway plays multiple roles during embryonic development and under pathological conditions. Although the core components of the Shh pathway are conserved, the regulation of signal transduction varies significantly among species and cell types. Protein kinases Ulk3 and Pka are involved in the Shh pathway as modulators of the activities of Gli transcription factors, which are the nuclear mediators of the signal. Here, we investigate the regulation and activities of two GLI1 isoforms, full-length GLI1 (GLI1FL) and GLI1ΔN. The latter protein lacks the first 128 amino acids including the conserved phosphorylation cluster and the binding motif for SUFU, the key regulator of GLI activity. Both GLI1 isoforms are co-expressed in all human cell lines analysed and possess similar DNA binding activity. ULK3 potentiates the transcriptional activity of both GLI1 proteins, whereas PKA inhibits the activity of GLI1ΔN, but not GLI1FL. In addition to its well-established role as a transcriptional activator, GLI1FL acts as a repressor by inhibiting transcription from the early promoters of human papillomavirus type 18 (HPV18). Additionally, compared to GLI1ΔN, GLI1FL is a more potent suppressor of replication of several HPV types. Altogether, our data show that the N-terminal part of GLI1FL is crucial for the realization of its full potential as a transcriptional regulator.

Activity of CK2α protein kinase is required for efficient replication of some HPV types

Inhibition of human papillomavirus (HPV) replication is a promising therapeutic approach for intervening with HPV-related pathologies. Primary targets for interference are two viral proteins, E1 and E2, which are required for HPV replication. Both E1 and E2 are phosphoproteins; thus, the protein kinases that phosphorylate them might represent secondary targets to achieve inhibition of HPV replication. In the present study, we show that CX4945, an ATP-competitive small molecule inhibitor of casein kinase 2 (CK2) catalytic activity, suppresses replication of different HPV types, including novel HPV5NLuc, HPV11NLuc and HPV18NLuc marker genomes, but enhances the replication of HPV16 and HPV31. We further corroborate our findings using short interfering RNA (siRNA)-mediated knockdown of CK2 α and α' subunits in U2OS and CIN612 cells; we show that while both subunits are expressed in these cell lines, CK2α is required for HPV replication, but CK2α' is not. Furthermore, we demonstrate that CK2α acts in a kinase activity-dependent manner and regulates the stability and nuclear retention of endogenous E1 proteins of HPV11 and HPV18. This unique feature of CK2α makes it an attractive target for developing antiviral agents.

Characterization of an HPV33 natural variant with enhanced transcriptional activity suggests a role for C/EBPβ in the regulation of the viral early promoter

The Long Control Region (LCR) of the human papillomavirus (HPV) genome encompasses the early promoter (EP) that drives expression of the viral oncogenes in infected cells and HPV-associated cancers. Here, we report on a natural variant of HPV33 that displays higher EP activity than the prototype in transfected C33A and HeLa cervical carcinoma cells, and in the osteosarcoma U2OS cell line which supports replication of HPV episomes. This increased promoter activity was ascribed to a single nucleotide variation in the LCR, T7791C, in a putative binding site for the transcription factor C/EBPβ. T7791C abrogated binding of recombinant C/EBPβ to this site in vitro and stimulated the EP in vivo, suggesting that it abrogates a negatively-acting regulatory element. A second C/EBPβ binding site was identified in vitro that activated the EP in vivo and whose function and location in the epithelial-specific enhancer is shown to be conserved in the highly prevalent HPV18. These results suggest that C/EBPβ is both an activator and a repressor of the HPV33 EP, acting via two distinct binding sites. Prediction of C/EBPβ sites in the LCR of 186 HPV types suggests that C/EBPβ regulation of the EP is common among high‐risk viruses from the α genus.

The molecular biology and HPV drug responsiveness of cynomolgus macaque papillomaviruses support their use in the development of a relevant in vivo model for antiviral drug testing

Due to the extreme tissue and species restriction of the papillomaviruses (PVs), there is a great need for animal models that accurately mimic PV infection in humans for testing therapeutic strategies against human papillomaviruses (HPVs). In this study, we present data that demonstrate that in terms of gene expression during initial viral DNA amplification, Macaca fascicularis PV (MfPV) types 5 and 8 appear to be similar to mucosal oncogenic HPVs, while MfPV1 (isolated from skin) resembles most high-risk cutaneous beta HPVs (HPV5). Similarities were also observed in replication properties during the initial amplification phase of the MfPV genomes. We demonstrate that high-risk mucosal HPV-specific inhibitors target the transient replication of the MfPV8 genomes, which indicates that similar pathways are used by the high-risk HPVs and MfPVs during their genome replication. Taking all into account, we propose that Macaca fascicularis may serve as a highly relevant model for preclinical tests designed to evaluate therapeutic strategies against HPV-associated lesions.

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.

Identification of several high-risk HPV inhibitors and drug targets with a novel high-throughput screening assay

Human papillomaviruses (HPVs) are oncogenic viruses that cause numerous different cancers as well as benign lesions in the epithelia. To date, there is no effective cure for an ongoing HPV infection. Here, we describe the generation process of a platform for the development of anti-HPV drugs. This system consists of engineered full-length HPV genomes that express reporter genes for evaluation of the viral copy number in all three HPV replication stages. We demonstrate the usefulness of this system by conducting high-throughput screens to identify novel high-risk HPV-specific inhibitors. At least five of the inhibitors block the function of Tdp1 and PARP1, which have been identified as essential cellular proteins for HPV replication and promising candidates for the development of antivirals against HPV and possibly against HPV-related cancers.

Human papillomaviruses are causative agents of many different cancers; they are most commonly associated with cervical cancer which leads to about quarter of a million deaths each year. Regardless of extensive studies for decades there is no specific cure against HPV infection. During this research, we have engineered modified HPV marker genomes that express Renilla luciferase reporter gene which expression level correlates directly with viral genome copy number. We have used such modified HPV genome in high-throughput screening of NCI Diversity Set IV chemical library and have identified a number of novel high-risk HPV-specific chemical compounds and drug targets. Such Renilla-expressing marker genomes could be used in various cell systems suitable for HPV replication studies to conduct high-throughput screens and quantify viral genome copy number quickly and effectively.

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.

Control of viral replication and transcription by the papillomavirus E8^E2 protein

Human papillomaviruses have adjusted their replication levels to the differentiation state of the infected keratinocyte. PV genomes replicate in undifferentiated cells at low levels and to high levels in differentiated cells. Genome replication requires the viral E1 helicase and the viral E2 transcription/replication activator. The limited replication in undifferentiated cells is predominantly due to the expression of the highly conserved E8^E2 viral repressor protein, which is a fusion between E8 and the C-terminal half of the E2 protein. E8^E2 is a sequence-specific DNA binding protein that inhibits viral gene expression and viral genome replication. The E8 domain is required for repression activities, which are mainly due to the interaction with cellular NCoR/SMRT corepressor complexes. In the case of HPV16, the most carcinogenic HPV type, E8^E2 not only limits genome replication in undifferentiated cells but also productive replication in differentiated epithelium. E8^E2 is expressed from a separate promoter that is controlled by unknown cellular factors and the viral transcription and replication regulators E1, E2 and E8^E2. In summary, E8^E2 is an important negative regulator whose levels may be critical for the outcome of HPV infections.

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.

Conservation of the E8 CDS of the E8^E2 protein among mammalian papillomaviruses

Papillomaviridae are small dsDNA viruses with a limited coding capacity. To fulfill all of the functional requirements for propagation and spreading, papillomaviruses use double coding and alternative protein isoforms. E8 ^ E2 is an alternative E2 protein isoform that is generated by fusing the short E8 CDS that completely overlaps E1 to the 'hinge' and the DNA-binding region of the E2 protein via alternative transcription/splicing. The papillomaviruses in which E8 ^ E2 mRNA sequences have been described exhibit a sparse phylogenomic distribution. Thus, it is not clear whether E8 ^ E2 is an ancestral protein that has not been described for other papillomavirus types or whether it randomly appears because of the conservation of the E1 protein and occurs only coincidentally. We searched for potential E8 coding sequences in a non-redundant set of papillomaviruses and applied SynPlot2 and an in-house-developed algorithm (cRegions) to determine the most plausible of the above-mentioned scenarios. Beginning with nine experimentally described E8 ^ E2 mRNAs, we predicted the potential E8 CDSs for more than 300 mammalian papillomavirus genomes. According to our analysis, E8 ^ E2 is not a result of E1 coding and represents a protein in its own right, and it most likely has an ancestral origin that precedes the divergence of major mammalian papillomavirus genera.

Full-length single-cell RNA-seq applied to a viral human cancer: applications to HPV expression and splicing analysis in HeLa S3 cells

Background

Viral infection causes multiple forms of human cancer, and HPV infection is the primary factor in cervical carcinomas. Recent single-cell RNA-seq studies highlight the tumor heterogeneity present in most cancers, but virally induced tumors have not been studied. HeLa is a well characterized HPV+ cervical cancer cell line.

Result

We developed a new high throughput platform to prepare single-cell RNA on a nanoliter scale based on a customized microwell chip. Using this method, we successfully amplified full-length transcripts of 669 single HeLa S3 cells and 40 of them were randomly selected to perform single-cell RNA sequencing. Based on these data, we obtained a comprehensive understanding of the heterogeneity of HeLa S3 cells in gene expression, alternative splicing and fusions. Furthermore, we identified a high diversity of HPV-18 expression and splicing at the single-cell level. By co-expression analysis we identified 283 E6, E7 co-regulated genes, including CDC25, PCNA, PLK4, BUB1B and IRF1 known to interact with HPV viral proteins.

Conclusion

Our results reveal the heterogeneity of a virus-infected cell line. It not only provides a transcriptome characterization of HeLa S3 cells at the single cell level, but is a demonstration of the power of single cell RNA-seq analysis of virally infected cells and cancers.

Electronic supplementary material

The online version of this article (doi:10.1186/s13742-015-0091-4) contains supplementary material, which is available to authorized users.

Initial amplification of the HPV18 genome proceeds via two distinct replication mechanisms

Determining the mechanism of HPV18 replication is paramount for identifying possible drug targets against HPV infection. We used two-dimensional and three-dimensional gel electrophoresis techniques to identify replication intermediates arising during the initial amplification of HPV18 episomal genomes. We determined that the first rounds of HPV18 replication proceed via bidirectional theta structures; however, a notable accumulation of almost fully replicated HPV18 genomes indicates difficulties with the completion of theta replication. We also observed intermediates that were created by a second replication mechanism during the initial amplification of HPV18 genomes. The second replication mechanism does not utilize specific initiation or termination sequences and proceeds via a unidirectional replication fork. We suggest a significant role for the second replication mechanism during the initial replication of the HPV18 genome and propose that the second replication mechanism is recombination-dependent replication.

Human Papillomavirus Type 18 cis-Elements Crucial for Segregation and Latency

Stable maintenance replication is characteristic of the latency phase of HPV infection, during which the viral genomes are actively maintained as extrachromosomal genetic elements in infected proliferating basal keratinocytes. Active replication in the S-phase and segregation of the genome into daughter cells in mitosis are required for stable maintenance replication. Most of our knowledge about papillomavirus genome segregation has come from studies of bovine papillomavirus type 1 (BPV-1), which have demonstrated that the E2 protein cooperates with cellular trans-factors and that E2 binding sites act as cis-regulatory elements in the viral genome that are essential for the segregation process. However, the genomic organization of the regulatory region in HPVs, and the properties of the viral proteins are different from those of their BPV-1 counterparts. We have designed a segregation assay for HPV-18 and used it to demonstrate that the E2 protein performs segregation in combination with at least two E2 binding sites. The cooperative binding of the E2 protein to two E2 binding sites is a major determinant of HPV-18 genome segregation, as demonstrated by the change in spacing between adjacent binding sites #1 and #2 in the HPV-18 Upstream Regulatory Region (URR). Duplication or triplication of the natural 4 bp 5’-CGGG-3’ spacer between the E2 binding sites increased the cooperative binding of the E2 molecules as well as E2-dependent segregation. Removal of any spacing between these sites eliminated cooperative binding of the E2 protein and disabled segregation of the URR and HPV-18 genome. Transfer of these configurations of the E2 binding sites into viral genomes confirmed the role of the E2 protein and binding sites #1 and #2 in the segregation process. Additional analysis demonstrated that these sites also play an important role in the transcriptional regulation of viral gene expression from different HPV-18 promoters.

DAXX modulates human papillomavirus early gene expression and genome replication in U2OS cells

Background

The human papillomavirus (HPV) genomes can replicate, and are maintained as autonomously replicating extrachromosomal plasmids in human U2OS cells. Previous studies have shown that HPV genomes are transcriptionally active in U2OS cells and can express the viral early proteins required for initiation and establishment of HPV replication. In the present work, we have examined the involvement of cellular DAXX protein in HPV replication in U2OS cells.

Methods

We have used indirect immunofluorescence and FISH analysis in order to study HPV replication compartments in U2OS cells. In addition, we have used siRNA knock-down for examining the effect of the DAXX protein on HPV replication and transcription in U2OS cells.

Results

We show that a portion of HPV replication foci are partially co-localized with components of ND10, cellular DAXX and PML proteins. In addition, we demonstrate that the knock-down of the cellular DAXX protein modulates the HPV genome replication and transcription in U2OS cells – papillomavirus replication is reduced in the absence of this component of ND10.

Conclusions

The DAXX protein modulates the early gene expression and the transient replication of HPV genomes in U2OS cells.

The Cell Cycle Timing of Human Papillomavirus DNA Replication

Viruses manipulate the cell cycle of the host cell to optimize conditions for more efficient viral genome replication. One strategy utilized by DNA viruses is to replicate their genomes non-concurrently with the host genome; in this case, the viral genome is amplified outside S phase. This phenomenon has also been described for human papillomavirus (HPV) vegetative genome replication, which occurs in G2-arrested cells; however, the precise timing of viral DNA replication during initial and stable replication phases has not been studied. We developed a new method to quantitate newly synthesized DNA levels and used this method in combination with cell cycle synchronization to show that viral DNA replication is initiated during S phase and is extended to G2 during initial amplification but follows the replication pattern of cellular DNA during S phase in the stable maintenance phase. E1 and E2 protein overexpression changes the replication time from S only to both the S and G2 phases in cells that stably maintain viral episomes. These data demonstrate that the active synthesis and replication of the HPV genome are extended into the G2 phase to amplify its copy number and the duration of HPV genome replication is controlled by the level of the viral replication proteins E1 and E2. Using the G2 phase for genome amplification may be an important adaptation that allows exploitation of changing cellular conditions during cell cycle progression. We also describe a new method to quantify newly synthesized viral DNA levels and discuss its benefits for HPV research.

Characterization of the Human Papillomavirus 16 E8 Promoter

Persistent infections with certain human papillomaviruses (HPV) such as HPV16 are a necessary risk factor for the development of anogenital and oropharyngeal cancers. HPV16 genomes replicate as low-copy-number plasmids in the nucleus of undifferentiated keratinocytes, which requires the viral E1 and E2 replication proteins. The HPV16 E8^E2C (or E8^E2) protein limits genome replication by repressing both viral transcription and the E1/E2-dependent DNA replication. How E8^E2C expression is regulated is not understood. Previous transcript analyses indicated that the spliced E8^E2C RNA is initiated at a promoter located in the E1 region upstream of the E8 gene. Deletion and mutational analyses of the E8 promoter region identify two conserved elements that are required for basal promoter activity in HPV-negative keratinocytes. In contrast, the transcriptional enhancer in the upstream regulatory region of HPV16 does not modulate basal E8 promoter activity. Cotransfection studies indicate that E8^E2C inhibits, whereas E2 weakly activates, the E8 promoter. Interestingly, the cotransfection of E1 and E2 induces the E8 promoter much more strongly than the major early promoter, and this is partially dependent upon binding of E2 to Brd4. Mutation of E8 promoter elements in the context of HPV16 genomes results in an increased genome copy number and elevated levels of viral early and late transcripts. In summary, the promoter responsible for the expression of E8^E2C is both positively and negatively regulated by viral and cellular factors, and this regulatory circuit may be crucial to maintain a low but constant copy number of HPV16 genomes in undifferentiated cells.

IMPORTANCE

HPV16 replicates in differentiating epithelia and can cause cancer. How HPV16 maintains its genome in undifferentiated cells at a low but constant level is not well understood but may be relevant for the immunological escape of HPV16 in the basal layers of the infected epithelium. This study demonstrates that the expression of the viral E8^E2C protein, which is a potent inhibitor of viral replication in undifferentiated cells, is driven by a separate promoter. The E8 promoter is both positively and negatively regulated by viral proteins and thus most likely acts as a sensor and modulator of viral copy number.

The transcription map of HPV11 in U2OS cells adequately reflects the initial and stable replication phases of the viral genome

BACKGROUND

Although prophylactic vaccines have been developed against HPV6, HPV11, HPV16 and HPV18 there is the clear unmet medical need in order to justify the development of drugs targeting human papillomavirus replication. The native host cells of HPVs are human primary keratinocytes which can be cultivated in raft cultures. However, this method is difficult to use in high-throughput screening assays and the need for a cost-effective cellular system for screening potential anti-HPV drug candidates during all stages of HPV genome replication remains.

METHODS

U2OS cells were transfected with HPV11 wt or E8- minicircle genomes and their gene expression was studied via 3' RACE, 5' RACE or via real time PCR methods. The DNA replication of these genomes was detected by Southern blot methods.

RESULTS

The analysis of HPV11 transcripts in U2OS cells showed that the patterns of promoter use, splice sites and polyadenylation cleavage sites are identical to those previously characterized in human HPV-related lesions, human squamous carcinoma cell lines (e.g., SSC-4) and laryngeal papillomas. Transcriptional initiation from the three previously described HPV11 promoters in the E6 and E7 ORFs (P90, P264, and P674-714) were functional, and these promoters were used together with two promoter regions in the E1 ORF (P1092 and P1372). Mutating the E8 ORF ATG start codon to ACG eliminated the translation of fusion proteins from the E8 ORF coupled to E1 and E2 proteins C-terminal sequences, leading to the de-repression of gene expression (particularly from the P1092 promoter) and to the activation of genome replication. These data suggested that the expression of the functional E8^E2 protein is used to control viral gene expression and copy number of the HPV11 genome. The analysis of HPV11 E1 expression plasmids showed that the E6/E7 region, together with the E1 coding region, is crucial for the production of functionally active E1 protein.

CONCLUSIONS

The data presented in this paper suggest that in human osteosarcoma cell line U2OS the gene expression pattern of the HPV11 truly reflect the expression profile of the replicating HPV genome and therefore this cellular system is suitable for drug development program targeting HPV replication.