Advertisement

  • Loading metrics

Open Access

Peer-reviewed

Research Article

Role of Ultraviolet Radiation in Papillomavirus-Induced Disease

  • Aayushi Uberoi,

    Affiliation McArdle Laboratory for Cancer Research, Department of Oncology, School of Medicine and Public Health, University of Wisconsin-Madison, Madison, Wisconsin, United States of America

    ⨯
  • Satoshi Yoshida,

    Affiliation McArdle Laboratory for Cancer Research, Department of Oncology, School of Medicine and Public Health, University of Wisconsin-Madison, Madison, Wisconsin, United States of America

    ⨯
  • Ian H. Frazer,

    Affiliation The University of Queensland Diamantina Institute, Translational Research Institute, Brisbane, Queensland, Australia

    ⨯
  • Henry C. Pitot,

    Affiliation McArdle Laboratory for Cancer Research, Department of Oncology, School of Medicine and Public Health, University of Wisconsin-Madison, Madison, Wisconsin, United States of America

    ⨯
  • Paul F. Lambert

    plambert@wisc.edu

    Affiliation McArdle Laboratory for Cancer Research, Department of Oncology, School of Medicine and Public Health, University of Wisconsin-Madison, Madison, Wisconsin, United States of America

    ⨯

Abstract

Human papillomaviruses are causally associated with 5% of human cancers. The recent discovery of a papillomavirus (MmuPV1) that infects laboratory mice provides unique opportunities to study the life cycle and pathogenesis of papillomaviruses in the context of a genetically manipulatable host organism. To date, MmuPV1-induced disease has been found largely to be restricted to severely immunodeficient strains of mice. In this study, we report that ultraviolet radiation (UVR), specifically UVB spectra, causes wild-type strains of mice to become highly susceptible to MmuPV1-induced disease. MmuPV1-infected mice treated with UVB develop warts that progress to squamous cell carcinoma. Our studies further indicate that UVB induces systemic immunosuppression in mice that correlates with susceptibility to MmuPV1-associated disease. These findings provide new insight into how MmuPV1 can be used to study the life cycle of papillomaviruses and their role in carcinogenesis, the role of host immunity in controlling papillomavirus-associated pathogenesis, and a basis for understanding in part the role of UVR in promoting HPV infection in humans.

Author Summary

Epidemiological studies have implicated that ultraviolet radiation (UVR) from sunlight drive papillomavirus-induced disease in healthy as well as immunocompromised humans. In this report we demonstrate that treatment of immunocompetent mice with UVR renders them susceptible to papillomas and associated squamous cell carcinoma when infected with the recently discovered murine papillomavirus (MmuPV1). Our data further suggest UVR increases susceptibility to virally induced disease by inducing immunosuppression.

Citation: Uberoi A, Yoshida S, Frazer IH, Pitot HC, Lambert PF (2016) Role of Ultraviolet Radiation in Papillomavirus-Induced Disease. PLoS Pathog 12(5): e1005664. https://doi.org/10.1371/journal.ppat.1005664

Editor: Denise A. Galloway, Fred Hutchinson Cancer Research Center, UNITED STATES

Received: February 11, 2016; Accepted: May 6, 2016; Published: May 31, 2016

Copyright: © 2016 Uberoi et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Data Availability: All relevant data are within the paper and its Supporting Information files.

Funding: This study was supported by grants from the National Cancer Institute (CA022443) and the National Institute of Arthritis, Musculoskelatal and Skin diseases (AR066524). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Competing interests: The authors have declared that no competing interests exist.

Introduction

Papillomaviruses are species-specific, epitheliotropic, double-stranded DNA viruses. There are over 200 strains or genotypes of human papillomaviruses (HPVs) [1]. Mucosotropic HPVs are the most common sexually transmitted pathogens, and a subset of these viruses cause 5% of human cancers, including cervical cancer, other anogenital cancers, and a growing fraction of head and neck cancers (reviewed in [2]). Other HPVs cause cutaneous warts, which are among the most common ailments treated by dermatologists [3–6]. They arise most frequently among children [7,8], and impose a significant burden in immunocompromised patients, particularly amongst organ transplant recipients [9–11]. They are ubiquitous in nature and can persist in the skin asymptomatically for years most clearly in context of immunosuppressed patients [9,12]. A subset of cutaneous HPVs also has been causally associated with skin cancer (reviewed in [9,13,14]).

The study of papillomavirus-induced disease has long been hindered by the absence of any identified strains of virus that infect laboratory mice. This limitation was overcome with the recent identification of the murine papillomavirus, MmuPV1, isolated from cutaneous warts arising on the T-cell deficient NMRI-FoxN1nu/nu strain of laboratory mice [15]. MmuPV1 belongs to the pi-papillomaviridae genus and is phylogenetically related to cutaneous HPVs and other animal PVs that cause cutaneous disease in exotic rodent species [16]. MmuPV1 causes warts in cutaneous epithelium as well as in mucosal epithelium lining the female reproductive tract and oral cavity [15,17–20], and in some cases these lesions show signs of neoplastic progression [20,21]. Multiple studies have shown that the ability of MmuPV1 to cause overt disease is largely restricted to immunodeficient strains of mice [18,20].

Epidemiological studies have suggested that there is a correlation between exposure to ultraviolet radiation (UVR) and the prevalence of cutaneous HPVs in healthy and immunosuppressed patients, respectively [22,23]. Cutaneous HPVs are more commonly found at anatomical sites exposed to sunlight, and a history of blistering sunburn is associated with prevalent and persistent cutaneous HPV infections [12,22–24]. UVR has also been shown to play a role in papillomavirus-associated disease caused by animal papillomaviruses that infect the African multimammate rat (MnPV) or the cottontail rabbit (CRPV) [25–27].

In this report we demonstrate that UVR, specifically UVB (280 to 315 nm) assists in development of MmuPV1 dependent papillomas and associated malignant progression to squamous cell carcinomas, in immunocompetent strains of mice. We further show that there is a correlation between UVR-induced susceptibility to MmuPV-1 associated disease and UVR-induced immunosuppression. These findings provide a potential explanation for the role of UVR-mediated immunosuppression in papillomavirus-associated disease in humans.

Results

MmuPV1 causes papillomas in immunocompetent FVB mice following exposure to UV irradiation

Several lines of evidence support a link between UVR exposure and papillomavirus infection [22,23,26]. We sought to determine if UVR facilitates MmuPV1-induced papillomatosis in immunocompetent strains of mice. We used the inbred FVB/NJ strain of mice for our initial studies because it has been classically used to study chemically induced skin tumorigenesis [28].

Ears and tails of 8–10 weeks old FVB/NJ mice were infected with 108 viral genome equivalents (VGE, a measure of the amount of encapsidated genomes in a stock of virus) of MmuPV1 virions following topical scarification of the epidermis. Twenty-four hours post-infection (h.p.i.), mice were or were not exposed to varying doses of UVB whole body irradiation (280-320nm range UVR). Mice were then observed weekly for papillomatosis (Fig 1A). By 3 months post-infection greater than 50% of infected ear sites developed papillomas in FVB/NJ mice treated with 300mJ/cm2 UVB (Fig 1B, Table 1, S1 Fig-Panel A). Sites on the tails of the same mice infected with an equal dose of virus failed to develop papillomas. We did not see any papillomas develop at infected sites on non-irradiated mice or mice treated with a lower dose (150mJ/cm2) of UVB (Table 1). Mock-infected mice exposed to the same doses of UVB did not develop warts (Table 1). Lateral transmission of MmuPV1 has been observed in immunodeficient strains of mice experimentally infected with MmuPV1 [15,20]; however, we did not see any warts arise at uninfected sites on the UVB-treated, MmuPV1-infected FVB/NJ mice.

thumbnail
Download:
Table 1. Papilloma incidence in UVB-irradiated FVB/NJ mice.

https://doi.org/10.1371/journal.ppat.1005664.t001

thumbnail
Download:
Fig 1. MmuPV1-induced papillomas in UVB-irradiated FVB/NJ mice.

