Tumourigenesis driven by the human papillomavirus type 16 Asian-American e6 variant in a three-dimensional keratinocyte model

Is High-Risk Sexual Behavior a Risk Factor for Oropharyngeal Cancer?

(1) Background: Several lines of evidence established a link between high-risk (HR) sexual behavior (SB), the persistence of human papillomavirus (HPV) DNA in saliva, and the presence of oncogenic HR-HPV subtypes in oropharyngeal squamous cell carcinoma (OPSCC). A highly influential case-control study by D'Souza et al. comparing OPSCC patients and ENT patients with benign diseases (hospital controls) established HR-SB as a putative etiological risk factor for OPSCC. Aiming to replicate their findings in a nested case-control study of OPSCC patients and propensity score (PS)-matched unaffected controls from a large population-based German cohort study, we here demonstrate discrepant findings regarding HR-SB in OPSCC. (2) Methods: According to the main risk factors for HNSCC (age, sex, tobacco smoking, and alcohol consumption) PS-matched healthy controls invited from the population-based cohort study LIFE and HNSCC (including OPSCC) patients underwent interviews, using AUDIT and Fagerström, as well as questionnaires asking for SB categories as published. Afterwards, by newly calculating PSs for the same four risk factors, we matched each OPSCC patient with two healthy controls and compared responses utilizing chi-squared tests and logistic regression. (3) Results: The HNSCC patients and controls showed significant differences in sex distribution, chronologic age, tobacco-smoking history (pack years), and alcohol dependence (based on AUDIT score). However, PS-matching decreased the differences between OPSCC patients and controls substantially. Despite confirming that OPSCC patients were more likely to self-report their first sexual intercourse before age 18, we found no association between OPSCC and HR-SB, neither for practicing oral-sex, having an increased number of oral- or vaginal-sex partners, nor for having casual sex or having any sexually transmitted disease. (4) Conclusions: Our data, by showing a low prevalence of HR-SB in OPSCC patients, confirm findings from other European studies that differ substantially from North American case-control studies. HR-SB alone may not add excess risk for developing OPSCC.

Single cell transcriptomic analysis of HPV16-infected epithelium identifies a keratinocyte subpopulation implicated in cancer

Persistent HPV16 infection is a major cause of the global cancer burden. The viral life cycle is dependent on the differentiation program of stratified squamous epithelium, but the landscape of keratinocyte subpopulations which support distinct phases of the viral life cycle has yet to be elucidated. Here, single cell RNA sequencing of HPV16 infected compared to uninfected organoids identifies twelve distinct keratinocyte populations, with a subset mapped to reconstruct their respective 3D geography in stratified squamous epithelium. Instead of conventional terminally differentiated cells, an HPV-reprogrammed keratinocyte subpopulation (HIDDEN cells) forms the surface compartment and requires overexpression of the ELF3/ESE-1 transcription factor. HIDDEN cells are detected throughout stages of human carcinogenesis including primary human cervical intraepithelial neoplasias and HPV positive head and neck cancers, and a possible role in promoting viral carcinogenesis is supported by TCGA analyses. Single cell transcriptome information on HPV-infected versus uninfected epithelium will enable broader studies of the role of individual keratinocyte subpopulations in tumor virus infection and cancer evolution.

The E6 Oncoprotein of HPV16 AA-c Variant Regulates Cell Migration through the MINCR/miR-28-5p/RAP1B Axis

The E6 oncoprotein of HPV16 variants differentially alters the transcription of the genes involved in migration and non-coding RNAs such as lncRNAs. The role of the lncRNA MINCR in cervical cancer and its relationship with variants of oncogenic HPV remain unknown. Therefore, the objective of this study was to analyze the effect of the E6 oncoprotein of the AA-c variant of HPV16 in cell migration through the MINCR/miR-28-5p/RAP1B axis. To explore the functional role of MINCR in CC, we used an in vitro model of C33-A cells with exogenous expression of the E6 oncoprotein of the AA-c variant of HPV16. Interfering RNAs performed MINCR silencing, and the expression of miR-28-5p and RAP1B mRNA was analyzed by RT-qPCR. We found that C33-A/AA-c cells expressed MINCR 8-fold higher compared to the control cells. There is an inverse correlation between the expression of miR-28-5p and RAP1B in C33-A/AA-c cells. Our results suggest that MINCR might regulate the expression of RAP1B through the inhibition of miR-28-5p in CC cells expressing the E6 oncoprotein of HPV16 AA-c. We report, for the first time, that the MINCR/miR-28-5p/RAP1B axis positively regulates cell migration in CC-derived cells that express the E6 oncoprotein of the AA-c variant of HPV16.

