Estrogen protects females from hepatocellular carcinoma (HCC). To determine whether this protection is mediated by classic estrogen receptors, we tested HCC susceptibility in estrogen receptor-deficient mice. In contrast to a previous study, we found that diethylnitrosamine induces hepatocarcinogenesis to a significantly greater extent when females lack Esr1, which encodes Estrogen Receptor-α. Relative to wild-type littermates, Esr1 knockout females developed 9-fold more tumors. Deficiency of Esr2, which encodes Estrogen Receptor-β, did not affect liver carcinogenesis in females. Using microarrays and QPCR to examine estrogen receptor effects on hepatic gene expression patterns, we found that germline Esr1 deficiency resulted in the masculinization of gene expression in the female liver. Six of the most dysregulated genes have previously been implicated in HCC. In contrast, Esr1 deletion specifically in hepatocytes of Esr1 conditional null female mice (in which Cre was expressed from the albumin promoter) resulted in the maintenance of female-specific liver gene expression. Wild-type adult females lacking ovarian estrogen due to ovariectomy, which is known to make females susceptible to HCC, also maintained female-specific expression in the liver of females. These studies indicate that Esr1 mediates liver cancer risk, and its control of sex-specific liver gene expression involves cells other than hepatocytes.

High-risk human papillomavirus strain 16 (HPV16) causes oral and anogenital cancers through the activities of two viral oncoproteins, E6 and E7, that dysregulate the host p53 and pRb tumor suppressor pathways, respectively. The maintenance of HPV16-positive cancers requires constitutive expression of E6 and E7. Therefore, inactivating these proteins could provide the basis for an anticancer therapy. Herein we demonstrate that a subset of aspartyl protease inhibitor drugs currently used to treat HIV/AIDS cause marked reductions in HPV16 E6 and E7 protein levels using two independent cell culture models: HPV16-transformed CaSki cervical cancer cells and NIKS16 organotypic raft cultures (a 3-D HPV16-positive model of epithelial pre-cancer). Treatment of CaSki cells with some (lopinavir, ritonavir, nelfinavir, and saquinavir) but not other (indinavir and atazanavir) protease inhibitors reduced E6 and E7 protein levels, correlating with increased p53 protein levels and decreased cell viability. Long-term (>7 day) treatment of HPV16-positive NIKS16 raft cultures with saquinavir caused epithelial atrophy with no discernible effects on HPV-negative rafts, demonstrating selectivity. Saquinavir also reduced HPV16's effects on markers of the cellular autophagy pathway in NIKS16 rafts, a hallmark of HPV-driven pre-cancers. Taken together, these data suggest HIV-1 protease inhibitors be studied further in the context of treating or preventing HPV16-positive cancers.

We investigated whether in human head and neck squamous cell carcinoma (HNSCC) high levels of expression of stress keratin 17 (K17) are associated with poor survival and resistance to immunotherapy.

