Dangerous Liaisons: Long-Term Replication with an Extrachromosomal HPV Genome

Roles of human papillomavirus in cancers: oncogenic mechanisms and clinical use

Human papillomaviruses, particularly high-risk human papillomaviruses, have been universally considered to be associated with the oncogenesis and progression of various cancers. The genome of human papillomaviruses is circular, double-stranded DNA that encodes early and late proteins. Each of the proteins is of crucial significance in infecting the epithelium of host cells persistently and supporting viral genome integrating into host cells. Notably, E6 and E7 proteins, classified as oncoproteins, trigger the incidence of cancers by fostering cell proliferation, hindering apoptosis, evading immune surveillance, promoting cell invasion, and disrupting the balance of cellular metabolism. Therefore, targeting human papillomaviruses and decoding molecular mechanisms by which human papillomaviruses drive carcinogenesis are of great necessity to better treat human papillomaviruses-related cancers. Human papillomaviruses have been applied clinically to different facets of human papillomavirus-related cancers, including prevention, screening, diagnosis, treatment, and prognosis. Several types of prophylactic vaccines have been publicly utilized worldwide and have greatly decreased the occurrence of human papillomavirus-related cancers, which have benefited numerous people. Although various therapeutic vaccines have been developed and tested clinically, none of them have been officially approved to date. Enhancing the efficacy of vaccines and searching for innovative technologies targeting human papillomaviruses remain critical challenges that warrant continuous research and attention in the future.

Fibroblasts regulate the transcriptional signature of human papillomavirus-positive keratinocytes

Persistent human papillomavirus (HPV) infection is necessary but insufficient for viral oncogenesis. Additional contributing co-factors, such as immune evasion and viral integration have been implicated in HPV-induced cancer progression. It is widely accepted that HPV + keratinocytes require co-culture with fibroblasts to maintain viral DNA as episomes. How fibroblasts regulate viral episome maintenance is a critical knowledge gap. Here we present comprehensive RNA sequencing and proteomic analysis demonstrating that coculture with fibroblasts is supportive of the viral life cycle, and is confirmatory of previous observations. Novel observations suggest that errors in "cross-talk" between fibroblasts and infected keratinocytes may regulate HPV integration and drive oncogenic progression. Our co-culture models offer new insights into HPV-related transformation mechanisms.

Fibroblast stromal support model for predicting human papillomavirus-associated cancer drug responses

Currently, there are no specific antiviral therapeutic approaches targeting Human papillomaviruses (HPVs), which cause around 5% of all human cancers. Specific antiviral reagents are particularly needed for HPV-related oropharyngeal cancers (HPV+OPCs) whose incidence is increasing and for which there are no early diagnostic tools available. We and others have demonstrated that the estrogen receptor alpha (ERα) is overexpressed in HPV+OPCs, compared to HPV-negative cancers in this region, and that these elevated levels are associated with an improved disease outcome. Utilizing this HPV+-specific overexpression profile, we previously demonstrated that estrogen attenuates the growth and cell viability of HPV+ keratinocytes and HPV+ cancer cells in vitro. Expansion of this work in vivo failed to replicate this sensitization. The role of stromal support from the tumor microenvironment (TME) has previously been tied to both the HPV lifecycle and in vivo therapeutic responses. Our investigations revealed that in vitro co-culture with fibroblasts attenuated HPV+-specific estrogen growth responses. Continuing to monopolize on the HPV+-specific overexpression of ERα, our co-culture models then assessed the suitability of the selective estrogen receptor modulators (SERMs), raloxifene and tamoxifen, and showed growth attenuation in a variety of our models to one or both of these drugs in vitro. Utilization of these SERMs in vivo closely resembled the sensitization predicted by our co-culture models. Therefore, the in vitro fibroblast co-culture model better predicts in vivo responses. We propose that utilization of our co-culture in vitro model can accelerate cancer therapeutic drug discovery.

IMPORTANCE

Human papillomavirus-related cancers (HPV+ cancers) remain a significant public health concern, and specific clinical approaches are desperately needed. In translating drug response data from in vitro to in vivo, the fibroblasts of the adjacent stromal support network play a key role. Our study presents the utilization of a fibroblast 2D co-culture system to better predict translational drug assessments for HPV+ cancers. We also suggest that this co-culture system should be considered for other translational approaches. Predicting even a portion of treatment paradigms that may fail in vivo with a co-culture model will yield significant time, effort, resource, and cost efficiencies.

