Abstract
Tumor heterogeneity is a major challenge for cancer treatment, especially due to the presence of various subpopulations with stem cell or progenitor cell properties. In mouse melanomas, both CD34+p75− (CD34+) and CD34−p75− (CD34−) tumor subpopulations were characterized as melanoma-propagating cells (MPC) that exhibit some of those key features. However, these two subpopulations differ from each other in tumorigenic potential, ability to recapitulate heterogeneity, and chemoresistance. In this study, we demonstrate that CD34+ and CD34− subpopulations carrying the BRAFV600E mutation confer differential sensitivity to targeted BRAF inhibition. Through elevated KDM5B expression, melanoma cells shift toward a more drug-tolerant, CD34− state upon exposure to BRAF inhibitor or combined BRAF inhibitor and MEK inhibitor treatment. KDM5B loss or inhibition shifts melanoma cells to the more BRAF inhibitor–sensitive CD34+ state. These results support that KDM5B is a critical epigenetic regulator that governs the transition of key MPC subpopulations with distinct drug sensitivity. This study also emphasizes the importance of continuing to advance our understanding of intratumor heterogeneity and ultimately develop novel therapeutics by altering the heterogeneous characteristics of melanoma.
Introduction
Intratumor heterogeneity serves as a significant clinical barrier to finding effective and long-lasting treatments for various types of cancer (1). Many reports suggest that melanoma does not follow the hierarchical cancer stem cell model because tumorigenic cells within a single melanoma tumor represent more than a rare subpopulation (2–5). Meanwhile, several other studies argue that there exist more subtleties to this conclusion than previously suggested (6, 7). Our laboratory demonstrated that mouse melanoma tumors driven by different genetic alterations are composed of three major subpopulations. These three subpopulations, CD34+p75− (CD34+), CD34−p75− (CD34−), and CD34−p75+ (p75+), exhibit distinguishing cellular behaviors (8). Specifically, CD34+ cells are highly tumorigenic in vivo, whereas p75+ cells are rarely tumorigenic; CD34− cells exhibit an intermediate level of tumorigenic potential but are capable of recapitulating tumor heterogeneity in vitro and in vivo. Thus, both CD34+ and CD34− cells are referred to as melanoma-propagating cells (MPC) and are speculated to play a critical role in tumor development. In addition, these three subpopulations demonstrate varying levels of resistance to traditional cytotoxic chemotherapeutics in vitro, suggesting that they may also respond differently to more targeted treatments (8). Overall, although we have made progress in understanding the cellular characteristics representing the different tumor cell subpopulations in melanoma, we have not yet identified key regulators mediating the conversion of CD34+ and CD34− MPCs. This is an essential step to provide effective therapeutic strategies for melanoma patients considering the critical cellular features unique to these subpopulations.
Several epigenetic regulators have been used to identify heterogeneous tumor subpopulations with distinct cellular phenotypes. Non–small cell lung cancer cell lines along with several other human cancer cell lines contained a reversibly drug-resistant subpopulation that required KDM5A to maintain the drug-tolerant state (9–11). A slow-cycling subpopulation of melanoma cells expressing high levels of histone demethylase KDM5B (also known as JARID1B) was shown to be important in maintaining melanoma growth and intrinsically resistant to various cancer drugs (12, 13). Moreover, although CD34+, CD34−, and p75+ cells have been identified in mouse melanoma tumors with different genetic alterations, subpopulations expressing the same distribution of markers exhibit similar progenitor cell properties (8). Taken together, these findings indicate that epigenetic variation likely exists between the two MPC and p75+ subpopulations and may also be the explanation for their diverse cellular characteristics.
In addition, cancer cells exploit intratumor heterogeneity as a defense mechanism in response to drug therapy (9, 13, 14). This highlights the clinical need to understand the regulation of intratumor heterogeneity during tumor development and in the context of therapeutic treatment. In melanoma, over 50% of patients harbor the V600 mutation in BRAF kinase. This particular genetic alteration leads to constitutive activation of the MAPK pathway and contributes to uncontrolled cell proliferation (15). Unfortunately, resistance to targeted BRAF inhibitors (BRAFi) such as vemurafenib presents a continuous challenge in the clinic. Multiple resistance mechanisms have been described, including reactivation of the MAPK pathway, upregulation of the compensatory PI3K signaling, and nongenetic modifications (16, 17). In this study, we examine how CD34+ and CD34− MPC subpopulations respond to targeted therapies to better understand the link between intratumor heterogeneity and drug response. We further demonstrate that histone demethylase KDM5B is a key regulator of the conversion between MPC subpopulations of mouse melanoma.
Materials and Methods
Cell culture
Yale University Mouse Melanoma cell lines (YUMM1.1, YUMM1.3, YUMM1.7, YUMM1.9, YUMM2.1, YUMM3.3, YUMM3.4, and YUMM5.1) used in this study were characterized previously (18). Cells were maintained in DMEM/F12 (1:1) media containing 10% heat-inactivated FBS, 1% penicillin streptomycin, and 1% MEM nonessential amino acids (Gibco by Life Technologies) at 37°C and 5% CO2. BRAFi-resistant YUMM1.7 (YUMM1.7R) cells were derived by continuously exposing YUMM1.7 cells to 2 μmol/L BRAFi (PLX4032, BioVision 2235-5; ref. 19) for 8 weeks as described previously (20) and maintained in 2 μmol/L BRAFi media unless indicated otherwise. Human melanoma cell lines 501mel and YUMAC (within 10 passages after being obtained from Dr. Ruth Halaban from Yale University; ref. 21) were grown in similar conditions except using OPTI-MEM (Invitrogen). Cell lines were authenticated using short tandem repeat profiling performed at the DNA Analysis Facility on Science Hill at Yale University. Sigma LookOut Mycoplasma PCR detection kit (Sigma-Aldrich, MP0035) was used for mycoplasma detection in cell culture.
For drug treatment, YUMM cells were treated with either 1.5 μmol/L BRAFi or with a combination of 1.5 μmol/L BRAFi and MEK inhibitor (MEKi; Trametinib, Santa Cruz Biotechnology sc-364639). Cells were harvested for Western blotting, flow cytometry, and reverse transcription and quantitative PCR (RT-qPCR) at end point. For KDM5 inhibitor (KDM5i) treatment, YUMM1.7 cells were pretreated in 1 μmol/L KDM5-C70 (Xcess Biosciences Inc. M60192-2s; refs. 22, 23), 1 μmol/L CPI-48 (NCGC00488278, prepared according to patent WO2015/135094; ref. 24), or 100 μmol/L YUKA1 (ChemBridge Hit2Lead 7870547; ref. 11) for 5 days before being cotreated with 1.5 μmol/L BRAFi for an additional 3 days. On the 8th day, cells were harvested for similar assays and analyzed as described above.
Flow cytometry analysis and cell sorting
At the end of each in vitro experiment, cells were washed with PBS and trypsinized with 0.25% Trypsin-EDTA (Gibco by Life Technologies). Cells were stained with rat anti-CD34 (ThermoFisher 14-0341-82) and rabbit anti-Nerve Growth Factor Receptor/p75 (Millipore Sigma AB1554) in PBS supplemented with 1% BSA (Sigma-Aldrich A9647) at 4°C for 30 minutes. Cells were washed 3 times with fluorescence-activated cell sorting (FACS) buffer and stained with Brilliant Violet 421 Donkey anti-rabbit IgG (BioLegend 406410), Alexa Fluor 647 Goat anti-rat IgG (ThermoFisher A-21247), and propidium iodide (BioLegend 421301) at 4°C for 15 minutes. Cells were washed 3 times, resuspended in PBS, and analyzed by BD LSRII Flow Cytometer. Plots were generated using FlowJo 9.9.6. To sort CD34+ cells, YUMM1.7 cells were stained using the same protocol and subjected to single-cell sorting using BD FACSAria into a 96-well plate. Individual cells were expanded and tested by flow cytometry to ensure CD34+ purity. CD34+ and CD34− cells used for the RT-qPCR experiment were sorted using the same staining protocol but collected using 15 mL conical tubes using the BD FACSAria or Sony SH800 Cell Sorter. To analyze subpopulation distribution in grafted tumors by flow cytometry, tumors were dissociated by finely mincing and incubating in DMEM/F12 (1:1) containing 1 mg/mL Collagenase/Dispase (Sigma-Aldrich 10269638001) and 1X DNase I (Qiagen 79254) at 37°C for 30 minutes. Dissociated tumors were titrated and filtered through a 70 μm filter to yield a single-cell suspension before FACS staining. Cells were stained using a similar protocol described above with an additional anti-CD45 APC-Cy7 antibody (BioLegend 103115).
Western blot analysis
Mouse cells were collected following digestion with 0.25% Trypsin-EDTA and lysed on ice with high salt lysis buffer [50 mmol/L Tris-HCl pH 7.9, 0.1 mmol/L EDTA pH 8.0, 320 mmol/L NaCl, 0.5% NP40, 10% Glycerol, and 1x Protease Inhibitor cocktail (Roche 11873580001)] to obtain whole-cell lysates. The remaining pellets were resuspended in Laemmli buffer and sonicated to extract histones. Protein concentration of whole-cell lysates was measured by Bradford assay (Bio-Rad) with known BSA standards. Samples in Laemmli buffer were heated for 10 minutes at 95°C and loaded onto 7% (whole-cell lysates) or 15% (histones) SDS-PAGE gels. Membranes were blocked in 4% nonfat milk in Tris-buffered saline (50 mmol/L Tris-HCl, 138 mmol/L NaCl, 2.7 mmol/L KCl, pH 7.4) with 0.05% Tween (TBS-T) and incubated with primary antibodies in the same buffer or 5% BSA in TBS-T overnight at 4°C. Membranes were incubated with secondary anti-rabbit or anti-mouse antibodies for 1 hour at room temperature. Blots were visualized by ECL Western Blotting Substrate (Thermo Scientific Pierce 32106) or Immobilon Western Chemiluminescent HRP Substrate (EMD Millipore WBKLS0100) on film or with KwikQuant Imager. Primary antibodies used included rabbit anti-KDM5B (Sigma HPA027179), mouse anti-tubulin (Sigma T5168), rabbit anti-histone H3 (Abcam ab1791), rabbit anti-H3K4me3 (C42D8, Cell Signaling Technology 9751), mouse anti-vinculin (Sigma V9131), and rabbit anti-GAPDH (14C10, Cell Signaling Technology 2118).
Proliferation assay and IC50 assays
To measure cell proliferation, cells were seeded in triplicate in 96-well tissue culture plates and grown at various experimental conditions: DMEM/F12 media, dimethyl sulfoxide (DMSO, Sigma-Aldrich 67-68-5), BRAFi, and BRAFi/MEKi. At the respective time points, cell proliferation was measured using CyQUANT NF Cell Proliferation Assay Kit (Invitrogen C35006). Values within the same group at each time point were averaged and compiled with t test statistical analysis using Prism 7. IC50 values were determined by incubating cells in different concentrations of BRAFi at 37°C for 3 days. The concentrations were 0.039, 0.078, 0.156, 0.313, 0.652, 1.25, 2.5, 5, and 10 μmol/L, in addition to DMSO controls. All samples were analyzed in triplicate, and all experiments were repeated 3 times.
Tumor graft
To study tumor growth and response to targeted inhibitors, 100,000 YUMM1.7 cells were injected subcutaneously into each C57BL/6J mouse (Jackson Laboratory). Mice were maintained in the animal facility of Yale School of Medicine. All animal studies were approved by the Institutional Animal Care and Use Committee of Yale University. Tumor volume was monitored and determined using the following equation: width × length × height × 0.5233 (mm). Once tumor volume reached 300 mm3, mice in treatment groups were switched to a special diet with BRAFi (Plexxikon PLX4720-200ppm; ref. 25) or combination BRAFi/MEKi (Plexxikon PLX4720/PLX2695-200/7ppm; ref. 26) chow, whereas mice in control groups continued with the regular chow diet. At end point, mice were euthanized and tumors were harvested for flow cytometry analysis and RT-qPCR analysis.
Generation of KDM5B knockdown cell lines
Stable knockdowns of KDM5B in mouse melanoma cell lines YUMM1.7 and YUMM3.4, and human melanoma cell lines YUMAC and 501mel were performed as described previously (27, 28) using lentiviral shRNA hairpins: pLKO.1- mouse-KDM5B sh1 (targeting GCCTACATCATGTGAAAGAAT), pLKO.1- mouse-KDM5B sh2 (targeting CCTGAAATTCAGGAGCTTTAT), pLKO.1- human-KDM5B sh1 (targeting CGAGATGGAATTAACAGTCTT), pLKO.1- human-KDM5B sh2 (targeting AGGGAGATGCACTTCGATATA), and pLKO.1-shGFP (control sh; a gift from Dr. William Hahn, Dana-Farber Cancer Institute, Boston, MA). To generate stable cell lines, cells were infected with lentivirus for 24 hours and refed with fresh media with 2 μg/mL puromycin for mouse lines and 1 μg/mL puromycin for human lines until uninfected control cells were completely eliminated.
RNA purification and RT-qPCR
Total RNA was extracted from cultured mouse or human cells, or mouse tumor tissues with the RNeasy Plus Mini Kit (Qiagen 74134). cDNA was obtained using the High Capacity cDNA Reverse Transcription Kit (Applied Biosystems 4368814) according to the manufacturer's instruction. RT-qPCR was performed with CFX96 (BioRad) using Fast SYBR Green Master Mix (Applied Biosystems 43856414). All values were normalized to the level of mouse Gapdh or human 18s abundance. Data were the average of triplicate experiments ± SEM. All primers used are described in Supplementary Table S1.
Results
CD34− cells were enriched in BRAFi-resistant melanoma cells
Several mouse melanoma cell lines known as YUMMs were recently derived from genetically engineered mouse models with relevant human driver mutations (18). By flow cytometry, we identified both CD34+ and CD34− MPC subpopulations at various percentage breakdowns in these cell lines across different genotypes (Supplementary Table S2). As a result, YUMMs served as an ideal experimental model to study potential epigenetic differences between subpopulations and to identify key factors that regulate intratumor heterogeneity. In YUMM1.7R cells, a BRAFi-resistant mouse melanoma line (BRAFi IC50 = 2.8 μmol/L) derived from BRAFi-sensitive YUMM1.7 cells (BRAFi IC50 = 0.3 μmol/L), we detected a major shift in CD34 positivity compared with parental YUMM1.7 cells by flow cytometry (Fig. 1A). Specifically, YUMM1.7R was composed of mostly CD34− cells, whereas CD34+ cells existed at a higher percentage in parental YUMM1.7. It is worth noting that YUMM1.7 cells do not have p75+ subpopulations as in some other cell lines (Supplementary Table S2). To further compare the properties of these cells, we examined the expression of BRAFV600E-responsive genes Gdf15, Ldlr, Map3k1, Ypel2, and Slc43a3, which are highly expressed in BRAFi-sensitive human melanoma cell lines and are downregulated by BRAFi treatment (29). We found that, in comparison with sorted CD34− cells, sorted CD34+ cells exhibited significantly higher mRNA levels of these BRAFV600E-responsive genes (Fig. 1B, left plot). Consistently, these BRAF-responsive genes are expressed at higher levels in YUMM1.7 cells than in YUMM1.7R cells (Fig. 1B, right plot), suggesting that the resistant cells adapted to BRAFi treatment by dampening the expression of these genes. Sorted CD34+ cells also demonstrated stronger growth inhibition by BRAFi treatment than CD34− or unsorted YUMM1.7 cells consisting of both CD34+ and CD34− cells (Fig. 1C). These data suggest that CD34+ and CD34− subpopulations are intrinsically distinct in response to targeted BRAF inhibition, in addition to their known differences in tumorigenicity, ability to recapitulate heterogeneity, and chemoresistance (8). Furthermore, these findings suggest that mouse melanoma cells undergo changes upon exposure to drug treatment that alter intratumor heterogeneity and shift subpopulations toward a more drug-tolerant state.
BRAFi induced reversible subpopulation shifting from CD34+ to CD34− cells
To evaluate the effect of BRAFi on subpopulation composition in YUMMs, we treated sorted CD34+ cells with either DMSO or 1.5 μmol/L BRAFi for 3 days and subjected the cells to flow cytometry analysis (Fig. 2A). As described previously (8), sorted CD34+ cells treated with DMSO vehicle mainly remained CD34-positive (Supplementary Fig. S1A). In contrast, BRAFi induced emergence of a small subpopulation of CD34− cells from sorted CD34+ cells, as shown by the reduced ratio of %CD34+ to %CD34− cells (Fig. 2A; Supplementary Fig. S1A). Consistently, in parental YUMM1.7 cells, in a treatment time course from 1 to 6 days, the degree of subpopulation shifting from CD34+ toward CD34− increased (Fig. 2B). To examine whether this was not just a genotype-specific phenomenon, we conducted similar experiment using YUMM3.3 cells. Unlike YUMM1.7 cells (BrafV600EPten−/−Cdkn2a−/−), YUMM3.3 cells (BrafV600ECdkn2a−/−) carry wild-type Pten. Despite the difference in their genotypes, YUMM3.3 cells showed a similar shifting from CD34+ to CD34− cells after 3 days of 1.5 μmol/L BRAFi treatment (Supplementary Fig. S1B). Furthermore, we treated YUMM1.7 cells with a combination of 1.5 μmol/L BRAFi and MEKi for 3 days and observed a comparable degree of CD34+ to CD34− cell shifting compared with BRAFi treatment alone (Fig. 2B and C). These data suggest that BRAFi alone is sufficient to induce changes in MPC subpopulations. To investigate whether subpopulation shifting is stable, we withdrew BRAFi from treated YUMM1.7 cells after 3 days of treatment. The percentage of CD34+ cells was restored in both YUMM1.7 and YUMM1.7R cells after withdrawal of BRAFi for 6 days (Fig. 2D). These data strongly indicate that BRAFi-induced subpopulation shifting between CD34+ and CD34− cells is reversible within a short period of time and therefore likely caused by epigenetic changes. In addition, propidium iodide staining revealed that BRAFi treatment or long-term exposure to BRAFi did not induce significant cell death in both YUMM1.7 and YUMM1.7R cells, respectively (Supplementary Fig. S1C and S1D), suggesting that the rise of CD34− subpopulation was not driven by selective pressure.
BRAFi induced reversible upregulation of KDM5B
A small subpopulation of slow-cycling, KDM5B (JARID1B)high human melanoma cells was shown to confer intrinsic multidrug resistance (12). To determine whether KDM5B level is associated with drug resistance and varies between CD34+ and CD34− cells, we treated YUMM1.7 cells for 3 days with 1.5 μmol/L BRAFi or DMSO control, sorted CD34+ and CD34− subpopulations from YUMM1.7 cells, and evaluated Kdm5b expression in these cells. After BRAFi treatment, Kdm5b mRNA level significantly increased in both subpopulations (Fig. 3A). Consistently, KDM5B protein level increased significantly under BRAFi treatment in CD34+ YUMM1.7 cells (Fig. 3B, left plot) and in YUMM3.3 cells (Fig. 3B, right plot). Moreover, KDM5B protein level is significantly higher in YUMM1.7R cells than that in YUMM1.7 cells (Fig. 3C). Although we also noticed that Kdm5b was expressed at a slightly higher level in CD34− cells than in CD34+ cells (Fig. 3A), this small difference is unlikely to be biologically important. These results suggested that BRAFi-induced KDM5B expression is not simply due to subpopulation switching. To investigate the effects of upregulated KDM5B, we compared the levels of KDM5B substrate H3K4me3 in DMSO control and 1.5 μmol/L BRAFi-treated cells. We found that BRAFi-treated cells had a lower level of H3K4me3 than DMSO control cells, consistent with the increased expression of KDM5B after BRAFi treatment (Fig. 3B). Similarly, the increased KDM5B level in YUMM1.7R cells compared with YUMM1.7 cells correlated with decreased H3K4me3 level (Fig. 3C). We further demonstrated that withdrawing BRAFi for 6 days from treated YUMM1.7 and YUMM3.3 (Fig. 3D) and BRAFi resistant YUMM1.7R (Fig. 3E) significantly lowered their Kdm5b mRNA levels by RT-qPCR, suggesting that induction of KDM5B expression by BRAFi is reversible. Taken together, these results suggest that KDM5B plays an important role in regulating BRAFi-induced subpopulation switching.
KDM5B regulated the transition between CD34+ and CD34− subpopulations
To further investigate whether KDM5B is important in changing MPC subpopulation switching in melanoma, we generated KDM5B knockdown derivatives in different YUMM cell lines. By flow cytometry, we observed an approximately 1.8-fold increase in the percentage of CD34+ subpopulation in YUMM1.7 cells with KDM5B shRNA–mediated knockdown (77.4% or 71.8% vs. 41.8%; Fig. 4A). The percentage of CD34− cells in the KDM5B knockdown cell lines decreased correspondingly, with no emergence of p75+ cells (Fig. 4A). Similarly, knocking down KDM5B in YUMM3.4 cells also led to an increase in the percentage of CD34+ cells (Supplementary Fig. S2A). Because CD34+ and CD34− cells exhibit different cellular features in terms of proliferation and sensitivity to BRAFi, we examined whether these characteristics changed in the KDM5B knockdown cells as a result of subpopulation shifting. Indeed, compared with control shRNA cells, YUMM1.7 cells with KDM5B knockdown exhibited significantly faster cell proliferation (Fig. 4B) and increased sensitivity to BRAFi, shown by with higher growth inhibition after 3 days of 1.5 μmol/L BRAFi treatment (Fig. 4C; Supplementary Fig. S2B).
We then asked whether the histone demethylase activity of KDM5B contributed to the subpopulation transition. As no KDM5B-specific inhibitor is current available, we examined the effects of pan-KDM5 enzymatic inhibitors (KDM5i) KDM5-C70 (now referred to as C70; refs. 22, 23), CPI-48 (24, 30), which inhibits KDM5A-D, and a KDM5A/C inhibitor YUKA1 (11). We pretreated YUMM1.7 cells with 1 μmol/L KDM5i C70 for 5 days before subjecting them to BRAFi or combined BRAFi and KDM5i treatment. On the 8th day, we harvested all cells and examined subpopulation changes by flow cytometry (Fig. 4D). Consistent with the phenotype induced by KDM5B knockdown on YUMM1.7 and YUMM3.4 cells, KDM5i C70 increased the percentage of CD34+ cells of YUMM1.7 cells. Furthermore, KDM5i C70 blocked BRAFi-induced decrease of CD34+ cells (Fig. 4D). Treatment of YUMM1.7 cells with 1 μmol/L CPI-48, another pan-KDM5i, delivered consistent results as KDM5i C70 (Supplementary Fig. S2C). In contrast, pretreatment with 100 μmol/L YUKA1, a KDM5A/C inhibitor, had no effect on YUMM1.7 CD34+ and CD34− subpopulations and did not affect the BRAFi-induced phenotype (Supplementary Fig. S2C). Western blotting analyses showed that C70 and CPI-48 induced dramatic induction of the H3K4me3 levels, whereas the effect of YUKA1 was more modest, suggesting that KDM5B is the dominant H3K4me3 demethylase in these melanoma cells (Fig. 4D; Supplementary Fig. S2C). Moreover, KDM5i C70 blocked the reduction of H3K4me3 level by BRAFi (Fig. 4D). These data suggest that the enzymatic activity of KDM5B is critical for the subpopulation transition between CD34+ and CD34− MPCs.
KDM5B modulates MPC-like features of human melanoma cells
To investigate whether KDM5B plays a similar role in regulating tumor cell behaviors in human melanoma, we examined the levels of KDM5B in YUMAC and 501mel human melanoma cells treated with DMSO control, 1.5 μmol/L BRAFi for 3 days, or 1.5 μmol/L BRAFi for 3 days followed by 6 days of BRAFi withdrawal. Similar to the YUMM cells, KDM5B levels increased significantly in both YUMAC and 501mel cells after BRAFi treatment and then decreased to DMSO control level after BRAFi removal (Fig. 5A). Furthermore, knockdown of KDM5B with two independent shRNA hairpins in 501mel and YUMAC significantly increased H3K4me3 level (Fig. 5B), cell proliferation (Fig. 5C), and growth inhibition by 1.5 μmol/L BRAFi (Fig. 5D). In mouse melanoma, enhanced cell proliferation and sensitivity to BRAFi are characteristics unique to CD34+ cells as opposed to CD34− cells. Here, we demonstrate that KDM5B loss in human melanoma cell lines also leads to phenotypic changes toward a cellular state that resembles that of the CD34+ subpopulation in mouse melanoma.
In vivo BRAFi therapy triggers reversible MPC subpopulation changes
To evaluate how heterogeneous subpopulations of melanoma cells respond to BRAFi in vivo, we engrafted 100,000 YUMM1.7 cells subcutaneously in C57BL/6J mice. When tumors reached a volume of 300 mm3 (day 25), these mice were provided with BRAFi or regular chow diet. Tumor growth curves revealed that grafted YUMM1.7 tumors responded to BRAFi chow (Fig. 6A). At day 33, replaced BRAFi chow with regular chow for some BRAFi treated mice, and we call this group of mice “BRAFi chow removed” group. Tumors regrew in these mice and reached 1,000 mm3 at day 42, which allowed us to harvest tumors of similar size to those in the “regular chow” group on day 33. We then evaluated the tumors for Kdm5b level and percentage of CD34+ cells in tumors for “regular chow” group on day 33, “BRAFi chow” group on day 33, and “BRAFi chow removed” group on day 42. Using RT-qPCR, we observed an increased Kdm5b level in tumors collected on day 33 from BRAFi-treated mice compared with that in tumors from mock-treated mice, whereas this increase was partially reversed in tumors harvested on day 42 from mice with BRAFi removed (Fig. 6B). This pattern of fluctuating Kdm5b expression in grafted tumors was consistent with our observation in vitro. We analyzed these tumors by flow cytometry for CD45− cells to exclude hematopoietic cells and enrich for tumor cells. We demonstrated a reduced percentage of CD45−CD34+ cells in BRAFi-treated tumors compared with tumors that received regular chow (Fig. 6C). Tumors that grew out after BRAFi chow removal exhibited a higher percentage of CD45−CD34+ cells than BRAFi-treated tumors on day 33 (Fig. 6C). In addition, we examined the subpopulation composition of several small “baseline” tumors from day 25 before providing any BRAFi chow to exclude the possibility that the difference observed between “regular chow” and “BRAFi chow” was due to changes in tumor volume. Indeed, there was no significant difference in the percentage of CD45−CD34+ cells between “baseline” and “regular chow” tumors, suggesting that tumor size does not affect the percentage of CD45−CD34+ cells (Supplementary Fig. S3A). Furthermore, combined treatment of BRAFi and MEKi chow in grafted YUMM1.7 tumors led to an increase in Kdm5b expression (Supplementary Fig. S4A and S4B) and a reduced percentage of CD45−CD34+ cells, when compared with tumors from the control group (Supplementary Fig. S4C). These results indicate that in vivo tumor MPC subpopulations respond to BRAFi or BRAFi/MEKi therapies similarly to in vitro–cultured cells by altering KDM5B expression in a reversible manner.
Discussion
In mouse melanoma, although both CD34+ and CD34− cells are referred to as MPCs, the two subpopulations differ significantly in their tumorigenic potential, in addition to several other critical cellular features, such as cell proliferation, resistance to chemotherapy, and ability to recapitulate heterogeneity (8). Specifically, CD34− cells are capable of re-establishing cellular heterogeneity both in vitro and in vivo, whereas CD34+ cells mainly produce CD34+ progeny. In this study, we identify a key regulator that modulates the CD34+ and CD34− MPC transition and shed light on novel therapeutic strategies that exploit the heterogeneous nature of melanoma.
YUMM1.7 cells consist of approximately 75% CD34+, 25% CD34− cells, and no p75+ cells (Supplementary Table S2). Despite long-term in vitro culturing and passaging, the subpopulation composition of YUMM1.7 remains relatively consistent. Although YUMM3.3 and several other genotypically different YUMM cell lines are composed of varying percentages of CD34+, CD34−, and p75+ cells, the size of each subpopulation does not change significantly over time in vitro. Nevertheless, in YUMM1.7R cells, which are derived from parental YUMM1.7 cells, we demonstrate a noticeable increase in the percentage of CD34− cells. These observations suggest that a well-controlled mechanism regulates the conversion of CD34+ and CD34− MPCs which is altered upon exposure to targeted therapies, shifting the subpopulation equilibrium in favor of a more drug-tolerant state. Indeed, compared with CD34+ cells, CD34− cells demonstrate lower expression of several BRAFV600E-responsive genes (ref. 29; Fig. 1B), indicating that CD34− cells are likely more resistant to BRAFi treatment. Consistent with these observations, sorted CD34+cells are indeed more sensitive to BRAFi treatment than sorted CD34− and parental YUMM1.7 (Fig. 1C).
Of the MPC subpopulations of CD34+ and CD34− cells, it is previously suggested that only the CD34− subpopulation is capable of generating all subpopulations (8). Here, we provide evidence that under certain circumstances, CD34+ cells can re-establish intratumor heterogeneity under drug treatment (Fig. 2A and B). Our findings are also consistent with several previous studies that have shown stressful exposures like drug treatment could trigger phenotype switching in tumor cells as a transient defense mechanism (20, 31, 32). We observed that subpopulation shifting from CD34+ to CD34− of YUMM1.7 cells occurs as early as one day after the initial BRAFi treatment (Fig. 2B). Moreover, the more drug-tolerant CD34− subpopulation expands under BRAFi exposure and the degree of subpopulation shifting intensifies as treatment prolongs (Fig. 2B). Subpopulation shifting from CD34+ to CD34− MPCs also occurs in BRAFi/MEKi cotreated YUMM1.7 cells (Fig. 2C), suggesting that the subpopulation shifting phenotype is not limited to BRAFi monotherapy and has broader clinical implications, as combined BRAFi and MEKi treatment is becoming the standard of care (33, 34). Similar to the findings in recent studies (35, 36), we observed reversibility in subpopulation shifting upon BRAFi withdrawal. The reversible phenotype applies to both BRAFi-sensitive YUMM1.7 and YUMM3.3 cells, and BRAFi-resistant YUMM1.7R cells, suggesting that BRAFi resistance of YUMM1.7R cells is reversibly regulated.
Furthermore, we provide evidence that increase of the percentage of CD34− cells is not driven by selective pressure due to BRAFi-induced cell death, as there is no significant cell death after BRAFi treatments. It is important to note that CD34+ cells actually exhibit greater survival in the presence of chemotherapeutic agents such as cisplatin and temozolomide than CD34− cells and are the main subpopulation present after treatment (8). In contrast to BRAFi, both drugs are cytotoxic and induce excessive cell death. Therefore, loss of CD34− cells after chemotherapy is likely a consequence of selection of CD34+ cells rather than subpopulation shifting. MPC subpopulations' distinct responses to different categories of therapeutic drugs highlight a potential benefit of combining cytotoxic chemotherapies and targeted therapies (37).
Moreover, because shifting of MPC subpopulations in YUMM cells occurs relatively quickly upon initial exposure to BRAFi and is a reversible phenotype, we suspect that epigenetic regulators mediate this transition. In fact, subpopulations enriched with epigenetic regulators such as KDM5A or KDM5B confer intrinsic drug tolerability in cancers including melanoma and non–small cell lung cancer (9, 13, 38). However, it is unclear how these epigenetic regulators are modulated during the development of drug resistance. It has been reported that BRAFi treatment upregulates MITF/PGC1α signaling in BRAF-mutated melanoma cells (39). It is possible that the induced MITF/PGC1α signaling set the stage for the increased KDM5B expression. Our study shows that BRAFi treatment leads to an increase in KDM5B and decrease in H3K4me3, suggesting that induced KDM5B is functional and reshapes the epigenetic landscape. In addition, we demonstrate that knockdown of KDM5B in YUMM1.7 and YUMM3.4 increases the percentage of CD34+ cells (Fig. 4A; Supplementary Fig. S2A). Compared with shRNA control cells, YUMM1.7 cells with KDM5B shRNAs exhibit increased cell proliferation and higher sensitivity to BRAF inhibition, both of which are the cellular features of CD34+ cells (8). KDM5B loss induced cell proliferation in YUMAC, and 501mel human melanoma cells were consistent with previous studies of WM3734, WM3899, and WM35 human melanoma cells (12). Taken together, our results show that KDM5B promotes the switching from CD34+ to CD34− MPCs of YUMM1.7 and YUMM3.4 cells. Further studies with the appropriate mouse melanoma cell lines are necessary to identify and characterize the regulators of transition between MPCs and the nontumorigenic p75+ subpopulation.
We show that pan-KDM5 inhibitors C70 or CPI-48 mitigate the degree of BRAFi-sensitive CD34+ cells shifting toward more BRAFi-tolerant CD34− cells (Fig. 4D; Supplementary Fig. S2C), whereas YUKA1, an enzymatic inhibitor specific to KDM5A and KDM5C (11, 40, 41), has minimal impact on subpopulation conversion (Supplementary Fig. S2C). These data suggest that the enzymatic activity of KDM5B plays major roles in subpopulation switching. However, the effects of pan-KDM5i are milder than that of KDM5B knockdown (compare Fig. 4D and Supplementary Fig. S2C with Fig. 4A). Although we cannot fully exclude the possibility that pan-KDM5 inhibitor has broader effects on multiple KDM5 family members than knockdown of only KDM5B, these differences suggest that both histone demethylase catalytic activity-dependent and -independent functions of KDM5B contribute to its regulation of subpopulation shifting. Consistent with this idea, in addition to the JmjC domain that mediates histone demethylation, KDM5 enzymes have several other domains involved in binding specific histone markers or interacting with other chromatin complex (42–44). Nonetheless, the KDM5i treatment tactic has the potential to prevent melanoma cells from transitioning to a more intrinsic drug-tolerant state, limiting the potential for developing acquired resistance. Taken together, these data emphasize the importance of understanding the precise cellular features and roles of different tumor subpopulations and identifying their key regulator(s) to develop novel and effective therapeutic paradigms.
The tumor microenvironment has been reported to influence tumor cell plasticity (45–47). In this study, by using immune-competent C57BL/6J mice grafted with YUMM1.7 cells, we provide evidence that a similar subpopulation shifting phenotype from CD34+ to CD34− tumor cells occurs with BRAFi or combined BRAFi and MEKi treatments in vivo. In addition, treated tumors are composed of a higher percentage of CD34− cells that express elevated levels of Kdm5b. Consistent with our in vitro findings, the subpopulation phenotype is also reversible upon treatment withdrawal, and Kdm5b levels fluctuate accordingly. Although we have not studied the interaction between tumor subpopulations and the immune/stromal cells in the microenvironment, it is possible that the more drug-tolerant CD34− populations are more immunosuppressive than CD34+ cells in immunocompetent mice. Consistently, studies have shown benefits to combining targeted therapies and immune checkpoint inhibitors for melanoma treatment (48, 49). Therefore, it will be interesting to examine whether the tumor microenvironment plays any role in shifting CD34+ and CD34− subpopulations in vivo. This could offer insight into designing more complex and effective treatment plans.
In summary, our study suggests that KDM5B is a key factor modulating the conversion of CD34+ and CD34− MPC subpopulations. It also demonstrates that by altering expression of KDM5B and therefore the level of H3K4me3, melanoma cells treated with BRAFi or BRAFi/MEKi are capable of shifting toward a more drug-tolerant CD34− state as a defense mechanism. Furthermore, this study suggests that novel combination treatment using pan-KDM5 or KDM5B-specific inhibitors with BRAFi could be used to increase drug sensitivity of tumor cells to BRAFi and thus improve overall therapeutic efficacy. It sheds light on the importance to continue advancing our understanding of intratumor heterogeneity in melanoma and designing more effective treatment paradigms by manipulating the tumor heterogeneity.
Disclosure of Potential Conflicts of Interest
M.W. Bosenberg is a consultant at Eli Lilly and Company. No potential conflicts of interest were disclosed by the other authors.
Disclaimer
The content of this publication does not necessarily reflect the views or policies of the Department of Health and Human Services, nor does mention of trade names, commercial products, or organizations imply endorsement by the U.S. Government.
Authors' Contributions
Conception and design: X. Liu, S.M. Zhang, A. Simeonov, M.W. Bosenberg, Q. Yan
Development of methodology: X. Liu, D.J. Jansen, A. Simeonov, M.W. Bosenberg, Q. Yan
Acquisition of data (provided animals, acquired and managed patients, provided facilities, etc.): X. Liu, S.M. Zhang, M.K. McGeary, I. Krykbaeva, A. Simeonov
Analysis and interpretation of data (e.g., statistical analysis, biostatistics, computational analysis): X. Liu, S.M. Zhang, M.K. McGeary, I. Krykbaeva, A. Simeonov, D.P. Kelly, M.W. Bosenberg, Q. Yan
Writing, review, and/or revision of the manuscript: X. Liu, S.M. Zhang, L. Lai, D.J. Jansen, A. Simeonov, M.D. Hall, D.P. Kelly, M.W. Bosenberg, Q. Yan
Administrative, technical, or material support (i.e., reporting or organizing data, constructing databases): X. Liu, S.C. Kales, M.D. Hall, M.W. Bosenberg, Q. Yan
Study supervision: A. Simeonov, M.D. Hall, M.W. Bosenberg, Q. Yan
Other (synthesized and formulated compounds that were used in study): D.J. Jansen
Acknowledgments
We would like to thank Dr. Ruth Halaban from the Department of Dermatology at Yale University for kindly providing the human melanoma cell lines used in this study. We thank Plexxikon, Inc. for providing PLX4720-200ppm and PLX4720/PLX2695-200/7ppm chow. We thank the members of the Bosenberg and Yan laboratories for their helpful discussions and support. We thank the whole KDM5 team from Emory Specialized Center and the National Center for Advancing Translational Sciences (NCATS) in the Chemical Biology Consortium and the University of Texas MD Anderson Cancer Center, especially Dr. Haian Fu, Dr. Xiaodong Cheng, Dr. Paula M. Vertino, Dr. John R. Horton, Ajit Jadhav, Dr. David J. Maloney, Dr. Margaret A. Johns, John R. Shanks, Dr. Philip J. Webber, and Dr. Melinda S. Hanes for their essential contribution to the selection and characterization of KDM5 inhibitors used in this article. We thank Drs. Anthony Welch and Barbara Mroczkowski (NCI) and Drs. Andrew Flint and Dane Liston (Frederick National Laboratory for Cancer Research) for administrative support.
These studies were supported in part by grants from by the NIH Yale SPORE in Skin Cancer, NCI P50 CA121974 (to M.W. Bosenberg and Q. Yan), Department of Defense Peer-Reviewed Cancer Research Program Award W81XWH-13-10235 (to Q. Yan) and W81XWH-16-1-0306 (to S.M. Zhang), the Melanoma Research Foundation (to M.W. Bosenberg and Q. Yan), R01 CA196660 (to M.W. Bosenberg), P01 CA128814 (to M.W. Bosenberg and D.P. Kelly), the Melanoma Research Alliance (M.W. Bosenberg), the American Skin Association (M.W. Bosenberg), the Hervey Family Foundation (M.W. Bosenberg), by pilot funding from the Yale Cancer Center, NCI P30 CA16359 (M.W. Bosenberg and Q. Yan), T32CA193200—Yale Cancer Biology Training Program (to M.K. McGeary), F30 CA228444 (to I. Krykbaeva), and by federal funds from the NCI, NIH, under Chemical Biology Consortium Contract No. HHSN261200800001E, and the intramural research program of the NCATS, NIH (to D.J. Jansen, S.C. Kales, A. Simeonov, and M.D. Hall).
The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.