Papillomaviruses exhibit species-specific tropism, thereby limiting understanding and research of several aspects of HPV infection and carcinogenesis. The discovery of a murine papillomavirus (MmuPV1) provides the opportunity to study papillomavirus infections in a tractable, in vivo laboratory model. MmuPV1 infects and causes disease in the cutaneous epithelium, as well as the mucosal epithelia of the oral cavity and anogenital tract. In this report, we describe a murine model of MmuPV1 infection and neoplastic disease in the female reproductive tracts of wild-type immunocompetent FVB mice. Low-grade dysplastic lesions developed in reproductive tracts of FVB mice infected with MmuPV1 for 4 months, and mice infected for 6 months developed significantly worse disease, including squamous cell carcinoma (SCC). We also tested the contribution of estrogen and/or UV radiation (UVR), two cofactors we previously identified as being involved in papillomavirus-mediated disease, to cervicovaginal disease. Similar to HPV16 transgenic mice, exogenous estrogen treatment induced high-grade precancerous lesions in the reproductive tracts of MmuPV1-infected mice by 4 months and together with MmuPV1 efficiently induced SCC by 6 months. UV radiation and exogenous estrogen cooperated to promote carcinogenesis in MmuPV1-infected mice. This murine infection model represents the first instance of de novo papillomavirus-mediated carcinogenesis in the female reproductive tract of wild-type mice resulting from active virus infection and is also the first report of the female hormone estrogen contributing to this process. This model will provide an additional platform for fundamental studies on papillomavirus infection, cervicovaginal disease, and the role of cellular cofactors during papillomavirus-induced carcinogenesis.IMPORTANCE Tractable and efficient models of papillomavirus-induced pathogenesis are limited due to the strict species-specific and tissue-specific tropism of these viruses. Here, we report a novel preclinical murine model of papillomavirus-induced cervicovaginal disease in wild-type, immunocompetent mice using the recently discovered murine papillomavirus, MmuPV1. In this model, MmuPV1 establishes persistent viral infections in the mucosal epithelia of the female reproductive tract, a necessary component needed to accurately mimic HPV-mediated neoplastic disease in humans. Persistent MmuPV1 infections were able to induce progressive neoplastic disease and carcinogenesis, either alone or in combination with previously identified cofactors of papillomavirus-induced disease. This new model will provide a much-needed platform for basic and translational studies on both papillomavirus infection and associated disease in immunocompetent mice.

The HPV16 E7 oncoprotein and 17β-estradiol are important factors for the induction of premalignant lesions and cervical cancer. The study of these factors is crucial for a better understanding of cervical tumorigenesis. Here, we assessed the global gene expression profiles induced by the HPV16 E7 oncoprotein and/or 17β-estradiol in cervical tissue of FvB and K14E7 transgenic mice. We found that the most dramatic changes in gene expression occurred in K14E7 and FvB groups treated with 17β-estradiol. A large number of differentially expressed genes involved in the immune response were observed in 17β-estradiol treated groups. The E7 oncoprotein mainly affected the expression of genes involved in cellular metabolism. Our microarray data also identified differentially expressed genes that have not previously been reported in cervical cancer. The identification of genes regulated by E7 and 17β-estradiol, provides the basis for further studies on their role in cervical carcinogenesis.

Head and neck squamous cell carcinoma (HNSCC) remains a challenging cancer to treat with overall 5-year survival on the order of 50-60%. Therefore, predictive biomarkers for this disease would be valuable to provide more effective and individualized therapeutic approaches for these patients. While prognostic biomarkers such as p16 expression correlate with outcome; to date, no predictive biomarkers have been clinically validated for HNSCC. We generated xenografts in immunocompromised mice from six established HNSCC cell lines and evaluated response to cisplatin, cetuximab, and radiation. Tissue microarrays were constructed from pre- and posttreatment tumor samples derived from each xenograft experiment. Quantitative immunohistochemistry was performed using a semiautomated imaging and analysis platform to determine the relative expression of five potential predictive biomarkers: epidermal growth factor receptor (EGFR), phospho-EGFR, phospho-Akt, phospho-ERK, and excision repair cross-complementation group 1 (ERCC1). Biomarker levels were compared between xenografts that were sensitive versus resistant to a specific therapy utilizing a two-sample t-test with equal standard deviations. Indeed the xenografts displayed heterogeneous responses to each treatment, and we linked a number of baseline biomarker levels to response. This included low ERCC1 being associated with cisplatin sensitivity, low phospho-Akt correlated with cetuximab sensitivity, and high total EGFR was related to radiation resistance. Overall, we developed a systematic approach to identifying predictive biomarkers and demonstrated several connections between biomarker levels and treatment response. Despite these promising initial results, this work requires additional preclinical validation, likely involving the use of patient-derived xenografts, prior to moving into the clinical realm for confirmation among patients with HNSCC.

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.

Autophagy is an intracellular catabolic process that removes and recycles unnecessary/dysfunctional cellular components, contributing to cellular health and survival. Autophagy is a highly regulated cellular process that responds to several intracellular signals, many of which are deregulated by human papillomavirus (HPV) infection through the expression of HPV-encoded oncoproteins. This adaptive inhibitory response helps prevent viral clearance. A strong correlation remains between HPV infection and the development of squamous cell carcinoma (SCC) of the anus, particularly in HIV positive and other immunosuppressed patients. We hypothesize that autophagy is inhibited by HPV-encoded oncoproteins thereby promoting anal carcinogenesis (Fig 1).

