Stress keratin 17 and estrogen support viral persistence and modulate the immune environment during cervicovaginal murine papillomavirus infection.

A murine papillomavirus, MmuPV1, infects both cutaneous and mucosal epithelia of laboratory mice and can be used to model high-risk human papillomavirus (HPV) infection and HPV-associated disease. We have shown that estrogen exacerbates papillomavirus-induced cervical disease in HPV-transgenic mice. We have also previously identified stress keratin 17 (K17) as a host factor that supports MmuPV1-induced cutaneous disease. Here, we sought to test the role of estrogen and K17 in MmuPV1 infection and associated disease in the female reproductive tract. We experimentally infected wild-type and K17 knockout (K17KO) mice with MmuPV1 in the female reproductive tract in the presence or absence of exogenous estrogen for 6 mon. We observed that a significantly higher percentage of K17KO mice cleared the virus as opposed to wild-type mice. In estrogen-treated wild-type mice, the MmuPV1 viral copy number was significantly higher compared to untreated mice by as early as 2 wk postinfection, suggesting that estrogen may help facilitate MmuPV1 infection and/or establishment. Consistent with this, viral clearance was not observed in either wild-type or K17KO mice when treated with estrogen. Furthermore, neoplastic disease progression and cervical carcinogenesis were supported by the presence of K17 and exacerbated by estrogen treatment. Subsequent analyses indicated that estrogen treatment induces a systemic immunosuppressive state in MmuPV1-infected animals and that both estrogen and K17 modulate the local intratumoral immune microenvironment within MmuPV1-induced neoplastic lesions. Collectively, these findings suggest that estrogen and K17 act at multiple stages of papillomavirus-induced disease at least in part via immunomodulatory mechanisms.

HPV16 E6 induces chromosomal instability due to polar chromosomes caused by E6AP-dependent degradation of the mitotic kinesin CENP-E.

Chromosome segregation during mitosis is highly regulated to ensure production of genetically identical progeny. Recurrent mitotic errors cause chromosomal instability (CIN), a hallmark of tumors. The E6 and E7 oncoproteins of high-risk human papillomavirus (HPV), which causes cervical, anal, and head and neck cancers (HNC), cause mitotic defects consistent with CIN in models of anogenital cancers, but this has not been studied in the context of HNC. Here, we show that HPV16 induces a specific type of CIN in patient HNC tumors, patient-derived xenografts, and cell lines, which is due to defects in chromosome congression. These defects are specifically induced by the HPV16 oncogene E6 rather than E7. We show that HPV16 E6 expression causes degradation of the mitotic kinesin CENP-E, whose depletion produces chromosomes that are chronically misaligned near spindle poles (polar chromosomes) and fail to congress. Though the canonical oncogenic role of E6 is the degradation of the tumor suppressor p53, CENP-E degradation and polar chromosomes occur independently of p53. Instead, E6 directs CENP-E degradation in a proteasome-dependent manner via the E6-associated ubiquitin protein ligase E6AP/UBE3A. This study reveals a mechanism by which HPV induces CIN, which may impact HPV-mediated tumor initiation, progression, and therapeutic response.

Deficiency in Ever2 does not increase susceptibility of mice to pathogenesis by the mouse papillomavirus, MmuPV1.

