Next Article in Journal
Multifunctional Peptides from Spanish Dry-Cured Pork Ham: Endothelial Responses and Molecular Modeling Studies
Next Article in Special Issue
Long-Term Vemurafenib Exposure Induced Alterations of Cell Phenotypes in Melanoma: Increased Cell Migration and Its Association with EGFR Expression
Previous Article in Journal
Hydrogen Sulfide Prevents Elastin Loss and Attenuates Calcification Induced by High Glucose in Smooth Muscle Cells through Suppression of Stat3/Cathepsin S Signaling Pathway
Previous Article in Special Issue
Gene-Specific Intron Retention Serves as Molecular Signature that Distinguishes Melanoma from Non-Melanoma Cancer Cells in Greek Patients
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Physiologically Relevant Oxygen Concentration (6% O2) as an Important Component of the Microenvironment Impacting Melanoma Phenotype and Melanoma Response to Targeted Therapeutics In Vitro

by
Marta Osrodek
,
Mariusz L. Hartman
and
Malgorzata Czyz
*
Department of Molecular Biology of Cancer, Medical University of Lodz, 6/8 Mazowiecka Street, 92-215 Lodz, Poland
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2019, 20(17), 4203; https://doi.org/10.3390/ijms20174203
Submission received: 26 July 2019 / Revised: 21 August 2019 / Accepted: 23 August 2019 / Published: 27 August 2019
(This article belongs to the Special Issue Molecular Biology of Melanoma)

Abstract

:
Cancer cell phenotype largely depends on oxygen availability. The atmospheric oxygen concentration (21%) used in in vitro studies is much higher than in any human tissue. Using well-characterized patient-derived melanoma cell lines, we compared: (i) activities of several signaling pathways, and (ii) the effects of vemurafenib and trametinib in hyperoxia (21% O2), normoxia (6% O2) and hypoxia (1% O2). A high plasticity of melanoma cells in response to changes in oxygen supplementation and drug treatment was observed, and the transcriptional reprograming and phenotypic changes varied between cell lines. Normoxia enhanced the expression of vascular endothelial growth factor (VEGF), glucose metabolism/transport-related genes, and changed percentages of NGFR- and MITF-positive cells in cell line-dependent manner. Increased protein stability might be responsible for high PGC1α level in MITFlow melanoma cells. Vemurafenib and trametinib while targeting the activity of MAPK/ERK pathway irrespective of oxygen concentration, were less effective in normoxia than hyperoxia in reducing levels of VEGF, PGC1α, SLC7A11 and Ki-67-positive cells in cell line-dependent manner. In conclusion, in vitro studies performed in atmospheric oxygen concentration provide different information on melanoma cell phenotype and response to drugs than performed in normoxia, which might partially explain the discrepancies between results obtained in vitro and in clinical settings.

1. Introduction

As our understanding of the role of tumor microenvironment (TME) in cancer biology and disease progression has expanded, numerous attempts to mimic TME in in vitro settings have been made. These include three-dimensional matrices [1,2], co-cultures with other cell types [3,4] and modifications of culture media supplementation [5,6,7,8]. However, in spite of oxygen being an essential biological component, its concentration characteristic for human tissues is misrepresented in standard in vitro research models and majority of studies in vitro are performed in an atmospheric oxygen concentration (21% O2) [9]. Mammalian cells are exposed to up to 14.5% oxygen in vivo, and distinct cell functions are attributed to tissue-specific oxygen levels [10]. Furthermore, there is accumulating evidence of the detrimental effect of the atmospheric oxygen level in cell culture, particularly affecting primary, non-transformed cells and stem cells [11,12,13,14]. Although cancer cells are more resilient and are easily propagated ex vivo, oxygen remains a crucial, regulatory factor in multiple aspects of cancer biology [9]. It is therefore not only intuitive but rational that ex vivo studies on mammalian cells should be conducted at oxygen concentrations that are more appropriate for the tissue they originated from.
Melanoma arises from malignant transformation of melanocytes residing in the epidermal layer of skin, where the oxygen concentration does not exceed 8% O2 [9,15]. Most melanomas harbor activating mutations in genes associated with the RAS/RAF/MEK/ERK pathway, with BRAFV600E being the most prevalent [16]. Efforts to target this pathway led to the development of selective therapeutics against BRAFV600E or MEK1/2 kinase, however, the vast majority of melanomas are either intrinsically resistant or develop drug resistance within the first year of treatment [17,18,19,20]. Despite being rather small in size, melanoma tumors, like other solid cancers, comprise poorly oxygenated areas of hypoxia, which prompted multiple studies focused on the role of hypoxic conditions in modifying melanoma phenotype and melanoma response to treatment [21,22,23,24]. Little is known, however, how higher but still physiological oxygen concentrations affect melanoma cell phenotype and the predominant research approach in all oxygen-related studies is a short-term exposure of melanoma cells to low (≤1%) oxygen concentrations following standard hyper-oxygenated cell culture.
In recent years, metabolic reprogramming has gained significant interest. In spite of previous reports showing that cancer cells are primarily glycolytic, it has become clearer that melanomas manifest great metabolic plasticity and heterogeneity. While enhanced glycolysis favors metastatic behavior in melanoma cells and several elements of the glycolytic pathway are considered as prognostic markers [25,26], some melanoma cells characterized by high expression of the peroxisome proliferator-activated receptor γ coactivator 1α (PGC1α), exhibit greater dependency on oxidative phosphorylation [27]. Moreover, results showing contradictory metabolic phenotypes have been obtained between in vitro and in vivo studies. For instance, while ovarian tumor-initiating cells are highly glycolytic in vivo they display a strong oxidative metabolism in vitro [28].
Three melanoma cell lines derived from surgical specimens of previously untreated patients, DMBC12, DMBC17, and DMBC28, were used in this study. They are genetically and phenotypically different, as published in detail elsewhere [29,30,31]. Briefly, two of the cell lines, DMBC12 and DMBC28, harbor V600E substitution in BRAF and one, DMBC17, harbors Q61R substitution in HRAS. As the microphthalmia-associated transcription factor (MITF) is one of the major regulators of mitochondrial respiration in the melanocytic lineage [27], we have chosen melanoma cell lines that differ in the expression and activity of MITF [31], with the highest in DMBC17 cells and almost undetectable in DMBC12 cells. When a multi-stage differentiation model [32] is applied to classify melanoma cells chosen for this study, and the expression of the sex determining region Y-box 10 (SOX10) and tyrosine-protein kinase receptor (AXL) is considered in addition, each cell line represents a different phenotype: DMBC12 cell line (MITFlow/no/SOX10high/AXLhigh) can be considered as invasive/neural crest-like, DMBC28 cell line (MITFhigh/SOX10high/AXLlow) as neural crest/pigmentated, and DMBC17 cell line (MITFhigh/SOX10high/AXLlow) as melanocyte-like [31].
We established an in vitro culture model of patient-derived melanoma cells in a physiological oxygen concentration corresponding to human peripheral tissues (6% O2). These conditions were defined in this study as normoxia, whereas standard laboratory conditions (21% O2) were designated as hyperoxia. We aimed to compare the phenotypes of investigated melanoma cell lines in different concentrations of oxygen. The influence of two drugs currently used in melanoma treatment, vemurafenib, a BRAFV600E inhibitor and trametinib, a MEK1/2 inhibitor, was also compared in these conditions.

2. Results

2.1. In Vitro Culture of Melanoma Cells in 6% O2 Alters their Metabolic Demands

First, we investigated proliferation of melanoma cells and observed that initial changes in cell confluence were similar in hyperoxia and normoxia, however, after approximately 30 h a strong cell density-dependent inhibition was detected in normoxia (Figure 1A). Moreover, there was a pronounced reduction of cell survival after three days in 6% O2, and Annexin V/propidium iodide staining revealed that apoptotic and necrotic cells could be counted, especially in DMBC17 and DMBC28 cell line cultures (Figure 1B). While the percentages of apoptotic and necrotic cells were low in hyperoxia, even at initially high density of cells reaching 3 × 105/mL, in normoxia both apoptosis and necrosis were induced already at lower seeding density of 2 × 105 cells/mL (Figure 1B). This might indicate that melanoma cells differ in metabolic demands when cultured in the presence of 6% and 21% of oxygen.

2.2. Oxygen Concentration-Dependent Changes in the Activity of Signaling Pathways are Diverse among Melanoma Cell Lines

Oxygen plays a fundamental role in tumor biology, and several signaling molecules respond to oxygen deficit. In three different oxygen concentrations, we assessed closely interconnected hypoxia-inducible factor (HIF)-1α and HIF-2α as the key regulators of the cellular response to oxygen deprivation, phosphorylation status of the pyruvate dehydrogenase (PDH) crucial for conversion of pyruvate into acetyl coenzyme A (acetyl-CoA), and the activity of the PI3K/AKT (phosphoinositide 3- kinase/protein kinase B) pathway promoting HIF-1α activity. HIF-1α and HIF-2α, whose protein levels were markedly increased in hypoxia in all cell lines, were also clearly detectable already in normoxia in DMBC28 cells (Figure 2A). While phosphorylation of PDH was constant in DMBC17 and DMBC28 cells, irrespective of oxygen concentration, DMBC12 cells exerted high PDH phosphorylation only in hypoxia, which indicates that PDH was active only in DMBC12 cells in normoxia and hyperoxia (Figure 2B). Activity of the PI3K/AKT pathway varied between cell lines in hyperoxia and was differently altered when cells were exposed to lower oxygen concentrations (Figure 2B). While activity of AKT was gradually increased in normoxia and hypoxia in phospho-AKTlow DMBC28 cells, in another phospho-AKTlow DMBC17 cell line, no increase was observed at lower oxygen concentrations (Figure 2B). In DMBC12 cells, a high level of phospho-AKT was detected in all oxygen conditions.

2.3. Oxygen Concentration-Dependent Changes in the Composition of Melanoma Cell Populations

In the following experiments, the percentages of nerve growth factor receptor (NGFR)- and MITF-positive cells were compared between cell populations grown in different oxygen concentrations (Figure 2C,D). In DMBC12 cell population, NGFR was expressed by 15.2 ± 1.5% cells in hyperoxia, and this percentage was significantly but only slightly higher in normoxia. In NGFRlow DMBC17 cell population (1.9 ± 0.4% in hyperoxia) it was significantly higher in both normoxia after 48 h and hypoxia already after 24 h. DMBC28 cell line, with 20.6 ± 4.3% NGFR-positive cells in hyperoxia, was exceptional as lowering concentration of oxygen to 6% significantly reduced the percentages of NGFR-positive cells after 48 h. Percentages of MITF-positive cells in MITFhigh cell lines were either significantly lower in normoxia and hypoxia than in hyperoxia (DMBC28) or remained unchanged (DMBC17). This suggests that melanoma cells cultured in vitro in the presence of 21% O2 may differ in their phenotypes from melanoma cells grown in vivo at much lower oxygen concentrations.