Epstein-Barr virus (EBV) encodes multiple miRNAs known to contribute to its pathogenicity. Previous studies have found that the levels of some EBV miRNAs are 10-100 fold higher in biopsies and in tumor xenografts than in cells grown in culture. We have asked if these increased levels reflect transcriptional enhancement resulting from the tumor microenvironment, selection for increased levels of the EBV genome, or both. We measured the levels of BART miRNAs and their DNA templates in tumor xenografts induced from EBV-positive gastric carcinoma cells and EBV-negative gastric carcinoma cells expressing plasmid replicons encoding these miRNAs. We focused on BART miRNAs which are expressed in all tumors and found that they provide tumors selective growth advantages as xenografts. Stem-loop PCR and real-time PCR revealed that the xenografts expressed both higher levels of some miRNAs and viral DNA templates than did the corresponding cells in culture.

The Hippo signaling pathway controls organ size and tumorigenesis through a kinase cascade that inactivates Yes-associated protein (YAP). Here, we show that YAP plays a central role in controlling the progression of cervical cancer. Our results suggest that YAP expression is associated with a poor prognosis for cervical cancer. TGF-α and amphiregulin (AREG), via EGFR, inhibit the Hippo signaling pathway and activate YAP to induce cervical cancer cell proliferation and migration. Activated YAP allows for up-regulation of TGF-α, AREG, and EGFR, forming a positive signaling loop to drive cervical cancer cell proliferation. HPV E6 protein, a major etiological molecule of cervical cancer, maintains high YAP protein levels in cervical cancer cells by preventing proteasome-dependent YAP degradation to drive cervical cancer cell proliferation. Results from human cervical cancer genomic databases and an accepted transgenic mouse model strongly support the clinical relevance of the discovered feed-forward signaling loop. Our study indicates that combined targeting of the Hippo and the ERBB signaling pathways represents a novel therapeutic strategy for prevention and treatment of cervical cancer.

High-risk types of human papillomavirus (HPV), such as HPV16, have been found in nearly all cases of cervical cancer. Therapies targeted at blocking the HPV16 E6 protein and its deleterious effects on the tumour suppressor pathways of the cell can reverse the malignant phenotype of affected keratinocytes while sparing uninfected cells. Through a strong interdisciplinary collaboration between engineering and biology, a novel, non-invasive intracellular delivery method for the HPV16 E6 antibody, F127-6G6, was developed. The method employs high intensity focused ultrasound (HIFU) in combination with microbubbles, in a process known as sonoporation. In this proof of principle study, it was first demonstrated that sonoporation antibody delivery into the HPV16 positive cervical carcinoma derived cell lines CaSki and SiHa was possible, using chemical transfection as a baseline for comparison. Delivery of the E6 antibody using sonoporation significantly restored p53 expression in these cells, indicating the antibody is able to enter the cells and remains active. This delivery method is targeted, non-cytotoxic, and non-invasive, making it more easily translatable for in vivo experiments than other transfection methods.

This work aimed to determine if cataractous changes associated with EMT occurring in the K14E6 mice lenses are associated with TGF-β and Wnt/β-catenin signaling activation.