Structure and transcription of integrated HPV DNA in vulvar carcinomas

HPV infections are associated with a fraction of vulvar cancers. Through hybridization capture and DNA sequencing, HPV DNA was detected in five of thirteen vulvar cancers. HPV16 DNA was integrated into human DNA in three of the five. The insertions were in introns of human NCKAP1, C5orf67, and LRP1B. Integrations in NCKAP1 and C5orf67 were flanked by short direct repeats in the human DNA, consistent with HPV DNA insertions at sites of abortive, staggered, endonucleolytic incisions. The insertion in C5orf67 was present as a 36 kbp, human-HPV-hetero-catemeric DNA as either an extrachromosomal circle or a tandem repeat within the human genome. The human circularization/repeat junction was defined at single nucleotide resolution. The integrated viral DNA segments all retained an intact upstream regulatory region and the adjacent viral E6 and E7 oncogenes. RNA sequencing revealed that the only HPV genes consistently transcribed from the integrated viral DNAs were E7 and E6*I. The other two HPV DNA+ tumors had coinfections, but no evidence for integration. HPV-positive and HPV-negative vulvar cancers exhibited contrasting human, global gene expression patterns partially overlapping with previously observed differences between HPV-positive and HPV-negative cervical and oropharyngeal cancers. A substantial fraction of the differentially expressed genes involved immune system function. Thus, transcription and HPV DNA integration in vulvar cancers resemble those in other HPV-positive cancers. This study emphasizes the power of hybridization capture coupled with DNA and RNA sequencing to identify a broad spectrum of HPV types, determine human genome integration status of viral DNAs, and elucidate their structures.

HPV, HBV, and HIV-1 Viral Integration Site Mapping: A Streamlined Workflow from NGS to Genomic Insights of Carcinogenesis

Viral integration within the host genome plays a pivotal role in carcinogenesis. Various disruptive mechanisms are involved, leading to genomic instability, mutations, and DNA damage. With next-generation sequencing (NGS), we can now precisely identify viral and host genomic breakpoints and chimeric sequences, which are useful for integration site analysis. In this study, we evaluated a commercial hybrid capture NGS panel specifically designed for detecting three key viruses: HPV, HBV, and HIV-1. We also tested workflows for Viral Hybrid Capture (VHC) and Viral Integration Site (VIS) analysis, leveraging customized viral databases in CLC Microbial Genomics. By analyzing sequenced data from virally infected cancer cell lines (including SiHa, HeLa, CaSki, C-33A, DoTc2, 2A3, SCC154 for HPV; 3B2, SNU-182 for HBV; and ACH-2 for HIV-1), we precisely pinpointed viral integration sites. The workflow also highlighted disrupted and neighboring human genes that may play a crucial role in tumor development. Our results included informative virus-host read mappings, genomic breakpoints, and integration circular plots. These visual representations enhance our understanding of the integration process. In conclusion, our seamless end-to-end workflow bridges the gap in understanding viral contributions to cancer development, paving the way for improved diagnostics and treatment strategies.

Knowledge on Human Papillomavirus Infections, Cancer Biology, Immune Interactions, Vaccination Coverage and Common Treatments: A Comprehensive Review

Human papillomavirus (HPV) is a circular, double-stranded DNA virus and recognized as the most prevalent sexually transmitted infectious agent worldwide. The HPV life cycle encompasses three primary stages. First, the virus infiltrates the basal cells of the stratified epidermis. Second, there is a low-level expression of viral genes and preservation of the viral genome in the basal layer. Lastly, productive replication of HPV occurs in differentiated cells. An effective immune response, involving various immune cells, including innate immunity, keratinocytes, dendritic cells, and natural killer T cells, is instrumental in clearing HPV infection and thwarting the development of HPV-associated tumors. Vaccines have demonstrated their efficacy in preventing genital warts, high-grade precancerous lesions, and cancers in females. In males, the vaccines can also aid in preventing genital warts, anal precancerous lesions, and cancer. This comprehensive review aims to provide a thorough and detailed exploration of HPV infections, delving into its genetic characteristics, life cycle, pathogenesis, and the role of high-risk and low-risk HPV strains. In addition, this review seeks to elucidate the intricate immune interactions that govern HPV infections, spanning from innate immunity to adaptive immune responses, as well as examining the evasion mechanisms used by the virus. Furthermore, the article discusses the current landscape of HPV vaccines and common treatments, contributing to a holistic understanding of HPV and its associated diseases.