(A) Sequence of manipulations to FVB/NJ mice in the MmuPV1-UVB infection model system. Ear sites were scarified and exposed to 108 VGE of MmuPV1 either 24 hours before (FVB+UVB 24h.b.i.+MmuPV) or 24 hours after (FVB+MmuPV+UVB 24h.p.i.) 300mJ/cm2 UVB whole-body irradiation. Sites were then scored weekly for presence of papillomas. (B) Kaplan-Meyer plots of the fraction of papilloma-free infected ear sites with respect to time. There is no significant difference in the temporal onset of disease between the two experimental groups (FVB+UVB24h.p.i.+MmuPV vs FVB+MmuPV+UVB 24h.b.i.; P = 0.814, Wilcoxon log-rank test, two-sided). (C) Percentage of sites with overt papillomas over a 6-month observation period. In both the cohorts of mice, a subset of papillomas completely regressed by 6 months post-infection. (D) Distribution of ear sites of UVB-treated FVBN/J mice (FVB+UVB24h.p.i.+MmuPV group) infected that developed papillomas that completely regressed (red), partially regressed (blue) or continued to grow (green) over the 6 month monitoring period. In grey is the fraction of sites that did not develop papillomas. (E) Growth profiles of individual papillomas arising in the FVB+UVB24h.p.i.+MmuPV group during the 6-month monitoring period. Individual lines represent the geometric mean diameter of each papilloma as a function of time. Red lines are papillomas that completed regressed. Blue lines are papillomas that partially regressed. Green lines are papillomas that continued to grow.

https://doi.org/10.1371/journal.ppat.1005664.g001

The presence and size of warts arising on the MmuPV1-infected ears sites of UVB-treated FVB/NJ mice were monitored for a total of 6 months. In this time frame 17% of the sites developed papillomas that completely regressed (Fig 1D and 1E: red lines in Fig 1E). 25% of sites developed papillomas that partially regressed (Fig 1D and 1E-blue lines). Another 17% of sites developed papillomas that continued to grow over the 6 month period (Fig 1D and 1E—green lines). 41% of the sites did not develop papillomas.

To address whether UVB is exerting its effect directly on the virus, we tested whether UVB treatment at 300mJ/cm2 24 hours before infection (h.b.i.) had the same impact on papilloma incidence. Over the initial 4-month period post-infection, papillomas arose with similar frequency as in groups of mice exposed to UVB 24 hours after infection (Fig 1B, P = 0.814, log-rank, two-sided). Again, we saw complete regression of a similar fraction of papillomas when mice were observed up to 6 months post-infection (Fig 1C). We also found that a subset of animals infected with MmuPV1 and treated with a UVB fourteen days post-infection also developed papillomas at a similar frequency (S1 Table).

To learn if the UVA spectra (315-400nm range UVR) also caused immunocompetent mice to become susceptible to MmuPV1-induced papillomatosis, we infected FVB/NJ mice with 108 VGE of MmuPV1 at sites on the ears and exposed them to 300 J/cm2 of UVA 24 hours post-infection. We failed to observe papillomas arise at any sites on these mice (S2 Table). Combining UVA (300 J/cm2) and UVB (300 mJ/cm2) treatment gave a similar incidence of papillomas as seen with UVB alone (S2 Table).

We also tested the susceptibility of other commonly used, wild-type, inbred strains of mice to MmuPV1-induced disease with varying doses of UVB (S3 Table): 30% of infected ear sites on C57/BL6 mice developed papillomas at the three-month time point when treated with UVB at a higher dose of 600mJ/cm2. Only one infected site on a BALB/c mouse developed a papilloma when exposed to 1200mJ/cm2 UVB (S3 Table). These results indicate that there are strain differences in susceptibility of mice to MmuPV1-induced disease following UVB irradiation. Because UVB-treated FVB/NJ mice were most susceptible to MmuPV1 disease (p = 0.041 FVB/NJ vs C57/Bl6 or 0.0044 FVB/NJ vs BalbC, at 300 mJ/cm2), we pursued all further studies of the role of UVB in inducing MmuPV1 infection using the FVB/NJ strain of mice.

To investigate the relationship between viral dose and papilloma incidence in FVB/NJ mice treated with UVB, we infected mice on their ears with 108, 107 or 106 VGE of MmuPV1 and exposed the mice to 300mJ/cm2 of UVB 24 hours post-infection (Fig 2). In parallel, T-cell deficient FoxN1nu/nu mice were also infected at ear sites with the same stock of virus at designated doses in the absence of UVB-treatment (Fig 2). By 3 months post-infection, 90% of ear sites in FoxN1nu/nu mice infected with 108 and 107 VGE MmuPV1 developed papillomas (Fig 2A and 2B –Right). At this time-point, papillomas developed in UVB-irradiated FVB/NJ mice infected with 108 and 107 VGE MmuPV1, but at lower penetrance (Fig 2A and 2B—Left) and no papillomas formed on UVB-treated FVB/NJ mice infected with 106 VGE of MmuPV1 up to 6 months post-infection (Fig 2B- Right). In contrast, 37.5% of FoxN1nu/nu mice ear sites infected with 106 VGE of MmuPV1 develop papillomas (Fig 2B-Left). Non-irradiated FVB/NJ mice infected with 108, 107 or 106 VGE of MmuPV1 did not develop any papillomas. These results indicate that a higher threshold in the amount of virus is required to see papilloma induction in immunocompetent animals after UVB irradiation. Furthermore, there was complete regression of a subset of papillomas in FVB/NJ mice infected with 107 and 108 VGE of MmuPV1 when monitored up to 6 months post-infection (Fig 2B). There was no regression of papillomas in immunodeficient FoxN1nu/nu mice.

thumbnail
Download:
Fig 2. Viral dose-dependent incidence of papillomas in UVB-treated FVB/NJ mice and in FoxN1nu/nu mice.

Immunocompetent FVB/NJ mice were infected at ear sites with 108, 107 or 106 VGE of MmuPV1 virions following scarification. Twenty-four hours post-infection (h.p.i.), mice were irradiated with 300mJ/cm2 UVB and scored weekly for presence of ear papillomas up to 6 months. In parallel, T-cell deficient FoxN1nu/nu mice not treated with UVB were infected at ear sites with the same stock of virus at designated doses. Graphs for each strain are shown separately. (A) Kaplan-Meyer plot of the fraction of papilloma-free infected sites over the initial four-month period following infection. There was no significant difference between incidences of disease between UVB-irradiated FVB mice treated with 107 versus 108 VGE of MmuPV1 (P = 0.21, Wilcoxon log-rank test, two-sided). FVB/NJ mice infected with 106 VGE of MmuPV1 did not develop any papillomas. (B) Percentage of sites with overt papillomas over a 6-month observation period. Note that in both groups, UVB-irradiated FVB mice infected with 107 and 108 VGE of MmuPV1, a subset of papillomas completely regressed by 6 months post-infection (left). There was no regression of papillomas in FoxN1nu/nu mice (right).

https://doi.org/10.1371/journal.ppat.1005664.g002

MmuPV1-induced papillomas in UVB-irradiated immunocompetent mice show productive infection and signs of neoplastic transformation

Lesions arising on the ear sites of MmuPV1-infected FVB mice treated with UVB were harvested at 6 months post-infection and analyzed histopathologically. They showed sessile papilloma-like morphology with hyperkeratosis (Fig 3A, left; S1 Fig-Panel B). Cells within the stratum granulosum showed presence of koilocytes consistent with productive papillomavirus infection (Fig 3A, middle). The MmuPV1 major viral capsid protein was expressed in these papillomas as indicated by L1-specific immunofluorescence (Fig 3A, right). Viral capsid protein was mostly detected in the koilocytes and other cells within the upper layers of the epithelia (suprabasal and terminally differentiated) with some L1-positive cells occasionally observed in the basal layer consistent with prior findings in immunodeficient mice [19]. The sessile papillomas were accompanied by multiple areas of atypical squamous cell hyperplasia, several of which were suggestive of early neoplastic transformation, consistent with other reports [17,19,20]. In our case, however, we also found focal areas of malignant progression consistent with squamous cell carcinoma with invasion extending into follicular structures deep within the dermis as evident from cytokeratin-14 staining (Fig 3B). These areas of malignant progression were observed in several histopathologically scored lesions (e.g. see S1 Fig- Panel B). Foci of chronic inflammation were also noted in the lesions, some of which extended into the atypical dermis (Fig 3B).

