Human papillomaviruses (HPVs) are the most common sexually transmitted infection in the United States and are a major etiological agent of cancers in the anogenital tract and oral cavity. Growing evidence suggests changes in the host microbiome are associated with the natural history and ultimate outcome of HPV infection. We sought to define changes in the host cervicovaginal microbiome during papillomavirus infection, persistence, and pathogenesis using the murine papillomavirus (MmuPV1) cervicovaginal infection model. Cervicovaginal lavages were performed over a time course of MmuPV1 infection in immunocompetent female FVB/N mice and extracted DNA was analyzed by qPCR to track MmuPV1 viral copy number. 16S ribosomal RNA (rRNA) gene sequencing was used to determine the composition and diversity of microbial communities throughout this time course. We also sought to determine whether specific microbial communities exist across the spectrum of MmuPV1-induced neoplastic disease. We, therefore, performed laser-capture microdissection to isolate regions of disease representing all stages of neoplastic disease progression (normal, low- and high-grade dysplasia, and cancer) from female reproductive tract tissue sections from MmuPV1-infected mice and performed 16S rRNA sequencing. Consistent with other studies, we found that the natural murine cervicovaginal microbiome is highly variable across different experiments. Despite these differences in initial microbiome composition between experiments, we observed that MmuPV1 persistence, viral load, and severity of disease influenced the composition of the cervicovaginal microbiome. These studies demonstrate that papillomavirus infection can alter the cervicovaginal microbiome.IMPORTANCEHuman papillomaviruses (HPVs) are the most common sexually transmitted infection in the United States. A subset of HPVs that infect the anogenital tract (cervix, vagina, anus) and oral cavity cause at least 5% of cancers worldwide. Recent evidence indicates that the community of microbial organisms present in the human cervix and vagina, known as the cervicovaginal microbiome, plays a role in HPV-induced cervical cancer. However, the mechanisms underlying this interplay are not well-defined. In this study, we infected the female reproductive tract of mice with a murine papillomavirus (MmuPV1) and found that key aspects of papillomavirus infection and disease influence the host cervicovaginal microbiome. This is the first study to define changes in the host microbiome associated with MmuPV1 infection in a preclinical animal model of HPV-induced cervical cancer. These results pave the way for using MmuPV1 infection models to further investigate the interactions between papillomaviruses and the host microbiome.

Human papillomavirus (HPV)-induced oropharyngeal cancer now exceeds HPV-induced cervical cancer, with a noticeable sex bias. Although it is well established that women have a more proficient immune system, it remains unclear whether immune control of oral papillomavirus infections differs between sexes. In the current study, we use genetically modified mice to target CCR2 and Stat1 pathways, with the aim of investigating the role of both innate and adaptive immune responses in clearing oral papillomavirus, using our established papillomavirus (MmuPV1) infection model. Persistent oral MmuPV1 infection was detected in Rag1ko mice with T and B cell deficiencies. Meanwhile, other tested mice were susceptible to MmuPV1 infections but were able to clear the virus. We found sex differences in key myeloid cells, including macrophages, neutrophils, and dendritic cells in the infected tongues of wild type and Stat1ko mice but these differences were not observed in CCR2ko mice. Intriguingly, we also observed a sex difference in anti-MmuPV1 E4 antibody levels, especially for two IgG isotypes: IgG2b and IgG3. However, we found comparable numbers of interferon-gamma-producing CD8 T cells stimulated by E6 and E7 in both sexes. These findings suggest that males and females may use different components of innate and adaptive immune responses to control papillomavirus infections in the MmuPV1 mouse model. The observed sex difference in immune responses, especially in myeloid cells including dendritic cell (DC) subsets, may have potential diagnostic and prognostic values for HPV-associated oropharyngeal cancer.