Cis-regulatory effect of HPV integration is constrained by host chromatin architecture in cervical cancers

Human papillomavirus (HPV) infections are the primary drivers of cervical cancers, and often HPV DNA gets integrated into the host genome. Although the oncogenic impact of HPV encoded genes is relatively well known, the cis-regulatory effect of integrated HPV DNA on host chromatin structure and gene regulation remains less understood. We investigated genome-wide patterns of HPV integrations and associated host gene expression changes in the context of host chromatin states and topologically associating domains (TADs). HPV integrations were significantly enriched in active chromatin regions and depleted in inactive ones. Interestingly, regardless of chromatin state, genomic regions flanking HPV integrations showed transcriptional upregulation. Nevertheless, upregulation (both local and long-range) was mostly confined to TADs with integration, but not affecting adjacent TADs. Few TADs showed recurrent integrations associated with overexpression of oncogenes within them (e.g. MYC, PVT1, TP63 and ERBB2) regardless of proximity. Hi-C and 4C-seq analyses in cervical cancer cell line (HeLa) demonstrated chromatin looping interactions between integrated HPV and MYC/PVT1 regions (~ 500 kb apart), leading to allele-specific overexpression. Based on these, we propose HPV integrations can trigger multimodal oncogenic activation to promote cancer progression.

Human papillomaviruses: Knowns, mysteries, and unchartered territories

There has been an explosion in the number of papillomaviruses that have been identified and fully sequenced. Yet only a minute fraction of these has been studied in any detail. Most of our molecular research efforts have focused on the E6 and E7 proteins of "high-risk," cancer-associated human papillomaviruses (HPVs). Interactions of the high-risk HPV E6 and E7 proteins with their respective cellular targets, the p53 and the retinoblastoma tumor suppressors, have been investigated in minute detail. Some have thus questioned if research on papillomaviruses remains an exciting and worthwhile area of investigation. However, fundamentally new insights on the biological activities and cellular targets of the high-risk HPV E6 and E7 proteins have been discovered and previously unstudied HPVs have been newly associated with human diseases. HPV infections continue to be an important cause of human morbidity and mortality and since there are no antivirals to combat HPV infections, research on HPVs should remain attractive to new investigators and biomedical funding agencies, alike.

Human Papillomavirus Genome Copy Number Is Maintained by S-Phase Amplification, Genome Loss to the Cytosol during Mitosis, and Degradation in G1 Phase

The current model of human papillomavirus (HPV) replication is comprised of three modes of replication. Following infectious delivery, the viral genome is amplified during the establishment phase to reach up to some hundred copies per cell. The HPV genome copy number remains constant during the maintenance stage. The differentiation of infected cells induces HPV genome amplification. Using highly sensitive in situ hybridization (DNAscope) and freshly HPV16-infected as well as established HPV16-positive cell lines, we observed that the viral genome is amplified in each S phase of undifferentiated keratinocytes cultured as monolayers. The nuclear viral genome copy number is reset to pre-S-phase levels during mitosis. The majority of the viral genome fails to tether to host chromosomes and is lost to the cytosol. Cytosolic viral genomes gradually decrease during cell cycle progression. The loss of cytosolic genomes is blocked in the presence of NH4Cl or other drugs that interfere with lysosomal acidification, suggesting the involvement of autophagy in viral genome degradation. These observations were also made with HPV31 cell lines obtained from patient samples. Cytosolic viral genomes were not detected in UMSCC47 cells carrying integrated HPV16 DNA. Analyses of organotypic raft cultures derived from keratinocytes harboring episomal HPV16 revealed the presence of cytosolic viral genomes as well. We conclude that HPV maintains viral genome copy numbers by balancing viral genome amplification during S phase with the loss of viral genomes to the cytosol during mitosis. It seems plausible that restrictions to viral genome tethering to mitotic chromosomes reset genome copy numbers in each cell cycle. IMPORTANCE HPV genome maintenance is currently thought to be achieved by regulating the expression and activity of the viral replication factors E1 and E2. In addition, the E8^E2 repressor has been shown to be important for restricting genome copy numbers by competing with E1 and E2 for binding to the viral origin of replication and by recruiting repressor complexes. Here, we demonstrate that the HPV genome is amplified in each S phase. The nuclear genome copy number is reset during mitosis by a failure of the majority of the genomes to tether to mitotic chromosomes. Rather, HPV genomes accumulate in the cytoplasm of freshly divided cells. Cytosolic viral DNA is degraded in G1 in a lysosome-dependent manner, contributing to the genome copy reset. Our data imply that the mode of replication during establishment and maintenance is the same and further suggest that restrictions to genome tethering significantly contribute to viral genome maintenance.

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.