The Potentials and Pitfalls of a Human Cervical Organoid Model Including Langerhans Cells

Three-dimensional cell culturing to capture a life-like experimental environment has become a versatile tool for basic and clinical research. Mucosal and skin tissues can be grown as “organoids” in a petri dish and serve a wide variety of research questions. Here, we report our experience with human cervical organoids which could also include an immune component, e.g., Langerhans cells. We employ commercially available human cervical keratinocytes and fibroblasts as well as a myeloid cell line matured and purified into langerin-positive Langerhans cells. These are then seeded on a layer of keratinocytes with underlying dermal equivalent. Using about 10-fold more than the reported number in healthy cervical tissue (1–3%), we obtain differentiated cervical epithelium after 14 days with ~1% being Langerhans cells. We provide a detailed protocol for interested researchers to apply the described “aseptic” organoid model for all sorts of investigations—with or without Langerhans cells.

Virus-Host Protein-Protein Interactions between Human Papillomavirus 16 E6 A1 and D2/D3 Sub-Lineages: Variances and Similarities

High-risk strains of human papillomavirus are causative agents for cervical and other mucosal cancers, with type 16 being the most frequent. Compared to the European Prototype (EP; A1), the Asian-American (AA; D2/D3) sub-lineage seems to have increased abilities to promote carcinogenesis. Here, we studied protein-protein interactions (PPIs) between host proteins and sub-lineages of the key transforming E6 protein. We transduced human keratinocyte with EP or AA E6 genes and co-immunoprecipitated E6 proteins along with interacting cellular proteins to detect virus-host binding partners. AAE6 and EPE6 may have unique PPIs with host cellular proteins, conferring gain or loss of function and resulting in varied abilities to promote carcinogenesis. Using liquid chromatography-mass spectrometry and stringent interactor selection criteria based on the number of peptides, we identified 25 candidates: 6 unique to AAE6 and EPE6, along with 13 E6 targets common to both. A novel approach based on pathway selection discovered 171 target proteins: 90 unique AAE6 and 61 unique EPE6 along with 20 common E6 targets. Interpretations were made using databases, such as UniProt, BioGRID, and Reactome. Detected E6 targets were differentially implicated in important hallmarks of cancer: deregulating Notch signaling, energetics and hypoxia, DNA replication and repair, and immune response.

An epithelial organoid model with Langerhans cells for assessing virus-host interactions

Persistent infection with oncogenic human papillomavirus (HPV) may lead to cancer in mucosal and skin tissue. Consequently, HPV must have developed strategies to escape host immune surveillance. Nevertheless, most HPV infections are cleared by the infected host. Our laboratory investigates Langerhans cells (LCs), acting at the interface between innate and adaptive immunity. We hypothesize that this first line of defence is vital for potential HPV elimination. As an alternative to animal models, we use smaller-scale epithelial organoids grown from human primary keratinocytes derived from various anatomical sites. This approach is amenable to large sample sizes-an essential aspect for scientific rigour and statistical power. To evaluate LCs phenotypically and molecularly during the viral life cycle and onset of carcinogenesis, we have included an engineered myeloid cell line with the ability to acquire an LC phenotype. This model is accurately tailored for the crucial time-window of early virus elimination in a complex organism and will shed more light on our long-standing research question of how naturally occurring HPV variants influence disease development. It may also be applied to other microorganism-host interaction research or enquiries of epithelium immunobiology. Finally, our continuously updated pathogen-host analysis tool enables state-of-the-art bioinformatics analyses of next-generation sequencing data. This article is part of the theme issue 'Silent cancer agents: multi-disciplinary modelling of human DNA oncoviruses'.

Epithelial stratification shapes infection dynamics

Infections of stratified epithelia contribute to a large group of common diseases, such as dermatological conditions and sexually transmitted diseases. To investigate how epithelial structure affects infection dynamics, we develop a general ecology-inspired model for stratified epithelia. Our model allows us to simulate infections, explore new hypotheses and estimate parameters that are difficult to measure with tissue cell cultures. We focus on two contrasting pathogens: Chlamydia trachomatis and Human papillomaviruses (HPV). Using cervicovaginal parameter estimates, we find that key infection symptoms can be explained by differential interactions with the layers, while clearance and pathogen burden appear to be bottom-up processes. Cell protective responses to infections (e.g. mucus trapping) generally lowered pathogen load but there were specific effects based on infection strategies. Our modeling approach opens new perspectives for 3D tissue culture experimental systems of infections and, more generally, for developing and testing hypotheses related to infections of stratified epithelia.

