Matrix metalloproteinases (MMP) degrade various components of the extracellular matrix, and their overexpression has been implicated in tumor progression. Nonsynonymous single nucleotide polymorphisms (SNPs) lead to amino acid substitutions that can alter the function of the encoded protein. We evaluated the associations of six nonsynonymous SNPs in the MMP3, MMP8, and MMP9 genes with skin cancer risk in a nested case-control study of Caucasians within the Nurses' Health Study among 218 melanoma cases, 285 squamous cell carcinoma (SCC) cases, 300 basal cell carcinoma (BCC) cases, and 870 normal controls. We observed that the MMP9 Arg668Gln polymorphism was significantly associated with a decreased risk of SCC. Compared with the Arg/Arg group, the multivariate odds ratio was 0.67 (95% confidence interval, 0.47-0.97) for the Arg/Gln group and 0.21 (95% confidence interval, 0.05-0.97) for the Gln/Gln group (Ptrend = 0.004). We did not observe any association of this SNP with the risks of melanoma and basal cell carcinoma. No associations were found for other SNPs with skin cancer risk. This study provides evidence for the contribution of the MMP9 Arg668Gln to SCC development. (Cancer Epidemiol Biomarkers Prev 2008;17(12):3551–7)

Skin cancer is the most common form of cancer in the United States, with ∼1 million new cases per year (1). Among three major types of skin cancer, melanoma is the most serious form. The most common type of skin cancer is basal cell carcinoma (BCC), followed by squamous cell carcinoma (SCC). In addition to the carcinogenic effects of UV radiation, genetic factors also contribute to the development of skin cancer (2-4).

Degradation and remodeling of the extracellular matrix (ECM) and basement membrane are essential steps in tumor progression. Matrix metalloproteinases (MMP) belong to a family of over 20 zinc-dependent enzymes that degrade various protein components of the ECM. Overexpression of these enzymes leads to increased proteolytic degradation of ECM, resulting in the invasion and metastasis of many types of tumors (5-9). Based on their structure, substrate specificity, and cellular localization, MMPs can be divided into five subgroups: collagenases (including MMP8), gelatinases (including MMP9), stromelysins (including MMP3), matrilysins, and membrane-type MMPs (10). Due to their crucial roles in tumor development, MMPs have become drug targets for cancer therapy in some clinical trials (11, 12).

Single nucleotide polymorphisms (SNPs) are the most common genetic variation in the human genome, and nonsynonymous SNPs resulting in amino acid substitutions can alter the function of the encoded protein. It has been reported that MMP genetic variants including nonsynonymous polymorphisms are associated with the development of different types of tumors, such as melanoma, breast, lung, bladder, and colorectal tumors (13-17).

We evaluated the associations of six nonsynonymous SNPs in the MMP3, MMP8, and MMP9 genes with the risks of three skin cancer types (melanoma, SCC, and BCC) in a nested case-control study of Caucasians within the Nurses' Health Study.

Study Population

The Nurses' Health Study was established in 1976, when 121,700 female registered nurses between the ages of 30 and 55 y residing in 11 larger U.S. states completed and returned the initial self-administered questionnaire on their medical histories and baseline health-related exposures, forming the basis for the Nurses' Health Study cohort. Updated information has been obtained by questionnaires every 2 y. From May 1989 through September 1990, we collected blood samples from 32,826 participants in the Nurses' Health Study cohort. The distribution of risk factors for skin cancer in the subcohort of those who donated blood samples was similar to that in the overall cohort (4). Eligible cases in this study consisted of women with incident skin cancer from the subcohort who had given a blood specimen, including SCC and BCC cases with a diagnosis any time after blood collection up to June 1, 1998, and melanoma cases up to June 1, 2000, who had no previously diagnosed skin cancer. A common control series was randomly selected from participants who gave a blood sample and were free of diagnosed skin cancer up to and including the questionnaire cycle during which the case was diagnosed. One or two controls were matched to each case by year of birth (±1 y). All subjects were the U.S. non-Hispanic Caucasian women in this study. The participation rates of cases and controls were 92% and 89%, respectively. The nested case-control study consisted of 218 incident melanoma cases, 285 incident SCC cases, a sample of 300 BCC cases from the large number of incident cases, and 870 age-matched controls. The study protocol was approved by the Committee on Use of Human Subjects of the Brigham and Women's Hospital, Boston, MA.

Exposure Data

We obtained information regarding skin cancer risk factors from the prospective biennial questionnaires and a retrospective supplementary questionnaire. Information on natural hair color at age 20 y and childhood and adolescent tanning tendency were collected in the 1982 prospective questionnaire. Ethnic group was ascertained in the 1992 questionnaire. In the skin cancer nested case-control study, natural skin color and other sun exposure–related information were collected by the retrospective supplementary questionnaire in 2002. In addition, the 11 states of residence of cohort members at baseline were grouped into three regions: Northeast (Connecticut, Massachusetts, Maryland, New Jersey, New York, and Pennsylvania), Northcentral (Michigan and Ohio), and West and South (California, Texas, and Florida). To estimate sunlight exposure for each subject, an UV database for 50 U.S. states was developed. The database used reports from the Climatic Atlas of the United States, which reported mean daily solar radiation (in Langleys) at the earth's surface for weather stations around the country (18). The records of average annual solar radiation for January and July were extracted to represent winter and summer radiation, respectively. The mean solar radiation for each individual's past (at different age categories) and current residences was derived from the UV values measured at the nearest weather station. Both summer and winter radiation indices were developed for the residence of each age category. A cumulative lifetime sun exposure was developed by combining the residence-linked UV value and hours spent outdoors at difference age categories obtained from the retrospective supplementary questionnaire. This derived variable was significantly associated with skin cancer risk (19).

Laboratory Assays

We genotyped 6 SNPs including MMP3 Lys45Glu (rs679620), MMP8 Lys87Glu (rs1940475), MMP9 Ala20Val (rs1805088), MMP9 Gln279Arg (rs17576), MMP9 Pro574Arg (rs2250889), and MMP9 Arg668Gln (rs2274756) by the 5 nuclease assay (TaqMan) in 384-well format, using the ABI PRISM 7900HT Sequence Detection System (Applied Biosystems). TaqMan primers and probes were designed using the Primer Express Oligo Design software v2.0 (ABI PRISM). Due to assay failure, we genotyped rs6032619 as a surrogate for the MMP9 rs2250889 (r2 = 1.0). We calculated the r2 value between the rs6032619 and the rs2250889 using the HapMap project 90 (30 trios) Caucasian samples from a U.S. Utah population with Northern and Western European ancestry collected in 1980 by the Centre d'Etude du Polymorphisme Humain (20). Laboratory personnel were blinded to case-control status, and 10% blinded quality control samples (duplicate samples) were inserted to validate genotyping procedures; concordance for the blinded quality control samples was 100%. Primers, probes, and conditions for genotyping assays are available upon request.