The First Human Vulvar Intraepithelial Neoplasia Cell Line with Naturally Infected Episomal HPV18 Genome

Persistent infection with high-risk HPV leads to cervical cancers and other anogenital cancers and head and neck carcinomas in both men and women. There is no effective drug fortreating HPV infection and HPV-associated carcinomas, largely due to a lack of models of natural HPV infection and the complexity of the HPV life cycle. There are no available cell lines from vulvar, anal, or penile lesions and cancers in the field. In this study, we established the first human cell line from vulvar intraepithelial neoplasia (VIN) with naturally infected HPV18 by conditional reprogramming (CR) method. Our data demonstrated that VIN cells possessed different biological characteristics and diploid karyotypes from HPV18-positive cancer cells (HeLa). Then, we determined that VIN cells contained episomal HPV18 using approaches including the ratio of HPV E2_copy/E7_copy, rolling cycle amplification, and sequencing. The VIN cells expressed squamous epithelium-specific markers that are different from HeLa cells, a cervical adenocarcinoma cell line. When cultured under 3D air-liquid interface (ALI) system, we observed the expression of both early and late differentiation markers involucrin and filaggrin. Most importantly, we were able to detect the expression of viral late gene L1 in the cornified layer of ALI 3D culture derived from VIN cells, suggesting quite different HPV genomic status from cancer cells. We also observed progeny viral particles under transmission electron microscopy (TEM) in ALI 3D cultures, confirming the episomal HPV18 genome and active viral life cycle in the new cell line. To our knowledge, this is the first human VIN cell line with naturally infected HPV18 genome and provides a valuable model for HPV biology studies, HPV-associated cancer initiation and progression, and drug-screening platforms.

HPV Integration Site Mapping: A Rapid Method of Viral Integration Site (VIS) Analysis and Visualization Using Automated Workflows in CLC Microbial Genomics

Human papillomavirus (HPV) integration within the host genome may contribute to carcinogenesis through various disruptive mechanisms. With next-generation sequencing (NGS), identification of viral and host genomic breakpoints and chimeric sequences are now possible. However, a simple, streamlined bioinformatics workflow has been non-existent until recently. Here, we tested two new, automated workflows in CLC Microbial Genomics, i.e., Viral Hybrid Capture (VHC) Data Analysis and Viral Integration Site (VIS) Identification for software performance and efficiency. The workflows embedded with HPV and human reference genomes were used to analyze a publicly available NGS dataset derived from pre- and cancerous HPV+ cervical cytology of 21 Gabonese women. The VHC and VIS workflow median runtimes were 19 and 7 min per sample, respectively. The VIS dynamic graphical outputs included read mappings, virus-host genomic breakpoints, and virus-host integration circular plots. Key findings, including disrupted and nearby genes, were summarized in an auto-generated report. Overall, the VHC and VIS workflows proved to be a rapid and accurate means of localizing viral-host integration site(s) and identifying disrupted and neighboring human genes. Applying HPV VIS-mapping to pre- or invasive tumors will advance our understanding of viral oncogenesis and facilitate the discovery of prognostic biomarkers and therapeutic targets.

ATM Pathway Is Essential for HPV-Positive Human Cervical Cancer-Derived Cell Lines Viability and Proliferation