2.4. Normoxia Promotes the Expression of Glucose Metabolism/Transport-Related Genes and to the Lower Extent Genes Associated with Glutamine Metabolism and Transport

The expression of pivotal glucose and glutamine metabolism/transport-related genes was assessed in melanoma cells exposed to 6% O2 and 1% O2. As the reference, the expression of these genes in 21% O2 was used. We analyzed the expression of genes encoding glucose transporter 1 (GLUT1), hexokinase 2 (HK2), the first enzyme of the glycolytic pathway, and pyruvate dehydrogenase kinase 1 (PDK1), a metabolic gatekeeper, which inhibits the activity of PDH and restrains pyruvate entry to the TCA cycle. All these genes are direct targets of HIF-1α. Accordingly, the expression of all three genes was significantly enhanced when cells were exposed to hypoxia for 24 h (Figure 3A).
PDK1 transcript levels were significantly increased also in normoxia, and this enhancement was especially high in DMBC28 cells. Normoxia induced a significant increase of GLUT1 mRNA levels in DMBC12 cells and DMBC28 cells, whereas the expression of HK2 was significantly increased only in DMBC12 cells. These results show that investigated melanoma cell lines vary in their reaction to a transition from hyperoxia to normoxia. We hypothesized that metabolic adaptation might be required, and indeed when melanoma cells were cultured in normoxic conditions for 3 weeks, the transcript levels of GLUT1, PDK1 and HK2 were markedly enhanced also in DMBC17 cells (Figure 3B). Although significantly changed, the expression of HK2 was the least responsive to hypoxia in all cell lines and proportionally, changes in HK2 expression were the least pronounced also in normoxia. In the following experiments, melanoma cells cultured in normoxia for at least 3 weeks were used.
First however, we investigated whether the expression of genes involved in glutamine uptake and metabolism is subjected to oxygen-dependent regulation after short and long exposure to low oxygen levels (Figure 3C,D). We assessed the transcript levels of solute carrier family 1 member 5 (SLC1A5, also known as ASCT2), which is the main transporter of glutamine into the cell, glutaminase (GLS), which catalyzes the first step of glutamine metabolism by converting glutamine to glutamate, and solute carrier family 7 member 11 transporter (SLC7A11), which exports glutamate in exchange for cystine. When the expression of glutamine-related genes was compared in normoxic and hyperoxic conditions, no significant differences were observed with the exception of SLC7A11 transcript in DMBC28 cells shortly exposed to 6% O2 (Figure 3C), and GLS transcript in DMBC28 cells exposed to normoxia for 3 weeks (Figure 3D).

2.5. Increased Protein Stability Might be Responsible for High PGC1α Level in MITFlow Melanoma Cells

Following the intriguing results of the analysis of PDH phosphorylation, which suggest that only DMBC12 cells might rely on mitochondrial metabolism in normoxic and hyperoxic conditions, we investigated the influence of oxygen level on the expression of PGC1α responsible for immediate cell adaptation for high energetic demands, and MITF, one of the key regulators of PGC1α expression in melanoma cells. Analysis of transcript levels of MITF and PGC1α showed a good correlation between them, as relative to DMBC17 cell line with the highest levels of MITF and PGC1α mRNAs, both transcripts were at markedly lower levels in DMBC28 cells and were almost undetectable in DMBC12 cells (Figure 4A).
Surprisingly, although MITF protein was almost undetectable in DMBC12 cells, a high level of PGC1α protein was found in these cells (Figure 4B). This altogether indicates that while there is a good correlation between transcript levels of MITF and PGC1α in melanoma cells, no correlation could be found at the protein level. To address this discrepancy, we carried out a protein decay assay. We found that while PGC1α half-life did not exceed three hours in DMBC17 and DMBC28 cells, PGC1α protein remained exceptionally stable in DMBC12 cells (Figure 4C,D).

2.6. The Response to Targeted Therapeutics Depends on Oxygen Availability and Original Phenotype of Melanoma Cells

Next, we investigated whether melanoma cell response to vemurafenib, an inhibitor of BRAFV600, and trametinib, an inhibitor of MEK1/2, depends on oxygen concentration. For that we assessed phosphorylation of ERK1/2, an effector protein of the RAS/RAF/MEK/ERK pathway, and MITF, a transcription factor determining the phenotypic state of melanoma cells. While the activity of ERK1/2 was inhibited by vemurafenib and trametinib to an almost undetectable level in all tested cell lines both in normoxia and hyperoxia, MITF level was reduced more substantially in normoxia than in hyperoxia in both MITFhigh cell lines, DMBC17 and DMBC28 (Figure 5A).
Having observed that drug-induced changes in MITF expression in MITFhigh cells depend on oxygen concentration when measured at the bulk protein level (Figure 5A), we asked the question whether MITF expression is downregulated uniformly in all cells or a subpopulation of MITF-positive cells is reduced. Flow cytometric analysis of MITF expression revealed that although the percentages of MITF-positive cells were not substantially changed by drug treatment in hyperoxia, they were significantly decreased in DMBC28 cell population cultured in normoxia (Figure 5B, Supplementary Materials Figure S1). Of note, the percentages of MITF-positive DMBC28 cells differed significantly even in control populations when cultured in hypoxic, normoxic and hyperoxic conditions (Figure 5B). The percentages of MITF-positive cells in DMBC17 cell line were similar in all oxygen concentrations, and remained unaltered in trametinib-treated cells (Figure 5B). MITF-positive cells were almost undetectable in DMBC12 cell line irrespective of oxygen concentration, also after drug treatment (Figure S1).
The percentages of Ki-67-positive cells were similar in control cells of all cell lines and regardless of oxygen concentration, except for DMBC28 cells grown in hypoxia, which were significantly lower than in hyperoxia (Figure 5C). Vemurafenib and trametinib reduced the frequency of Ki-67-positive cells to a variable extent in investigated cell lines after 48 h of treatment. Interestingly, drugs reduced the percentage of Ki-67-positive cells to significantly lower level in hyperoxia than in normoxia, which might indicate that melanoma cells in normoxia could be less responsive to drug treatment. The majority of DMBC28 cells that survived drug treatment in hyperoxia were MITF-positive and Ki-67-negative, in contrast to drug-treated melanoma cells in normoxia, which were predominantly MITF-negative and Ki-67-positive (Figure 5D).

2.7. Drug-Triggered Reduction of PGC1α Level is Accompanied with Induction of Apoptosis

PGC1α has been shown to regulate the survival of PGC-1α-positive melanoma cells, and knockdown of PGC1α significantly induced apoptosis [33]. Our results support and extend these findings as a tendency of drug-triggered PGC1α downregulation in samples with greater PARP cleavage was observed, which might indicate that the stability of this protein can be connected with induction of apoptosis, especially in melanoma cells with lower expression of PGC1α, DMBC28 and DMBC17 (Figure 6).

2.8. Vemurafenib and Trametinib Reduce the Expression of VEGF in Hyperoxia but not in Normoxia

Although melanoma cell lines used in this study did not differ in the basal expression level of VEGFA in hyperoxia (Figure 7A), lowering the oxygen level was associated with significant VEGF transcript level enhancement that varied markedly between cell lines with DMBC12 showing the lowest increase (Figure 7B). Interestingly, they were significantly reduced by vemurafenib and trametinib, however, only under hyperoxic conditions (Figure 7C).

2.9. Vemurafenib and Trametinib Reduce the Expression of Genes Related to Glutamine and Glucose Metabolism and Transport to a Different Extent and in Cell Line-Dependent Manner

We analyzed changes in the expression of glutamine and glucose metabolism-related genes in response to trametinib and vemurafenib in melanoma cells grown in hyperoxia and normoxia for at least 3 weeks (Figure 8). Transcript levels of SLC1A5 were uniformly and significantly decreased in all melanoma cell lines in normoxia and hyperoxia, which suggests that glutamine uptake is downregulated by vemurafenib and trametinib irrespective of the oxygen availability. SLC7A11 expression was downregulated in both oxygen concentrations in two out of three drug-treated cell lines, whereas GLS expression was markedly and significantly decreased only in DMBC28 cells, which suggests that glutamine export and metabolism might be also downregulated by targeted therapeutics, however, not in all melanomas. The effects of vemurafenib and trametinib on the expression of glucose metabolism-related genes were less pronounced. It is worth to note that effects of drugs on transcript level of PDK1 in DMBC12 cells and SLC7A11 in DMBC12 cells and DMBC28 cells were significantly weaker in normoxia than in hyperoxia.

3. Discussion

