Abstract
Background: We previously reported a significant association between higher UV radiation exposure before diagnosis and greater survival with melanoma in a population-based study in Connecticut. We sought to evaluate the hypothesis that sun exposure before diagnosis was associated with greater survival in a larger, international population-based study with more detailed exposure information.
Methods: We conducted a multicenter, international population-based study in four countries—Australia, Italy, Canada, and the United States—with 3,578 cases of melanoma with an average of 7.4 years of follow-up. Measures of sun exposure included sunburn, intermittent exposure, hours of holiday sun exposure, hours of water-related outdoor activities, ambient ultraviolet B (280–320 nm) dose, histologic solar elastosis, and season of diagnosis.
Results: Results were not strongly supportive of the earlier hypothesis. Having had any sunburn in 1 year within 10 years of diagnosis was inversely associated with survival; solar elastosis—a measure of lifetime cumulative exposure—was not. In addition, none of the intermittent exposure measures—water-related activities and sunny holidays—were associated with melanoma-specific survival. Estimated ambient UVB dose was not associated with survival.
Conclusion: Although there was an apparent protective effect of sunburns within 10 years of diagnosis, there was only weak evidence in this large, international, population-based study of melanoma that sun exposure before diagnosis is associated with greater melanoma-specific survival.
Impact: This study adds to the evidence that sun exposure before melanoma diagnosis has little effect on survival with melanoma. Cancer Epidemiol Biomarkers Prev; 23(10); 2145–52. ©2014 AACR.
Introduction
UV radiation (UVR) exposure is the major environmental risk factor for the development of melanoma (1) with intermittent UVR exposure, including sunburn, generally the measure of sun exposure most strongly associated with the development of melanoma (2, 3). In a Connecticut population-based study of 650 melanoma cases followed for an average of 5 years, Berwick and colleagues (4) reported that several measures of UVR before the diagnosis of melanoma were inversely associated with mortality from melanoma, suggesting that something about sun exposure, possibly its role in vitamin D production, was limiting cancer progression. Subsequently, Newton-Bishop and colleagues in a UK study of 872 patients with melanoma (5) reported that serum vitamin D levels were higher among those with better overall survival, and Rosso and colleagues in a European study of 260 patients with melanoma (6) found that patients with melanoma with more sunny vacations before diagnosis had better melanoma-specific survival. Laboratory studies have shown that vitamin D suppresses tumor proliferation (7) and suggest that increased vitamin D levels might keep a melanoma “in check.” To test the hypothesis that increased sun exposure before diagnosis is associated with improved survival from melanoma, we evaluated measures of solar UVR exposure before diagnosis in 3,578 incident patients with melanoma in the Genes, Environment and Melanoma study (GEM), an international, population-based study (8).
Materials and Methods
Subjects
A detailed description of the methods used in this study is available elsewhere (9). Briefly, this multicenter, international population-based study was conducted in four countries through eight population-based tumor registries—in Australia in the states of New South Wales and Tasmania, in Italy in the province of Piedmont, in Canada in the provinces of British Columbia and Ontario, and in the United States in the state of New Jersey, a 39-county area of North Carolina, two Southern California cancer registry populations (the Orange County Registry and the San Diego/Imperial Organization for Cancer Control), and through a hospital-based registry in the state of Michigan.
Institutional review board approval was obtained from all centers and written informed consent was obtained before interview. We interviewed 2,372 patients with incident first primary melanoma cases and 1,206 with incident multiple primary cases. Of the 1,206 with multiple primary cases, 96 had been first ascertained with single primaries. Single primary melanoma cases were diagnosed in 2000 and multiple primary cases from 1998 (British Columbia, California, New Jersey, and Tasmania) or 2000 (New South Wales, North Carolina, and Ontario) to 2003.
The overall participation rate was 54% for individuals completing all aspects of the study and submitting a DNA sample.
Data collection
A structured questionnaire administered by telephone assessed basic demographics, phenotypic characteristics, family history of cancer, recreational and occupational sun exposure at each decade of life, sunbed use, changes in sun-related behavior after a melanoma diagnosis, and a lifetime residential history. Nevi on the back were self-assessed using a set of photos and by reference to charts showing different patterns of nevi and freckles as previously described (2, 9).
UVR exposure measures
We evaluated effects on survival of measures of UVR exposure in various periods before diagnosis.
Sunburns.
Individuals reported whether they had been burned severely enough to have pain or blisters for two or more days in a specified year in the 10 years before diagnosis. This was coded as “once or more” or “never.”
Solar elastosis.
Solar elastosis, an indicator of sun exposure accumulated over a lifetime (10), was evaluated on histopathologic slide review as absent or present. Slides from 2,781 (78%) subjects were reviewed by expert dermatopathologists (L. From, K. Busam, and P.A. Groben) to standardize pathologic criteria and add parameters that community pathology laboratories often do not report, such as solar elastosis. Inter-reviewer reliability for solar elastosis was assessed as very good (κ = 0.65).
Intermittent sun exposure.
In a previous GEM analysis, two variables were considered to represent intermittent sun exposure—hours of holiday sun exposure in a place sunnier than usual residence and hours of water-related outdoor activities (2). These measures for 1 year in the most recent decade were categorized into quartiles based on the distribution among the entire population and ranked from low (quartile one) to high (quartile four).
UVB radiation dose.
Individual residential histories were coded for latitude, longitude, and altitude from birth to age at diagnosis, and then ambient UVB irradiances were calculated for each decade of age from records of satellite measurements of irradiance at the earth's surface as unweighted wavelength integrated spectral irradiance between 280 and 320 nm. UVB was used in analyses as this wavelength is thought to be the most effective in inducing serum vitamin D levels. Details of the calculations are available in Thomas and colleagues (10). Ambient UVB levels in the decade of life that included the melanoma diagnosis, at age 10 and over the lifetime (at each decade), were multiplied by the reported time spent outdoors on weekends and weekdays in the same period and categorized into quartiles based on the distribution among the entire population.
Season of diagnosis.
Diagnoses were classified by season, with data pooled for summer (December to February in the Southern hemisphere and June to August in the Northern), autumn (March to May in the Southern hemisphere and September to November in the Northern), winter (June to August in the Southern hemisphere and December to February in the Northern), and spring (September to November in the Southern hemisphere and March to May in the Northern).
Follow-up for survival
Patient follow-up for vital status was complete through 2007 except in British Columbia and Turin, where vital status was complete through 2008. Date and cause of death 7 years after diagnosis were obtained from National Death Indexes, cancer registries, and municipal records. We analyzed an average of 7.4 years of melanoma-specific survival. Individuals were classified as “died of melanoma,” “died of other cause,” and “alive at the end of follow-up.” An event was considered death due to melanoma. Among patients with multiple primaries, Breslow thickness (see Supplementary Tables S1–S3) and anatomic site for the thickest of their lesions were used in statistical models.
Data analysis
Cox proportional hazards models were used to calculate HRs and 95% confidence intervals (CI) for associations of categories of each exposure variable with melanoma outcome. Time to death from melanoma from diagnosis for those with single primaries or the most recent melanoma for those with multiple primaries was the outcome. Those who died of other causes or who were still alive at follow-up were censored in this analysis.
Age at diagnosis, sex, recruitment center, education level, and anatomic site were potential confounders of the association of sun exposure measures and melanoma survival. We found that there was no difference in effects of sun exposure measures and survival by primary status and therefore included both single and multiple primary melanomas in analyses to improve precision and included an indicator variable for primary status in all models. Kricker and colleagues (11) previously reported that there was no survival difference between multiple and single primaries in GEM. A time-dependent covariate was used for the 96 patients who developed a second primary during the study follow-up period. Pigmentary characteristics, prior history of nonmelanoma skin cancer, and family history of melanoma were assessed but found not to be potential confounders of sun exposure measures in relation to survival. Stratified analyses were conducted to determine if any effect of sun exposure measures on risk of death from melanoma was modified by MC1R status (with or without “red hair color” variants D84E, R151C, R160W, and D294H), ability to tan (good and poor), and propensity to sunburn (high and low). Likelihood ratio tests for heterogeneity were used to evaluate significance of any apparent effect modification. Tests for linear trend were performed for ordered categorical variables. All tests were two-sided, and P < 0.05 was considered statistically significant. All data were analyzed using SAS 9.3 (Cary, NC).
Results
Of the 3,578 eligible individuals diagnosed with melanoma in this study (2,007 males and 1,571 females), 563 died by the end of follow-up (15.7%): 255 (7.1%) from melanoma and 308 (8.6%) from other causes.
Survival analyses are presented as baseline models, with HRs adjusted for center, age, sex, primary status, and the time-dependent covariate, and as fully adjusted models, which included the above variables as well as others significantly associated with survival, educational level, and anatomic site.
Clinical and host characteristics and melanoma-specific survival
Anticipated associations for host and clinical characteristics were seen (Table 1). Primary status was not associated with hazard of death from melanoma in the fully adjusted model. Women had a lower risk of dying from melanoma in both the baseline model (P < 0.001) and the fully adjusted model (P = 0.0002). The hazard of death increased with increasing age (fully adjusted HR 1.02 for each year of age; 95% CI, 1.01–1.03, P < 0.0001). Melanomas on the arms were at lowest risk for poor survival relative to melanoma of the head and neck (fully adjusted HR, 0.47; 95% CI, 0.31–0.71, P = 0.003). Relative to superficial spreading melanoma, the fully adjusted HR for lentigo maligna melanoma was decreased (HR, 0.57; 95% CI, 0.33–0.98, P = 0.04). Breslow thickness (fully adjusted HR, 13.79; 95% CI, 9.12–20.84, for thickness of 4.00 mm or higher relative to thickness of less than 1.00 mm) was strongly and significantly associated with poor prognosis (P < 0.001). Similar to most other studies, those with more education had a significantly reduced hazard of dying from melanoma (fully adjusted HR, 0.56; 95% CI, 0.40–0.78, P = 0.0005). Having a family history of melanoma (fully adjusted HR, 0.85; 95% CI, 0.58–1.24, P = 0.39) or a prior history of nonmelanoma skin cancer (fully adjusted HR, 0.93; 95% CI, 0.71–1.23, P = 0.63) did not affect the hazard of dying from melanoma.
Variable . | Level . | Subjects in study (N = 3,578) . | Subjects died from melanoma (N = 255) . | Baseline modela HR (95% CI) . | Fully adjusted modelb HR (95% CI) . |
---|---|---|---|---|---|
Primary status | Single | 2,372 | 152 | 1.00 | 1.00 |
Multiple | 1,206 | 103 | 0.99 (0.75–1.32) | 1.03 (0.78–1.37) | |
P value | 0.98 | 0.83 | |||
Sex | Male | 2,007 | 184 | 1.00 | 1.00 |
Female | 1,571 | 71 | 0.56 (0.43–0.75) | 0.56 (0.42–0.76) | |
P value | <0.001 | 0.002 | |||
Age at diagnosis | Per year | 1.03 (1.02–1.04) | 1.02 (1.01–1.03) | ||
P value | <0.001 | < 0.001 | |||
Anatomic site | Head and neck | 578 | 77 | 1.00 | 1.00 |
Trunk | 1,585 | 107 | 0.54 (0.40–0.73) | 0.53 (0.39–0.73) | |
Arms | 666 | 34 | 0.47 (0.31–0.71) | 0.47 (0.31–0.71) | |
Legs | 749 | 37 | 0.51 (0.33–0.77) | 0.51 (0.34–0.78) | |
Histology | SSM | 2,302 | 106 | 1.00 | 1.00 |
NM | 333 | 70 | 4.27 (3.13–5.81) | 3.74 (2.72–5.14) | |
LMM | 366 | 18 | 0.85 (0.51–1.41) | 0.57 (0.33–0.98) | |
ALM | 16 | 3 | 8.99 (3.62–22.36) | 9.90 (3.87–25.38) | |
NOS | 496 | 40 | 1.95 (1.33–2.86) | 1.85 (1.26–2.73) | |
Other | 65 | 18 | 4.51 (2.59–8.15) | 3.04 (1.65–5.61) | |
Breslow thickness | 0.01–1.00 | 2,228 | 45 | 1.00 | 1.00 |
1.01–2.00 | 727 | 79 | 5.33 (3.69–7.70) | 5.13 (3.53–7.41) | |
2.01–4.00 | 361 | 75 | 10.06 (6.92–14.60) | 9.65 (0.62–14.07) | |
>4.00 | 175 | 52 | 15.03 (10.02–22.53) | 13.81 (9.13–20.88) | |
Missing | 87 | 4 | |||
P value for trend | <0.001 | <0.001 | |||
Education | <College | 2,415 | 203 | 1.00 | 1.00 |
College+ | 1,133 | 47 | 0.56 (0.40–0.78) | 0.69 (0.49–0.97) | |
P value | 0.0006 | 0.03 | |||
Family history of melanoma | |||||
None | 2,953 | 212 | 1.00 | 1.00 | |
Present | 551 | 31 | 0.82 (0.56–1.21) | 0.85 (0.58–1.24) | |
Do not know | 74 | 12 | |||
P value | |||||
History of NMSC | None | 2,449 | 187 | 1.00 | 1.00 |
Yes | 1,081 | 86 | 0.91 (0.69–1.20) | 0.93 (0.70–1.23) | |
Do not know | 48 | 2 | |||
P value | 0.51 | 0.59 |
Variable . | Level . | Subjects in study (N = 3,578) . | Subjects died from melanoma (N = 255) . | Baseline modela HR (95% CI) . | Fully adjusted modelb HR (95% CI) . |
---|---|---|---|---|---|
Primary status | Single | 2,372 | 152 | 1.00 | 1.00 |
Multiple | 1,206 | 103 | 0.99 (0.75–1.32) | 1.03 (0.78–1.37) | |
P value | 0.98 | 0.83 | |||
Sex | Male | 2,007 | 184 | 1.00 | 1.00 |
Female | 1,571 | 71 | 0.56 (0.43–0.75) | 0.56 (0.42–0.76) | |
P value | <0.001 | 0.002 | |||
Age at diagnosis | Per year | 1.03 (1.02–1.04) | 1.02 (1.01–1.03) | ||
P value | <0.001 | < 0.001 | |||
Anatomic site | Head and neck | 578 | 77 | 1.00 | 1.00 |
Trunk | 1,585 | 107 | 0.54 (0.40–0.73) | 0.53 (0.39–0.73) | |
Arms | 666 | 34 | 0.47 (0.31–0.71) | 0.47 (0.31–0.71) | |
Legs | 749 | 37 | 0.51 (0.33–0.77) | 0.51 (0.34–0.78) | |
Histology | SSM | 2,302 | 106 | 1.00 | 1.00 |
NM | 333 | 70 | 4.27 (3.13–5.81) | 3.74 (2.72–5.14) | |
LMM | 366 | 18 | 0.85 (0.51–1.41) | 0.57 (0.33–0.98) | |
ALM | 16 | 3 | 8.99 (3.62–22.36) | 9.90 (3.87–25.38) | |
NOS | 496 | 40 | 1.95 (1.33–2.86) | 1.85 (1.26–2.73) | |
Other | 65 | 18 | 4.51 (2.59–8.15) | 3.04 (1.65–5.61) | |
Breslow thickness | 0.01–1.00 | 2,228 | 45 | 1.00 | 1.00 |
1.01–2.00 | 727 | 79 | 5.33 (3.69–7.70) | 5.13 (3.53–7.41) | |
2.01–4.00 | 361 | 75 | 10.06 (6.92–14.60) | 9.65 (0.62–14.07) | |
>4.00 | 175 | 52 | 15.03 (10.02–22.53) | 13.81 (9.13–20.88) | |
Missing | 87 | 4 | |||
P value for trend | <0.001 | <0.001 | |||
Education | <College | 2,415 | 203 | 1.00 | 1.00 |
College+ | 1,133 | 47 | 0.56 (0.40–0.78) | 0.69 (0.49–0.97) | |
P value | 0.0006 | 0.03 | |||
Family history of melanoma | |||||
None | 2,953 | 212 | 1.00 | 1.00 | |
Present | 551 | 31 | 0.82 (0.56–1.21) | 0.85 (0.58–1.24) | |
Do not know | 74 | 12 | |||
P value | |||||
History of NMSC | None | 2,449 | 187 | 1.00 | 1.00 |
Yes | 1,081 | 86 | 0.91 (0.69–1.20) | 0.93 (0.70–1.23) | |
Do not know | 48 | 2 | |||
P value | 0.51 | 0.59 |
Abbreviations: ALM, acral lentiginous melanoma; LMM, lentigo maligna melanoma; NMSC, nonmelanoma skin cancer; NM, nodular melanoma; NOS, not otherwise specified; SSM, superficial spreading melanoma.
aAdjusted for center, primary status, crossover time-dependent status, age at diagnosis, and sex.
bAdjusted for center, primary status, crossover time-dependent status, age at diagnosis, sex, anatomic site, and education.
Recent sun exposure
We found a reduced HR of melanoma death with one or more sunburns in a year in the decade before diagnosis (fully adjusted HR, 0.27; 95% CI, 0.09–0.85, P = 0.03; Table 2). Other sun exposure variables in the decade before diagnosis, including holiday sun hours in a place sunnier than usual residence and hours of water-related activities and estimated UVB dose, and season of diagnosis, were not significantly associated with survival from melanoma in either the baseline or the fully adjusted models.
Variable . | Level . | Number of subjects in study . | Number of subjects died from melanoma . | Baseline modela HR (95% CI) . | Fully adjusted modelb HR (95% CI) . |
---|---|---|---|---|---|
Sunburns within 10 years of diagnosis | 0 | 3,246 | 240 | 1.00 | 1.00 |
1+ | 252 | 4 | 0.36 (0.13–0.98) | 0.27 (0.09–0.85) | |
Missing | 80 | 11 | |||
P value | 0.05 | 0.03 | |||
Holiday sun hours within 10 years of diagnosis | 0 | 1,852 | 145 | 1.00 | 1.00 |
>0–<56.5 | 740 | 37 | 0.75 (0.52–1.08) | 0.77 (0.53–1.11) | |
56.5+ | 739 | 55 | 0.87 (0.63–1.21) | 0.90 (0.65–1.25) | |
Missing | 247 | 18 | |||
P for trend | 0.28 | 0.38 | |||
Water-related activities within 10 years of diagnosis | 0–<1,314 | 848 | 60 | 1.00 | 1.00 |
1,314–<3,120 | 830 | 55 | 0.89 (0.61–1.28) | 0.86 (0.59–1.26) | |
3,120–<6,140 | 868 | 63 | 0.87 (0.60–1.25) | 0.87 (0.60–1.26) | |
6,140+ | 848 | 64 | 0.88 (0.61–1.27) | 0.84 (0.48–1.22) | |
Missing | 183 | 13 | |||
P for trend | 0.51 | 0.41 | |||
UVB dose within 10 years of diagnosis (kJ/m2) | 0–<2,134 | 836 | 58 | 1.00 | 1.00 |
2,134–<3,757 | 838 | 45 | 0.70 (0.47–1.05) | 0.68 (0.46–1.01) | |
3,757–<6,413 | 837 | 67 | 0.95 (0.66–1.37) | 0.86 (0.59–1.24) | |
6,413+ | 837 | 69 | 0.78 (0.53–1.15) | 0.69 (0.47–1.02) | |
Missing | |||||
P for trend | 0.51 | 0.18 | |||
Season of diagnosis | Winter | 741 | 52 | 1.00 | 1.00 |
Fall | 962 | 71 | 0.82 (0.64–1.04) | 0.90 (0.62–1.30) | |
Spring | 803 | 50 | 0.87 (0.69–1.11) | 0.97 (0.65–1.45) | |
Summer | 1,060 | 81 | 0.89 (0.79–1.16) | 1.09 (0.77–1.55) | |
P for trend | 0.49 | 0.46 |
Variable . | Level . | Number of subjects in study . | Number of subjects died from melanoma . | Baseline modela HR (95% CI) . | Fully adjusted modelb HR (95% CI) . |
---|---|---|---|---|---|
Sunburns within 10 years of diagnosis | 0 | 3,246 | 240 | 1.00 | 1.00 |
1+ | 252 | 4 | 0.36 (0.13–0.98) | 0.27 (0.09–0.85) | |
Missing | 80 | 11 | |||
P value | 0.05 | 0.03 | |||
Holiday sun hours within 10 years of diagnosis | 0 | 1,852 | 145 | 1.00 | 1.00 |
>0–<56.5 | 740 | 37 | 0.75 (0.52–1.08) | 0.77 (0.53–1.11) | |
56.5+ | 739 | 55 | 0.87 (0.63–1.21) | 0.90 (0.65–1.25) | |
Missing | 247 | 18 | |||
P for trend | 0.28 | 0.38 | |||
Water-related activities within 10 years of diagnosis | 0–<1,314 | 848 | 60 | 1.00 | 1.00 |
1,314–<3,120 | 830 | 55 | 0.89 (0.61–1.28) | 0.86 (0.59–1.26) | |
3,120–<6,140 | 868 | 63 | 0.87 (0.60–1.25) | 0.87 (0.60–1.26) | |
6,140+ | 848 | 64 | 0.88 (0.61–1.27) | 0.84 (0.48–1.22) | |
Missing | 183 | 13 | |||
P for trend | 0.51 | 0.41 | |||
UVB dose within 10 years of diagnosis (kJ/m2) | 0–<2,134 | 836 | 58 | 1.00 | 1.00 |
2,134–<3,757 | 838 | 45 | 0.70 (0.47–1.05) | 0.68 (0.46–1.01) | |
3,757–<6,413 | 837 | 67 | 0.95 (0.66–1.37) | 0.86 (0.59–1.24) | |
6,413+ | 837 | 69 | 0.78 (0.53–1.15) | 0.69 (0.47–1.02) | |
Missing | |||||
P for trend | 0.51 | 0.18 | |||
Season of diagnosis | Winter | 741 | 52 | 1.00 | 1.00 |
Fall | 962 | 71 | 0.82 (0.64–1.04) | 0.90 (0.62–1.30) | |
Spring | 803 | 50 | 0.87 (0.69–1.11) | 0.97 (0.65–1.45) | |
Summer | 1,060 | 81 | 0.89 (0.79–1.16) | 1.09 (0.77–1.55) | |
P for trend | 0.49 | 0.46 |
aAdjusted for center, primary status, crossover time-dependent status, age at diagnosis, and sex.
bAdjusted for center, primary status, crossover time-dependent status, age at diagnosis, sex, education, and anatomic site.
Early-life sun exposure
We found a significant trend for increasing melanoma mortality with increasing UVB dose at age 10, (fully adjusted HR, 1.49; 95% CI, 0.97–2.30; P = 0.03) for the highest quartile compared with the lowest. Other sun exposure variables in early life were not significantly associated with survival from melanoma (Table 3).
Variable . | Level . | Number of subjects in study . | Number of subjects died from melanoma . | Baseline modelHR (95% CI)a . | Fully adjusted modelHR (95% CI)b . |
---|---|---|---|---|---|
Sunburns—early life | 0 | 1,584 | 114 | 1.00 | 1.00 |
1+ | 1,496 | 104 | 1.03 (0.78–1.36) | 1.08 (0.81–1.42) | |
Missing | 498 | 37 | |||
P value | 0.82 | 0.61 | |||
Holiday sun hours—early life | 0 | 2,726 | 197 | 1.00 | 1.00 |
1+ | 769 | 52 | 1.10 (0.80–1.52) | 1.19 (0.86–1.67) | |
Missing | 83 | 6 | |||
P value | 0.56 | 0.29 | |||
Water-related activities—early life | 0–<386 | 849 | 57 | 1.00 | 1.00 |
386–<1,404 | 848 | 57 | 1.02 (0.70–1.49) | 1.02 (0.70–1.49) | |
1,404–<3,414 | 852 | 57 | 0.96 (0.66–1.41) | 0.91 (0.61–1.35) | |
3,414+ | 850 | 71 | 1.18 (0.82–1.70) | 1.17 (0.81–1.70) | |
Missing | 179 | 13 | |||
P for trend | 0.42 | 0.49 | |||
UVB dose—early life (kJ/m2) | 0–<3,333 | 839 | 47 | 1.00 | 1.00 |
3,333–<4,916.5 | 838 | 43 | 0.98 (0.64–1.50) | 0.93 (0.60–1.43) | |
4,916.5–<6,796 | 838 | 69 | 1.46 (0.97–2.20) | 1.35 (0.89–2.05) | |
6,796+ | 839 | 77 | 1.65 (1.07–2.52) | 1.49 (0.97–2.31) | |
Missing | 224 | 19 | |||
P for trend | 0.009 | 0.03 |
Variable . | Level . | Number of subjects in study . | Number of subjects died from melanoma . | Baseline modelHR (95% CI)a . | Fully adjusted modelHR (95% CI)b . |
---|---|---|---|---|---|
Sunburns—early life | 0 | 1,584 | 114 | 1.00 | 1.00 |
1+ | 1,496 | 104 | 1.03 (0.78–1.36) | 1.08 (0.81–1.42) | |
Missing | 498 | 37 | |||
P value | 0.82 | 0.61 | |||
Holiday sun hours—early life | 0 | 2,726 | 197 | 1.00 | 1.00 |
1+ | 769 | 52 | 1.10 (0.80–1.52) | 1.19 (0.86–1.67) | |
Missing | 83 | 6 | |||
P value | 0.56 | 0.29 | |||
Water-related activities—early life | 0–<386 | 849 | 57 | 1.00 | 1.00 |
386–<1,404 | 848 | 57 | 1.02 (0.70–1.49) | 1.02 (0.70–1.49) | |
1,404–<3,414 | 852 | 57 | 0.96 (0.66–1.41) | 0.91 (0.61–1.35) | |
3,414+ | 850 | 71 | 1.18 (0.82–1.70) | 1.17 (0.81–1.70) | |
Missing | 179 | 13 | |||
P for trend | 0.42 | 0.49 | |||
UVB dose—early life (kJ/m2) | 0–<3,333 | 839 | 47 | 1.00 | 1.00 |
3,333–<4,916.5 | 838 | 43 | 0.98 (0.64–1.50) | 0.93 (0.60–1.43) | |
4,916.5–<6,796 | 838 | 69 | 1.46 (0.97–2.20) | 1.35 (0.89–2.05) | |
6,796+ | 839 | 77 | 1.65 (1.07–2.52) | 1.49 (0.97–2.31) | |
Missing | 224 | 19 | |||
P for trend | 0.009 | 0.03 |
aAdjusted for center, primary status, crossover time-dependent status, age at diagnosis, and sex.
bAdjusted for center, primary status, crossover time-dependent status, age at diagnosis, sex, and anatomic site.
Lifetime average annual sun exposure
None of the lifetime cumulative or annual average sun exposure measures were associated either positively or negatively with melanoma-specific survival (Table 4). Solar elastosis was not associated with an increased risk of dying from melanoma in the baseline or the fully adjusted model (HR, 0.74; 95% CI, 0.52–1.07; P = 0.11). Lifetime annual average levels of holiday sun hours in a place sunnier than usual residence, water-related activities, and estimated solar UVB dose were also not significantly associated with melanoma-specific survival (Table 4).
Variable . | Level . | Number of subjects in study . | Number of melanoma deaths . | Baseline model HR (95% CI)a . | Fully adjusted model HR (95% CI)b . |
---|---|---|---|---|---|
Solar elastosis | |||||
Absent | 889 | 55 | 1.00 | 1.00 | |
Present | 1,892 | 141 | 0.88 (0.63–1.24) | 0.74 (0.51–1.06) | |
Missing | 797 | 59 | |||
P value | 0.47 | 0.10 | |||
Ever sunburned | No | 1,139 | 86 | 1.00 | 1.00 |
Yes | 2,174 | 167 | 1.05 (0.80–1.36) | 1.05 (0.80–1.37) | |
Missing | 12 | 2 | |||
P value | 0.75 | 0.75 | |||
Holiday sun hours (average annual) | |||||
0–<1.02 | 716 | 54 | 1.00 | 1.00 | |
1.02–<19.7 | 718 | 40 | 0.71 (0.47–1.07) | 0.76 (0.50–1.14) | |
19.7–<44.9 | 717 | 51 | 0.85 (0.58–1.25) | 0.88 (0.59–1.30) | |
44.9+ | 717 | 53 | 0.83 (0.56–1.23) | 0.88 (0.59–1.31) | |
Missing | 710 | 57 | |||
P for trend | 0.52 | 0.67 | |||
Water-related activities (average annual) | |||||
0–<0.8 | 886 | 58 | 1.00 | 1.00 | |
>0.8–<25.39 | 887 | 72 | 1.16 (0.81–1.65) | 1.23 (0.86–1.77) | |
25.39–<76.5 | 887 | 66 | 1.15 (0.80–1.65) | 1.25 (0.86–1.81) | |
76.5+ | 887 | 57 | 1.00 (0.69–1.45) | 1.08 (0.74–1.58) | |
Missing | 31 | 2 | |||
P for trend | 0.95 | 0.72 | |||
Average annual UVB dose (kJ/m2) | 0–<2,8572,857–<4,106.84,106.8–<5,8885,888+Missing | 822823823823287 | 5243578023 | 1.000.85 (0.55–1.31)1.08 (0.71–1.65)1.23 (0.80–1.89) | 1.000.79 (0.51–1.22)0.95 (0.62–1.45)1.10 (0.71–1.70) |
P for trend | 0.16 | 0.39 |
Variable . | Level . | Number of subjects in study . | Number of melanoma deaths . | Baseline model HR (95% CI)a . | Fully adjusted model HR (95% CI)b . |
---|---|---|---|---|---|
Solar elastosis | |||||
Absent | 889 | 55 | 1.00 | 1.00 | |
Present | 1,892 | 141 | 0.88 (0.63–1.24) | 0.74 (0.51–1.06) | |
Missing | 797 | 59 | |||
P value | 0.47 | 0.10 | |||
Ever sunburned | No | 1,139 | 86 | 1.00 | 1.00 |
Yes | 2,174 | 167 | 1.05 (0.80–1.36) | 1.05 (0.80–1.37) | |
Missing | 12 | 2 | |||
P value | 0.75 | 0.75 | |||
Holiday sun hours (average annual) | |||||
0–<1.02 | 716 | 54 | 1.00 | 1.00 | |
1.02–<19.7 | 718 | 40 | 0.71 (0.47–1.07) | 0.76 (0.50–1.14) | |
19.7–<44.9 | 717 | 51 | 0.85 (0.58–1.25) | 0.88 (0.59–1.30) | |
44.9+ | 717 | 53 | 0.83 (0.56–1.23) | 0.88 (0.59–1.31) | |
Missing | 710 | 57 | |||
P for trend | 0.52 | 0.67 | |||
Water-related activities (average annual) | |||||
0–<0.8 | 886 | 58 | 1.00 | 1.00 | |
>0.8–<25.39 | 887 | 72 | 1.16 (0.81–1.65) | 1.23 (0.86–1.77) | |
25.39–<76.5 | 887 | 66 | 1.15 (0.80–1.65) | 1.25 (0.86–1.81) | |
76.5+ | 887 | 57 | 1.00 (0.69–1.45) | 1.08 (0.74–1.58) | |
Missing | 31 | 2 | |||
P for trend | 0.95 | 0.72 | |||
Average annual UVB dose (kJ/m2) | 0–<2,8572,857–<4,106.84,106.8–<5,8885,888+Missing | 822823823823287 | 5243578023 | 1.000.85 (0.55–1.31)1.08 (0.71–1.65)1.23 (0.80–1.89) | 1.000.79 (0.51–1.22)0.95 (0.62–1.45)1.10 (0.71–1.70) |
P for trend | 0.16 | 0.39 |
aAdjusted for center, primary status, crossover time-dependent status, age at diagnosis, and sex.
bAdjusted for center, primary status, crossover time-dependent status, age at diagnosis, sex, anatomic site, and education.
Stratified analyses
There was little evidence that any association of sun exposure variables and hazard of death from melanoma varied among categories of MC1R status, ability to tan, and propensity to burn in relationship to melanoma survival (data not shown).
Discussion
This study of 3,578 highly annotated patients with melanoma shows the expected associations of host characteristics and clinical variables with survival, but provides only a little support for our previous study in Connecticut where sun exposure before diagnosis was inversely associated with melanoma survival, such that individuals with higher levels of intermittent sun exposure, presence of solar elastosis, and any sunburns before diagnosis had better survival. The present study found only an inverse association of sunburns within the 10 years before diagnosis with survival from melanoma. Lifetime sunburn history was not associated with survival with melanoma, which is opposite to the finding in the Connecticut study.
Analytic studies of sun exposure and melanoma survival are few. There are differences of study design and study population among the several studies that show an inverse association with either solar UVB or circulating serum vitamin D and survival compared with the present study. Lesions were generally somewhat deeper in the Connecticut study with a mean thickness of 1.81 mm (median, 0.81 mm) versus 1.30 mm (median, 0.78 mm) in this study. This difference is indicative of a general trend to diagnose thinner lesions over time (12). The inclusion of Breslow thickness in the fully adjusted model did not materially modify associations in models without its inclusion (Supplementary Tables S1–S3). It is important to note that because this study is population-based, it includes many individuals with very thin melanomas and hence high overall survival. Such population-based studies are critical for public health recommendations, but any particular effects of lifestyle on survival would be most relevant for the more selected group of people whose melanoma characteristics place them at a higher likelihood of mortality from melanoma.
In the Rosso and colleagues (6) study, the population from Turin, Italy, was quite small. The major variable associated with improved survival with melanoma was number of holidays to sunny places; it is possible that this variable is confounded with socioeconomic status, which has been found to be inversely associated with hazard of death from melanoma in three studies (13–15).
In the Newton-Bishop and colleagues (5) study, measures of circulating serum vitamin D were positively associated with relapse-free survival and lower Breslow thickness at diagnosis. This study did not look at melanoma-specific survival, but rather overall survival. In addition, only individuals with tumors greater than 0.75 mm were included. These results differ from our studies in Connecticut and the present GEM study that both focus on melanoma-specific survival and inclusion of all tumors unrestricted by Breslow thickness. We have evaluated overall survival, however, and found that several measures of intermittent sun exposure before diagnosis—UVB dose in quartiles (P for trend = 0.004), hours spent in water-related activities (P for trend = 0.01), and hours of holiday sun exposure (P for trend = 0.03)—are significantly and inversely associated with survival (Supplementary Table S4). Our data indicate a possible impact of sun exposure on overall survival; however, this study was not designed to evaluate deaths other than melanoma.
Several limitations deserve note, particularly the potential for misclassification in recalled sun exposure. Because the “dose” information relies on reported hours of sun exposure multiplied by the ambient exposure, there is the potential for misclassification that is likely nondirectional and would bias results to the null. In addition, although sunburn is likely subject to recall bias (16), the fact that sunburn represents overexposure to the sun, whereas exposure to high ambient levels of UV is modified by behaviors and phenotype, may make the single finding that sunburn before diagnosis is “protective” and more salient. Caution is necessary in interpreting that finding due to the very small number of deaths in the group experiencing sunburn (n = 4). Misclassification could also result from differences among centers in non–UV sun-related behaviors that might affect mortality in comparison with previous single-center studies where more uniform non-UV behaviors factors might be more uniform.
Another concern lies with the use of death certificates for verification of mortality as death certificates are sometimes misclassified (17). Each of the centers in this study had high-quality identification of deaths, using death certification, such as the National Death Index in the United States and Australia and the Provincial Cancer Registries in Canada. In Italy, deaths were verified by linking to the municipal rosters. If in fact a patient died from a metastasis from his melanoma but was classified as dying from another cancer, such as lung cancer, then our statistical power will have been reduced. Furthermore, it is noted that deaths from melanoma continue to occur over a relatively long period of time, and we have survival information for 7.4 years, so that a longer follow-up period may produce somewhat different results.
Many studies have demonstrated positive associations between solar UV exposure at season of diagnosis and survival from different cancers. Results are mixed although the majority of studies demonstrate that those cancers diagnosed in the fall, when circulating serum vitamin D levels are generally the highest, have better prognosis than those diagnosed in other seasons. For melanoma, one study found higher survival in patients diagnosed in summer or fall (18) and one did not (19); both were from Australia.
Our study's strengths include the large number of participants, the variety of latitudes, the relatively long follow-up, a reliable sun exposure questionnaire (20, 21), the ability to control for confounders, and the extensive pathologic review of cases.
In conclusion, this study provides only weak evidence that high levels of sun exposure before diagnosis have a benefit for melanoma survival.
Disclosure of Potential Conflicts of Interest
B.K. Armstrong is a Chair of Skin Cancer Prevention Advisory Committee at Cancer Institute NSW. No potential conflicts of interest were disclosed by the other authors.
Authors' Contributions
Conception and design: M. Berwick, A.S. Reiner, B.K. Armstrong, A. Kricker, K. Busam, I. Orlow, L.D. Marrett, R.P. Gallagher, H. Anton-Culver, S. Rosso, A. Venn, C.B. Begg, S.B. Gruber
Development of methodology: B.K. Armstrong, A. Kricker, N.E. Thomas, K. Busam, R.P. Gallagher, S. Rosso, T. Dwyer, C.B. Begg, S.B. Gruber, M. Berwick
Acquisition of data (provided animals, acquired and managed patients, provided facilities, etc.): M. Berwick, S. Paine, B.K. Armstrong, A. Kricker, N.E. Thomas, P.A. Groben, L. From, K. Busam, L.D. Marrett, R.P. Gallagher, H. Anton-Culver, S. Rosso, R. Zanetti, T. Dwyer, A. Venn, S.B. Gruber
Analysis and interpretation of data (e.g., statistical analysis, biostatistics, computational analysis): M. Berwick, A.S. Reiner, S. Paine, B.K. Armstrong, C. Goumas, A.E. Cust, N.E. Thomas, P.A. Groben, K. Busam, R.P. Gallagher, S. Rosso, J. Lee-Taylor, C.B. Begg, S.B. Gruber
Writing, review, and/or revision of the manuscript: M. Berwick, A.S. Reiner, B.K. Armstrong, A. Kricker, C. Goumas, A.E. Cust, N.E. Thomas, K. Busam, I. Orlow, L.D. Marrett, R.P. Gallagher, H. Anton-Culver, S. Rosso, R. Zanetti, P.A. Kanetsky, T. Dwyer, A. Venn, C.B. Begg, S.B. Gruber
Administrative, technical, or material support (i.e., reporting or organizing data, constructing databases): M. Berwick, A. Kricker, C. Goumas, N.E. Thomas, T. Dwyer, S.B. Gruber
Study supervision: M. Berwick, B.K. Armstrong, N.E. Thomas, L.D. Marrett, R.P. Gallagher, H. Anton-Culver, T. Dwyer, S.B. Gruber
Acknowledgments
The study was conducted by the GEM Study Group:
Coordinating Center, Memorial Sloan-Kettering Cancer Center, New York, NY, USA: Marianne Berwick [Principal Investigator (PI), currently at the University of New Mexico], Colin B. Begg (Co-PI), Irene Orlow (co-investigator), Klaus Busam (Dermatopathologist), Anne S. Reiner (Biostatistician), Pampa Roy (Laboratory Technician), Ajay Sharma (Laboratory Technician), Jaipreet Rayar (Laboratory Technician). The University of New Mexico, Albuquerque: Marianne Berwick (PI), Li Luo (Biostatistician), Kirsten White (Laboratory Manager), Susan Paine (Data Manager), Harold Nelson (Data Manager). Study centers included the following: The University of Sydney and The Cancer Council New South Wales, Sydney, Australia: Bruce K. Armstrong (PI), Anne Kricker (co-PI), Melisa Litchfield (Study Coordinator); Menzies Research Institute Tasmania, University of Tasmania, Hobart, Australia: Alison Venn (current PI), Terry Dwyer (PI, currently at Oxford University), Paul Tucker (Dermatopathologist); British Columbia Cancer Agency, Vancouver, Canada: Richard P. Gallagher (PI), Donna Kan (coordinator); Cancer Care Ontario, Toronto, Canada: Loraine D. Marrett (PI), Elizabeth Theis (co-investigator), Lynn From (Dermatopathologist); Center for Cancer Prevention, Torino, Italy: Roberto Zanetti (PI), Stefano Rosso (co-PI); University of California, Irvine, CA: Hoda Anton-Culver (PI), Argyrios Ziogas (Statistician); University of Michigan, Ann Arbor, MI: Stephen B. Gruber (PI, currently at the University of Southern California, Los Angeles, CA), Timothy Johnson (Director of Melanoma Program), Shu-Chen Huang (co-investigator, joint at USC-University of Michigan); New Jersey Department of Health and Senior Services, Trenton, NJ: Judith Klotz (PI, currently retired), Homer Wilcox (Co-PI, currently retired); University of North Carolina, Chapel Hill, NC: Nancy E. Thomas (PI), Robert C. Millikan (previous PI, deceased), David Ollila (co-investigator), Kathleen Conway (co-investigator), Pamela A. Groben (Dermatopathologist), Sharon N. Edmiston (Research Analyst), Honglin Hao (Laboratory Specialist), Elois Parrish (Laboratory Specialist); University of Pennsylvania, Philadelphia, PA: Timothy Rebbeck (PI), Peter A. Kanetsky (co-investigator). UV data consultants: Julia Lee-Taylor and Sasha Madronich, National Centre for Atmospheric Research, Boulder, CO.
Grant Support
This work was supported by the NIH (U01 CA 83101, R01 CA112524, R01 CA112524-05S2, and K05 CA13165, to M. Berwick; R01 CA112243, R01 CA112243-05S1, and P30 ES010126, to N.E. Thomas) and a Michael Smith Foundation for Health Research Infrastructure Award (to R.P. Gallagher). The National Center for Atmospheric Research is sponsored by the National Science Foundation
The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.