thumbnail
Download:
Fig 3. Histopathological analysis of MmuPV1-induced papillomas in UVB –irradiated mice.

(A) H&E staining of ear tissue indicating presence of a sessile ear papilloma caused by MmuPV1 in FVB mice following UVB irradiation (left). Presence of koilocytes is indicated (middle). L1-immunofluorescent staining of papillomas (right). (B) Sessile papilloma with areas indicating frank malignant transformation (inset;i) and foci of inflammation (inset;ii). K14 staining of papilloma and associated epithelial invasion into underlying dermis. Papillomas on the tail of an immunodeficient BALB/c FoxN1nu mice (right, top). These papillomas were formed as a result of infection with viral extracts prepared from papillomas that arose in immunocompetent FVB mice infected with MmuPV1 followed by UVB irradiation. (C) L1 (green)-MmuPV1 FISH (red) co-staining of MmuPV1-induced ear papillomas in UVB irradiated FVB/NJ mice. The nuclei were counterstained with DAPI (blue). Scale bars represent 100μm.

https://doi.org/10.1371/journal.ppat.1005664.g003

Southern blot quantification showed that viral extracts harvested from papillomas from immunocompetent mice had a 1000-fold reduction in the amount of VGE/mg of wart harvested when compared to the papillomas arising in immunocompromised FoxN1nu mice (S2 Fig). Regardless, the viral extracts from immunocompetent mice were infectious and caused papillomas in FoxN1nu mice (Fig 3B). MmuPV1 DNA-specific in situ hybridization coupled with L1 immunohistochemistry analysis of MmuPV1-induced ear papillomas in UVB irradiated FVB/NJ mice (Fig 3C) showed presence of amplified viral DNA and L1 capsid. Robust L1 expression was seen throughout the papilloma most frequently in the suprabasal layers of the epithelia. The FISH positive cells were comparatively less frequent and predominantly seen in the spinous epithelium occasionally showing co-localization with L1 positive nuclei. We observed that while regions of papillomatosis showed presence of amplified MmuPV1 DNA areas of malignancy showed little to no presence of MmuPV1 amplified DNA (S3 Fig). Presence of the viral major capsid protein (L1) and amplified viral DNA, which are markers for the productive stage of the viral life cycle [reviewed in [29]], coupled with findings from the transmission experiments (Fig 3B), confirm that MmuPV1 establishes a productive infection in immunocompetent mice following UVB treatment.

UVB irradiation has a systemic effect on host biology

There are several lines of evidence in the field of photoimmunology indicating that UVB impairs a variety of immune responses in humans and laboratory animals both locally, within UV-irradiated skin, and systemically, at distant sites [30–32]. To test whether UVB-assisted pathogenesis in MmuPV1-infected FVB/NJ mice is due to a systemic or a local effect of UVB-irradiation, we infected FVB/NJ mice at ear sites as described previously. Twenty-four hours after infection, mice were anesthetized and the infected ear sites were shielded with tin foil leaving the rest of the mouse exposed to 300mJ/cm2 UVB irradiation (Fig 4A). Mice were then observed weekly to score for papillomatosis up to six months post-infection. At the end of 16 weeks approximately 55% of infected mice developed papillomas on ear sites (Fig 4B). There was no significant difference in the temporal onset of papillomas between shielded and unshielded mice (P = 0.938, log-rank, two sided). The animals were kept under observation for 6 months and scored weekly for papilloma incidence. We found that by six months post-infection 27% of the papillomas completely regressed in the mouse cohort whose infection sites were shielded, similar to that seen in the unshielded cohort, where 28.5% of the papillomas had completely regressed. These results indicate that systemic effects of UVB on host biology must contribute to the ability of UVB to induce MmuPV1-dependent disease.

thumbnail
Download:
Fig 4. Systemic effect of UVB assists in MmuPV-dependent papillomatosis.

(A) Schematic illustrating the experimental design. Twenty-four hours post infection with 108 VGE, groups of mice were either whole body UVB irradiated (FVB+MmuPV+UVB 24h.p.i.) or infected sites were shielded from UVB (FVB+MmuPV+UVB 24h.p.i. shielded) and scored weekly for presence of ear papillomas up to 6 months. UVB dose was 300mJ/cm2. Scoring data is compared to the FVB+MmuPV+UVB 24h.p.i. group also shown in Fig 2, as both experiments were performed at the same time. (B) Kaplan-Meyer plots of the fraction of papilloma-free infected sites over the first 12 weeks post-infection. There was no significant difference between the two experimental groups (shielded versus unshielded ears) (P = 0.938) as assessed by Wilcoxon log-rank analysis. (C) Percentage of sites with overt papillomas over a 6-month observation period. Some papillomas completely regressed in both groups. Control animals i.e. non-UVB irradiated FVB mice infected with MmuPV1 (FVB+MmuPV-UVB) did not develop any papillomas over the 6 month observation period.

https://doi.org/10.1371/journal.ppat.1005664.g004

UVB-mediated immunosuppression correlates with development of MmuPV1- dependent cutaneous disease

Through the study of different strains of immunodeficient mice, it has been found that T-cell deficiency is necessary for the development of MmuPV1-dependent papillomas [18,20]. Several studies have demonstrated that UVB causes cell-mediated immunosuppression in mice [33,34], reviewed in [30,32,35,36]. Cell-mediated immunosuppression by UVB has traditionally been measured by monitoring delayed type hypersensitivity (DTH) responses [30,32,35,37,38]. A single exposure to UVB irradiation is sufficient to inhibit DTH responses in some strains of mice [39,40]. To determine if the single dose of UVB irradiation (300mJ/cm2) that makes FVB/NJ mice susceptible to MmuPV1-induced papillomatosis (Table 1) is sufficient to cause immunosuppression in this strain of mice, we measured DTH responses in these and control mice not treated with UVB [41,42]. Ten days post-UVB treatment, FVB/NJ mice were sensitized to an antigen by topically applying 0.5 mg 1-Chloro-2, 4-Di-Nitrobenzene (DNCB) on the shaved backs of the mice. Five days later we determined the level of immune response to this antigen by applying 0.2 mg DNCB to a distant site (the ears) and monitoring DTH responses (Fig 5). DTH was assessed by measuring ear swelling every 24 hours for four days. A single exposure to UVB (300mJ/cm2) was capable of causing immunosuppression in FVB/NJ mice as evidenced by a marked decrease in swelling in response to DNCB challenge in the UVB treated mice compared to non-UVB treated mice (Fig 5). This difference was statistically significant (p<0.005, T-test). To assess whether this immunosuppression was systemic or local, we repeated the DTH assay on mice in which we shielded their ears from UVB as described previously. These mice also displayed reduced swelling in response to challenge with DNCB on their ears, indicating that UVB-induced immunosuppression in these mice is systemic (Fig 5). This observation suggests that UVB causes systemic immunosuppression that may assist MmuPV1-dependent disease (Figs 1–4).

thumbnail
Download:
Fig 5. UVB causes systemic immunosuppression of host.

Three groups of mice were treated as follows: a control group that was not UVB irradiated (FVB-UVB, White bars, n = 3), UVB-irradiated group (FVB+UVB, Striped bars, n = 6) and UVB-irradiated group in which ears were shielded from UVB exposure [FVB+UVB (ears shielded), Black bars]. Ear thickness reported as the average of the difference between DNCB-challenged (left) and unchallenged (right) ears four days post challenge. SEM refers to standard error in measurement computed by measuring standard deviation. DTH measured by the change in thickness of the ear (ΔEar) was found to be statistically significant between UVB-treated mice and non-irradiated control mice (student t-test, p<0.005, two-sided) for both UVB-treated groups.

https://doi.org/10.1371/journal.ppat.1005664.g005

Mice displaying long-term UVB-induced immune suppression preferentially develop MmuPV1-induced warts

