Quantitative proteome analysis of Merkel cell carcinoma cell lines using SILAC

Background Merkel cell carcinoma (MCC) is an aggressive neuroendocrine tumour of the skin with growing incidence. To better understand the biology of this malignant disease, immortalized cell lines are used in research for in vitro experiments. However, a comprehensive quantitative proteome analysis of these cell lines has not been performed so far. Methods Stable isotope labelling by amino acids in cell culture (SILAC) was applied to six MCC cell lines (BroLi, MKL-1, MKL-2, PeTa, WaGa, and MCC13). Following tryptic digest of labelled proteins, peptides were analysed by mass spectrometry. Proteome patterns of MCC cell lines were compared to the proteome profile of an immortalized keratinocyte cell line (HaCaT). Results In total, 142 proteins were upregulated and 43 proteins were downregulated. Altered proteins included mitoferrin-1, histone H2A type 1-H, protein-arginine deiminase type-6, heterogeneous nuclear ribonucleoproteins A2/B1, protein SLX4IP and clathrin light chain B. Furthermore, several proteins of the histone family and their variants were highly abundant in MCC cell lines. Conclusions The results of this study present a new protein map of MCC and provide deeper insights in the biology of MCC. Data are available via ProteomeXchange with identifier PXD008181.


Background
Merkel cell carcinoma (MCC) is a rare malignant tumour of the skin with neuroendocrine differentiation [1,2] and growing incidence rates ranging from 2 to 4 cases per million per year in Europe and the US, to 8 cases per million per year in Australia [3]. MCC shows a very aggressive behaviour with significant potential to build metastases and a high locoregional recurrence rate [4,5]. The overall 10-year survival is reported to be 57.3% [6]. The main risk factors are UV radiation as MCC mainly appears in sun-exposed areas, and immunosuppression since the incidence is higher in HIV-infected patients, transplant recipients and patients with chronic lymphocytic leukaemia [7]. A further factor that plays an important role in the development of MCC is the Merkel cell polyomavirus [8]. Firstly described by Feng et al. in 2008 it can be found in up to 80% of the cases [9,10]. However other biological processes are likely to be involved in the development of MCC.
A genomic profiling study of patients with MCC revealed that the most frequent abnormalities are related to the TP53 gene and the cell cycle pathway. Further abnormalities were found in the PI3K/AKT/mTOR pathway and DNA repair genes [11]. However, despite of changes at the chromosomal level, the dysfunction of biochemical pathways is expressed at the protein level. Therefore this proteomic study was conducted to gain deeper insights into the biology of MCC and possibly

Open Access
Clinical Proteomics *Correspondence: goran.mitulovic@meduniwien.ac.at 4 Clinical Department of Laboratory Medicine, Medical University of Vienna, Waehringer Guertel 18-20, 1090 Vienna, Austria Full list of author information is available at the end of the article to find new molecular targets for therapy. Stable isotope labelling by amino acids in cell culture (SILAC) was applied to investigate the expression patterns of six MCC cell lines.

Cell culture
The human Merkel cell carcinoma cell lines BroLi, MKL-1, MKL-2, PeTa, WaGa, and MCC13 were a kind gift of Prof. Houben [12]. HaCaT, a human skin keratinocyte cell line was obtained from AddexBio (San Diego, CA, USA). Cells were cultured in RPMI buffer (Thermo Fisher Scientific, Waltham, MA, USA) supplemented with 10% fetal bovine serum (Thermo Fisher Scientific, Waltham, MA, USA) at 37 °C and 5% CO 2 in a humidified incubator. For SILAC labelling cells were grown in SILAC Media supplemented either with 13 C 6 l-Lysine-2HCl (heavy) and l-Arginine-HCl (light) or with l-Lysine-2HCl (light) and l-Arginine-HCl (light) (Thermo Fisher Scientific, Waltham, MA, USA). Cells were cultured for at least ten cell doublings. Three biological replicates were cultured per each cell line, and these were merged into one sample upon cell lysis for further treatment. The labelling efficiency was estimated using the method described by Rappsilber et al. [13]. Briefly, incorporation efficiency of the heavy labeled amino acids into proteins was assessed in a pilot experiment, where a small aliquot of cells was lysed, and proteins were reduced, alkylated, and tryptically digested. The resulting peptides were subjected to MS analysis as described below. Heavy label incorporation into proteins obtained from cells was assessed to be more than 95%.

