
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
Human Cancer Biology |
Authors' Affiliation: Veridex, LLC, a Johnson & Johnson Company, San Diego, California
Requests for reprints: Yixin Wang, Veridex, LLC, a Johnson & Johnson Company, 3210 Merryfield Row, San Diego, CA 92121. Phone: 858-784-3295; Fax: 858-450-2070; E-mail: ywang{at}vrxus.jnj.com.
| Abstract |
|---|
|
|
|---|
Experimental Design: Total RNA isolated from 45 primary melanoma, 18 benign skin nevi, and 7 normal skin tissue specimens were analyzed on an Affymetrix Hu133A microarray containing 22,000 probe sets.
Results: Hierarchical clustering revealed a distinct separation of the melanoma samples from the benign and normal specimens. Novel genes associated with malignant melanoma were identified. Differential gene expression of two melanoma-specific genes, PLAB and L1CAM, were tested by a one-step quantitative reverse transcription-PCR assay on primary malignant melanoma, benign nevi, and normal skin samples, as well as on malignant melanoma lymph node metastasis and melanoma-free lymph nodes. The performance of the markers was compared with conventional melanoma markers such as tyrosinase, gp100, and MART1.
Conclusion: Our study systematically identified novel melanoma-specific genes and showed the feasibility of using a combination of PLAB and L1CAM in a reverse transcription-PCR assay to differentiate clinically relevant samples containing benign or malignant melanocytes.
High-density microarrays have been applied to simultaneously monitor the expression of thousands of genes in biological samples. Studies have resulted in the identification of genes differentially expressed in benign and malignant lesions (12), as well as genes that might be of prognostic value (13). In their pioneering work, Bittner et al. (14) did gene expression profiling of malignant melanoma, using a microarray containing probes for 8,150 cDNAs. These researchers identified several genes that might be associated with aggressive tumor behavior. In the recent work of Hoek et al. (15), comparison of gene expression profiles of a few melanoma and normal melanocyte cell lines led to the identification of differentially expressed genes and pathways modulated in melanoma. In this study, we report the gene expression profiling of an extensive set of clinically relevant tissue samples. Forty-five primary malignant melanomas, 18 benign skin nevi, and 7 normal skin tissues were hybridized on an Affymetrix Hu133A microarray (Santa Clara, CA) containing 22,000 probe sets. Differentially expressed genes in malignant melanoma, as compared with benign tissue, were identified. A one-step quantitative RT-PCR assay was used to test a combination of two melanoma-specific genes, PLAB and L1CAM, in a panel of clinically relevant samples that included primary malignant melanoma, benign nevi, melanoma lymph node metastasis and melanoma-free lymph node samples.
| Materials and Methods |
|---|
|
|
|---|
RNA isolation and microarray hybridization. Qiagen RNeasy mini kit (Qiagen Inc., Valencia, CA) was used with a modified protocol to minimize the residue melanin in the RNA sample. For melanocyte-containing tissues, four replicate tissue samples derived from individual patients, each weighed
5 mg, were used and processed separately. Tissue samples were homogenized in 1.0 mL RLT buffer (Qiagen) containing 10 µL of ß-mercaptoethanol (Sigma Chemical Co., St. Louis, MO) by a mechanical homogenizer (UltraTurrex T8, IKA-Werke, Staufen, Germany). Following homogenization, samples were loaded on Qiagen RNeasy columns and followed by centrifugation. After discarding the flow-through, 700 mL of RW1 buffer was added. The column was kept at room temperature for 5 minutes and then centrifuged; this step was repeated thrice. Then we followed the standard Qiagen RNeasy mini kit protocol. To remove RNA from the silica gel membrane, a two-step elution was done. The total RNA derived from the same individual patient tissue was pooled and used for further analysis.
Standard Trizol protocol was used for RNA isolation from tissues that do not contain a significant proportion of melanocytes. Tissue was homogenized in Trizol reagent (Invitrogen, Carlsbad, CA). After centrifugation, the top liquid phase was collected and total RNA was precipitated with isopropyl alcohol at 20°C. RNA pellets were washed with 75% ethanol, resolved in water and stored at 80°C until use. The RNA quality was examined with an Agilent 2100 Bioanalyzer RNA 6000 Nano Assay (Agilent Technologies, Palo Alto, CA).
Labeled cRNA was prepared and hybridized with the high-density oligonucleotide array Hu133A Gene Chip (Affymetrix) containing a total of 22,000 probe sets according to standard manufacturer protocols. Arrays were scanned using Affymetrix protocols and scanners. For subsequent analysis, each probe set was considered as a separate gene. Expression values for each gene were calculated by using Affymetrix Gene Chip analysis software MAS 5.0. Arrays that met three quality control standards were used: "present" call was >35%, scale factor was <12 when scaled to a target intensity of 600, and background level was <150. A lower than usual percentage of "present" calls cutoff was chosen because it is difficult to isolate undegraded RNA from the skin samples (16), resulting in lower overall gene expression levels.
Data analysis. Gene expression data were filtered to include only genes that are called "present" in two or more samples. This filter was used to remove genes that did not change expression in the samples. Of the 22,000 genes presented on the array, 15,795 passed this filter and were used for hierarchical clustering. Prior to clustering, each gene expression signal was divided by the median expression in all the samples in the data set. This standardization step minimized the effect of the magnitude of gene expression and grouped together genes with similar expression patterns in the clustering analysis. Average linkage hierarchical clustering using Pearson correlation was done on both the genes and the samples using GeneSpring 6.1 software.
In order to identify differentially expressed genes, we compared the melanoma samples to the benign nevi and the normal skin samples separately. The first analysis consisted of the 45 melanoma and 7 normal skin samples; the second analysis consisted of 45 melanoma and 18 nevi samples. Significance analysis of microarrays (17), percentile analysis, and t test were used in gene selection (Supplementary Fig. S1). We selected a short list of genes with at least 10-fold overexpression in melanoma as compared with the benign specimens. The complete array data set has been submitted to the National Center for Biotechnology Information/Genbank GEO database (series entry GSE3189).
Reverse transcription-PCR validation of microarray results. Ten micrograms of total RNA from each sample were treated with DNase I and reverse-transcribed with oligo (dT) primer using Superscript II reverse transcriptase according to the manufacturer's instructions (Invitrogen). A control gene PBGD was previously tested and reported as a housekeeping gene (18). Primers and MGB-probes for me20m (gp100), L1CAM, NTRK3, and the control gene PBGD were designed using Primer Express software (Applied Biosystems, Foster City, CA). The PLAB (MIC1) gene probe was FAM-TAMRA-based because sequences were inadequate to design MGB based probes. Primer/probe sequences were as follows: me20m (gp100) forward, 5'-TGTGTCTCTGGCTGATACCAACA-3'; me20m (gp100) reverse, 5'-TTCTTGACCAGGCATGATAAGCT-3'; me20m (gp100) probe, 5'-(6-FAM) CTGGCAGTGGTCAGC-3'; L1CAM forward, 5'-GCTGGGACTGGGAACAGAACT-3'; L1CAM reverse, 5'-GGAGCAGAGATGGCAAAGAAA-3'; L1CAM probe, 5'-(6-FAM) TCCCCACCATCTGCTGT-3'; NTRK3 forward, 5'-GCCCCGGCACCCTTTA-3'; NTRK3 reverse, 5'-AACCCTGCCAGTGGTGGAT-3'; NTRK3 probe, 5'-(6-FAM) CAGATGGGTGTTTTC-3'; PLAB (MIC1) forward, 5'-GGCAGAATCTTCGTCCGCA-3'; PLAB (MIC1) reverse, 5'-GGACAGTGGTCCCCGTTG-3'; PLAB (MIC1) probe, 5'-(6-FAM) CCCAGCTGGAGTTGCACTTGCGGCC(TAMRA)-3'; PBGD forward, 5'-CTGCTTCGCTGCATCGCTGAAA-3'; PBGD reverse, 5'-CAGACTCCTCCAGTCAGGTACA-3'; PBGD probe, 5'-(6-FAM) CCTGAGGCACCTGGAAGGAGGCTGCAGTGT(TAMRA)-3'. All primers and probes were tested for optimal amplification efficiency above 90%. The standard curve was composed of six 10-fold dilutions of target gene PCR product with copy numbers ranging from 10 to 106. RT-PCR amplification was carried out in a 20 µL reaction mix containing 50 ng template cDNA, 2x TaqMan universal PCR master mix (12.5 µL; Applied Biosystems), 500 nmol/L forward and reverse primers, and 250 nmol/L probe. Reactions were run on an ABI PRISM 7900HT Sequence Detection System (Applied Biosystems). The cycling conditions were: 2 minutes of AmpErase UNG activation at 50°C, 10 minutes of polymerase activation at 95°C and 50 cycles at 95°C for 15 seconds and annealing temperature (60°C) for 60 seconds. In each assay, a standard curve and a no-template control along with template cDNA were included in duplicate for both the gene of interest and the control gene. The relative quantity of each target gene was represented as
Ct, which is equal to Ct of the target gene subtracted by Ct of the control gene.
One-step quantitative reverse transcription-PCR assay on primary and lymph node samples. Evaluation of expression of selected genes was carried out with one-step RT-PCR with RNA from primary melanoma, benign nevi, normal skin, melanoma lymph node metastasis, and melanoma-free lymph nodes. ß-Actin was used as a housekeeping gene to control for the input quantity and quality of RNA in the reactions. DNase treatment was not used. Instead, primers or probes were designed to span an intron so they would not report on genomic DNA. Eight nanograms of total RNA were used for the RT-PCR. The total RNA was reverse-transcribed using 40x Multiscribe and RNase inhibitor mix contained in the TaqMan one-step PCR Master Mix reagents kit (Applied Biosystems). The cDNA was then subjected to the 2x Master Mix without UNG, and PCR amplification was carried out on the ABI 7900HT sequence detection system (Applied Biosystems) in the 384-well block format using a 10 µL reaction size. The primer and probe concentrations were 4 and 2.5 µmol/L, respectively. The reaction mixture was incubated at 48°C for 30 minutes for the reverse transcription, followed by a Amplitaq activation step of 95°C for 10 minutes and finally 40 cycles of 95°C for 15 seconds denaturing and 60°C for 1 minute annealing and extension. On each plate, a standard curve is generated from 8 pg to 80 ng, and when the R2 value was >0.99, the cycle threshold (Ct) values were accepted. Sequences used in the reactions were as follows, each written in the 5' to 3' direction: L1CAM forward, CCACAGATGACATCAGCCTCAA; L1CAM reverse, GGTCACACCCAGCTCTTCCTT; L1CAM probe, TGGCAAGCCCGAAGTGCAGTTCC; tyrosinase forward, CTTTAGAAATACACTGGAAGGATTTGCTA; tyrosinase reverse, CATTGTGCATGCTGCTTTGA; tyrosinase probe, TCCACTTACTGGGATAGCGGATGCCTC; MART1 forward, ACTTCATCTATGGTTACCCCAAGAA; MART1 reverse, TCCCAGCGGCCTCTTCA; MART1 probe, CACGGCCACTCTTACACCACGGC; me20m forward, CTTAAGGCTGGTGAAGAGACAAGTC; me20m reverse, CAGGATCTCGGCACTTTCAATAC; me20m probe, TCGATATGGTTCCTTTTCCGTCACCCTG; PLAB forward, ATTCGAACACCGACCTCGTC; PLAB reverse, CGCAGGTGCAGGTGGC; PLAB probe, GATACTCACGCCAGAAGTGCGGCT. All primers and probes were optimized towards the same amplification efficiency.
For each sample
Ct = Ct (target gene) Ct ß-actin was calculated.
Ct normalization has been widely used in clinical RT-PCR assays and was chosen as a straightforward method (19). For each gene marker, t test was done on
Ct between the melanoma and nonmelanoma samples, including both primary and lymph node samples. We then used the
Ct to construct two scores for each patient. One score was derived from a combination of two melanoma-specific genes, PLAB and L1CAM; and the other score was derived from a combination of three conventional melanoma markers, tyrosinase, me20m, and MART1. The score was defined as the weighted sum of
Ct values of the tested genes with the corresponding t statistics as the weight. The two scores were normalized to have the same mean in order to compare them on the same scale.
| Results |
|---|
|
|
|---|
|
Unsupervised hierarchical clustering results revealed a distinct separation of the melanoma, benign nevi, and normal skin samples (Fig. 1). We observed four clusters, including two clusters consisting of the majority of the melanoma samples (43 out of 45), the third cluster contained all 7 normal skin, 3 benign nevi, and 2 melanoma specimens and the fourth cluster, that included 14 of the 18 benign nevi samples. The source of the samples didn't affect clustering. Specimens originating from different sources were clustered together according to the sample type (melanoma, benign, or normal). To further test the stability of the clustering patterns, we used an alternative cutoff on gene filtering prior to the cluster analysis. Specifically, we retained genes that have at least 10% "present" calls in each of the melanoma, benign nevi, and skin samples. With this cutoff, we obtained 15,306 genes and repeated hierarchical clustering. The cluster pattern on the patient samples was the same as the one from the 15,795 from the two "present" calls, confirming clustering stability.
|
Distinct gene clusters were found in association with melanoma. This can be characterized by up-regulated (Fig. 1A, B, and C) and down-regulated (Fig. 1E) genes in the melanoma samples. At the same time, melanoma and benign nevi samples showed high expression of known melanocyte markers, such as MART-1 (Fig. 2D), confirming the comparable content of melanocytes in these samples and the inability of melanocyte-specific markers to differentiate them. Our data indicates that melanoma, benign nevi, and normal skin samples have distinct gene expression profiles and can be separated on a molecular basis.
|
|
One-step quantitative reverse transcription-PCR assay for malignant melanoma. We examined a combination of two genes highly overexpressed in melanoma, PLAB and L1CAM, in a variety of clinical tissue samples containing malignant melanocytes (primary melanoma and melanoma lymph node metastasis), benign melanocytes (benign skin nevi), and normal samples (normal skin and melanoma-free lymph node) by RT-PCR. The primary tissues were the same as those used for the microarray study whereas all the lymph node specimens were derived from independent patients. Conventional melanoma markers, such as tyrosinase, gp100, and MART1 were also tested on the same samples as the controls because they are the most commonly used markers for the melanoma molecular assays in current clinical studies (2527). Calculated scores were presented for PLAB and L1CAM (Fig. 3A) and for tyrosinase, me20m, and MART1 (Fig. 3B). The results showed a significant difference in the expression of PLAB and L1CAM between malignant melanoma samples (primary and lymph node metastasis) and benign nevi and normal lymph node. In contrast, three conventional markers showed similar expression levels in benign and melanoma samples. To further show the ability of these gene markers to separate benign and malignant tissues, we tested two cutoffs; the first was set-up as the highest score in primary normal samples, and the second as the highest score in benign nevi samples. For each cutoff, we estimated the sensitivity and specificity of the assay in the lymph node samples. With the cutoff determined on the normal samples, the new markers and the conventional markers gave sensitivities of 90% and 83%, respectively. Using of the cutoff determined on the benign samples, the sensitivities for the new and conventional markers were 88% and 42%. The results indicated that the new markers potentially have much better abilities to differentiate tissues containing benign and malignant melanocytes.
|
| Discussion |
|---|
|
|
|---|
We have compared our results to the recent study of Haqq et al. (30). In their work, a cDNA microarray containing 20,862 probes was used to profile benign nevi, primary melanoma, and metastatic melanoma specimens. The sample set included metastatic and primary melanoma and benign nevi. Similar clustering results that separated the benign nevi and the primary malignant melanoma tissues were found in their study. Common genes were reported in both studies which could discriminate melanoma from benign nevi including kinesin-like 5 (KNSL5), prostate differentiation factor (PLAB), CITED1, osteopontin (SPP1), cathepsin B (CSTB), cadherin 3 (CDH3), presenilin 2 (PSEN2).
The results of our one-step RT-PCR assay showed that novel melanoma-specific genes, PLAB and L1CAM, were expressed not only in primary melanoma tissues but also in melanoma lymph node metastasis. Moreover, their ability to differentiate malignant melanoma from benign nevi made them better candidates than the conventional markers for the molecular test of melanoma diagnostics. With further validation in clinical studies, these genes could be developed as specific markers for a molecular staging assay to detect melanoma micrometastasis during a sentinel lymph node biopsy procedure. Another potential application of the genes is for the diagnosis of melanocyte lesions with uncertain pathologic features. We are evaluating the markers for use in sentinel lymph node biopsies determining the presence of melanoma disease and the prognosis of tumors.
| Acknowledgments |
|---|
| Footnotes |
|---|
Note: Supplementary data for this article are available at Clinical Cancer Research (http://clincancerres.aacrjournals.org/).
Received 3/28/05; revised 5/31/05; accepted 7/14/05.
| References |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
J. W.J. van Kilsdonk, R. H. Wilting, M. Bergers, G. N.P. van Muijen, J. Schalkwijk, L. C.L.T. van Kempen, and G. W.M. Swart Attenuation of Melanoma Invasion by a Secreted Variant of Activated Leukocyte Cell Adhesion Molecule Cancer Res., May 15, 2008; 68(10): 3671 - 3679. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. R. Moore, D. L. Persons, J. A. Sosman, D. Bobadilla, V. Bedell, D. D. Smith, S. R. Wolman, R. J. Tuthill, J. Moon, V. K. Sondak, et al. Detection of Copy Number Alterations in Metastatic Melanoma by a DNA Fluorescence In situ Hybridization Probe Panel and Array Comparative Genomic Hybridization: A Southwest Oncology Group Study (S9431) Clin. Cancer Res., May 15, 2008; 14(10): 2927 - 2935. [Abstract] [Full Text] [PDF] |
||||
![]() |
C. Rolny, L. Capparuccia, A. Casazza, M. Mazzone, A. Vallario, A. Cignetti, E. Medico, P. Carmeliet, P. M. Comoglio, and L. Tamagnone The tumor suppressor semaphorin 3B triggers a prometastatic program mediated by interleukin 8 and the tumor microenvironment J. Exp. Med., May 12, 2008; 205(5): 1155 - 1171. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. Bazzaro, A. Santillan, Z. Lin, T. Tang, M. K. Lee, R. E. Bristow, I.-M. Shih, and R. B.S. Roden Myosin II Co-Chaperone General Cell UNC-45 Overexpression Is Associated with Ovarian Cancer, Rapid Proliferation, and Motility Am. J. Pathol., November 1, 2007; 171(5): 1640 - 1649. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. B. Axelsen, J. Lotem, L. Sachs, and E. Domany Genes overexpressed in different human solid cancers exhibit different tissue-specific expression profiles PNAS, August 7, 2007; 104(32): 13122 - 13127. [Abstract] [Full Text] [PDF] |
||||
![]() |
H. Jin, X. Wang, J. Ying, A. H. Y. Wong, Y. Cui, G. Srivastava, Z.-Y. Shen, E.-M. Li, Q. Zhang, J. Jin, et al. Epigenetic silencing of a Ca2+-regulated Ras GTPase-activating protein RASAL defines a new mechanism of Ras activation in human cancers PNAS, July 24, 2007; 104(30): 12353 - 12358. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. R. Alonso, L. Tracey, P. Ortiz, B. Perez-Gomez, J. Palacios, M. Pollan, J. Linares, S. Serrano, A. I. Saez-Castillo, L. Sanchez, et al. A High-Throughput Study in Melanoma Identifies Epithelial-Mesenchymal Transition as a Major Determinant of Metastasis Cancer Res., April 1, 2007; 67(7): 3450 - 3460. [Abstract] [Full Text] [PDF] |
||||
![]() |
L. Packer, S. Pavey, A. Parker, M. Stark, P. Johansson, B. Clarke, P. Pollock, M. Ringner, and N. Hayward Osteopontin is a downstream effector of the PI3-kinase pathway in melanomas that is inversely correlated with functional PTEN Carcinogenesis, September 1, 2006; 27(9): 1778 - 1786. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| Cancer Research | Clinical Cancer Research |
| Cancer Epidemiology Biomarkers & Prevention | Molecular Cancer Therapeutics |
| Molecular Cancer Research | Cancer Prevention Research |
| Cancer Prevention Journals Portal | Cancer Reviews Online |
| Annual Meeting Education Book | Meeting Abstracts Online |