|Home | About | Journals | Submit | Contact Us | Français|
A giant congenital melanocytic nevus (GCMN) is a malformation of the pigment cells. It is a distress to the patients for two reasons: one is disfigurement, and the other is the possibility of malignant changes. However, the underlying mechanisms of the development of GCMN and melanotumorigenesis in GCMN are unknown. Hence, the aim of this study was to identify the proteomic alterations and associated functional pathways in GCMN.
Proteomic differences between GCMN (n=3) and normal skin samples (n=3) were analyzed by one-dimensional-liquid chromatography-tandem mass spectrometry Relative levels of the selected proteins were validated using western blot analysis. The biological processes associated with the abundance modified proteins were analyzed using bioinformatic tools. Among the 46 abundance modified proteins, expression of 4 proteins was significantly downregulated and expression of 42 proteins was significantly upregulated in GCMN compared to normal skin samples (p<0.05). More importantly, 31% of the upregulated proteins were implicated in various cancers, with five proteins being specifically related with melanoma. The abundance modified proteins in GCMN were involved in the biological processes of neurotrophin signaling, melanosome, and downregulated of MTA-3 in ER-negative breast tumors. In particular, an increase in the expression of the 14-3-3 protein family members appeared to be associated with key cellular biological functions in GCMN. Western blot analysis confirmed the upregulation of 14-3-3epsilon, 14-3-3 tau, and prohibitin in GCMN.
These findings suggest that GCMN exhibits potential proteomic alterations, which may play a role in melanotumorigenesis, and the significant alteration of 14-3-3 family proteins could be a key regulator of the biological pathway remodeling in GCMN.
Congenital melanocytic nevi (CMN) are pigment cell malformations that are visible at birth, or are nevi showing congenital features that become clinically obvious shortly after birth . CMN is caused by abnormal melanocyte differentiation, migration, and deposition in the dermis during the early stages of embryogenesis [2,3]. It is a distress to patients for two reasons: one is disfigurement, and the other is the increased risk of developing malignant melanoma [3,4], especially in individuals with giant congenital melanocytic nevi (GCMN; over 20 cm in diameter) . Several genomic and proteomic studies have been performed to elucidate the mechanism of melanotumorigenesis arising from CMN. Gene-based analyses have revealed that the oncogenic BRAF and NRAS mutations are frequently seen in CMN. Additionally, increased Bcl-2 expression in CMN has been suggested to suppress apoptosis, which otherwise plays an important role in the maintenance of nevocytes . In spite of these findings, the major biological processes and pathways of melanotumorigenesis remain unclear. Therefore, the complete proteomic characterization of GCMN and related biological pathways through comparative protein profiling is essential to understand the underlying mechanism of the origin of GCMN and the subsequent process of melanotumorigenesis.
In this study, the proteomic alteration and systemic properties of GCMN were assessed, with the aim of gaining an insight into the functional association between GCMN and melanotumorigenesis. We used label-free liquid chromatography-mass spectrometry (LC-MS) and established bioinformatic tools to identify the proteins that may play a key role in the malignant transformation of GCMN. We found that the 46 proteins were changed in protein abundance levels between normal skin and GCMN samples, and these proteins belonged to tightly organized functional clusters. Moreover, we found that five of the identified proteins were implicated in melanoma. The results of this study will improve our understanding in the biological identification of GCMN and the possible mechanisms that give rise to GCMN-associated melanotumorigenesis.
The clinical characteristics of patients who provided GCMN and normal skin samples are listed in Table Table1.1. All GCMN were over 20 cm in diameter. A representative GCMN lesion in a patient is shown in Figure Figure1A.1A. To minimize environmental bias between individuals, we collected three paired normal and GCMN skin samples from three GCMN donors for one-dimensional-liquid chromatography-tandem mass spectrometry (1D-LC-MS/MS) analysis (Table (Table1).1). Seven unpaired samples from normal skin and GCMN were used for western blot analysis to validate the proteomic results. The average age of GCMN patients who donated samples for western blot analysis was lower (5.57±1.90) than that of the normal skin donors (22.71±11.01, p=0.0058, n=7).
The same amount (20 μg) of soluble proteins obtained from GCMN and normal skin samples were analyzed on sodium dodecyl sulfate polyacrylamide gel electrophoresis (SDS-PAGE; Figure Figure1B).1B). The Coomassie blue-stained gels were sliced into 15 uniform slices. The total ion current chromatogram of each slice was acquired (Additional file 1: Figure S1), which led to the identification of 438 non-redundant proteins (Additional file 2: Table S2). Among these, 343 proteins were commonly identified in both the groups, 62 proteins were detected only in GCMN samples, and 33 proteins were detected only in normal skin samples (Figure (Figure1C).1C). Forty-six abundance modified proteins were identified whose protein abundance level was significantly (Student’s t-test with Bonferroni correction, p<0.05) different between GCMN and normal skin (Figure (Figure1D).1D). The large majority of these abundance modified proteins were upregulated (42 of 46, 91.3%), and very few (4 proteins, 8.7%) were downregulated in GCMN (Figure (Figure2,2, Table Table2).2). The most upregulated protein in GCMN was cathepsin D (CTSD) [GenBank Gene ID: 1509], which was expressed at a 6.8-fold higher level than in normal skin. In addition, protein disulfide-isomerase A3 precursor (PDIA3, 5.4-fold) [GenBank Gene ID: 2923], prohibitin (PHB, 5.3-fold) [GenBank Gene ID: 5245], heat-shock protein 70 (5.2-fold) [GenBank Gene ID: 3312], D-3-phosphoglycerate dehydrogenase (5-fold) [GenBank Gene ID: 26227], and adenosine-5'-triphosphate synthase subunit beta (5-fold) [GenBank Gene ID: 506] were highly (>4-fold) upregulated in GCMN. The 4 proteins downregulated in GCMN were cytokeratin-1 (−1.4-fold) [GenBank Gene ID: 3848], filaggrin 2 (−4.7-fold) [GenBank Gene ID: 388698], hornerin (−6.0-fold) [GenBank Gene ID: 388697], and alcohol dehydrogenase 1B (−6.2-fold) [GenBank Gene ID: 125]. Interestingly, we found that 33% (14 of total 42) of the significantly upregulated proteins in GCMN samples were implicated in various kinds of cancers (Table (Table2).2). Five of the 14 proteins, namely CTSD [GenBank Gene ID: 1509][9,10], PHB [GenBank Gene ID: 5245][11,12], chondroitin sulfate proteoglycan 4 [GenBank Gene ID: 1464] , phosphatidylethanolamine-binding protein 1 [GenBank Gene ID: 5037] , and ribosomal protein SA [GenBank Gene ID: 3921]  were related with melanoma.
Proteomic profiling of certain diseases has been increasingly used to acquire data on a large number of proteins to identify potential biomarkers and to gain insight into the underlying mechanisms of a variety of diseases. However, the massive amount of generated information often only increases the complexity of these analyses. To address this issue, we developed a systemic approach to analyze the proteomic data generated in the present study. To gain an insight into the systemic properties of the GCMN proteome, the 46 abundance modified proteins were categorized based on molecular function, biological process, and Clusters of Orthologous Groups of proteins (COG) analysis. The proteins were classified into the following classes according to their molecular functions: chaperones (16% of total), oxidoreductases (8%), and cytoskeletal elements (7%) (Figure (Figure3A).3A). Proteins associated with the biological processes of signal transduction (15%), cell structure and motility (10%), protein metabolism and modification (10%), and cell cycle (10%) were markedly altered in their expression. The COG analysis showed that the altered proteins in GCMN samples were involved in posttranslational modification, turnover, chaperones (30%), and the cytoskeleton (12%; Figure Figure3C),3C), which was consistent with the molecular function and biological process classification.
Further network analysis was carried out using the software tool, ClueGO, to identify enriched functional groups in the GCMN samples and obtain detailed information on each one. Gene Ontology (GO)_biological process, GO_cellular component, GO_molecular function, KEGG pathway, and Reactome_biocarta were selected as ontology sources. The specific functional groups identified as significantly enriched in GCMN were melanosome (GO_cellular component), neurotrophin signaling pathway (KEGG pathway), downregulated of MTA-3 in ER-negative breast tumors (Biocarta), cell cycle (KEGG pathway), phospholipase inhibitor activity (GO_molecular function), and glycolysis/gluconeogenesis (KEGG pathway; Table Table3,3, Additional file 3: Figure S2).
Because proteins rarely act alone, but rather interact as a group, thus comprising a functional cluster, a protein-protein interacting network was constructed using STRING 8.0. This network consisted of 45 of the 46 abundance modified proteins, with 52 interactions between them. In the GCMN network, the upregulated proteins were tightly linked with each other, forming a large cluster. The biological meaning of the cluster was determined by overlaying the COG of each protein on the network, which revealed shared COGs between close neighboring proteins (Figure (Figure4).4). Moreover, the GCMN network contained a highly specific protein cluster, identified as the 14-3-3 protein family, whose members were expressed at significantly higher levels than in normal skin as follows: 14-3-3 sigma (or SFN, 1.7-fold) [GenBank Gene ID: 2810], 14-3-3 protein beta/alpha (YWHAB, 2.2-fold) [GenBank Gene ID: 7529], 14-3-3 protein epsilon (YWHAE, 4.2-fold) [GenBank Gene ID: 7531], 14-3-3 protein gamma (YWHAG, 2.1-fold) [GenBank Gene ID: 7532], 14-3-3 protein theta (YWHAQ, 3.3-fold) [GenBank Gene ID: 10971], and 14-3-3 protein zeta/delta (YWHAZ, 2.6-fold) [GenBank Gene ID: 7534].
To investigate the inter-relationship between the proteomic characteristics of GCMN and melanoma, we compared our GCMN proteome network with a recently reported proteome of metastatic melanoma cell line that contained 110 non-redundant proteins . We selected 63 of the total melanoma proteins that had been included in the two high-score biological networks and analyzed protein-protein interactions between the GCMN and melanoma proteomes (Figure (Figure5A).5A). Five proteins including PHB, CTSD, elongation factor 1-alpha 1 (EEF1A1), Annexin A2 (ANXA2), and heat-shock protein beta-1 (HSPB1) were similarly upregulated in both proteome sets. Subsequently, we found that 217 interactions existed between the GCMN and melanoma proteomes. To estimate the biological importance of the 14-3-3 family proteins in the process of melanotumorigenesis, we next analyzed the protein-protein interactions of 14-3-3 family proteins with the GCMN and melanoma proteomes. We identified a total of 23 proteins, out of which 6 proteins in the GCNM and 17 proteins in the melanoma proteome possibly interacted with the 14-3-3 family proteins (Figure (Figure5B).5B). Two proteins interacting with 14-3-3, PHB and EEF1A1, were co-detected in both GCMN and melanoma proteomes. Additionally, the average interaction number of the 14-3-3 family proteins (7.66) was 2-fold higher than the average interaction number of other proteins (3.80). These results suggest that there are similar proteomic alterations in both GCMN and metastatic melanoma and the 14-3-3 family proteins may play an important role in melanotumorigenesis and tumor progression.
To further validate the proteomic results, we evaluated the altered expression of 14-3-3 alpha/beta, sigma, epsilon, zeta, tau, and of prohibitin in GCMN and unpaired normal skin samples (n=7) using semi-quantitative western blots. The integrated values of each protein band are reported in supplementary Additional files 4: Table S2 and Additional files 5.: Table S3 Because it was difficult to obtain skin specimens from normal children younger than 10 years of age, we could not use age-matched normal skin samples. However, we used comparatively young skin samples obtained from subjects whose average age was validated as 22.7 years. As seen in the representative western blots (Figures (Figures6A6A and B), the protein expression of 14-3-3 epsilon, 14-3-3 tau, and prohibitin in GCMN samples was 1.55-fold, 2-fold, and 2.53-fold, respectively, which was significantly higher than that in normal skin samples. However, the expression of 14-3-3alpha/beta, sigma, and zeta were not significantly different in GCMN and normal skin samples. One possible explanation for the unmatched results between the results of the proteomic and western blot analysis might be the diversity of samples used in the western blot analysis, which unlike the paired samples utilized in the proteomic analysis, were collected from normal skin samples of seven non-GCMN patients.
We further analyzed the expression levels of 14-3-3 epsilon, 14-3-3 tau, and prohibitin in normal skin fibroblast cell line (Detroit 551) and three kinds of melanoma cell lines (SK-MEL-2, SK-MEL-5, and SK-MEL-28) to validate whether our proteomic findings are truly relevant to clinical melanoma. Our results showed that the protein levels of 14-3-3 epsilon were significantly increased in all melanoma cell lines, and the levels of 14-3-3 tau were significantly increased in the SK-MEL-2 and SK-MEL-28 cell lines. Additionally, the protein levels of prohibitin were increased in the SK-MEL-2 and SK-MEL-5 cells, but were decreased in the SK-MEL-28 cell line (Figures (Figures7A7A and B).
In the present study, the proteomic composition of GCMN was compared with that of normal skin. A major aim of the study was the identification of proteins whose expression is altered in GCMN, which will help understand the altered biological processes in GCMN and help gain an insight into the mechanism of melanotumorigenesis in these malformations. LC-MS/MS analysis showed that 46 of the 438 identified proteins changed in their abundance levels between the normal skin and GCMN samples. In the GCMN samples, 92% of the abundance modified proteins were upregulated, but only 8% were downregulated (Figure (Figure22 and Table Table2).2). The use of different bioinformatic tools showed that GCMN clearly differed from normal skin in terms of protein expression patterns, which suggested that specific biological processes are altered in GCMN. As derived from the GO categories, KEGG pathways, and Reactome_biocarta, these processes were shown to encompass several major biological functions, namely the neurotrophin signaling pathway, downregulated of MTA-3 in ER-negative breast tumors, the cell cycle, phospholipase inhibitor activity, and glycolysis/gluconeogenesis. Strikingly, among these, neurotrophin signaling [17,18], MTA-3 downregulation (Table (Table3)3) , cell cycle deregulation , and glycolysis/gluconeogenesis  have been implicated in the development and progression of melanoma and other cancers.
Comparison of systemic properties of the GCMN and metastatic melanoma proteomes revealed that these two different disease proteomes shared at least five proteomic alterations in common and their abundance modified proteins closely interacted with each other (Figure (Figure5A).5A). Because closely related diseases are known to share common proteins or common interactions , our results suggested the close relationship between GCMN and melanoma.
Our proteomic analysis also revealed the significantly increased expression of 14 cancer-related proteins in GCMN compared to normal skin samples. Among them, PHB is a molecular maker of malignant cancers, and overexpression of PHB has been reported in melanoma [11,12] and various kinds of cancers, including gastric carcinoma , thyroid cancer , and hepatocellular carcinoma . This significant upregulation of cancer-related proteins in GCMN, specifically which of melanoma-implicated proteins, strengthened the possible risk of melanotumorigenesis in GCMN.
The 14-3-3 proteins comprise a highly conserved family of proteins whose members are found in both plants and mammals. They mediate signal transduction by binding to phosphoserine-containing proteins and are involved in many biological cellular processes, such as metabolism, protein trafficking, signal transduction, apoptosis, and cell cycle regulation, through interaction with various phosphoserine-containing proteins, such as CDC25 phosphatases, RAF1, and IRS1 proteins. In the present study, 14-3-3 family proteins were estimated to interact with 23 proteins in GCMN and melanoma (Figure (Figure5B),5B), and their average number of interactions was about 2-fold higher than the average number of interactions of other abundance modified proteins. These results suggested that 14-3-3 family proteins could play an important role in the alteration of biological processes in GCMN and melanoma.
The 14-3-3 family proteins consist of seven isoforms: beta, gamma, epsilon, sigma, zeta, tau, and eta. The alpha and sigma isoforms are the phosphoforms of 14-3-3 beta and zeta, respectively. All 14-3-3 proteins are ubiquitously expressed, with the exception of 14-3-3 sigma, which is exclusively expressed in epithelial cells . Among the 14-3-3 family members, the overexpression of the 14-3-3 sigma gene or its respective protein is frequently found in cancers such as ovarian carcinomas , pancreatic cancer , papillary thyroid carcinoma , hepatocellular carcinoma , and breast cancer . In our proteomic analysis, the expression of 14-3-3 proteins was significantly higher in GCMN than in normal skin samples, which strongly supports the greater tendency toward melanotumorigenesis in GCMN. In particular, the enhanced expression of 14-3-3 epsilon and tau proteins was clearly shown in western blot analysis (Figure (Figure66).
Compared to the other isoforms, little is known about the molecular and biological role of 14-3-3 tau and epsilon proteins. Like other isoforms, 14-3-3 tau is also involved in cell death and survival processes. For example, 14-3-3 tau binds to ataxia telangiectasia-mutated (ATM)-phosphorylated E2F1 during DNA damage and promotes E2F1 stability, leading to the induction of apoptosis , and the deletion of 14-3-3 tau leads to embryonic lethality in a mouse model . Interestingly, a recent study suggested that 14-3-3 tau exhibits an oncogenic role by downregulating p21 in breast cancer .
14-3-3 epsilon has been shown to play an essential role in cell development. Studies in Drosophila showed that 14-3-3 epsilon is required for the correct timing of mitosis in undisturbed post-blastoderm cell cycle . More recently, defects in neuronal migration during the development of 14-3-3 epsilon-knockout mice were reported .
The phosphorylation-induced binding of 14-3-3 epsilon to the pro-apoptotic transcription factor forkhead transcription factor-like 1 (FKHRL1 or FOXO3a) leads to structural changes in 14-3-3 epsilon and inhibits its pro-apoptotic activity . In inflammation and carcinogenesis, 14-3-3 epsilon interacts with key molecules of the mitogen-activated protein kinase signaling module to selectively modulate tumor necrosis factor-alpha-induced nuclear factor-kappa-beta activity . The function and regulatory mechanism of 14-3-3 epsilon in carcinogenesis is controversial and appears to be tumor-specific. Expression of the protein is higher in renal cell carcinoma than that in normal kidney . Moreover, on the basis of their involvement in the tumorigenesis of meningioma, 14-3-3 epsilon, zeta, and theta are thought to be efficient markers for predicting the degree of malignancy of these tumors . In contrast, mRNA and protein expression of 14-3-3 epsilon in laryngeal squamous cell carcinoma tissues was shown to be significantly lower than that in normal tissues . An early role of 14-3-3epsilon in tumorigenesis is suggested by the observation that 14-3-3 epsilon expression is increased in intrinsically aged and photoaged human skin . Interestingly, we found even higher protein levels of 14-3-3 epsilon, 14-3-3 tau, and PHB in GCMN than those in aged skin samples. This result suggested that GCMN may have a higher risk of tumorigenesis than aged skin. Because of the limitation in sample availability, we could not directly determine the expression level of 14-3-3 proteins and PHB in malignant melanoma tissue; however, we demonstrated significantly increased protein expression of 14-3-3 epsilon and tau in two different melanoma cell lines, SK-MEL-2 and SK-MEL-28, compared to normal skin cell line (Detroit 551). This result might support the association of 14-3-3 epsilon and tau upregulation with clinical melanotumorigenesis (Figures (Figures7A7A and B).
Nevertheless, further studies are needed to validate the functional role of 14-3-3 proteins in melanotumorigenesis through the proteomic comparison of different malignant melanoma patients with giant congenital melanocytic nevi. Furthermore, it is also necessary to carefully validate the biological meaning of the upregulation of melanoma-implicated proteins in GCMN and their role in melanotumorigenesis.
Taken together, our data suggest that proteomic modifications with tumorigenic potential are present in GCMN, and these proteomic alterations possibly modify six important biological processes or pathways that include melanosome, neurotrophin signaling pathway, downregulated of MTA-3 in ER-negative breast tumors, cell cycle, phospholipase inhibitor activity, and glycolysis/gluconeogenesis These pathways may be significantly altered in GCMN skins. The intensive alteration of 14-3-3 family proteins and PHB possibly acts as a central regulator of GCMN biological pathway remodeling, which may have an important role in the development of GCMN and could be associated with melanotumorigenesis.
A total of 10 normal and GCMN skin samples, which were defatted, were obtained from patients who underwent excision procedures at the Department of Plastic Surgery, Inje University Ilsan Paik Hospital, Korea. The collection and use of the samples were approved by the Institutional Review Board of Inje University Ilsan Paik Hospital (IRB No. IB-0902-015). The present study was carried out in accordance with The Code of Ethics of the World Medical Association (Declaration of Helsinki) for experiments involving humans.
The human embryo skin cell line Detroit 551 and the human malignant melanoma cell lines SK-MEL-2, SK-MEL-5, and SK-MEL-28 were obtained from the American Type Culture Collection (ATCC; Rockville, MD). The culture medium used throughout these experiments was RPMI-1640 (Lonza, Verviers, Belgium) containing 10% fetal bovine serum (PAA, Pasching, Austria) and 100 μg/ml penicillin-streptomycin (Lonza). The cells were incubated at 37°C in a humidified atmosphere of 5% CO2.
Three paired normal and GCMN skin samples were selected for 1D-LC-MS/MS proteomic analysis to exclude environmental bias. Excised skin samples were ground to a powder in liquid nitrogen, dissolved in lysis buffer (9 M urea, 2 M thiourea, 4% CHAPS (3-[(3-cholamidopropyl)dimethylammonio]-1-propanesulfonate), 40 mM dithiothreitol (DTT), and 1% protease inhibitor cocktail), vortexed, and incubated on ice for 1 h. The mixture was then centrifuged (10,000×g, 30 min, 4°C), and the total proteins contained in the supernatant were used for the experiments. The total protein content of the solution was determined using the 2D Quant kit (GE Healthcare, Milwaukee, WI), with bovine serum albumin (0–50 mg/ml) as the standard.
Protein separation and LC-MS analysis were performed as previously described . Briefly, dissolved skin proteins were separated on a 12% polyacrylamide gel by SDS-PAGE. The gels were washed three times with ddH2O for 5 min each and stained with Bio-Safe Coomassie stain solution (Coomassie G250 stain; Bio-Rad, Hercules, CA) for 1 h, with gentle shaking at room temperature. The Coomassie-stained gels were evenly sliced into 15 slices and then destained by incubation in 75 mM ammonium bicarbonate/40% ethanol (1:1). Disulfides were reduced by treatment with 5 mM DTT/25 mM ammonium bicarbonate at 60°C for 30 min, followed by alkylation with 55 mM iodoacetoamide at room temperature for 30 min. The gel pieces were then dehydrated in 100% acetonitrile (ACN), dried, and swollen overnight at 37°C in 10 μl 25 mM ammonium bicarbonate buffer containing 20 μg modified sequencing-grade trypsin (Roche Applied Science, Indianapolis, IN) per ml. The tryptic peptide mixture was eluted from the gel using 0.1% formic acid. LC-MS/MS analysis was performed using a ThermoFinnigan ProteomeX workstation LTQ linear ion trap MS (Thermo Electron, San Jose, CA) equipped with a nanospray ionization (NSI) source (Thermo Electron). Briefly, 12 μl peptide sample obtained from the in-gel digestion was injected and loaded onto a peptide trap cartridge (Agilent, Palo Alto, CA). Trapped peptides were eluted onto a 10-cm reversed-phase PicoFrit column packed in-house with 5-μm, 300-Å pore size C18 and separated by gradient elution. The mobile phases consisted of H2O and ACN, both containing 0.1% v/v formic acid. The flow rate was maintained at 200 nl/min. The gradient started at 2% ACN, then reached 60% ACN in 50 min, 80% ACN in the next 5 min, and 100% H2O in the final 15 min. Data-dependent acquisition (m/z 400–1800) was enabled, and each MS survey scan was followed by five MS/MS scans within 30 s, with the dynamic exclusion option enabled. The spray voltage was 1.9 kV, the temperature of the ion transfer tube was 195°C, and the normalized collision energy was 35% .
Data-analyzed tandem mass spectra were extracted, and the charge state was deconvoluted and deisotoped using the Sorcerer 3.4 beta2 platform (Sorcerer software 3.1.4, Sorcerer Web interface 2.2.0 r334, and Trans-, Proteomic Pipeline 2.9.5). All MS/MS samples were analyzed using SEQUEST (version v.27, rev. 11; ThermoFinnigan, San Jose, CA), which was set to search the ipiHuman 3.29 database (IPI ver.3.29, 40131 entries), with semitrypsin as the digestion enzyme. The search used a fragment-ion mass tolerance of 1.00 Da and a parent-ion mass tolerance of 1.5 Da. Iodoacetamide-derivatized cysteine was specified as a fixed modification. Methionine oxidation, iodoacetamide derivatizion of cysteine, and phosphorylation of serine, threonine, and tyrosine were specified as variable modifications. The Scaffold software (version Scaffold-2.0; Proteome Software Inc., Portland, OR) was used to validate MS/MS-based peptide and protein identifications. Peptide identifications were accepted if their probability was >95.0%, as specified by the Peptide Prophet algorithm, and if they contained at least one identified peptide. Protein probabilities were assigned by the Protein Prophet algorithm. Proteins containing similar peptides such that they could not be differentiated based on MS/MS analysis alone were grouped to satisfy the principles of parsimony. After identifying the proteins, each dataset was used for a subtractive analysis by semi-quantitative normalized spectral counts, which were normalized by total spectral counts in the Scaffold program .
A systemic bioinformatics analysis of the GCMN proteome was conducted using the Search Tool for the Retrieval of Interacting Genes/Proteins (STRING 8.3) , the Protein Analysis Through Evolutionary Relationships classification system (PANTHER 7.0) , the National Center for Biotechnology Information (NCBI) COG database , Cytoscape, and ClueGO .
Protein expressions of 14-3-3 alpha+beta (30 kDa), 14-3-3 epsilon (29 kDa), 14-3-3 zeta (28 kDa), 14-3-3 sigma (25 kDa), 14-3-3 tau (31 kDa), and prohibitin (30 kDa) in normal (n=7) and GCMN (n=7) samples were analyzed by western blots to confirm the proteomic results. In addition, protein expressions of 14-3-3 epsilon, 14-3-3 tau, and prohibitin in a normal cell line (Detroit 551) and melanoma cell lines (SK-MEL-2, SK-MEL-5, and SK-MEL-28) were analyzed by western blots. Relative expression of each protein was normalized to an internal standard protein, β-tubulin (55 kDa) or glyceraldehyde-3-phosphate dehydrogenase (37 kDa) (Abcam, Cambridge, MA). The values were expressed as mean±standard error.
Two-tailed unpaired Student’s t-test in the Scaffold (version 2.06.02) software was used to compare abundance of each protein in the normal and GCMN skin samples. Using the Origin software (version 7.0220, OriginLabs, MA, USA), Bonferroni correction was applied to control the rate of false positives in the comparison of means of each protein’s abundance. A two-tailed unpaired Student’s t-test was used to compare the results of western blot analysis between the normal and GCMN skin samples. p<0.05 was considered statistically significant.
GCMN: Giant congenital melanocytic nevus; CTSD: cathepsin D; PDIA3: Protein disulfide-isomerase A3 precursor; PHB: Prohibitin; HSPA8: Heat shock protein 70; PHGDH: D-3-phosphoglycerate dehydrogenase; ATP5B: ATP synthase subunit beta; KRT1: Cytokeratin-1; FLG2: Filaggrin; HRNR: Hornerin; ADH1B: Alcohol dehydrogenase 1B; CSPG4: Chondroitin sulfate proteoglycan 4; PEBP1: Phosphatidylethanolamine-binding protein; RPSA: Ribosomal protein SA; COG: Clusters of Orthologous Groups; STRING: Search Tool for the Retrieval of Interacting Genes/Proteins; PANTHER: Protein Analysis Through Evolutionary Relationships; KEGG: Kyoto Encyclopedia of Genes and Genomes; GO: Gene Ontology.
The authors declare that they have no competing interests.
HKK performed the proteomic analysis including its design, coordination, analysis of the data, and drafted the manuscript. YKK provided GCMN and normal skin samples and was largely involved in the design of the study and in writing the manuscript. ISS, SRL, SHJ, MHK, and DYS performed sample preparation and western blot analysis. NK, BDR, KSK, KCT, and CGP critically revised the clinical aspect of the manuscript. JYC carried out the MS/MS proteomic analysis. JH conceived the overall experimental design and manuscript preparation. All authors read and approved the final manuscript.
Figure S1. Peptide mass peak of 15 slices dissected from 1D gels of normal skin and GCMN.
Table S1. Identified proteins in skin samples of Normal and GCMN patients.
Figure S2. The functional network and enriched functional group of GCMN was analyzed using ClueGO, a biological term enrichment analyzer. The following proteins were significantly enriched in GCMN: melanosome (GO_cellular component), neurotrophin signaling pathway (KEGG pathway), downregulated of MTA-3 in ER-negative breast tumors (Biocarta), cell cycle (KEGG pathway), phospholipase inhibitor activity (GO_molecular function), and glycolysis/gluconeogenesis (KEGG pathway).
Table S2. Integrated densitometry value of Western blot band.
Table S3. Statistical analysis result of western blot.
This work was supported by Priority Research Centers Program and Basic Research Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Education, Science and Technology (353-2009-1-E00002 and 2010–0020224).