To determine if UVB-mediated immune suppression correlates with papilloma incidence we performed long-term DTH assays in FVB/NJ mice following infection with MmuPV1 and UVB-treatment. Mice were infected with MmuPV1 in their right ear and exposed to a single dose of UVB (300 mJ/cm2) 24 hrs post-infection as described previously (Fig 1A). Ten days post UVB-treatment animals were sensitized to antigen by painting 0.5 mg of DNCB on their shaved backs. Mice were monitored for papilloma formation. Three months post infection, animals were challenged with 0.2 mg DNCB on the uninfected, left ear. We chose the three-month time period because this is the time by which FVB/NJ mice treated with UVB develop the maximal number of MmuPV1-induced warts (Figs 1C and 2B). Post-challenge we measured ear swelling up to 96 hrs (Fig 6). We found that in the control group (animals not treated with UVB) a long-lived DTH response was established as evidenced by swelling of the ears. UVB-treated, uninfected control mice displayed a range in levels of DTH response, indicative of variable levels of long-term immune suppression. This was also seen in UVB-treated, MmuPV1-infected mice. Of specific note, the mice in this latter group that retained immune-suppression at the end of 3 months were the same ones that had developed MmuPV1-induced warts. The difference in ear swelling between the UVB-irradiated animals infected with MmuPV1 that developed warts and those that did not develop warts was statistically significant (P = 0.016, Wilcoxon rank-sum, two-sided). These data demonstrate a strong correlation between long-term, UVB-induced immunosuppression and MmuPV1-dependent pathogenesis.

thumbnail
Download:
Fig 6. Long-lived immune suppression in UVB-irradiated mice correlates with papilloma incidence.

Mice were either infected with MmuPV1 or vehicle followed by UVB irradiation (300mJ/cm2) in the right ear. Mice were sensitized with DNCB 10 days post UVB exposure and were challenged DNCB in the left ear 3 months following infection. Ear swelling was measured by means of a Vernier calipers. Ear thickness is reported as the average of the difference between ear thickness 0 hrs post challenge and 72 hrs post challenge. Ear thickness for each mouse is shown as a dot plot. The black lines represent the mean reading for each group. Wilcoxon rank-sum test was used to analyze difference between several groups. There was significant difference between UVB-irradiated (UV) and control (no UV) groups (*p = 0.021, two-sided). There was significant difference between the UVB irradiated animals infected with MmuPV1 that developed warts and those that did not develop warts (**p = 0.016, two-sided). There was no significant difference between control (no UV) and UVB-irradiated mice infected with MmuPV1 that did not develop warts (MmuPV1+ UV without warts).

https://doi.org/10.1371/journal.ppat.1005664.g006

Discussion

The lack of infection models in a tractable laboratory animal has limited our ability to study the pathogenesis of papillomaviruses in their natural hosts. The murine papillomavirus, MmuPV1, isolated from cutaneous warts arising on immunodeficient NCR-FoxN1nu/nu laboratory mice [15,16] is a valuable animal papillomavirus because it provides us, for the first time, the opportunity to study papillomavirus infections in the context of a genetically manipulatable host.

Prior to this study, MmuPV1-associated papillomatosis has been described primarily in the context of immunodeficient strains of mice [17,18,20,21], though there have been reports that MmuPV1 can cause warts on hairless strains of mice, which are thought to be immunocompetent [15,43]. Consistent with a role of the host immune system in limiting MmuPV1-induced disease, Handisurya, et al. (2014) found that Cyclosporin A treatment is required for induction and maintenance of MmuPV1-induced papillomas in immunocompetent mice [18]. These authors did find that, at very high doses of virus (1012 VGEs), the SENCAR strain of mice, selectively bred for high susceptibility to skin tumor induction by chemical carcinogens [44], developed papillomas, but these papillomas regressed within two weeks of appearing.

In this study we investigated whether the UVB spectra impacts susceptibility of immunocompetent mice to MmuPV1. There was limited prior evidence suggesting the role of UV radiation in other animal models for papillomavirus infection. Mastomys natalensis Papillomavirus (MnPV) is a rodent papillomavirus that is shown to cause papillomatosis in the African multimammate rat [25]. MnPV DNA was found in UV induced tumors in HRA/Skh mice that are hairless but immunocompetent [26]. Further, studies showed that cell-free extracts containing MnPV enhanced UV-induced tumorigenesis [45]. We found that when FVB/NJ mice were exposed to high doses of UVB (Table 1) greater than 50% of infected sites developed papillomas by 3 months, some of which persisted for 6 months (Figs 1–4). UVB also induced MmuPV1-dependent papillomatosis in other strains of immunocompetent mice (S3 Table).

Histopathological analysis of the papillomas in the MmuPV1/UVB infection model also indicated progression to squamous cell carcinoma (Fig 3, S1 Fig-Panel B). Interestingly, we observed that while regions of papillomatosis showed many cells that had amplified viral DNA, we only observed a few cells harboring amplified viral DNA in papilloma-associated malignant regions (S3 Fig). More sensitive in situ hybridization techniques will be required to determine whether MmuPV1 viral genomes are lost as malignant progression arises, as is thought to occur in beta-HPV associated non-melanoma skin cancers [46].

Several lines of evidence supports the hypothesis that UVB is having an indirect effect on increasing the susceptibility to MmuPV-1 induced papillomatosis by inducing systemic immunosuppression. First, we did not find any significant difference in susceptibility to wart formation in mice that were treated with UVB 24 hours before or 24 hours after infection with MmuPV1 (Fig 1). Second, UVB's ability to increase susceptibility to MmuPV1-induced papillomatosis was equally efficient whether or not the site of infection was exposed to UVB (Fig 4). And third, long-lived immunosuppression correlated with wart formation (Fig 6).

A second but not mutually exclusive hypothesis is that UVB also directly influences MmuPV1-induced pathogenesis. Supportive of this hypothesis, expression of the viral genes in K14HPV8 transgenic mice were increased following exposure to UVR [47]. The physiological relevance of this observation remains unclear, however, because the HPV8 genes are under the control of a heterologous, keratin 14 promoter in the context of this transgenic mouse model. The MmuPV1 infection model should provide a valuable experimental platform for further testing this second hypothesis.

The ability of papillomaviruses to persist in their host in the absence of causing overt disease, i.e. latency, has long been suspected. Compelling evidence for latency comes from the observation that cyclosporine-induced immunosuppression led to elevation of the viral DNA copy number at sites of wart regression in cottontail rabbits infected with rabbit oral papillomavirus, consistent with reactivation of virus from latency [48]. In the MmuPV1-infection model, ELISA analysis of serum from C57BL/6 mice infected with MmuPV1 showed seroreactivity to MmuPV1 virus particles at 70 days post-infection even though they did not develop warts [49]. Likewise, seroconversion was seen in all mice (n = 20) in a group of SKH-1 mice infected with MmuPV1, of which 3 mice actually developed papillomas [43]. Together, these observations indicate that the virus is presented to the immune system even in the absence of causing overt disease. Whether this results from the original exposure of the animals to the virus or the consequence of a latent infection remains unclear. Our own observation that a single exposure to UVB 14 days post-infection led to papillomatosis at sites infected with MmuPV1 (S1 Table) indicates either that infectious virus is stably retained at the site of infection for that period of time and then initiates infection post UVB exposure or that latent infections arose that were then activated by UVB-induced immunosuppression. It remains to be determined whether latency arises in the MmuPV1 infection model, and, if so, the nature of this latency.

In this study we observed an unexpected tissue specificity in terms of susceptibility to MmuPV1-induced papillomatosis; while infection of the ear led to efficient formation of papillomas in UVB-treated FVB/NJ mice, papillomas did not arise at tail sites that were infected in the same mice. Others have reported that there is some site specificity for MmuPV1-induced disease in immunodeficient mice [19,20]. For example, while tail and muzzle of immunodeficient mice are susceptible to MmuPV1-induced papillomatosis, the torso skin is not. There are limited studies that have been directed towards looking at difference in susceptibility of different cutaneous sites of mice to cancer. Therefore the molecular differences between the sites are not very clear. In our studies with HPV 16 transgenic (K14E6, K14E7) mice we have consistently observed epithelial hyperplasia most extensively in the ear skin [50–52]. Likewise, ear skin of HPV 38 E6/E7 transgenic mice display patches of hyperproliferation [53]. Recently, there has been one report that suggests that difference in miR-155 expression of ear versus chest/torso skin of K14HPV16 mice could explain the difference in susceptibility of different tissues to develop HPV-mediated carcinogenesis [54]. It is also possible that grooming behaviors of mice, such as scratching of ear, promotes wounding [55] thereby allowing virus to access lower layers of epithelia better.