Statistical Methods

We used the χ2 test to assess whether the genotypes for all six SNPs were in Hardy-Weinberg equilibrium among the controls. We evaluated the association between each genotype and skin cancer risk using unconditional logistic regression. An additive model was used to calculate the P value for trend on skin cancer risk according to an ordinal coding for genotype (0, 1, or 2 copies of SNP minor allele). We compared each type of skin cancer with the common control series to increase the statistical power. For the four MMP9 SNPs, haplotype frequencies and expected haplotype counts for each individual were estimated using a simple expectation-maximization algorithm, as implemented in SAS PROC HAPLOTYPE. The analyses of possible associations between haplotypes and skin cancer risk were done using the expectation-substitution technique (21). All statistical analyses were two-sided and carried out using SAS V9.1 (SAS Institute).

Computational Prediction for Functional Significance of Nonsynonymous SNPs

We used five bioinformatic tools, (a) BLOSUM62 (22), (b) PMut (23), (c) PolyPhen (24), (d) SIFT (25), and (e) SNPs3D (26), to predict the functional significance of the six nonsynonymous polymorphisms. These methods are described in the Supplementary Text.

Descriptive Characteristics of Skin Cancer Cases and Controls

Basic characteristics of cases and controls in this study are presented in Table 1. Detailed description was published previously (27). In brief, at the beginning of the follow-up of this nested case-control study, the nurses were between ages 43 and 68 y (mean age, 58.7 y). The mean ages at diagnosis for incident melanoma, SCC, and BCC cases were 63.4, 64.7, and 64.0 y, respectively. Skin cancer cases were more likely to possess red hair color and fair skin color. The childhood tanning ability of cases was less than that of controls. Women in the West and South regions were more likely to be diagnosed with SCC or BCC compared with those in Northeast. A family history of skin cancer was a risk factor for all three types of skin cancer. Those with skin cancers were more likely to have used sunlamps or attended tanning salons. Those with skin cancers had higher cumulative sun exposure while wearing a bathing suit and more lifetime severe sunburns that blistered. For melanoma, there were 20.4% cases on the head and neck, 29.1% on the trunk, 23.0% on the arm and hand, and 27.6% on the leg and foot. For SCC, there were 44.5% cases on the head and neck, 17.4% on the trunk, 20.0% on the arm and hand, and 18.1% on the leg and foot.

Table 1.

Characteristics of skin cancer cases and controls in the nested case-control study

CharacteristicControls (n = 870)Melanoma cases (n = 218)SCC cases (n = 285)BCC cases (n = 300)
Age at diagnosis (mean, y) 64.5 63.4 64.7 64.0 
Natural hair color at age 20 (%)     
    Black or dark brown 43.9 31.5 41.3 30.3 
    Light brown 40.0 42.5 34.6 45.7 
    Blonde 12.0 15.5 16.8 18.0 
    Red 2.9 10.5 5.2 4.7 
Natural skin color (%)     
    Fair 40.0 57.1 54.6 53.0 
    Medium 36.7 25.6 32.2 31.3 
    Olive 4.8 0.9 1.8 1.3 
Tanning ability (%)     
    Practically none 8.3 13.8 13.9 11.7 
    Light tan 21.1 29.5 24.1 25.5 
    Average tan 47.3 40.0 47.4 46.5 
    Tan 23.3 16.7 14.7 16.3 
Geographic region at baseline (%)     
    Northeast 55.2 58.0 51.7 49.3 
    Northcentral 23.4 16.9 17.1 20.3 
    West and South 21.4 25.1 31.1 30.3 
Family history of skin cancer (%) 25.1 36.5 35.7 42.7 
Sunlamp use or tanning salon attendance (%) 10.0 19.2 14.3 14.7 
Highest tertile of cumulative sun exposure with a bathing suit (%) 33.4 53.3 46.1 42.6 
Number of lifetime severe sunburns (mean) 5.4 9.6 7.8 8.2 
CharacteristicControls (n = 870)Melanoma cases (n = 218)SCC cases (n = 285)BCC cases (n = 300)
Age at diagnosis (mean, y) 64.5 63.4 64.7 64.0 
Natural hair color at age 20 (%)     
    Black or dark brown 43.9 31.5 41.3 30.3 
    Light brown 40.0 42.5 34.6 45.7 
    Blonde 12.0 15.5 16.8 18.0 
    Red 2.9 10.5 5.2 4.7 
Natural skin color (%)     
    Fair 40.0 57.1 54.6 53.0 
    Medium 36.7 25.6 32.2 31.3 
    Olive 4.8 0.9 1.8 1.3 
Tanning ability (%)     
    Practically none 8.3 13.8 13.9 11.7 
    Light tan 21.1 29.5 24.1 25.5 
    Average tan 47.3 40.0 47.4 46.5 
    Tan 23.3 16.7 14.7 16.3 
Geographic region at baseline (%)     
    Northeast 55.2 58.0 51.7 49.3 
    Northcentral 23.4 16.9 17.1 20.3 
    West and South 21.4 25.1 31.1 30.3 
Family history of skin cancer (%) 25.1 36.5 35.7 42.7 
Sunlamp use or tanning salon attendance (%) 10.0 19.2 14.3 14.7 
Highest tertile of cumulative sun exposure with a bathing suit (%) 33.4 53.3 46.1 42.6 
Number of lifetime severe sunburns (mean) 5.4 9.6 7.8 8.2 

NOTE: The percentages may not sum to 100 due to rounding.

Associations of the Six SNPs in the MMP3, MMP8, and MMP9 Genes with Skin Cancer Risk

The distributions of genotypes for 6 SNPs evaluated in this study were in Hardy-Weinberg equilibrium among controls (P > 0.05).

We evaluated the main effect of each polymorphism across the three types of skin cancer (Table 2). We observed that the Gln668 allele of the MMP9 Arg668Gln polymorphism was significantly associated with a decreased risk of SCC. Compared with the Arg/Arg group, the multivariate odds ratio (OR) was 0.67 [95% confidence interval (95% CI), 0.47-0.97] for the Arg/Gln group and 0.21 (95% CI, 0.05-0.97) for the Gln/Gln group (Ptrend = 0.004). The multivariate estimates were essentially the same as the age-adjusted ones.

Table 2.

Six SNPs in the MMP3, MMP8, and MMP9 genes and skin cancer risk