In vitro cell cultures are an inherent part of cancer research. Much of our knowledge on the molecular mechanisms crucial for cancer maintenance and response to treatment derives from in vitro studies performed in the atmospheric oxygen concentration (21%), which is much higher than in any human tissue. Therefore, despite the undeniable value of cell line-based research, the model is rather flawed and its direct contribution to clinical cancer research is questionable [34]. The main aim of this study was to compare the phenotype of melanoma cells grown in different oxygen concentrations: 21% designated as hyperoxia, 6% chosen as mimicking the physiological oxygen level present in many tissues and regions of tumors and called in this study normoxia, and 1% as hypoxia. In addition, we investigated the oxygen concentration-dependent melanoma cell response to targeted therapeutics, vemurafenib and trametinib.
The patient-derived melanoma cell lines chosen for this study differ in oncogenic mutations in the MAPK/ERK pathway and activity of transcription factor MITF, one of the critical regulators of melanoma phenotype [30,31]. As the oxygen concentration is critical for cancer cell energy metabolism and both, the MAPK/ERK pathway and MITF signaling, have been associated with the control of metabolic activities of melanoma cells [27], the selection of DMBC12 cells (BRAFV600E/MITFlow), DMBC28 cells (BRAFV600E/MITFhigh) and DMBC17 (HRASQ61R/MITFhigh) [30,31] for this study was well justified. And indeed, we were able to demonstrate: (i) cell type-dependent phenotypic changes in response to different oxygen concentrations; (ii) oxygen level-dependent effects of targeted therapeutics in melanoma cells.
It has been demonstrated that melanoma development and metastasis are associated with changes in cell metabolism involving a shift from oxidative phosphorylation (OXPHOS) to aerobic glycolysis [35]. The enhanced glycolytic activity was connected to BRAFV600E, constitutively active in majority of melanomas [36]. BRAF-dependent metabolic rewiring includes activation of HIF-1α [37], which is normally induced in response to low oxygen levels and leads to increased expression of glucose transporters and glycolytic enzymes. Many melanomas display the activity of HIF-1α also under non-hypoxic conditions [38,39]. In our study, HIF-1α was expressed in BRAFV600E melanoma cells in hypoxia and to a lesser extent in normoxia, whereas in HRASQ61R melanoma cells oxygen concentration had to be reduced to 1% to stimulate HIF-1α expression. We showed that the expression of glucose transporter/enzyme-related genes was also oxygen-dependent in both BRAF mutant and HRAS mutant melanoma cells, however, while their expression was significantly higher in BRAF mutant DMBC28 cells already after short exposure to normoxia, HRAS mutant DMBC17 cells needed a longer adaptation period to promote the expression of these genes. In addition to maintaining a state of aerobic glycolysis by enhancing the expression of glucose transporters and glycolytic enzymes, HIF-1α participates in controlling PDH enzymatic activity [40]. PDH is a key enzyme of a multienzyme complex catalyzing pyruvate conversion into acetyl-CoA and linking glycolysis to the tricarboxylic acid (TCA) cycle, thus promoting OXPHOS [41]. Under hypoxic conditions, phosphorylation that inactivates PDH follows an upregulation of its kinase, PDK1, which is HIF-1-dependent [42,43]. However, PDH phosphorylation can be also HIF-1-independent, as it can precede the upregulation of PDK1 mediated by HIF-1α [44]. In this study, protein levels of HIF-1α and phosphorylated PDH suggest that DMBC12 cells display an OXPHOS phenotype since similar and low level of HIF-1α could be detected in normoxia and hypoxia, and a strong inactivating PDH phosphorylation occurred in these cells only in hypoxia, whereas DMBC28 and DMBC17 cell lines might be defined as having a glycolytic phenotype with a high level of inactive PDH irrespective of oxygen availability, especially DMBC28 cell line, in which HIF-1α could be detected also in normoxic conditions. An OXPHOShigh subset of melanomas shows a high expression rate of PGC1α [33], which is an important transcriptional coactivator regulating oxidative metabolism in many tissues [45,46,47,48,49]. The role of PGC1α is controversial as it can act as a tumor promoter and a tumor suppressor, which can be partially explained by its cell type specificity, including expression and diverse interacting proteins [49,50,51,52]. Multiple promoters of PPARGC1A, a gene encoding PGC1α, are selectively responsive in distinct tissues, and MITF was found to directly stimulate the expression of PGC1α by binding to its promoter region in melanoma cells [27]. MITF is a master regulator of melanocytes and a lineage-specific melanoma oncogene [53], and it is important for cell differentiation, proliferation and survival [54,55]. MITF activity is tightly regulated in melanocytes and melanoma [56,57], including transcriptional repression by HIF-1α [58]. A correlation between transcript levels of MITF and PGC1α was previously demonstrated in melanoma cells [27], and this correlation was also confirmed in our study. Importantly, for the first time we have shown that while MITFlow/no melanoma cells display a low level of PGC1α transcript, they may accumulate PGC1α protein reaching high level, probably due to mechanism(s) that increase its stability. According to literature, PGC1α protein has an approximate half-life of 2–3 h and multiple factors are responsible for its stability [59]. The mechanism(s) preserving a high level of PGC1α protein in melanoma cells for longer than 10 h needs to be elucidated.
In addition to altering cancer cell metabolism, PGC1α has been implicated in cell survival, and knockdown of PGC1α results in induction of apoptosis in melanoma cells [33]. Our results are in agreement with this notion as we have shown that drug-triggered reduction of PGC1α level is accompanied with induction of apoptosis. Therefore, the paradigm that PGC1α is increased in response to drug-triggered inhibition of ERK1/2 activity [27] might be true unless the inhibition of ERK1/2 is accompanied with induction of apoptosis. Interestingly, in PGC1αhigh DMBC12 cell line both drug-induced reduction of PGC1α and induction of apoptosis were less pronounced than in PGC1αlow cell lines, DMBC28 and DMBC17.
Cell line-specific effects of adjusting oxygen level to this more closely resembling tumor natural one can be also exemplified by variable increase in AKT activity, expression of VEGF and MITF, and cell response to targeted therapeutics in terms of percentages of MITF- and Ki-67-positive cells at the cell population level. Our finding that the activity of the PI3K/AKT pathway is not universally activated in all patient-derived melanoma cell lines exposed to decreased oxygen level, has also been reported in other cancer types [60]. Notably, neither of the cell lines used in our study harbored mutations in PTEN [30], therefore, the mechanism behind high AKT activity in DMBC12 is likely not of genetic background. It has been recently shown that AMPK, a major energy sensor, phosphorylates mTORC2, which leads to phosphorylation of AKT at Ser473 [61]. This could partially explain the increase in phospho-AKT in DMBC28 in hypoxia and high levels of it in the fast proliferating DMBC12 cell line, however, it would not explain lack of AKT activity in DMBC17. Moreover, we observed an overlapping activity of AKT and expression of PGC1α, which is also activated by AMPK [62].
The effect of targeted therapeutics on the expression of glutaminase and glutamine transporters in drug-sensitive, patient-derived melanoma cells have not been described before. However, several reports indicate the role of mitochondrial metabolism of glutamine in drug-resistance [63,64]. In those studies, melanoma cells with acquired resistance to BRAFV600E or MEK inhibitors depended on glutamine metabolism for survival and were susceptible to glutaminase inhibition. The relevance of glutamine transport in melanoma is supported by data indicating that SLC1A5 and SLC7A11 are overexpressed in melanoma cells when compared to melanocytes [65,66]. Whether the expression of these genes would be restored in resistant cells remains to be clarified. The selective increase in SLC7A11 expression in hypoxic DMBC28 cells is unclear, however, such increase has been previously described in breast cancer [67] and renal cancer cell lines [68], due to high HIF activity or glucose starvation, respectively.
The MITF/PGC1a axis has also been implicated in the melanoma cell response to targeted therapy. Some studies have shown that selective inhibition of BRAFV600E leads to increased MITF expression followed by upregulated PGC1α levels that facilitate drug resistance [17,27,69]. On the other hand, other studies, including our own, have shown that melanoma cells may lose MITF expression during acquisition of drug resistance [31,70,71]. In this context, both upregulation and downregulation of MITF expression may confer an advantage for melanoma cells treated with targeted therapeutics. Whether in this setting MITF downregulation in physiological normoxia provided direct advantage is unclear, however, it was concomitant with limited effects of drugs on the frequency of Ki-67-positive cells in normoxia (6% O2) when compared with drug effects in the atmospheric oxygen concentration (21%). Clinically, high proliferative index of recurrent melanomas, displayed as high Ki-67 expression, is an independent predictor of worse overall survival [72,73]. Our study indicates that the response of melanoma cells to treatment can be greatly altered by artificial levels of oxygen in cell culture. Downregulation of VEGF expression by vemurafenib and trametinib in hyperoxia falls into the same category; it is likely an artefact of 21% O2 in vitro culture and would presumably be disproved in a mouse model. It serves, however, as a proof of concept that an excessive oxygen concentration in cell culture may contribute data that cannot be reproduced in vivo. It is reasonable to assume that an opposite effect may also be observed where a tested compound is less effective in 21% O2 than it would be in physiologically relevant oxygen levels and becomes screened out before being investigated in vivo.
Concentration of oxygen plays an important role in maintenance of stemness, in normal tissue but also in cancer [74,75]. While in some studies, low oxygen levels promote stem like properties shown as increased NGFR level, in other studies low oxygen concentrations predominantly promote invasive phenotype with downregulated NGFR expression [76,77,78]. Both stem-like and invasive properties are connected with increased glycolysis and hypoxia [26,77,79]. All differences in the composition of subpopulations, namely NGFR-positive, MITF-positive and Ki-67-positive cells observed in this study between normoxia and hyperoxia, were cell-line specific but also point at oxygen as the causative factor of melanoma heterogeneity.
This is a preliminary study limited to a small sample number. Nevertheless, within just three patient-derived cell lines we observed distinct responses to normoxia, which demonstrate the diversity and plasticity of melanoma cells. A large scale comparison of metabolic phenotypes in varying oxygen concentrations would be required. We have shown that patient-derived melanoma cells in normoxia is a valid research model and that although additional adjustments in cell culture have to be made, melanoma cells can be successfully propagated in physiologically relevant oxygen levels in a regular manner. The application of normoxia research model could potentially reach beyond in vitro experiments as well. A comparative study of patient-derived xenografts (PDX) and cell line-derived xenografts (CDX) found transcriptome-wide differences between the two xenograft models that were attributed to divergent regulation of hypoxic response. Bhadury and colleagues concluded that cell lines became ‘pseudo-hypoxic’ when injected into mice, owing it to prior adaptation of these cells to standard laboratory culture conditions [80]. Arguably, tumor cells pre-adapted in vitro to an oxygen concentration relevant to in vivo settings might also prove advantageous in animal research models. Thus, this study, although preliminary, produced encouraging data that validates normoxia (6% O2) as a research model in vitro.

4. Materials and Methods

4.1. Tumor Tissues and Cell Culture

Patient-derived melanoma cells were obtained during surgical interventions [31,81]. The study was approved by Ethical Commission of Medical University of Lodz (identification code: RNN/84/09/KE). Each patient signed an informed consent before tissue acquisition. Tumor tissues were processed immediately after surgical procurement and suspensions of melanoma cells for culturing were generated within 2 h. Cells were named DMBC12, DMBC17, and DMBC28, (for Department of Molecular Biology of Cancer, DMBC). Melanoma cells were maintained in stem cell medium (SCM), consisting of DMEM/F12 medium, B-27 supplement (Gibco, Paisley, UK), insulin (10 μg/mL), heparin (1 ng/mL), 10 ng/mL bFGF, 20 ng/mL EGF (BD Biosciences, San Jose, CA, USA) and antibiotics (100 IU/mL penicillin, 100 μg/mL streptomycin). For initial part of the study involving proliferation assays, optimization of cell density, preliminary investigation of oxygen-dependent changes in gene expression as well as all experiments in hypoxia, cells were exposed to indicated oxygen concentration for the length of the assay. For all further analysis on physiological normoxia, melanoma cells were maintained in 6% for at least 3 weeks prior to experiments and all throughout the study. For experiments, cells were left to grow overnight before being treated with trametinib at the concentrations of 30 nM or 50 nM and vemurafenib at 5 μM or 10 μM. Then, cells were collected for RNA isolation (after 24 h), cell lysates (after 26 h or 44 h), viability analysis (after 72 h) and immunophenotype analysis (after 24 h or 48 h).