Epidermodysplasia verruciformis (EV) is a rare genetic skin disorder that is characterized by the development of papillomavirus-induced skin lesions that can progress to squamous cell carcinoma (SCC). Certain high-risk, cutaneous ß-genus human papillomaviruses (ß-HPVs), in particular HPV5 and HPV8, are associated with inducing EV in individuals who have a homozygous mutation in one of three genes tied to this disease: EVER1, EVER2, or CIB1. EVER1 and EVER2 are also known as TMC6 and TMC8, respectively. Little is known about the biochemical activities of EVER gene products or their roles in facilitating EV in conjunction with ß-HPV infection. To investigate the potential effect of EVER genes on papillomavirus infection, we pursued in vivo infection studies by infecting Ever2-null mice with mouse papillomavirus (MmuPV1). MmuPV1 shares characteristics with ß-HPVs including similar genome organization, shared molecular activities of their early, E6 and E7, oncoproteins, the lack of a viral E5 gene, and the capacity to cause skin lesions that can progress to SCC. MmuPV1 infections were conducted both in the presence and absence of UVB irradiation, which is known to increase the risk of MmuPV1-induced pathogenesis. Infection with MmuPV1 induced skin lesions in both wild-type and Ever2-null mice with and without UVB. Many lesions in both genotypes progressed to malignancy, and the disease severity did not differ between Ever2-null and wild-type mice. However, somewhat surprisingly, lesion growth and viral transcription was decreased, and lesion regression was increased in Ever2-null mice compared with wild-type mice. These studies demonstrate that Ever2-null mice infected with MmuPV1 do not exhibit the same phenotype as human EV patients infected with ß-HPVs.IMPORTANCEHumans with homozygous mutations in the EVER2 gene develop epidermodysplasia verruciformis (EV), a disease characterized by predisposition to persistent ß-genus human papillomavirus (ß-HPV) skin infections, which can progress to skin cancer. To investigate how EVER2 confers protection from papillomaviruses, we infected the skin of homozygous Ever2-null mice with mouse papillomavirus MmuPV1. Like in humans with EV, infected Ever2-null mice developed skin lesions that could progress to cancer. Unlike in humans with EV, lesions in these Ever2-null mice grew more slowly and regressed more frequently than in wild-type mice. MmuPV1 transcription was higher in wild-type mice than in Ever2-null mice, indicating that mouse EVER2 does not confer protection from papillomaviruses. These findings suggest that there are functional differences between MmuPV1 and ß-HPVs and/or between mouse and human EVER2.

MmuPV1 E7's interaction with PTPN14 delays Epithelial differentiation and contributes to virus-induced skin disease.

Human papillomaviruses (HPVs) contribute to approximately 5% of all human cancers. Species-specific barriers limit the ability to study HPV pathogenesis in animal models. Murine papillomavirus (MmuPV1) provides a powerful tool to study the roles of papillomavirus genes in pathogenesis arising from a natural infection. We previously identified Protein Tyrosine Phosphatase Non-Receptor Type 14 (PTPN14), a tumor suppressor targeted by HPV E7 proteins, as a putative cellular target of MmuPV1 E7. Here, we confirmed the MmuPV1 E7-PTPN14 interaction. Based on the published structure of the HPV18 E7/PTPN14 complex, we generated a MmuPV1 E7 mutant, E7K81S, that was defective for binding PTPN14. Wild-type (WT) and E7K81S mutant viral genomes replicated as extrachromosomal circular DNAs to comparable levels in mouse keratinocytes. E7K81S mutant virus (E7K81S MmuPV1) was generated and used to infect FoxN/Nude mice. E7K81S MmuPV1 caused neoplastic lesions at a frequency similar to that of WT MmuPV1, but the lesions arose later and were smaller than WT-induced lesions. The E7K81S MmuPV1-induced lesions also had a trend towards a less severe grade of neoplastic disease. In the lesions, E7K81S MmuPV1 supported the late (productive) stage of the viral life cycle and promoted E2F activity and cellular DNA synthesis in suprabasal epithelial cells to similar degrees as WT MmuPV1. There was a similar frequency of lateral spread of infections among mice infected with E7K81S or WT MmuPV1. Compared to WT MmuPV1-induced lesions, E7K81S MmuPV1-induced lesions had a significant expansion of cells expressing differentiation markers, Keratin 10 and Involucrin. We conclude that an intact PTPN14 binding site is necessary for MmuPV1 E7's ability to contribute to papillomavirus-induced pathogenesis and this correlates with MmuPV1 E7 causing a delay in epithelial differentiation, which is a hallmark of papillomavirus-induced neoplasia.

Merkel cell polyomavirus large T antigen binding to pRb promotes skin hyperplasia and tumor development.