Isolation of Biopsy-Derived, Human Cervical Keratinocytes Propagated as Monolayer and Organoid Cultures

The successful isolation and propagation of patient-derived keratinocytes from cervical lesions constitute a more appropriate model of cervical disease than traditional cervical cancer-derived cell lines such as SiHa and CaSki. Our aim was to streamline the growth of patient-obtained, cervical keratinocytes into a reproducible process. We performed an observational case series study with 60 women referred to colposcopy for a diagnostic biopsy. Main outcome measures were how many samples could be passaged at least once (n = 11), and where enough cells could be established, to precisely define their proliferation profile over time (n = 3). Altering cell culture conditions over those reported by other groups markedly improved outcomes. We were also successful in making freeze backs which could be resuscitated to successfully propagate multi-layered, organoids from cervical keratinocytes (n = 3). For best results, biopsy-intrinsic factors such as size and tissue digestion appear to be major variables. This seems to be the first systematic report with a well characterized and defined sample size, detailed protocol, and carefully assessed cell yield and performance. This research is particularly impactful for constituting a sample repository-on-demand for appropriate disease modelling and drug screening under the umbrella of personalized health.

Differences in the viral genome between HPV-positive cervical and oropharyngeal cancer

Human papillomavirus (HPV)-driven oropharyngeal cancer incidence in the United States has steadily increased in the past decades and has now become the most frequently diagnosed HPV-associated cancer type, surpassing cervical cancer. Variations in the HPV genome correlate with tumorigenic risk, and the distribution of genetic variants is extensively studied in cervical cancer, but very little is known about new mutations or the distribution of HPV types and variants in oropharyngeal cancer. Here we present an archival tissue cohort study that compares genomic characteristics of HPV associated with cervical versus oropharyngeal tumors using DNA sequence analysis. We found HPV16 to be more prevalent in oropharyngeal samples than in cervical samples (91.2% versus 52.9%), while HPV18 (1.5% versus 18.2%) and HPV45 (0.7% versus 9.9%) were much less prevalent. Differences between cervix and oropharynx in HPV16 variants distribution were more subtle, but the combined European + Asian (EUR+AS) variant group was more prevalent (90.2% versus 71.4%), while the American Asian 1 + American Asian 2 (AA1+AA2) variant group was much less prevalent (4.4% versus 22.5%) in oropharyngeal cancers. HPV prevalence in oropharyngeal cancers showed an increasing trend from 60% in 2003 to 80% in 2016. We also identified over nine times more nonsynonymous mutations in the HPV E6 gene amplified from oropharyngeal samples, but for E7 the difference in mutation rates between the two anatomical locations was not significant. Overall, we showed that HPV genome in oropharyngeal cancer presents important differences when compared to cervical cancer and this may explain the distinct pathomechanisms and susceptibility to treatment of HPV-associated oropharyngeal cancer.

Synthetic siRNA targeting human papillomavirus 16 E6: a perspective on in vitro nanotherapeutic approaches

High-risk human papillomaviruses infect skin and mucosa, causing approximately 5% of cancers worldwide. In the search for targeted nanotherapeutic approaches, siRNAs against the viral E6 transcript have been molecules of interest but have not yet seen successful translation into the clinic. By reviewing the past approximately 15 years of in vitro literature, we identify the need for siRNA validation protocols which concurrently evaluate ranges of key treatment parameters as well as characterize downstream process restoration in a methodical, quantitative manner and demonstrate their implementation using our own data. We also reflect on the future need for more appropriate cell culture models to represent patient lesions as well as the application of personalized approaches to identify optimal treatment strategies.

E6/E7 oncogenes in epithelial suprabasal layers and estradiol promote cervical growth and ear regeneration

Tissue growth is a common characteristic of carcinogenesis and regeneration. Here we show that suprabasal expression of human papillomavirus (HPV)16 E6/E7 oncogenes in Tg(K6b-E6/E7) mice, similar to that observed in HPV-infected human tissue, and estradiol increased cervical epithelium growth and ear-hole closure efficiency. Oncogenes in combination with estradiol had a significant contribution to the proliferation of suprabasal cells of cervical epithelium that correlated with an increased expression of keratin genes. Remarkably, long-term treatments with estradiol resulted in evident cellular and tissue abnormalities indicative of a precancerous phenotype. Regenerating ear epithelium of transgenic mice also showed increased suprabasal cell proliferation and expression of keratin genes. Unexpectedly, we observed higher ear regeneration efficiency in adult than in young female mice, which was further increased by E6/E7 oncogenes. Supporting a role of estradiol in this phenomenon, ovariectomy and treatment with an estrogen receptor inhibitor caused a significant reduction in regenerative capacity. Our data suggest that Tg(K6b-E6/E7) mice are unique to mimic the initial stages of HPV-mediated cervical carcinogenesis, and ear regeneration could facilitate the elucidation of mechanisms involved.

