Human papillomavirus (HPV) is the most common sexually transmitted pathogen, and high-risk HPVs contribute to 5% of human cancers, including 25% of head and neck squamous cell carcinomas (HNSCCs). Despite the significant role played by HPVs in HNSCC, there is currently no available in vivo system to model the process from papillomavirus infection to virus-induced HNSCC. In this paper, we describe an infection-based HNSCC model, utilizing a mouse papillomavirus (MmuPV1), which naturally infects laboratory mice. Infections of the tongue epithelium of two immunodeficient strains with MmuPV1 caused high-grade squamous dysplasia with early signs of invasive carcinoma over the course of 4 months. When combined with the oral carcinogen 4-nitroquinoline-1-oxide (4NQO), MmuPV1 caused invasive squamous cell carcinoma (SCC) on the tongue of both immunodeficient and immunocompetent mice. These tumors expressed markers of papillomavirus infection and HPV-associated carcinogenesis. This novel preclinical model provides a valuable new means to study how natural papillomavirus infections contribute to HNSCC.IMPORTANCE The species specificity of papillomavirus has limited the development of an infection-based animal model to study HPV-associated head and neck carcinogenesis. Our study presents a novel in vivo model using the mouse papillomavirus MmuPV1 to study papillomavirus-associated head and neck cancer. In our model, MmuPV1 infects and causes lesions in both immunodeficient and genetically immunocompetent strains of mice. These virally induced lesions carry features associated with both HPV infections and HPV-associated carcinogenesis. Combined with previously identified cancer cofactors, MmuPV1 causes invasive squamous cell carcinomas in mice. This model provides opportunities for basic and translational studies of papillomavirus infection-based head and neck disease.

The E6 protein encoded by the murine papillomavirus (MmuPV1) is essential for MmuPV1-induced skin disease. Our previous work has identified a number of cellular interacting partners of MmuPV1 E6 and E7 through affinity purification/mass spectrometry analysis. These studies revealed that MmuPV1 E6 potently inhibits keratinocyte differentiation through multiple molecular mechanisms including inhibition of NOTCH and TGF-β signaling. Here, we report that MmuPV1 E6 has additional important oncogenic activities when expressed in its natural host cells, mouse keratinocytes, including increasing proliferation, overcoming density-mediated growth arrest, and proliferation under conditions of limited supply of growth factors. Unbiased proteomic/transcriptomic analyses of mouse keratinocytes expressing MmuPV1 E6 substantiated its effect on these cellular processes and divulged that some of these effects may be mediated in part through it upregulating E2F activity. Our analyses also revealed that MmuPV1 E6 may alter other cancer hallmarks including evasion of growth suppressors, inhibition of immune response, resistance to cell death, and alterations in DNA damage response. Collectively, our results suggest that MmuPV1 E6 is a major driver of multiple hallmarks of cancer in MmuPV1's natural host cells, mouse keratinocytes.IMPORTANCEThe Mus musculus papillomavirus 1 (MmuPV1) E6 and E7 proteins are required for MmuPV1-induced disease. Our understanding of the activities of MmuPV1 E6 has been based on affinity purification/mass spectrometry studies where cellular interacting partners of MmuPV1 E6 were identified, and these studies revealed that MmuPV1 E6 can inhibit keratinocyte differentiation through multiple mechanisms. We report that MmuPV1 E6 encodes additional activities including the induction of proliferation, resistance to density-mediated growth arrest, and decreased dependence on exogenous growth factors. Proteomic and transcriptomic analyses provided evidence that MmuPV1 E6 increases the expression and steady state levels of a number of cellular proteins that promote cellular proliferation and other hallmarks of cancer. These results indicate that MmuPV1 E6 is a major driver of MmuPV1-induced pathogenesis.

High-risk human papillomaviruses (HPVs) are causally associated with multiple human cancers. Previous studies have shown that the HPV oncoprotein E7 induces immune suppression; however, the underlying mechanisms remain unknown. To understand the mechanisms by which HPV deregulates host immune responses in the tumor microenvironment, we analyzed gene expression changes of all known chemokines and their receptors using our global gene expression data sets from human HPV-positive and -negative head/neck cancer and cervical tissue specimens in different disease stages. We report that, while many proinflammatory chemokines increase expression throughout cancer progression, CXCL14 is dramatically downregulated in HPV-positive cancers. HPV suppression of CXCL14 is dependent on E7 and associated with DNA hypermethylation in the CXCL14 promoter. Using in vivo mouse models, we revealed that restoration of Cxcl14 expression in HPV-positive mouse oropharyngeal carcinoma cells clears tumors in immunocompetent syngeneic mice, but not in Rag1-deficient mice. Further, Cxcl14 reexpression significantly increases natural killer (NK), CD4(+) T, and CD8(+) T cell infiltration into the tumor-draining lymph nodes in vivo In vitro transwell migration assays show that Cxcl14 reexpression induces chemotaxis of NK, CD4(+) T, and CD8(+) T cells. These results suggest that CXCL14 downregulation by HPV plays an important role in suppression of antitumor immune responses. Our findings provide a new mechanistic understanding of virus-induced immune evasion that contributes to cancer progression.

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 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.