Clear evidence supports a causal link between Merkel cell polyomavirus (MCPyV) and the highly aggressive human skin cancer called Merkel cell carcinoma (MCC). Integration of viral DNA into the human genome facilitates continued expression of the MCPyV small tumor (ST) and large tumor (LT) antigens in virus-positive MCCs. In MCC tumors, MCPyV LT is truncated in a manner that renders the virus unable to replicate yet preserves the LXCXE motif that facilitates its binding to and inactivation of the retinoblastoma tumor suppressor protein (pRb). We previously developed a MCPyV transgenic mouse model in which MCC tumor-derived ST and truncated LT expression were targeted to the stratified epithelium of the skin, causing epithelial hyperplasia, increased proliferation, and spontaneous tumorigenesis. We sought to determine if any of these phenotypes required the association between the truncated MCPyV LT and pRb. Mice were generated in which K14-driven MCPyV ST/LT were expressed in the context of a homozygous RbΔLXCXE knock-in allele that attenuates LT-pRb interactions through LT's LXCXE motif. We found that many of the phenotypes including tumorigenesis that develop in the K14-driven MCPyV transgenic mice were dependent upon LT's LXCXE-dependent interaction with pRb. These findings highlight the importance of the MCPyV LT-pRb interaction in an in vivo model for MCPyV-induced tumorigenesis.

Type 1 and Type 2 Epstein-Barr viruses induce proliferation, and inhibit differentiation, in infected telomerase-immortalized normal oral keratinocytes.

Differentiated epithelial cells are an important source of infectious EBV virions in human saliva, and latent Epstein-Barr virus (EBV) infection is strongly associated with the epithelial cell tumor, nasopharyngeal carcinoma (NPC). However, it has been difficult to model how EBV contributes to NPC, since EBV has not been shown to enhance proliferation of epithelial cells in monolayer culture in vitro and is not stably maintained in epithelial cells without antibiotic selection. In addition, although there are two major types of EBV (type 1 (T1) and type 2 (T2)), it is currently unknown whether T1 and T2 EBV behave differently in epithelial cells. Here we inserted a G418 resistance gene into the T2 EBV strain, AG876, allowing us to compare the phenotypes of T1 Akata virus versus T2 AG876 virus in a telomerase-immortalized normal oral keratinocyte cell line (NOKs) using a variety of different methods, including RNA-seq analysis, proliferation assays, immunoblot analyses, and air-liquid interface culture. We show that both T1 Akata virus infection and T2 AG876 virus infection of NOKs induce cellular proliferation, and inhibit spontaneous differentiation, in comparison to the uninfected cells when cells are grown without supplemental growth factors in monolayer culture. T1 EBV and T2 EBV also have a similar ability to induce epithelial-to-mesenchymal (EMT) transition and activate canonical and non-canonical NF-κB signaling in infected NOKs. In contrast to our recent results in EBV-infected lymphoblastoid cells (in which T2 EBV infection is much more lytic than T1 EBV infection), we find that NOKs infected with T1 and T2 EBV respond similarly to lytic inducing agents such as TPA treatment or differentiation. These results suggest that T1 and T2 EBV have similar phenotypes in infected epithelial cells, with both EBV types enhancing cellular proliferation and inhibiting differentiation when growth factors are limiting.

Mouse papillomavirus type 1 (MmuPV1) DNA is frequently integrated in benign tumors by microhomology-mediated end-joining.

MmuPV1 is a useful model for studying papillomavirus-induced tumorigenesis. We used RNA-seq to look for chimeric RNAs that map to both MmuPV1 and host genomes. In tumor tissues, a higher proportion of total viral reads were virus-host chimeric junction reads (CJRs) (1.9‰ - 7‰) than in tumor-free tissues (0.6‰ - 1.3‰): most CJRs mapped to the viral E2/E4 region. Although most of the MmuPV1 integration sites were mapped to intergenic regions and introns throughout the mouse genome, integrations were seen more than once in several genes: Malat1, Krt1, Krt10, Fabp5, Pard3, and Grip1; these data were confirmed by rapid amplification of cDNA ends (RACE)-Single Molecule Real-Time (SMRT)-seq or targeted DNA-seq. Microhomology sequences were frequently seen at host-virus DNA junctions. MmuPV1 infection and integration affected the expression of host genes. We found that factors for DNA double-stranded break repair and microhomology-mediated end-joining (MMEJ), such as H2ax, Fen1, DNA polymerase Polθ, Cdk1, and Plk1, exhibited a step-wise increase and Mdc1 a decrease in expression in MmuPV1-infected tissues and MmuPV1 tumors relative to normal tissues. Increased expression of mitotic kinases CDK1 and PLK1 appears to be correlated with CtIP phosphorylation in MmuPV1 tumors, suggesting a role for MMEJ-mediated DNA joining in the MmuPV1 integration events that are associated with MmuPV1-induced progression of tumors.

ΔNp63α promotes Epstein-Barr virus latency in undifferentiated epithelial cells.