In our studies we observed a clear viral-dose dependence in causing disease. Interestingly, we found that a viral dose of 107 VGE or greater was required to see papillomas in UVB-irradiated FVB/NJ mice whereas 106 VGE of MmuPV1 were sufficient to induce papillomas on BALB/c-Foxn1nu/nu mice, albeit at low incidence over the same 6-month observation period (Fig 2). This suggests that there is a higher threshold of virus required to see MmuPV1-dependent disease in UVB treated immunocompetent animals. We posit that this observation reflects, at least in part, the fact that long-term immunosuppression induced by UVB in FVB/NJ mice only occurs in a subset of the UVB-treated mice (Fig 6).

In our testing of different genetic backgrounds, only FVB/NJ was susceptible to MmuPV1-associated pathogenesis at 300mJ/cm2, whereas C57/BL6 and BALB/c mice were susceptible to MmuPV1-associated pathogenesis at higher doses of UVB (S3 Table). This observation is consistent with a previous study that suggests that different genetic backgrounds respond differently to UVB in terms of levels of immunosuppression based upon DTH assays [56]. While that study did not test FVB/NJ mice, they found that C57/BL6 mice were more sensitive to UVB-induced immunosuppression than BALB/c mice. This correlates with the level of MmuPV1-dependent papillomatosis observed in our studies (S3 Table).

In our study we observed that UVB, but not UVA alone makes immunocompetent mice susceptible to MmuPV1-dependent papillomatosis (S2 Table). UVB has been shown to play a key role in initiating and mediating immunosuppression, whereas the mechanisms and roles of UVA in immunosuppression are not well understood [57,58]. UVB can cause both short-term as well as long-term defects in cell-mediated DTH responses [35,38,40], at least in part by inhibiting development of memory T-cells and by causing an overall reduction in T-cell subpopulations in the skin [59]. We observed that UVB caused variable levels of long-term immunosuppression in FVB/NJ animals and that MmuPV1-induced papillomas developed preferentially in those animals retaining long-term immunosuppression (Fig 6). This can explain the observation that papillomas arise in approximately 50% of animals. Based upon these findings we also posit that the observed regression of papillomas (Figs 1–2) correlates with a loss of immunosuppression.

The observation that there is a strong correlation between long-term immunosuppression induced by UVB and MmuPV1-dependent pathogenesis is perhaps the most notable finding of this study. This correlation supports the hypothesis that UVB-induced immunosuppression can help drive papillomavirus-induced disease. There is correlative epidemiological data from human studies in which anatomical sites on individuals that are exposed to sunlight, or at which sunburn has occurred are increased in their susceptibility to HPV-induced warts. One difference however, is that in our studies with FVB/NJ mice, the effect of UVB was found to be systemic; UVB irradiation did not have to be applied to the infection site. This raises the interesting question: is there a difference in the role of UVB in HPV-driven pathogenesis compared to MmupV1-driven pathogenesis? Further studies are needed to assess whether a local effect of UVB on MmuPV1 can be identified in mice. Recently, it has been shown that MmuPV1 also infects the mucosal epithelium of the female reproductive tract and oral cavity [21]. In this regards MmuPV1 infection model is truly unique and the systemic immunosuppression by UVB can further be tested in the context of mucosal disease.

In conclusion, we have reported the novel finding that UVR makes immunocompetent mice susceptible to development of MmuPV1-induced cutaneous papillomas, and this correlates with UVB-induced systemic immunosuppression. This observation opens the door to pursuing studies using genetically engineered mice to study molecular pathways that mediate the role of UVB in making mice susceptible to papillomavirus-induced pathogenesis, as well as identifying cellular targets of MmuPV-1 encoded factors that mediate their role in pathogenesis.

Materials and Methods

Animals

The following mice were obtained and bred for the purpose of this study (vendor in parenthesis): immunocompetent FVB/NJ (Taconic), BALB/c (Charles River), C57BL/6 (Jacksons Lab); immunodeficient athymic BALB/c FoxN1nu/nu (Harlan). All infected mice were housed in aseptic conditions in micro-isolator cages. Animals were handled only by designated personnel and personal protection gear was changed between cages to prevent any cross contamination from virus.

Ethics statement

Mice were housed at McArdle Laboratory Animal Care Unit in strict accordance with guidelines approved by the Association for Assessment of Laboratory Animal Care, at the University of Wisconsin Medical School. All protocols for animal work were approved by the University of Wisconsin Medical School Institutional Animal Care and Use Committee (Protocol number: M02478).

MmuPV1 virus stock

MmuPV1 virus stock was generated by isolating MmuPV1 virions from papillomas in nude mice as described previously [18,19,60]. The MmuPV1 infection model was established in nude mice using quasivirions generated (as described previously [19,61]) using a clone of MmuPV1 obtained in the Lambert lab. Briefly, this clone of MmuPV1 was made by performing rolling circle amplification on the virus extracts (generously provided by Dr. Aravind Ingle, ACTREC, India) from the original colony of nude mice infected with MmuPV1 [15] followed by cloning into the pUC19 vector [15,16]. To encapsidate MmuPV1 genome we used pMusSHELL, a Mammalian expression plasmid with codon modified L1 and L2 genes of MmuPV1 (generously provided by Dr. Chris Buck, NIH) [60]. To confirm that the virus stock was infectious, nude mice were infected in parallel during each experiment as positive controls.

Isolation of infectious MmuPV1 virions

This protocol has been modified from previously described methods of isolating non-enveloped viruses from human skin or tissue samples [62,63]. Animals with warts were euthanized and 10 mg of excised wart was homogenized in 700 μl PBS containing Triton-X-100 (1%). Benzonase was added to the homogenized wart sample and incubated at 37°C for 30 minutes. Collagenase H (2 mg) was added, and the sample was vortexed and then incubated at 4°C overnight. Sodium chloride concentration was adjusted to 0.8M, the sample was centrifuged for 5 minutes at 5000g, and the supernatant was clarified by ultracentrifugation through an Optiprep (iodixanol) step gradient followed by fractionation as described in detail on the following website: http://home.ccr.cancer.gov/LCO/pseudovirusproduction.htm. Viral genome equivalence was estimated by comparing the amount of encapsidated viral DNA in the viral stock, liberated by treatment with proteinase K, to known standards of cloned MmuPV-1 genome by Southern analysis using MmuPV1-specific probes, followed by quantification using ImageJ software (S1 Text).

MmuPV1 infection model

In vivo infection with purified MmuPV1 virions was performed on scarified skin of the animals' ears. Animals were anesthetized and tails or inner ears were scarified using a 27-gauge syringe needle to scrape the epithelia (not sufficient to cause bleeding) followed by pipette delivery of virus solution using a siliconized pipette tip. This method is modified from a previously reported infection model [60]. As controls mice were mock infected with vehicle i.e. optiprep as described above.

UV treatment

Animals were exposed to a single dose of UVB using a custom designed Research Irradiation Unit (Daavlin, Bryan, OH) [64–66]. This irradiation unit consists of an exposure unit mounted on fixed legs. Within the exposure unit there are four UVA and UVB lamps controlled using Daavlin Flex Control Integrating Dosimeters. In this system, dose units can be entered in milli-Joules per Centimeter Square for UVB (mJ/cm2) and Joules per Centimeter Square for UVA (J/cm2); variations in energy output are automatically compensated to deliver the desired dose. This enables us to expose the animals to an accurate dosimetery of UVB radiation. For accuracy, the machine is periodically calibrated using International Light IL 1400, digital light meter (Daavlin Company). For ear shielding experiments, ear sites were shielded from UV exposure by covering the head of anesthetized mice with tin foil during UV exposure.

Papilloma scoring and size determination

