- Open Access
Interleukin-12 elicits a non-canonical response in B16 melanoma cells to enhance survival
Cell Communication and Signaling volume 18, Article number: 78 (2020)
Oncogenesis rewires signaling networks to confer a fitness advantage to malignant cells. For instance, the B16F0 melanoma cell model creates a cytokine sink for Interleukin-12 (IL-12) to deprive neighboring cells of this important anti-tumor immune signal. While a cytokine sink provides an indirect fitness advantage, does IL-12 provide an intrinsic advantage to B16F0 cells?
Acute in vitro viability assays were used to compare the cytotoxic effect of imatinib on a melanoma cell line of spontaneous origin (B16F0) with a normal melanocyte cell line (Melan-A) in the presence of IL-12. The results were analyzed using a mathematical model coupled with a Markov Chain Monte Carlo approach to obtain a posterior distribution in the parameters that quantified the biological effect of imatinib and IL-12. Intracellular signaling responses to IL-12 were compared using flow cytometry in 2D6 cells, a cell model for canonical signaling, and B16F0 cells, where potential non-canonical signaling occurs. Bayes Factors were used to select among competing signaling mechanisms that were formulated as mathematical models. Analysis of single cell RNAseq data from human melanoma patients was used to explore generalizability.
Functionally, IL-12 enhanced the survival of B16F0 cells but not normal Melan-A melanocytes that were challenged with a cytotoxic agent. Interestingly, the ratio of IL-12 receptor components (IL12RB2:IL12RB1) was increased in B16F0 cells. A similar pattern was observed in human melanoma. To identify a mechanism, we assayed the phosphorylation of proteins involved in canonical IL-12 signaling, STAT4, and cell survival, Akt. In contrast to T cells that exhibited a canonical response to IL-12 by phosphorylating STAT4, IL-12 stimulation of B16F0 cells predominantly phosphorylated Akt. Mechanistically, the differential response in B16F0 cells is explained by both ligand-dependent and ligand-independent aspects to initiate PI3K-AKT signaling upon IL12RB2 homodimerization. Namely, IL-12 promotes IL12RB2 homodimerization with low affinity and IL12RB2 overexpression promotes homodimerization via molecular crowding on the plasma membrane.
The data suggest that B16F0 cells shifted the intracellular response to IL-12 from engaging immune surveillance to favoring cell survival. Identifying how signaling networks are rewired in model systems of spontaneous origin can inspire therapeutic strategies in humans.
Plain english summary
Interleukin-12 is a key cytokine that promotes anti-tumor immunity, as it is secreted by antigen presenting cells to activate Natural Killer cells and T cells present within the tumor microenvironment. Thinking of cancer as an evolutionary process implies that an immunosuppressive tumor microenvironment could arise during oncogenesis by interfering with endogenous anti-tumor immune signals, like IL-12. Previously, we found that B16F0 cells, a cell line derived from a spontaneous melanoma, interrupts this secreted heterocellular signal by sequestering IL-12, which provides an indirect fitness advantage. Normally, IL-12 signals via a receptor comprised of two components, IL12RB1 and IL12RB2, that are expressed in a 1:1 ratio and activates STAT4 as a downstream effector. Here, we report that B16F0 cells gain an intrinsic advantage by rewiring the canonical response to IL-12 to instead initiate PI3K-AKT signaling, which promotes cell survival. The data suggest a model where overexpressing one component of the IL-12 receptor, IL12RB2, enables melanoma cells to shift the functional response via both IL-12-mediated and molecular crowding-based IL12RB2 homodimerization. To explore the generalizability of these results, we also found that the expression of IL12RB2:IL12RB1 is similarly skewed in human melanoma based on transcriptional profiles of melanoma cells and tumor-infiltrating lymphocytes.
Responses of the different cell types within tissues, such as T cells and Langerhans cells within the epidermal layer of skin, are coordinated by the extracellular release of protein signals. One such extracellular signal is the cytokine Interleukin-12 (IL-12) . IL-12 is a heterodimer comprised of p35 and p40 subunits that is produced by antigen presenting cells to boost natural killer (NK) and cytotoxic T lymphocyte (CTL) activity and to polarize T helper cells towards a type 1 phenotype . Activated NK, CTLs and type 1 T helper cells play important roles in controlling viral infections and defending against malignancy. Towards this aim, local delivery of IL-12 to the tumor microenvironment has been shown to promote tumor regression in both mouse models and human melanoma [3–8] while knocking out the IL-12 pathway increases malignancy .
Inside the cell, these extracellular signals are processed through a series of protein-protein interactions that transmit this information from the cell membrane to the nucleus to orchestrate a cellular response. In the case of IL-12, these extracellular signals are transduced by the canonical Janus-Kinase (JAK)-Signal Transducer and Activator of Transcription (STAT) signaling pathway . The IL-12 receptor, composed of a β1 (IL12RB1) and β2 (IL12RB2) subunit, lacks intrinsic enzymatic activity and forms a complex with two JAKs, JAK2 and TYK2. Autophosphorylation of the receptor complex results in binding and activation of STAT4. Phosphorylated STAT4 dimerizes and then translocates to the nucleus to function as a transcription factor . As a transcription factor, STAT4 promotes the expression and release of Interferon-gamma (IFN- γ), which upregulates MHC class I antigen presentation . An increase in antigen presentation enables CTLs to target and kill virally infected and malignant cells present within the tissue .
While most of our understanding of how information is transmitted between cells is based on normal biology, cancer cells arise through a process of mutation and selection [15, 16]. This implies that cancer cells develop ways during oncogenesis to block these extracellular signals or rewire intracellular signaling pathways to gain a fitness advantage [17, 18]. Of note, B16F0 cells correspond to the parental cell line isolated from a spontaneous tumor that developed in C57BL/6J mice and are used as a syngeneic transplantable mouse model for human melanoma [19, 20]. In contrast to genetically engineered mouse models where genetic alterations that promote oncogenic transformation are invoked a priori , tumor cell lines derived from spontaneous tumors are likely to have a broader collection of mechanisms that promote oncogenic transformation. Following from this context, the mouse melanoma cell line B16F0 overexpresses one component of the IL-12 receptor, IL12RB2, to form a cytokine sink for IL-12, which can deprive neighboring immune cells of this important anti-tumor cytokine . As rewiring of intracellular cell signaling networks is another way that tumor cells can gain advantage, we asked here whether IL-12 provides an intrinsic advantage to B16F0 tumor cells.
B16F0 mouse melanoma (ATCC, Manassas, VA), Melan A mouse melanocytes (V. Hearing, National Cancer Institute, Bethesda, MD) and 2D6 cells (M. Grusby, Harvard University, Cambridge, MA) were cultured as described previously . All cell lines were revived from frozen stock and used within 10-15 passages that did not exceed a period of 6 months. Specifically, complete medium for B16F0 cells consisted of Dulbecco’s Modified Eagle Medium (DMEM, Corning, Manassas, VA, USA) supplemented with 10% heat-inactivated fetal bovine serum (FBS, Hyclone, Logan UT, USA), 1% penicillin-streptomycin (Pen-Strep, Hyclone), and L-Glutamine (2 mM final concentration, Medtech Inc., Herndon, VA, USA). Melan A culture medium was DMEM with 10% FBS, 20 mM HCl, 10mM HEPES, 1% Pen-Strep, and fresh 12-O-tetradecanoyl phorbol-13-acetate (200 nM) and phenylthiourea (200 μM). Cells were detached using either Trypsin-Versene or magnesium and calcium-free Dulbecco’s Phosphate-Buffered Saline (DPBS, Corning). 2D6 cells were maintained in RPMI 1640 supplemented with 10% FBS-HI, 1% Pen-Strep, 2 mM L-Glutamine, 50 mM HEPES, 49 mM b-mercaptoethanol (Sigma Chemical, St Louis, MO), 100 mM sodium pyruvate (Fisher Scientific, Pittsburgh, PA, USA), 1% MEM non-essential amino-acid solution (100x Fisher Scientific) and 6.7 ng/ml recombinant mouse IL12p70 (IL-12, eBioscience, San Diego, CA, USA).
Apoptosis and necrosis assay
Cells were plated in a 12-well tissue culture plate (BD Falcon, BD Biosciences, San Jose, CA) for 12 hours, then pre-treated with complete culture media with a final concentration of 40, 100, or 200 ng/ml IL-12. After 18 hours of pretreatment, cells were treated with 0, 5, 10, or 20 μM Imatinib (Cayman Chemical, Ann Arbor Michigan). A 0.1% dimethyl sulfoxide (DMSO) vehicle control in complete media was used as a negative control, while cells heated to 50 oC for 5-10 minutes or exposed to 10% DMSO served as positive controls. After 24 hours of incubation with imatinib or other treatment, cell viability was assessed using either a flow cytometry-based assay or an ATPlite assay (PerkinElmer, Shelton, CT), which was performed according to the manufacturer’s instructions. For the flow cytometry-based assay, cell media and floating cells were pooled with cells detached with DPBS without magnesium or calcium, centrifuged in 15 mL conical tubes at 250xg, and cells transferred to a 96 well round bottom plate (BD Falcon) for staining. Cells were then washed with DPBS with magnesium and calcium. After washing, cells were stained with an apoptosis and necrosis staining kit (ABD Serotech, Bio-Rad, Raleigh, NC) according to manufacturer’s instructions and analyzed by flow cytometry on a BD LSRFortessa. Data was analyzed in the BD FACSDiva 3.0 software and in Microsoft Excel. Statistical significance in the difference in viability upon IL-12 treatment for a given imatinib concentration was assessed using a two-tailed homoscedastic t-Test, where a p-value less than 0.05 was considered significant.
Cells were grown in a 12-well tissue culture plate if adherent (B16F0 and Melan-A) or in a 96-well round bottom plate if non-adherent (2D6) (BD Falcon). At each time point, adherent cells were detached from the plate with Trypsin-EDTA or DPBS without calcium/magnesium. For IL12RB1 and IL12RB2 copy number, cells were stained for 30 minutes with Yellow or Violet Live/Dead Fixable Stain (BD Biosciences), fixed with Lyse/Fix Buffer (BD Biosciences), then stained for PE-IL12RB1 (BD Biosciences) and AlexaFluor 488-IL12RB2 (R&D Systems, Minneapolis, MN) antibodies. Quantum Simply Cellular microspheres (Bangs Laboratories, Fishers, IN) were stained with IL12RB1 or IL12RB2 antibody within the same experiment to assess receptor copy number. Single stained controls were used to assess fluorescent spillover. To assay STAT4 and Akt activity, B16F0 and 2D6 cells were starved of FBS and IL12p70 for 12-24 hours and then treated for 1hr or 3hr with the indicated concentration of IL12p70 or vehicle. Positive control cells were treated with 10% FBS and 40 ng/mL IL12p70. Cells were stained using Violet Fixable Live/Dead (Invitrogen) to assess cell viability, fixed (BD Biosciences Lyse/Fix), permeabilized with Perm Buffer III (BD Biosciences), blocked with Mouse IgG (1:100, concentration, Jackson Labs, West Grove, PA), stained with FITC-IL12RB2 (R&D Systems), PE-pAKT and AlexaFluor 647-pSTAT4 (BD Biosciences) and stored at 4 oC in PBSaz before analysis. Cells were analyzed using the BD LSRFortessa flow cytometry analyzer, FCSExpress, and BD FACSDiva 3.0 software. Following acquiring 20,000 events on average, forward and side scatter were used to gate cell populations, where live cells were identified using Live-Dead staining and used for subsequent analysis.
CRISPR/Cas9 constructs, transfection and cloning validation
Two vectors were prepared by GenScript (Piscataway, NJ) with puromycin resistance genes and the flanking guide RNA (gRNA) sequences TCCGCATACGTTCACGTTCT and AGAATTTCCAAGATCGTTGA and transfected into wild-type B16F0 cells using Lipofectamine 2000 (Thermo Fisher Scientific, Waltham, MA). Transfection was confirmed via microscopy visually using a green fluorescent protein control. Derivatives of the B16F0 cell line were obtained from single clones isolated using puromycin selection and single cell plating. About 25% of clones exhibited some level of IL12RB2 knockdown, which was confirmed by flow cytometry.
B16F0 and 2D6 cells were grown in complete media in 6-well plates 24 hours before the beginning of the experiment. At 12 hours prior, cells were starved of FBS and IL-12. One hour before start of experiment, cells were treated with control media, PI3K inhibitor LY294002 (50 μM, Cell Signaling Technology, Danvers, MA), or with Akt inhibitor MK-2206 (2.4 ug/mL or 5 μMol/L, Cayman Chemical, Ann Arbor, MI). Cells were then treated with 100 ng/mL IL-12 at t=0 or stimulated with FBS as a positive control. After 30 minutes or 3 hours, cells were collected for cell lysates for western blot analysis for total Akt (1:1000, #9272, Cell Signaling Technology), pAkt S472/473 (1:1000, #9271, Cell Signaling Technology), pAkt Thr308 (1:1000, #9275, Cell Signaling Technology), and β-actin (1:1000, Sigma Aldrich, St. Louis, MO). The cells were lysed in radioimmunoprecipitation assay buffer (50 mM Tris-HCl, pH 8.0, 150 mM NaCl, 1% Nonidet P-40, 0.5% sodium deoxycholate, 0.1% SDS, 2 mM EDTA, 1 mM NaVO4, 10 mM NaF, and protease inhibitors), and the protein concentration in the lysates was determined by a spectrophotometer. Equal amounts of the lysates were subjected to SDS-polyacrylamide gel electrophoresis and the immunoblotting was performed with use of the primary antibodies listed above and secondary antibodies conjugated to horseradish peroxidase (Jackson ImmunoResearch, West Grove, PA). Proteins were visualized by chemiluminescence.
Data analysis and statistics
Single-cell gene expression (scRNA-seq) analysis of 31 human melanoma tumors were downloaded from the Gene Expression Omnibus entry GSE115978 . All data was analyzed in R (V3.5.1) using the ‘stats’ package (V3.5.1). A Pearson’s Chi-squared test was used to assess statistical differences between distributions.
The impact of IL-12 on decreasing the cytotoxic efficacy of imatinib was inferred from the data using a mathematical model:
where CIm is the concentration of imatinib, CIL12 is the concentration of IL-12, EC50 is the effective concentration required to elicit a 50% drop in viability, Kd is the equilibrium binding constant of IL-12 and set equal to 40 ng/mL, and α is a parameter that quantifies the impact of IL-12 in shifting the EC50 of imatinib. As described previously , a Markov Chain Monte Carlo approach based on a Metropolis-Hasting algorithm was used to obtain converged Markov Chains that represent samples from the distributions in the model parameters that were consistent with the data, that is posterior distributions in the model parameters. In brief, three Markov Chains were generated with over-disperse starting points. New steps in the chains were proposed using a normally distributed random number generator with a mean of zero and adjusted standard deviation such that the acceptance fraction was 0.2. The Gelman-Rubin potential scale reduction factor applied to the parameters was used to assess convergence of the chains to sampling of the posterior distribution. Performance and diagnostic metrics of the chains are shown as Supplementary Figures S2 and S4. The posterior distribution in parameter α was used to infer whether the effect was greater than zero and significant. A p-value (P(α<0|Model,data)) of less 0.05 was considered as significant, that is the data support that IL-12 provides a rescue effect to imatinib. The significance associated with differences in STAT4 and Akt phosphorylation was assessed using a two-tailed homoscedastic t-Test, where a difference with a p-value less than 0.05 was considered significant.
Mathematical modeling and model selection
Generally, the IL-12 signaling pathway transduces signals by activating Akt and STAT4 in three steps. Here, we postulate that the first step is one of three competing mechanisms that create an activated receptor complex: a canonical ligand-dependent mechanism, a non-canonical receptor-dependent mechanism, and a hybrid model comprised of both canonical and non-canonical mechanisms. The canonical ligand-dependent mechanism corresponds to the reversible binding of a ligand (IL12) to a receptor multi-protein complex (RC) resulting in a conformation change that increases kinase activity of the complex (ARC):
with the forward and reverse rate constants kf1 and kr1, respectively. Assuming that the total number of receptor complexes (TOTRC) are conserved, the amount of activated receptor in equilibrium with the ligand concentration is given by:
where KC is equal to kr1/kf1. Similarly, a non-canonical receptor-dependent mechanism represents the spontaneous change in activity of a receptor complex:
At equilibrium, the amount of activated receptor is given by:
where KN is equal to kr2/kf2. The hybrid model combines both mechanisms:
with the equilibrium relationships:
Based on the relative abundance of IL-12 receptor proteins, we assumed that the total receptor complex (TOTRC) is comprised primarily of IL12RB1 and IL12RB2 heterodimers in 2D6 cells while IL12RB2 homodimers are primarily present in B16F0 cells and correspond to the abundances measured by flow cytometry. The concentration of IL12 is an experimental variable and the parameters KC and KN were fit to the experimental data.
An activated receptor complex catalyzes the second step by phosphorylating (p-) the signaling intermediate (S), such as Akt or STAT4:
with the rate constant for phosphorylation, kA. While two different mechanisms activate the receptor complex (e.g., ARC1 versus ARC2), we assumed that the phosphorylation kinetics of the signaling intermediate are the same. Activation is coupled with a deactivation step:
which is catalyzed by a generic phosphatase (PASE) with a rate constant of kD and captures the decline in phosphorylation upon removal of the ligand. Assuming the total signaling intermediate is conserved (TOTS), the amount of phosphorylated signaling intermediate for a given abundance of activated receptor complex (ARC) at steady state is:
where KD is equal to kD·PASE/kA. TOTS and KD were free parameters used in calibrating the mathematical models to the data. The three competing models are expressed as:
Except for KC, the three mathematical models were fit separately to Akt and STAT4 phosphorylation (expressed in mean fluorescent intensity (MFI)) versus IL12RB2 data (expressed in copies) for each cell line stratified by IL12RB2 expression and observed under unstimulated and IL-12 stimulated conditions. The ligand binding constant, KC, was determined for each cell line, as extracellular binding affinity should not influence the specificity of intracellular signal transduction. As the number of experimental observations (Nobs≥69) exceeded the number of free parameters (Nparameters≤4), the model parameters were theoretically identifiable from the data. A summed squares of the error (SSE) metric was used as a goodness-of-fit for model i:
which quantifies the deviation between specific observations, Yu, and their respective prediction by model i: YMi. The maximum likelihood estimate for the parameters (ΘML) of the corresponding mathematical models were obtained at the minimum SSE. To quantify the strength of evidence favoring model i over model j, we used the Bayes Factor :
where P(Y|Mi,ΘML) is the likelihood for observing data (Y), given the mathematical model i and parameter values ΘML, and can be estimated by:
where Nobs is the number of observations of the experimental variable. To interpret the Bayes Factor, values for the Bij between 1 and 3 are considered weak, between 3 and 20 are positive, between 20 and 150 are strong, and >150 are very strong evidence favoring model i over model j .
To investigate whether IL12RB2 provides an intrinsic advantage to the B16F0 cell line, we first asked whether the subunits of the IL-12 receptor are similarly upregulated in Melan-A cells, a model of ‘normal’ melanocytes. The 2D6 T cells, a mouse type 1 T helper cell line, were used as a positive control for IL-12 receptor expression . IL12RB1 and IL12RB2 abundance was assayed by flow cytometry and calibrated to copy number using antibody calibration beads (Fig. 1).
In 2D6 cells, IL12RB1 and IL12RB2 were present in almost a 2:1 copy ratio, with 43,000 copies of IL12RB1 and 19,500 copies of IL12RB2 and consistent with previous findings . Interestingly, the melanoma-related cell lines stained brightly for IL12RB2. The copy number of IL12RB2 was 1.5 times higher in the melanoma cell line versus the ‘normal’ melanocyte line, where the B16F0 cell line expressed over 330,000 copies of IL12RB2 and Melan-A cells had 214,000 copies. In contrast, IL12RB1 abundance was barely detectable in both cell lines, as staining was similar to the unstained controls. Only the B16F0 cell line was quantified at 1000 copies, while the Melan-A cell line was below the detection level of 200 copies per cell. These data suggest that expression of components of the IL-12 receptor, IL12RB1 and IL12RB2, are differentially regulated in mouse cell models of melanoma compared with T cells. This differential regulation is manifested as a high copy number of IL12RB2 and an absence of IL12RB1, such that copies of the heterodimeric IL-12 receptor are limited. Of note, we previously reported that the apparent molecular weight of IL12RB2 is similar between B16F0 and 2D6 cells when assayed by Western blot using an antibody recognizing the C-terminus (see Fig. 1F in ). As the antibody used for flow cytometry recognizes the extracellular domain of IL12RB2 (N-terminus), the Western blot and flow cytometry results suggest that a similar isoform of IL12RB2 is expressed in B16F0 and 2D6 cells.
As B16F0 cells do not make IL-12 (see Fig. 2D and S2A in ), creating a cytokine sink by overexpressing IL12RB2 provides an indirect survival advantage to malignant melanocytes. We also wondered whether overexpression of IL12RB2 provides a direct survival advantage to malignant melanocytes. To test for an intrinsic survival advantage, we used an apoptosis and necrosis assay to assess cell viability during a chemotherapy challenge with imatinib in the B16F0 melanoma (Fig. 2) in combination with an IL-12 rescue. We used imatinib, a small molecule tyrosine kinase inhibitor that was originally developed to inhibit the BCR-ABL fusion protein in leukemia  but also inhibits c-Kit, PDGF receptors, and CSF1R and is being investigated as a possible chemotherapy for metastatic melanoma [28, 29]. To determine whether IL-12 could rescue the melanoma cells from the cytotoxic effects of chemotherapy, we incubated the cells with three different concentration of IL-12 (40, 100, or 200 ng/mL) for 18 hours, which correspond to concentrations of IL-12 in the serum of various mouse models at approximately 3 – 134 ng/mL . We then treated with three different concentrations of imatinib (5, 10, or 20 μM) for 24 hours, with 0.1% DMSO in DPBS serving as a vehicle control. Cell viability was assayed by flow cytometry, whereby cells were identified based on forward and side scatter characteristics. Gating on “cells”, the presence of Annexin-V staining on the cell surface and no PI staining indicates a cell undergoing apoptosis, while cells that are positive for PI and negative for Annexin-V indicates a cell that underwent necrotic cell death.
As summarized in Fig. 2, the viability of B16F0 cells was greater than 90% at baseline. Upon treatment with imatinib, the viability was decreased in a dose-dependent manner, as expected. At 24 hours, very few of the cells were observed during the early stages of apoptosis (Annexin-V+, PI-) and the majority of non-viable cells either underwent necrotic cell death (Annexin-V-, PI+) or were in the late stages of apoptotic cell death (Annexin-V+, PI+). The addition of IL-12 reduced the cytotoxic efficacy of imatinib such that the maximum shift in viability of cells pre-treated with IL-12 before imatinib challenge was about 15%. To quantify the effects of imatinib and IL-12 on the viability of the B16F0 cells, we analyzed the data using a mathematical model and inferred the model parameters using a Markov Chain Monte Carlo (MCMC) approach. In theory, the model is identifiable as there are more data points (n=16) than model parameters: EC50 and α. However given the non-linear nature of the model, the MCMC approach was used to assess whether the model parameters could be practically identified from the available data in addition to estimating parameter values. Traces from three independent chains initiated with different starting values and the corresponding Gelman-Rubin convergence diagnostics are shown in Figure S2. Using the converged segments of the chains, the MCMC results suggested that the effective concentration of imatinib (EC50) and the rescue effect of IL-12 could be practically identified from the data as they exhibited symmetric bounded distributions. The posterior distributions provided a median value of the EC50 of imatinib equal to 40 nM for imatinib and of α equal to 1 (Fig. 2d). Moreover, the posterior distributions indicated that the rescue effect of IL-12 from the cytotoxic actions of imatinib was significant (p = 0.0095 that α<0). Quantitatively, a saturating concentration of IL-12 would increase the EC50 for imatinib by 100% (maximum likelihood value of α=1), thereby shifting the imatinib dose-response curve to the right. We repeated this experiment using a plate-based viability assay (ATPlite) using both B16F0 and Melan-A cells (Fig. 2e). We used a calibration curve to confirm that these results were within the dynamic range of the assay (see Figure S3a). Of note, higher concentrations of imatinib resulted in near complete cell death, although the results were outside of the dynamic range and thus not quantitative (see Figure S3b). While the variability of the ATPlite assay was higher than the flow cytometry-based assay, IL-12 did not provide a rescue effect when Melan-A cells, an immortalized normal melanocyte cell line, were exposed to imatinib. Specifically, MCMC-based analysis of the ATPlite-generated data with Melan-A cells provided a posterior distribution in α with a maximum likelihood value of 0 (see bottom panel of Fig. 2e and Figure S4 for MCMC diagnostics). While cell viability was assayed after only 24 hours, the observed difference observed with B16F0 cells implies that, over a 1 week period, tumor load would increase by about 10-fold in the IL-12 pre-treated cells versus chemotherapy alone.
To confirm the role of IL12RB2 in the rescue effect of IL-12 to cytotoxic stress with imatinib, we knocked down IL12RB2 in the B16F0 cell line using a CRISPR-Cas9 construct (Fig. 3a). Generally, knockdown clones grew at similar rates, showed minor differences in color, and displayed similar morphology to the wild-type cell line. Knockdown of IL12RB2 was confirmed by flow cytometry. To compare the functional response, wild type (wt) and IL12RB2 knockdown (KD) B16F0 cells were treated with increasing concentrations of imatinib for 12 hours and pre-stimulation with 100 ng/mL IL-12 overnight was included as a second independent variable. Cell viability was assessed using an ATPlite assay. Consistent with data shown in Fig. 3, imatinib dose-dependently reduced viability in the wild-type B16F0 cell line that was prevented with pre-treatment with IL-12 (Fig. 3b). The rescue effect of IL-12 was most pronounced at the higher concentrations of imatinib. In comparison, imatinib also dose-dependently reduced viability of the IL12RB2 knockdown cells but IL-12 provided no benefit (Fig. 3c). Visually, we also noticed that IL-12 had no effect on preventing cell detachment or inhibiting growth in IL12RB2 KD cells. Collectively, the results suggested that IL-12 enhances the survival of B16F0 cells in the presence of cytotoxic stressors, like imatinib, and its effect is dependent on IL12RB2 expression.
Given the imbalance of IL-12 receptor expression in B16F0 cells, it was unclear what signals were transduced in B16F0 cells in response to IL-12 stimulation. Moreover, activating STAT4, which plays a role in Interferon- γ gene expression and not survival, is thought to be the canonical signaling response to IL-12. To investigate, we assayed the phosphorylation of STAT4, a canonical signal transducer, and Akt, a non-canonical signal transducer involved in promoting cell survival and therapeutic resistance , in response to IL-12 stimulation in both live B16F0 and 2D6 cells by flow cytometry (Fig. 4). The 2D6 T helper cell line was used as a model for canonical IL-12 signal transduction. Following 12-24 hours of preconditioning in media free of IL-12 (2D6 only) and serum (2D6 and B16), we stimulated B16F0 and 2D6 cells with increasing concentrations of IL-12 (0, 7, 100, and 300 ng/ml). Pre-conditioning ensured that activation of the IL-12 signaling pathway was due to IL-12, and not due to factors in the serum or basal IL-12. PBS was used as a negative vehicle control for both cell lines.
Only 2D6 cells phosphorylated STAT4 in response to IL-12 stimulation (p <0.0015 for 2D6 and p >0.05 for the B16), as expected. Basal STAT4 phosphorylation was minimal in both cell lines for the vehicle only-treated controls. In terms of Akt phosphorylation, the B16 melanoma showed a marked activation of Akt within 3 hours at all concentrations of IL-12 tested, with a progressive increase in pAkt with increasing concentrations of IL-12 treatment (from p <0.05 to p <0.0001). In contrast, 2D6 cells exhibited a less pronounced shift in phosphorylation of Akt above baseline in response to IL-12 treatment (p <0.0015). Pre-conditioning cells in the absence of serum was critical to observe differences in Akt phosphorylation.
To determine whether Akt activation in the B16 melanoma line was specific, we repeated the experiment and used a Western Blot assay to determine basal levels of Akt and test several inhibitors against the IL-12-induced effect on Akt phosphorylation (Fig. 4 panels e and f). To serve as a negative control and to determine whether the PI3K pathway may be cross-activated in response to IL-12 stimulation, we used a PI3K inhibitor LY294002. In addition, we used a selective pAkt inhibitor, MK-2206, to investigate whether the observed effect was solely through possible direct activation of Akt by the IL-12 receptor. Two different phosphorylation sites on Akt were assessed, threonine 308 (Thr308) and serine 472/473 (Ser472/473). We also measured total Akt, and β-actin was used as a loading control. We found no change in phosphorylation of Thr308 in response to IL-12 stimulation, whereas baseline phosphorylation was turned off by both LY294002 and MK-2206. However, IL-12 stimulation produced a marked increase in phosphorylation of the Ser472/473 site, and again, phosphorylation was blocked by both the PI3K and Akt-selective inhibitors.
Conventionally, signal transduction within the IL-12 pathway is comprised of three steps. The first step involves a change in conformation of the multi-protein receptor complex upon ligand binding that promotes receptor kinase activity. Upon gaining kinase activity, the receptor complex phosphorylates signaling intermediates, like STAT4. In the third step, phosphorylated signaling intermediates localize to the nucleus and promote gene transcription. Yet the different responses to IL-12 in B16F0 versus 2D6 cells suggest that the first step uses one of three competing mechanisms to create an activated receptor complex: a canonical ligand-dependent mechanism, a non-canonical receptor-dependent mechanism, and a hybrid model comprised of both canonical and non-canonical mechanisms. To select among these competing mechanisms, we used the data presented in Fig. 4a-d coupled with a Bayes Factor analysis. First we stratified the flow cytometry results into different subsets based on IL12RB2 expression (see Figure S1, panel a) and quantified the phosphorylation of Akt and STAT4 in each IL12RB2 subset under different treatment conditions and in each cell line (Fig. 5). We also confirmed that the observed increases in Akt and STAT4 phosphorylation with increasing IL12RB2 abundance were not attributed to fluorescent spillover (see Figure S1, panels b and c). Mathematical relationships were developed based on the three competing mechanisms and fit to the experimental data (see Methods). Bayes Factors were used to compare how well each competing mechanism captured the data and tabulated at the bottom of Fig. 5.
A Bayes Factor (Bij) is a pairwise comparator that provides a rational basis for selecting among competing mathematical models by quantifying the strength of evidence supporting one mathematical model (model i) over another mathematical model (model j). The strength of evidence, that is the probability of generating model predictions that are similar to the observed data (i.e., P(Y|Model,ΘML)), is estimated by the maximum likelihood, that is the sum of squared errors between the experimental observations and the corresponding model predictions obtained using optimal parameter values (ΘML). Since three models are considered, we have displayed each possible pairwise comparison as a matrix where model i corresponds to the row and model j corresponds to the column. Numerical entries in the matrix are colored black when there is very strong evidence favoring model i over model j, that is when Bij is greater than 150. The Bayes Factors suggest that the hybrid model best described Akt phosphorylation in B16F0 and STAT4 phosphorylation in 2D6 cells while the non-canonical model best described STAT4 phosphorylation in B16F0 cells and Akt phosphorylation in 2D6 cells.
In regressing the mathematical models to the data, the ligand binding constant, KC, was determined to be 0.320 and 60.7 ng/ml in 2D6 and B16F0 cells, respectively. These values suggest that IL-12 has a 190-times higher affinity for IL12RB1:IL12RB2 compared to a receptor comprised of IL12RB2 only, as similarly reported . Interestingly, the canonical ligand-dependent mechanism provided the worst explanatory power in three out of four combinations. Specifically, the canonical model was unable to capture the increasing phosphorylation of Akt in both B16F0 and 2D6 cells and of STAT4 in B16F0 cells when cells with increasing copies of IL12RB2 were left unstimulated. In contrast, the non-canonical receptor-dependent mechanism provided the best explanation for phosphorylation of STAT4 in B16F0 cells and of Akt in 2D6 cells, although the extent of phosphorylation was low in both cases. In B16F0 cells, the non-canonical model captured the central trend whereby Akt phosphorylation increased with increasing IL12RB2 abundance but failed to capture IL-12-dependent variation that increased with IL12RB2 copies. For STAT4 phosphorylation in 2D6 cells, the non-canonical model provided the worst explanatory power. The hybrid model provided the best explanation of Akt phosphorylation in B16F0 cells. Moreover, knocking down IL12RB2 to below 100,000 copies in B16F0 cells is sufficient to eliminate receptor-dependent and IL-12-induced phosphorylation of Akt. For STAT4 in 2D6 cells, the hybrid model slightly edged out the canonical model as it was able to capture the slight increase in STAT4 phosphorylation observed in unstimulated cells. The high variability of STAT4 phosphorylation among replicates in 2D6 with low IL12RB2 copies is attributed to the low number of flow cytometric events included in these subsets.
Given the general trend in 2D6 cells where STAT4 and Akt phosphorylation increased with increasing IL12RB2 abundance, we next asked whether the observed responses in 2D6 cells and B16F0 cells represent different regions within a common signaling mechanism. Signals transduced by a receptor depend on assembly of the receptor complex, which can be facilitated by ligand binding. However, molecular crowding of receptors on the plasma membrane can increase ligand-independent signaling by promoting complex assembly . To explore the effect of molecular crowding, we plotted STAT4 and Akt phosphorylation against the estimated surface density of IL12 receptor complex (IL12R, Fig. 6). For 2D6 cells, the receptor complex was assumed to be a 1:1 ratio of IL12RB1 to IL12RB2 and the complex copy number corresponded to copies of IL12RB2. For B16F0 cells, the receptor complex was assumed to be a homodimer of IL12RB2 and the complex copy number corresponded to half of the measured copies of IL12RB2. The available membrane surface area for 2D6 and B16F0 cells were estimated to be 104 and 344 μm2, respectively , and adjusted based on a normalized forward scatter area obtained by flow cytometry as a measure of cell size .
As expected, IL-12 stimulation increased the level of STAT4 phosphorylation in 2D6 cells, with the level almost constant with increasing IL12R density (Fig. 6a). Even in the absence of IL-12, STAT4 phosphorylation remained constant as the IL12R density increased, with the levels in B16F0 cells similar than those observed in unstimulated 2D6 cells. The STAT4 response suggests that a different signaling mechanism applies to STAT4 phosphorylation in 2D6 cells versus B16F0 cells. For Akt phosphorylation (Fig. 6b), 2D6 and B16F0 cells provided identical levels of activation for the lower densities of IL12R, which suggests a common signaling mechanism. Only at higher densities of IL12R were differences observed upon IL-12 simulation. To rule out confounding effects that total Akt and STAT4 may increase in parallel with IL12RB2, we stimulated B16F0 cells with serum-containing media for 12 hours and assayed Akt and STAT4 phosphorylation by flow. We make the assumption that serum induced a maximal activation of Akt and STAT4 such that the levels of phosphorylated Akt and STAT4 represent the total Akt and STAT4 expression (see Figure S4, panels a and b). In subgroups stratified by IL12RB2 expression, STAT4 phosphorylation was independent of IL12R complex density. In contrast, phosphorylated Akt seemed to increase slightly with IL12R complex density. From these data, we developed a linear relationship between total Akt and IL12RB2 density (Total Akt (MFI) = 0.352 * IL12RB2 (in copies/ μm2) + 382). The total Akt was then estimated for each IL12RB2 subgroup and used for the value of TOTS in equation (12). The ratio pS/TOTS still increased both as a function of IL12R complex density and IL12 (see Figure S4c), which is consistent with a hybrid model.
In comparing the responses in Fig. 6, the similar response in STAT4 for all three concentrations of IL-12 in 2D6 cells while Akt phosphorylation progressively increased with increasing concentrations of IL-12 suggests that the IL12 receptor complex in 2D6 cells (i.e., IL12RB1:IL12RB2) is more efficient in phosphorylating STAT4 than either of the complexes (i.e., IL12RB1:IL12RB2 or IL12RB2:IL12RB2) in phosphorylating Akt. More importantly, the overlapping response in Akt phosphorylation between 2D6 and B16F0 cells for similar densities of IL12R suggests that the non-canonical response observed in B16F0 cells can be simply attributed to an increase in IL12RB2 copies. In addition, a gain in Akt phosphorylation is a continuous function of IL12RB2 expression. Collectively, the results suggest that IL-12 provides a survival advantage to B16F0 cells and that IL-12 induced a different cellular response in B16F0 cells compared to 2D6 T cells, whereby phosphorylation of Akt, a non-canonical signal transducer that is associated with enhanced cell survival is favored over STAT4, which is the canonical signal transducer in the IL-12 pathway.
To extend our focus beyond mouse cell models, we also wondered whether similar patterns of gene expression are observed in human melanoma based on single cell RNA-seq analysis. In a study of 31 human melanoma tumors (Fig. 7), malignant melanocytes, CD4+ T cells, and NK cells expressed measurable mRNA reads for both IL12RB1 and IL12RB2 (20.7%, 21.8%, and 23.9%, respectively) while only 10.3% of CD8+ T cells had non-zero reads for both IL12RB1 and IL12RB2. The fraction of cells with non-zero reads for IL12RB1 and not IL12RB2 increased progressively in malignant melanocytes, NK cells, CD4+ T cells, and CD8+ T cells, where the proportion being 55.9%, 65.2%, 72.7%, and 81.6%, respectively. Conversely, the fraction of cells with non-zero reads for IL12RB2 and not IL12RB1 decreased from 5.2% in malignant melanocytes to 0.3% in CD8+ T cells (5.2, 2.2, 0.8, and 0.3). The distribution in cells with non-zero reads for only IL12RB1, only IL12RB2, or both were statistically different among cell types (Fisher exact Test p-value <0.0005). It follows then that the ratio of IL12RB2 to IL12RB1 expression was higher for malignant melanocytes compared to both CD4+ and CD8+ T cells (p-value = 1e-5 for both). Interestingly, CD8+ T cells seemed to predominantly have an IL12RB2:IL12RB1 ratio of less than 1 (median = 0.085). The distribution in the IL12RB2:IL12RB1 ratio for NK cells seemed similar to CD4+ T cells, although there were too few cells (n=22) to assess significance. Collectively, these data suggest that human melanoma cells have an increased expression of IL12RB2 relative to IL12RB1 that is different than observed in immune cells that canonically respond to IL-12, such as NK cells, CD4+ T cells, and CD8+ T cells and qualitatively mirrors what we observed with mouse cell lines. Emerging proteogenomic methods (e.g., [36, 37]) that can simultaneously measure protein and mRNA abundance at the single cell level would help strengthen the clinical relevance of these observations.
Given the importance of IL-12 in promoting type 1 anti-tumor immune response, we had previously observed that, by overexpressing one component of the IL-12 receptor IL12RB2, B16 melanoma cells could create a “cytokine sink" and deprive tumor-infiltrating lymphocytes of this local signal to sustain an anti-tumor response . Here, our results indicate that IL-12 provided an intrinsic survival advantage to B16 melanoma cells and acted directly on B16 melanoma cells to phosphorylate Akt (Fig. 4a and c), which is in contrast to the canonical phosphorylation of STAT4 in T helper cells (Fig. 4b and d). The survival benefit of IL-12 was mediated through IL12RB2, as knock down of IL12RB2 abrogated the effect. As the IL-12 rescue was not observed with Melan-A cells, the collective results suggest that melanoma cells rewire signal transduction to engage the non-canonical signal transducer Akt to promote survival in response to cytokines within the tumor environment.
In Fig. 8, we propose a mechanism by which rewiring of the canonical IL-12 pathway could occur that is initiated by homodimerization of IL12RB2 and is enhanced by IL-12 stimulation and molecular crowding. IL-12 has been shown to promote a low but significant level of homodimerization of IL12RB2 subunits . Overexpression of IL12RB2 makes this association more likely to occur due to molecular crowding, as an increased number of receptors on the cell membrane increases the probability that a given receptor can get close enough to another receptor to phosphorylate the cytoplasmic domain. Following IL12RB2 homodimerization, the receptor associated Janus kinase, JAK2, interacts to promote cross phosphorylation and activation . Activated JAK2 can then recruit and activate PI3K via the p85a subunit [40, 41]. While other intermediates may play a role, activated PI3K can then directly activate Thr308-primed Akt via phosphorylating the Ser 472/473 loci [42, 43]. Our findings support the direct role of PI3K in the observed phosphorylation of Akt to affect cell survival, as inhibitors directed against PI3K pathway like LY294002 ameliorated this effect. While the data suggests that phosphorylation of Akt can occur with the IL12RB1:IL12RB2 complex, the copies of IL12RB1:IL12RB2 in 2D6 cells are too low to have a functional impact.
Rewiring of signaling pathways during oncogenesis is an important concept in cancer biology . For instance, a recent study across all known human kinases showed there are six types of network-attacking mutations that perturb signaling networks to promote cancer, including mutations that dysregulate phosphorylation sites to cause constitutive activation or deactivation, rewire networks through upstream or downstream interactions, and dysregulate network dynamics by activated or deactivating nodes of response . As for precedence of non-canonical signaling within the JAK-STAT pathway, STAT3 has been reported as a non-canonical activator of the NF- κB pathway to promote myeloid-derived suppressor cells in the tumor microenvironment . Beside oncogenic fusion proteins, increased abundance of intracellular signaling proteins through either copy number amplifications or epigenetic changes can also create non-canonical protein interactions that could redirect the flow of intracellular information . Given that the Akt response between the 2D6 and B16F0 cell lines are the same where the IL12 receptor complex densities overlap, favorable changes in coding sequence are not required but merely changes that increase the abundance of IL12RB2 are sufficient to impart a fitness advantage. Moreover, the improvement in fitness, as suggested by increased Akt phosphorylation, is a continuous function of IL12RB2 copies. Knowing the basis for cancer cell rewiring can be used for therapeutic advantage. For instance, the non-canonical use of signaling cross-talk downstream of the EGF receptor sensitized resistant triple negative breast cancers to DNA-damaging chemotherapies by inhibition of EGFR . Similar approaches could be applied to circumvent rewiring of the IL-12 pathway.
In summary, IL-12 provides an important function in coordinating an anti-tumor immune response. In melanoma, overexpression of one component of the IL-12 receptor, IL12RB2, by malignant melanocytes may not only serve to deprive immune cells of this important cytokine, but malignant melanocytes may also use IL-12 as a means of activating cell survival pathways. Supported by prior studies, we present evidence showing that IL-12 increased the dose of imatinib chemotherapy needed to induce cell death in B16 melanoma cells and that IL-12 activated a non-canonical cell survival protein, Akt, in a melanoma cell line, but not in the 2D6 T helper cell line. The differential response in B16F0 cells is explained by a mechanism involving both ligand-dependent and ligand-independent aspects to initiate PI3K-AKT signaling upon IL12RB2 homodimerization. Namely, IL-12 promotes IL12RB2 homodimerization with low affinity and IL12RB2 overexpression promotes homodimerization via molecular crowding on the plasma membrane. While subsequent studies will help clarify the mechanistic details to improve the translational impact, these studies highlight the importance of non-canonical pathway activation in cancer.
Availability of data and materials
The single cell RNAseq datasets analyzed during the current study are available in the Gene Expression Omnibus entry GSE115978. Flow cytometry datasets generated during the current study are available from the corresponding author on reasonable request. All other data generated or analyzed during this study are included in this published article and its supplementary information files.
cytotoxic T lymphocyte
Dulbecco’s Modified Eagle Medium
fetal bovine serum
- IFN-γ :
Major Histocompatibility Complex
Markov Chain Monte Carlo
mean fluorescent intensity
phosphate buffered saline
Signal Transducer and Activator of Transcription
Colombo MP, Trinchieri G. Interleukin-12 in anti-tumor immunity and immunotherapy. Cytokine Growth Factor Rev. 2002; 13(2):155–68.
Athie-Morales V, Smits HH, Cantrell DA, Hilkens CMU. Sustained IL-12 Signaling Is Required for Th1 Development. J Immunol. 2004; 172(1):61–9. https://doi.org/10.4049/jimmunol.172.1.61.
Zhang L, Morgan RA, Beane JD, Zheng Z, Dudley ME, Kassim SH, Nahvi AV, Ngo LT, Sherry RM, Phan GQ, Hughes MS, Kammula US, Feldman SA, Toomey MA, Kerkar SP, Restifo NP, Yang JC, Rosenberg SA. Tumor-infiltrating lymphocytes genetically engineered with an inducible gene encoding interleukin-12 for the immunotherapy of metastatic melanoma. Clin Cancer Res. 2015; 21(10):2278–88. https://doi.org/10.1158/1078-0432.CCR-14-2085.
Zhao X, Bose A, Komita H, Taylor JL, Kawabe M, Chi N, Spokas L, Lowe DB, Goldbach C, Alber S, Watkins SC, Butterfield LH, Kalinski P, Kirkwood JM, Storkus WJ. Intratumoral IL-12 gene therapy results in the crosspriming of Tc1 cells reactive against tumor-associated stromal antigens. Mol Ther. 2011; 19(4):805–14. https://doi.org/10.1038/mt.2010.295.
Koneru M, O’Cearbhaill R, Pendharkar S, Spriggs DR, Brentjens RJ. A phase I clinical trial of adoptive T cell therapy using IL-12 secreting MUC-16(ecto) directed chimeric antigen receptors for recurrent ovarian cancer. J Trans Med. 2015; 13(1):102. https://doi.org/10.1186/s12967-015-0460-x.
Chinnasamy Z, Yu D, Kerkar SP, Zhang L, Morgan RA, Restifo NP, Rosenberg SA. Local delivery of interleukin-12 using T cells targeting VEGF receptor-2 eradicates multiple vascularized tumors in mice. Clin Cancer Res. 2012; 18(6):1672–83. https://doi.org/10.1158/1078-0432.CCR-11-3050.
Choi I-K, Lee J-S, Zhang S-N, Park J, Lee K-M, Sonn CH, Yun C-O. Oncolytic adenovirus co-expressing IL-12 and IL-18 improves tumor-specific immunity via differentiation of T cells expressing IL-12R β2 or IL-18R α. Gene Ther. 2011; 18(9):898–909. https://doi.org/10.1038/gt.2011.37.
Seo SH, Kim KS, Park SH, Suh YS, Kim SJ, Jeun S-S, Sung YC. The effects of mesenchymal stem cells injected via different routes on modified IL-12-mediated antitumor activity. Gene Ther. 2011; 18(5):488–95. https://doi.org/10.1038/gt.2010.170.
Airoldi I, Di Carlo E, Cocco C, Sorrentino C, Fais F, Cilli M, D’Antuono T, Colombo MP, Colombo MP, Pistoia V. Lack of Il12rb2 signaling predisposes to spontaneous autoimmunity and malignancy. Blood. 2005; 106(12):3846–53.
Murray PJ. The JAK-STAT signaling pathway: input and output integration. J Immunol. 2007; 178(5):2623–9. https://doi.org/10.4049/jimmunol.178.5.2623.
Watford WT, Hissong BD, Bream JH, Kanno Y, Muul L, O’Shea JJ. Signaling by IL-12 and IL-23 and the immunoregulatory roles of STAT4. Immunol Rev. 2004; 202(1):139–56. https://doi.org/10.1111/j.0105-2896.2004.00211.x.
Yamamoto K, Quelle FW, Thierfelder WE, Kreider BL, Gilbert DJ, Jenkins NA, Copeland NG, Silvennoinen O, Ihle JN. Stat4, a novel gamma interferon activation site-binding protein expressed in early myeloid differentiation. Mol Cell Biol. 1994; 14(7):4342–9. https://doi.org/10.1128/MCB.14.7.4342.
Gaczynska M, Rock KL, Goldberg AL. Gamma-interferon and expression of MHC genes regulate peptide hydrolysis by proteasomes. Nature. 1993; 365(6443):264–7. https://doi.org/10.1038/365264a0.
Restifo NP, Bacík I, Irvine KR, Yewdell JW, McCabe BJ, Anderson RW, Eisenlohr LC, Rosenberg SA, Bennink JR. Antigen processing in vivo and the elicitation of primary CTL responses. J Immunol. 1995; 154(9):4414–22.
Bignell GR, Greenman CD, Davies H, Butler AP, Edkins S, Andrews JM, Buck G, Chen L, Beare D, Latimer C, Widaa S, Hinton J, Fahey C, Fu B, Swamy S, Dalgliesh GL, Teh BT, Deloukas P, Yang F, Campbell PJ, Futreal PA, Stratton MR. Signatures of mutation and selection in the cancer genome. Nature. 2010; 463(7283):893–8. https://doi.org/10.1038/nature08768.
Klinke DJ. An evolutionary perspective on anti-tumor immunity. Front Oncol. 2012; 2:202.
Heidel FH, Arreba-Tutusaus P, Armstrong SA, Fischer T. Evolutionarily conserved signaling pathways: acting in the shadows of acute myelogenous leukemia’s genetic diversity. Clin Cancer Res. 2015; 21(2):240–8. https://doi.org/10.1158/1078-0432.CCR-14-1436.
Kvitek DJ, Sherlock G. Whole genome, whole population sequencing reveals that loss of signaling networks is the major adaptive strategy in a constant environment. PLoS Genet. 2013; 9(11):1003972. https://doi.org/10.1371/journal.pgen.1003972.
Almeida JL, Dakic A, Kindig K, Kone M, Letham DLD, Langdon S, Peat R, Holding-Pillai J, Hall EM, Ladd M, Shaffer MD, Berg H, Li J, Wigger G, Lund S, Steffen CR, Fransway BB, Geraghty B, Natoli M, Bauer B, Gollin SM, Lewis DW, Reid Y. Interlaboratory study to validate a STR profiling method for intraspecies identification of mouse cell lines. PLoS ONE. 2019; 14(6):0218412.
Overwijk WW, Restifo NP. B16 as a mouse model for human melanoma. Curr Protoc Immunol. 2001; Chapter 20:20–1.
DuPage M, Jacks T. Genetically engineered mouse models of cancer reveal new insights about the antitumor immune response. Curr Opin Immunol. 2013; 25(2):192–9.
Kulkarni YM, Chambers E, McGray AJR, Ware JS, Bramson JL, Klinke DJn. A quantitative systems approach to identify paracrine mechanisms that locally suppress immune response to Interleukin-12 in the B16 melanoma model. Integr Biol (Cambridge). 2012; 4(8):925–36. https://doi.org/10.1039/c2ib20053h.
Wu Y, Deng W, McGinley EC, Klinke DJ. Melanoma exosomes deliver a complex biological payload that upregulates PTPN11 to suppress T lymphocyte function. Pigment Cell Melanoma Res. 2017; 30(2):203–18.
Jerby-Arnon L, Shah P, Cuoco MS, Rodman C, Su MJ, Melms JC, Leeson R, Kanodia A, Mei S, Lin JR, Wang S, Rabasha B, Liu D, Zhang G, Margolais C, Ashenberg O, Ott PA, Buchbinder EI, Haq R, Hodi FS, Boland GM, Sullivan RJ, Frederick DT, Miao B, Moll T, Flaherty KT, Herlyn M, Jenkins RW, Thummalapalli R, Kowalczyk MS, Canadas I, Schilling B, Cartwright ANR, Luoma AM, Malu S, Hwu P, Bernatchez C, Forget MA, Barbie DA, Shalek AK, Tirosh I, Sorger PK, Wucherpfennig K, Van Allen EM, Schadendorf D, Johnson BE, Rotem A, Rozenblatt-Rosen O, Garraway LA, Yoon CH, Izar B, Regev A. A Cancer Cell Program Promotes T Cell Exclusion and Resistance to Checkpoint Blockade. Cell. 2018; 175(4):984–97.
Klinke DJ. In silico model-based inference: a contemporary approach for hypothesis testing in network biology. Biotechnol Prog. 2014; 30(6):1247–61.
Kass RE, Raftery AE. Bayes factors. J Am Stat Assoc. 1995; 90:773–95.
Druker BJ, Talpaz M, Resta DJ, Peng B, Buchdunger E, Ford JM, Lydon NB, Kantarjian H, Capdeville R, Ohno-Jones S, Sawyers CL. Efficacy and safety of a specific inhibitor of the BCR-ABL tyrosine kinase in chronic myeloid leukemia. N Engl J Med. 2001; 344(14):1031–7.
Carvajal RD, Antonescu CR, Wolchok JD, Chapman PB, Roman R-A, Teitcher J, Panageas KS, Busam KJ, Chmielowski B, Lutzky J, Pavlick AC, Fusco A, Cane L, Takebe N, Vemula S, Bouvier N, Bastian BC, Schwartz GK. KIT as a Therapeutic Target in Metastatic Melanoma. JAMA. 2011; 305(22):2327. https://doi.org/10.1001/jama.2011.746.
Hodi FS, Corless CL, Giobbie-Hurder A, Fletcher JA, Zhu M, Marino-Enriquez A, Friedlander P, Gonzalez R, Weber JS, Gajewski TF, O’Day SJ, Kim KB, Lawrence D, Flaherty KT, Luke JJ, Collichio FA, Ernstoff MS, Heinrich MC, Beadling C, Zukotynski KA, Yap JT, Van den Abbeele AD, Demetri GD, Fisher DE. Imatinib for melanomas harboring mutationally activated or amplified KIT arising on mucosal, acral, and chronically sun-damaged skin. J Clin Oncol. 2013; 31(26):3182–90. https://doi.org/10.1200/JCO.2012.47.7836.
Schulz O, Edwards AD, Schito M, Aliberti J, Manickasingham S, Sher A, Reis e Sousa C. CD40 triggering of heterodimeric IL-12 p70 production by dendritic cells in vivo requires a microbial priming signal. Immunity. 2000; 13(4):453–62.
Davies MA. The Role of the PI3K-AKT Pathway in Melanoma. Cancer J. 2012; 18(2):142–7. https://doi.org/10.1097/PPO.0b013e31824d448c.
Gately MK, Renzetti LM, Magram J, Stern AS, Adorini L, Gubler U, Presky DH. The interleukin-12/interleukin-12-receptor system: role in normal and pathologic immune responses. Annu Rev Immunol. 1998; 16:495–521. https://doi.org/10.1146/annurev.immunol.16.1.495.
Hsieh MY, Yang S, Raymond-Stinz MA, Steinberg S, Vlachos DG, Shu W, Wilson B, Edwards JS. Stochastic simulations of ErbB homo and heterodimerisation: potential impacts of receptor conformational state and spatial segregation. IET Syst Biol. 2008; 2(5):256–72.
Wang Q, Klinke DJ, Wang Z. CD8(+) T cell response to adenovirus vaccination and subsequent suppression of tumor growth: modeling, simulation and analysis. BMC Syst Biol. 2015; 9:27.
Tzur A, Moore JK, Jorgensen P, Shapiro HM, Kirschner MW. Optimizing optical flow cytometry for cell volume-based sorting and analysis. PLoS ONE. 2011; 6(1):16053.
Stoeckius M, Hafemeister C, Stephenson W, Houck-Loomis B, Chattopadhyay PK, Swerdlow H, Satija R, Smibert P. Simultaneous epitope and transcriptome measurement in single cells. Nat Methods. 2017; 14(9):865–8.
Lin J, Jordi C, Son M, Van Phan H, Drayman N, Abasiyanik MF, Vistain L, Tu HL, Tay S. Ultra-sensitive digital quantification of proteins and mRNA in single cells. Nat Commun. 2019; 10(1):3544.
Collison LW, Delgoffe GM, Guy CS, Vignali KM, Chaturvedi V, Fairweather D, Satoskar AR, Garcia KC, Hunter CA, Drake CG, Murray PJ, Vignali DAA. The composition and signaling of the IL-35 receptor are unconventional. Nat Immunol. 2012; 13(3):290–9. https://doi.org/10.1038/ni.2227.
Kershaw NJ, Murphy JM, Liau NPD, Varghese LN, Laktyushin A, Whitlock EL, Lucet IS, Nicola NA, Babon JJ. SOCS3 binds specific receptor-JAK complexes to control cytokine signaling by direct kinase inhibition. Nat Struct Mol Biol. 2013; 20(4):469–76. https://doi.org/10.1038/nsmb.2519.
Fuhrer DK, Yang YC. Complex formation of JAK2 with PP2A, P13K, and Yes in response to the hematopoietic cytokine interleukin-11. Biochem Biophys Res Commun. 1996; 224(2):289–96.
Yamada O, Ozaki K, Akiyama M, Kawauchi K. JAK-STAT and JAK-PI3K-mTORC1 pathways regulate telomerase transcriptionally and posttranslationally in ATL cells. Mol Cancer Ther. 2012; 11(5):1112–21. https://doi.org/10.1158/1535-7163.MCT-11-0850.
Tsuchiya A, Kanno T, Nishizaki T. PI3 kinase directly phosphorylates Akt1/2 at Ser473/474 in the insulin signal transduction pathway. J Endocrinol. 2014; 220(1):49–59.
Yoo JK, Cho JH, Lee SW, Sung YC. IL-12 provides proliferation and survival signals to murine CD4+ T cells through phosphatidylinositol 3-kinase/Akt signaling pathway. J Immunol. 2002; 169(7):3637–43.
Pawson T, Warner N. Oncogenic re-wiring of cellular signaling pathways. Oncogene. 2007; 26(9):1268–75. https://doi.org/10.1038/sj.onc.1210255.
Creixell P, Schoof EM, Simpson CD, Longden J, Miller CJ, Lou HJ, Perryman L, Cox TR, Zivanovic N, Palmeri A, Wesolowska-Andersen A, Helmer-Citterich Ma, Ferkinghoff-Borg J, Itamochi H, Bodenmiller B, Erler JT, Turk BE, Linding R. Kinome-wide decoding of network-attacking mutations rewiring cancer signaling. Cell. 2015; 163(1):202–17. https://doi.org/10.1016/j.cell.2015.08.056.
Yu J, Wang Y, Yan F, Zhang P, Li H, Zhao H, Yan C, Yan F, Ren X. Noncanonical NF- κB activation mediates STAT3-stimulated IDO upregulation in myeloid-derived suppressor cells in breast cancer. J Immunol. 2014; 193(5):2574–86. https://doi.org/10.4049/jimmunol.1400833.
Klinke DJ. Signal transduction networks in cancer: quantitative parameters influence network topology. Cancer Res. 2010; 70(5):1773–82.
Lee MJ, Ye AS, Gardino AK, Heijink AM, Sorger PK, MacBeath G, Yaffe MB. Sequential application of anticancer drugs enhances cell death by rewiring apoptotic signaling networks. Cell. 2012; 149(4):780–94. https://doi.org/10.1016/j.cell.2012.03.031.
The authors thank Cassidy Bland and Adam Palmer for assisting with manuscript preparation and Vanessa Cuppett and Ning Cheng for technical assistance.
Sources of funding that contributed to this work include NIH/NCI R01CA193473 (DJK), NSF 1644932 (DJK), WVU Flow Cytometry Core funding: MBRCC CoBRE grant GM103488/RR032138, ARIA S10 grant RR020866, Fortessa S10 grant OD016165, and WV InBRE grant GM103434. The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Cancer Institute, the National Institutes of Health, or the National Science Foundation. These funding agencies played no role in the design of the study; the collection, analysis, and interpretation of data; and the writing of the manuscript.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional file 1 Supplemental Figure S1: cell subsets were established that vary in expression of iL12RB2 and with minimal fluorescent spillover into pAkt and pSTAT4 measures. (a) Flow cytometric events were subdivided into groups based on IL12RB2 expression, where solid vertical lines denote limits of IL12RB2 expression associated with each subgroup. Red curve denotes distribution in IL12RB2 expression by B16F0 while gray shaded indicates unstained controls. Single-stained controls for IL12RB2 in B16F0 (b) and 2D6 (c) cells were used to establish that fluorescence associated with measuring Akt and STAT4 phosphorylation was not due to fluorescent spillover. Results for two independent experiments are shown for each subgroup based on IL12RB2 expression.
Additional file 2 Supplemental Figure S2: Diagnostics for Markov Chain Monte Carlo estimates of the posterior distribution in the Imatinib EC50 and IL12 effect (α) estimated from flow cytometry viability assay using B16F0 cells. The traces for the initial 10000 steps (A and B) and full length (C and D) of three independent chains (colored in black, blue, and red, respectively) with over disperse starting points. (E) The Gelman-Rubin potential scale reduction factor (PSRF) was used to assess convergence of the Markov Chains to the posterior distribution in the Imatinib EC50 (left panel) and IL12 effect α (right panel), where a value of less than 1.2 indicated that the chains have converged to sampling the posterior distribution. The MCMC chains converged after less than 50,000 steps. (F) New steps in the Markov Chain were proposed using a normally distributed random number generator with a mean of zero and adjusted standard deviation such that the acceptance fraction was 0.2.
Additional file 3 Supplemental Figure S3: Calibration curves for quantifying cell viability using the ATPlite assay. (A) Increasing concentrations of B16F0 cells were plated just prior to reading viability using the ATPlite assay to establish the dynamic range of the assay (left panel). Results for experimental conditions that were acquired within the dynamic range of the assay are indicated by green overlay. (B) In a separate dose-finding experiment, increasing concentrations of B16F0 cells were plated just prior to reading viability using the ATPlite assay (left panel). While higher doses of imatinib appeared to reduce cell viability to near zero, the experimental conditions were acquired outside of the dynamic range of the assay (green overlay).
Additional file 4 Supplemental Figure S4: Diagnostics for Markov Chain Monte Carlo estimates of the posterior distribution in the Imatinib EC50 and IL12 effect (α) estimated from ATPlite assay results using Melan-A cells. The traces for the initial 10000 steps (A and B) and full length (C and D) of three independent chains (colored in black, blue, and red, respectively) with over disperse starting points. (E) The Gelman-Rubin potential scale reduction factor (PSRF) was used to assess convergence of the Markov Chains to the posterior distribution in the Imatinib EC50 (left panel) and IL12 effect α (right panel), where a value of less than 1.2 indicated that the chains have converged to sampling the posterior distribution. The MCMC chains converged after less than 50,000 steps. (F) New steps in the Markov Chain were proposed using a normally distributed random number generator with a mean of zero and adjusted standard deviation such that the acceptance fraction was 0.2.
Additional file 5 Supplemental Figure S5: Estimating total Akt and STAT4 values. Media contained serum was used to elicit a near maximal phosphorylation of STAT4 (a) and Akt (b) in B16F0 cells following a 12 hour stimulation (red squares). Phosphorylation of STAT4 and Akt was assayed by flow cytometry, where results are shown for each subgroup based on IL12RB2 expression. Single-stained controls for IL12RB2 in B16F0 (black x’s) cells were used to establish that fluorescence associated with measuring Akt and STAT4 phosphorylation was not due to fluorescent spillover. From these data, we developed a linear relationship between total Akt and IL12RB2 density (Total Akt (MFI) = 0.352 * IL12RB2 (in copies/ μm2) + 382). The total Akt was then estimated for each IL12RB2 subgroup and used for the value of TOTS in equation (12). (c) The ratio of phosphorylated Akt to total Akt, which corresponds to pS/TOTS in equations 12, 14, 16, and 19, was calculated for the different experimental conditions. Total Akt was assumed to follow the same dependence on IL12RB2 in 2D6 and B16F0 cells.
About this article
Cite this article
Byrne-Hoffman, C.N., Deng, W., McGrath, O. et al. Interleukin-12 elicits a non-canonical response in B16 melanoma cells to enhance survival. Cell Commun Signal 18, 78 (2020). https://doi.org/10.1186/s12964-020-00547-4
- Intracellular signaling
- JAK-STAT signaling
- PI3K-AKT signaling