4.2. A Time-Lapse Fluorescence Microscopy

For proliferation analysis, a time-lapse fluorescence microscope system IncuCyte ZOOM (Essen Bioscience, Ann Arbor, MI, USA) was used. The data were analyzed using the IncuCyte ZOOM software. Proliferation was determined as the image area occupied by melanoma cells. Changes in the confluence are presented as fold-change relative to time T0.

4.3. Flow Cytometry

Flow cytometric data were acquired with FACSVerse (BD Biosciences, San Jose, CA, USA), and analyzed using BD FACSuite. For viability assay, cells were stained with Annexin V-FLUOS Staining Kit (Roche, Mannheim, Germany) for 15 min. For immunophenotype analysis, cells were left to grow overnight before being treated with drugs for 24 h or 48 h. LIVE/DEAD Fixable Violet Dead Cell Stain Kit (Life Technologies, Eugene, OR, USA) was used to exclude dead cells from analysis. For intracellular staining, cells were fixed with 4% paraformaldehyde and permeabilized with 0.1% Triton X-100 in PBS for 20 min. Antibodies against Ki-67 (Alexa Fluor647-conjugated) and NGFR (PE-conjugated) were from BD Biosciences, and antibodies against MITF (Alexa Fluor488-conjugated) were from Abcam (Cambridge, UK). Appropriate isotype controls were included in each experiment.

4.4. RNA Isolation, cDNA Synthesis, and Quantitative RT-PCR (qRT-PCR)

Total RNA was isolated and purified using Total RNA Isolation kit with mini column system (A&A Biotechnology, Gdynia, Poland). Total RNA (1 μg) was transcribed into cDNA using 300 ng of random primers and SuperScript II Reverse Transcriptase (Invitrogen Life Technologies, Carlsbad, CA, USA). The evaluation of mRNA expression of selected genes was performed by quantitative real-time polymerase chain reaction using the Rotor-Gene 3000 Real-Time DNA analysis system (Corbett Research, Mortlake, Australia). Amplification was performed using KAPA SYBR FAST qPCR Kit Universal 2X qPCR Master Mix (Kapa Biosystems, Cape Town, South Africa), 200 nM of each primer and 25 ng cDNA template per reaction. Primer sequences were as follows: GLS forward 5′-TGC ATT CCT GTG GCA TGT AT-3′ and reverse 5′-TTG CCC ATC TTA TCC AGA GG-3′, GLUT1 forward 5′-CGG GCC AAG AGT GTG CTA AA-3′ and reverse 5′-TGA CGA TAC CGG AGC CAA TG-3′, HK2 forward 5′-TAG GGC TTG AGA GCA CCT GT-3′ and reverse 5′-CCA CAC CCA CTG TCA CTT TG-3′, MITF-M forward 5′-ACC GTC TCT CAC TGG ATT GG-3′ and reverse 5′-TAC TTG GTG GGG TTT TCG AG-3′, PDK1 forward 5′-CTA TGA AAA TGC TAG GCG TCT-3′ and reverse 5′-AAC CAC TTG TAT TGG CTG TCC-3′, PGC1α forward 5′-AAC AGC AGC AGA GAC AAA TGC ACC-3′ and reverse 5′-TGC AGT TCC AGA GAG TTC CAC ACT-3′, RPS17 forward 5′-AAT CTC CTG ATC CAA GGC TG-3′ and reverse 5′-CAA GAT AGC AGG TTA TGT CAC G-3′, SLC1A5 forward 5′-GAG CTG CTT ATC CGC TTC TTC-3′ and reverse 5′-GGG GCG TAC CAC ATG ATC C-3′, SLC7A11 forward 5′-TGC TGG GCT GAT TTA TCT TCG-3′ and reverse 5′-GAA AGG GCA ACC ATG AAG AGG-3′, VEGF forward 5′-TTG CCT TGC TGC TCT ACC TCC A-3′ and reverse 5′-GAT GGC AGT AGC TGC GCT GAT A-3′. To calculate the relative expression of target genes vs. a reference gene RPS17, a mathematical model including an efficiency correction was used.

4.5. Western Blot

Melanoma cells were lysed in RIPA buffer containing 50 mmol/l Tris-HCl pH 8.0, 150 mmol/l NaCl, 1% Triton X-100, 0.5% sodium deoxycholate, 0.1% SDS supplemented with freshly added protease and phosphatase inhibitors (Merck, Darmstadt, Germany). Protein concentration was determined by Bradford assay. Cell lysates were diluted in 2 × Laemmli buffer (125 mM Tris-HCl pH 6.8, 0.004% bromophenol blue, 20% glycerol, 4% SDS and 10% β-mercaptoethanol) and protein samples were loaded on 7% SDS-polyacrylamide gel followed by electrophoresis at constant voltage 25 V/cm. The proteins were transferred onto an Immobilon-P PVDF membrane (Millipore, Billerica, MA, USA). The membrane was incubated in a blocking solution: 5% nonfat milk in PBS-Tween 0.05% or 5% phospho- BLOCKER (Cell Biolabs, San Diego, CA, USA) in PBS-Tween 0.05% for 45 min. Primary antibodies detecting ERK1/2, phospho-ERK1/2 (Thr202/Tyr204), AKT, phospho-AKT (Ser473), HIF-2α, MITF, PGC1α, MCL-1 and PARP were from Cell Signaling Technology (Danvers, MA, USA), PDH and GAPDH from Santa Cruz Biotechnology (Santa Cruz, CA, USA), phospho-PDH (Ser293) and β-actin from Merck (Darmstadt, Germany), and HIF-1α from BD Biosciences. Secondary HRP-conjugated anti-mouse antibodies were from Santa Cruz Biotechnology and anti-rabbit antibodies were from Cell Signaling. The membrane was incubated with Pierce® ECL Western Blotting Substrate (Pierce, Rockford, IL, USA) for 1 min and the proteins were visualized using ChemiDoc Imaging System (BioRad, Hercules, CA, USA).

4.6. Protein Decay Assay

Melanoma cells were incubated with 100 µg/mL cycloheximide (Sigma-Aldrich, Saint Louis, MO, USA) to inhibit de novo protein synthesis. Cells were collected at 1 h, 2 h, 3 h, 4 h, 6 h and 10 h intervals after treatment with cycloheximide for cell lysate preparation and Western blotting.

4.7. Drugs

Vemurafenib and trametinib were purchased from Selleck Chemicals LLC (Houston, TX, USA).

4.8. Statistical Analysis

Graphs are presented as mean ± SD. Student’s t-test was performed to determine significant differences between the mean values of the tested parameters. The difference was considered significant if p < 0.05.

Supplementary Materials

Supplementary Materials can be found at https://www.mdpi.com/1422-0067/20/17/4203/s1.

Author Contributions

Conceptualization, M.C.; investigation, M.O. and M.H.; data curation, M.O. and M.C.; writing—original draft preparation, M.O.; writing—review and editing, M.C.; supervision, M.C.; funding acquisition, M.C.

Funding

This research was funded by the Grant 2014/15/B/NZ7/00947 from National Science Center (Poland).

Acknowledgments

We thank Ewa Lewandowska and Karolina Niewinna for excellent technical support.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

AKTProtein Kinase B
AMPKAMP-Activated Protein Kinase
AXLTyrosine-protein kinase receptor
bFGFBasic fibroblast growth factor
CDXCell line-derived xenografts
CHXCycloheximide
DMBCDepartment of Molecular Biology of Cancer
EGFEpidermal growth factor
GAPDHGlyceraldehyde 3-phosphate dehydrogenase
GLSGlutaminase
GLUT1Glucose Transporter Type 1
HIF-1αHypoxia Inducible Factor 1 Subunit Alpha
HIF-2αHypoxia Inducible Factor 2 Subunit Alpha
HK2Hexokinase 2
MAPKMitogen-Activated Protein Kinase
MCL-1Myeloid Cell Leukemia 1
MITFMicrophthalmia-Associated Transcription Factor
mTORC2Mammalian target of rapamycin complex 2
NGFRNerve Growth Factor Receptor
OXPHOSOxidative phosphorylation
PDHPyruvate Dehydrogenase
PDK1Pyruvate Dehydrogenase Kinase 1
PDXPatient-derived xenografts
PGC1αPeroxisome Proliferator-Activated Receptor Gamma Coactivator 1-Alpha
PI3KPhosphatidylinositol-4,5-Bisphosphate 3-Kinase
RPS17Ribosomal Protein S17
SLC1A5Solute Carrier Family 1 Member 5
SLC7A11Solute Carrier Family 7 Member 11
SOX10Sex Determining Region Y-Box 10
TCAtricarboxylic acid cycle
TMETumor microenvironment
TRATrametinib
VEGFVascular Endothelial Growth Factor
VEMVemurafenib