Variation of HPV Subtypes with Focus on HPV-Infection and Cancer in the Head and Neck Region

The human papillomavirus (HPV) comprises a heterogeneous group of double-strand DNA viruses with variable potential to infect human epithelial cells and trigger neoplastic transformation. Its 8 kb genome encodes proteins required for virus replication and self-organized formation of infectious particles but also for early proteins E6 and E7 able to trigger neoplastic transformation. E6 and E7 of high-risk (HR) HPV subtypes can bind to p53 or release E2F and abrogate replication control. Due to variable amino acid sequence (AAS) in the binding sites of E6 and E7 particular HR-HPV variants within subtypes are essentially heterogeneous in efficacy triggering neoplastic transformation and cancer development. This could explain differences in the clinical course of HPV-driven head and neck cancer.

Different Association of Human Papillomavirus 16 Variants with Early and Late Presentation of Cervical Cancer

The median age of cervical cancer (CC) presentation coincides with the mean age of menopause presentation (49 years) in Mexico. Here, we investigated the association between different HPV16 variants and early (≤ 49 years) or delayed (≥ 50 years) CC presentation. We conducted a case-case study that included 462 CCs, 386 squamous cell carcinomas (SCC), 63 adenocarcinomas (ACC), and 13 additional cell types. Variants were identified by PCR and DNA sequencing. The risk conferred by each variant for developing CC earlier than 50 years was analyzed using a univariate logistic regression model considering old-aged patients (≥ 50 years) and non-HPV16 cases as the reference variables. Overall, the frequency of HPV16 was 50.9%, and the only identified variants were the European A1/2 (31.2%) and the Asian-American D2 (10.8%), and D3 (8.9%). D2 was mainly associated with ≤ 49-year-old patients (15.9%); A1/2 was uniformly distributed between the two age groups (~31%), whereas D3 increased with age to a frequency of 11.8% in the older group. Only the D2 variant conferred a 3.3-fold increase in the risk of developing CC before 50 years of age (OR = 3.3, 95% CI = 1.7-6.6, p < 0.001) in relation with non-HPV16 cases. Remarkably, this risk was higher for ACC (OR = 6.0, 95% CI = 1.1-33, p < 0.05) than for SCC (OR = 2.8, 95% CI = 1.3-5.9, p < 0.01). Interestingly, when analyzing only the HPV16-positive CC, D2 increases (OR = 2.5, 95% CI = 1.2-5, p < 0.05) and D3 decreases (OR = 0.45, 95% CI 0.2-0.9, p < 0.05) the risk to develop CC before 50 years old in relation with A1/2 variant. These results indicated that D2 variant is associated with early and D3 with delayed CC presentation, whereas A1/2 variant was uniformly distributed between the two age groups.

Molecular mechanisms of HPV mediated neoplastic progression

Human Papillomavirus is the major etiological agent in the development of cervical cancer but not a sufficient cause. Despite significant research, the underlying mechanisms of progression from a low-grade squamous intraepithelial lesion to high grade squamous intraepithelial lesion are yet to be understood. Deregulation of viral gene expression and host genomic instability play a central role in virus-mediated carcinogenesis. Key events such as viral integration and epigenetic modifications may lead to the deregulation of viral and host gene expression. This review has summarized the available literature to describe the possible mechanism and role of viral integration in mediating carcinogenesis. HPV integration begins with DNA damage or double strand break induced either by oxidative stress or HPV proteins and the subsequent steps are driven by the DNA damage responses. Inflammation and oxidative stress could be considered as cofactors in stimulating viral integration and deregulation of cellular and viral oncogenes during the progression of cervical carcinoma. All these events together with the host and viral genetic and epigenetic modifications in neoplastic progression have also been reviewed which may be relevant in identifying a new preventive therapeutic strategy. In the absence of therapeutic intervention for HPV-infected individuals, future research focus should be directed towards preventing and reversing of HPV integration. DNA damage response, knocking out integrated HPV sequences, siRNA approach, modulating the selection mechanism of cells harboring integrated genomes and epigenetic modifiers are the possible therapeutic targets.