Melanoma
SCC
BCC
MMP3 Lys45GluControls (%)Cases (%)Multivariate OR*Multivariate ORCases (%)Multivariate OR*Multivariate ORCases (%)Multivariate OR*Multivariate OR
AA 211 (25.2) 45 (21.3) 1.00 1.00 74 (27.2) 1.00 1.00 77 (25.9) 1.00 1.00 
AG 428 (51.1) 120 (56.9) 1.31 (0.89-1.92) 1.25 (0.83-1.88) 141 (51.8) 0.94 (0.68-1.31) 0.92 (0.65-1.30) 161 (54.2) 1.04 (0.75-1.43) 1.02 (0.73-1.44) 
GG 198 (23.7) 46 (21.8) 1.08 (0.68-1.70) 0.95 (0.58-1.55) 57 (21.0) 0.82 (0.55-1.22) 0.81 (0.53-1.23) 59 (19.9) 0.83 (0.56-1.22) 0.79 (0.52-1.19) 
Additive OR — — 1.04 (0.83-1.29) 0.97 (0.77-1.23) — 0.91 (0.75-1.10) 0.90 (0.73-1.11) — 0.92 (0.76-1.11) 0.89 (0.73-1.10) 
Ptrend — — 0.74 0.83 — 0.33 0.32 — 0.37 0.28 
AG or GG 626 (74.8) 166 (78.7) 1.24 (0.86-1.78) 1.15 (0.78-1.70) 198 (72.8) 0.90 (0.66-1.23) 0.88 (0.64-1.23) 220 (74.1) 0.97 (0.72-1.31) 0.95 (0.69-1.31) 
           
MMP8 Lys87Glu
 
Controls (%)
 
Cases (%)
 
Multivariate OR*
 
Multivariate OR
 
Cases (%)
 
Multivariate OR*
 
Multivariate OR
 
Cases (%)
 
Multivariate OR*
 
Multivariate OR
 
TT 222 (26.8) 56 (27.2) 1.00 1.00 74 (28.0) 1.00 1.00 64 (22.1) 1.00 1.00 
TC 409 (49.5) 104 (50.5) 1.00 (0.69-1.43) 1.08 (0.73-1.60) 117 (44.3) 0.85 (0.61-1.19) 0.87 (0.61-1.24) 146 (50.3) 1.25 (0.89-1.75) 1.34 (0.94-1.91) 
CC 196 (23.7) 46 (22.3) 0.93 (0.60-1.43) 1.01 (0.63-1.62) 73 (27.7) 1.12 (0.77-1.63) 1.22 (0.82-1.82) 80 (27.6) 1.43 (0.98-2.10) 1.52 (1.01-2.29) 
Additive OR — — 0.96 (0.78-1.20) 1.01 (0.80-1.27) — 1.06 (0.87-1.28) 1.10 (0.90-1.35) — 1.19 (0.99-1.44) 1.23 (1.01-1.51) 
Ptrend — — 0.74 0.94 — 0.58 0.36 — 0.07 0.04 
TT or CC 605 (73.2) 150 (72.8) 0.97 (0.69-1.37) 1.06 (0.73-1.54) 190 (72.0) 0.94 (0.69-1.28) 0.98 (0.71-1.36) 226 (77.9) 1.31 (0.95-1.80) 1.40 (1.00-1.96) 
           
MMP9 Ala20Val
 
Controls (%)
 
Cases (%)
 
Multivariate OR*
 
Multivariate OR
 
Cases (%)
 
Multivariate OR*
 
Multivariate OR
 
Cases (%)
 
Multivariate OR*
 
Multivariate OR
 
CC 832 (99.0) 210 (99.5) 1.00 1.00 274 (100.0) 1.00 1.00 289 (98.3) 1.00 1.00 
CT 8 (1.0) 1 (0.5) 0.50 (0.06-4.08) 0.24 (0.03-2.04) 0 (0.0) — — 5 (1.7) 1.85 (0.60-5.74) 1.24 (0.38-4.02) 
TT 0 (0.0) 0 (0.0) — — 0 (0.0) — — 0 (0.0) — — 
           
MMP9 Gln279Arg
 
Controls (%)
 
Cases (%)
 
Multivariate OR*
 
Multivariate OR
 
Cases (%)
 
Multivariate OR*
 
Multivariate OR
 
Cases (%)
 
Multivariate OR*
 
Multivariate OR
 
AA 336 (40.4) 87 (42.2) 1.00 1.00 117 (42.9) 1.00 1.00 115 (39.2) 1.00 1.00 
AG 389 (46.8) 93 (45.1) 0.91 (0.66-1.27) 0.83 (0.58-1.18) 131 (48.0) 0.97 (0.73-1.30) 0.88 (0.65-1.19) 141 (48.1) 1.07 (0.80-1.42) 1.03 (0.76-1.40) 
GG 107 (12.9) 26 (12.6) 0.94 (0.57-1.53) 0.89 (0.52-1.50) 25 (9.2) 0.67 (0.41-1.09) 0.57 (0.34-0.95) 37 (12.6) 1.02 (0.66-1.57) 0.95 (0.60-1.49) 
Additive OR — — 0.95 (0.76-1.19) 0.91 (0.71-1.16) — 0.87 (0.71-1.07) 0.80 (0.64-0.99) — 1.03 (0.84-1.25) 0.99 (0.80-1.22) 
Ptrend — — 0.66 0.44 — 0.20 0.04 — 0.79 0.92 
AG or GG 496 (59.6) 119 (57.8) 0.92 (0.67-1.25) 0.84 (0.60-1.18) 156 (57.1) 0.91 (0.69-1.20) 0.81 (0.60-1.08) 178 (60.8) 1.06 (0.81-1.39) 1.01 (0.76-1.35) 
           
MMP9 Pro574Arg
 
Controls (%)
 
Cases (%)
 
Multivariate OR*
 
Multivariate OR
 
Cases (%)
 
Multivariate OR*
 
Multivariate OR
 
Cases (%)
 
Multivariate OR*
 
Multivariate OR
 
CC 755 (93.7) 189 (96.9) 1.00 1.00 250 (93.6) 1.00 1.00 263 (94.6) 1.00 1.00 
CG 51 (6.3) 6 (3.1) 0.50 (0.21-1.19) 0.60 (0.24-1.50) 17 (6.4) 0.99 (0.56-1.76) 1.13 (0.62-2.05) 15 (5.4) 0.83 (0.46-1.51) 0.95 (0.51-1.77) 
GG 0 (0.0) 0 (0.0) — — 0 (0.0) — — 0 (0.0) — — 
           
MMP9 Arg668Gln
 
Controls (%)
 
Cases (%)
 
Multivariate OR*
 
Multivariate OR
 
Cases (%)
 
Multivariate OR*
 
Multivariate OR
 
Cases (%)
 
Multivariate OR*
 
Multivariate OR
 