References

  1. Smalley, K.S.; Lioni, M.; Noma, K.; Haass, N.K.; Herlyn, M. In vitro three-dimensional tumor microenvironment models for anticancer drug discovery. Expert Opin. Drug Discov. 2008, 3, 1–10. [Google Scholar] [CrossRef] [PubMed]
  2. Vörsmann, H.; Groeber, F.; Walles, H.; Busch, S.; Beissert, S.; Walczak, H.; Kulms, D. Development of a human three-dimensional organotypic skin-melanoma spheroid model for in vitro drug testing. Cell Death Dis. 2013, 4, e719. [Google Scholar] [CrossRef] [PubMed]
  3. Flach, E.H.; Rebecca, V.W.; Herlyn, M.; Smalley, K.S.; Anderson, A.R. Fibroblasts contribute to melanoma tumour growth and drug resistance. Mol. Pharm. 2011, 8, 2039–2049. [Google Scholar] [CrossRef] [PubMed]
  4. Somasundaram, R.; Zhang, G.; Fukunaga-Kalabis, M.; Perego, M.; Krepler, C.; Xu, X.; Wagner, C.; Hristova, D.; Zhang, J.; Tian, T.; et al. Tumor-associated B-cells induce tumor heterogeneity and therapy resistance. Nat. Commun. 2017, 8, 607. [Google Scholar] [CrossRef] [PubMed]
  5. Sztiller-Sikorska, M.; Hartman, M.L.; Talar, B.; Jakubowska, J.; Zalesna, I.; Czyz, M. Phenotypic diversity of patient-derived melanoma populations in stem cell medium. Lab. Invest. 2015, 95, 672–683. [Google Scholar] [CrossRef] [Green Version]
  6. Zalesna, I.; Osrodek, M.; Hartman, M.L.; Rozanski, M.; Sztiller-Sikorska, M.; Niewinna, K.; Nejc, D.; Czyz, M. Exogenous growth factors bFGF, EGF and HGF do not influence viability and phenotype of V600EBRAF melanoma cells and their response to vemurafenib and trametinib in vitro. PLoS ONE 2017, 12, e0183498. [Google Scholar] [CrossRef] [PubMed]
  7. Lee, J.; Kotliarova, S.; Kotliarov, Y.; Li, A.; Su, Q.; Donin, N.M.; Pastorino, S.; Purow, B.W.; Christopher, N.; Zhang, W.; et al. Tumor stem cells derived from glioblastomas cultured in bFGF and EGF more closely mirror the phenotype and genotype of primary tumors than do serum-cultured cell lines. Cancer Cell 2006, 9, 391–403. [Google Scholar] [CrossRef] [Green Version]
  8. Hartman, M.L.; Talar, B.; Noman, M.Z.; Gajos-Michniewicz, A.; Chouaib, S.; Czyz, M. Gene expression profiling identifies microphthalmia-associated transcription factor (MITF) and Dickkopf-1 (DKK1) as regulators of microenvironment-driven alterations in melanoma phenotype. PLoS ONE 2014, 9, e95157. [Google Scholar] [CrossRef]
  9. McKeown, S.R. Defining normoxia, physoxia and hypoxia in tumours-implications for treatment response. Br. J. Radiol. 2014, 87, 20130676. [Google Scholar] [CrossRef]
  10. Carreau, A.; El Hafny-Rahbi, B.; Matejuk, A.; Grillon, C.; Kieda, C. Why is the partial oxygen pressure of human tissues a crucial parameter? Small molecules and hypoxia. J. Cell Mol. Med. 2011, 15, 1239–1253. [Google Scholar] [CrossRef] [Green Version]
  11. Forsyth, N.R.; Musio, A.; Vezzoni, P.; Simpson, A.H.; Noble, B.S.; McWhir, J. Physiologic oxygen enhances human embryonic stem cell clonal recovery and reduces chromosomal abnormalities. Cloning Stem Cells 2006, 8, 16–23. [Google Scholar] [CrossRef] [PubMed]
  12. Atkuri, K.R.; Herzenberg, L.A.; Niemi, A.K.; Cowan, T.; Herzenberg, L.A. Importance of culturing primary lymphocytes at physiological oxygen levels. Proc. Natl. Acad. Sci. USA 2007, 104, 4547–4552. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Guo, C.W.; Kawakatsu, M.; Idemitsu, M.; Urata, Y.; Goto, S.; Ono, Y.; Hamano, K.; Li, T.S. Culture under low physiological oxygen conditions improves the stemness and quality of induced pluripotent stem cells. J. Cell Physiol. 2013, 228, 2159–2166. [Google Scholar] [CrossRef] [PubMed]
  14. Mas-Bargues, C.; Sanz-Ros, J.; Román-Domínguez, A.; Inglés, M.; Gimeno-Mallench, L.; El Alami, M.; Viña-Almunia, J.; Gambini, J.; Viña, J.; Borrás, C. Relevance of Oxygen Concentration in Stem Cell Culture for Regenerative Medicine. Int. J. Mol. Sci. 2019, 20, E1195. [Google Scholar] [CrossRef] [PubMed]
  15. Rezvani, H.R.; Ali, N.; Nissen, L.J.; Harfouche, G.; de Verneuil, H.; Taïeb, A.; Mazurier, F. HIF-1α in epidermis: Oxygen sensing, cutaneous angiogenesis, cancer, and non-cancer disorders. J. Invest. Dermatol. 2011, 131, 1793–1805. [Google Scholar] [CrossRef] [PubMed]
  16. Savoia, P.; Fava, P.; Casoni, F.; Cremona, O. Targeting the ERK Signaling Pathway in Melanoma. Int. J. Mol. Sci. 2019, 20, E1483. [Google Scholar] [CrossRef] [PubMed]
  17. Gopal, Y.N.; Rizos, H.; Chen, G.; Deng, W.; Frederick, D.T.; Cooper, Z.A.; Scolyer, R.A.; Pupo, G.; Komurov, K.; Sehgal, V.; et al. Inhibition of mTORC1/2 overcomes resistance to MAPK pathway inhibitors mediated by PGC1alpha and oxidative phosphorylation in melanoma. Cancer Res. 2014, 74, 7037–7047. [Google Scholar] [CrossRef]
  18. Amaral, T.; Sinnberg, T.; Meier, F.; Krepler, C.; Levesque, M.; Niessner, H.; Garbe, C. MAPK pathway in melanoma part II-secondary and adaptive resistance mechanisms to BRAF inhibition. Eur. J. Cancer 2017, 73, 93–101. [Google Scholar] [CrossRef]
  19. Hugo, W.; Shi, H.; Sun, L.; Piva, M.; Song, C.; Kong, X.; Moriceau, G.; Hong, A.; Dahlman, K.B.; Johnson, D.B.; et al. Non-genomic and Immune Evolution of Melanoma Acquiring MAPKi Resistance. Cell 2015, 162, 1271–1285. [Google Scholar] [CrossRef] [Green Version]
  20. Arozarena, I.; Wellbrock, C. Phenotype plasticity as enabler of melanoma progression and therapy resistance. Nat. Rev. Cancer 2019, 19, 377–391. [Google Scholar] [CrossRef] [Green Version]
  21. Bedogni, B.; Powell, M.B. Hypoxia, melanocytes and melanoma – survival and tumor development in the permissive microenvironment of the skin. Pigment. Cell Melanoma Res. 2009, 22, 166–174. [Google Scholar] [CrossRef]
  22. Pucciarelli, D.; Lengger, N.; Takáčová, M.; Csaderova, L.; Bartosova, M.; Breiteneder, H.; Pastorekova, S.; Hafner, C. Hypoxia increases the heterogeneity of melanoma cell populations and affects the response to vemurafenib. Mol. Med. Rep. 2016, 13, 3281–3288. [Google Scholar] [CrossRef]
  23. Widmer, D.S.; Hoek, K.S.; Cheng, P.F.; Eichhoff, O.M.; Biedermann, T.; Raaijmakers, M.I.G.; Hemmi, S.; Dummer, R.; Levesque, M.P. Hypoxia contributes to melanoma heterogeneity by triggering HIF1α-dependent phenotype switching. J. Invest. Dermatol. 2013, 133, 2436–2443. [Google Scholar] [CrossRef]
  24. Qin, Y.; Roszik, J.; Chattopadhyay, C.; Hashimoto, Y.; Liu, C.; Cooper, Z.A.; Wargo, J.A.; Hwu, P.; Ekmekcioglu, S.; Grimm, E.A. Hypoxia-Driven Mechanism of Vemurafenib Resistance in Melanoma. Mol. Cancer Ther. 2016, 15, 2442–2454. [Google Scholar] [CrossRef]
  25. Koch, A.; Lang, S.A.; Wild, P.J.; Gantner, S.; Mahli, A.; Spanier, G.; Berneburg, M.; Müller, M.; Bosserhoff, A.K.; Hellerbrand, C. Glucose transporter isoform 1 expression enhances metastasis of malignant melanoma cells. Oncotarget 2015, 6, 32748–32760. [Google Scholar] [CrossRef] [Green Version]
  26. Pinheiro, C.; Miranda-Gonçalves, V.; Longatto-Filho, A.; Vicente, A.L.; Berardinelli, G.N.; Scapulatempo-Neto, C.; Costa, R.F.; Viana, C.R.; Reis, R.M.; Baltazar, F.; et al. The metabolic microenvironment of melanomas: Prognostic value of MCT1 and MCT4. Cell Cycle 2016, 15, 1462–1470. [Google Scholar] [CrossRef] [Green Version]
  27. Haq, R.; Shoag, J.; Andreu-Perez, P.; Yokoyama, S.; Edelman, H.; Rowe, G.C.; Frederick, D.T.; Hurley, A.D.; Nellore, A.; Kung, A.L.; et al. Oncogenic BRAF regulates oxidative metabolism via PGC1α and MITF. Cancer Cell 2013, 23, 302–315. [Google Scholar] [CrossRef]
  28. Anderson, A.S.; Roberts, P.C.; Frisard, M.I.; Hulver, M.W.; Schmelz, E.M. Ovarian tumor-initiating cells display a flexible metabolism. Exp. Cell. Res. 2014, 328, 44–57. [Google Scholar] [CrossRef] [Green Version]
  29. Hartman, M.L.; Rozanski, M.; Osrodek, M.; Zalesna, I.; Czyz, M. Vemurafenib and trametinib reduce expression of CTGF and IL-8 in V600EBRAF melanoma cells. Lab. Invest. 2017, 97, 217–227. [Google Scholar] [CrossRef] [Green Version]
  30. Hartman, M.L.; Sztiller-Sikorska, M.; Czyz, M. Whole-exome sequencing reveals novel genetic variants associated with diverse phenotypes of melanoma cells. Mol. Carcinog. 2019, 58, 588–602. [Google Scholar] [CrossRef]
  31. Czyz, M.; Sztiller-Sikorska, M.; Gajos-Michniewicz, A.; Osrodek, M.; Hartman, M.L. Plasticity of Drug-Naïve and Vemurafenib- or Trametinib-Resistant Melanoma Cells in Execution of Differentiation/Pigmentation Program. J. Oncol. 2019, 2019, 1–15. [Google Scholar] [CrossRef]
  32. Tsoi, J.; Robert, L.; Paraiso, K.; Galvan, C.; Sheu, K.M.; Lay, J.; Wong, D.J.L.; Atefi, M.; Shirazi, R.; Wang, X.; et al. Multi-stage differentiation defines melanoma subtypes with differential vulnerability to drug-induced iron-dependent oxidative stress. Cancer Cell 2018, 33, 890–904. [Google Scholar] [CrossRef]
  33. Vazquez, F.; Lim, J.H.; Chim, H.; Bhalla, K.; Girnun, G.; Pierce, K.; Clish, C.B.; Granter, S.R.; Widlund, H.R.; Spiegelman, B.M.; et al. PGC1alpha expression defines a subset of human melanoma tumors with increased mitochondrial capacity and resistance to oxidative stress. Cancer Cell 2013, 23, 287–301. [Google Scholar] [CrossRef]
  34. Gillet, J.P.; Calcagno, A.M.; Varma, S.; Marino, M.; Green, L.J.; Vora, M.I.; Patel, C.; Orina, J.N.; Eliseeva, T.A.; Singal, V.; et al. Redefining the relevance of established cancer cell lines to the study of mechanisms of clinical anti-cancer drug resistance. Proc. Natl. Acad. Sci. USA 2011, 108, 18708–18713. [Google Scholar] [CrossRef] [Green Version]
  35. Abildgaard, C.; Guldberg, P. Molecular drivers of cellular metabolic reprogramming in melanoma. Trends. Mol. Med. 2015, 21, 164–171. [Google Scholar] [CrossRef]
  36. Hall, A.; Meyle, K.D.; Lange, M.K.; Klima, M.; Sanderhoff, M.; Dahl, C.; Abildgaard, C.; Thorup, K.; Moghimi, S.M.; Jensen, P.B.; et al. Dysfunctional oxidative phosphorylation makes malignant melanoma cells addicted to glycolysis driven by the (V600E)BRAF oncogene. Oncotarget 2013, 4, 584–599. [Google Scholar] [CrossRef]
  37. Kumar, S.M.; Yu, H.; Edwards, R.; Chen, L.; Kazianis, S.; Brafford, P.; Acs, G.; Herlyn, M.; Xu, X. Mutant V600EBRAF increases hypoxia inducible factor-1alpha expression in melanoma. Cancer Res. 2007, 67, 3177–3184. [Google Scholar] [CrossRef]
  38. Kuphal, S.; Winklmeier, A.; Warnecke, C.; Bosserhoff, A.K. Constitutive HIF-1 activity in malignant melanoma. Eur. J. Cancer 2010, 46, 1159–1169. [Google Scholar] [CrossRef]
  39. Mills, C.N.; Joshi, S.S.; Niles, R.M. Expression and function of hypoxia inducible factor-1 alpha in human melanoma under non-hypoxic conditions. Mol. Cancer 2009, 8, 104. [Google Scholar] [CrossRef]
  40. Kluza, J.; Corazao-Rozas, P.; Touil, Y.; Jendoubi, M.; Maire, C.; Guerreschi, P.; Jonneaux, A.; Ballot, C.; Balayssac, S.; Valable, S.; et al. Inactivation of the HIF-1alpha/PDK3 signaling axis drives melanoma toward mitochondrial oxidative metabolism and potentiates the therapeutic activity of pro-oxidants. Cancer Res. 2012, 72, 5035–5047. [Google Scholar] [CrossRef]
  41. Patel, M.S.; Nemeria, N.S.; Furey, W.; Jordan, F. The pyruvate dehydrogenase complexes: Structure-based function and regulation. J. Biol. Chem. 2014, 289, 16615–16623. [Google Scholar] [CrossRef]
  42. Kim, J.W.; Tchernyshyov, I.; Semenza, G.L.; Dang, C.V. HIF-1- mediated expression of pyruvate dehydrogenase kinase: A metabolic switch required for cellular adaptation to hypoxia. Cell Metab. 2006, 3, 177–185. [Google Scholar] [CrossRef]
  43. Semenza, G.L. Hypoxia-inducible factors: Mediators of cancer progression and targets for cancer therapy. Trends Pharmacol. Sci. 2012, 33, 207–214. [Google Scholar] [CrossRef]
  44. Zimmer, A.D.; Walbrecq, G.; Kozar, I.; Behrmann, I.; Haan, C. Phosphorylation of the pyruvate dehydrogenase complex precedes HIF-1-mediated effects and pyruvate dehydrogenase kinase 1 upregulation during the first hours of hypoxic treatment in hepatocellular carcinoma cells. Hypoxia 2016, 4, 135–145. [Google Scholar] [CrossRef]
  45. Martínez-Redondo, V.; Pettersson, A.T.; Ruas, J.L. The hitchhiker’s guide to PGC-1α isoform structure and biological functions. Diabetologia 2015, 58, 1969–1977. [Google Scholar] [CrossRef]
  46. St-Pierre, J.; Lin, J.; Krauss, S.; Tarr, P.T.; Yang, R.; Newgard, C.B.; Spiegelman, B.M. Bioenergetic analysis of peroxisome proliferator-activated receptor gamma coactivators 1alpha and 1beta (PGC-1alpha and PGC-1beta) in muscle cells. J. Biol. Chem. 2003, 278, 26597–26603. [Google Scholar] [CrossRef]
  47. Scarpulla, R.C. Metabolic control of mitochondrial biogenesis through the PGC-1 family regulatory network. Biochim. Biophys. Acta 2011, 1813, 1269–1278. [Google Scholar] [CrossRef] [Green Version]
  48. Lin, J.; Handschin, C.; Spiegelman, B.M. Metabolic control through the PGC-1 family of transcription coactivators. Cell Metab. 2005, 1, 361–370. [Google Scholar] [CrossRef] [Green Version]
  49. Cheng, C.F.; Ku, H.C.; Lin, H. PGC-1α as a Pivotal Factor in Lipid and Metabolic Regulation. Int. J. Mol. Sci. 2018, 19, E3447. [Google Scholar] [CrossRef]
  50. Luo, C.; Lim, J.H.; Lee, Y.; Granter, S.R.; Thomas, A.; Vazquez, F.; Widlund, H.R.; Puigserver, P. A PGC1α-mediated transcriptional axis suppresses melanoma metastasis. Nature 2016, 537, 422–426. [Google Scholar] [CrossRef]
  51. Andrzejewski, S.; Klimcakova, E.; Johnson, R.M.; Tabariès, S.; Annis, M.G.; McGuirk, S.; Northey, J.J.; Chénard, V.; Sriram, U.; Papadopoli, D.J.; et al. PGC-1alpha Promotes Breast Cancer Metastasis and Confers Bioenergetic Flexibility against Metabolic Drugs. Cell Metab. 2017, 26, 778–787. [Google Scholar] [CrossRef]
  52. LeBleu, V.S.; O’Connell, J.T.; Gonzalez Herrera, K.N.; Wikman, H.; Pantel, K.; Haigis, M.C.; de Carvalho, F.M.; Damascena, A.; Domingos Chinen, L.T.; Rocha, R.M.; et al. PGC-1α mediates mitochondrial biogenesis and oxidative phosphorylation in cancer cells to promote metastasis. Nat. Cell Biol. 2014, 16, 992–1003. [Google Scholar] [CrossRef]
  53. Garraway, L.A.; Widlund, H.R.; Rubin, M.A.; Getz, G.; Berger, A.J.; Ramaswamy, S.; Beroukhim, R.; Milner, D.A.; Granter, S.R.; Du, J.; et al. Integrative genomic analyses identify MITF as a lineage survival oncogene amplified in malignant melanoma. Nature 2005, 436, 117–122. [Google Scholar] [CrossRef]
  54. Levy, C.; Khaled, M.; Fisher, D.E. MITF: Master regulator of melanocyte development and melanoma oncogene. Trends Mol. Med. 2006, 12, 406–414. [Google Scholar] [CrossRef]
  55. Hartman, M.L.; Czyz, M. Pro-survival role of MITF in melanoma. J. Invest. Dermatol. 2015, 135, 352–358. [Google Scholar] [CrossRef]
  56. Hartman, M.L.; Czyz, M. MITF in melanoma: Mechanisms behind its expression and activity. Cell Mol. Life Sci. 2015, 72, 1249–1260. [Google Scholar] [CrossRef]
  57. Kawakami, A.; Fisher, D.E. The master role of microphthalmia-associated transcription factor in melanocyte and melanoma biology. Lab. Invest. 2017, 97, 649–656. [Google Scholar] [CrossRef] [Green Version]
  58. Feige, E.; Yokoyama, S.; Levy, C.; Khaled, M.; Igras, V.; Lin, R.J.; Lee, S.; Widlund, H.R.; Granter, S.R.; Kung, A.L.; et al. Hypoxia-induced transcriptional repression of the melanoma-associated oncogene MITF. Proc. Natl. Acad. Sci. USA 2011, 108, E924–E933. [Google Scholar] [CrossRef]
  59. Puigserver, P.; Rhee, J.; Lin, J.; Wu, Z.; Yoon, J.C.; Zhang, C.Y.; Krauss, S.; Mootha, V.K.; Lowell, B.B.; Spiegelman, B.M. Cytokine stimulation of energy expenditure through p38 MAP kinase activation of PPARgamma coactivator-1. Mol. Cell 2001, 8, 971–982. [Google Scholar] [CrossRef]
  60. Stegeman, H.; Span, P.N.; Peeters, W.J.; Verheijen, M.M.; Grénman, R.; Meijer, T.W.; Kaanders, J.H.; Bussink, J. Interaction between hypoxia, AKT and HIF-1 signaling in HNSCC and NSCLC: Implications for future treatment strategies. Future Sci. OA. 2016, 2, FSO84. [Google Scholar] [CrossRef]
  61. Kazyken, D.; Magnuson, B.; Bodur, C.; Acosta-Jaquez, H.A.; Zhang, D.; Tong, X.; Barnes, T.M.; Steinl, G.K.; Patterson, N.E.; Altheim, C.H.; et al. AMPK directly activates mTORC2 to promote cell survival during acute energetic stress. Sci. Signal. 2019, 12, eaav3249. [Google Scholar] [CrossRef]
  62. Jäger, S.; Handschin, C.; St-Pierre, J.; Spiegelman, B.M. AMP-activated protein kinase (AMPK) action in skeletal muscle via direct phosphorylation of PGC-1alpha. Proc. Natl. Acad. Sci. USA 2007, 104, 12017–12022. [Google Scholar] [CrossRef]
  63. Baenke, F.; Chaneton, B.; Smith, M.; Van Den Broek, N.; Hogan, K.; Tang, H.; Viros, A.; Martin, M.; Galbraith, L.; Girotti, M.R.; et al. Resistance to BRAF inhibitors induces glutamine dependency in melanoma cells. Mol. Oncol. 2016, 10, 73–84. [Google Scholar] [CrossRef]
  64. Hernandez-Davies, J.E.; Tran, T.Q.; Reid, M.A.; Rosales, K.R.; Lowman, X.H.; Pan, M.; Moriceau, G.; Yang, Y.; Wu, J.; Lo, R.S.; et al. Vemurafenib resistance reprograms melanoma cells towards glutamine dependence. J. Transl. Med. 2015, 13, 210. [Google Scholar] [CrossRef]
  65. Wang, Q.; Beaumont, K.A.; Otte, N.J.; Font, J.; Bailey, C.G.; van Geldermalsen, M.; Sharp, D.M.; Tiffen, J.C.; Ryan, R.M.; Jormakka, M.; et al. Targeting glutamine transport to suppress melanoma cell growth. Int. J. Cancer 2014, 135, 1060–1071. [Google Scholar] [CrossRef]
  66. Shin, S.S.; Jeong, B.S.; Wall, B.A.; Li, J.; Shan, N.L.; Wen, Y.; Goydos, J.S.; Chen, S. Participation of xCT in melanoma cell proliferation in vitro and tumorigenesis in vivo. Oncogenesis 2018, 7, 86. [Google Scholar] [CrossRef]
  67. Samanta, D.; Semenza, G.L. Maintenance of redox homeostasis by hypoxia-inducible factors. Redox Biol. 2017, 13, 331–335. [Google Scholar] [CrossRef]
  68. Koppula, P.; Zhang, Y.; Shi, J.; Li, W.; Gan, B. The glutamate/cystine antiporter SLC7A11/xCT enhances cancer cell dependency on glucose by exporting glutamate. J. Biol. Chem. 2017, 292, 14240–14249. [Google Scholar] [CrossRef] [Green Version]
  69. Zhang, G.; Frederick, D.T.; Wu, L.; Wei, Z.; Krepler, C.; Srinivasan, S.; Chae, Y.C.; Xu, X.; Choi, H.; Dimwamwa, E.; et al. Targeting mitochondrial biogenesis to overcome drug resistance to MAPK inhibitors. J. Clin. Invest. 2016, 126, 1834–1856. [Google Scholar] [CrossRef]
  70. Müller, J.; Krijgsman, O.; Tsoi, J.; Robert, L.; Hugo, W.; Song, C.; Kong, X.; Possik, P.A.; Cornelissen-Steijger, P.D.; Geukes Foppen, M.H.; et al. Low MITF/AXL ratio predicts early resistance to multiple targeted drugs in melanoma. Nat. Commun. 2014, 5, 5712. [Google Scholar] [CrossRef]
  71. Ji, Z.; Erin Chen, Y.; Kumar, R.; Taylor, M.; Jenny Njauw, C.N.; Miao, B.; Frederick, D.T.; Wargo, J.A.; Flaherty, K.T.; Jönsson, G.; et al. MITF Modulates Therapeutic Resistance through EGFR Signaling. J. Invest. Dermatol. 2015, 135, 1863–1872. [Google Scholar] [CrossRef] [Green Version]
  72. Tu, T.J.; Ma, M.W.; Monni, S.; Rose, A.E.; Yee, H.; Darvishian, F.; Polsky, D.; Berman, R.S.; Shapiro, R.L.; Pavlick, A.C.; et al. A high proliferative index of recurrent melanoma is associated with worse survival. Oncology 2011, 80, 181–187. [Google Scholar] [CrossRef]
  73. Gimotty, P.A.; Van Belle, P.; Elder, D.E.; Murry, T.; Montone, K.T.; Xu, X.; Hotz, S.; Raines, S.; Ming, M.E.; Wahl, P.; et al. Biologic and prognostic significance of dermal Ki67 expression, mitoses, and tumorigenicity in thin invasive cutaneous melanoma. J. Clin. Oncol. 2005, 23, 8048–8056. [Google Scholar] [CrossRef]
  74. Tong, W.W.; Tong, G.H.; Liu, Y. Cancer stem cells and hypoxia-inducible factors (Review). Int. J. Oncol. 2018, 53, 469–476. [Google Scholar] [CrossRef] [Green Version]
  75. Huang, X.; Trinh, T.; Aljoufi, A.; Broxmeyer, H.E. Hypoxia Signaling Pathway in Stem Cell Regulation: Good and Evil. Curr. Stem Cell Rep. 2018, 4, 149–157. [Google Scholar] [CrossRef]
  76. Cheli, Y.; Bonnazi, V.F.; Jacquel, A.; Allegra, M.; De Donatis, G.M.; Bahadoran, P.; Bertolotto, C.; Ballotti, R. CD271 is an imperfect marker for melanoma initiating cells. Oncotarget 2014, 5, 5272–5283. [Google Scholar] [CrossRef] [Green Version]
  77. Marconi, A.; Borroni, R.G.; Truzzi, F.; Longo, C.; Pistoni, F.; Pellacani, G.; Pincelli, C. Hypoxia-Inducible Factor-1α and CD271 inversely correlate with melanoma invasiveness. Exp. Dermatol. 2015, 24, 396–398. [Google Scholar] [CrossRef]
  78. Saltari, A.; Truzzi, F.; Quadri, M.; Lotti, R.; Palazzo, E.; Grisendi, G.; Tiso, N.; Marconi, A.; Pincelli, C. CD271 Down-Regulation Promotes Melanoma Progression and Invasion in Three-Dimensional Models and in Zebrafish. J. Invest. Dermatol. 2016, 136, 2049–2058. [Google Scholar] [CrossRef] [Green Version]
  79. Intlekofer, A.M.; Finley, L.W.S. Metabolic signatures of cancer cells and stem cells. Nat. Metab. 2019, 1, 177–188. [Google Scholar] [CrossRef] [Green Version]
  80. Bhadury, J.; Einarsdottir, B.O.; Podraza, A.; Bagge, R.O.; Stierner, U.; Ny, L.; Dávila López, M.; Nilsson, J.A. Hypoxia-regulated gene expression explains differences between melanoma cell line-derived xenografts and patient-derived xenografts. Oncotarget 2016, 7, 23801–23811. [Google Scholar] [CrossRef] [Green Version]
  81. Sztiller-Sikorska, M.; Koprowska, K.; Jakubowska, J.; Zalesna, I.; Stasiak, M.; Duechler, M.; Czyz, M.E. Sphere formation and self-renewal capacity of melanoma cells is affected by the microenvironment. Melanoma Res. 2012, 22, 215–224. [Google Scholar] [CrossRef]