Select protease inhibitors (PI) have been found to be effective in decreasing human papillomavirus oncoprotein expression. This study evaluated whether the topical PI, Saquinavir (SQV), promotes viral clearance in an infectious mouse model with Mus musculus papillomavirus 1 (MmuPV1). NOD scid gamma (NSG) mice were anally infected with ∼4 × 108 viral genome equivalents of MmuPV1 and 120 days post-infection (when majority have high-grade anal dysplasia), began topical treatments: control (mock), 7,12-dimethylbenz(a)anthracene (DMBA) only, once weekly to promote carcinogenesis, 1% SQV only, daily (Monday - Friday), and SQV + DMBA. Viral MmuPV1 load was analyzed from anal lavages pre and post-treatment. Anal tissue was harvested, processed, and evaluated for drug absorption, grade of anal disease, and anal viral RNA. Results suggest that topical SQV promotes decreased viral shedding in female mice treated with SQV.

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.

Human papillomavirus (HPV) is a significant health burden and leading cause of virus-induced cancers. However, studies have been hampered due to restricted tropism that makes production and purification of high titer virus problematic. This issue has been overcome by developing alternative HPV production methods such as virus-like particles (VLPs), which are devoid of a native viral genome. Structural studies have been limited in resolution due to the heterogeneity, fragility, and stability of the VLP capsids. The mouse papillomavirus (MmuPV1) presented here has provided the opportunity to study a native papillomavirus in the context of a common laboratory animal. Using cryo EM to solve the structure of MmuPV1, we achieved 3.3 Å resolution with a local symmetry refinement method that defined smaller, symmetry related subparticles. The resulting high-resolution structure allowed us to build the MmuPV1 asymmetric unit for the first time and identify putative L2 density. We also used our program ISECC to quantify capsid flexibility, which revealed that capsomers move as rigid bodies connected by flexible linkers. The MmuPV1 flexibility was comparable to that of a HPV VLP previously characterized. The resulting MmuPV1 structure is a promising step forward in the study of papillomavirus and will provide a framework for continuing biochemical, genetic, and biophysical research for papillomaviruses.

Papillomavirus L1 and L2, the major and minor capsid proteins, play significant roles in viral assembly, entry, and propagation. In the current study, we investigate the impact of L1 and L2 on viral life cycle and tumor growth with a newly established mouse papillomavirus (MmuPV1) infection model. MmuPV1 L1 knockout, L2 knockout, and L1 plus L2 knockout mutant genomes (designated as L1ATGko-4m, L2ATGko, and L1-L2ATGko respectively) were generated. The mutants were examined for their ability to generate lesions in athymic nude mice. Viral activities were examined by qPCR, immunohistochemistry (IHC), in situ hybridization (ISH), and transmission electron microscopy (TEM) analyses. We demonstrated that viral DNA replication and tumor growth occurred at both cutaneous and mucosal sites infected with each of the mutants. Infections involving L1ATGko-4m, L2ATGko, and L1-L2ATGko mutant genomes generally resulted in smaller tumor sizes compared to infection with the wild type. The L1 protein was absent in L1ATGko-4m and L1-L2ATGko mutant-treated tissues, even though viral transcripts and E4 protein expression were robust. Therefore, L1 is not essential for MmuPV1-induced tumor growth, and this finding parallels our previous observations in the rabbit papillomavirus model. Very few viral particles were detected in L2ATGko mutant-infected tissues. Interestingly, the localization of L1 in lesions induced by L2ATGko was primarily cytoplasmic rather than nuclear. The findings support the hypothesis that the L2 gene influences the expression, location, transport, and assembly of the L1 protein in vivo.

HPV infections in the oral cavity that progress to cancer are on the increase in the USA. Model systems to study co-factors for progression of these infections are lacking as HPVs are species-restricted and cannot grow in preclinical animal models. We have recently developed a mouse papillomavirus (MmuPV1) oral mucosal infection model that provides opportunities to test, for the first time, the hypothesis that tobacco carcinogens are co-factors that can impact the progression of oral papillomas to squamous cell carcinoma (SCC). Four cohorts of mice per sex were included: (1) infected with MmuPV1 and treated orally with DMSO-saline; (2) infected with MmuPV1 and treated orally with the tobacco carcinogen, dibenzo[def,p]chrysene (DBP); (3) uninfected and treated orally with DMSO-saline, and (4) uninfected and treated orally with DBP. Oral swabs were collected monthly for subsequent assessment of viral load. Oral tissues were collected for in situ viral DNA/RNA detection, viral protein staining, and pathological assessment for hyperplasia, papillomas and SCC at study termination. We observed increased rates of SCC in oral tissue infected with MmuPV1 and treated with DBP when compared to mice treated with DBP or virus individually, each of which showed minimal disease. Virally-infected epithelium showed strong levels of viral DNA/RNA and viral protein E4/L1 staining. In contrast, areas of SCC showed reduced viral DNA staining indicative of lower viral copy per nucleus but strong RNA signals. Several host markers (p120 ctn, p53, S100A9) were also examined in the mouse oral tissues; of particular significance, p120 ctn discriminated normal un-infected epithelium from SCC or papilloma epithelium. In summary, we have confirmed that our infection model is an excellent platform to assess the impact of co-factors including tobacco carcinogens for oral PV cancerous progression. Our findings can assist in the design of novel prevention/treatment strategies for HPV positive vs. HPV negative disease.