GG 623 (74.1) 144 (68.6) 1.00 1.00 218 (80.1) 1.00 1.00 219 (75.0) 1.00 1.00 
GA 196 (23.3) 59 (28.1) 1.32 (0.93-1.86) 1.28 (0.88-1.85) 52 (19.1) 0.75 (0.54-1.06) 0.67 (0.47-0.97) 61 (20.9) 0.88 (0.64-1.23) 0.84 (0.59-1.18) 
AA 22 (2.6) 7 (3.3) 1.44 (0.60-3.45) 1.27 (0.49-3.28) 2 (0.7) 0.26 (0.06-1.11) 0.21 (0.05-0.97) 12 (4.1) 1.55 (0.75-3.19) 1.61 (0.74-3.46) 
Additive OR — — 1.27 (0.96-1.69) 1.22 (0.90-1.65) — 0.69 (0.51-0.94) 0.62 (0.45-0.86) — 1.02 (0.79-1.32) 1.00 (0.76-1.31) 
Ptrend — — 0.10 0.20 — 0.019 0.004 — 0.88 0.97 
GA or AA 218 (25.9) 66 (31.4) 1.33 (0.95-1.85) 1.28 (0.89-1.82) 54 (19.9) 0.71 (0.50-0.99) 0.63 (0.44-0.89) 73 (25.0) 0.95 (0.70-1.29) 0.91 (0.66-1.26) 
Melanoma
SCC
BCC
MMP3 Lys45GluControls (%)Cases (%)Multivariate OR*Multivariate ORCases (%)Multivariate OR*Multivariate ORCases (%)Multivariate OR*Multivariate OR
AA 211 (25.2) 45 (21.3) 1.00 1.00 74 (27.2) 1.00 1.00 77 (25.9) 1.00 1.00 
AG 428 (51.1) 120 (56.9) 1.31 (0.89-1.92) 1.25 (0.83-1.88) 141 (51.8) 0.94 (0.68-1.31) 0.92 (0.65-1.30) 161 (54.2) 1.04 (0.75-1.43) 1.02 (0.73-1.44) 
GG 198 (23.7) 46 (21.8) 1.08 (0.68-1.70) 0.95 (0.58-1.55) 57 (21.0) 0.82 (0.55-1.22) 0.81 (0.53-1.23) 59 (19.9) 0.83 (0.56-1.22) 0.79 (0.52-1.19) 
Additive OR — — 1.04 (0.83-1.29) 0.97 (0.77-1.23) — 0.91 (0.75-1.10) 0.90 (0.73-1.11) — 0.92 (0.76-1.11) 0.89 (0.73-1.10) 
Ptrend — — 0.74 0.83 — 0.33 0.32 — 0.37 0.28 
AG or GG 626 (74.8) 166 (78.7) 1.24 (0.86-1.78) 1.15 (0.78-1.70) 198 (72.8) 0.90 (0.66-1.23) 0.88 (0.64-1.23) 220 (74.1) 0.97 (0.72-1.31) 0.95 (0.69-1.31) 
           
MMP8 Lys87Glu
 
Controls (%)
 
Cases (%)
 
Multivariate OR*
 
Multivariate OR
 
Cases (%)
 
Multivariate OR*
 
Multivariate OR
 
Cases (%)
 
Multivariate OR*
 
Multivariate OR
 
TT 222 (26.8) 56 (27.2) 1.00 1.00 74 (28.0) 1.00 1.00 64 (22.1) 1.00 1.00 
TC 409 (49.5) 104 (50.5) 1.00 (0.69-1.43) 1.08 (0.73-1.60) 117 (44.3) 0.85 (0.61-1.19) 0.87 (0.61-1.24) 146 (50.3) 1.25 (0.89-1.75) 1.34 (0.94-1.91) 
CC 196 (23.7) 46 (22.3) 0.93 (0.60-1.43) 1.01 (0.63-1.62) 73 (27.7) 1.12 (0.77-1.63) 1.22 (0.82-1.82) 80 (27.6) 1.43 (0.98-2.10) 1.52 (1.01-2.29) 
Additive OR — — 0.96 (0.78-1.20) 1.01 (0.80-1.27) — 1.06 (0.87-1.28) 1.10 (0.90-1.35) — 1.19 (0.99-1.44) 1.23 (1.01-1.51) 
Ptrend — — 0.74 0.94 — 0.58 0.36 — 0.07 0.04 
TT or CC 605 (73.2) 150 (72.8) 0.97 (0.69-1.37) 1.06 (0.73-1.54) 190 (72.0) 0.94 (0.69-1.28) 0.98 (0.71-1.36) 226 (77.9) 1.31 (0.95-1.80) 1.40 (1.00-1.96) 
           
MMP9 Ala20Val
 
Controls (%)
 
Cases (%)
 
Multivariate OR*
 
Multivariate OR
 
Cases (%)
 
Multivariate OR*
 
Multivariate OR
 
Cases (%)
 
Multivariate OR*
 
Multivariate OR
 
CC 832 (99.0) 210 (99.5) 1.00 1.00 274 (100.0) 1.00 1.00 289 (98.3) 1.00 1.00 
CT 8 (1.0) 1 (0.5) 0.50 (0.06-4.08) 0.24 (0.03-2.04) 0 (0.0) — — 5 (1.7) 1.85 (0.60-5.74) 1.24 (0.38-4.02) 
TT 0 (0.0) 0 (0.0) — — 0 (0.0) — — 0 (0.0) — — 
           
MMP9 Gln279Arg
 
Controls (%)
 
Cases (%)
 
Multivariate OR*
 
Multivariate OR
 
Cases (%)
 
Multivariate OR*
 
Multivariate OR
 
Cases (%)
 
Multivariate OR*
 
Multivariate OR
 
AA 336 (40.4) 87 (42.2) 1.00 1.00 117 (42.9) 1.00 1.00 115 (39.2) 1.00 1.00 
AG 389 (46.8) 93 (45.1) 0.91 (0.66-1.27) 0.83 (0.58-1.18) 131 (48.0) 0.97 (0.73-1.30) 0.88 (0.65-1.19) 141 (48.1) 1.07 (0.80-1.42) 1.03 (0.76-1.40) 
GG 107 (12.9) 26 (12.6) 0.94 (0.57-1.53) 0.89 (0.52-1.50) 25 (9.2) 0.67 (0.41-1.09) 0.57 (0.34-0.95) 37 (12.6) 1.02 (0.66-1.57) 0.95 (0.60-1.49) 
Additive OR — — 0.95 (0.76-1.19) 0.91 (0.71-1.16) — 0.87 (0.71-1.07) 0.80 (0.64-0.99) — 1.03 (0.84-1.25) 0.99 (0.80-1.22) 
Ptrend — — 0.66 0.44 — 0.20 0.04 — 0.79 0.92 
AG or GG 496 (59.6) 119 (57.8) 0.92 (0.67-1.25) 0.84 (0.60-1.18) 156 (57.1) 0.91 (0.69-1.20) 0.81 (0.60-1.08) 178 (60.8) 1.06 (0.81-1.39) 1.01 (0.76-1.35) 
           