Figure 1. Proliferation and survival of melanoma cells cultured in 6% O2 greatly depend on melanoma cell density. (A) Proliferation time-courses of melanoma cells of varying seeding densities in 21% O2 and 6% O2. Cell proliferation was measured as changes in the occupied area over time using IncuCyte, and shown as fold-change relative to the confluence at time T0. (B) Representative dot plots of Annexin V/propidium iodine (PI) stained cells after 72 h in culture, showing early apoptotic (blue), late apoptotic (red), and necrotic (green) cells. Numbers in each quadrant are percentages of cells (mean values of three biological replicates ± SD). DMBC, melanoma cell populations obtained at the Department of Molecular Biology of Cancer from surgical specimens.
Figure 1. Proliferation and survival of melanoma cells cultured in 6% O2 greatly depend on melanoma cell density. (A) Proliferation time-courses of melanoma cells of varying seeding densities in 21% O2 and 6% O2. Cell proliferation was measured as changes in the occupied area over time using IncuCyte, and shown as fold-change relative to the confluence at time T0. (B) Representative dot plots of Annexin V/propidium iodine (PI) stained cells after 72 h in culture, showing early apoptotic (blue), late apoptotic (red), and necrotic (green) cells. Numbers in each quadrant are percentages of cells (mean values of three biological replicates ± SD). DMBC, melanoma cell populations obtained at the Department of Molecular Biology of Cancer from surgical specimens.
Ijms 20 04203 g001
Figure 2. Phenotypic alterations in melanoma cells in different oxygen concentrations depend on initial phenotype of melanoma cells. Melanoma cells were incubated at different oxygen concentrations, 21% O2 (hyperoxia) 6% O2 (normoxia) and 1% O2 (hypoxia). Protein levels of (A) HIF (hypoxia-inducible factor)-1α and HIF-2α and (B) p-PDH (phosphorylated pyruvate dehydrogenase, inactive form), PDH (total), p-AKT (phosphorylated protein kinase B) and total AKT were assessed by immunoblotting with representative images shown. β-actin was used as a loading control. The percentages of (C) NGFR (nerve growth factor receptor)-positive cells and (D) MITF (microphthalmia-associated transcription factor)-positive cells in normoxia and hypoxia relative to the standard culture conditions (21% O2). MITF-positive cells were almost undetectable in DMBC12 cell population (see Figure S1). n = 3, except for hypoxia (n = 2). Differences are considered significant at * p < 0.05, ** p < 0.01, *** p < 0.001.
Figure 2. Phenotypic alterations in melanoma cells in different oxygen concentrations depend on initial phenotype of melanoma cells. Melanoma cells were incubated at different oxygen concentrations, 21% O2 (hyperoxia) 6% O2 (normoxia) and 1% O2 (hypoxia). Protein levels of (A) HIF (hypoxia-inducible factor)-1α and HIF-2α and (B) p-PDH (phosphorylated pyruvate dehydrogenase, inactive form), PDH (total), p-AKT (phosphorylated protein kinase B) and total AKT were assessed by immunoblotting with representative images shown. β-actin was used as a loading control. The percentages of (C) NGFR (nerve growth factor receptor)-positive cells and (D) MITF (microphthalmia-associated transcription factor)-positive cells in normoxia and hypoxia relative to the standard culture conditions (21% O2). MITF-positive cells were almost undetectable in DMBC12 cell population (see Figure S1). n = 3, except for hypoxia (n = 2). Differences are considered significant at * p < 0.05, ** p < 0.01, *** p < 0.001.
Ijms 20 04203 g002
Figure 3. Normoxia stimulates the expression of genes associated with glucose metabolism and to the lower extent with glutamine metabolism in cell line-dependent manner. (A) Transcript levels of GLUT1 (glucose transporter 1), PDK1 (pyruvate dehydrogenase kinase 1) and HK2 (hexokinase 2) in melanoma cells incubated in the presence of 21% O2, 6% O2 or 1% O2 for 24 h were determined by qRT-PCR and normalized to the expression of a reference gene RPS17. Gene expression in 6% O2 and 1% O2 is presented relative to the expression in 21% O2. (B) Transcript levels of GLUT1, PDK1 and HK2 in melanoma cells cultured in the presence of 6% O2 for at least 3 weeks (established 6% O2 culture) relative to their levels in cells cultured in 21% O2. (C) Transcript levels of GLS (glutaminase), SLC1A5 (solute carrier family 1 member 5) and SLC7A11 (solute carrier family 7 member 11 transporter) in melanoma cells after 24 h incubation in 21% O2, 6% O2 and 1% O2, or (D) in the established 6% O2 culture, relative to their levels in 21% O2. Bars represent mean values of 3-4 biological replicates ± SD. Differences are considered significant at * p < 0.05, ** p < 0.01 or *** p < 0.001.
Figure 3. Normoxia stimulates the expression of genes associated with glucose metabolism and to the lower extent with glutamine metabolism in cell line-dependent manner. (A) Transcript levels of GLUT1 (glucose transporter 1), PDK1 (pyruvate dehydrogenase kinase 1) and HK2 (hexokinase 2) in melanoma cells incubated in the presence of 21% O2, 6% O2 or 1% O2 for 24 h were determined by qRT-PCR and normalized to the expression of a reference gene RPS17. Gene expression in 6% O2 and 1% O2 is presented relative to the expression in 21% O2. (B) Transcript levels of GLUT1, PDK1 and HK2 in melanoma cells cultured in the presence of 6% O2 for at least 3 weeks (established 6% O2 culture) relative to their levels in cells cultured in 21% O2. (C) Transcript levels of GLS (glutaminase), SLC1A5 (solute carrier family 1 member 5) and SLC7A11 (solute carrier family 7 member 11 transporter) in melanoma cells after 24 h incubation in 21% O2, 6% O2 and 1% O2, or (D) in the established 6% O2 culture, relative to their levels in 21% O2. Bars represent mean values of 3-4 biological replicates ± SD. Differences are considered significant at * p < 0.05, ** p < 0.01 or *** p < 0.001.
Ijms 20 04203 g003
Figure 4. PGC1α (peroxisome proliferator-activated receptor γ coactivator 1 alpha) protein accumulates in MITFlow melanoma cells. (A) Transcript levels of MITF-M and PGC1α were determined by qRT-PCR and normalized to the expression of a reference gene RPS17. Gene expression in DMBC12 and DMBC28 cells is shown relative to the expression in DMBC17 cells. Bars represent mean values of 3 biological replicates ± SD. Differences are considered significant at ** p < 0.01, *** p < 0.001. (B) Protein levels of MITF and PGC1α in different oxygen concentrations were assessed by western blot with representative images shown (n = 3). β-actin was used as a loading control. (C) For protein decay assay, melanoma cells were treated with 100 μg/mL cycloheximide (CHX) and protein levels of PGC1α and MCL-1 (myeloid cell leukemia 1) were determined at indicated time points. MCL-1, having short half-life, served as control of protein stability. GAPDH was used as a loading control. (D) Quantification of average protein levels of PGC1α calculated from two independent experiments.
Figure 4. PGC1α (peroxisome proliferator-activated receptor γ coactivator 1 alpha) protein accumulates in MITFlow melanoma cells. (A) Transcript levels of MITF-M and PGC1α were determined by qRT-PCR and normalized to the expression of a reference gene RPS17. Gene expression in DMBC12 and DMBC28 cells is shown relative to the expression in DMBC17 cells. Bars represent mean values of 3 biological replicates ± SD. Differences are considered significant at ** p < 0.01, *** p < 0.001. (B) Protein levels of MITF and PGC1α in different oxygen concentrations were assessed by western blot with representative images shown (n = 3). β-actin was used as a loading control. (C) For protein decay assay, melanoma cells were treated with 100 μg/mL cycloheximide (CHX) and protein levels of PGC1α and MCL-1 (myeloid cell leukemia 1) were determined at indicated time points. MCL-1, having short half-life, served as control of protein stability. GAPDH was used as a loading control. (D) Quantification of average protein levels of PGC1α calculated from two independent experiments.
Ijms 20 04203 g004
Figure 5. Vemurafenib and trametinib while targeting the MAPK/ERK pathway irrespective of oxygen concentration, affect MITF-positive and Ki-67-positive melanoma cells in oxygen- and cell line-dependent manner. (A) The effects of trametinib (TRA) and vemurafenib (VEM) on ERK1/2 activity and MITF level in 21% and 6% O2 were determined by Western blotting with representative images shown (n = 3). β-actin was used as a loading control. (B) The percentage of MITF-positive cells in control cultures and after treatment with either 50 nM trametinib (TRA) or 10 µM vemurafenib (VEM) are shown as mean values ± SD (n = 3). Results for DMBC12 were not quantified due to almost undetectable MITF level in control and lack of changes in drug-treated cells. Differences are considered significant at * p < 0.05 ** p < 0.01 *** p < 0.001. See Figure S1 for representative dot plots. (C) The percentage of Ki-67-positive cells in control cultures and after treatment with 50 nM trametinib or 10 µM vemurafenib are shown as mean values of 3 biological replicates ± SD. Differences are considered significant at * p < 0.05 ** p < 0.01 *** p < 0.001. See Figure S1 for representative dot plots. (D) Representative density plots of dual-stained DMBC28 cells with anti-MITF and anti-Ki-67 antibodies after 48 h of drug treatment. Average percentages of cells positive for one or both markers are shown (n = 2).
Figure 5. Vemurafenib and trametinib while targeting the MAPK/ERK pathway irrespective of oxygen concentration, affect MITF-positive and Ki-67-positive melanoma cells in oxygen- and cell line-dependent manner. (A) The effects of trametinib (TRA) and vemurafenib (VEM) on ERK1/2 activity and MITF level in 21% and 6% O2 were determined by Western blotting with representative images shown (n = 3). β-actin was used as a loading control. (B) The percentage of MITF-positive cells in control cultures and after treatment with either 50 nM trametinib (TRA) or 10 µM vemurafenib (VEM) are shown as mean values ± SD (n = 3). Results for DMBC12 were not quantified due to almost undetectable MITF level in control and lack of changes in drug-treated cells. Differences are considered significant at * p < 0.05 ** p < 0.01 *** p < 0.001. See Figure S1 for representative dot plots. (C) The percentage of Ki-67-positive cells in control cultures and after treatment with 50 nM trametinib or 10 µM vemurafenib are shown as mean values of 3 biological replicates ± SD. Differences are considered significant at * p < 0.05 ** p < 0.01 *** p < 0.001. See Figure S1 for representative dot plots. (D) Representative density plots of dual-stained DMBC28 cells with anti-MITF and anti-Ki-67 antibodies after 48 h of drug treatment. Average percentages of cells positive for one or both markers are shown (n = 2).
Ijms 20 04203 g005
Figure 6. Changes in the protein level of PGC1α during drug-induced apoptosis. Protein levels of PARP and PGC1α after 44 h of treatment with vemurafenib (VEM) and trametinib (TRA) at indicated concentrations in normoxia (6% O2) and hyperoxia (21% O2) were assessed by western blotting with representative images shown (n = 3). Upper and lower arrows indicate full-length and cleaved PARP (cPARP), respectively. β-actin was used as a loading control.
Figure 6. Changes in the protein level of PGC1α during drug-induced apoptosis. Protein levels of PARP and PGC1α after 44 h of treatment with vemurafenib (VEM) and trametinib (TRA) at indicated concentrations in normoxia (6% O2) and hyperoxia (21% O2) were assessed by western blotting with representative images shown (n = 3). Upper and lower arrows indicate full-length and cleaved PARP (cPARP), respectively. β-actin was used as a loading control.
Ijms 20 04203 g006
Figure 7. Effects of trametinib (TRA) and vemurafenib (VEM) on VEGF transcript level assessed at different oxygen concentrations. VEGF transcript levels were determined by qRT-PCR and normalized to the expression of a reference gene RPS17. (A) Transcript levels of VEGF were determined in cell cultures grown in hyperoxia by qRT-PCR and their expression is compared between cell lines with transcript levels in DMBC17 cells set as one. (B) Transcript level of VEGF in 6% O2 and 1% O2 is presented relative to its level in 21% O2. (C) Changes in VEGF transcript level after drug treatment is presented relative to its expression in control. Bars represent mean values of 3-4 independent experiments ± SD. Differences are considered significant at * p < 0.05, ** p < 0.01.
Figure 7. Effects of trametinib (TRA) and vemurafenib (VEM) on VEGF transcript level assessed at different oxygen concentrations. VEGF transcript levels were determined by qRT-PCR and normalized to the expression of a reference gene RPS17. (A) Transcript levels of VEGF were determined in cell cultures grown in hyperoxia by qRT-PCR and their expression is compared between cell lines with transcript levels in DMBC17 cells set as one. (B) Transcript level of VEGF in 6% O2 and 1% O2 is presented relative to its level in 21% O2. (C) Changes in VEGF transcript level after drug treatment is presented relative to its expression in control. Bars represent mean values of 3-4 independent experiments ± SD. Differences are considered significant at * p < 0.05, ** p < 0.01.
Ijms 20 04203 g007
Figure 8. Effects of trametinib (TRA) and vemurafenib (VEM) on the expression of genes related to glutamine and glucose metabolism assessed at different oxygen concentrations. Drug-induced changes in transcript levels of glutamine (A) and glucose (B) metabolism-related genes after 24 h of treatment with either vemurafenib or trametinib. qRT-PCR data were normalized to the expression of a reference gene RPS17. Bars represent mean values of 3-4 biological replicates ± SD. Statistically significant differences are indicated: * p < 0.05, ** p < 0.01 or *** p < 0.001.
Figure 8. Effects of trametinib (TRA) and vemurafenib (VEM) on the expression of genes related to glutamine and glucose metabolism assessed at different oxygen concentrations. Drug-induced changes in transcript levels of glutamine (A) and glucose (B) metabolism-related genes after 24 h of treatment with either vemurafenib or trametinib. qRT-PCR data were normalized to the expression of a reference gene RPS17. Bars represent mean values of 3-4 biological replicates ± SD. Statistically significant differences are indicated: * p < 0.05, ** p < 0.01 or *** p < 0.001.
Ijms 20 04203 g008

