- Research article
- Open Access
- Open Peer Review
Genomics of NSCLC patients both affirm PD-L1 expression and predict their clinical responses to anti-PD-1 immunotherapy
BMC Cancer volume 18, Article number: 225 (2018)
The Correction to this article has been published in BMC Cancer 2018 18:413
Programmed Death Ligand 1 (PD-L1) is a co-stimulatory and immune checkpoint protein. PD-L1 expression in non-small cell lung cancers (NSCLC) is a hallmark of adaptive resistance and its expression is often used to predict the outcome of Programmed Death 1 (PD-1) and PD-L1 immunotherapy treatments. However, clinical benefits do not occur in all patients and new approaches are needed to assist in selecting patients for PD-1 or PD-L1 immunotherapies. Here, we hypothesized that patient tumor cell genomics influenced cell signaling and expression of PD-L1, chemokines, and immunosuppressive molecules and these profiles could be used to predict patient clinical responses.
We used a recent dataset from NSCLC patients treated with pembrolizumab. Deleterious gene mutational profiles in patient exomes were identified and annotated into a cancer network to create NSCLC patient-specific predictive computational simulation models. Validation checks were performed on the cancer network, simulation model predictions, and PD-1 match rates between patient-specific predicted and clinical responses.
Expression profiles of these 24 chemokines and immunosuppressive molecules were used to identify patients who would or would not respond to PD-1 immunotherapy. PD-L1 expression alone was not sufficient to predict which patients would or would not respond to PD-1 immunotherapy. Adding chemokine and immunosuppressive molecule expression profiles allowed patient models to achieve a greater than 85.0% predictive correlation among predicted and reported patient clinical responses.
Our results suggested that chemokine and immunosuppressive molecule expression profiles can be used to accurately predict clinical responses thus differentiating among patients who would and would not benefit from PD-1 or PD-L1 immunotherapies.
In clinical trials of PD-1 or PD-L1 checkpoint immunotherapies, patients with NSCLC separate into groups that respond or do not respond to immunotherapy treatment [1,2,3,4]. Objective responses for NSCLC in these studies range from 19.0–23.0%. Patients are selected for immunotherapy based on immunohistochemistry (IHC) detection of PD-L1 reactivity. Positive PD-L1 reactivity in tumors is considered to be important to predicting the success of PD-1 and PD-L1 immunotherapy treatments [2, 5]. However, IHC results to detect PD-L1 reactivity can vary depending upon different IHC platforms, differences in anti-PD-L1 antibodies, differences in scoring systems, and differences in positivity cut-off values [4, 6,7,8,9]. In the Blueprint PD-L1 IHC Assay Comparison Project , similar antibody-specific differences were seen. In all, this variability presents challenges for using PD-L1 reactivity as a sole marker for diagnosis and as a marker to predict the success of PD-1 and PD-L1 immunotherapy treatments.
More likely, there is a complex profile of molecules that contributes to the regulation of PD-L1 and to the subsequent immunosuppressive effects that NSCLC cells have on immune cells . Chae et al. suggested that reliable predictive molecules need to be identified that can be used to select patients who would benefit from immunotherapy, yet limit the exposure of patients who would not benefit or have adverse reactions . Multifaceted predictive biomarker systems have also been proposed that contain input on PD-L1 expression, tumor mutations, and the roles of inflammatory cells to identify patients that would respond or not respond to immunotherapy treatment [13, 14]. For better treatment outcomes, there is a recognized need to develop additional methods that can identify a profile of molecules that contributes to the regulation of PD-L1 expression and affirms IHC PD-L1 positive reactivity. One approach is to use the influence of patient cell genomics on tumor cell signaling to identify the downstream effects on PD-L1 expression.
In this study, we hypothesized that patient tumor cell genomics influences cell signaling and the expression of PD-L1, chemokines, and immunosuppressive molecules. We also hypothesized that these profiles can be used to predict patient clinical responses. Rizvi et al. assessed the mutational profiles that determined sensitivity to PD-L1 blockade from patients with NSCLC treated with pembrolizumab  and we used the Rizvi et al. dataset to test our hypothesis. We first assessed the effect of patient genomics on the expression profile of 24 molecules: PD-L1, 9 chemokines, and 14 immunosuppressive molecules. Differences among patient-specific models reflected the input that their deleterious gene mutation profiles had on modeled signaling pathways and the expression of PD-L1, chemokines, and immunosuppressive molecules. Second, we used the expression profiles of these 24 chemokines and immunosuppressive molecules to sort patients into those that would or would not respond to PD-1 immunotherapy. The 9 chemokines were used to generate an index to predict dendritic cell infiltration and PD-L1 and the 14 immunosuppressive molecules were selected as tumor-derived molecules with a long list of reported immunosuppressive functions (Additional file 1: Table S1). Our results suggest that patient-specific chemokine and immunosuppressive molecule expression profiles can be used to accurately predict clinical responses thus differentiate among patients who would or would not respond to PD-1 immunotherapy.
Patient clinical characteristics and mutation profiles
This was a retrospective study and patient data, clinical characteristics, and exome sequencing information for each of 34 patients were obtained directly from Supplement Table 3 of the Rizvi et al study. study . To maintain anonymity, a random string generator was used to create a new random, 6-character uppercase alpha numeric string for each patient. This blinded both the identities of the patients in this study and their link to the prior published dataset we modeled.
All patients had stage IV NSCLC and were treated at Memorial Sloan Kettering Cancer Center (n = 29) or the University of California at Los Angeles (n = 5) on protocol NCT01295827. All patients had consented to Institutional Review Board-approved protocols permitting tissue collection and sequencing by the co-authors in this study (Naiyer A. Rizvi and Timothy A. Chan). All patients initiated therapy in 2012–2013 and were treated at 10 mg/kg every 2–3 weeks. Five patients were treated at 2 mg/kg every 3 weeks. The overall response rate and progression-free survival were reported to be similar across dose and schedules. PD-L1 expression on NSCLC tumor cells and immune cells by IHC was reported and scored semi-quantitatively: ≥50.0% membranous staining was considered strong, 1–49.0% was considered weak, and < 1.0% was considered negative .
Exomes from each NSCLC patient were examined using FannsDB , FATHMM , Mutation Assessor , Polyphen , PROVEAN , and SIFT . Gene mutations deleterious to gene function were identified (Additional file 2: Table S2). For example, there were 1192 gene mutations listed for patient SA97V5 and 36 mutations were deleterious to gene function (Fig. 1).
A validated cancer network containing a database of proteins involved in cell signal transduction, metabolism, and epigenetics obtained from manual review of new and published research (Additional file 3: Figure S1) was used to create patient NSCLC-specific predictive computational simulation models. This approach modeled protein-protein interactions at each step in a signaling pathway using ordinary differential equations (ODE)  and to predict specific pathway output . Pathway protein-protein interactions at each specific node were modeled as Michaelis-Menten equations that contained the reaction, enzyme, initial concentrations of protein intermediate reactants, and parameters of the reaction like Ka, Km, kcat, Vmax, etc. ODE were solved at each step by the Radau method . To demonstrate this modeling approach, an annexure section of the PD-L1 pathway is illustrated showing the step-by-step details of the protein-protein interactions at each node in the pathway as an example of the modeling process that also occured in all of the other pathways (Additional file 4: Supplementary Materials and Methods and Supplement Table 3 of the Rizvi et al study). The cancer network and the schema for creating these simulation models, predicting molecule responses, and identifying those patients who would or would not respond to PD-L1 immunotherapy is shown in Fig. 1.
NSCLC models in the cancer network were created for each patient. At the initial step, models did not contain patient-specific deleterious gene mutation profiles and were simulated to reach a homeostatic steady state, which served as the control baseline for the molecules of interest. Then patient-specific deleterious gene mutation profiles were converted into a computational format and annotated into the NSCLC cancer network, simulated to induce the patient-specific cancer disease states, and used to predict the expression of PD-L1, chemokines, and immunosuppressive molecules. At the network level, mutations of oncogenes were represented as gain of function at the activity level and mutations of tumor suppressor genes were represented as a loss of function at the activity level unless explicit functionality of the mutation was known from published studies. Copy number variations such as amplifications and deletions were represented as over-expression or deletion of gene function at the expression level. The time required to achieve a patient-specific network varied depending upon on the complexity of the patient-specific deleterious gene mutation profile.
The modeled output contained the expression profiles of 24 molecules (e.g., PD-L1, 9 chemokines, and 14 immunosuppressive molecules). PD-L1 expression was reported as percent change calculated as ((D/C)-1)*100. C was the absolute value of the non-tumorigenic baseline control (μM) and D was the absolute value of PD-L1 obtained from the patient-specific cancer state network (μM) . CCL2 , CCL3 , CCL4 , CCL5 , CCL11 , CCL20 , and CX3CL1  expression were determined similarly. These chemokines are capable of trafficking dendritic cells into the tumor microenvironment. Individual chemokine percent expression values were given weightage and normalized to sum to 1. A dendritic cell infiltration index was then calculated to be the sum of each prediction % change * weightage (Additional file 6: Table S4). Finally, the expression of 14 immunosuppressive molecules thought to facilitate the ability of cancer cells to escape normal tumor surveillance was determined (Additional file 1: Table S1).
Patient-specific simulation model predictions were also assessed using Weka 3, a data mining software program in Java . Weka 3 contained machine learning algorithms for data pre-processing; data classification, regression, clustering, and association rules; and data visualization. Using the predicted responses in Table 1, several machine-learning algorithms were implemented to learn prediction models (Additional file 7: Table S5).
Clinical response projections
Differences among the expression of 14 molecules were used in a 3-step process to sort patients into those that would or would not respond to PD-1 immunotherapy (Fig. 2). Patients were sorted by their PD-L1 expression (Step 1), their dendritic cell infiltration index (Steps 2a and b), and their immunosuppressive molecule expression (Steps 3a and b).
Graphic representations of the simulation model networks for each patient-specific model and the underlying network relationships were created as previously described  to identify similarities in patient-specific signaling pathways and to identify the influence of the pathway intermediates altered by the patient deleterious gene mutation profiles.
Simulation model validations
A series of internal control check analyses were used to validate the cancer network input and output data. These control checks monitored a) the effects of select pathway molecule over-expression or knockdown on pathway predictions, b) the effects of select drugs on pathway predictions, and c) the effects of activation, regulation, and cross-talk interactions among pathway intermediates on pathway predictions.
A cross-validation approach was used to assess the match scores of the PD-1 predicted responses against the PD-1 clinical responses in the Rizvi et al. 2015 Discovery dataset vs. the Validation dataset . The datasets were then pooled and re-partitioned into two new Training and Test datasets. A similar cross-validation approach was then used to assess the match scores of the PD-1 predicted responses vs. the PD-1 clinical responses. Differences between the match rates of the PD-1 predicted responses and the PD-1 clinical responses were performed via chi-square test or Fisher’s exact as previously described . All statistical tests utilized a 0.05 level of significance.
In this retrospective study, we created 29 of 34 separate and patient-specific simulation models from the exome sequencing information for each of the 34 patients listed in Supplement Table 3 of the Rizvi et al study . In the Discovery dataset, 13 of 16 patients had sufficient information to create simulation models and patients RYRJFL, IYXPLI, and GOFKQI did not (Additional file 2: Table S2). In the Validation dataset, 16 of 18 patients had sufficient information to create computational simulation models and patients 6NLFT5 and 32I5VC did not (Additional file 2: Table S2). The 5 patients with insufficient information were omitted from this study since their deleterious gene mutation profiles lacked driver genes and we were unable to achieve an increase in tumor phenotypes of proliferation and viability with the subset of gene aberrations reported. The remaining patients in the Discovery dataset (n = 13) contained 5 clinical responders and 8 clinical non-responders and the patients in the Validation dataset (n = 16) contained 6 clinical responders and 10 clinical non-responders (Table 1). Objective responses to PD-1/PD-L1 immunotherapies are known to vary. For example, objective responses to PD-1 immunotherapies for NSCLC were reported to range from 19.0–21.0% and objective responses to PD-L1 immunotherapies for NSCLC were reported to range from 10.0–23.0% [1, 2]. Hence the proportion of more non-responders than responders in this sample size was representative of the responses previously reported in larger patient populations.
Patient molecule responses
Results reported for patient SA97V5 were used as an example to clearly illustrate the process of model creation, model prediction, and model validation.
Modeled PD-L1 expression ranged from − 8.3% (patient 67K46M) to 185.5% (patient M9GYO4) (Table 1). Patient SA97V5 had a PD-L1 expression value of 67.0%.
Modeled chemokine expression was used to create a dendritic cell infiltration index. This index was a weighted function of the percentage change of each of the 9 individual chemokines (Table 1) and ranged from 4.20 (patient 67K46M) to 79.85 (patient DFZLO2). Patient SA97V5 chemokine expression for CCL2 (28.7%), CCL3 (14.3%), CCL4 (27.2%), CCL5 (13.4%), CCL7 (38.0%), CCL11 (36.9%), CCL20 (30.6%), CX3CL (32.9%), and CXCL14 (− 3.3%); formula details; and calculations for creating the index value of 23.9% are shown in Additional file 6: Table S4.
Modeled expression profiles for 14 immunosuppressive molecules, including those for patient SA97V5, are also shown in Table 1.
Clinical response projections
The expression of PD-L1, dendritic cell infiltration index, and immunosuppressive molecules were used in a 3-step process to sort patients into those that would or would not respond to PD-1 immunotherapy (Fig. 2).
At step 1, 9 patients with PD-L1 expression below 29.0% were identified as PD-1 drug non-responders (Figs. 2 and 3a). The remaining 16 patients with PD-L1 expression equal to or greater than 29.0% proceeded to step 2. Patient SA97V5 had a PD-L1 expression value of 67.0% and proceeded to step 2.
At Step 2, 2 patients with dendritic cell infiltration index values below 20.0% were identified as non-responders (Figs. 2 and 3b) and 2 patients with index values greater than 60.0% were identified as PD-1 drug responders (Figs. 2 and 3b). One mismatch occurred at this step. Clinical non-responder GI7AGZ with a dendritic cell index of 64.9% was misidentified as a PD-1 responder (Fig. 3b). Twelve patients with index values greater than 20.0%, but less than 60.0% proceeded to step 3. Patient SA97V5 had a dendritic cell infiltration index of 23.9% and proceeded to step 3.
At Step 3, 4 patients with immunosuppressive molecule values higher than that of their PD-L1 expression with a margin of greater than 5.0% were considered to be non-responders (Fig. 2, Step 3a) and 8 patient-specific models with values lower than that of their PD-L1 expression with a margin of greater than 5.0% were considered to be responders (Fig. 2, Step 3b). Three mismatches occurred at this step. Clinical responder patient 2FCOH7 had an immunosuppressive molecule expression profile of a non-responder: vascular endothelial growth factor (VEGF), cytotoxic T-lymphocyte-associated protein 4 (CTLA4), ganglioside GM3 (GM3), and ganglioside GD2 (GD2) were all higher than that of PD-L1 with a margin of greater than 5.0%. Clinical non-responder patients F3FK2W and 6QFSVV had immunosuppressive molecule expression profiles of responders. Patient SA97V5 had all 14 immunosuppressive molecules below the threshold of PD-L1 and was identified as a PD-1 drug responder (Figs. 2 and 4a) and patient QIA43T had molecules TGFB1 and IL6 above the threshold of PD-L1 and was identified as a PD-1 drug non-responder (Figs. 2 and 4b).
Patient-specific model predictions in the Discovery and Validation datasets were also checked using Weka 3  and SMO support vector machine with a normalized polynomial kernel had the best performance (Additional file 7: Table S5). The relationship between PD-L1 expression and predicted TGFB1 expression using Weka 3 algorithms for all patients in the dataset is shown in Additional file 8: Figure S2 and similar trends were seen when comparing the PD-L1 expression level to the other 13 predicted molecules. Weka 3 correctly identified 24 out of 29 patients whereas the computational simulation models correctly identified 25 of 29 patients.
In the cross-validation analysis of the match scores between the Rizvi et al. 2015 Discovery and Validation datasets, there were no significant differences between the match scores of non-responders and responders in the PD-1 clinical response group (38.5% vs. 37.5%; p = 0.9577, Additional file 9: Table S6) and PD-1 predicted response group (30.8% vs. 56.3%; p = 0.2642, Additional file 9: Table S6). Even though the Discovery dataset had a higher match score rate among the PD-1 clinical response group and the PD-1 predicted response group than the Validation dataset (92.3% vs. 81.2%, respectively), there was no significant difference between the two datasets (p = 0.6059).
Similarly, in the cross-validation analysis of the match scores between the Training and Test datasets, there were no significant differences between the match scores of non-responders and responders in the PD-1 clinical response group (38.9% vs. 36.4%; p = 0.9999, Additional file 9: Table S6) and PD-1 predicted response group (44.4% vs. 45.5%; p = 0.9577, Additional file 9: Table S6). In the Training dataset, the PD-1 predicted responses had an 83.3% match score with the PD-1 clinical responses and in the Test dataset the PD-1 predicted responses had a 90.9% match score with the PD-1 clinical responses. Again, there was no significant difference between the two datasets (p = 0.9999).
Predicted pathway comparisons
Deleterious gene mutations in patients were mapped to unique and common signaling pathways involved in PD-L1 expression (Fig. 5). Common pathways were utilized among a number of patient-specific models. Mutations in patient C9TGAJ (kirsten rat sarcoma viral oncogene homolog, KRAS mutation), patient RDD2UW (KRAS mutation), patient M9GYO4 (mitogen-activated protein kinase kinase 2, MAP2K2 mutation), and patient DFZLO2 (mitogen-activated protein kinase kinase kinase 1, MAP3K1 mutation) altered the extracellular signal-regulated kinase (ERK) activation pathway. Mutations in patient P90A0O (B-Raf proto-oncogene 1, BRAF1 mutation and tumor protein p53, TP53 mutation) and patient L8MTGU (KRAS mutation) altered ERK activation and apoptotic pathways. Mutations in patient SA97V5 (breast cancer anti-estrogen resistance protein 1, BCAR1 mutation; ankyrin-2, ANK2 mutation; insulin receptor substrate 1, IRS1 mutation; and cAMP response element-binding binding protein, CREBBP mutation), patient 26YMUF (rapamycin-insensitive companion of TOR, RICTOR mutation and ERK mutation), and patient L6ADEL (TP53 mutation and TNF receptor-associated factor, TRAF3 mutation) all had non-KRAS and non-B-Raf proto-oncogene (BRAF) driven activation pathways of PD-1 drug responder status.
Models also reinforced the association between PD-L1 expression and the presence of a KRAS mutation (or high ERK activation). Thus if a patient responder was found to have a KRAS mutation or positive regulator around mitogen-activated protein kinase kinase (MEK) pathway, this may identify a means to regulate PD-L1 by the MEK mediated pathway. A KRAS/BRAF/MEK related mutation in the profile leads to stronger expression of PD-L1 in the profile. Patients 195P5D, J0T9TJ, and 6QFSVV had KRAS mutations but were non-responders to PD-L1 inhibitor indicating complex and additional factors and pathways driving PD-L1 expression and response to the checkpoint inhibitor.
KRAS mutations and PD-L1 expression
To further affirm the association between the presence of KRAS mutations or KRAS co-mutations and PD-L1 expression, 2 additional datasets [33, 34] were modelled beyond the Supplement Table 3 of the Rizvi et al study .
KRAS mutations in lung adenocarcinoma were reported to be associated with co-mutations in TP53. In modeled simulations of the Dong, et al. dataset , the KRAS+TP53 co-mutation (KP subgroup) was predicted to increase PD-L1 expression. The KRAS+TP53 co-mutation had higher levels of predicted PD-L1 expression than the KRAS mutation and TP53 mutation alone.
KRAS mutations in lung adenocarcinoma were also reported to be associated with co-mutations in STK11/LKB1 (the KL subgroup) . In modeled simulations of the Skoulidis, et al. dataset , KRAS+STK11+KEAP1 co-mutation was predicted to reduce PD-L1 expression. The KRAS+STK11+KEAP1 co-mutation had lower levels of predicted PD-L1 expression than the KRAS+TP53 co-mutation.
KRAS mutations in lung adenocarcinoma were reported to be associated with co-mutations in TP53 (KP subgroup) and CDKN2A/B . In modeled simulations of the Skoulidis, et al. dataset , KRAS+CDKN2A/B co-mutations (KC subgroup) were predicted to reduce PD-L1 expression. KRAS+CDKN2A/B co-mutation had lower levels of predicted PD-L1 expression than the KRAS+STK11+KEAP1 co-mutation and the KRAS+TP53 co-mutation.
In this retrospective study, we used a recent dataset from NSCLC patients treated with pembrolizumab and identified deleterious gene mutational profiles in patient exomes. We annotated the deleterious gene mutational profiles into a cancer network to create NSCLC patient-specific predictive computational simulation models. We used these models as a tool to identify and validate a profile of 24 chemokines and immunosuppressive molecules that could accurately affirm expression of PD-L1 and predict patient clinical responses to PD-1 immunotherapy. We found that patient tumor cell genomics influenced cell signaling and altered the expression of PD-L1, 9 chemokines, and 14 immunosuppressive molecules. We also found that expression profiles of these 24 chemokines and immunosuppressive molecules could be used to identify patients who would or would not respond to PD-1 immunotherapy. Adding chemokine and immunosuppressive molecule expression profiles to a predicted PD-L1 profile allowed models to achieve a greater than 85.0% predictive correlation among predicted and reported patient clinical responses. This differentiated patients who would and would not benefit from PD-1 or PD-L1 immunotherapies. To validate our results, we used retrospective correlation of our simulation models against patient genomic signatures and clinical outcome data that was available in the NSCLC cohort of the Rizvi et al. study . We also used Weka 3 to validate predictions determined using the predictive computational simulation models. The Weka 3 results were similar to that generated via machine-learning methods and the Chi-square test was used to show no differences among the match rate results in these datasets. It is important to note that expanding PD-L1 expression profiles to include 23 additional chemokine and immunosuppressive molecule expression responses allowed models to achieve a greater than 85.0% correlation among predicted and reported patient clinical responses.
The 24 molecules used in this study have immunosuppressive properties. The role of PD-L1 in tumor pathogenesis is well known. Increased expression of PD-L1 on tumor cells inhibits T-cell proliferation, reduces T-cell survival, inhibits cytokine release, and promotes T-cell apoptosis [3, 35,36,37,38]. This leads to T-cell exhaustion and adaptive tumor immunosuppression .
Cytokines also have a role. Cytokine scores were recently found to be associated with overall survival in CheckMate 017 and 057 (both for nivolumab and docetaxel treated patients) . In this present study, 9 chemokines were selected that chemoattractant dendritic cells [41,42,43]. In other studies, dendritic cells were present in NSCLC tumors [42, 44] and dendritic cell infiltration was reported to be an independent prognostic factor for NSCLC . The simulation models here captured the role of dendritic cells by predicting chemokine expression in the form of a functional dendritic cell index (Additional file 6: Table S4).
Fourteen molecules were included that have pleiotropic functions including immunosuppressive properties (Additional file 1: Table S1). IL6, can prevent dendritic cell maturation, prime tumor-specific T-cells via signal transducer and activator of transcription 3 (STAT3) signaling, inhibit NF-κB binding activity, and inhibit C-C chemokine receptor type 7 (CCR7) expression [45,46,47]. IL10 can impair dendritic cell function and protect tumor cells from cytotoxic T-cell-mediated cytotoxicity by downregulating transporter-associated with antigen processing (TAP)1 and TAP2 [48, 49]. TGFβ can alter immune surveillance of regulatory T-cells. It represses CTL-mediated tumor cytotoxicity by altering the expression of perforin, granzyme A, granzyme B, Fas ligand (FASL), and IFNγ [49,50,51]. VEGF is a marker of tumor invasion and metastasis and can inhibit maturation of dendritic cells [44, 52]. IDO, a tryptophan-metabolizing enzyme that limits tryptophan, inhibits the proliferation of lymphocytes, and contributes to peripheral immunologic tolerance [53,54,55,56]. Increased IDO production by cancer cells down regulates natural killer (NK) receptors and induces NK cell apoptosis. It induces cell cycle arrest, decreases activation, and increases apoptosis in cytotoxic T-cells. Tryptophan 2, 3-dioxygenase 2 (TDO2) inhibits tryptophan 2,3-dioxygenase . Prostaglandin E2 (PGE2) suppresses NK cell function through the E2 prostaglandin receptor 4 (EP4) . Lectin, galactoside-binding, soluble, 9 (LGALS9) mediates T-cell dysfunction and T-cell senescence . Cluster of Differentiation 47 (CD47) is a negative regulator of dendritic cells binding to signal regulatory protein (SIRP) on dendritic cells and directly repressing dendritic cell phagocytosis, maturation, and production of IFNγ . CTLA4 restrains the adaptive immune response of T-cells towards tumor-associated antigens [60,61,62]. Gangliosides GM3 and GD2 induce monocyte apoptosis and impair differentiation to dendritic cells .
Using a 24 molecule expression profile allowed computational models to achieve a greater than 85.0% predictive correlation. However, this list was not exclusive and incorporating additional molecules into patient NSCLC-specific expression profiles may have merit and improve computational model accuracy. These included Lymphocyte activation gene-3 (LAG-3), T cell immunoglobulin-3 (TIM-3), and T cell immunoglobulin and ITIM domain (TIGIT) co-inhibitory receptors . LAG-3 is a co-inhibitory receptor upregulated on activated CD4+ T cells, CD8+ T cells, and subsets of natural killer (NK) cells [64,65,66]. It impairs T cell proliferation and cytokine production and alters NK cell cytotoxicity and cytokine production. TIM-3 is a cell surface molecule expressed on IFNγ-producing CD4+ T helper 1 cells, CD8+ T cytotoxic 1 T cells, NK cells, monocytes, and dendritic cells . TIM-3 dampens the development of protective immunity and TIM-3 blockade improves cell function. In patients with NSCLC, co-blockade of the TIM-3 and PD-1 pathways suppresses tumor growth. TIGIT is another co-inhibitory receptor expressed on NK cells, T cells, and Treg cells . CD155, CD112, and TIGIT ligands suppress immune responses through CD155 on dendritic cells. TIGIT is thought to work with PD-1 and TIM-3 to attenuate T cell responses and promote T cell dysfunction.
After the NSCLC patient-specific predictive computational simulation models were created and the profiles of 24 chemokines and immunosuppressive molecules were predicted, we created a decision tree to identify patients who would or would not respond to PD-1 immunotherapy. Decision cutoffs were established at 29.0% PD-L1 expression (Step 1), < 20.0% dendritic cell infiltration (Step 2a), > 60.0% dendritic cell infiltration (Step 2b), and immunosuppressive molecule expression as < PD-L1 with a margin of greater than 5.0% (Step 3) (Fig. 2). The decision tree was robust and had built-in redundancy. Basing the PD-L1 drug responder status on 3 separate predicted criteria allowed a responder/non-responder not identified at one step to be identified at a later step. Also the thresholds were specific. At 29.0% PD-L1 expression (Step 1), 9 non-responder patients were identified. Decreasing the PD-L1 expression cutoff from 29.0% to 25.0% identified only 6 non-responders. Increasing the PD-L1 expression cutoff from 29.0% to 35.0% identified a number of false negatives and setting the PD-L1 expression cutoff at 35.0% identified up to 13 non-responder patients: the three additional patient responders L8MTGU, P90A0O, and 26YMUF would be falsely identified as non-responders.
A diversity of signaling pathways are reported to be involved in the expression and regulation of PD-L1 [67,68,69] and these pathways were observed in expression and regulation of PD-L1 in this study (Fig. 5). Responder patients had mutations around the rapidly accelerated fibrosarcoma (RAF)-rat sarcoma (RAS)-ERK pathway including KRAS/BRAF and MEK-related mutations that predicted the profiles to have stronger expression of PD-L1. However, the presence of KRAS cannot be the only criteria for predicting strong expression of PD-L1 and thus a likely PD-1 drug responder, since there were non-responder profiles that also had KRAS mutations. We observed that matched predicted and clinical non-responder patients 195P5D and J0T9TJ and mismatched predicted responder and clinical non-responder patient 6QFSVV all had KRAS mutations (Additional file 2: Table S2, Fig. 5).
Recent studies support the concept that NSCLC is not a homogeneous disease and at least 3 subtypes of KRAS mutations involving LKB1 or TP53 can be identified. The tumors with these mutations have different PD-L1 expression patterns (higher in KRAS mutations and TP53 mutations) and different sensitivities to immune checkpoint blockade. Thus the effects of KRAS mutations and KRAS co-mutations on PD-L1 expression was further assessed using 2 additional datasets [33, 34] beyond the Supplement Table 3 of the Rizvi et al study .
Dong, et al.  reported that TP53 and KRAS mutations may predict which patients would or would not respond to PD-1 immunotherapy. Modeling the dataset in their study, we predicted that KRAS+TP53 co-mutation (KP Subgroup) would lead to increased PD-L1 expression. Skoulidis, et al.  reported that KRAS mutations in lung adenocarcinoma were associated with co-mutations in STK11/LKB1 (the KL subgroup) . KL tumors had high rates of KEAP1 mutations with lower PD-L1 expression. Modeling the dataset in their study, we predicted that KRAS+STK11+KEAP1 co-mutations (KL Subgroup) also would lead to reduced PD-L1 expression. We predicted that KRAS+CDKN2A/B co-mutation (KC Subgroup) would lead to reduced PD-L1 expression. There was a reduction in positive regulation due to reduction in AMPK, mTOR pathway and also due to KEAP1 loss of function. There was an increase in the WT TP53 mediated inhibitory regulation of PD-L1 expression. This was a novel finding based on network analysis.
The techniques described in this retrospective study have application. Although the techniques were complicated and need more extensive validation with larger datasets, their utility in clinical practice is possible. Profiling of tumors is becoming more main stream for precision personalized medicine. The approach may not necessarily be expensive, but in fact provides more utility to the generated profiling data for most tumor samples.
Patient tumor cell genomics were found to influence cell signaling with downstream effects on the expression of 24 chemokines and immunosuppressive molecules. This allowed us to establish patient-specific profiles of these molecules that could be used to predict patient clinical responses with greater than 85.0% correlation among predicted and reported patient clinical responses. Developing a workflow incorporating immunosuppressive molecules could a) be used as a potential complementary assay to affirm IHC results or used as an alternate assay where IHC in unfeasible, b) affirm patient PD-1 and PD-L1 drug responder status, c) as a method to determine influencing factors on PD-L1 expression, and d) as a potential clinical decision support system facilitating selection of therapies based on individual patient mutational profiles. The latter application used shortly after cancer diagnosis and just before cancer treatment could generate important patient-specific treatment options that could assist clinicians in selecting appropriate mono-therapies or combination therapies.
V-akt murine thymoma viral oncogene homolog
Activator protein 1
Breast cancer anti-estrogen resistance protein 1
B-Raf proto-oncogene 1
B-Raf proto-oncogene, serine/threonine kinase
C-C chemokine receptor type 7
Cluster of Differentiation 47
Cyclin-dependent kinase Inhibitor 2A
cAMP response element-binding binding protein
Cytotoxic T-lymphocyte-associated protein 4
E2 prostaglandin receptor 4
Extracellular signal-regulated kinase
Interferon gamma receptor 1
Interferon regulatory factor 1
Insulin receptor substrate 1
Kelch-like ECH-associated protein 1
Kirsten rat sarcoma viral oncogene homolog
Lymphocyte-activation protein 3
Lectin, galactoside-binding, soluble, 9
Liver kinase B1
- MAP K1:
Mitogen-activated protein kinase 1
- MAP K3:
Mitogen-activated protein kinase 3
Mitogen-activated protein kinase kinase 1
Mitogen-activated protein kinase kinase 2
Mitogen-activated protein kinase kinase kinase 1
Mitogen-activated protein kinase kinase
Mitogen-activated protein kinase kinase 1/2
Mechanistic target of rapamycin
MYC proto-oncogene, bHLH transcription factor
NK2 Homeobox 1
Neuroblastoma RAS viral oncogene homolog
Non-small cell lung cancer(s)
Ordinary differential equation(s)
Phosphatidylinositol-4,5-bisphosphate 3-kinase catalytic subunit alpha
Phosphatase and tensin homolog
Rapidly accelerated fibrosarcoma
Rapamycin-insensitive companion of TOR
Signal transducer and activator of transcription 1
Signal transducer and activator of transcription 3
Serine/threonine kinase 11
Transporter-associated with antigen processing
Tryptophan 2, 3-dioxygenase 2
T-Cell Immunoreceptor With Ig And ITIM Domains
T-cell immunoglobulin and mucin-domain containing-3
Tumor protein p53
TNF receptor-associated factor
Vascular endothelial growth factor
Shih K, Arkenau HT, Infante JR. Clinical impact of checkpoint inhibitors as novel cancer therapies. Drugs. 2014;74(17):1993–2013.
Sunshine J, Taube JM. PD-1/PD-L1 inhibitors. Curr Opin Pharmacol. 2015;23:32–8.
He J, Hu Y, Hu M, Li B. Development of PD-1/PD-L1 pathway in tumor immune microenvironment and treatment for non-small cell lung cancer. Sci Rep. 2015;5:13110.
Shukuya T, Carbone DP. Predictive markers for the efficacy of anti-PD-1/PD-L1 antibodies in lung cancer. J Thorac Oncol. 2016;11(7):976–88.
Carbognin L, Pilotto S, Milella M, Vaccaro V, Brunelli M, Calio A, Cuppone F, Sperduti I, Giannarelli D, Chilosi M, et al. Differential activity of Nivolumab, Pembrolizumab and MPDL3280A according to the tumor expression of programmed death-ligand-1 (PD-L1): sensitivity analysis of trials in melanoma, lung and genitourinary cancers. PLoS One. 2015;10(6):e0130142.
McLaughlin J, Han G, Schalper KA, Carvajal-Hausdorf D, Pelekanou V, Rehman J, Velcheti V, Herbst R, LoRusso P, Rimm DL. Quantitative assessment of the heterogeneity of PD-L1 expression in non-small-cell lung cancer. JAMA Oncology. 2016;2(1):46–54.
Kerr KM, Hirsch FR. Programmed death ligand 1 immunohistochemistry: friend or foe? Arch Pathol Lab Med. 2016;140(4):326–31.
Yu H, Boyle TA, Zhou C, Rimm DL, Hirsch FR. PD-L1 expression in lung cancer. J Thorac Oncol. 2016;11(7):964–75.
Brunnstrom H, Johansson A, Westbom-Fremer S, Backman M, Djureinovic D, Patthey A, Isaksson-Mettavainio M, Gulyas M, Micke P. PD-L1 immunohistochemistry in clinical diagnostics of lung cancer: inter-pathologist variability is higher than assay variability. Mod Pathol. 2017;30(10):1411–21.
Hirsch FR, McElhinny A, Stanforth D, Ranger-Moore J, Jansson M, Kulangara K, Richardson W, Towne P, Hanks D, Vennapusa B, et al. PD-L1 immunohistochemistry assays for lung cancer: results from phase 1 of the blueprint PD-L1 IHC assay comparison project. J Thorac Oncol. 2017;12(2):208–22.
Herbst RS, Soria JC, Kowanetz M, Fine GD, Hamid O, Gordon MS, Sosman JA, McDermott DF, Powderly JD, Gettinger SN, et al. Predictive correlates of response to the anti-PD-L1 antibody MPDL3280A in cancer patients. Nature. 2014;515(7528):563–7.
Chae YK, Pan A, Davis AA, Raparia K, Mohindra NA, Matsangou M, Giles FJ. Biomarkers for PD-1/PD-L1 blockade therapy in non-small-cell lung cancer: is PD-L1 expression a good marker for patient selection? Clinical Lung Cancer. 2016;17(5):350–61.
Bhaijee F, Anders RA. PD-L1 expression as a predictive biomarker: is absence of proof the same as proof of absence? JAMA Oncology. 2016;2(1):54–5.
Blank CU, Haanen JB, Ribas A, Schumacher TN. The “cancer immunogram”. Science. 2016;352(6286):658–60.
Rizvi NA, Hellmann MD, Snyder A, Kvistborg P, Makarov V, Havel JJ, Lee W, Yuan J, Wong P, Ho TS, et al. Mutational landscape determines sensitivity to PD-1 blockade in non-small cell lung cancer. Science. 2015;348(6230):124–8.
Martelotto LG, Ng CK, De Filippo MR, Zhang Y, Piscuoglio S, Lim RS, Shen R, Norton L, Reis-Filho JS, Weigelt B. Benchmarking mutation effect prediction algorithms using functionally validated cancer-related missense mutations. Genome Biol. 2014;15(10):484.
Reva B, Antipin Y, Sander C. Predicting the functional impact of protein mutations: application to cancer genomics. Nucleic Acids Res. 2011;39(17):e118.
Adzhubei IA, Schmidt S, Peshkin L, Ramensky VE, Gerasimova A, Bork P, Kondrashov AS, Sunyaev SR. A method and server for predicting damaging missense mutations. Nat Methods. 2010;7(4):248–9.
Choi Y, Sims GE, Murphy S, Miller JR, Chan AP. Predicting the functional effect of amino acid substitutions and indels. PLoS One. 2012;7(10):e46688.
Sim NL, Kumar P, Hu J, Henikoff S, Schneider G, Ng PC: SIFT web server: predicting effects of amino acid substitutions on proteins. In: Nucleic Acids Res. vol. 40 (Web Server issue); 2012: W452–7. doi: https://doi.org/10.1093/nar/gks539.
Kholodenko BN, Demin OV, Moehren G, Hoek JB. Quantification of short term signaling by the epidermal growth factor receptor. J Biol Chem. 1999;274(42):30169–81.
Kaushik P, Gorin F, Vali S. Dynamics of tyrosine hydroxylase mediated regulation of dopamine synthesis. J Comput Neurosci. 2007;22(2):147–60.
Hairer E, Wanner G. Stiff differential equations solved by Radau methods. J Comput Appl Math. 1999;111(1–2):93–111.
Lanzel EA, Hernandez MPG, Bates AM, Treinen CN, Starman EE, Fischer CL, Parashar D, Guthmiller JM, Johnson GK, Abbasi T, et al. Predicting PD-L1 expression on human cancer cells using next-generation sequencing information in computational simulation models. Cancer Immunol Immunother. 2016;65(12):1511–22.
Sagar D, Lamontagne A, Foss CA, Khan ZK, Pomper MG, Jain P. Dendritic cell CNS recruitment correlates with disease severity in EAE via CCL2 chemotaxis at the blood-brain barrier through paracellular transmigration and ERK activation. J Neuroinflammation. 2012;9:245.
Charmoy M, Brunner-Agten S, Aebischer D, Auderset F, Launois P, Milon G, Proudfoot AE, Tacchini-Cottier F. Neutrophil-derived CCL3 is essential for the rapid recruitment of dendritic cells to the site of Leishmania major inoculation in resistant mice. PLoS Pathog. 2010;6(2):e1000755.
Nobile C, Lind M, Miro F, Chemin K, Tourret M, Occhipinti G, Dogniaux S, Amigorena S, Hivroz C. Cognate CD4+ T-cell-dendritic cell interactions induce migration of immature dendritic cells through dissolution of their podosomes. Blood. 2008;111(7):3579–90.
Chabot V, Reverdiau P, Iochmann S, Rico A, Senecal D, Goupille C, Sizaret PY, Sensebe L. CCL5-enhanced human immature dendritic cell migration through the basement membrane in vitro depends on matrix metalloproteinase-9. J Leukoc Biol. 2006;79(4):767–78.
Beaulieu S, Robbiani DF, Du X, Rodrigues E, Ignatius R, Wei Y, Ponath P, Young JW, Pope M, Steinman RM, et al. Expression of a functional eotaxin (CC chemokine ligand 11) receptor CCR3 by human dendritic cells. J Immunol. 2002;169(6):2925–36.
Demedts IK, Bracke KR, Van Pottelberge G, Testelmans D, Verleden GM, Vermassen FE, Joos GF, Brusselle GG. Accumulation of dendritic cells and increased CCL20 levels in the airways of patients with chronic obstructive pulmonary disease. Am J Respir Crit Care Med. 2007;175(10):998–1005.
Johnson LA, Jackson DG. The chemokine CX3CL1 promotes trafficking of dendritic cells through inflamed lymphatics. J Cell Sci. 2013;126(Pt 22):5259–70.
Hall M, Frank E, Holmes G, Pfahringer B, Reutemann P, Witten IH: The WEKA data mining software: an update. In: ACM SIGKDD Explorations Newsletter. 2009;11(1):10–18. https://dl.acm.org/citation.cfm?id=1656278
Dong ZY, Zhong WZ, Zhang XC, Su J, Xie Z, Liu SY, Tu HY, Chen HJ, Sun YL, Zhou Q, et al. Potential predictive value of TP53 and KRAS mutation status for response to PD-1 blockade immunotherapy in lung adenocarcinoma. Clin Cancer Res. 2017;23(12):3012–24.
Skoulidis F, Byers LA, Diao L, Papadimitrakopoulou VA, Tong P, Izzo J, Behrens C, Kadara H, Parra ER, Canales JR, et al. Co-occurring genomic alterations define major subsets of KRAS-mutant lung adenocarcinoma with distinct biology, immune profiles, and therapeutic vulnerabilities. Cancer Discov. 2015;5(8):860–77.
Zitvogel L, Kroemer G. Targeting PD-1/PD-L1 interactions for cancer immunotherapy. OncoImmunology. 2012;1(8):1223–5.
Malm IJ, Bruno TC, Fu J, Zeng Q, Taube JM, Westra W, Pardoll D, Drake CG, Kim YJ. Expression profile and in vitro blockade of programmed death-1 in human papillomavirus-negative head and neck squamous cell carcinoma. Head Neck. 2015;37(8):1088–95.
Dong H, Strome SE, Salomao DR, Tamura H, Hirano F, Flies DB, Roche PC, Lu J, Zhu G, Tamada K, et al. Tumor-associated B7-H1 promotes T-cell apoptosis: a potential mechanism of immune evasion. Nat Med. 2002;8(8):793–800.
Tseng SY, Otsuji M, Gorski K, Huang X, Slansky JE, Pai SI, Shalabi A, Shin T, Pardoll DM, Tsuchiya H. B7-DC, a new dendritic cell molecule with potent costimulatory properties for T cells. J Exp Med. 2001;193(7):839–46.
Jiang Y, Li Y, Zhu B. T-cell exhaustion in the tumor microenvironment. Cell Death Dis. 2015;6:e1792.
Borghaei H, Brahmer JR, Horn L, Ready N, Steins M, Felip E, Paz-Ares LG, Arrieta O, Barlesi F, Antonia SJ, et al. Nivolumab (nivo) vs docetaxel (doc) in patients (pts) with advanced NSCLC: CheckMate 017/057 2-y update and exploratory cytokine profile analyses. J Clin Oncol. 2016;34(15_suppl):9025.
Balkwill F. Cancer and the chemokine network. Nat Rev Cancer. 2004;4(7):540–50.
Perrot I, Blanchard D, Freymond N, Isaac S, Guibert B, Pacheco Y, Lebecque S. Dendritic cells infiltrating human non-small cell lung cancer are blocked at immature stage. J Immunol. 2007;178(5):2763–9.
Mukaida N, Sasaki S, Baba T. Chemokines in cancer development and progression and their potential as targeting molecules for cancer treatment. Mediat Inflamm. 2014;2014:170381.
Inoshima N, Nakanishi Y, Minami T, Izumi M, Takayama K, Yoshino I, Hara N. The influence of dendritic cell infiltration and vascular endothelial growth factor expression on the prognosis of non-small cell lung cancer. Clin Cancer Res. 2002;8(11):3480–6.
Menetrier-Caux C, Montmain G, Dieu MC, Bain C, Favrot MC, Caux C, Blay JY. Inhibition of the differentiation of dendritic cells from CD34(+) progenitors by tumor cells: role of interleukin-6 and macrophage colony-stimulating factor. Blood. 1998;92(12):4778–91.
Hegde S, Pahne J, Smola-Hess S. Novel immunosuppressive properties of interleukin-6 in dendritic cells: inhibition of NF-kappaB binding activity and CCR7 expression. FASEB J. 2004;18(12):1439–41.
Ishii H, Tanaka S, Masuyama K. Therapeutic strategy for cancer immunotherapy in head and neck cancer. Adv Cell Mol Otolaryngol. 2015;3:27690.
Pinzon-Charry A, Maxwell T, Lopez JA. Dendritic cell dysfunction in cancer: a mechanism for immunosuppression. Immunol Cell Biol. 2005;83(5):451–61.
Rabinovich GA, Gabrilovich D, Sotomayor EM. Immunosuppressive strategies that are mediated by tumor cells. Annu Rev Immunol. 2007;25:267–96.
Sauter ER, Nesbit M, Watson JC, Klein-Szanto A, Litwin S, Herlyn M. Vascular endothelial growth factor is a marker of tumor invasion and metastasis in squamous cell carcinomas of the head and neck. Clin Cancer Res. 1999;5(4):775–82.
Moutsopoulos NM, Wen J, Wahl SM. TGF-beta and tumors--an ill-fated alliance. Curr Opin Immunol. 2008;20(2):234–40.
Gabrilovich DI, Chen HL, Girgis KR, Cunningham HT, Meny GM, Nadaf S, Kavanaugh D, Carbone DP. Production of vascular endothelial growth factor by human tumors inhibits the functional maturation of dendritic cells. Nat Med. 1996;2(10):1096–103.
Munn DH, Zhou M, Attwood JT, Bondarev I, Conway SJ, Marshall B, Brown C, Mellor AL. Prevention of allogeneic fetal rejection by tryptophan catabolism. Science. 1998;281(5380):1191–3.
Lob S, Konigsrainer A, Rammensee HG, Opelz G, Terness P. Inhibitors of indoleamine-2,3-dioxygenase for cancer therapy: can we see the wood for the trees? Nat Rev Cancer. 2009;9(6):445–52.
Spranger S, Spaapen RM, Zha Y, Williams J, Meng Y, Ha TT, Gajewski TF. Up-regulation of PD-L1, IDO, and T(regs) in the melanoma tumor microenvironment is driven by CD8(+) T cells. Sci Transl Med. 2013;5(200):200ra116.
Iversen TZ, Andersen MH, Svane IM. The targeting of indoleamine 2,3 dioxygenase -mediated immune escape in cancer. Basic Clin Pharmacol Toxicol. 2015;116(1):19–24.
Pilotte L, Larrieu P, Stroobant V, Colau D, Dolusic E, Frederick R, De Plaen E, Uyttenhove C, Wouters J, Masereel B, et al. Reversal of tumoral immune resistance by inhibition of tryptophan 2,3-dioxygenase. Proc Natl Acad Sci U S A. 2012;109(7):2497–502.
Holt D, Ma X, Kundu N, Fulton A. Prostaglandin E(2) (PGE (2)) suppresses natural killer cell function primarily through the PGE(2) receptor EP4. Cancer Immunol Immunother. 2011;60(11):1577–86.
Li H, Wu K, Tao K, Chen L, Zheng Q, Lu X, Liu J, Shi L, Liu C, Wang G, et al. Tim-3/galectin-9 signaling pathway mediates T-cell dysfunction and predicts poor prognosis in patients with hepatitis B virus-associated hepatocellular carcinoma. Hepatology. 2012;56(4):1342–51.
Avogadri F, Yuan J, Yang A, Schaer D, Wolchok JD. Modulation of CTLA-4 and GITR for cancer immunotherapy. Curr Top Microbiol Immunol. 2011;344:211–44.
Mocellin S, Nitti D. CTLA-4 blockade and the renaissance of cancer immunotherapy. Biochim Biophys Acta. 2013;1836(2):187–96.
Grosso JF, Jure-Kunkel MN. CTLA-4 blockade in tumor models: an overview of preclinical and translational research. Cancer Immun. 2013;13:5.
Peguet-Navarro J, Sportouch M, Popa I, Berthier O, Schmitt D, Portoukalian J. Gangliosides from human melanoma tumors impair dendritic cell differentiation from monocytes and induce their apoptosis. J Immunol. 2003;170(7):3488–94.
Anderson AC, Joller N, Kuchroo VK. Lag-3, Tim-3, and TIGIT: co-inhibitory receptors with specialized functions in immune regulation. Immunity. 2016;44(5):989–1004.
Monney L, Sabatos CA, Gaglia JL, Ryu A, Waldner H, Chernova T, Manning S, Greenfield EA, Coyle AJ, Sobel RA, et al. Th1-specific cell surface protein Tim-3 regulates macrophage activation and severity of an autoimmune disease. Nature. 2002;415(6871):536–41.
Anderson AC. Tim-3: an emerging target in the cancer immunotherapy landscape. Cancer Immunol Res. 2014;2(5):393–8.
Pedoeem A, Azoulay-Alfaguter I, Strazza M, Silverman GJ, Mor A. Programmed death-1 pathway in cancer and autoimmunity. Clin Immunol. 2014;153(1):145–52.
Chen J, Jiang CC, Jin L, Zhang XD. Regulation of PD-L1: a novel role of pro-survival signalling in cancer. Ann Oncol. 2016;27(3):409–16.
Ritprajak P, Azuma M. Intrinsic and extrinsic control of expression of the immunoregulatory molecule PD-L1 in epithelial cells and squamous cell carcinoma. Oral Oncol. 2015;51(3):221–8.
Parsa AT, Waldron JS, Panner A, Crane CA, Parney IF, Barry JJ, Cachola KE, Murray JC, Tihan T, Jensen MC, et al. Loss of tumor suppressor PTEN function increases B7-H1 expression and immunoresistance in glioma. Nat Med. 2007;13(1):84–8.
Lastwika KJ, Wilson W 3rd, Li QK, Norris J, Xu H, Ghazarian SR, Kitagawa H, Kawabata S, Taube JM, Yao S, et al. Control of PD-L1 expression by oncogenic activation of the AKT-mTOR pathway in non-small cell lung cancer. Cancer Res. 2016;76(2):227–38.
Casey SC, Tong L, Li Y, Do R, Walz S, Fitzgerald KN, Gouw AM, Baylot V, Gutgemann I, Eilers M, et al. MYC regulates the antitumor immune response through CD47 and PD-L1. Science. 2016;352(6282):227–31.
The authors thank Neeraj Kumar Singh and Upasana Mitra of Cellworks Research India Ltd., Whitefield, Bangalore, India for their predictive computational modeling work to show the relationship among KRAS mutations and PD-L1 expression.
This research was supported by a grant from the National Institute of Dental and Craniofacial Research (NIDCR) of the National Institutes of Health (R01 DE014390). NIDCR did not have a role in the study design, analysis, or interpretation of the data, writing of the manuscript, or the decision to submit the manuscript for publication.
Availability of data and materials
All data generated or analyzed during this study are included in this published article, figures, table, and additional files. Patient data, clinical characteristics, and exome sequencing information for each of 34 patients were obtained directly from Additional file 5: Table S3 of the Rizvi et al. study . To maintain anonymity, a random string generator was used to create a new random, 6-character uppercase alpha numeric string for each patient. This blinded both the identities of the patients in this study and their link to the prior published dataset we modeled.
Ethics approval and consent to participate
All patients had stage IV non-small cell lung cancer (NSCLC) and were treated at Memorial Sloan Kettering Cancer Center (n = 29) or the University of California at Los Angeles (n = 5) on protocol NCT01295827. All patients had consented to the Memorial Sloan Kettering Cancer Center Institutional Review Board-approved protocols permitting tissue collection and sequencing.
All patient related research was Memorial Sloan Kettering Cancer Center Institutional Review Board-approved and treated under protocol NCT0129827. Written informed consent was obtained from all patients.
Consent for publication
KAB has had a Cooperative Research and Development Agreement with Cellworks Group Inc., San Jose, CA. TA and SV work for Cellworks Group Inc., San Jose, California. DP, NKS, and UM work for Cellworks Research India Ltd., Whitefield, Bangalore, India. TAC is a co-founder of Gritstone Oncology and holds equity in the company. All other authors declare no competing financial interests in the findings of this study or with Cellworks Group Inc., San Jose, CA or Cellworks Research India Ltd., Whitefield, Bangalore, India.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
The original article has been corrected.
Table S1. Molecules with immunosuppressive functions used in simulation models to predict PD-1 drug responder status. (DOCX 55 kb)
Table S2. Individual mutational profiles of patient drug responders (n = 11) and nonresponders (n = 18). (DOCX 19 kb)
Figure S1. A schematic pyramid showing the levels of information used to develop the validated, cancer network. This network was created from published reports on cell receptors, signaling pathways, pathway signaling intermediates, activation factors, transcription factors, and enzyme kinetics. Information on each pathway node, its functionality, and its links with other genes, proteins, and pathways was manually researched, analyzed, curated, and aggregated to construct the integrated network maze. Every process or reaction was modeled mathematically using Michaelis Menton kinetics, mass action kinetics, and variations of these representations using ordinary differential equations (ODEs). Modeled events included but were not limited to interactions at the cell surface (e.g., binding of ligands to receptors, etc.), metabolic and cell signaling (e.g., signal pathway events, cross talk interactions among pathways, feedback control, etc.), activation and regulation of genes (e.g., activation links of transcription factors, etc.), intracellular processes such as proteasomal degradation, endoplasmic reticulum (ER) stress, oxidative stress, DNA damage and repair pathways, and cell cycle pathways. Time-dependent changes in signaling pathway fluxes of every biological reaction modeled utilizing modified ODEs were solved with a proprietary solver. Models were validated with a series of internal control analysis checks on predictions. These checks included assessing the effects of pathway molecule over-expression or knockdown on pathway predictions; effects of drugs on pathway predictions; and activation, regulation, and cross-talk interactions among pathway intermediates on pathway predictions. (DOCX 53 kb)
Supplementary Materials and Methods. (DOCX 109 kb)
Table S3. An example of the predictive computational modeling process. Specific details on an annexure section of the PD-L1 pathway show the step-by-step reactions, mechanisms, and reaction equations that occur. Such reactions also occurred in all of the other pathways. (DOCX 102 kb)
Table S4. Creation of the dendritic cell infiltration index for the patient SA97V5-specific simulation model. Chemokines CCL11, CCL20, CCL2, CCL3, CCL4, CCL5, CCL7, CX3CL1, and CXCL14, capable of trafficking of dendritic cells into the tumor microenvironment, were used to create the index. Individual chemokine percent expression (with respect to non-tumorigenic baseline controls) was predicted and given weightage so as to normalize the total to 1. The index was then calculated to be the sum of each prediction % change * weightage. (DOCX 16 kb)
Table S5. Analysis of the Discovery and Validation datasets was performed using Weka 3. The first number in each column represented the number of patient treatment responses correctly classified by the model. The second number represented the number of incorrectly classified patient treatment responses. The GOAL row at the bottom of each column described the number of correctly and incorrectly classified patients in the simulation models. The Test Set columns described the output from applying the model trained on the Discovery set to the Validation set. The “Test and Train” columns described test set accuracy (test set column) plus the training error (results obtained by applying the model to the training set, i.e. training error). (DOCX 19 kb)
Figure S2. An example of the relationship between PD-L1 expression and predicted TGFB1 expression using Weka 3 algorithms for all patients in the dataset. Similar trends were seen when comparing the PD-L1 expression level to the other 13 predicted molecules. For this, the number of gene mutations identified for each patient ranged from 2 to 36 with a total of 264 unique genes between all patients. This categorical data was preprocessed and expanded into a gene vector of length 264 to represent each of the unique genes. For each gene in the vector, the data was represented in binary; a 1 was assigned if the patient had a mutation in this gene, a 0 otherwise. Two datasets, one including gene mutations (Molecules and Gene Mutations) and one without (Molecules), were both used to learn prediction models. The Discovery and Validation datasets were determined based on the split provided to allow for comparable results. The performance of a subset of these models on the testing and training sets for both Molecules and Molecules and Gene Mutations datasets are shown. The SMO support vector machine with a normalized polynomial kernel had the best performance when applied to the molecule dataset. This model correctly identified 24 out of 29 patients whereas the simulation models correctly identified 25 of 29. This was only a difference of one match between the two prediction methods. Still, several other methods, while not performing as well overall, were able to identify 9 patients in the test dataset accurately. This was near the computational simulation model prediction capability in which 10 patients were successfully identified in the test dataset. In general, adding the gene mutation data to the molecule data either maintained or decreased the performance of a model. (DOCX 4114 kb)
Table S6. Comparisons of clinical and predicated responses and match scores. We used a cross-validation approach to assess the match scores in Table 1 of the PD-1 predicted responses against the PD-1 clinical responses in the Rizvi et al. 2015 Discovery dataset vs. the Validation dataset. We then pooled and re-partitioned the dataset into two new Training and Test datasets. We then used a similar cross-validation approach to assess the match scores of the PD-1 predicted responses vs. the PD-1 clinical responses. (DOCX 17 kb)