MMP9 Pro574Arg
 
Controls (%)
 
Cases (%)
 
Multivariate OR*
 
Multivariate OR
 
Cases (%)
 
Multivariate OR*
 
Multivariate OR
 
Cases (%)
 
Multivariate OR*
 
Multivariate OR
 
CC 755 (93.7) 189 (96.9) 1.00 1.00 250 (93.6) 1.00 1.00 263 (94.6) 1.00 1.00 
CG 51 (6.3) 6 (3.1) 0.50 (0.21-1.19) 0.60 (0.24-1.50) 17 (6.4) 0.99 (0.56-1.76) 1.13 (0.62-2.05) 15 (5.4) 0.83 (0.46-1.51) 0.95 (0.51-1.77) 
GG 0 (0.0) 0 (0.0) — — 0 (0.0) — — 0 (0.0) — — 
           
MMP9 Arg668Gln
 
Controls (%)
 
Cases (%)
 
Multivariate OR*
 
Multivariate OR
 
Cases (%)
 
Multivariate OR*
 
Multivariate OR
 
Cases (%)
 
Multivariate OR*
 
Multivariate OR
 
GG 623 (74.1) 144 (68.6) 1.00 1.00 218 (80.1) 1.00 1.00 219 (75.0) 1.00 1.00 
GA 196 (23.3) 59 (28.1) 1.32 (0.93-1.86) 1.28 (0.88-1.85) 52 (19.1) 0.75 (0.54-1.06) 0.67 (0.47-0.97) 61 (20.9) 0.88 (0.64-1.23) 0.84 (0.59-1.18) 
AA 22 (2.6) 7 (3.3) 1.44 (0.60-3.45) 1.27 (0.49-3.28) 2 (0.7) 0.26 (0.06-1.11) 0.21 (0.05-0.97) 12 (4.1) 1.55 (0.75-3.19) 1.61 (0.74-3.46) 
Additive OR — — 1.27 (0.96-1.69) 1.22 (0.90-1.65) — 0.69 (0.51-0.94) 0.62 (0.45-0.86) — 1.02 (0.79-1.32) 1.00 (0.76-1.31) 
Ptrend — — 0.10 0.20 — 0.019 0.004 — 0.88 0.97 
GA or AA 218 (25.9) 66 (31.4) 1.33 (0.95-1.85) 1.28 (0.89-1.82) 54 (19.9) 0.71 (0.50-0.99) 0.63 (0.44-0.89) 73 (25.0) 0.95 (0.70-1.29) 0.91 (0.66-1.26) 
*

Unconditional logistic regression adjusted for the age.

Unconditional logistic regression adjusted for age, constitutional susceptibility score (tertiles), family history of skin cancer (yes/no), the number of lifetime severe sunburns that blistered (none, 1-5, 6-11, >11), sunlamp use or tanning salon attendance (yes/no), cumulative sun exposure while wearing a bathing suit (tertiles), and geographic region.

In the multivariate analyses controlling for age and skin cancer risk factors, we observed that MMP8 Lys87Glu was associated with an increased risk of BCC (OR, 1.23; 95% CI, 1.01-1.51), and MMP9 Gln279Arg was associated with a decreased risk of SCC (OR, 0.80; 95% CI, 0.64-0.99). However, these associations were not statistically significant in the age-adjusted analyses. The ORs (95% CI) of MMP8 Lys87Glu for BCC risk and MMP9 Gln279Arg for SCC risk were 1.19 (0.99-1.44) and 0.87 (0.71-1.07), respectively. We did not observe any significant associations between any of these six SNPs and melanoma risk.

Haplotypes for the MMP9 Gene and Skin Cancer Risk

We did a global test to evaluate the difference in haplotype frequencies between cases and controls (Table 3). We found a significant difference in haplotype frequency for SCC (Pglobal test = 0.02). The haplotype that carried the Gln279Arg and Arg668Gln variant alleles was significantly associated with a decreased risk of SCC after adjusting for age (OR, 0.62; 95% CI, 0.44-0.87). This association remained significant in the multivariate model (Pglobal test = 0.004; OR, 0.56; 95% CI, 0.40-0.80). This association was consistent with the result of single SNP analysis presented in Table 2.

Table 3.

Haplotypes for four SNPs in the MMP9 gene and skin cancer risk

MMP9
Controls
Melanoma
SCC
BCC
Cases
Cases
Cases
ABCDn (%)n (%)n (%)n (%)
925 (60.9) 229 (63.3) 323 (65.1) 318 (59.6) 
    Multivariate OR 1.00 1.00 1.00 
318 (20.9) 64 (17.7) 109 (22.0) 115 (21.5) 
    Multivariate OR 0.80 (0.58-1.09) 0.97 (0.75-1.26) 1.08 (0.83-1.39) 
222 (14.6) 63 (17.4) 47 (9.5) 82 (15.4) 
    Multivariate OR 1.15 (0.84-1.57) 0.62 (0.44-0.87) 1.08 (0.82-1.42) 
49 (3.2) 5 (1.4) 17 (3.4) 14 (2.6) 
    Multivariate OR 0.43 (0.17-1.10) 0.99 (0.55-1.76) 0.82 (0.44-1.52) 
Rare < 1%    6 (0.4) 1 (0.3) 0 (0.0) 5 (0.9) 
    Multivariate OR 0.76 (0.09-6.48) — 2.60 (0.77-8.74) 
A: Ala20Val; B: Gln279Arg; C: Pro574Arg; D: Arg668Gln        
MMP9
Controls
Melanoma
SCC
BCC
Cases
Cases
Cases
ABCDn (%)n (%)n (%)n (%)
925 (60.9) 229 (63.3) 323 (65.1) 318 (59.6) 
    Multivariate OR 1.00 1.00 1.00 
318 (20.9) 64 (17.7) 109 (22.0) 115 (21.5) 
    Multivariate OR 0.80 (0.58-1.09) 0.97 (0.75-1.26) 1.08 (0.83-1.39) 
222 (14.6) 63 (17.4) 47 (9.5) 82 (15.4) 
    Multivariate OR 1.15 (0.84-1.57) 0.62 (0.44-0.87) 1.08 (0.82-1.42) 
49 (3.2) 5 (1.4) 17 (3.4) 14 (2.6) 
    Multivariate OR 0.43 (0.17-1.10) 0.99 (0.55-1.76) 0.82 (0.44-1.52) 
Rare < 1%    6 (0.4) 1 (0.3) 0 (0.0) 5 (0.9) 
    Multivariate OR 0.76 (0.09-6.48) — 2.60 (0.77-8.74) 
A: Ala20Val; B: Gln279Arg; C: Pro574Arg; D: Arg668Gln        

NOTE: 0, common allele; 1, rare allele. Logistic regression adjusted for age. P value for global test for melanoma is 0.12; for SCC is 0.02; and for BCC is 0.53.

Computational Prediction of the Functional Significance for the Six SNPs

The results of the computational prediction of the functional significance of the six polymorphisms are summarized in Table 4. Results based on BLOSUM62 show that only the MMP9 Pro574Arg polymorphism has a negative score, which is indicative of a nonconservative change (28). This suggests that MMP9 Pro574Arg is likely to be deleterious. Results based on PMut show that MMP9 Pro574Arg and MMP9 Arg668Gln had pathogenicity prediction scores of >0.5, suggesting that these 2 nonsynonymous SNPs may be pathologic, whereas the other SNPs were predicted to be neutral. Only the MMP9 Pro574Arg polymorphism displayed a reliability index of >5. Results based on PolyPhen show that, except for the MMP9 Gln279Arg polymorphism, which had a PSIC score Δ of 1.537 (“possibly damaging”), all other SNPs had PSIC score Δ values <1.0, (“benign”). SIFT showed that the 6 SNPs had TI scores ranging from 0.35 to 1.00, with MSCS values below 3.25 across the board; therefore, all the six SNPs are considered “tolerated.” SNPs3D prediction scores, such as svm profile score and svm structure score, were calculated for each SNP. All six SNPs had positive svm profile scores, showing that they are not highly conserved in the respective protein families. Only the MMP9 Arg668Gln polymorphism received a negative svm structure score and was predicted to cause losses of both hydrogen bond and salt bridge; therefore, this SNP may have a deleterious effect on protein structural stability.

Table 4.

Computational prediction of the functional significance of six SNPs in the MMP3, MMP8, and MMP9 genes

SNPBLOSUM62 score*PMut pathogenicity prediction scorePolyPhen PSIC score DSIFT TI score§SNPs3D prediction scores (molecular effects)
MMP3 Lys45Glu 1.00 (evolutionarily more acceptable) 0.4299 [1] (neutral) 0.143 (benign) 1.00 [1.99] (tolerated) Svm profile = 1.79, svm structure = 0.79 (on the protein surface; electrostatic repulsion) 
MMP8 Lys87Glu 1.00 (evolutionarily more acceptable) 0.4159 [1] (neutral) 0.793 (benign) 0.80 [1.98] (tolerated) Svm profile = 0.01, svm structure = 1.26 (on the protein surface) 
MMP9 Ala20Val 0.00 (evolutionarily more acceptable) 0.3748 [2] (neutral) 0.421 (benign) 0.56 [2.55] (tolerated) Svm profile = 1.67 
MMP9 Gln279Arg 1.00 (evolutionarily more acceptable) 0.2577 [4] (neutral) 1.537 (possibly damaging) 0.38 [2.61] (Tolerated) Svm profile = 2.28, svm structure = 1.24 
MMP9 Pro574Arg −2.00 (evolutionarily less acceptable) 0.8465 [6] (pathologic) 0.764 (benign) 0.47 [2.29] (tolerated) Svm profile = 2.40, svm structure = 0.79 (saltbridge lost) 
MMP9 Arg668Gln 1.00 (evolutionarily more acceptable) 0.5389 [0] (pathologic) 0.574 (benign) 0.52 [2.28] (tolerated) Svm profile = 1.90, svm structure = −0.68 (hydrogen bond lost; saltbridge lost) 
SNPBLOSUM62 score*PMut pathogenicity prediction scorePolyPhen PSIC score DSIFT TI score§SNPs3D prediction scores (molecular effects)
MMP3 Lys45Glu 1.00 (evolutionarily more acceptable) 0.4299 [1] (neutral) 0.143 (benign) 1.00 [1.99] (tolerated) Svm profile = 1.79, svm structure = 0.79 (on the protein surface; electrostatic repulsion) 
MMP8 Lys87Glu 1.00 (evolutionarily more acceptable) 0.4159 [1] (neutral) 0.793 (benign) 0.80 [1.98] (tolerated) Svm profile = 0.01, svm structure = 1.26 (on the protein surface) 
MMP9 Ala20Val 0.00 (evolutionarily more acceptable) 0.3748 [2] (neutral) 0.421 (benign) 0.56 [2.55] (tolerated) Svm profile = 1.67 
MMP9 Gln279Arg 1.00 (evolutionarily more acceptable) 0.2577 [4] (neutral) 1.537 (possibly damaging) 0.38 [2.61] (Tolerated) Svm profile = 2.28, svm structure = 1.24 
MMP9 Pro574Arg −2.00 (evolutionarily less acceptable) 0.8465 [6] (pathologic) 0.764 (benign) 0.47 [2.29] (tolerated) Svm profile = 2.40, svm structure = 0.79 (saltbridge lost) 
MMP9 Arg668Gln 1.00 (evolutionarily more acceptable) 0.5389 [0] (pathologic) 0.574 (benign) 0.52 [2.28] (tolerated) Svm profile = 1.90, svm structure = −0.68 (hydrogen bond lost; saltbridge lost) 
*

A Score < 0 signals an evolutionarily less acceptable substitution, and a Score ≥ 0 signals an evolutionarily more acceptable substitution.

A Pathogenicity Prediction Score of >0.5 signals a pathologic mutation. The Reliability Index is shown in the square brackets for each SNP. The prediction becomes more reliable as the Reliability Index becomes higher.

The functional significance prediction PSIC score “0.00-0.99” and “1.50-1.99” are indicative of “Benign” and “Possibly damaging,” respectively.

§

TI, tolerance index; the median sequence conservation score is shown in the square brackets for each SNP. An median sequence conservation score of >3.25 is indicative of a low confidence because of an insufficiency of homologous sequences.

For both svm profile and svm structure scores, a negative value is indicative of a deleterious substitution, whereas a positive value is indicative of a neutral change, with greater absolute values indicating stronger confidence in the prediction.

The MMPs play an important role in proteolytic degradation and remodeling of various ECM components in both physiologic and pathologic situations, including tissue repair, angiogenesis, and tumor cell invasion and metastasis. Under normal physiologic conditions, MMPs are expressed at very low levels, whereas overexpression and activation have been found in many types of tumors (8, 29-31). Several previous studies examined the associations between the MMP polymorphisms and the development of tumors, such as lung cancer, breast cancer, and melanoma metastasis (11, 13-17).

Here, we assessed the associations between six nonsynonymous SNPs in the MMP3, MMP8, and MMP9 genes and the risk of three types of skin cancer. We observed a significant inverse association between the MMP9 Gln668 allele and SCC risk, mirrored in the haplotype analysis as well. MMP9 (gelatinase B) is produced by multiple cells including keratinocytes (32). Expression of the MMP9 gene was detected in malignantly transformed keratinocytes in SCC but not in normal keratinocytes (32). Furthermore, overexpression of MMP9 has been linked to the dyskeratotic foci of Bowen's disease, which is cutaneous SCC in situ (32-34). In addition, MMP9 protein expression was up-regulated in SCC of the skin when compared with BCC (35). In a transgenic mouse model, the absence of MMP9 was shown to delay angiogenesis and to reduce the frequency of SCC (33).

The results of the computational prediction of the functional significance for the SNPs in this study indicated the possible functional effects of the MMP9 nonsynonymous SNPs, such as Pro574Arg and Arg668Gln, on the MMP9 protein. We showed that the MMP9 Pro574Arg and Arg668Gln had two predictions of being less evolutionarily conservative or deleterious. The MMP9 Pro574Arg is predicted by BLOSUM62 to be evolutionarily less acceptable and by PMut to be pathologic to protein function; the MMP9 Arg668Gln is predicted by PMut to be pathologic and by SNPs3D to be deleterious to the protein structural stability. Similarly, a synthetic three-dimensional structure of MMP9 derived by Cotignola et al. (13) reflects the potential alteration of MMP9 protein function due to an amino acid change of nonsynonymous SNPs in the MMP9 gene including Arg668Gln, Pro574Arg, and Gln279Arg evaluated in this study. Presumably, the Gln allele of the MMP9 Arg668Gln is associated with reduced activity of the MMP9 protein, reducing the degradation of ECM. This could help explain our finding that women carrying the variant allele of MMP9 Arg668Gln have lower risk of SCC than those carrying the wild-type allele. On the other hand, compared with melanocytes and basal cells, keratinocytes are less tolerant of DNA damage and have a lower apoptotic threshold, which makes easier to enter the apoptosis pathway as a protective mechanism for SCC (36-38). This further supports the effect of MMP9 protein on SCC risk rather than melanoma or BCC risks because it has been proposed that MMP9 contributes to apoptosis (39-41). However, further biochemical studies are needed to confirm this because the role of MMPs in tumorigenesis is much more complex. Individual MMPs may play different roles in different cell types or in different states of transformation. For example, it has been shown that MMP9 can induce tumorigenesis and also inhibit it by enhancing angiostatin production (42-44).

For melanoma and BCC risks, we did not observe any suggestive associations with the six nonsynonymous SNPs we evaluated. Previous studies focused on the association between the SNPs in the MMP2, MMP3, and MMP9 genes and melanoma progression, showing significant inverse associations of the three MMP9 SNPs (Gln279Arg, Pro574Arg, and Arg668Gln) with metastasis of melanoma (13, 45). No studies have assessed the effects of SNPs in the MMP genes on the risks of SCC and BCC. To our knowledge, this is the first simultaneous evaluation of the associations between nonsynonymous SNPs in the MMP genes and the risks of three most common types of skin cancers.

In summary, we found an inverse association between the MMP9 Arg668Gln polymorphism and SCC risk. The nested case-control design, high follow-up rate, and high response rate for the retrospective supplementary questionnaire are among strengths of this study. The limitations of the study include the misclassification of the self-reported assessment of pigmentary phenotypes and sun exposure–related information. Our previous report showed that the retrospective assessment of risk factors was not likely to substantially bias the estimate of risk in this study (19). Because we selected incident cases after blood collection at the mean age of 58.7 years, the cases had an older age at diagnosis than that in the general population. The sample size in each cancer type in this study was modest, and additional studies are warranted to confirm these associations. Further biochemical studies are warranted to explore the biological characteristics of the MMP genetic variants.

No potential conflicts of interest were disclosed.

Grant support: NIH grant number CA122838 and CA132175.

Note: Supplementary data for this article are available at Cancer Epidemiology Biomakers and Prevention Online (http://cebp.aacrjournals.org/).

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.

We thank Dr. Hardeep Ranu and Pati Soule of the Dana-Farber/Harvard Cancer Center High-Throughput Polymorphism Detection Core for their laboratory assistance, and Carolyn Guo for her programming support, and the participants in the Nurses' Health Study for their dedication and commitment.

1
Howe HL, Wingo PA, Thun MJ, et al. Annual report to the nation on the status of cancer (1973 through 1998), featuring cancers with recent increasing trends.
J Natl Cancer Inst
2001
;
93
:
824
–42.
2
Brash DE. Sunlight and the onset of skin cancer.
Trends Genet
1997
;
13
:
410
–4.
3
de Gruijl FR, van Kranen HJ, Mullenders LH. UV-induced DNA damage, repair, mutations and oncogenic pathways in skin cancer.
J Photochem Photobiol B
2001
;
63
:
19
–27.
4
Han J, Kraft P, Colditz GA, Wong J, Hunter DJ. Melanocortin 1 receptor variants and skin cancer risk.
Int J Cancer
2006
;
119
:
1976
–84.
5
Bergers G, Coussens LM. Extrinsic regulators of epithelial tumor progression: metalloproteinases.
Curr Opin Genet Dev
2000
;
10
:
120
–7.
6
Coussens LM, Werb Z. Matrix metalloproteinases and the development of cancer.
Chem Biol
1996
;
3
:
895
–904.
7
Egeblad M, Werb Z. New functions for the matrix metalloproteinases in cancer progression.
Nat Rev Cancer
2002
;
2
:
161
–74.
8
Nagase H, Woessner JF, Jr. Matrix metalloproteinases.
J Biol Chem
1999
;
274
:
21491
–4.
9
Stetler-Stevenson WG, Aznavoorian S, Liotta LA. Tumor cell interactions with the extracellular matrix during invasion and metastasis.
Annu Rev Cell Biol
1993
;
9
:
541
–73.
10
Visse R, Nagase H. Matrix metalloproteinases and tissue inhibitors of metalloproteinases: structure, function, and biochemistry.
Circ Res
2003
;
92
:
827
–39.
11
Decock J, Paridaens R, Ye S. Genetic polymorphisms of matrix metalloproteinases in lung, breast and colorectal cancer.
Clin Genet
2008
;
73
:
197
–211.
12
Overall CM, Kleifeld O. Tumour microenvironment - opinion: validating matrix metalloproteinases as drug targets and anti-targets for cancer therapy.
Nat Rev Cancer
2006
;
6
:
227
–39.
13
Cotignola J, Reva B, Mitra N, et al. Matrix metalloproteinase-9 (MMP-9) polymorphisms in patients with cutaneous malignant melanoma.
BMC Med Genet
2007
;
8
:
10
.
14
Elander N, Soderkvist P, Fransen K. Matrix metalloproteinase (MMP) -1, -2, -3 and -9 promoter polymorphisms in colorectal cancer.
Anticancer Res
2006
;
26
:
791
–5.
15
Heist RS, Marshall AL, Liu G, et al. Matrix metalloproteinase polymorphisms and survival in stage I non-small cell lung cancer.
Clin Cancer Res
2006
;
12
:
5448
–53.
16
Kader AK, Liu J, Shao L, et al. Matrix metalloproteinase polymorphisms are associated with bladder cancer invasiveness.
Clin Cancer Res
2007
;
13
:
2614
–20.
17
Lei H, Hemminki K, Altieri A, et al. Promoter polymorphisms in matrix metalloproteinases and their inhibitors: few associations with breast cancer susceptibility and progression.
Breast Cancer Res Treat
2007
;
103
:
61
–9.
18
National Oceanic and Atmospheric Administration. Climatic Atlas of the United States. 1983.
19
Han J, Colditz GA, Hunter DJ. Risk factors for skin cancers: a nested case-control study within the Nurses' Health Study.
Int J Epidemiol
2006
;
35
:
1514
–21.
20
Dausset J, Cann H, Cohen D, Lathrop M, Lalouel JM, White R. Centre d'etude du polymorphisme humain (CEPH): collaborative genetic mapping of the human genome.
Genomics
1990
;
6
:
575
–7.
21
Kraft P, Cox DG, Paynter RA, Hunter D, De Vivo I. Accounting for haplotype uncertainty in matched association studies: a comparison of simple and flexible techniques.
Genet Epidemiol
2005
;
28
:
261
–72.
22
Henikoff S, Henikoff JG. Amino acid substitution matrices from protein blocks.
Proc Natl Acad Sci U S A
1992
;
89
:
10915
–9.
23
Ferrer-Costa C, Gelpi JL, Zamakola L, Parraga I, de la Cruz X, Orozco M. PMUT: a web-based tool for the annotation of pathological mutations on proteins.
Bioinformatics
2005
;
21
:
3176
–8.
24
Sunyaev SR, Eisenhaber F, Rodchenkov IV, Eisenhaber B, Tumanyan VG, Kuznetsov EN. PSIC: profile extraction from sequence alignments with position-specific counts of independent observations.
Protein Eng
1999
;
12
:
387
–94.
25
Ng PC, Henikoff S. SIFT: Predicting amino acid changes that affect protein function.
Nucleic Acids Res
2003
;
31
:
3812
–4.
26
Yue P, Melamud E, Moult J. SNPs3D: candidate gene and SNP selection for association studies.
BMC Bioinformatics
2006
;
7
:
166
.
27
Han J, Colditz GA, Liu JS, Hunter DJ. Genetic variation in XPD, sun exposure, and risk of skin cancer.
Cancer Epidemiol Biomarkers Prev
2005
;
14
:
1539
–44.
28
Cargill M, Altshuler D, Ireland J, et al. Characterization of single-nucleotide polymorphisms in coding regions of human genes.
Nat Genet
1999
;
22
:
231
–8.
29
Kahari VM, Saarialho-Kere U. Matrix metalloproteinases in skin.
Exp Dermatol
1997
;
6
:
199
–213.
30
McCawley LJ, Matrisian LM. Matrix metalloproteinases: multifunctional contributors to tumor progression.
Mol Med Today
2000
;
6
:
149
–56.
31
Vu TH, Werb Z. Matrix metalloproteinases: effectors of development and normal physiology.
Genes Dev
2000
;
14
:
2123
–33.
32
Kerkela E, Saarialho-Kere U. Matrix metalloproteinases in tumor progression: focus on basal and squamous cell skin cancer.
Exp Dermatol
2003
;
12
:
109
–25.
33
Coussens LM, Tinkle CL, Hanahan D, Werb Z. MMP-9 supplied by bone marrow-derived cells contributes to skin carcinogenesis.
Cell
2000
;
103
:
481
–90.
34
Kobayashi T, Onoda N, Takagi T, et al. Immunolocalizations of human gelatinase (type IV collagenase, MMP-9) and TIMP (tissue inhibitor of metalloproteinases) in normal epidermis and some epidermal tumors.
Arch Dermatol Res
1996
;
288
:
239
–44.
35
Dumas V, Kanitakis J, Charvat S, Euvrard S, Faure M, Claudy A. Expression of basement membrane antigens and matrix metalloproteinases 2 and 9 in cutaneous basal and squamous cell carcinomas.
Anticancer Res
1999
;
19
:
2929
–38.
36
Bowen AR, Hanks AN, Allen SM, Alexander A, Diedrich MJ, Grossman D. Apoptosis regulators and responses in human melanocytic and keratinocytic cells.
J Invest Dermatol
2003
;
120
:
48
–55.
37
Danno K, Horio T. Sunburn cell: factors involved in its formation.
Photochem Photobiol
1987
;
45
:
683
–90.
38
Gilchrest BA, Eller MS, Geller AC, Yaar M. The pathogenesis of melanoma induced by ultraviolet radiation.
N Engl J Med
1999
;
340
:
1341
–8.
39
Ovechkin AV, Tyagi N, Rodriguez WE, Hayden MR, Moshal KS, Tyagi SC. Role of matrix metalloproteinase-9 in endothelial apoptosis in chronic heart failure in mice.
J Appl Physiol
2005
;
99
:
2398
–405.
40
Stintzing S, Heuschmann P, Barbera L, et al. Overexpression of MMP9 and tissue factor in unstable carotid plaques associated with Chlamydia pneumoniae, inflammation, and apoptosis.
Ann Vasc Surg
2005
;
19
:
310
–9.
41
Yu F, Kamada H, Niizuma K, Endo H, Chan PH. Induction of mmp-9 expression and endothelial injury by oxidative stress after spinal cord injury.
J Neurotrauma
2008
;
25
:
184
–95.
42
Coussens LM, Raymond WW, Bergers G, et al. Inflammatory mast cells up-regulate angiogenesis during squamous epithelial carcinogenesis.
Genes Dev
1999
;
13
:
1382
–97.
43
Nishizuka I, Ichikawa Y, Ishikawa T, et al. Matrilysin stimulates DNA synthesis of cultured vascular endothelial cells and induces angiogenesis in vivo.
Cancer Lett
2001
;
173
:
175
–82.
44
Pozzi A, Moberg PE, Miles LA, Wagner S, Soloway P, Gardner HA. Elevated matrix metalloprotease and angiostatin levels in integrin α 1 knockout mice cause reduced tumor vascularization.
Proc Natl Acad Sci U S A
2000
;
97
:
2202
–7.
45
Cotignola J, Roy P, Patel A, et al. Functional polymorphisms in the promoter regions of MMP2 and MMP3 are not associated with melanoma progression.
J Negat Results Biomed
2007
;
6
:
9
.

Supplementary data