The human papillomavirus DNA genome undergoes three distinct stages of replication: establishment, maintenance and amplification. We show that the HPV16 E6 protein is required for the maintenance of the HPV16 DNA genome as an extrachromosomal, nuclear plasmid in its natural host cell, the human keratinocyte. Based upon mutational analyses, inactivation of p53 by E6, but not necessarily E6-mediated degradation of p53, was found to correlate with the ability of E6 to support maintenance of the HPV16 genome as a nuclear plasmid. Inactivation of p53 with dominant negative p53 rescued the ability of HPV16 E6STOP and E6SAT mutant genomes to replicate as extrachromosomal genomes, though not to the same degree as observed for the HPV16 E6 wild-type (WT) genome. Inactivation of p53 also rescued the ability of HPV18 and HPV31 E6-deficient genomes to be maintained at copy numbers comparable to that of HPV18 and HPV31 E6WT genomes at early passages, though upon further passaging copy numbers for the HPV18 and 31 E6-deficient genomes lessened compared to that of the WT genomes. We conclude that inactivation of p53 is necessary for maintenance of HPV16 and for HPV18 and 31 to replicate at WT copy number, but that additional functions of E6 independent of inactivating p53 must also contribute to the maintenance of these genomes. Together these results suggest that re-activation of p53 may be a possible means for eradicating extrachromosomal HPV16, 18 or 31 genomes in the context of persistent infections.

Mouse papillomavirus type 1 (MmuPV1) provides, for the first time, the opportunity to study infection and pathogenesis of papillomaviruses in the context of laboratory mice. In this report, we define the transcriptome of MmuPV1 genome present in papillomas arising in experimentally infected mice using a combination of RNA-seq, PacBio Iso-seq, 5' RACE, 3' RACE, primer-walking RT-PCR, RNase protection, Northern blot and in situ hybridization analyses. We demonstrate that the MmuPV1 genome is transcribed unidirectionally from five major promoters (P) or transcription start sites (TSS) and polyadenylates its transcripts at two major polyadenylation (pA) sites. We designate the P7503, P360 and P859 as "early" promoters because they give rise to transcripts mostly utilizing the polyadenylation signal at nt 3844 and therefore can only encode early genes, and P7107 and P533 as "late" promoters because they give rise to transcripts utilizing polyadenylation signals at either nt 3844 or nt 7047, the latter being able to encode late, capsid proteins. MmuPV1 genome contains five splice donor sites and three acceptor sites that produce thirty-six RNA isoforms deduced to express seven predicted early gene products (E6, E7, E1, E1^M1, E1^M2, E2 and E8^E2) and three predicted late gene products (E1^E4, L2 and L1). The majority of the viral early transcripts are spliced once from nt 757 to 3139, while viral late transcripts, which are predicted to encode L1, are spliced twice, first from nt 7243 to either nt 3139 (P7107) or nt 757 to 3139 (P533) and second from nt 3431 to nt 5372. Thirteen of these viral transcripts were detectable by Northern blot analysis, with the P533-derived late E1^E4 transcripts being the most abundant. The late transcripts could be detected in highly differentiated keratinocytes of MmuPV1-infected tissues as early as ten days after MmuPV1 inoculation and correlated with detection of L1 protein and viral DNA amplification. In mature warts, detection of L1 was also found in more poorly differentiated cells, as previously reported. Subclinical infections were also observed. The comprehensive transcription map of MmuPV1 generated in this study provides further evidence that MmuPV1 is similar to high-risk cutaneous beta human papillomaviruses. The knowledge revealed will facilitate the use of MmuPV1 as an animal virus model for understanding of human papillomavirus gene expression, pathogenesis and immunology.

L83V-related variants of human papillomavirus (HPV) 16 E6, exemplified by the Asian-American variant Q14H/H78Y/L83V, were shown to be more prevalent than E6 prototype in progressing lesions and cervical cancer. We evaluated functions relevant to carcinogenesis for the E6 variants L83V, R10/L83V and Q14H/H78Y/L83V as well as the prototype in a model of human normal immortalized keratinocytes (NIKS). All E6 expressing NIKS equally abrogated growth arrest and DNA damage responses. Organotypic cultures derived from these keratinocytes demonstrated hyperplasia and aberrantly expressed keratin 5 in the suprabasal compartment. In contrast, differentiation and induction of apoptosis varied. The E6 variant rafts expressed keratin 10 in nearly all suprabasal cells while the prototype raft showed keratin 10 staining in a subset of suprabasal cells only. In addition, E6 variant NIKS expressing R10G/L83V and Q14H/H78Y/L83V were more prone to undergo cell-detachment-induced apoptosis (anoikis) than NIKS expressing E6 prototype. The combined differentiation and apoptosis pattern of high-risk E6 variants, especially of Q14H/H78Y/L83V, may reflect a phenotype beneficial to carcinogenesis and viral life cycle.