Epstein-Barr virus (EBV) is a human herpesvirus that causes infectious mononucleosis and contributes to both B-cell and epithelial-cell malignancies. EBV-infected epithelial cell tumors, including nasopharyngeal carcinoma (NPC), are largely composed of latently infected cells, but the mechanism(s) maintaining viral latency are poorly understood. Expression of the EBV BZLF1 (Z) and BRLF1 (R) encoded immediate-early (IE) proteins induces lytic infection, and these IE proteins activate each other's promoters. ΔNp63α (a p53 family member) is required for proliferation and survival of basal epithelial cells and is over-expressed in NPC tumors. Here we show that ΔNp63α promotes EBV latency by inhibiting activation of the BZLF1 IE promoter (Zp). Furthermore, we find that another p63 gene splice variant, TAp63α, which is expressed in some Burkitt and diffuse large B cell lymphomas, also represses EBV lytic reactivation. We demonstrate that ΔNp63α inhibits the Z promoter indirectly by preventing the ability of other transcription factors, including the viral IE R protein and the cellular KLF4 protein, to activate Zp. Mechanistically, we show that ΔNp63α promotes viral latency in undifferentiated epithelial cells both by enhancing expression of a known Zp repressor protein, c-myc, and by decreasing cellular p38 kinase activity. Furthermore, we find that the ability of cis-platinum chemotherapy to degrade ΔNp63α contributes to the lytic-inducing effect of this agent in EBV-infected epithelial cells. Together these findings demonstrate that the loss of ΔNp63α expression, in conjunction with enhanced expression of differentiation-dependent transcription factors such as BLIMP1 and KLF4, induces lytic EBV reactivation during normal epithelial cell differentiation. Conversely, expression of ΔNp63α in undifferentiated nasopharyngeal carcinoma cells and TAp63α in Burkitt lymphoma promotes EBV latency in these malignancies.

Cell-penetrating peptide inhibits retromer-mediated human papillomavirus trafficking during virus entry.

Virus replication requires critical interactions between viral proteins and cellular proteins that mediate many aspects of infection, including the transport of viral genomes to the site of replication. In human papillomavirus (HPV) infection, the cellular protein complex known as retromer binds to the L2 capsid protein and sorts incoming virions into the retrograde transport pathway for trafficking to the nucleus. Here, we show that short synthetic peptides containing the HPV16 L2 retromer-binding site and a cell-penetrating sequence enter cells, sequester retromer from the incoming HPV pseudovirus, and inhibit HPV exit from the endosome, resulting in loss of viral components from cells and in a profound, dose-dependent block to infection. The peptide also inhibits cervicovaginal HPV16 pseudovirus infection in a mouse model. These results confirm the retromer-mediated model of retrograde HPV entry and validate intracellular virus trafficking as an antiviral target. More generally, inhibiting virus replication with agents that can enter cells and disrupt essential protein-protein interactions may be applicable in broad outline to many viruses.

IQGAP1 and RNA Splicing in the Context of Head and Neck via Phosphoproteomics.

IQGAP1 (IQ motif-containing GTPase-activating protein 1) scaffolds several signaling pathways in mammalian cells that are implicated in carcinogenesis, including the RAS and PI3K pathways that involve multiple protein kinases. IQGAP1 has been shown to promote head and neck squamous cell carcinoma (HNSCC); however, the underlying mechanism(s) remains unclear. Here, we report a mass spectrometry-based analysis identifying differences in phosphorylation of cellular proteins in vivo and in vitro in the presence or absence of IQGAP1. By comparing the esophageal phosphoproteome profiles between Iqgap1+/+ and Iqgap1-/- mice, we identified RNA splicing as one of the most altered cellular processes. Serine/arginine-rich splicing factor 6 (SRSF6) was the protein with the most downregulated levels of phosphorylation in Iqgap1-/- tissue. We confirmed that the absence of IQGAP1 reduced SRSF6 phosphorylation both in vivo and in vitro. We then expanded our analysis to human normal oral keratinocytes. Again, we found factors involved in RNA splicing to be highly altered in the phosphoproteome profile upon genetic disruption of IQGAP1. Both the Clinical Proteomic Tumor Analysis Consortium (CPTAC) and the Cancer Genome Atlas (TCGA) data sets indicate that phosphorylation of splicing-related proteins is important in HNSCC prognosis. The Biological General Repository for Interaction Datasets (BioGRID) repository also suggested multiple interactions between IQGAP1 and splicing-related proteins. Based on these collective observations, we propose that IQGAP1 regulates the phosphorylation of splicing proteins, which potentially affects their splicing activities and, therefore, contributes to HNSCC. Raw data are available from the MassIVE database with identifier MSV000087770.