Papillomas were measured bi-weekly. Fraction of papillomas that completely regressed was computed by expressing number of papillomas at end-point of study (i.e. 6-months post infection) compared to maximum number of papillomas formed (i.e. papillomas at 3-months post infection). Papilloma size was determined by calculating the cubic root of the product of length × width × height to obtain a geometric mean diameter (GMD) as described previously [67]. Data were represented as the means (± standard errors of the mean [SEM] of the GMDs for each test group.

Delayed Type Hypersensitivity assays

DTH assays were performed by topically applying 1-Chloro-2,4-Di-Nitrobenzene (DNCB) (Sigma) using a modification of a method described previously [41]. For short-term immunosuppression studies, mice irradiated with 300mJ/cm2 UVB. Ten days post UVB-irradiation mice were shaved on their backs and sensitized by topically applying 0.5% DNCB in 50μl vehicle (4:1 acetone:olive oil). Five days post-sensitization they were challenged with DNCB in the left ear and ear thickness was measured. To assess correlation between wart incidence and long-term immunosuppression mice were either infected with MmuPV1 or mock infected in the right ear followed by treatment with 300mJ/cm2 UVB. Mice were sensitized with DNCB (0.3% in 50μl vehicle) 10 days post UVB exposure on their backs. Three months post-sensitization, mice were challenged with DNCB in the left ear and ear thickness was measured. Ear thickness was measured by means of Vernier calipers up to 96 hrs. The average of three readings taken at different points across the ear swelling was considered as a single measurement. Ear thickness was reported as the average of the difference between challenged and control ears 24–96 hours post challenge. The control (right) ear was the ear that was challenged with the vehicle only. The standard error in measurement was computed by measuring standard deviation.

Histology

Skin was harvested, fixed in 4% paraformaldehyde, and embedded in paraffin. Serial sections (5 μm thick) were stained with hematoxylin and eosin (H&E) and evaluated for histopathological features. Immunofluorescent staining was performed on sections after deparaffinizing with xylenes and rehydrating with graded ethanol, respectively. For cytokeratin14 staining, sections were blocked for 1hour at room temperature with goat serum followed by incubation with K14 antibody (Covance) at 1:1000 dilution for 1hour at room temperature. K14 signals were detected with Alexa-fluor 594 against rabbit. For L1-FISH dual immunofluorescent staining, first L1 immunofluorescent staining was performed. Antigen retrieval was performed using Proteinase K (20 μg/ml) for 15 minutes at 37°C. Samples were blocked for 1hour at room temperature with 5% goat serum and incubated overnight at 4°C with rabbit polyclonal immune serum directed against MusPV1 L1 [18,19] at 1:5000 dilution (Gift from Dr. Chris Buck, NCI). To proceed with MmuPV1 fluorescent in situ hybridization (FISH), the L1 stained tissue was dehydrated using a series of ice-cold ethanols (70%, 80%, 95%) for 2 minutes each. Slides were dried by placing them in an empty container at 50°C for 5 minutes and then placed in denaturation solution (28 mL formamide, 4 mL 20X SSC pH 5.3, 8 mL water) at 72°C for 2 minutes. The ethanol series was repeated, sections were dried, and denatured digoxigenin (DIG-11-dUTP, Roche)-labeled probe was hybridized to cells overnight at 37°C in a humidified chamber. To make the probe, nick translation was used to label MmuPV1-plasmid DNA with digoxigenin. After washing with 2X SSC and 50% formamide at 50°C (for 30 minutes twice) and 2X SSC at 50°C (for 30 minutes twice), MmuPV1 FISH DNA signals were detected with a digoxigenin-specific antibody conjugated to fluorescein isothiocyanate (Sigma, F3523) at 2% by volume in STM solution (4X SSC, 5% non-fat dried milk, 0.05% Tween-20, 0.002% sodium azide) for 30 minutes at 37°C. L1 protein signals were detected using Alexa-fluor 488 antibody against rabbit. Nuclei were counterstained with DAPI. All images were captured using a Zeiss AxioImager M2 microscope and AxioVision software version 4.8.2 (Jena, Germany).

Statistical analysis

All statistical analyses were performed using MSTAT statistical software version 6.1.4 (http://www.mcardle.wisc.edu/mstat).

Supporting Information

S1 Fig.

(A) Papillomas on ears of immunocompetent FVB/NJ mice. Examples of papillomas that arose on ears of FVB/NJ mice infected with 108 VGE MmuPV1 followed by irradiation with 300mJ/cm2 UVB. These images represent papillomas at 6 months post infection. (B) Histopathology of MmuPV1-induced ear papillomas in UVB irradiated FVB/NJ mice. Top panel shows H&E images (taken using a 2.5X objective) of lesions. Bottom panel with insets (taken using a 20X objective) indicate corresponding areas of lesions showing focal regions invasivity within the underlying dermis. Scale bars denote 100μm.

https://doi.org/10.1371/journal.ppat.1005664.s001

(TIF)

S2 Fig. Quantification of Viral Genome Equivalence by Southern Hybridization.

Southern analysis used to quantify virus stocks obtained from MmuPV1 induced ear warts in immunodeficient nude mice and UVB-irradiated FVB/NJ mice.

https://doi.org/10.1371/journal.ppat.1005664.s002

(TIF)

S3 Fig. Fluorescent in situ hybridization to detect MmuPV1 in papilloma and malignant regions.

MmuPV1-FISH was performed on papillomas that arose on ears of FVB/NJ mice infected with 108 VGE MmuPV1 followed by irradiation with 300mJ/cm2 UVB. Areas of papillomatosis show presence of amplified viral DNA (red) whereas cancerous regions show little to no presence of amplivied viral DNA. Nuclei were counterstained with DAPI (blue). Scale bars denote 100μm.

https://doi.org/10.1371/journal.ppat.1005664.s003

(TIF)

S1 Table. Papilloma incidence in UVB-irradiated FVB/NJ mice infected with MmuPV1 at different time points.

FVB/NJ mice were either treated with UVB twenty-four hours prior to infection (Group 1), twenty-four hours post-infection (Group 2) or fourteen days post-infection (Group 3) with MmuPV1. All animals were infected with 108 VGE MmuPV1. UVB dose was 300mJ/cm2. Data shown here represents sites scored at 3 months post-infection.

https://doi.org/10.1371/journal.ppat.1005664.s004

(TIF)

S2 Table. Papilloma incidence in immunocompetent mice infected with MmuPV1 followed by UVA or UVB irradiation.

FVB/NJ mice were either treated with UVA alone (Group 1), UVB alone (Group 2), or both UVA and UVB (Group 3) twenty-four hours post-infection with 108 VGE MmuPV1. UVB dose was 300mJ/cm2 and UVA dose was 300J/ cm2. Data shown here represents sites scored at 3 months post-infection.

https://doi.org/10.1371/journal.ppat.1005664.s005

(TIF)

S3 Table. Papilloma incidence in different strains of immunocompetent mice infected with MmuPV1 followed by UVB irradiation.

Mice of different genetic backgrounds were infected with 108 VGE MmuPV1 and irradiated with designated dose of UVB twenty-four hours post-infection. Data shown here represents sites scored at 3 months post-infection.

https://doi.org/10.1371/journal.ppat.1005664.s006

(TIF)

S1 Text. Detailed protocol for Southern analysis.

Text accompanying materials and methods.

https://doi.org/10.1371/journal.ppat.1005664.s007

(PDF)

Acknowledgments

We thank Dr. Chris Buck (National Cancer Institute, Bethesda, MD) for providing us with antibody against L1 and MusPVSheLL plasmid. We thank Dr. Aravind Ingle (Advanced Centre for Treatment and Research in Cancer, Tata Memorial Hospital, India) for providing DNA extracts from the originally observed warts arising in their colony of Nude mice. We thank Amy Liem and Dr. Haoshun Huang (University of Wisconsin-Madison) for assisting in animal experiments. We thank the UWCC histology facility for embedding tissue, and Harlene Edwards (UW-Madison) for sectioning tissue. We thank Drs. Jack Longley and Ruth Sullivan (UW-Madison) for providing valuable guidance for histopathology. We thank Norman Drinkwater (UW-Madison) for helpful advice in performing statistical analysis. We thank Dr. Karl Munger and Jordan Meyers (Tufts University) for providing valuable guidance and discussions during this study. Finally, we thank Dr. Megan Spurgeon, Kathleen Makielski and Patrick Nyman (University of Wisconsin-Madison) for helpful discussions and critical reading of the manuscript.

Author Contributions

Conceived and designed the experiments: AU IHF PFL. Performed the experiments: AU SY HCP. Analyzed the data: AU IHF PFL. Wrote the paper: AU PFL.

References

  1. 1. de Villiers EM, Fauquet C, Broker TR, Bernard HU, zur Hausen H (2004) Classification of papillomaviruses. Virology 324: 17–27. pmid:15183049
  2. 2. McLaughlin-Drubin ME, Meyers J, Munger K (2012) Cancer associated human papillomaviruses. Curr Opin Virol 2: 459–466. pmid:22658985
  3. 3. Jablonska S, Majewski S, Obalek S, Orth G (1997) Cutaneous warts. Clin Dermatol 15: 309–319. pmid:9255438
  4. 4. Jablonska S, Orth G, Obalek S, Croissant O (1985) Cutaneous warts. Clinical, histologic, and virologic correlations. Clin Dermatol 3: 71–82. pmid:2850861
  5. 5. Tschandl P, Rosendahl C, Kittler H (2014) Cutaneous human papillomavirus infection: manifestations and diagnosis. Curr Probl Dermatol 45: 92–97. pmid:24643180
  6. 6. Lebwohl MG, Rosen T, Stockfleth E (2010) The role of human papillomavirus in common skin conditions: current viewpoints and therapeutic options. Cutis 86: suppl 1–11; quiz suppl 12.
  7. 7. Mammas IN, Spandidos DA, Sourvinos G (2014) Genomic diversity of human papillomaviruses (HPV) and clinical implications: an overview in adulthood and childhood. Infect Genet Evol 21: 220–226. pmid:24291228
  8. 8. van Haalen FM, Bruggink SC, Gussekloo J, Assendelft WJ, Eekhof JA (2009) Warts in primary schoolchildren: prevalence and relation with environmental factors. Br J Dermatol 161: 148–152. pmid:19438464
  9. 9. Akgul B, Cooke JC, Storey A (2006) HPV-associated skin disease. J Pathol 208: 165–175. pmid:16362995
  10. 10. Astori G, Lavergne D, Benton C, Hockmayr B, Egawa K, et al. (1998) Human papillomaviruses are commonly found in normal skin of immunocompetent hosts. J Invest Dermatol 110: 752–755. pmid:9579540
  11. 11. Gassenmaier A, Fuchs P, Schell H, Pfister H (1986) Papillomavirus DNA in warts of immunosuppressed renal allograft recipients. Arch Dermatol Res 278: 219–223. pmid:3015050
  12. 12. Antonsson A, Forslund O, Ekberg H, Sterner G, Hansson BG (2000) The ubiquity and impressive genomic diversity of human skin papillomaviruses suggest a commensalic nature of these viruses. J Virol 74: 11636–11641. pmid:11090162
  13. 13. Howley PM, Pfister HJ (2015) Beta genus papillomaviruses and skin cancer. Virology 479–480: 290–296.
  14. 14. Connolly K, Manders P, Earls P, Epstein RJ (2014) Papillomavirus-associated squamous skin cancers following transplant immunosuppression: one Notch closer to control. Cancer Treat Rev 40: 205–214. pmid:24051018
  15. 15. Ingle A, Ghim S, Joh J, Chepkoech I, Bennett Jenson A, et al. (2011) Novel laboratory mouse papillomavirus (MusPV) infection. Vet Pathol 48: 500–505. pmid:20685915
  16. 16. Joh J, Jenson AB, King W, Proctor M, Ingle A, et al. (2011) Genomic analysis of the first laboratory-mouse papillomavirus. J Gen Virol 92: 692–698. pmid:21084500
  17. 17. Cladel NM, Budgeon LR, Cooper TK, Balogh KK, Hu J, et al. (2013) Secondary infections, expanded tissue tropism, and evidence for malignant potential in immunocompromised mice infected with Mus musculus papillomavirus 1 DNA and virus. J Virol 87: 9391–9395. pmid:23785210
  18. 18. Handisurya A, Day PM, Thompson CD, Bonelli M, Lowy DR, et al. (2014) Strain-specific properties and T cells regulate the susceptibility to papilloma induction by Mus musculus papillomavirus 1. PLoS Pathog 10: e1004314. pmid:25121947
  19. 19. Handisurya A, Day PM, Thompson CD, Buck CB, Pang YY, et al. (2013) Characterization of Mus musculus papillomavirus 1 infection in situ reveals an unusual pattern of late gene expression and capsid protein localization. J Virol 87: 13214–13225. pmid:24067981
  20. 20. Sundberg JP, Stearns TM, Joh J, Proctor M, Ingle A, et al. (2014) Immune status, strain background, and anatomic site of inoculation affect mouse papillomavirus (MmuPV1) induction of exophytic papillomas or endophytic trichoblastomas. PLoS One 9: e113582. pmid:25474466
  21. 21. Cladel NM, Budgeon LR, Balogh KK, Cooper TK, Hu J, et al. (2015) A novel pre-clinical murine model to study the life cycle and progression of cervical and anal papillomavirus infections. PLoS One 10: e0120128. pmid:25803616
  22. 22. Chen AC, McMillan NA, Antonsson A (2008) Human papillomavirus type spectrum in normal skin of individuals with or without a history of frequent sun exposure. J Gen Virol 89: 2891–2897. pmid:18931088
  23. 23. Hampras SS, Giuliano AR, Lin HY, Fisher KJ, Abrahamsen ME, et al. (2014) Natural history of cutaneous human papillomavirus (HPV) infection in men: the HIM study. PLoS One 9: e104843. pmid:25198694
  24. 24. Iannacone MR, Michael KM, Giuliano AR, Waterboer T, Pawlita M, et al. (2010) Risk factors for cutaneous human papillomavirus seroreactivity among patients undergoing skin cancer screening in Florida. J Infect Dis 201: 760–769. pmid:20105078
  25. 25. Amtmann E, Volm M, Wayss K (1984) Tumour induction in the rodent Mastomys natalensis by activation of endogenous papilloma virus genomes. Nature 308: 291–292. pmid:6322013
  26. 26. Tilbrook PA, Greenoak GE, Reeve VE, Canfield PJ, Gissmann L, et al. (1989) Identification of papillomaviral DNA sequences in hairless mouse tumours induced by ultraviolet irradiation. J Gen Virol 70 (Pt 4): 1005–1009. pmid:2543760
  27. 27. Zhang P, Nouri M, Brandsma JL, Iftner T, Steinberg BM (1999) Induction of E6/E7 expression in cottontail rabbit papillomavirus latency following UV activation. Virology 263: 388–394. pmid:10544111
  28. 28. Hennings H, Glick AB, Lowry DT, Krsmanovic LS, Sly LM, et al. (1993) FVB/N mice: an inbred strain sensitive to the chemical induction of squamous cell carcinomas in the skin. Carcinogenesis 14: 2353–2358. pmid:8242866
  29. 29. Egawa N, Egawa K, Griffin H, Doorbar J (2015) Human Papillomaviruses; Epithelial Tropisms, and the Development of Neoplasia. Viruses 7: 3863–3890. pmid:26193301
  30. 30. Kripke ML (2013) Reflections on the field of photoimmunology. J Invest Dermatol 133: 27–30. pmid:22854621
  31. 31. Fisher MS, Kripke ML (1977) Systemic alteration induced in mice by ultraviolet light irradiation and its relationship to ultraviolet carcinogenesis. Proc Natl Acad Sci U S A 74: 1688–1692. pmid:300876
  32. 32. Ullrich SE, Byrne SN (2012) The immunologic revolution: photoimmunology. J Invest Dermatol 132: 896–905. pmid:22170491
  33. 33. DeFabo EC, Kripke ML (1979) Dose-response characteristics of immunologic unresponsiveness to UV-induced tumors produced by UV irradiation of mice. Photochem Photobiol 30: 385–390. pmid:515194
  34. 34. De Fabo EC, Kripke ML (1980) Wavelength dependence and dose-rate independence of UV radiation-induced immunologic unresponsiveness of mice to a UV-induced fibrosarcoma. Photochem Photobiol 32: 183–188. pmid:7433530
  35. 35. Kripke ML (1984) Immunological unresponsiveness induced by ultraviolet radiation. Immunol Rev 80: 87–102. pmid:6237978
  36. 36. Schwarz A, Schwarz T (2002) Molecular determinants of UV-induced immunosuppression. Exp Dermatol 11 Suppl 1: 9–12. pmid:12444952
  37. 37. Kripke ML, Thorn RM, Lill PH, Civin CI, Pazmino NH, et al. (1979) Further characterization of immunological unresponsiveness induced in mice by ultraviolet radiation. Growth and induction of nonultraviolet-induced tumors in ultraviolet-irradiated mice. Transplantation 28: 212–217. pmid:91250
  38. 38. Noonan FP, De Fabo EC, Kripke ML (1981) Suppression of contact hypersensitivity by UV radiation and its relationship to UV-induced suppression of tumor immunity. Photochem Photobiol 34: 683–689. pmid:6458828
  39. 39. Kripke ML, Cox PA, Alas LG, Yarosh DB (1992) Pyrimidine dimers in DNA initiate systemic immunosuppression in UV-irradiated mice. Proc Natl Acad Sci U S A 89: 7516–7520. pmid:1502162
  40. 40. Mottram PL, Mirisklavos A, Clunie GJ, Noonan FP (1988) A single dose of UV radiation suppresses delayed type hypersensitivity responses to alloantigens and prolongs heart allograft survival in mice. Immunol Cell Biol 66 (Pt 5–6): 377–385. pmid:2976024
  41. 41. Dunn LA, Evander M, Tindle RW, Bulloch AL, de Kluyver RL, et al. (1997) Presentation of the HPV16E7 protein by skin grafts is insufficient to allow graft rejection in an E7-primed animal. Virology 235: 94–103. pmid:9300040
  42. 42. Frazer IH, Leippe DM, Dunn LA, Liem A, Tindle RW, et al. (1995) Immunological responses in human papillomavirus 16 E6/E7-transgenic mice to E7 protein correlate with the presence of skin disease. Cancer Res 55: 2635–2639. pmid:7540107
  43. 43. Wang JW, Jiang R, Peng S, Chang YN, Hung CF, et al. (2015) Immunologic Control of Mus musculus Papillomavirus Type 1. PLoS Pathog 11: e1005243. pmid:26495972
  44. 44. Nesnow S, Bergman H, Slaga TJ (1986) Comparison of the tumorigenic response of SENCAR and C57BL/6 mice to benzo(a)pyrene and the inter-experimental variability over a three-year period. Environ Health Perspect 68: 19–25. pmid:3780628
  45. 45. Reeve VE, Greenoak GE, Canfield PJ, Boehm-Wilcox C, Tilbrook PA, et al. (1989) Enhancement of u.v.-induced skin carcinogenesis in the hairless mouse by inoculation with cell-free extracts of skin tumours. Immunol Cell Biol 67 (Pt 6): 421–427. pmid:2560465
  46. 46. Forslund O, Lindelof B, Hradil E, Nordin P, Stenquist B, et al. (2004) High prevalence of cutaneous human papillomavirus DNA on the top of skin tumors but not in "Stripped" biopsies from the same tumors. J Invest Dermatol 123: 388–394. pmid:15245440
  47. 47. Hufbauer M, Lazic D, Akgul B, Brandsma JL, Pfister H, et al. (2010) Enhanced human papillomavirus type 8 oncogene expression levels are crucial for skin tumorigenesis in transgenic mice. Virology 403: 128–136. pmid:20451232
  48. 48. Maglennon GA, McIntosh PB, Doorbar J (2014) Immunosuppression facilitates the reactivation of latent papillomavirus infections. J Virol 88: 710–716. pmid:24173230
  49. 49. Joh J, Jenson AB, Proctor M, Ingle A, Silva KA, et al. (2012) Molecular diagnosis of a laboratory mouse papillomavirus (MusPV). Exp Mol Pathol 93: 416–421. pmid:22796029
  50. 50. Herber R, Liem A, Pitot H, Lambert PF (1996) Squamous epithelial hyperplasia and carcinoma in mice transgenic for the human papillomavirus type 16 E7 oncogene. J Virol 70: 1873–1881. pmid:8627712
  51. 51. Riley RR, Duensing S, Brake T, Munger K, Lambert PF, et al. (2003) Dissection of human papillomavirus E6 and E7 function in transgenic mouse models of cervical carcinogenesis. Cancer Res 63: 4862–4871. pmid:12941807
  52. 52. Song S, Pitot HC, Lambert PF (1999) The human papillomavirus type 16 E6 gene alone is sufficient to induce carcinomas in transgenic animals. J Virol 73: 5887–5893. pmid:10364340
  53. 53. Dong W, Kloz U, Accardi R, Caldeira S, Tong WM, et al. (2005) Skin hyperproliferation and susceptibility to chemical carcinogenesis in transgenic mice expressing E6 and E7 of human papillomavirus type 38. J Virol 79: 14899–14908. pmid:16282489
  54. 54. Paiva I, Gil da Costa RM, Ribeiro J, Sousa H, Bastos M, et al. (2015) A role for microRNA-155 expression in microenvironment associated to HPV-induced carcinogenesis in K14-HPV16 transgenic mice. PLoS One 10: e0116868. pmid:25625305
  55. 55. Hashimoto Y, Arai I, Nakanishi Y, Sakurai T, Nakamura A, et al. (2004) Scratching of their skin by NC/Nga mice leads to development of dermatitis. Life Sci 76: 783–794. pmid:15581910
  56. 56. Streilein JW, Bergstresser PR (1988) Genetic basis of ultraviolet-B effects on contact hypersensitivity. Immunogenetics 27: 252–258. pmid:3346042
  57. 57. Poon TS, Barnetson RS, Halliday GM (2005) Sunlight-induced immunosuppression in humans is initially because of UVB, then UVA, followed by interactive effects. J Invest Dermatol 125: 840–846. pmid:16185286
  58. 58. Byrne SN, Spinks N, Halliday GM (2002) Ultraviolet a irradiation of C57BL/6 mice suppresses systemic contact hypersensitivity or enhances secondary immunity depending on dose. J Invest Dermatol 119: 858–864. pmid:12406331
  59. 59. Rana S, Byrne SN, MacDonald LJ, Chan CY, Halliday GM (2008) Ultraviolet B suppresses immunity by inhibiting effector and memory T cells. Am J Pathol 172: 993–1004. pmid:18292235
  60. 60. Handisurya A, Day PM, Thompson CD, Buck CB, Kwak K, et al. (2012) Murine skin and vaginal mucosa are similarly susceptible to infection by pseudovirions of different papillomavirus classifications and species. Virology 433: 385–394. pmid:22985477
  61. 61. Pyeon D, Lambert PF, Ahlquist P (2005) Production of infectious human papillomavirus independently of viral replication and epithelial cell differentiation. Proc Natl Acad Sci U S A 102: 9311–9316. pmid:15958530
  62. 62. Schowalter RM, Pastrana DV, Pumphrey KA, Moyer AL, Buck CB (2010) Merkel cell polyomavirus and two previously unknown polyomaviruses are chronically shed from human skin. Cell Host Microbe 7: 509–515. pmid:20542254
  63. 63. Peretti A, FitzGerald PC, Bliskovsky V, Buck CB, Pastrana DV (2015) Hamburger polyomaviruses. J Gen Virol 96: 833–839. pmid:25568187
  64. 64. Ahsan H, Reagan-Shaw S, Eggert DM, Tan TC, Afaq F, et al. (2007) Protective effect of sanguinarine on ultraviolet B-mediated damages in SKH-1 hairless mouse skin: implications for prevention of skin cancer. Photochem Photobiol 83: 986–993. pmid:17645675
  65. 65. Calvo-Castro L, Syed DN, Chamcheu JC, Vilela FM, Perez AM, et al. (2013) Protective effect of tropical highland blackberry juice (Rubus adenotrichos Schltdl.) against UVB-mediated damage in human epidermal keratinocytes and in a reconstituted skin equivalent model. Photochem Photobiol 89: 1199–1207. pmid:23711186
  66. 66. Afaq F, Syed DN, Malik A, Hadi N, Sarfaraz S, et al. (2007) Delphinidin, an anthocyanidin in pigmented fruits and vegetables, protects human HaCaT keratinocytes and mouse skin against UVB-mediated oxidative stress and apoptosis. J Invest Dermatol 127: 222–232. pmid:16902416
  67. 67. Cladel NM, Hu J, Balogh K, Mejia A, Christensen ND (2008) Wounding prior to challenge substantially improves infectivity of cottontail rabbit papillomavirus and allows for standardization of infection. J Virol Methods 148: 34–39. pmid:18061687