The DEK oncogene is overexpressed in many human malignancies including at early tumor stages. Our reported in vitro and in vivo models of squamous cell carcinoma have demonstrated that DEK contributes functionally to cellular and tumor survival and to proliferation. However, the underlying molecular mechanisms remain poorly understood. Based on recent RNA sequencing experiments, DEK expression was necessary for the transcription of several metabolic enzymes involved in anabolic pathways. This identified a possible mechanism whereby DEK may drive cellular metabolism to enable cell proliferation. Functional metabolic Seahorse analysis demonstrated increased baseline and maximum extracellular acidification rates, a readout of glycolysis, in DEK-overexpressing keratinocytes and squamous cell carcinoma cells. DEK overexpression also increased the maximum rate of oxygen consumption and therefore increased the potential for oxidative phosphorylation (OxPhos). To detect small metabolites that participate in glycolysis and the tricarboxylic acid cycle (TCA) that supplies substrate for OxPhos, we carried out NMR-based metabolomics studies. We found that high levels of DEK significantly reprogrammed cellular metabolism and altered the abundances of amino acids, TCA cycle intermediates and the glycolytic end products lactate, alanine and NAD+. Taken together, these data support a scenario whereby overexpression of the human DEK oncogene reprograms keratinocyte metabolism to fulfill energy and macromolecule demands required to enable and sustain cancer cell growth.

High-risk human papillomaviruses (HPVs) cause 5% of human cancers. Despite the availability of HPV vaccines, there remains a strong urgency to find ways to treat persistent HPV infections, as current HPV vaccines are not therapeutic for individuals already infected. We used a mouse papillomavirus infection model to characterize virus-host interactions. We found that mouse papillomavirus (MmuPV1) suppresses host immune responses via overexpression of stress keratins. In mice deficient for stress keratin K17 (K17KO), we observed rapid regression of papillomas dependent on T cells. Cellular genes involved in immune response were differentially expressed in the papillomas arising on the K17KO mice correlating with increased numbers of infiltrating CD8+ T cells and upregulation of IFNγ-related genes, including CXCL9 and CXCL10, prior to complete regression. Blocking the receptor for CXCL9/CXCL10 prevented early regression. Our data provide a novel mechanism by which papillomavirus-infected cells evade host immunity and defines new therapeutic targets for treating persistent papillomavirus infections.

The artemisinin family of compounds is cytopathic in certain cancer cell lines that are positive for human papillomaviruses (HPV) and can potentially drive the regression of dysplastic lesions. We evaluated the efficacy of topical dihydroartemisinin (DHA) on cervical dysplasia and anal dysplasia in two papillomavirus mouse models: K14E6/E7 transgenic mice, which express HPV16 oncogenes; and immunodeficient NOD/SCID gamma (NSG) mice infected with Mus musculus papillomavirus (MmuPV1). Mice started treatment with DHA at 25 weeks of age (K14E6/E7) or 20 weeks post infection (MmuPV1-infected), when the majority of mice are known to have papillomavirus-induced low- to high-grade dysplasia. Mice were treated with or without topical DHA at the cervix or anus and with or without topical treatment with the chemical carcinogen 7,12 dimethylbenz(a)anthracene (DMBA) at the anus of in transgenic mice to induce neoplastic progression. Mice were monitored for overt tumor growth, and tissue was harvested after 20 weeks of treatment and scored for severity of histological disease. For MmuPV1-infected mice, anogenital lavages were taken to monitor for viral clearance. Tissues were also evaluated for viral gene expression at the RNA and/or protein levels. Treatment with topical DHA did not reduce dysplasia in the anogenital tract in either papillomavirus-induced mouse model and did not prevent progression to anal cancer in the DMBA-treated K14E6/E7 mice.