Stress keratin 17 enhances papillomavirus infection-induced disease by downregulating T cell recruitment.

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.

Mus musculus Papillomavirus 1: a New Frontier in Animal Models of Papillomavirus Pathogenesis.

Animal models of viral pathogenesis are essential tools in human disease research. Human papillomaviruses (HPVs) are a significant public health issue due to their widespread sexual transmission and oncogenic potential. Infection-based models of papillomavirus pathogenesis have been complicated by their strict species and tissue specificity. In this Gem, we discuss the discovery of a murine papillomavirus, Mus musculus papillomavirus 1 (MmuPV1), and how its experimental use represents a major advancement in models of papillomavirus-induced pathogenesis/carcinogenesis, and their transmission.

Reactivation of Epstein-Barr Virus by HIF-1α Requires p53.

We previously reported that the cellular transcription factor hypoxia-inducible factor 1α (HIF-1α) binds a hypoxia response element (HRE) located within the promoter of Epstein-Barr virus's (EBV's) latent-lytic switch BZLF1 gene, Zp, inducing viral reactivation. In this study, EBV-infected cell lines derived from gastric cancers and Burkitt lymphomas were incubated with HIF-1α-stabilizing drugs: the iron chelator deferoxamine (Desferal [DFO]), a neddylation inhibitor (pevonedistat [MLN-4924]), and a prolyl hydroxylase inhibitor (roxadustat [FG-4592]). DFO and MLN-4924, but not FG-4592, induced accumulation of both lytic EBV proteins and phosphorylated p53 in cell lines that contain a wild-type p53 gene. FG-4592 also failed to activate transcription from Zp in a reporter assay despite inducing accumulation of HIF-1α and transcription from another HRE-containing promoter. Unexpectedly, DFO failed to induce EBV reactivation in cell lines that express mutant or no p53 or when p53 expression was knocked down with short hairpin RNAs (shRNAs). Likewise, HIF-1α failed to activate transcription from Zp when p53 was knocked out by CRISPR-Cas9. Importantly, DFO induced binding of p53 as well as HIF-1α to Zp in chromatin immunoprecipitation (ChIP) assays, but only when the HRE was present. Nutlin-3, a drug known to induce accumulation of phosphorylated p53, synergized with DFO and MLN-4924 in inducing EBV reactivation. Conversely, KU-55933, a drug that inhibits ataxia telangiectasia mutated, thereby preventing p53 phosphorylation, inhibited DFO-induced EBV reactivation. Lastly, activation of Zp transcription by DFO and MLN-4924 mapped to its HRE. Thus, we conclude that induction of BZLF1 gene expression by HIF-1α requires phosphorylated, wild-type p53 as a coactivator, with HIF-1α binding recruiting p53 to Zp.IMPORTANCE EBV, a human herpesvirus, is latently present in most nasopharyngeal carcinomas, Burkitt lymphomas, and some gastric cancers. To develop a lytic-induction therapy for treating patients with EBV-associated cancers, we need a way to efficiently reactivate EBV into lytic replication. EBV's BZLF1 gene product, Zta, usually controls this reactivation switch. We previously showed that HIF-1α binds the BZLF1 gene promoter, inducing Zta synthesis, and HIF-1α-stabilizing drugs can induce EBV reactivation. In this study, we determined which EBV-positive cell lines are reactivated by classes of HIF-1α-stabilizing drugs. We found, unexpectedly, that HIF-1α-stabilizing drugs only induce reactivation when they also induce accumulation of phosphorylated, wild-type p53. Fortunately, p53 phosphorylation can also be provided by drugs such as nutlin-3, leading to synergistic reactivation of EBV. These findings indicate that some HIF-1α-stabilizing drugs may be helpful as part of a lytic-induction therapy for treating patients with EBV-positive malignancies that contain wild-type p53.

The vaginal and fecal microbiota of a murine cervical carcinoma model under synergistic effect of 17ß-Estradiol and E7 oncogene expression.