The papillomavirus E5 gene contributes to transformation and tumorigenesis; however, its exact function in these processes and viral pathogenesis is unclear. While E5 is present in high-risk mucosotropic HPVs that cause anogenital and head and neck cancers, it is absent in cutaneous HPVs and the recently discovered mouse papillomavirus (MmuPV1), which causes papillomas and squamous cell carcinomas of the skin and mucosal epithelia in laboratory mice. We infected K14E5 transgenic mice, which express the high-risk mucosotropic HPV16 E5 gene in stratified epithelia, with MmuPV1 to investigate the effects of E5 on papillomavirus-induced pathogenesis. Skin lesions in MmuPV1-infected K14E5 mice had earlier onset, higher incidence, and reduced frequency of spontaneous regression compared to those in non-transgenic mice. K14E5 mice were also more susceptible to cervicovaginal cancers when infected with MmuPV1 and treated with estrogen compared to non-transgenic mice. Our studies support the hypothesis that E5 contributes to papillomavirus-induced pathogenesis.

Preclinical model systems to study multiple features of the papillomavirus life cycle are extremely valuable tools to aid our understanding of Human Papillomavirus (HPV) biology, disease progression and treatments. Mouse papillomavirus (MmuPV1) is the first ever rodent papillomavirus that can infect the laboratory strain of mice and was discovered recently in 2011. This model is an attractive model to study papillomavirus pathogenesis due to the ubiquitous availability of lab mice and the fact that this mouse species is easily genetically modifiable. Several other groups, including ours, have reported that MmuPV1-induced papillomas are restricted to T-cell deficient immunosuppressed mice. In our lab we showed for the first time that MmuPV1 causes skin cancers in UVB-irradiated immunocompetent animals. In this report we describe in detail the MmuPV1-UV infection model that can be adapted to study MmuPV1 biology in immunocompetent animals.

We report secondary cutaneous infections in the mouse papillomavirus (MmuPV1)/mouse model. Our previous study demonstrated that cutaneous MmuPV1 infection could spread to mucosal sites. Recently, we observed that mucosal infections could also spread to various cutaneous sites including the back, tail, muzzle and mammary tissues. The secondary site lesions were positive for viral DNA, viral capsid protein and viral particles as determined by in situ hybridization, immunohistochemistry and transmission electron microscopy analyses, respectively. We also demonstrated differential viral production and tumour growth at different secondarily infected skin sites. For example, fewer viral particles were detected in the least susceptible back tissues when compared with those in the infected muzzle and tail, although similar amounts of viral DNA were detected. Follow-up studies demonstrated that significantly lower amounts of viral DNA were packaged in the back lesions. Lavages harvested from the oral cavity and lower genital tracts were equally infectious at both cutaneous and mucosal sites, supporting the broad tissue tropism of this papillomavirus. Importantly, two secondary skin lesions on the forearms of two mice displayed a malignant phenotype at about 9.5 months post-primary infection. Therefore, MmuPV1 induces not only dysplasia at mucosal sites such as the vagina, anus and oral cavity but also skin carcinoma at cutaneous sites. These findings demonstrate that MmuPV1 mucosal infection can be spread to cutaneous sites and suggest that the model could serve a useful role in the study of the viral life cycle and pathogenesis of papillomavirus.