Several recent studies have converged upon the innate immune DNA cytosine deaminase APOBEC3B (A3B) as a significant source of genomic uracil lesions and mutagenesis in multiple human cancers, including those of the breast, head/neck, cervix, bladder, lung, ovary, and other tissues. A3B is upregulated in these tumor types relative to normal tissues, but the mechanism is unclear. Because A3B also has antiviral activity in multiple systems and is a member of the broader innate immune response, we tested the hypothesis that human papillomavirus (HPV) infection causes A3B upregulation. We found that A3B mRNA expression and enzymatic activity were upregulated following transfection of a high-risk HPV genome and that this effect was abrogated by inactivation of E6. Transduction experiments showed that the E6 oncoprotein alone was sufficient to cause A3B upregulation, and a panel of high-risk E6 proteins triggered higher A3B levels than did a panel of low-risk or noncancer E6 proteins. Knockdown experiments in HPV-positive cell lines showed that endogenous E6 is required for A3B upregulation. Analyses of publicly available head/neck cancer data further support this relationship, as A3B levels are higher in HPV-positive cancers than in HPV-negative cancers. Taken together with the established role for high-risk E6 in functional inactivation of TP53 and published positive correlations in breast cancer between A3B upregulation and genetic inactivation of TP53, our studies suggest a model in which high-risk HPV E6, possibly through functional inactivation of TP53, causes derepression of A3B gene transcription. This would lead to a mutator phenotype that explains the observed cytosine mutation biases in HPV-positive head/neck and cervical cancers.

Epstein-Barr virus (EBV) is a human herpesvirus that is associated with lymphomas as well as nasopharyngeal and gastric carcinomas. Although carcinomas account for almost 90% of EBV-associated cancers, progress in examining EBV's role in their pathogenesis has been limited by difficulty in establishing latent infection in nontransformed epithelial cells. Recently, EBV infection of human telomerase reverse transcriptase (hTERT)-immortalized normal oral keratinocytes (NOKs) has emerged as a model that recapitulates aspects of EBV infection in vivo, such as differentiation-associated viral replication. Using uninfected NOKs and NOKs infected with the Akata strain of EBV (NOKs-Akata), we examined changes in gene expression due to EBV infection and differentiation. Latent EBV infection produced very few significant gene expression changes in undifferentiated NOKs but significantly reduced the extent of differentiation-induced gene expression changes. Gene set enrichment analysis revealed that differentiation-induced downregulation of the cell cycle and metabolism pathways was markedly attenuated in NOKs-Akata relative to that in uninfected NOKs. We also observed that pathways induced by differentiation were less upregulated in NOKs-Akata. We observed decreased differentiation markers and increased suprabasal MCM7 expression in NOKs-Akata versus NOKs when both were grown in raft cultures, consistent with our transcriptome sequencing (RNA-seq) results. These effects were also observed in NOKs infected with a replication-defective EBV mutant (AkataΔRZ), implicating mechanisms other than lytic-gene-induced host shutoff. Our results help to define the mechanisms by which EBV infection alters keratinocyte differentiation and provide a basis for understanding the role of EBV in epithelial cancers.IMPORTANCE Latent infection by Epstein-Barr virus (EBV) is an early event in the development of EBV-associated carcinomas. In oral epithelial tissues, EBV establishes a lytic infection of differentiated epithelial cells to facilitate the spread of the virus to new hosts. Because of limitations in existing model systems, the effects of latent EBV infection on undifferentiated and differentiating epithelial cells are poorly understood. Here, we characterize latent infection of an hTERT-immortalized oral epithelial cell line (NOKs). We find that although EBV expresses a latency pattern similar to that seen in EBV-associated carcinomas, infection of undifferentiated NOKs results in differential expression of a small number of host genes. In differentiating NOKs, however, EBV has a more substantial effect, reducing the extent of differentiation and delaying the exit from the cell cycle. This effect may synergize with preexisting cellular abnormalities to prevent exit from the cell cycle, representing a critical step in the development of cancer.

Up to 95% of all anal cancers are associated with infection by human papillomavirus (HPV); however, no established preclinical model exists for high-grade anal disease and cancer mediated by a natural papillomavirus infection. To establish an infection-mediated model, we infected both immunocompromised NSG and immunocompetent FVB/NJ mice with the recently discovered murine papillomavirus MmuPV1, with and without the additional cofactors of UV B radiation (UVB) and/or the chemical carcinogen 7,12-dimethylbenz(a)anthracene (DMBA). Infections were tracked via lavages and swabs for MmuPV1 DNA, and pathology was assessed at the endpoint. Tissues were analyzed for biomarkers of viral infection and papillomavirus-mediated disease, and the localization of viral infection was investigated using biomarkers to characterize the anal microanatomical zones. IMPORTANCE We show, for the first time, that MmuPV1 infection is sufficient to efficiently mediate high-grade squamous intraepithelial lesions in the anal tract of mice using the NSG immunocompromised strain and that MmuPV1, in combination with the chemical carcinogen DMBA, has carcinogenic potential. We further show that MmuPV1 is able to persist for up to 6 months in the anal tract of FVB/NJ mice irradiated with UVB and contributes to high-grade disease and cancer in an immunocompetent strain. We demonstrate that MmuPV1 preferentially localizes to the anal transition zone and that this localization is not an artifact of infection methodology. This study presents a valuable new preclinical model for studying papillomavirus-mediated anal disease driven by a natural infection.