Sample preparation and protein identification, quantification and analysis Cell lysis and protein digests
MCC cells were lysed using the "Chemicon ® -Total Protein Extraction Kit" containing TM buffer (HEPES, pH7.9, MgCl2, KCl, EDTA, Sucrose, Glycerol, Sodium deoxycholate, NP-40, Sodium Ortho Vanadate, Merck Millipore Vienna, Austria) according to manufacturer's manual. Protein precipitation was performed using modified Wessel-Fluegge method as described elsewhere [14][15][16][17]. Lysed cell content was treated with methanol and dichloromethane and the interphase was collected. Proteins were pelleted by addition of methanol, air dried and dissolved in 50 mM triethylammonium bicarbonate (Sigma-Aldrich, Vienna, Austria). Protein content was measured using the Direct Detect FT-IR spectrophotometer (Merck Millipore Vienna, Austria) and equimolar amounts of proteins, 1 µg total proteins, from control and the MCC cells were mixed and submitted for tryptic digestion. Protein digest was performed using sequencing grade trypsin (Promega, Mannheim, Germany) as described earlier [17]. Briefly, proteins were reduced with 5 mM DTT for 30 min at 60 °C, and alkylated for 30 min with 15 mM iodoacetamide in the dark. Finally, porcine trypsin was added in a ratio 1:50 (w/w). After 16 h of incubation at 37 °C, aliquots of 20 µl were prepared and stored in 0.5 ml protein low-bind vials at − 20 °C until further usage. Iodoacetamide and DTT were purchased from Sigma-Aldrich (Vienna, Austria). All extractions and digestion steps were performed in protein low-bind vials of different volumes (Eppendorf, Vienna, Austria).

Nano-HPLC separation and MS detection
Methanol, was purchased from Merck (Vienna, Austria), 98% formic acid, acetonitrile, trifluoroacetic acid were purchased from Sigma-Aldrich (Vienna, Austria). HPLC grade water was prepared using an in-house Milli-Q plus device from Millipore (Vienna, Austria), and trifluoroethanol was purchased from Alfa Aesar (Karlsruhe, Germany). Reverse phase separation of tryptic peptides was conducted on a nanoRSLC UltiMate3000 (Thermo Fischer Scientific, Vienna, Austria) HPLC system consisting of an autosampler, thermal compartment and pumping module. PepMap C18 trap-column (300 µm ID × 5 mm length, 5 µm particle size, 100 Å pore size, Thermo Fisher Scientific, Vienna, Austria) was used for sample loading and desalting. The analytical column used for peptide separation was a 75 µm ID × 50 cm length Acclaim ® PepMap100 (C18, 3 µm particle size, 100 Å pore size, Thermo Fisher Scientific, Vienna, Austria). Both columns were operated at 60 °C in the column compartment as described earlier [15]. A total of 500 ng digested protein was injected onto the trap column at 30 µl/min and was loaded using the loading pump flow rate of 30 µl/ min. After 10 min the switching valve changed the position and the trap column was switched into the flow path of the nano pump. Peptides were eluted from the trap column onto the separation column using the gradient described in Supplementary information and in Table 1. HPLC was hyphenated with the maXis Impact Q-Time-of-Flight mass spectrometer (Bruker, Bremen, Germany) equipped with a nano ESI captive spray source. Peptides were ionized using positive electrospray and Data-Dependent collision-induced-dissociation was used for peptide fragmentation (MS/MS data). MS data was acquired using the data-dependent mode with positive ionization. Capillary was set to 1.8 kV and 20 most intense ions were fragmented using the collisional induced fragmentation by ramping the collisional energy from 15 to 35 eV. Fragmented ions were excluded from further fragmentation for 60 s.
All raw data were converted into Mascot "mgf " files by using Bruker's Data Analysis and these files were  All samples were analysed as technical triplicates to ensure statistical sound data and avoid artefacts due to variations in ionization efficiency.
Detailed information on separation gradient, the MS settings, and the data search and quantitation can be found in supplemental information (Additional file 1).