Share and Cite

MDPI and ACS Style

Osrodek, M.; Hartman, M.L.; Czyz, M. Physiologically Relevant Oxygen Concentration (6% O2) as an Important Component of the Microenvironment Impacting Melanoma Phenotype and Melanoma Response to Targeted Therapeutics In Vitro. Int. J. Mol. Sci. 2019, 20, 4203. https://doi.org/10.3390/ijms20174203

AMA Style

Osrodek M, Hartman ML, Czyz M. Physiologically Relevant Oxygen Concentration (6% O2) as an Important Component of the Microenvironment Impacting Melanoma Phenotype and Melanoma Response to Targeted Therapeutics In Vitro. International Journal of Molecular Sciences. 2019; 20(17):4203. https://doi.org/10.3390/ijms20174203

Chicago/Turabian Style

Osrodek, Marta, Mariusz L. Hartman, and Malgorzata Czyz. 2019. "Physiologically Relevant Oxygen Concentration (6% O2) as an Important Component of the Microenvironment Impacting Melanoma Phenotype and Melanoma Response to Targeted Therapeutics In Vitro" International Journal of Molecular Sciences 20, no. 17: 4203. https://doi.org/10.3390/ijms20174203

APA Style

Osrodek, M., Hartman, M. L., & Czyz, M. (2019). Physiologically Relevant Oxygen Concentration (6% O2) as an Important Component of the Microenvironment Impacting Melanoma Phenotype and Melanoma Response to Targeted Therapeutics In Vitro. International Journal of Molecular Sciences, 20(17), 4203. https://doi.org/10.3390/ijms20174203

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop