High ALDH1 expression correlates with better prognosis in tumorigenic malignant melanoma.

Aldehyde dehydrogenase 1 (ALDH1) has been proposed as biomarker of stem cells for certain human cancers. ALDH1 expression has been correlated with poor patient outcomes in a variety of malignancies but better patient outcomes in others, and its prognostic significance in malignant melanoma is unclear. Thus, 68 melanoma patients with comprehensive clinical and pathologic follow-up data were used to construct a tissue microarray. A modified histological score (H-score) with a maximum score of 300 was used to quantify immunohistochemical staining for ALDH1. Survival time was defined as the time between diagnosis and melanoma-specific death. Using univariate logistic regression, a low (<80 H-score) ALDH1 score showed 3.7-fold increase in risk for melanoma-specific death within 10 years when compared with high (>80) ALDH1 levels (P=0.017). Odds of MSD were lower by a factor of ~0.9 for each 10-point increase in H-Score. Median survival time was 44.1 months and 180.9 months for patients with low and high ALDH1 expression, respectively. Using multivariate analysis, ALDH1 H-score was found to be an independent prognostic factor. These findings suggest that ALDH1 expression in malignant melanoma has a favorable effect on patient survival. Further study is needed elucidate the function of this enzymatic protein in melanoma progression.

Introduction

Reported incidence of malignant melanoma has been increasing in the U.S. since the 1930s.[1,2] Although it accounts for less than one percent of skin cancer it is responsible for approximately 75% of skin-cancer-related deaths.[3] Due to its aggressive nature and its ability to acquire resistance to many standard therapies there is significant interest in understanding the targetable molecular pathways involved in its tumorigenesis in hopes of discovering new and effective treatment alternatives.[4] It has been hypothesized that tumor initiation and therapy resistance may be associated with cells exhibiting so-called cancer-stem-cell[5] properties. Levels of activity of the hypothesized cancer-stem-cell marker Aldehyde dehydrogenase (ALDH), and in particular its isoform 1, have been shown to be consistent with cancer-stem-cell activity in a variety of cancers.[6-14] Unfortunately studies of melanoma have yielded conflicting results.[15-19]

The level of activity of ALDH1, which belongs to a family of ubiquitous detoxifying enzymes that function to convert aldehydes to their respective carboxylic acids,[20,21] is associated with poor clinical outcomes in breast cancer[6] and leukemia,[13] but correlates with better clinical outcomes in ovarian cancer[22] and non-small-cell lung cancer.[23] Recently, a study testing chemotherapeutics in melanoma xenograft mice models has shown promising results for ALDH-targeted therapies demonstrating the need for increased understanding of this pathway and its role in disease progression.[24] However, to date no investigations have directly assessed the prognostic significance of ALDH1 in melanoma. In this paper, we investigate the relationship between clinical outcomes and ALDH1 activity in melanoma.

Materials and Methods

Study patients and design

The 68 patients eligible for this retrospective study of a prospective cohort had vertical growth phase tumorigenic primary melanomas, no regional nodal procedure or apparent metastases at the time of definitive treatment between 1972 and 1991 at the University of Pennsylvania's Pigmented Lesion Clinic, and at least ten years of follow up. They had at least 10 years of protocol-driven, prospective follow-up and had paraffin blocks available for tissue microarray. The protocol was approved by the Institutional Review Board of the University of Pennsylvania.

Tissue microarray construction and immunohistochemical staining

Tissues from 68 patients were included in the tissue microarray. These patients had vertical growth phase primary melanomas and no regional nodal procedure or apparent metastases at the time of definitive treatment between 1972 and 1991 at the University of Pennsylvania’s Pigmented Lesion Clinic. The tissue microarray was composed of 2 mm cores. When we designed the tissue microarray, we considered the possibility of heterogeneity within melanomas. Sixty three percent of cases had three cores from different areas of primary melanomas and thirty four percent of the cases had four cores. Two cases had one or two cores due to smaller primary tumor volumes. Immunohistochemical staining was performed for ALDH1 antibody (BD biosciences) using a 1:800 dilution with citrate retrieval. ALDH1 expression was analyzed by two independent pathologists (RMA, XX) using a modified histological score (H-score) based both on the percentage of positively stained cells and on the intensity of staining, with a maximum score of 300. For statistical analysis, the scores for each core reported by the two readers were averaged. The tissue microarray was constructed at the Penn tissue microarray facility in the Department of Pathology and Laboratory Medicine.

Pathology definitions

All original hematoxylin and eosin (H&E)-stained slides were read at the time of definitive treatment by pathologists associated with the Pigmented Lesion Clinic of the University of Pennsylvania (chiefly David E. Elder and Wallace H. Clark, Jr.). Attributes recorded at the time of diagnosis include thickness and anatomical level; dermal mitotic rate, expressed in terms of mitoses per square millimeter; vertical growth phase tumor-infiltrating lymphocytes, classified as either brisk, nonbrisk, or absent; regression, if present in the radial growth phase that often accompanies the vertical growth phase; microscopic satellites; ulceration; and vascular (blood or lymphatic) invasion. For statistical analysis, some characteristics were converted to categorical variables: thickness (1.01 to 2.00 mm, 2.01 to 4.00 mm or 4.01 mm and more); mitotic rate (<1 and ≥1 mitoses/mm2) and Clark’s level (III, IV/V).

Clinical definitions

Clinical characteristics available for all patients included age, sex, and the anatomic site of the primary lesion. For statistical analysis, binary variables were created for age at diagnosis (≥ 60 years, < 60 years) and anatomic site (extremity, axial). Survival time was defined as the time between definitive treatment and melanoma-specific death. The 10-year melanoma-specific death rate was defined as the proportion of patients who died of melanoma within 10 years of definitive treatment.

Statistical Methods

Univariate logistic regression analyses were used to estimate unadjusted odds ratios for 10-year melanoma-specific death and H-score as well as each of the other 10 available factors. The three factors that were significantly associated with 10-year melanoma-specific death in univariate models were included in a multivariate model. One factor that that did not make a significant contribution to the multivariate model was removed based on a likelihood ratio test to obtain a parsimonious, reduced multivariate model. Based on the number of events in this sample, 2–3 variables in the regression model are recommended to avoid overfitting (See Harrell, 2015, pp. 72–74).[26] X-tile software[25] was used to compute optimal H-score subgroup thresholds (2 thresholds, and thus “low,” “middle,” and “high” groups). These thresholds were determined by computing the two Kaplan-Meier survival curves defined by all possible H-score cut points, and the associated log-rank statistics to identify the cut-point where the two curves were most different based on the log-rank statistics. Exact chi-square tests (SAS Proc FREQ) were used to compare distributions for pathological and clinical factors among the H-Score subgroups. All statistical analyses were performed using SAS/STAT software, Version 9.4 of the SAS System for PCs.

Results

The ALDH1 immunohistochemical staining pattern was primarily cytoplasmic. We noted heterogeneity of ALDH1 staining in individual tumor cells. However, the variability within a tumor is relatively uniform which is reflected by similar H-scores of different cores from the same case. We stained some whole sections of melanomas with ALDH1 and confirmed that the variability in staining was relatively even within a melanoma. There was no specific distribution of positive cells in relationship to tumor necrosis or perivascular localization.

Table 1 presents the clinical and histopathologic characteristics for the 68 study patients. The median age at diagnosis was approximately 59 years (range, 21–82 years). The median thickness was 3.25 mm (range, 1.04–8.86mm). The majority of patients were Clark’s level IV/V (56%), 46% exhibited ulceration, and satellites were present in 72% of cases. Median mitotic rate was 6.3 mitoses/mm2 (range, 0–38.9 mitoses/mm2). The mean H-Score was computed for each patient using either 4 cores (34%), 3 cores (63%) or 1 or 2 cores (2.7%). The median H-score was 27 (range, 0–241). Within 10 years of definitive excision of their primary lesions, 42 patients died of causes related to their disease, resulting in a 10-year melanoma specific death rate of 62%.

Table 1

Patient and Tumor Characteristics, Overall and by H-Score Subgroup

	All Patients	H-Score Subgroups
	n=68Percent	Low(0 – 2)n=10Percent	Intermediate(3 – 120)n=43Percent	High(121 – 300)n=15Percent	p-value*
Age at Diagnosis					0.430
< 60 years	55.9	40	55.8	66.7
≥ 60 years	44.1	60	44.2	33.3
Sex					0.564
Male	67.6	80	62.8	73.3
Female	32.4	20	37.2	26.7
Anatomic Site					0.863
Axial	72.1	80	69.8	73.3
Extremity	27.9	20	30.2	26.7
Thickness					0.770
1.01–2.00mm	14.7	20	14.0	13.3
2.01–3.00mm	27.9	10	30.2	33.3
3.01–9.00mm	57.4	70	55.8	53.3
Regression					0.055^
Present	44.1	50	44.1	84.6
Absent	36.8	50	55.9	15.4
Unknown	19.1	---	---	---
Mitotic Rate**					0.645^
< 1.0	8.8	10	4.8	13.3
1.0–40.0	89.7	90	95.2	86.7
Unknown	1.5	---	---	---
TIL					0.912
Present	77.9	80	79.1	73.3
Absent	22.1	20	20.9	26.7
Ulceration					0.349^
Present	45.6	70	45.0	42.9
Absent	48.5	30	55.0	57.1
Unknown	5.9	---	---	---
Clark’s Level					0.016
III	44.1	20	39.5	73.3
IV & V	55.9	80	60.5	26.6
Satellites					0.331
Present	27.9	40	30.2	13.3
Absent	72.1	60	69.8	86.7

Abreviation: TIL, tumor infiltrating lymphocytes.

Exact chi-square to test for equality of proportions across the three H-score subgroups.

Mitoses/mm2

Unknown values excluded from the analysis