Cervical cancer is an important health issue worldwide. Many factors are related to this condition as the persistence of human papillomavirus (HPV) infection (e.g. type 16 and 18), the use of hormonal contraceptives for long periods of time, pH changes and bacterial vaginosis. The association between the microbiota and cervical human cancer is an interesting issue to be explored; given that environmental and hormonal factors may change the vaginal microbiota contributing to this condition. Our hypothesis was that changes in the microbiota diversity is associated with the development of cervical cancer. We evaluated the microbiota diversity in vaginal lavages and fecal samples at different stages of cervical cancer development in a mice model (K14HPV16E7) with type 16 E7 oncogene expression (E7), under continuous or not continuous stimulus of 17ß-estradiol (E2) and compared it with a non-transgenic isogenic control (FVB) under same conditions. Our results indicate that continuous E2 administration during 6 months in the model with type 16 E7 expression causing development of cancer, is associated with significant changes in the microbiota diversity of the cervicovaginal lavages. Similar results were not observed in the same model when no E2 was administered to the mice. The FVB mice with no E7 expression which do not develop cervical cancer, did not show comparable changes in the microbiota diversity when E2 was administered during the same period. Normal evolution of the cervical epithelium and microbiota diversity were observed for the FVB mice with no E2 administration. Large changes in the microbiota diversity in fecal samples were not observed suggesting a specific organ effect of E7 expression associated to E2 on the vaginal microbiota.

Directed differentiation of human pluripotent stem cells into epidermal stem and progenitor cells.

BACKGROUND: Pluripotent stem cells (PSCs) produced by somatic cell reprogramming self-renew in culture and can differentiate into any cell type, representing a powerful tool for disease modeling, drug screening, regenerative medicine, and the discovery of personalized therapies to treat tissue-specific pathologies. We previously reported the directed differentiation of human PSCs into epidermal stem and progenitor cells (ESPCs) and 3D epidermis to model the inherited syndrome Fanconi anemia (FA), wherein epidermal cell-junctional defects discovered using this system were validated in patient populations. Here, we describe in detail the corresponding protocol for generating PSC-derived keratinocytes using a distinct, normal PSC line (209.2 PSC). METHODS AND RESULTS: Our approach modifies previous protocols to minimize spontaneous cell death and terminal differentiation, eliminate cell stress-inducing keratinocyte selection steps, and reduce total protocol duration and cost. Independent donor-derived PSC lines were converted into ESPCs through the addition of relevant morphogens and a ROCK inhibitor. Results for the 209.2 PSC line highlight consistencies in 2D and also variable features in 3D epidermis compared to the previously published FA-PSC lines. 209.2 PSC-derived ESPCs exhibited a basal cell phenotype while maintaining the capacity to form epidermal organotypic rafts with morphology consistent with fetal epidermis. Transcriptional analyses demonstrated 209.2 ESPCs express epidermis-selective markers and not early endoderm markers, thus supporting an immature stage of p63+ epidermal development. CONCLUSIONS: This protocol provides an accelerated path for the generation of human ESPCs and 3D epidermal models to study normal epidermal development and homeostasis, elucidate mechanisms of epidermal disease pathogenesis, and provides a platform for developing personalized therapies.

Cooperation of genes in HPV16 E6/E7-dependent cervicovaginal carcinogenesis trackable by endoscopy and independent of exogenous estrogens or carcinogens.

Human papillomavirus (HPV) infection is necessary but insufficient for progression of epithelial cells from dysplasia to carcinoma-in situ (CIS) to invasive cancer. The combination of mutant cellular and viral oncogenes that regulate progression of cervical cancer (CC) remains unclear. Using combinations of HPV16 E6/E7 (E+), mutant Kras (mKras) (K+) and/or loss of Pten (P-/-), we generated autochthonous models of CC without exogenous estrogen, carcinogen or promoters. Furthermore, intravaginal instillation of adenoCre virus enabled focal activation of the oncogenes/inactivation of the tumor suppressor gene. In P+/+ mice, E6/E7 alone (P+/+E+K-) failed to cause premalignant changes, while mKras alone (P+/+E-K+) caused persistent mucosal abnormalities in about one-third of mice, but no cancers. To develop cancer, P+/+ mice needed both E6/E7 and mKras expression. Longitudinal endoscopies of P+/+E+K+ mice predicted carcinoma development by detection of mucosal lesions, found on an average of 23 weeks prior to death, unlike longitudinal quantitative PCRs of vaginal lavage samples from the same mice. Endoscopy revealed that individual mice differed widely in the time required for mucosal lesions to appear after adenoCre and in the time required for these lesions to progress to cancer. These cancers developed in the transition zone that extends, unlike in women, from the murine cervix to the distal vagina. The P-/-E+K+ genotype led to precipitous cancer development within a few weeks and E6/E7-independent cancer development occurred in the P-/-E-K+ genotype. In the P-/-E+K- genotype, mice only developed CIS. Thus, distinct combinations of viral and cellular oncogenes are involved in distinct steps in cervical carcinogenesis.