Term enrichment analysis
Heavy/light ratios where calculated using WARP-LC v. 1.3 (Bruker, Bremen, Germany). Proteins in MCC cell lines with H/L ratios of > 1.5 or < 0.5 where considered as significantly differently abundant. To put differential protein abundance into biologic context, Cytoscape (Seattle, WA, USA) in combination with ClueGO/CluePedia (a Cytoscape plug-in) was used with default parameters except for following: Database Gene Ontology Biological Process, levels between 4 and 13 and GO Fusion set on true [18]. p-values where corrected for multitesting according to Benjamini Hochberg.

Differentially expressed proteins in MCC cell lines and the control cell line
The SILAC method was used in six MCC cell lines to determine quantitative changes of proteins at the proteome level. Proteins detected in MCC cells were compared to the reference cell line HaCaT (Fig. 1). We chose the keratinocyte cell line HaCaT as reference cell line since, to our best knowledge, there is no commercially or otherwise available cell line with healthy Merkel cells.
In order to visualize the similarity and the difference between the particular cell lines, a heat map was created showing all quantified proteins. A list of proteins (overlapping and specific) for each sample is provided as Additional file 2: Table S1. As seen in Fig. 2 every cell line has its own distinct protein abundance pattern. The most similar cell line compared to the control cell line HaCaT was BroLi, whereas WaGa differed significantly from the other cell lines. All cell lines originated from different old patients and different anatomic locations. While WaGa was derived from ascites of a 67 years old man, MKL-1 was derived from a nodal metastasis of a 26 years old man. MKL-2 stem from a 72 years old man and the localization is unknown. BroLi was obtained from pleural effusion of a 55 year old man [19]. MCC13 was gained from a nodal metastasis of a 80 year old female patient and is called in literature also "variant" MCC cell line since unlike BroLi, MKL-1, MKL-2, PeTa and WaGa, it is a Merkel cell polyomavirus negative cell line and lacks some typical markers in immunohistochemical staining [20]. Nevertheless we decided to include this cell line into our study since a number of studies in the field of MCC research are still performed using this particular cell line.
Next, the protein profile of MCC cell lines was compared with the reference cell line HaCaT. Proteins present specifically in the MCC cells were determined. Then a Venn diagram was constructed. Since BroLi was the cell line with the least number of proteins and the difference between HaCaT and BroLi was small, the cell line BroLi has been omitted. Figure 3 shows a Venn diagram with the specific proteins for the cell lines MKL-1, MKL-2, PeTa, WaGa, and MCC13. Remarkably, only 10 proteins were found in all five cell lines at the same time: alpha 2-HS glycoprotein, inter-alpha-trypsin inhibitor heavy chain 2, FUS RNA binding protein, mechanistic target of rapamycin, SUB1 homolog transcriptional regulator, Y-box binding protein 1, serine and arginine rich splicing factor 2, testis specific 10 interacting protein, sperm associated antigen 5 and heterogeneous nuclear ribonucleoprotein A/B. A complete list of all specific proteins is provided as Additional file 3: Table S2.

Term enrichment analysis of proteins
It is of crucial importance and of highest interest to identify and quantify biological processes involved in the biology of cancer. Term Enrichment Analysis using ClueGO showed that multiple pathways where affected by differentially represented proteins (Fig. 4). Cellular processes like metabolic processes, protein folding, and signal transductions were affected. In particular, viral transcription was present in all cell lines but mostly in MKL-2. This can be explained by the fact that the Merkel cell polyomavirus has an important function in the pathogenesis of the development of MCC [21]. In the cell line MKL-2 also several mRNA and rRNA processes were more prevalent. In MCC13 the spliceosomal complex assembly was very active together with filament cytoskeleton organization and regulation of cell death. Further processes that play a role in cancer cell motility, like regulation of actin filament depolymerization [22] were enriched in several cell lines.

