Advertisement
Original Article| Volume 135, ISSUE 5, P1275-1282, May 2015

The Presence of Betapapillomavirus Antibodies around Transplantation Predicts the Development of Keratinocyte Carcinoma in Organ Transplant Recipients: A Cohort Study

      Organ transplant recipients (OTRs) have an increased risk of developing keratinocyte carcinomas (KCs). The aim of this study was to correlate infection with human papillomaviruses (HPVs) belonging to the beta genus (Beta-papillomavirus (Beta-PV)) at transplantation with later development of KCs. In a cohort study, sera collected between 1 year before and 1 year after transplantation of OTRs transplanted between 1990 and 2006 were tested for antibody responses against the L1 capsid antigen of Beta-PV and other HPV genera (Gamma-, Mu-, Nu-, and Alpha-PV) using multiplex serology. The OTRs were followed for a maximum of 22 years. Cox regression models with KC, squamous cell carcinoma (SCC), and basal cell carcinoma (BCC) as outcome variables were used. Out of 445 OTRs, 60 had developed KC: 14 developed only SCC, 24 only BCC, and 22 both types of KC. The time-dependent hazard ratio (HR) to develop either or both types of KC, adjusted for age, sex, and transplanted organ, in tested Beta-PV-seropositive OTR around the time of transplantation compared with Beta-PV-seronegative OTR was 2.9 (95% confidence interval (CI) 1.3–6.4). The HR for SCC was 2.9 (95% CI 0.99–8.5) and for BCC it was 3.1 (95% CI 1.2–8.0). There was also an association between Mu-PV seropositivity and KC, but there were no significant associations between other HPV genera tested and KC. A positive seroresponse for Beta-PV around transplantation significantly predicted the development of KC in OTRs up to 22 years later, providing additional evidence that infection with Beta-PV has a role in KC carcinogenesis.

      Abbreviations

      BCC
      basal cell carcinoma
      CI
      confidence interval
      HPV
      human papillomavirus
      HR
      hazard ratio
      KC
      keratinocyte carcinoma
      LUMC
      Leiden University Medical Center
      OR
      odds ratio
      OTR
      organ transplant recipient
      PV
      papillomavirus
      SCC
      squamous cell carcinoma

      INTRODUCTION

      Keratinocyte carcinomas (KCs) are the most prevalent malignancies seen in organ transplant recipients (OTRs). The incidence of squamous cell carcinoma (SCC) is 60–250 times increased compared with the general population, and for basal cell carcinomas (BCCs) this is 10–40 times (
      • Hartevelt M.M.
      • Bouwes Bavinck J.N.
      • Kootte A.M.
      • et al.
      Incidence of skin cancer after renal transplantation in The Netherlands.
      ;
      • Euvrard S.
      • Kanitakis J.
      • Claudy A.
      Skin cancers after organ transplantation.
      ;
      • Moloney F.J.
      • Comber H.
      • O’Lorcain P.
      • et al.
      A population-based study of skin cancer incidence and prevalence in renal transplant recipients.
      ;
      • Tessari G.
      • Naldi L.
      • Boschiero L.
      • et al.
      Incidence and clinical predictors of a subsequent nonmelanoma skin cancer in solid organ transplant recipients with a first nonmelanoma skin cancer: a multicenter cohort study.
      ). In OTRs, SCCs are more frequently observed than BCC, with a SCC/BCC ratio of 5:1 compared with the general population where a SCC/BCC ratio of 1:4 is found (
      • Wisgerhof H.C.
      • van der Boog P.J.
      • de Fijter J.W.
      • et al.
      Increased risk of squamous-cell carcinoma in simultaneous pancreas kidney transplant recipients compared with kidney transplant recipients.
      ).
      Life-long immunosuppressive therapy is the most important risk factor for developing KCs in OTRs. Other important risk factors include sun exposure, male gender, older age, smoking, and fair skin with susceptibility to sunburn. These are similar risk factors as in the immunocompetent population (
      • Terhorst D.
      • Drecoll U.
      • Stockfleth E.
      • et al.
      Organ transplant recipients and skin cancer: assessment of risk factors with focus on sun exposure.
      ). A role of human papillomaviruses (HPVs) in the development of KC has also been frequently suggested (
      • Proby C.M.
      • Harwood C.A.
      • Neale R.E.
      • et al.
      A case-control study of betapapillomavirus infection and cutaneous squamous cell carcinoma in organ transplant recipients.
      ;
      • Andersson K.
      • Michael K.M.
      • Luostarinen T.
      • et al.
      Prospective study of human papillomavirus seropositivity and risk of nonmelanoma skin cancer.
      ;
      • Neale R.E.
      • Weissenborn S.
      • Abeni D.
      • et al.
      Human papillomavirus load in eyebrow hair follicles and risk of cutaneous squamous cell carcinoma.
      ).
      HPVs are small double-stranded DNA viruses that infect epithelia of the skin or mucosa. Over 150 different types have been described (
      • de Villiers E.M.
      • Fauquet C.
      • Broker T.R.
      • et al.
      Classification of papillomaviruses.
      ;
      • Bernard H.U.
      • Burk R.D.
      • Chen Z.
      • et al.
      Classification of papillomaviruses (PVs) based on 189 PV types and proposal of taxonomic amendments.
      ). HPV has been estimated to have a causative role in 5.2% of all human cancers (
      • Parkin D.M.
      The global health burden of infection-associated cancers in the year 2002.
      ). HPV types are classified in different classes (1–4) by the IARC (International Agency for Research on Cancer) based on their carcinogenic potential (
      • IARC
      Human papillomaviruses.
      ;
      • Bouvard V.
      • Baan R.
      • Straif K.
      • et al.
      A review of human carcinogens—Part B: biological agents.
      ).
      The first time that HPV infection was linked with KC was in patients with a rare autosomal recessive disease called epidermodysplasia verruciformis (
      • Jablonska S.
      • Dabrowski J.
      • Jakubowicz K.
      Epidermodysplasia verruciformis as a model in studies on the role of papovaviruses in oncogenesis.
      ). Patients with epidermodysplasia verruciformis have an increased susceptibility to widespread Beta-papillomavirus (Beta-PV) infections of the skin that progress to SCC in one-third of the patients, mainly on sun-exposed sites (
      • Lutzner M.A.
      Epidermodysplasia verruciformis. An autosomal recessive disease characterized by viral warts and skin cancer. A model for viral oncogenesis.
      ). In the years after transplantation, OTRs start to develop a clinical picture that resembles that of patients with epidermodysplasia verruciformis. Epidemiological and experimental data have advocated a potential carcinogenic role of cutaneous HPV infection in skin cancer, especially of SCC in the immunosuppressed population (
      • Bouwes Bavinck J.N.
      • Plasmeijer E.I.
      • Feltkamp M.C.
      Beta-papillomavirus infection and skin cancer.
      ).
      There are different methods to detect Beta-PV infection, based on measurement of viral protein, viral DNA, or serum antibody responses (
      • Aldabagh B.
      • Angeles J.G.
      • Cardones A.R.
      • et al.
      Cutaneous squamous cell carcinoma and human papillomavirus: is there an association?.
      ). Serological responses are considered as the most consistent marker of biologically relevant Beta-PV infection (
      • Andersson K.
      • Michael K.M.
      • Luostarinen T.
      • et al.
      Prospective study of human papillomavirus seropositivity and risk of nonmelanoma skin cancer.
      ), reflecting a past or present infection strong enough to evoke an immune response.
      Most studies investigating the association between serological responses to HPV and SCC are cross-sectional or case–control studies that were carried out in the immunocompetent population, and showed an association between detection of Beta-PV antibodies and SCC (
      • Feltkamp M.C.
      • Broer R.
      • di Summa F.M.
      • et al.
      Seroreactivity to epidermodysplasia verruciformis-related human papillomavirus types is associated with nonmelanoma skin cancer.
      ;
      • Karagas M.R.
      • Nelson H.H.
      • Sehr P.
      • et al.
      Human papillomavirus infection and incidence of squamous cell and basal cell carcinomas of the skin.
      ;
      • Waterboer T.
      • Abeni D.
      • Sampogna F.
      • et al.
      Serological association of beta and gamma human papillomaviruses with squamous cell carcinoma of the skin.
      ;
      • Karagas M.R.
      • Waterboer T.
      • Li Z.
      • et al.
      Genus beta human papillomaviruses and incidence of basal cell and squamous cell carcinomas of skin: population based case-control study.
      ;
      • Bzhalava D.
      • Guan P.
      • Franceschi S.
      • et al.
      A systematic review of the prevalence of mucosal and cutaneous human papillomavirus types.
      ). As far as we know, cohort studies investigating the association between serological responses to HPV and later development of a first KC were not carried out in OTR (
      • Aldabagh B.
      • Angeles J.G.
      • Cardones A.R.
      • et al.
      Cutaneous squamous cell carcinoma and human papillomavirus: is there an association?.
      ). Therefore, we designed a retrospective follow-up study to establish whether there is a relationship between the presence of serological responses to HPV, in particular Beta-PV, around transplantation and the development of KC (SCC and BCC) in the years after the organ transplantation.

      RESULTS

      The baseline characteristics of the 445 OTRs included in the study are shown in Table 1. The majority of patients were male (65%) and most patients (79%) underwent kidney transplantation. The mean age of the patients at transplantation was 47.3 years and the mean follow-up after transplantation was 11.8 year. A mean of 4.2 sera (range 1–17) was available from 1 year before to 1 year after transplantation. Among the patients who developed KC, the mean time from transplantation until the first KC was 7.7 years.
      Table 1Baseline characteristics of the organ transplant recipients who were included in this study
      All patients N(%)Keratinocyte carcinomaSquamous cell carcinoma, irrespective of basal cell carcinomaBasal cell carcinoma, irrespective of squamous cell carcinoma
      +-+-+-
      N(%)N(%)N(%)N(%)N(%)N(%)
      Transplantation type
       Kidney345 (78.0)41 (68.3)304 (79.0)22 (61.1)323 (79.0)30 (65.2)315 (78.9)
       Kidney and pancreas100 (22.0)19 (31.7)81 (21.0)14 (38.9)86 (21.0)16 (34.8)84 (21.1)
      Gender
       Male291 (65.4)40 (66.7)251 (65.2)28 (77.8)263 (64.3)29 (63.0)262 (65.7)
       Female154 (34.6)20 (33.3)134 (34.8)8 (22.2)146 (35.7)17 (37.0)137 (34.3)
      Age (years) at transplantation
       Mean–median (range)47.3–46.8 (19.9–76.4)47.5–47.7 (19.9–68.7)47.3–46.7 (21.2–76.4)50.0–50.0 (31.8–68.7)47.1–46.5 (19.9–76.4)45.4–46.1 (19.9–64.2)47.5–47.0 (21.2–76.4)
      Number sera
      Sera available and tested for human papillomavirus (HPV) antibody presence between 1 year before and until 1 year after transplantation.
       Mean–median (range)4.22–4 (1–17)5.87–7 (1–17)3.97–4 (1–10)6.92–8 (1–17)3.99–4 (1–10)5.76–6 (1–17)4.05–4 (1–10)
       Time from Tx to first KC
       Mean—median (range)7.7–6.9 (0.28–19.8)9.0–8.3 (0.78–18.6)8.0–6.6 (0.28–19.8)
      Abbreviations: KC, keratinocyte carcinoma; Tx, transplantation.
      1 Sera available and tested for human papillomavirus (HPV) antibody presence between 1 year before and until 1 year after transplantation.
      Supplementary Table S1 online shows the percentages of OTRs who were seropositive for any HPV type tested between 1 year before and 1 year after the transplantation, for all OTRs together and stratified for those with KC, SCC, and BCC. Table 2 summarizes Supplementary Table S1 online and provides the percentages of OTRs who were seropositive for at least one HPV type per HPV genus (Beta-, Gamma-, Mu-, Nu, or Alpha-PV). Seroresponses against the beta genus were more frequently observed in the patients with KC compared with patients without KC, reaching statistical significance (P=0.005; Table 2). Patients who were seropositive for at least one HPV type of a specific genus had a greater probability for being seropositive for another HPV genus (Supplementary Table S2 online).
      Table 2Numbers and percentages of organ transplant recipients with positive antibody responses against one or more human papillomavirus (HPV) types of different HPV genera
      HPV types tested were: 7 Alpha-HPV (HPVs 2, 3, 6, 7, 13, 16, and 27b), 16 Beta-HPV (HPVs 5, 8, 9, 15, 17, 20, 23, 24, 36, 38, 49, 75, 76, 92, 93, and 96), 8 Gamma-HPV (HPVs 4, 48, 50, 60, 65, 95, 101, and 103), 2 Mu-HPV (HPV1 and HPV63), and 1 Nu-HPV (HPV41).
      All patientsN(%)Keratinocyte carcinomaSquamous cell carcinoma, irrespective of basal cell carcinomaBasal cell carcinoma, irrespective of squamous cell carcinoma
      +-+-+-
      N(%)N(%)N(%)N(%)N(%)N(%)
      Any Beta-HPV
       Negative119 (26.7)7 (11.7)112 (29.1)4 (11.1)115 (28.1)5 (10.9)114 (28.6)
       Positive326 (73.3)53 (88.3)273 (70.9)32 (88.9)294 (71.9)41 (89.1)285 (71.4)
      Any Gamma-HPV
       Negative113 (25.4)12 (20.0)101 (26.2)6 (16.7)107 (26.2)10 (21.7)103 (25.8)
       Positive332 (74.6)48 (80.0)284 (73.8)30 (83.3)302 (73.8)36 (78.3)296 (74.2)
      Any Mu-HPV
       Negative151 (33.9)12 (20.0)139 (36.1)6 (16.7)145 (35.5)9 (19.6)142 (35.6)
       Positive294 (66.1)48 (80.0)246 (63.9)30 (83.3)264 (64.5)37 (80.4)257 (64.4)
      Any Nu-HPV
       Negative367 (82.5)46 (76.7)321 (83.4)27 (75.0)340 (83.1)37 (80.4)330 (82.7)
       Positive78 (17.5)14 (23.3)64 (16.6)9 (25.0)69 (16.9)9 (19.6)69 (17.3)
      Any Alpha-HPV
       Negative106 (23.8)13 (21.7)93 (24.2)8 (22.2)98 (24.0)10 (21.7)96 (24.1)
       Positive339 (76.2)47 (78.3)292 (75.8)28 (77.8)311 (76.0)36 (78.3)303 (75.9)
      1 HPV types tested were: 7 Alpha-HPV (HPVs 2, 3, 6, 7, 13, 16, and 27b), 16 Beta-HPV (HPVs 5, 8, 9, 15, 17, 20, 23, 24, 36, 38, 49, 75, 76, 92, 93, and 96), 8 Gamma-HPV (HPVs 4, 48, 50, 60, 65, 95, 101, and 103), 2 Mu-HPV (HPV1 and HPV63), and 1 Nu-HPV (HPV41).
      The cumulative incidence of KC in Beta-PV-seronegative and Beta-PV-seropositive OTR is depicted in Figure 1 and Supplementary Figure S1 online. Log-rank tests were performed and showed a significant difference between Beta-PV-seropositive and -seronegative OTRs (P=0.043). The Kaplan-Meier plots for SCC and BCC showed similar curves as the plots for KC (Supplementary Figures S2 and S3 online). Table 3 shows the time-dependent nonadjusted hazard ratios (HRs) and the HRs adjusted for age at transplantation, sex, and type of transplantation for the development of KC, SCC, and BCC. A positive association was found between seropositivity for at least one Beta-PV type detected between 1 year before and 1 year after the organ transplantation and the development of KC, SCC, and/or BCC. After adjustment, Mu-PV seropositivity was also associated with KC (Table 3). The Cox proportional hazard analyses (Table 3) and the Kaplan-Meier plots of cumulative incidence of KC, SCC, and/or BCC in HPV-seronegative and -seropositive OTRs for the other HPV genera tested (Gamma-, Nu-, and Alpha-PV) showed no significant association between seropositivity against any of these HPV genera (Supplementary Figures S1–3 online).
      Figure thumbnail gr1
      Figure 1Kaplan–Meier plot for development of keratinocyte carcinomas (KCs) in relation to human papillomavirus (HPV) serology. The cumulative incidence of keratinocyte carcinomas (KCs) is statistically significantly higher in organ transplant recipients (OTRs) who were seropositive against one or more HPV types of the Beta-papillomavirus (Beta-PV) genus between 1 year before and 1 year after the transplantation (P=0.043).
      Table 3Time-dependent Cox proportional hazard ratios for keratinocyte carcinoma, squamous cell carcinoma, and basal cell carcinoma in organ transplant recipients related to human papillomavirus (HPV) genus serostatus between 1 year before and 1 year after the transplantation
      Keratinocyte carcinomaSquamous cell carcinoma, irrespective of basal cell carcinomaBasal cell carcinoma, irrespective of squamous cell carcinoma
      Unadjusted hazard ratioAdjusted
      Adjusted for sex, age at transplantation, and type of transplantation.
      hazard ratio
      Unadjusted hazard ratioAdjusted
      Adjusted for sex, age at transplantation, and type of transplantation.
      hazard ratio
      Unadjusted hazard ratioAdjusted
      Adjusted for sex, age at transplantation, and type of transplantation.
      hazard ratio
      Beta-HPV2.2 (1.0–4.9)2.8 (1.3–6.4)2.1 (0.75–6.1)2.9 (0.99–8.5)2.3 (0.92–6.0)3.1 (1.2–8.0)
      Gamma-HPV1.1 (0.59–2.1)1.3 (0.66–2.4)1.2 (0.51–3.0)1.7 (0.68–4.1)0.97 (0.48–2.0)1.1 (0.54–2.2)
      Mu-HPV1.8 (0.95–3.4)2.0 (1.0–3.7)2.1 (0.86–5.0)2.4 (0.98–5.7)1.8 (0.86–3.7)1.9 (0.93–4.0)
      Nu-HPV1.2 (0.68–2.2)1.5 (0.80–2.7)1.3 (0.61–2.7)1.8 (0.84–4.0)1.0 (0.46–2.0)1.1 (0.53–2.3)
      Alpha-HPV0.93 (0.50–1.7)0.99 (0.53–1.8)0.84 (0.38–1.9)0.97 (0.44–2.2)0.92 (0.46–1.9)0.97 (0.48–2.0)
      1 Adjusted for sex, age at transplantation, and type of transplantation.
      Beta-PV species 1 (specifically HPV8 and HPV20) and Beta-PV species 3 (specifically HPV75 and HPV76) were most relevant for the positive association between Beta-PV seropositivity and the development of KC (Supplementary Figure S4 online). The potential importance of these Beta-PV types was confirmed for HPV8 and HPV75 by conditional step forward regression analyses, and a statistically significant negative association was observed for HPV96 (Table 4).
      Beta-PV species 1 (HPV8 and HPV20) and Beta-PV species 3 (HPV76) were also associated with the development of SCC (Supplementary Figure S5 online), confirmed by conditional step forward regression analyses (Table 4). A statistically significant negative association was found for HPV24 and HPV38 (Table 4).
      Table 4Association between 16 different Beta-papillomavirus (Beta-PV) types and skin cancer
      Keratinocyte carcinomaSquamous cell carcinomaBasal cell carcinoma
      All variables in the regression model
      Logistic regression model.
      Forward regression final modelAll variables in the regression modelForward regression final modelAll variables in the regression modelForward regression final model
      Species 1
       HPV51.3 (0.57; 3.1)1.7 (0.54; 5.1)1.2 (0.46; 3.1)
       HPV82.6 (1.2; 5.6)2.4 (1.3; 4.6)3.4 (1.3; 8.8)3.4 (1.4; 8.3)2.2 (0.96; 5.1)
       HPV202.6 (1.1; 5.9)2.9 (1.1; 8.0)2.6 (1.0; 6.6)2.6 (1.0; 6.6)
       HPV240.43 (0.15; 1.2)0.19 (0.05; 0.71)0.20 (0.07; 0.63)0.68 (0.22; 2.1)
       HPV360.50 (0.18; 1.4)0.61 (0.18; 2.1)0.40 (0.13; 1.3)
       HPV930.48 (0.17; 1.4)0.35 (0.09; 1.5)0.27 (0.07; 1.1)
      Species 2
       HPV90.37 (0.13; 0.99)0.50 (0.15; 1.7)0.25 (0.08; 0.81)0.29 (0.11; 0.76)
       HPV151.5 (0.60; 3.7)1.3 (0.38; 4.2)2.3 (0.86; 6.2)2.7 (1.2; 6.3)
       HPV170.82 (0.35; 1.9)0.72 (0.24)0.95 (0.37; 2.4)
       HPV231.3 (0.53; 3.3)3.1 (1.0; 9.8)2.4 (0.92; 6.4)1.4 (0.53; 4.0)
       HPV380.94 (0.38; 2.3)0.23 (0.06; 0.83)0.21 (0.07; 0.66)0.86 (0.31; 2.4)
      Species 3
       HPV491.4 (0.56; 3.6)1.9 (0.62; 5.9)1.4 (0.48; 3.9)
       HPV752.8 (1.1; 7.4)3.1 (1.5; 6.6)1.6 (0.47; 5.6)2.2 (0.74; 6.7)
       HPV761.8 (0.67; 4.8)2.3 (0.70; 7.8)2.6 (1.1; 6.5)1.9 (0.64; 5.6)2.9 (1.3; 6.5)
      Species 4
       HPV921.2 (0.49; 3.1)1.5 (0.45; 4.9)0.82 (0.28; 2.4)
      Species 5
       HPV960.41 (0.17; 1.0)0.38 (0.18; 0.83)0.56 (0.19; 1.7)0.63 (0.24; 1.7)
       Age
      Age: per year.
      1.0 (0.98; 1.0)1.0 (1.0; 1.1)1.1 (1.0; 1.1)0.99 (0.96; 1.0)
       Sex
      Sex: men compared with women.
      0.87 (0.47; 1.6)1.4 (0.56; 3.3)0.81 (0.40; 1.6)
       Organ
      Organ: combined pancreas and kidney compared with kidney.
      2.2 (1.1; 4.6)2.1 (1.1; 4.0)4.7 (1.8; 12.3)4.9 (2.0; 12.1)2.2 (1.0; 4.9)2.4 (1.2; 4.7)
      The bold numbers are statistically significant.
      1 Logistic regression model.
      2 Age: per year.
      3 Sex: men compared with women.
      4 Organ: combined pancreas and kidney compared with kidney.
      Beta-PV species 2 (HPV15) and species 3 (HPV75 and HPV76) were relevant for the development of BCC (Supplementary Figure S6 online), which was confirmed in the conditional step forward regression analyses for HPV15 and HPV26, whereas there was a statistically significant negative association with HPV9 (Table 4).

      DISCUSSION

      We have found a statistically significant association between the development of KC, SCC, and BCC after transplantation and antibody responses to one or more Beta-PV types, measured in serum drawn in a window between 1 year before and 1 year after organ transplantation. The antibody response against Beta-PV types was specific for the beta genus, as there were no significant associations for positive seroresponses against the gamma, nu, and alpha genera with the development of KC, SCC, or BCC. Interestingly, we also found a positive association between KC and Mu-PV seropositivity, an HPV genus not recognized previously as potentially associated with KC development.
      So far, most studies investigating HPV seroresponses in association with KC development concerned cross-sectional and case–control studies in the immunocompetent population (
      • Aldabagh B.
      • Angeles J.G.
      • Cardones A.R.
      • et al.
      Cutaneous squamous cell carcinoma and human papillomavirus: is there an association?.
      ). We are aware of only a few published cohort studies that were performed in the immunocompetent population. A Scandinavian study described a weak association for overall Beta-PV seropositivity before SCC diagnoses, whereas a significant association was observed for Beta-PV species 2 (odds ratio (OR)=1.3, 95% confidence interval (CI) 1.1–1.7), and for serum samples taken >18 years before diagnosis (OR=1.8, 95% CI 1.1–2.8) (
      • Andersson K.
      • Michael K.M.
      • Luostarinen T.
      • et al.
      Prospective study of human papillomavirus seropositivity and risk of nonmelanoma skin cancer.
      ). A population-based cohort study in Australia revealed no associations between the presence of any of the Beta-PV antibodies and the occurrence of SCC (
      • Plasmeijer E.I.
      • Pandeya N.
      • O’Rourke P.
      • et al.
      The association between cutaneous squamous cell carcinoma and betapapillomavirus seropositivity: a cohort study.
      ). However, among people who were <50 years old in 1992, the presence of Beta-PV antibodies was associated with a 2-fold increased risk of SCC (
      • Plasmeijer E.I.
      • Pandeya N.
      • O’Rourke P.
      • et al.
      The association between cutaneous squamous cell carcinoma and betapapillomavirus seropositivity: a cohort study.
      ). In the United Kingdom, a prospective pilot study found no significant differences in seropositivity of any of the 38 HPV types (
      • Casabonne D.
      • Michael K.M.
      • Waterboer T.
      • et al.
      A prospective pilot study of antibodies against human papillomaviruses and cutaneous squamous cell carcinoma nested in the Oxford component of the European Prospective Investigation into Cancer and Nutrition.
      ). One study investigated the possible role of seropositivity for different HPVs in the incidence of a subsequent KC in 107 patients with a follow-up of 5 years who were enrolled in a previous case–control study. HPV seropositivity at baseline was strongly associated with the risk of developing a second SCC after 5 years for a number of beta and gamma HPV types, but no association was found for developing a second BCC (
      • Paradisi A.
      • Waterboer T.
      • Sampogna F.
      • et al.
      Seropositivity for human papillomavirus and incidence of subsequent squamous cell and basal cell carcinomas of the skin in patients with a previous nonmelanoma skin cancer.
      ). As far as we know, our study is the first cohort study with a prolonged follow-up period in OTRs in which the relationship between HPV serology before and just after the organ transplantation was studied in relation to later development of a first KC.
      • Proby C.M.
      • Harwood C.A.
      • Neale R.E.
      • et al.
      A case-control study of betapapillomavirus infection and cutaneous squamous cell carcinoma in organ transplant recipients.
      conducted a case–control study in OTRs from The Netherlands, United Kingdom, France, and Italy to investigate the association between SCC development and seroresponses against Beta-PV combined with detection of Beta-PV DNA in plucked eyebrow hairs. A positive association was found for concordant DNA and seropositivity for HPV36 (OR 2.4; CI 1.0–5.4), with similar, but not statistically significant, associations for HPV5, HPV9, and HPV24. Concordant Beta-PV DNA presence in hair and a serological antibody response for at least one HPV type was significantly associated with SCC risk (OR 1.6; CI 1.1–2.5). Another case–control study conducted in immunocompetent patients in Florida, USA, revealed a positive association between seropositivity to increasing numbers of Beta-PV types and SCC (OR 1.9; 95% CI 1.2–3.0) (
      • Iannacone M.R.
      • Gheit T.
      • Waterboer T.
      • et al.
      Case-control study of cutaneous human papillomaviruses in squamous cell carcinoma of the skin.
      ). A similar association for Beta-PV types DNA in eyebrow hairs was found in another clinic-based case–control study (4 types vs. HPV-negative: OR 2.0, 95% CI 1.1–3.8) (
      • Iannacone M.R.
      • Gheit T.
      • Pfister H.
      • et al.
      Case-control study of genus-beta human papillomaviruses in plucked eyebrow hairs and cutaneous squamous cell carcinoma.
      ).
      • Waterboer T.
      • Abeni D.
      • Sampogna F.
      • et al.
      Serological association of beta and gamma human papillomaviruses with squamous cell carcinoma of the skin.
      performed a case–control study in an Italian immunocompetent population and found a significant 2- to 3-fold risk of SCC with seropositivity for any species 2 Beta-PV type (OR 3.3, 95% CI 1.2–8.7) and for the Gamma-PV types (OR 3.1, 95% CI 1.1–8.6). In a multinational case–control study in immunocompetent subjects, a positive antibody response against four or more Beta-PV types was associated with a doubled risk for SCC in Australia and the Netherlands (
      • Bouwes Bavinck J.N.
      • Neale R.E.
      • Abeni D.
      • et al.
      Multicenter study of the association between betapapillomavirus infection and cutaneous squamous cell carcinoma.
      ).
      Little data are published about the association between HPV infection and the development of BCC and the findings are contradictory. A relationship between seroresponses to Beta-PV and the development of BCC has been described previously (
      • Steger G.
      • Olszewsky M.
      • Stockfleth E.
      • et al.
      Prevalence of antibodies to human papillomavirus type 8 in human sera.
      ;
      • Stark S.
      • Petridis A.K.
      • Ghim S.J.
      • et al.
      Prevalence of antibodies against virus-like particles of Epidermodysplasia verruciformis-associated HPV8 in patients at risk of skin cancer.
      ;
      • Feltkamp M.C.
      • Broer R.
      • di Summa F.M.
      • et al.
      Seroreactivity to epidermodysplasia verruciformis-related human papillomavirus types is associated with nonmelanoma skin cancer.
      ;
      • Andersson K.
      • Michael K.M.
      • Luostarinen T.
      • et al.
      Prospective study of human papillomavirus seropositivity and risk of nonmelanoma skin cancer.
      ). Interestingly,
      • Paolini F.
      • Carbone A.
      • Benevolo M.
      • et al.
      Human papillomaviruses, p16INK4a and Akt expression in basal cell carcinoma.
      reported a high frequency (11.5%) of HPV15 DNA in a small cohort of BCC patients, in agreement with the observed positive association between HPV15 and BCC in our study. However, some other studies did not report a positive association between Beta-PV seroreactivity and BCC (
      • Karagas M.R.
      • Nelson H.H.
      • Sehr P.
      • et al.
      Human papillomavirus infection and incidence of squamous cell and basal cell carcinomas of the skin.
      ;
      • Andersson K.
      • Waterboer T.
      • Kirnbauer R.
      • et al.
      Seroreactivity to cutaneous human papillomaviruses among patients with nonmelanoma skin cancer or benign skin lesions.
      ;
      • Karagas M.R.
      • Waterboer T.
      • Li Z.
      • et al.
      Genus beta human papillomaviruses and incidence of basal cell and squamous cell carcinomas of skin: population based case-control study.
      ;
      • Iannacone M.R.
      • Gheit T.
      • Waterboer T.
      • et al.
      Case-control study of cutaneous human papillomavirus infection in basal cell carcinoma of the skin.
      ). A recent systematic review showed no association between Beta-PV serological responses and BCC (
      • Bzhalava D.
      • Guan P.
      • Franceschi S.
      • et al.
      A systematic review of the prevalence of mucosal and cutaneous human papillomavirus types.
      ). In a subgroup of patients in Florida, USA, with Beta-PV DNA-positive BCC, however, an association between a serological response against Beta-PV and BCC was observed (
      • Iannacone M.R.
      • Gheit T.
      • Waterboer T.
      • et al.
      Case-control study of cutaneous human papillomavirus infection in basal cell carcinoma of the skin.
      ). In a population-based case–control study of patients with SCC, BCC, and controls from New Hampshire, USA, it was found that seropositivity to Beta-PV (OR 1.5, 95% CI 1.0–2.1) was significantly associated with SCC. This risk increased with positivity for multiple HPV types. No associations were found with BCC risk (
      • Karagas M.R.
      • Nelson H.H.
      • Sehr P.
      • et al.
      Human papillomavirus infection and incidence of squamous cell and basal cell carcinomas of the skin.
      ).
      • Karagas M.R.
      • Waterboer T.
      • Li Z.
      • et al.
      Genus beta human papillomaviruses and incidence of basal cell and squamous cell carcinomas of skin: population based case-control study.
      confirmed these findings in a larger group of 2366 skin cancer cases and controls. We cannot exclude, however, that the OTRs who developed BCC in our study are at increased risk of developing SCC later on, because BCCs often precede the development of SCC in these patients.
      Besides the advantages, the retrospective design of our study also has some disadvantages. For instance, the sera we tested were not systematically collected with the purpose to study a relationship with KC but were randomly selected from routinely archived sera collected for HLA antibody screening and other diagnostic purposes. The HRs we calculated were adjusted for sex, age at transplantation, and type of transplantation (
      • Wisgerhof H.C.
      • van der Boog P.J.
      • de Fijter J.W.
      • et al.
      Increased risk of squamous-cell carcinoma in simultaneous pancreas kidney transplant recipients compared with kidney transplant recipients.
      ). We could not adjust for other risk factors for KC, like sun exposure, sunburns, skin type, and type of immunosuppression, because these data were not systematically collected for all patients. These risk factors are associated with an increased prevalence of HPV DNA in eyebrow hairs, but it is not known whether this is also the case for Beta-PV seroresponses (
      • Alotaibi L.
      • Provost N.
      • Gagnon S.
      • et al.
      Diversity of cutaneous human papillomavirus types in individuals with and without skin lesion.
      ;
      • Chen A.C.
      • McMillan N.A.
      • Antonsson A.
      Human papillomavirus type spectrum in normal skin of individuals with or without a history of frequent sun exposure.
      ;
      • de Koning M.N.
      • Weissenborn S.J.
      • Abeni D.
      • et al.
      Prevalence and associated factors of betapapillomavirus infections in individuals without cutaneous squamous cell carcinoma.
      ).
      A drawback of our study was that data about seroconversion were not systematically available. It can be speculated that seroconversion in the years following transplantation in OTRs who were seronegative around transplantation might have increased the risk of tumor development in the seronegative group that, if true, could have led to a stronger association between Beta-PV-type seropositivity and KC.
      In general, a drawback of studies reporting the association between serological responses and skin cancer is that because of differences in the employed methodologies (initially starting with western blotting, followed by HPV-type-specific ELISAs using L1 virus-like particles, and more recently by high-throughput multiplex fluorescent bead–based assays enabling the simultaneous detection of antibody responses in large series of serologic samples against a variety of HPVs) the outcomes between the different studies cannot always be reliably compared. In our study, however, this is only true for the earlier studies referred to (western blot by
      • Steger G.
      • Olszewsky M.
      • Stockfleth E.
      • et al.
      Prevalence of antibodies to human papillomavirus type 8 in human sera.
      and ELISA technique by
      • Stark S.
      • Petridis A.K.
      • Ghim S.J.
      • et al.
      Prevalence of antibodies against virus-like particles of Epidermodysplasia verruciformis-associated HPV8 in patients at risk of skin cancer.
      ;
      • Feltkamp M.C.
      • Broer R.
      • di Summa F.M.
      • et al.
      Seroreactivity to epidermodysplasia verruciformis-related human papillomavirus types is associated with nonmelanoma skin cancer.
      ), as the other studies all used the same multiplex serology technique that was performed in the same laboratory in Heidelberg. Other factors that may explain the different associations between HPV infection and KC are differences in immunosuppressive regimens between centers that may exert different effects on the present HPV types, different exposure to sunlight reflected by the different latitudes of the centers, and differences in exposure to different HPV types.
      The exact mechanism by which Beta-PV infection predisposes for the development of KC remains unclear. It has been hypothesized that the effects of the virus on UV-induced DNA damage withholds apoptosis and DNA repair by accumulation of mutations (
      • Bouwes Bavinck J.N.
      • Feltkamp M.C.
      Milk of human kindness?—HAMLET, human papillomavirus, and warts.
      ). This is supported by the fact that the E6 and E7 proteins from Beta-PV types are capable of inhibiting UV-induced cell cycle checkpoints and DNA repair mechanisms (
      • Jackson S.
      • Harwood C.
      • Thomas M.
      • et al.
      Role of Bak in UV-induced apoptosis in skin cancer and abrogation by HPV E6 proteins.
      ;
      • Jackson S.
      • Storey A.
      E6 proteins from diverse cutaneous HPV types inhibit apoptosis in response to UV damage.
      ;
      • Struijk L.
      • van der Meijden E.
      • Kazem S.
      • et al.
      Specific betapapillomaviruses associated with squamous cell carcinoma of the skin inhibit UVB-induced apoptosis of primary human keratinocytes.
      ;
      • Underbrink M.P.
      • Howie H.L.
      • Bedard K.M.
      • et al.
      E6 proteins from multiple human betapapillomavirus types degrade Bak and protect keratinocytes from apoptosis after UVB irradiation.
      ). This potentially deregulates cell cycle control and growth in keratinocytes (
      • zur Hausen H.
      Immortalization of human cells and their malignant conversion by high risk human papillomavirus genotypes.
      ;
      • Iannacone M.R.
      • Gheit T.
      • Waterboer T.
      • et al.
      Case-control study of cutaneous human papillomaviruses in squamous cell carcinoma of the skin.
      ). Some Beta-PV types could upregulate the p16INK4a and Akt/P13K pathway and might have a role in the carcinogenesis of BCC (
      • Paolini F.
      • Carbone A.
      • Benevolo M.
      • et al.
      Human papillomaviruses, p16INK4a and Akt expression in basal cell carcinoma.
      ). The epidemiological evidence that HPV is involved in BCC carcinogenesis, however, is much weaker compared with the possible involvement of HPV in cutaneous SCC carcinogenesis. In conclusion, a detectable seroresponse for one or more Beta-PV types around the time of transplantation predicts an increased risk for the development of KC in OTR. This study provides additional evidence that infection with Beta-PV types has a role in keratinocyte carcinogenesis. Additional studies will be necessary to confirm the observed association between Mu-PV infection and KC development.

      MATERIALS AND METHODS

      Study design

      Between 1966 and 2006, a total of 2,136 patients were transplanted in the Leiden University Medical Center (LUMC) (Supplementary Figure S2 online). OTRs consisted of both kidney transplant recipients and simultaneous pancreas–kidney transplant recipients.
      During the past 25 years, sera were collected from a considerable number of OTRs. The most extensive source of sera consisted of samples collected from 1989 for clinical care–related screening and diagnostic purposes that were systematically stored in the Department of Medical Microbiology of the LUMC. In the scope of earlier case–control studies, additional sera had been collected and stored (
      • Bouwes Bavinck J.N.
      • Vermeer B.J.
      • van der Woude F.J.
      • et al.
      Relation between skin cancer and HLA antigens in renal-transplant recipients.
      ;
      • Bouwes Bavinck J.N.
      • Gissmann L.
      • Claas F.H.
      • et al.
      Relation between skin cancer, humoral responses to human papillomaviruses, and HLA class II molecules in renal transplant recipients.
      ). From the 101 patients who participated in a prospective study (
      • Antonsson A.
      • Waterboer T.
      • Bouwes Bavinck J.N.
      • et al.
      Longitudinal study of seroprevalence and serostability of 34 human papillomavirus types in European organ transplant recipients.
      ), pretransplant sera were retraced from the Eurotransplant serumbank. Altogether, 7,912 sera were available from 1,269 OTRs (Supplementary Figure S2 online). Because of practical reasons we had to limit the number of sera to be tested to a maximum of 3,000, the number thought to provide enough statistical power to answer our questions.
      We selected 648 representative sera from all 164 OTRs who developed KC during follow-up for whom at least one serum sample was available. For each of these OTRs, we selected between 3 and 5 control OTRs without skin cancer. These patients were matched for type of transplantation (kidney transplant recipients or simultaneous pancreas–kidney transplant recipients), sex, and age at transplantation, resulting in the selection of 2,112 sera from 534 OTRs without skin cancer (Supplementary Figure S2 online). These sera had been collected between 7 years before and 42 years after the transplantation.
      As we were specifically interested in the association between seroresponses to HPV around the transplantation and the later development of KC, we further restricted our selection to OTRs who had sera collected between 1 year before and 1 year after the transplantation. This resulted in a final cohort consisting of 1,880 sera from 445 OTRs, of whom 60 had developed KC (14 only SCC, 24 only BCC, and 22 both types of skin cancer) during the 22-year follow-up period. The number of sera collected in these patients ranged from 1 to 11 (mean: 4.2, median: 4.0). The remaining sera will be analyzed in a separate study.
      Data about development of KC were retracted from the LUMC oncologic database, the LUMC pathology database, and by hand searching of the medical files. The study was approved by the medical ethical committee of the LUMC.

      Serum processing and analysis

      Serum samples were obtained from -20 °C freezers at the department of Medical Microbiology. They were defrosted and pipetted in 96-well plates. Some selected serum samples were excluded from further analysis, because either they could not be found, were not properly centrifuged initially, or contained too little volume (Supplementary Figure S7 online). The 96-well plates were shipped on dry ice to the German Cancer Research Center (DKFZ; Heidelberg, Germany) for analysis using Luminex technology (Austin, TX).
      This multiplex serology technique is an antibody detection method based on glutathione S-transferase capture ELISA in combination with fluorescent bead technology. Positive serology cutoff points are standardized at 200 mean fluorescence intensity targeted on antibodies against major capsid antigen L1 of the various HPV types (
      • Waterboer T.
      • Sehr P.
      • Michael K.M.
      • et al.
      Multiplex human papillomavirus serology based on in situ-purified glutathione s-transferase fusion proteins.
      ;
      • Michael K.M.
      • Waterboer T.
      • Sehr P.
      • et al.
      Seroprevalence of 34 human papillomavirus types in the German general population.
      ). Sera were tested for 7 Alpha-PV (HPVs 2, 3, 6, 7, 13, 16, and 27b), 16 Beta-PV (HPVs 5, 8, 9, 15, 17, 20, 23, 24, 36, 38, 49, 75, 76, 92, 93, and 96), 8 Gamma-PV (HPVs 4, 48, 50, 60, 65, 95, 101, and 103), 2 Mu-PV (HPV1 and HPV63), and one Nu-PV (HPV41).
      The OTRs were considered seropositive for specific HPV genera or specific HPV types if at least one of the sera tested between 1 year before and 1 year after the transplantation had a positive antibody response against these specific HPV genera or specific HPV types.

      Statistical analyses

      Baseline characteristics are given in percentages and simple mean calculation. Kaplan–Meier survival analyses were used in conjunction with multivariable Cox proportional hazard methods to calculate corresponding HRs with 95% CIs for development of KC in relation to HPV serology. The date of transplantation was used as opening date. As censoring dates we used the date of diagnosis of first KC, the date of the patient’s death, the date of last follow-up, the date lost to follow-up, or the date of end of study (29 October 2012). Factors considered as potential confounders were sex, type of transplantation, and age at transplantation. The Cox models were adjusted for these variables. As the variable of HPV seropositivity had a time scope of 1year before and after the transplantation (starting point), time-dependent Cox regression analyses were performed. The analyses of SCC were performed regardless of the presence of BCC among the cases and controls and the analyses of BCC were performed regardless of the presence of SCC among the cases and controls.
      To investigate whether there were Beta-PV types that were the driving force causing the positive association between Beta-PV seropositivity and the development of KC, we formed three subgroups: (1) patients who were seronegative for any Beta-PV type, (2) patients who were seropositive for the specific Beta-PV type regardless of the other Beta-PV types, and (3) patients who were seronegative for the specific Beta-PV type regardless of the other Beta-PV types. The cumulative incidence of KC, SCC, and BCC was calculated with Kaplan–Meier survival analyses. In addition, we performed conditional step forward logistic regression to get a rough estimation of which individual Beta-PV may possibly contribute to the association between seropositivity and the development of KC, SCC, and/or BCC. Starting with no variables in the model, all Beta-PV types, age, sex, and type of organ were introduced into the model, starting with the most relevant one until the addition of new variables did not improve the model any longer.
      All analyses were performed with SPSS 20.0 software (IBM SPSS Statistics, IBM, Chicago, IL).

      SUPPLEMENTARY MATERIAL

      Supplementary material is linked to the online version of the paper at http://www.nature.com/jid

      REFERENCES

        • Aldabagh B.
        • Angeles J.G.
        • Cardones A.R.
        • et al.
        Cutaneous squamous cell carcinoma and human papillomavirus: is there an association?.
        Dermatol Surg. 2013; 39: 1-23
        • Alotaibi L.
        • Provost N.
        • Gagnon S.
        • et al.
        Diversity of cutaneous human papillomavirus types in individuals with and without skin lesion.
        J Clin Virol. 2006; 36: 133-140
        • Andersson K.
        • Michael K.M.
        • Luostarinen T.
        • et al.
        Prospective study of human papillomavirus seropositivity and risk of nonmelanoma skin cancer.
        Am J Epidemiol. 2012; 175: 685-695
        • Andersson K.
        • Waterboer T.
        • Kirnbauer R.
        • et al.
        Seroreactivity to cutaneous human papillomaviruses among patients with nonmelanoma skin cancer or benign skin lesions.
        Cancer Epidemiol Biomarkers Prev. 2008; 17: 189-195
        • Antonsson A.
        • Waterboer T.
        • Bouwes Bavinck J.N.
        • et al.
        Longitudinal study of seroprevalence and serostability of 34 human papillomavirus types in European organ transplant recipients.
        Virology. 2013; 436: 91-99
        • Bernard H.U.
        • Burk R.D.
        • Chen Z.
        • et al.
        Classification of papillomaviruses (PVs) based on 189 PV types and proposal of taxonomic amendments.
        Virology. 2010; 401: 70-79
        • Bouvard V.
        • Baan R.
        • Straif K.
        • et al.
        A review of human carcinogens—Part B: biological agents.
        Lancet Oncol. 2009; 10: 321-322
        • Bouwes Bavinck J.N.
        • Feltkamp M.C.
        Milk of human kindness?—HAMLET, human papillomavirus, and warts.
        N Engl J Med. 2004; 350: 2639-2642
        • Bouwes Bavinck J.N.
        • Gissmann L.
        • Claas F.H.
        • et al.
        Relation between skin cancer, humoral responses to human papillomaviruses, and HLA class II molecules in renal transplant recipients.
        J Immunol. 1993; 151: 1579-1586
        • Bouwes Bavinck J.N.
        • Neale R.E.
        • Abeni D.
        • et al.
        Multicenter study of the association between betapapillomavirus infection and cutaneous squamous cell carcinoma.
        Cancer Res. 2010; 70: 9777-9786
        • Bouwes Bavinck J.N.
        • Plasmeijer E.I.
        • Feltkamp M.C.
        Beta-papillomavirus infection and skin cancer.
        J Invest Dermatol. 2008; 128: 1355-1358
        • Bouwes Bavinck J.N.
        • Vermeer B.J.
        • van der Woude F.J.
        • et al.
        Relation between skin cancer and HLA antigens in renal-transplant recipients.
        N Engl J Med. 1991; 325: 843-848
        • Bzhalava D.
        • Guan P.
        • Franceschi S.
        • et al.
        A systematic review of the prevalence of mucosal and cutaneous human papillomavirus types.
        Virology. 2013; 445: 224-231
        • Casabonne D.
        • Michael K.M.
        • Waterboer T.
        • et al.
        A prospective pilot study of antibodies against human papillomaviruses and cutaneous squamous cell carcinoma nested in the Oxford component of the European Prospective Investigation into Cancer and Nutrition.
        Int J Cancer. 2007; 121: 1862-1868
        • Chen A.C.
        • McMillan N.A.
        • Antonsson A.
        Human papillomavirus type spectrum in normal skin of individuals with or without a history of frequent sun exposure.
        J Gen Virol. 2008; 89: 2891-2897
        • de Koning M.N.
        • Weissenborn S.J.
        • Abeni D.
        • et al.
        Prevalence and associated factors of betapapillomavirus infections in individuals without cutaneous squamous cell carcinoma.
        J Gen Virol. 2009; 90: 1611-1621
        • de Villiers E.M.
        • Fauquet C.
        • Broker T.R.
        • et al.
        Classification of papillomaviruses.
        Virology. 2004; 324: 17-27
        • Euvrard S.
        • Kanitakis J.
        • Claudy A.
        Skin cancers after organ transplantation.
        N Engl J Med. 2003; 348: 1681-1691
        • Feltkamp M.C.
        • Broer R.
        • di Summa F.M.
        • et al.
        Seroreactivity to epidermodysplasia verruciformis-related human papillomavirus types is associated with nonmelanoma skin cancer.
        Cancer Res. 2003; 63: 2695-2700
        • Hartevelt M.M.
        • Bouwes Bavinck J.N.
        • Kootte A.M.
        • et al.
        Incidence of skin cancer after renal transplantation in The Netherlands.
        Transplantation. 1990; 49: 506-509
        • Iannacone M.R.
        • Gheit T.
        • Pfister H.
        • et al.
        Case-control study of genus-beta human papillomaviruses in plucked eyebrow hairs and cutaneous squamous cell carcinoma.
        Int J Cancer. 2014; 134: 2231-2244
        • Iannacone M.R.
        • Gheit T.
        • Waterboer T.
        • et al.
        Case-control study of cutaneous human papillomaviruses in squamous cell carcinoma of the skin.
        Cancer Epidemiol Biomarkers Prev. 2012; 21: 1303-1313
        • Iannacone M.R.
        • Gheit T.
        • Waterboer T.
        • et al.
        Case-control study of cutaneous human papillomavirus infection in basal cell carcinoma of the skin.
        J Invest Dermatol. 2013; 133: 1512-1520
        • IARC
        Human papillomaviruses.
        IARC Monogr Eval Carcinog Risks Hum. 2007; 90: 1-636
        • Jablonska S.
        • Dabrowski J.
        • Jakubowicz K.
        Epidermodysplasia verruciformis as a model in studies on the role of papovaviruses in oncogenesis.
        Cancer Res. 1972; 32: 583-589
        • Jackson S.
        • Harwood C.
        • Thomas M.
        • et al.
        Role of Bak in UV-induced apoptosis in skin cancer and abrogation by HPV E6 proteins.
        Genes Dev. 2000; 14: 3065-3073
        • Jackson S.
        • Storey A.
        E6 proteins from diverse cutaneous HPV types inhibit apoptosis in response to UV damage.
        Oncogene. 2000; 19: 592-598
        • Karagas M.R.
        • Nelson H.H.
        • Sehr P.
        • et al.
        Human papillomavirus infection and incidence of squamous cell and basal cell carcinomas of the skin.
        J Natl Cancer Inst. 2006; 98: 389-395
        • Karagas M.R.
        • Waterboer T.
        • Li Z.
        • et al.
        Genus beta human papillomaviruses and incidence of basal cell and squamous cell carcinomas of skin: population based case-control study.
        BMJ. 2010; 341: c2986
        • Lutzner M.A.
        Epidermodysplasia verruciformis. An autosomal recessive disease characterized by viral warts and skin cancer. A model for viral oncogenesis.
        Bull Cancer. 1978; 65: 169-182
        • Michael K.M.
        • Waterboer T.
        • Sehr P.
        • et al.
        Seroprevalence of 34 human papillomavirus types in the German general population.
        PLoS Pathog. 2008; 4: e1000091
        • Moloney F.J.
        • Comber H.
        • O’Lorcain P.
        • et al.
        A population-based study of skin cancer incidence and prevalence in renal transplant recipients.
        Br J Dermatol. 2006; 154: 498-504
        • Neale R.E.
        • Weissenborn S.
        • Abeni D.
        • et al.
        Human papillomavirus load in eyebrow hair follicles and risk of cutaneous squamous cell carcinoma.
        Cancer Epidemiol Biomarkers Prev. 2013; 22: 719-727
        • Paolini F.
        • Carbone A.
        • Benevolo M.
        • et al.
        Human papillomaviruses, p16INK4a and Akt expression in basal cell carcinoma.
        J Exp Clin Cancer Res. 2011; 30: 108
        • Paradisi A.
        • Waterboer T.
        • Sampogna F.
        • et al.
        Seropositivity for human papillomavirus and incidence of subsequent squamous cell and basal cell carcinomas of the skin in patients with a previous nonmelanoma skin cancer.
        Br J Dermatol. 2011; 165: 782-791
        • Parkin D.M.
        The global health burden of infection-associated cancers in the year 2002.
        Int J Cancer. 2006; 118: 3030-3044
        • Plasmeijer E.I.
        • Pandeya N.
        • O’Rourke P.
        • et al.
        The association between cutaneous squamous cell carcinoma and betapapillomavirus seropositivity: a cohort study.
        Cancer Epidemiol Biomarkers Prev. 2011; 20: 1171-1177
        • Proby C.M.
        • Harwood C.A.
        • Neale R.E.
        • et al.
        A case-control study of betapapillomavirus infection and cutaneous squamous cell carcinoma in organ transplant recipients.
        Am J Transplant. 2011; 11: 1498-1508
        • Stark S.
        • Petridis A.K.
        • Ghim S.J.
        • et al.
        Prevalence of antibodies against virus-like particles of Epidermodysplasia verruciformis-associated HPV8 in patients at risk of skin cancer.
        J Invest Dermatol. 1998; 111: 696-701
        • Steger G.
        • Olszewsky M.
        • Stockfleth E.
        • et al.
        Prevalence of antibodies to human papillomavirus type 8 in human sera.
        J Virol. 1990; 64: 4399-4406
        • Struijk L.
        • van der Meijden E.
        • Kazem S.
        • et al.
        Specific betapapillomaviruses associated with squamous cell carcinoma of the skin inhibit UVB-induced apoptosis of primary human keratinocytes.
        J Gen Virol. 2008; 89: 2303-2314
        • Terhorst D.
        • Drecoll U.
        • Stockfleth E.
        • et al.
        Organ transplant recipients and skin cancer: assessment of risk factors with focus on sun exposure.
        Br J Dermatol. 2009; 161: 85-89
        • Tessari G.
        • Naldi L.
        • Boschiero L.
        • et al.
        Incidence and clinical predictors of a subsequent nonmelanoma skin cancer in solid organ transplant recipients with a first nonmelanoma skin cancer: a multicenter cohort study.
        Arch Dermatol. 2010; 146: 294-299
        • Underbrink M.P.
        • Howie H.L.
        • Bedard K.M.
        • et al.
        E6 proteins from multiple human betapapillomavirus types degrade Bak and protect keratinocytes from apoptosis after UVB irradiation.
        J Virol. 2008; 82: 10408-10417
        • Waterboer T.
        • Abeni D.
        • Sampogna F.
        • et al.
        Serological association of beta and gamma human papillomaviruses with squamous cell carcinoma of the skin.
        Br J Dermatol. 2008; 159: 457-459
        • Waterboer T.
        • Sehr P.
        • Michael K.M.
        • et al.
        Multiplex human papillomavirus serology based on in situ-purified glutathione s-transferase fusion proteins.
        Clin Chem. 2005; 51: 1845-1853
        • Wisgerhof H.C.
        • van der Boog P.J.
        • de Fijter J.W.
        • et al.
        Increased risk of squamous-cell carcinoma in simultaneous pancreas kidney transplant recipients compared with kidney transplant recipients.
        J Invest Dermatol. 2009; 129: 2886-2894
        • zur Hausen H.
        Immortalization of human cells and their malignant conversion by high risk human papillomavirus genotypes.
        Semin Cancer Biol. 1999; 9: 405-411