Functional variants of human papillomavirus type 16 demonstrate host genome integration and transcriptional alterations corresponding to their unique cancer epidemiology

BACKGROUND

Human papillomaviruses (HPVs) are a worldwide burden as they are a widespread group of tumour viruses in humans. Having a tropism for mucosal tissues, high-risk HPVs are detected in nearly all cervical cancers. HPV16 is the most common high-risk type but not all women infected with high-risk HPV develop a malignant tumour. Likely relevant, HPV genomes are polymorphic and some HPV16 single nucleotide polymorphisms (SNPs) are under evolutionary constraint instigating variable oncogenicity and immunogenicity in the infected host.

RESULTS

To investigate the tumourigenicity of two common HPV16 variants, we used our recently developed, three-dimensional organotypic model reminiscent of the natural HPV infectious cycle and conducted various "omics" and bioinformatics approaches. Based on epidemiological studies we chose to examine the HPV16 Asian-American (AA) and HPV16 European Prototype (EP) variants. They differ by three non-synonymous SNPs in the transforming and virus-encoded E6 oncogene where AAE6 is classified as a high- and EPE6 as a low-risk variant. Remarkably, the high-risk AAE6 variant genome integrated into the host DNA, while the low-risk EPE6 variant genome remained episomal as evidenced by highly sensitive Capt-HPV sequencing. RNA-seq experiments showed that the truncated form of AAE6, integrated in chromosome 5q32, produced a local gene over-expression and a large variety of viral-human fusion transcripts, including long distance spliced transcripts. In addition, differential enrichment of host cell pathways was observed between both HPV16 E6 variant-containing epithelia. Finally, in the high-risk variant, we detected a molecular signature of host chromosomal instability, a common property of cancer cells.

CONCLUSIONS

We show how naturally occurring SNPs in the HPV16 E6 oncogene cause significant changes in the outcome of HPV infections and subsequent viral and host transcriptome alterations prone to drive carcinogenesis. Host genome instability is closely linked to viral integration into the host genome of HPV-infected cells, which is a key phenomenon for malignant cellular transformation and the reason for uncontrolled E6 oncogene expression. In particular, the finding of variant-specific integration potential represents a new paradigm in HPV variant biology.

Changes in global gene expression profiles induced by HPV 16 E6 oncoprotein variants in cervical carcinoma C33-A cells

We analyzed the effects of the expression of HPV 16 E6 oncoprotein variants (AA-a, AA-c, E-A176/G350, E-C188/G350, E-G350), and the E-Prototype in global gene expression profiles in an in vitro model. E6 gene was cloned into an expression vector fused to GFP and was transfected in C33-A cells. Affymetrix GeneChip Human Transcriptome Array 2.0 platform was used to analyze the expression of over 245,000 coding transcripts. We found that HPV16 E6 variants altered the expression of 387 different genes in comparison with E-Prototype. The altered genes are involved in cellular processes related to the development of cervical carcinoma, such as adhesion, angiogenesis, apoptosis, differentiation, cell cycle, proliferation, transcription and protein translation. Our results show that polymorphic changes in HPV16 E6 natural variants are sufficient to alter the overall gene expression profile in C33-A cells, explaining in part the observed differences in oncogenic potential of HPV16 variants.

Genetic variations of human papillomavirus type 16: implications for cervical carcinogenesis

Human papillomaviruses (HPVs) are the causative agent of cervical cancer, and among approximately 15 high-risk genotypes, HPV16 accounts for more than half the cases of cervical cancer worldwide. Recent progress in determining HPV genomic sequences from clinical samples has revealed a wide variety in HPV16 genome sequences, and has allowed for comprehensive classification of intratype HPV16 variants. These consist of four variant lineages containing nucleotide variations in 1.0%-10.0% of the complete viral genome sequence. Epidemiological data suggest that the non-European-Asian lineages of HPV16 entail a higher risk of progression to invasive cervical cancer than the European-Asian lineage. Deep sequencing analysis has recently demonstrated that HPV16 genome sequences are highly homogeneous in individual clinical specimens compared with those of RNA viruses. However, an extremely sensitive PCR method, differential DNA denaturation PCR, has detected hypermutations from C to T or G to A in the E2 gene and the long control region of the HPV16 genome, which suggests the involvement of cellular apolipoprotein B mRNA editing enzyme, catalytic polypeptide-like (APOBEC) proteins in this hypermutation. The quasispecies status of the HPV16 genome in the infected cervix may affect the development of cervical cancer and warrants further investigation.