Overexpression of multiple proteins in different MCC cell lines
In total, 317 dysregulated (i.e. up-and downregulated) proteins with significance threshold of p < 0.05 were identified. Proteins altered > 1.5-fold were considered as upregulated and proteins altered < 0.5-fold were considered as downregulated. Based on these criteria, 142 proteins were identified as upregulated and 43 proteins were  were identified in the BroLi cell line. BroLi cell line is a very slowly growing cell line with a doubling time of 5 days [19] and this could be the reason why only a limited number of proteins were identified. For the BroLi cell line, mitoferrin-1 was found to be 9.64-fold upregulated compared to HaCaT cell line. Mitoferrin-1 is a protein involved in the mitochondrial iron transport and storage [23]. As iron is an important co-factor in DNA synthesis, dysregulated iron metabolism in cells is believed to play a role in tumorigenesis. The disturbance in iron transport between cytosol and mitochondrion is thought to lead to mitochondrial dysfunction and it therefore may contribute to tumour formation and propagation [24].
In MKL-2 cells, protein-arginine deiminase type-6 was upregulated 9.47-fold compared to the control cell line making it the most differently regulated protein for this cell line. This protein is an enzyme involved in post-translational modifications, which can have substantial effects on the structure and function of proteins. Citrullination is one such post-translational modification being catalysed by the family of protein arginine deiminase (PADs) enzymes. Five isoenzymes (PAD1-4 and 6) are known and they were identified in different types of tissue [25]. An overexpression of PADs has been detected in diseases like rheumatoid arthritis, neurologic diseases and cancer. In particular, the overexpression of PAD4 is associated with cancer since it plays a role in histone citrullination [26]. We identified PAD6, an isoenzyme mainly found in oocytes and embryos, to be the most abundant protein in the cell line MKL-2. Although the relation of PAD6 and cancer has not been described in the literature so far, we assume that it can be of interest due to its high occurrence. Furthermore, we identified the heterogeneous nuclear ribonucleoprotein A2/B1 (hnRNPA2/B1) to be the most upregulated protein in the PeTa cell line. The hnRNPs are a group of proteins binding to RNA and playing a role in mRNA processing [27]. So far, hnRNPA2/B1 was found to be overexpressed in lung cancer were it promotes tumour growth by activation of COX-2 signalling [28,29]. Furthermore, it was also found to be up regulated in hepatoma cell lines, gastric cancer, breast cancer and glioblastoma [30][31][32] but it was not described for Merkel cell carcinoma yet.
In the WaGa cell line, SLX4IP (SLX4 interacting protein) was the protein showing the highest upregulation. SLX4 is a DNA repair protein and it coordinates structure-specific endonucleases [33] but its role in cancer has not been determined up to now.
Finally, clathrin light chain B was the protein with the highest upregulation in the MCC13 cell line. Clathrin light chain B is a part of the clathrin protein, which is the main component of vesicles involved in intracellular transport. Recently, it was reported that clathrin light chains promote cell migration and therefore may play a role in cancer metastasis [34].
In the PeTa cell line various histones and their variants were found to be dysregulated compared to the control cell line HaCaT. Histones are substantial components for the packaging of the DNA in the chromosomes. The smallest packaging units are nucleosomes consisting of DNA wrapped around a histone octamer. A histone octamer in turn consists of two copies of each of the core histones: H2A, H2B, H3, and H4, being the smallest units. The linker histone H1 holds the nucleosome together and is the fifth member of the histone protein family [35]. Beside their structural function, histones play an important role in DNA replication and transcription regulation. Recently, it became evident that changes in histone expression are associated with cancer since an altered nucleosome structure can lead to instability and accessibility for different transcription factors [36]. So far, most of the histone variants were found in the histone H1, H2A, H2B, and H3 family. Some variants have been studied more detailed, but for many variants the function is still not known [37]. Furthermore, some histones serve as markers for cellular proliferation. In case of MCC, Henderson et al. used H3KT (histone-associated mitotic marker H3K79me3T80ph) and PHH3 (phosphohistone H3) as surrogates for detecting mitotic figures. Detection of H3KT and PHH3 correlated with a worse overall survival [38].
In the current study, proteins from all five major histone families with 15 different subfamily members were differently abundant in MCC cell lines compared to control samples. In particular, H2A1H and H2B1O were found to be overexpressed in MKL-1 and H2B1H and H2B2E were upregulated in MKL-2. Most of the highly abundant histone variants were found in the cell line PeTa: H2A1B, H2AJ, H1.4, H3.3, H1.5, H2B1N, and H2B2E. Furthermore, H2B1O and H2B1D were identified to be upregulated in WaGa, and H2B1C, H2B1K, H2B2F and H2AJ in MCC13. The role of histone variants in the development of carcinomas has been discussed and described in a number of publications [39][40][41][42][43][44][45][46]. As described in a recent review, canonical histones can be replaced with variant histones after environmentalstress-induced DNA damage repair, which subsequently results in a change in chromatin structure and stability [47]. A well-known environmental-stress factor is UV radiation, which in turn is a recognized risk factor for the development of Merkel cell carcinomas. This study shows for the first time that histone variants play an important role in the biology of Merkel cell carcinomas. A large number of histone variants was identified in all examined cell lines, except for the BroLi cell line. Due to the small difference between HaCaT and BroLi, BroLi has been omitted. Few proteins were found in several cell lines at the same time and only 10 proteins were found to be common in all cell lines. These 10 proteins are shown in the middle of the Venn diagram (alpha 2-HS glycoprotein, inter-alpha-trypsin inhibitor heavy chain 2, FUS RNA binding protein, mechanistic target of rapamycin, SUB1 homolog transcriptional regulator, Y-box binding protein 1, serine and arginine rich splicing factor 2, testis specific 10 interacting protein, sperm associated antigen 5 and heterogeneous nuclear ribonucleoprotein A/B). A complete list of all proteins identified is provided in Additional file 3: Table S2 Another interesting group of proteins that were identified as dyregulated in several MCC cell lines were the heat shock proteins (HSPs). It is a group of proteins that inhibit the unfolding or denaturation of cellular proteins and therefore being known as molecular chaperones whose expression is induced by stress. The major groups are classified according to their sizes and imply HSP10, HSP27, HSP40, HSP60, HSP70, and HSP90. Recent studies have shown that HSPs are highly expressed in many malignant tumours and due to their important role in cell proliferation and differentiation they are involved in carcinogenesis and metastasis [48,49]. In case of MCC, presence of HSP70 seems to be necessary for the interaction of large T antigen and the tumour-suppressing retinoblastoma protein. In detail, the large T antigen is an oncoprotein expressed by polyomavirus affected cells and Merkel cell carcinoma in turn is highly associated with polyomavirus [9]. Binding of large T antigen to retinoblastoma protein leads to inactivation of retinoblastoma protein [50] and subsequently to cell proliferation via activation of cell cycle progression associated genes [51].
Beside HSP70, we found HSP60 and co-chaperone HSP10 to be overexpressed in the tested MCC cell lines. Actively produced by cancer cells, HSP60 exhibits a protective effect against cell stressors like chemotherapeutics. In particular, HSP60 stabilizes the anti-apoptotic protein survivin, a protein over-expressed in most human tumours, and therefore it inhibits apoptosis. Furthermore, HSP60 builds a complex with p53, which leads to the loss of the pro-apoptotic function of p53 and this process again results in inhibition of apoptosis [52].