Three factors were significantly associated with 10-year melanoma-specific death in the univariate logistic regression analyses (Table 2, Univariate Models). ALDH1 H-score was inversely associated with the likelihood of 10-year melanoma-specific death (odds ratio=0.91, p=0.020), indicating that for each 10-point increase in ALDH1 H-score, the odds of 10-year melanoma-specific death decreases by a factor of 0.91. Both microscopic satellites (odds ratio=4.72, p=0.025) and Clark level IV/V (odds ratio=3.20, p=0.025) were associated with a higher odds of 10-year melanoma-specific death compared to patients without microscopic satellites and for patients with Clark level III melanomas, respectively.

Table 2

Logistic Regression Analysis for 10-year Melanoma-specific Death (n=68)

	Univariate Models		Multivariate Modela		Reduced Model
	UnadjustedOR	p-value	AdjustedOR	p-value	AdjustedOR	p-value
Modified H-score*	0.91	0.020	0.93	0.089	0.92	0.034
Satellites
Present	4.72	0.025	3.78	0.070	4.35	0.040
Absent	1.00	---	1.00	---	1.00	---
Clark’s Level
III	1.00	0.025	1.00	0.283	---	---
IV & V	3.2	---	1.86	---	---	---
Thickness (mm)	1.34	0.060	---	---	---	---
Mitotic Rate**	1.05	0.273^	---	---	---	---
Age at Diagnosis (years)	1.01	0.425	---	---	---	---
Sex
Male	0.89	0.826	---	---	---	---
Female	1.00	---	---	---	---	---
Anatomic Site
Axial	2.29	0.133	---	---	---	---
Extremity	1.00	---	---	---	---	---
Regression
Present	0.76	0.540^	---	---	---	---
Absent	1.00	---	---	---	---	---
TIL
Present	0.76	0.659	---	---	---	---
Absent	1.00	---	---	---	---	---
Ulceration
Present	1.34	0.570^	---	---	---	---
Absent/Unknown	1.00	---	---	---	---	---

Abbreviation: OR, odds ratio

Modified H-score (0 to 30) is equal to the H-score divided by 10.

Mitoses per high powered field

Unknown data excluded from the analysis

Multivariate model included the 3 variables significant in the univariate analyses

The multivariate logistic regression model included the three factors: ALDH1 H-score, the presence of satellites, and Clark’s level, (see Table 2 Multivariate model). The reduced multivariate model included ALDH1 H-score (odds ratio=0.92, p =0.034) and microscopic satellites (odds ratio=4.35, p=0.040) as independent prognostic factors for 10-year melanoma-specific death (Table 2, Reduced model). The interpretations and magnitudes of these coefficients are similar to those from the univariate models.

To further investigate the relationship between survival and ALDH1 expression, two optimal cut-points for H-score were determined to examine characteristics of those with low, intermediate and high risk of 10-year melanoma-specific death. The two optimal cut points were estimated to be 2.83 and 120 and three subgroups were defined: low ALDH1 (H-score between 0 and 2), intermediate ALDH1 (H-score between 3 and 119) and high ALDH1 (H-score between 120 and 300). The “high” group had significantly more lesions having Clark’s level III (73% in the high group versus 20% and 40% in low and medium groups, p=0.016) and had marginally significantly more lesions with radial growth phase regression (73% in the high group versus 40% and 35% in low and medium groups, p=0.055).

Discussion

In this study of 68 patients with tumorigenic malignant melanoma, comprehensive staging and long-term clinical follow-up, high levels of ALDH1 activity as interpreted by immunohistochemical staining was associated with better patient outcomes. Odds of 10-year melanoma-specific death was lower by a factor of approximately 0.9 for each 10-point increase in H-Score. Using multivariate analysis, ALDH1 H-score was also found to be an independent prognostic marker (p-value = 0.034). Overall, the data support that higher ALDH1 activity in this cohort correlated with better prognosis.

As noted above, studies have linked ALDH1 levels with both positive and negative prognostic significance depending on the type of cancer studied. Our finding of higher levels of ALDH1’s association with better prognosis in melanoma mirrors similar findings in non-small cell lung cancer[23] but contrasts with data on breast cancer[6,11] and leukemia,[13] where higher ALDH1 levels correlated with poor prognosis. These apparent discrepancies between ALDH1 expression in different tumors and clinical outcomes could reflect the notion that ALDH1, a ubiquitous metabolic enzyme, may serve different functions in each organ system, and thus may play a specific role in each organ’s tumor biology. Its positive correlation with prognosis suggests that ALDH1 expression is unlikely to serve as a cancer-stem-cell marker for melanoma in the way it may in breast cancer, where ALDH1 expression is believed to correlate with cancer-stem-cell status. Moreover, it may be a marker of differentiation rather than undifferentiation.

ALDH1 is important for tumor cell proliferation, survival, and resistance to chemotherapeutic agents. Small molecules that inhibit ALDH1 functions have been developed, such as diethylaminobenzaldehyde and disulfiram. Therefore, targeted therapy to ALDH1 may be tested in preclinical models and future clinical trials. However, as noted above, the effects of ALDH1 expression on tumor growth and metastasis appear to be tumor-type specific. To better target ALDH1 for therapy, additional studies are needed to understand the mechanisms underlying functional role of ALDH1 in different tumor types.