Murine papillomavirus, MmuPV1, causes natural infections in laboratory mice that can progress to squamous cell carcinoma (SCC) making it a useful preclinical model to study the role of papillomaviruses in cancer. Papillomavirus can infect cells within hair follicles, which contain multiple epithelial progenitor cell populations, including Lgr5+ progenitors, and transgenic mice expressing human papillomavirus oncogenes develop tumors derived from Lgr5 progenitors. We therefore tested the hypothesis that Lgr5+ progenitors contribute to neoplastic lesions arising in skins infected with MmuPV1 by performing lineage tracing experiments. Ears of 6-8-week-old Lgr5-eGFP-IRES-CreERT2/Rosa26LSLtdTomato mice were treated topically with 4-OH Tamoxifen to label Lgr5+ progenitor cells and their progeny with tdTomato and, 72 h later, infected with MmuPV1. Four months post-infection, tissue at the infection site was harvested for histopathological analysis and immunofluorescence to determine the percentage of tdTomato+ cells within the epithelial lesions caused by MmuPV1. Squamous cell dysplasia showed a low percentage of tdTomato+ cells (7%), indicating that it arises primarily from non-Lgr5 progenitor cells. In contrast, cutaneous SCC (cSCC) was substantially more positive for tdTomato+ cells (42%), indicating that cSCCs preferentially arise from Lgr5+ progenitors. Biomarker analyses of dysplasia vs. cSCC revealed further differences consistent with cSCC arising from LGR5+ progenitor cells.

The species specificity of papillomaviruses has been a significant roadblock for performing in vivo pathogenesis studies in common model organisms. The Mus musculus papillomavirus type 1 (MmuPV1) causes cutaneous papillomas that can progress to squamous cell carcinomas in laboratory mice. The papillomavirus E6 and E7 genes encode proteins that establish and maintain a cellular milieu that allows for viral genome synthesis and viral progeny synthesis in growth-arrested, terminally differentiated keratinocytes. The E6 and E7 proteins provide this activity by binding to and functionally reprogramming key cellular regulatory proteins. The MmuPV1 E7 protein lacks the canonical LXCXE motif that mediates the binding of multiple viral oncoproteins to the cellular retinoblastoma tumor suppressor protein, RB1. Our proteomic experiments, however, revealed that MmuPV1 E7 still interacts with RB1. We show that MmuPV1 E7 interacts through its C terminus with the C-terminal domain of RB1. Binding of MmuPV1 E7 to RB1 did not cause significant activation of E2F-regulated cellular genes. MmuPV1 E7 expression was shown to be essential for papilloma formation. Experimental infection of mice with MmuPV1 expressing an E7 mutant that is defective for binding to RB1 caused delayed onset, lower incidence, and smaller sizes of papillomas. Our results demonstrate that the MmuPV1 E7 gene is essential and that targeting noncanonical activities of RB1, which are independent of RB1's ability to modulate the expression of E2F-regulated genes, contribute to papillomavirus-mediated pathogenesis. IMPORTANCE Papillomavirus infections cause a variety of epithelial hyperplastic lesions, or warts. While most warts are benign, some papillomaviruses cause lesions that can progress to squamous cell carcinomas, and approximately 5% of all human cancers are caused by human papillomavirus (HPV) infections. The papillomavirus E6 and E7 proteins are thought to function to reprogram host epithelial cells to enable viral genome replication in terminally differentiated, normally growth-arrested cells. E6 and E7 lack enzymatic activities and function by interacting and functionally altering host cell regulatory proteins. Many cellular proteins that can interact with E6 and E7 have been identified, but the biological relevance of these interactions for viral pathogenesis has not been determined. This is because papillomaviruses are species specific and do not infect heterologous hosts. Here, we use a recently established mouse papillomavirus (MmuPV1) model to investigate the role of the E7 protein in viral pathogenesis. We show that MmuPV1 E7 is necessary for papilloma formation. The retinoblastoma tumor suppressor protein (RB1) is targeted by many papillomaviral E7 proteins, including cancer-associated HPVs. We show that MmuPV1 E7 can bind RB1 and that infection with a mutant MmuPV1 virus that expresses an RB1 binding-defective E7 mutant caused smaller and fewer papillomas that arise with delayed kinetics.