Infection with some mucosal human papillomavirus (HPV) types is the etiological cause of cervical cancer and of a significant fraction of vaginal, vulvar, anal, penile, and head and neck carcinomas. DNA repair machinery is essential for both HPV replication and tumor cells survival suggesting that cellular DNA repair machinery may play a dual role in HPV biology and pathogenesis. Here, we silenced genes involved in DNA Repair pathways to identify genes that are essential for the survival of HPV-transformed cells. We identified that inhibition of the ATM/CHK2/BRCA1 axis selectively affects the proliferation of cervical cancer-derived cell lines, without altering normal primary human keratinocytes (PHK) growth. Silencing or chemical inhibition of ATM/CHK2 reduced the clonogenic and proliferative capacity of cervical cancer-derived cells. Using PHK transduced with HPV16 oncogenes we observed that the effect of ATM/CHK2 silencing depends on the expression of the oncogene E6 and on its ability to induce p53 degradation. Our results show that inhibition of components of the ATM/CHK2 signaling axis reduces p53-deficient cells proliferation potential, suggesting the existence of a synthetic lethal association between CHK2 and p53. Altogether, we present evidence that synthetic lethality using ATM/CHK2 inhibitors can be exploited to treat cervical cancer and other HPV-associated tumors.

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.

Diabetes associated with HPV infection in women aged over 50 years: A cross-sectional study from China's largest academic woman's hospital

BACKGROUND

Metabolic disturbances and immune alterations caused by diabetes are not just bystanders of HPV infection, but the conclusion that diabetes increases the risk of HPV infection requires more clinical epidemiological evidence to confirm. Our aim was to evaluate the association of diabetes with HPV infection risk in female patients aged over 50 years in the cervical clinic.

METHODS

We conducted a cross-sectional study of 6402 women aged over 50 years in the cervical clinic between May 2019 and March 2022 from China's largest academic woman's hospital. The quantitative-effect relationship between diabetes and HPV infection was observed by dose-response graph. Segmented multivariate logistic regression analysis was conducted to estimate the relative risk of HPV infection in diabetes patients. Multivariable predicted marginal proportions from logistic regression models were used to compute adjusted risk ratios.

RESULTS

There is a nonlinear relationship between HbA1c and the risk of HPV infection. When the HbA1c exceeds 5.7%, there is a saturation effect. After adjustment for confounders, the risk ratio for HPV infection in women with prediabetes was 1.09 (95% CI: 1.00-1.18) compared with women with HbA1c <5.7%, and the risk ratio for HPV infection in women with diabetes was 1.18 (95%). CI: 1.04-1.33). Sensitivity analysis showed that the risk ratio for HPV infection was 1.47 (95% CL: 1.07-1.91) when diabetes was associated with vaginitis. E-value analysis suggested robustness to unmeasured confounding.

CONCLUSIONS

Diabetes and prediabetes are at increased risk of coinfection with HPV in female patients aged over 50 years in the cervical clinic.

Recurrent integration of human papillomavirus genomes at transcriptional regulatory hubs

Oncogenic human papillomavirus (HPV) genomes are often integrated into host chromosomes in HPV-associated cancers. HPV genomes are integrated either as a single copy or as tandem repeats of viral DNA interspersed with, or without, host DNA. Integration occurs frequently in common fragile sites susceptible to tandem repeat formation and the flanking or interspersed host DNA often contains transcriptional enhancer elements. When co-amplified with the viral genome, these enhancers can form super-enhancer-like elements that drive high viral oncogene expression. Here we compiled highly curated datasets of HPV integration sites in cervical (CESC) and head and neck squamous cell carcinoma (HNSCC) cancers, and assessed the number of breakpoints, viral transcriptional activity, and host genome copy number at each insertion site. Tumors frequently contained multiple distinct HPV integration sites but often only one “driver” site that expressed viral RNA. As common fragile sites and active enhancer elements are cell-type-specific, we mapped these regions in cervical cell lines using FANCD2 and Brd4/H3K27ac ChIP-seq, respectively. Large enhancer clusters, or super-enhancers, were also defined using the Brd4/H3K27ac ChIP-seq dataset. HPV integration breakpoints were enriched at both FANCD2-associated fragile sites and enhancer-rich regions, and frequently showed adjacent focal DNA amplification in CESC samples. We identified recurrent integration “hotspots” that were enriched for super-enhancers, some of which function as regulatory hubs for cell-identity genes. We propose that during persistent infection, extrachromosomal HPV minichromosomes associate with these transcriptional epicenters and accidental integration could promote viral oncogene expression and carcinogenesis.