Human papillomavirus oncogenes reprogram the cervical cancer microenvironment independently of and synergistically with estrogen.

High-risk human papillomaviruses (HPVs) infect epithelial cells and are causally associated with cervical cancer, but HPV infection is not sufficient for carcinogenesis. Previously, we reported that estrogen signaling in the stromal tumor microenvironment is associated with cervical cancer maintenance and progression. We have now determined how HPV oncogenes and estrogen treatment affect genome-wide host gene expression in laser-captured regions of the cervical epithelium and stroma of untreated or estrogen-treated nontransgenic and HPV-transgenic mice. HPV oncogene expression in the cervical epithelium elicited significant gene-expression changes in the proximal stromal compartment, and estrogen treatment uniquely affected gene expression in the cervical microenvironment of HPV-transgenic mice compared with nontransgenic mice. Several potential estrogen-induced paracrine-acting factors were identified in the expression profile of the cervical tumor microenvironment. The microenvironment of estrogen-treated HPV-transgenic mice was significantly enriched for chemokine/cytokine activity and inflammatory and immune functions associated with carcinogenesis. This inflammatory signature included several proangiogenic CXCR2 receptor ligands. A subset of the same CXCR2 ligands was likewise increased in cocultures of early-passage cells from human cervical samples, with levels highest in cocultures of cervical fibroblasts and cancer-derived epithelial cells. Our studies demonstrate that high-risk HPV oncogenes profoundly reprogram the tumor microenvironment independently of and synergistically with estrogen. These observations illuminate important means by which HPVs can cause cancer through alterations in the tumor microenvironment.

Cutaneous HPV8 and MmuPV1 E6 Proteins Target the NOTCH and TGF-ß Tumor Suppressors to Inhibit Differentiation and Sustain Keratinocyte Proliferation.

Cutaneous beta-papillomaviruses are associated with non-melanoma skin cancers that arise in patients who suffer from a rare genetic disorder, Epidermodysplasia verruciformis (EV) or after immunosuppression following organ transplantation. Recent studies have shown that the E6 proteins of the cancer associated beta human papillomavirus (HPV) 5 and HPV8 inhibit NOTCH and TGF-ß signaling. However, it is unclear whether disruption of these pathways may contribute to cutaneous HPV pathogenesis and carcinogenesis. A recently identified papillomavirus, MmuPV1, infects laboratory mouse strains and causes cutaneous skin warts that can progress to squamous cell carcinoma. To determine whether MmuPV1 may be an appropriate model to mechanistically dissect the molecular contributions of cutaneous HPV infections to skin carcinogenesis, we investigated whether MmuPV1 E6 shares biological and biochemical activities with HPV8 E6. We report that the HPV8 and MmuPV1 E6 proteins share the ability to bind to the MAML1 and SMAD2/SMAD3 transcriptional cofactors of NOTCH and TGF-beta signaling, respectively. Moreover, we demonstrate that these cutaneous papillomavirus E6 proteins inhibit these two tumor suppressor pathways and that this ability is linked to delayed differentiation and sustained proliferation of differentiating keratinocytes. Furthermore, we demonstrate that the ability of MmuPV1 E6 to bind MAML1 is necessary for papilloma formation in experimentally infected mice. Our results, therefore, suggest that experimental MmuPV1 infection in mice will be a robust and useful experimental system to model key aspects of cutaneous HPV infection, pathogenesis and carcinogenesis.

The full transcription map of mouse papillomavirus type 1 (MmuPV1) in mouse wart tissues.

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.