Predicting patient response

Abstract

Methods of predicting response of a subject suffering from cancer to a therapy, comprising calculating a resistance score for factors expressed by the subject, summing the resistance score to produce a total resistance score, wherein a total resistance score beyond a predetermined threshold indicates a subject is predicted to be resistant to the therapy, are provided.

Description

PREDICTING PATIENT RESPONSECROSS REFERENCE TO RELATED APPLICATIONS

[0001] This application claims the benefit of priority of U.S. Provisional Patent Application No. 63 / 465,026 filed on May 9, 2023 and of International Patent Application No. PCT / IL2023 / 050841 filed on August 10, 2023, the contents of which are all incorporated herein by reference in their entirety.FIELD OF THE INVENTION

[0002] The present invention is in the field of patient- specific diagnostics.BACKGROUND OF THE INVENTION

[0003] Immunotherapy based on immune checkpoint inhibitors (ICIs) represents a significant breakthrough in clinical oncology. ICIs augment an anti-tumor immune response by targeting checkpoint proteins such as PD-1, PD-L1 and CTLA-4 expressed on tumor and immune cells. Although ICI therapies can achieve unprecedented long-term disease control across multiple tumor types, efficacy varies widely between patients, with the majority exhibiting primary or subsequent acquired resistance to therapy. In metastatic non-small cell lung cancer (NSCLC) for which ICI regimens are standard of care, response rates range between 10%-50% depending on tumor PD-L1 expression as well as type and line of treatment. Identifying patients who are likely to benefit from ICI therapy is still a major clinical challenge because available predictive biomarkers are not sufficiently accurate.

[0004] To date, tumor PD-L1 expression and tumor mutational burden (TMB) are the most prominent biomarkers for predicting ICI response. While immunohistochemical tests for assessing PD-L1 expression in tumor tissue are used as companion diagnostics for informing treatment decisions in NSCLC, TMB is still not used routinely. According to current guidelines for NSCLC patients lacking oncogenic driver mutations, patients with high tumor PD-L1 expression (defined as PD-L1 expression on at least 50% of tumor cells) are eligible for first-line ICI monotherapy whereas ICI in combination with chemotherapy is the preferred choice for patients with PD-L1 expression <50%. However, clinical evidencedemonstrates limitations of the PD-L1 biomarker in predicting ICI response. For example, in the KEYNOTE-024 trial, approximately half of the PD-Ll-high patient cohort did not respond to pembrolizumab monotherapy. Moreover, several clinical trials report clinical benefit from ICI therapies in some patients with low tumor PD-L1 expression. Notably, although both PD-L1 and TMB are related to the mechanism of action of ICIs, these biomarkers do not account for the complexity of tumor-immune system interactions and the heterogenous mechanisms underpinning response and resistance to ICI therapy. In addition, the PD-L1 test requires tumor tissues, which are sometimes not available.

[0005] To address this, a more comprehensive characterization of the tumor, the tumor microenvironment (TME), peripheral immune cells and other host factors is needed. Indeed, a growing number of emerging predictive biomarkers for ICI outcome are based on tumor genomic features and expression patterns, the abundance and phenotype of tumorinfiltrating lymphocytes in the TME, peripheral T cell dynamics, and properties of other immune cell types. Importantly, integrative models combining several biomarkers show promise for improving predictive performance, presumably by better capturing the multifaceted nature of therapeutic benefit. For example, combining the PD-L1 and TMB biomarkers improves prediction of response to ICI therapy in lung cancer patients. In addition, several studies have demonstrated improved prediction of ICI outcomes using integrated genomic, transcriptomic, and immune repertoire data. Although such models are promising, they are limited in that they are based on multiple assays and usually require tumor tissue specimens.

[0006] Plasma proteomics represents a promising strategy for predictive biomarker discovery. Circulating blood contains thousands of proteins derived from the developing tumor, TME, peripheral immune cells and other host cells. As such, the plasma proteome reflects tumor-intrinsic properties, immune cell dynamics, angiogenesis, extracellular matrix remodeling and metabolic changes, making it a rich source of potential biomarkers that can be sampled in a minimally invasive manner and measured with a single assay. A method of determining patient-specific response to immunotherapy that integrates PD-L1 levels and plasma proteomic data, is greatly needed.SUMMARY OF THE INVENTION

[0007] The present invention provides methods of predicting response of a subject suffering from a cancer to a therapy, comprising calculating a resistance score for factors expressed by the subject, combining the resistance score to produce a total resistance score, wherein a total resistance score beyond a predetermined threshold indicates a subject is predicted to be resistant to the monotherapy or combination therapy. Response scores that are 1-the resistance score are also used.

[0008] According to a first aspect, there is provided a method of predicting response of a subject suffering from cancer to an anticancer therapy, the method comprising: a. receiving factor expression levels for a plurality of factors i. in a population of subjects suffering from cancer and known to respond to the anticancer therapy (responders); ii. in a population of subjects suffering from cancer and known to not respond to the anticancer therapy (non-responders); and iii. in the subject; b. calculate for factors of the plurality of factors a resistance score, wherein the calculating comprises applying a machine learning algorithm trained on a training set comprising the received factor expression levels in responders and non-responders to individual received factor expression levels from the subject and wherein the machine learning algorithm outputs the resistance score; and c. combine the calculated resistance scores to produce a total resistance score or determining the number of factors with a resistance score above a predetermined threshold to produce a total number of resistance- associated factors in the subject and convert the total number to a total resistance score; wherein a subject with a total resistance score beyond a predetermined threshold is predicted to not respond to the anticancer therapy and a subject with a total resistance score within the predetermined threshold is predicted to respond to the anticancer therapy; wherein the plurality of factors is selected from the factors provided in Tables 4 and 6; thereby predicting response of a subject to a monotherapy.

[0009] According to some embodiments, the total resistance score is converted to a total response score and wherein a total response score above a predetermined threshold indicates the subject is responsive to the anticancer therapy and a total response score below a predetermined threshold indicates the subject is not responsive to the anticancer therapy.

[0010] According to some embodiments, the anticancer therapy is a monotherapy comprising chemotherapy, targeted therapy or immunotherapy.[Oi l] According to some embodiments, the immunotherapy is anti-PD-l / PD-Ll immunotherapy .

[0012] According to some embodiments, the cancer is a PD-L1 high cancer and the monotherapy is an anti-PD-l / PD-Ll immunotherapy.

[0013] According to some embodiments, the monotherapy is chemotherapy.

[0014] According to some embodiments, the training set further comprises received factor expression levels in subjects suffering from cancer and known to respond to a combination therapy comprising an anti-PD-l / PD-Ll immunotherapy and chemotherapy (comboresponders) and received factor expression levels in subject suffering from cancer and known to not respond to the combination therapy (combo-non-responders).

[0015] According to some embodiments, the cancer is a PD-L1 low or negative cancer and the anticancer therapy is a combination therapy comprising an anti-PD-l / PD-Ll immunotherapy and chemotherapy.

[0016] According to some embodiments, the plurality of factors comprises at least two factors selected from the factors provided in Table 4.

[0017] According to some embodiments, the plurality of factors consists of factors selected from Table 4.

[0018] According to some embodiments, the plurality of factors comprises at least two factors selected from the factors provided in Table 5.

[0019] According to some embodiments, the plurality of factors consists of factors selected from Table 5.

[0020] According to some embodiments, the plurality of factors comprises at least two factors selected from the factors provided in Table 8.

[0021] According to some embodiments, the plurality of factors consists of factors selected from Table 8.

[0022] According to some embodiments, the plurality of factors comprises at least two factors selected from the factors provided in Table 9.

[0023] According to some embodiments, the plurality of factors consists of factors selected from Table 9.

[0024] According to some embodiments, the plurality of factors comprises at least two factors selected from the factors provided in Table 10.

[0025] According to some embodiments, the plurality of factors consists of factors selected from Table 10.

[0026] According to some embodiments, the plurality of factors comprises at least two factors selected from the factors provided in Table 11.

[0027] According to some embodiments, the plurality of factors consists of factors selected from Table 11.

[0028] According to some embodiments, the plurality of factors comprises at least two factors selected from the factors provided in Table 12.

[0029] According to some embodiments, the plurality of factors consists of factors selected from Table 12.

[0030] According to some embodiments, the plurality of factors comprises at least two factors selected from the factors provided in Table 6.

[0031] According to some embodiments, the plurality of factors consists of factors selected from Table 6.

[0032] According to some embodiments, the plurality of factors comprises at least two factors selected from the factors provided in Table 7.

[0033] According to some embodiments, the plurality of factors consists of factors selected from Table 7.

[0034] According to some embodiments, the responders and non-responders are determined based on progression free survival (PFS) at 1 year after initiation of the monotherapy or combination therapy.

[0035] According to some embodiments, the method comprises before (b) selecting a subset of the plurality of factors, wherein the subset comprises factors that best differentiate between the responders and non-responders, and wherein the calculating is for each factor of the subset.

[0036] According to some embodiments, the selecting comprises applying a statistical test to the received factor expression levels, optionally wherein the statistical test is a Kolmogorov-Smirnov test.

[0037] According to some embodiments, the subset consists of at least 50 factors.

[0038] According to some embodiments, the factor expression level is from a time point before administration of the anticancer therapy to the subject.

[0039] According to some embodiments, the combining is averaging.

[0040] According to some embodiments, the combining comprises determining the total number of factors with a resistance score above a predetermined threshold and producing a total resistance score proportional to the total number.

[0041] According to some embodiments, the converting comprises transformation by linear regression.

[0042] According to some embodiments, the cancer is selected from hepato-biliary cancer, cervical cancer, urogenital cancer, anogenital cancer, prostate cancer, thyroid cancer, ovarian cancer, nervous system cancer, ocular cancer, lung cancer, soft tissue cancer, bone cancer, pancreatic cancer, bladder cancer, skin cancer, intestinal cancer, hepatic cancer, rectal cancer, colorectal cancer, esophageal cancer, gastric cancer, gastroesophageal cancer, breast cancer, renal cancer, skin cancer, head and neck cancer, leukemia and lymphoma.

[0043] According to some embodiments, the cancer is selected from lung cancer, skin cancer, anogenital cancer, cervical cancer, renal cancer and head and neck cancer.

[0044] According to some embodiments, the cancer is non-small cell lung cancer (NSCLC).

[0045] According to some embodiments, the cancer is a tyrosine kinase inhibitor resistant cancer.

[0046] According to some embodiments, the predetermined threshold is determined by performing a cross-validation within the training set or is the median score of the training set.

[0047] According to some embodiments, the plurality of factors is at least 200 factors.

[0048] According to some embodiments, the factors expression levels are factors expression levels in a biological sample provided by the subjects.

[0049] According to some embodiments, the biological sample is selected from blood plasma, whole blood, blood serum or peripheral blood mononuclear cells.

[0050] According to some embodiments, the biological sample is blood plasma or blood serum.

[0051] According to some embodiments, the method further comprises administering the anticancer therapy to the subject predicted to respond to the anticancer therapy or administering an alternative therapy to the subject predicted to not respond to the anticancer therapy.

[0052] According to some embodiments, the method further comprises administering the monotherapy to the subject predicted to respond to the monotherapy or administering a combined therapy comprising the anti-PD-l / PD-Ll immunotherapy and chemotherapy to the subject predicted to not respond to the monotherapy.

[0053] According to some embodiments, the method further comprises administering the combination therapy to the subject predicted to respond to the combination therapy or administering an alternative therapy to the subject predicted to not respond to the combination therapy.

[0054] According to some embodiments, the anti-PD-l / PD-Ll immunotherapy is selected from Pembrolizumab, Nivolumab, Durvalumab and Atezolizumab.

[0055] According to some embodiments, the chemotherapy is selected from Carboplatin, Paclitaxel, Nab-Paclitaxel, Pemetrexed, Vinorelbine, and Cisplatin.

[0056] According to some embodiments, the combination therapy is selected from: a. Carboplatin, Durvalumab, and Paclitaxel; b. Atezolizumab, Bevacizumab, Carboplatin, and Paclitaxel; c. Carboplatin, Nab-Paclitaxel, and Pembrolizumab; d. Carboplatin, Nivolumab, and Paclitaxel;e. Carboplatin, Nivolumab, Pemetrexed; f. Carboplatin, Paclitaxel, Pembrolizumab; g. Carboplatin, Paclitaxel, Pembrolizumab, and radiation; h. Carboplatin, and Pembrolizumab; i. Carboplatin, Pembrolizumab, and Pemetrexed; j. Carboplatin, Pembrolizumab, and Vinorelbine; and k. Cisplatin, Pembrolizumab, and Pemetrexed.

[0057] According to some embodiments, predicting response comprises predicting overall survival.

[0058] According to some embodiments, predicting response comprises predicting progression free survival.

[0059] According to some embodiments, progression free survival is at 1 year after initiation of the monotherapy or combination therapy.

[0060] According to some embodiments, the subject suffers from a negative PD-L1 cancer.

[0061] According to some embodiments, PD-L1 high cancer comprises at least 50% of cancer cells being positive for surface expression of PD-L1 and PD-L1 low or negative cancer comprises fewer than 50% of cancer cells being positive for surface expression of PD-L1.

[0062] According to some embodiments, the PD-L1 low or negative cancer is PD-L1 negative cancer comprising less than 1% of cells being positive for surface expression of PD-L1.

[0063] According to some embodiments, the trained machine learning algorithm is trained by a method comprising: at a training stage, training a machine learning algorithm on a training set comprising:(i) factor expression levels of resistance-associated factors in samples from subjects suffering from cancer and known to be responsive to the anticancer therapy and factor expression levels of resistance-associated factors insamples from subjects suffering from the cancer and known to be non- responsive to the anticancer therapy; and(ii) labels associated with the responsiveness of the subjects suffering from the cancer; to produce a trained machine learning algorithm, wherein the trained machine learning algorithm is trained to output the resistance score and wherein the resistance-associated factors are selected from those provided in Tables 4 and 6.

[0064] According to some embodiments, the expression levels of resistance-associated factors are labeled with the labels.

[0065] According to some embodiments, the total resistance score predetermined threshold is 5 and a resistance score above 5 indicates the subject is resistant to the therapy or the total resistance score is converted to a total response score by the equation (10-total resistance score) and wherein a total response score above a predetermined threshold indicates the subject is responsive to therapy, optionally wherein the total response score predetermined threshold is 5.

[0066] Further embodiments and the full scope of applicability of the present invention will become apparent from the detailed description given hereinafter. However, it should be understood that the detailed description and specific examples, while indicating preferred embodiments of the invention, are given by way of illustration only, since various changes and modifications within the spirit and scope of the invention will become apparent to those skilled in the art from this detailed description.BRIEF DESCRIPTION OF THE DRAWINGS

[0067] Figures 1A-1C: Illustration of protein expression distributions in responders and non-responders populations at the single protein level. (1A-1C) Computer-generated examples of distributions of protein expression for responder and non-responder populations. Example of protein expression levels that may be considered as RAPs (lightgray dashed line) or are not RAPs (dark-gray dashed line) based on the population expression distribution data are shown.

[0068] Figure 2: Illustration of the RAP score of Equation 2 as implemented in Algorithm 1. The RAP score was calculated using synthetic data, where the responder and non-responder populations were generated by sampling from a normal distribution. The synthetic populations expression levels are shown in histograms, the responder population in darkgrey and the non-responder population in light-grey. Given these distributions, the RAP score was calculated for each expression level. The resulting RAP score is plotted in a blue curve, the values are indicated in the secondary Y-axis on the right.

[0069] Figures 3A-3C: RAP score threshold determination based on AUC as a function of RAP score. AUC at each RAP score was calculated and the peak of the obtained curve was determined as the threshold (dotted line) for determining a certain protein as a RAP or not. (3A-3B) Graph describing determination of the RAP score threshold using the mathematical approach for protein measures at (3A) T1 and (3B) TO. (3C) Graph describing determination of the RAP score threshold using the machine learning approach.

[0070] Figures 4A-4B: (4A) Bar chart showing the number of RAPs for each patient in the test cohort (n=30). Responders - light-grey; non-responders - dark-grey. (4B) ROC curve for the RAP analysis.

[0071] Figure 5: Heat map of hallmarks of cancer significantly enriched in six patients. The enrichment analysis was based on Fisher exact test (FDR < 0.05). Next to each patient identifier the number of the patient’s RAPs is indicated in brackets.

[0072] Figure 6: Protein-protein network of the key RAPs in the current cohort. The network is based on STRING database. Each node (protein) is colored based on the hallmarks of cancer to which it is associated. A black circular frame indicates a targetable RAP. The size of each node correlates with the number of patients that had the examined RAP. The compartment / s of each node is indicated in the middle (based on Human Protein Atlas). I, intracellular. M, membranal. S, soluble. A protein can have more than one compartment.

[0073] Figure 7: Chart of the significantly enriched hallmarks of cancer among the 19 RAPs. The analysis was done using Fisher exact test. Enrichment factor above 1 indicates enrichment.

[0074] Figure 8: Heat map of the protein expression levels of the 19 RAPs in healthy tissues. The expression data is based on Human Protein Atlas (HPA) database.

[0075] Figure 9: Heat map of the percentage of medium-high level staining in patients with different cancer types, including NSCLC. The expression data is based on Human Protein Atlas (HPA) database.

[0076] Figures 10A-10F: Clinical description of the 184 patients included in the analysis. (10A) Heatmap representing patient clinical characteristics: response to treatment (ORR1, ORR2, 1-year DCB); percent of cells expressing PD-L1 in biopsy immuno staining, a prognostic marker of response to treatment; treatment type: ICI only or combined treatment of ICI and chemotherapy; line of treatment: first line indicates ICI treatment was given as the first systemic treatment for NSCLC, advanced line indicates a previous non-ICI treatment was given before the current ICI treatment was administered. Sex indicates patient sex at birth. Histology indicates the lung cancer histological type (ADC-adenocarcinoma, SCC-squamous cell carcinoma). (10B-10C) Violin plots of the correlation of the patient age with response in each time point: (10B) ORR1 and (10C) ORR2. ORR1 and ORR2 are overall response determined 3 months and 6 months following treatment initiation, respectively. (10D-10E) Graphical display of the response groups in (10D) ORR1 and (10E) ORR2. NR = non -responders. R = responders (partial responders or complete responders). SD = stable disease (in the model they are included in the responder group). (10F) Graphical display of the division of the population into the development and the validation sets.

[0077] Figures 11A-11B: Performance of the classification model. ROC AUC was calculated using the total resistance score together with actual overall response evaluation at 3-month ORR, -6-month ORR and 1 year duration of clinical benefit (DCB) for both TO and Tl. Results at TO for the (11A, upper panel) development set and for the (11A, lower panel and 11B, upper panel) validation set are shown. (11B, lower panel) A similar classification model was generated based on Tl.

[0078] Figures 12A-12B: (12A) Patients sorted by their response probability score (calculated as 1 -resistance score) based on protein levels at TO. Actual observed response at 3 months ORR is indicated by color for each patient. (12B) Dot plot of the agreement between the predicted response probability based on TO protein expression and observed response probability at either 3 months, 6 months or 1 year. Each point on the graph indicates a specific patient, and the different time points are indicated by different hues and marker types. The black diagonal line indicates the line y = x, the red diagonal line indicates the fitted regression line for all the points and the goodness of fit of the regression (R2) is indicated. The horizontal lines indicate the average observed response probability for the 3 timepoints (color coded) across the entire validation set.

[0079] Figures 13A-13B: Survival analysis based on prediction results for ORR at 3 -months based on TO protein measurements for (13A) PFS and (13B) OS.

[0080] Figures 14A-14B: (14A) A functional network of all potential RAPs from this analysis. Each node represents a RAP, and the edge between nodes indicates a functional relation. Nodes with a larger size, and protein name provided, indicate investigational new drugs (INDs) in combination with immunotherapy. The nodes are colored based on the protein function. (14B) Functional network for two patients a predicted non-responder (upper) and a predicted responder (lower). RAPs detected for each patient are outlined in black. The non-responder patient had 44 RAPs detected and the responder had 10 RAPs detected.

[0081] Figure 15: Functional differences between RAPs higher in each response group. Each polygon in the Voronoi plot represents a RAP, and the size correlates with the difference between responders and non-responders. While non-responder RAPs are involved in splicing, signaling and cytoskeleton-related processes, the responder RAPs mainly involved proteolysis and cell adhesion. Each color indicates a different overall function.

[0082] Figure 16: A table describing the clinical parameters of the 339 patients included in the analysis.

[0083] Figure 17: Line graphs of patient number at each time point are indicated per response group (NR, non-responder; R, responder) and in total. The patient cohort was divided into development and validation sets.

[0084] Figure 18: Association of clinical parameters with CB at 3, 6 and 12 months. The examined clinical parameters are age, sex, histology type, treatment type, PD-L1 status and ECOG performance status. NSCC, Non-squamous cell carcinoma; SCC, Squamous cell carcinoma; CB, Clinical benefit; NCB, No clinical benefit; ICI, immune checkpoint inhibitor.

[0085] Figures 19A-19B: Performance of clinical parameter-based predictive models. (19A) Receiver operating characteristics (ROC) plot of the PD-Ll-based predictive model. (19B) ROC plot of the predictive clinical model based on PD-L1, sex, ECOG and treatment line. The area under the curve (AUC) values are indicated for each time point. CI, confidence interval.

[0086] Figures 20A-20B: (20A) Development of the RAP prediction model. A cohort of advanced stage NSCLC patients receiving ICI-based therapy was assembled. Pre-treatment blood samples were obtained, and plasma proteomes were profiled. Clinical benefit (CB) was assessed at 3, 6 and 12 months after starting treatment, and patients were followed upfor 2 years. A predictive model for CB was developed for each time point as follows: Proteins displaying differential plasma levels in CB and NCB patient populations were selected for model training using a statistical test. Such proteins are collectively termed Resistance Associated Proteins (RAPs). A predictive model for CB was developed per RAP using a machine learning algorithm. CB predictions inferred from each RAP were summed up to yield a RAP score per patient. RAP scores (total number of active RAPs) were linearly scaled to values between 0 and 1, enabling the conversion of a given patient’s RAP score into a CB probability. (20B) Development and validation of the RAP model. The cohort was divided into development and validation sets (75% and 25%, respectively). The development set was randomly divided into train and test sets (75% and 25%, respectively). The train set was used for RAP selection followed by model training resulting in a predictive model per RAP. Clinical benefit (CB) predictions were then generated per RAP for each patient in the test set. CB predictions from all selected RAPs were summed up to yield a RAP score per patient in the test set. The process was repeated 80 times, each time with a random division of development set patients into train and test sets. RAP scores were averaged per patient in the development set and linearly scaled. Model output is CB probability (a value between 0 and 1). The model was then locked and tested on the independent validation set.

[0087] Figure 21: The effect of RAP number on model performance per time point. Different numbers of RAPs (ranging from 1-400) were selected. For each number, the model was run 10 times. Model performance was assessed by ROC analysis. AUC is indicated. Based on this analysis, 50 was set as the cut-off for the number of selected RAPs.

[0088] Figures 22A-22G: RAP identification during model development. (22A) Histograms showing the number of identified RAPs grouped according to the number of times they were selected over 80 iterations. The top, middle and bottom histograms are for the 3-, 6- and 12- month time points, respectively. (22B) The total number of RAPs identified per time point. Some proteins measured by the SomaScan assay are redundant due to different aptamers binding to the same protein. The numbers of overall and non -redundant RAPs are indicated by blank and dotted bars, respectively. The numbers of non-redundant RAPs identified at least 40 times in a total of 80 iterations are indicated by lined bars. (22C) A Venn diagram showing the number of RAPs identified per time point. (22D) Hierarchical clustering-based heatmap showing the number of iterations in which a given protein was classified as a RAP. (22E) RAP cellular localization and potential cellular origin. Data were obtained from Human Protein Atlas. The same protein may be assigned more than onecellular localization. (22F) Voronoi plots displaying the main biological functions of the RAPs per time point. Each polygon represents a RAP, and the size correlates with the number of times that the protein was selected as a RAP. Proteins from the same KEGG biological process are grouped together (using default settings of the Proteomaps tool). (22G) Enrichment analysis of RAPs per time point. The enrichment analysis was performed for RAPs that were selected in at least 10 iterations. Fisher exact test (FDR < 0.1) was used.

[0089] Figures 23A-23E: Performance of the RAP predictive model. (23A) Bar plot showing predicted clinical benefit (CB) probabilities sorted from lowest to highest. Observed CB and no CB (NCB) patients are shown as light and dark blue bars, respectively. (23B) Overall survival analysis of patients stratified to high and low CB probability groups. The median CB probability per time point was used as the stratification threshold. HR, hazard ratio. CI, confidence interval. (23C) Progression-free survival analysis of patients stratified to high and low CB probability groups. The median CB probability per time point was used as the stratification threshold. HR, hazard ratio. CI, confidence interval. (23D) Predicted CB probability as a function of observed CB rate. Each dot represents a patient. The observed CB rate for each predicted CB probability datapoint refers to the proportion of observed CB patients within a patient group assigned the CB probability ±0.15. X=Y is indicated by a black line. The goodness of fit is indicated. (23E) Receiver operating characteristics (ROC) plot for the RAP model per time point. The area under the curve (AUC) is indicated. The dashed line indicates AUC = 0.5. CI, confidence interval.

[0090] Figure 24: Enrichment analysis for CB probabilities and observed CB rates at each time point. The enrichment analysis was done using 2D-enrichment test. The X-axis indicates the enrichment score for predicted CB probability. The Y-axis indicates the enrichment score for observed rates (as defined by the proportion of observed CB patients within a patient group assigned the CB probability ±0.15). The enrichment score is a value between 1 and -1. Positive and negative enrichment scores indicate enrichment in high and low CB probabilities, respectively and in high and low observed CB rates, respectively. The solid line indicates the X=Y line.

[0091] Figures 25A-25D: (25A-25C) Comparison between CB probabilities at sequential time points. Each dot represents a patient in the cohort. CB probability at one time point is plotted against CB probability at a subsequent time point. The colors indicate patient CB labels per time point, and whether the clinical benefit label changed between time points. (25A) Comparison between 3 and 6 months. (25B) Comparison between 3 and 12 months.(25C) Comparison between 6 and 12 months. (25D) Sankey plot displaying the flow of CB labeling over time. CB, Clinical benefit; NCB, No clinical benefit; NA, not available.

[0092] Figures 26A-26B: The RAP model outperforms PD-L1 and clinical parameter-based models. Predictive performance was compared across five models: RAP model; PD-L1- based model (PD-L1); Clinical model (CM); Integrated RAP + PD-L1; Integrated RAP + CM. (26A) Receiver operating characteristics (ROC) plots of the five models at each time point. The area under the curve (AUC) is indicated. The dashed line indicates AUC=0.5. CI, confidence interval. (26B) Forest plot comparing the five models. Top, Cox regression analysis based on overall survival (OS) data. Bottom, Cox regression analysis based on progression-free survival (PFS) data.

[0093] Figure 27: RAP model performance in different patient subsets. NSCC, non- squamous cell carcinoma; SCC, squamous cell carcinoma.

[0094] Figure 28: Kaplan-Meier plots of PD-Ll-high, PD-Ll-low and PD-L1 -negative patients in the overall cohort. Left, overall survival (OS); right, progression-free survival (PFS). Dashed line indicates median survival.

[0095] Figure 29: The RAP model predicts differential survival outcomes in patients with PD-L1 >50%. PD-Ll-high patients were stratified to high (left) and low (right) CB probability groups using the cohort median CB probability as the stratification threshold. Overall survival (OS; lower panel) and progression-free survival (PFS; upper panel) were evaluated in patients treated with ICI monotherapy vs combination Id-chemotherapy. Dashed line indicates median survival.

[0096] Figure 30: The RAP model predicts differential survival outcomes in patients with PD-L1 <50%. PD-Ll-low and PD-L1 -negative patients (PD-L1 -low-negative) were stratified to high (left) and low (right) CB probability groups using the cohort median CB probability as the stratification threshold. Overall survival (OS; lower panel) and progression-free survival (PFS; upper panel) were evaluated in patients treated with ICI monotherapy vs combination Id-chemotherapy. Dashed line indicates median survival.

[0097] Figure 31A-31C: (31A) Patient clinical data. (31B) Development of the PROphet prediction model. A cohort of advanced stage NSCLC patients receiving Id-based therapy was assembled. Pre-treatment blood samples were obtained, and plasma proteomes were profiled using SomaScan technology. Clinical benefit (CB) was assessed at 12 months after starting treatment, and patients were followed up for 2 years. A predictive model for CB wasdeveloped as follows: Proteins displaying differential plasma levels in CB and NCB patient populations were selected for model training using a statistical test. Such proteins are collectively termed Resistance Associated Proteins (RAPs). A predictive model for CB was developed per RAP using a machine learning algorithm. CB predictions inferred from each RAP were summed up to yield a RAP score per patient. RAP scores were linearly scaled to values between 0 and 1, enabling the conversion of a given patient’s RAP score into a CB probability, which determines the PROphet result, negative or positive on the scale of 0 and 10. (31C) Development and validation of the RAP model. The cohort was divided into development and validation sets. The development set was randomly divided into train and test sets (75% and 25%, respectively). The train set was used for RAP selection followed by model training resulting in a predictive model per RAP. Clinical benefit (CB) predictions were then generated per RAP for each patient in the test set. CB predictions from all selected RAPs were summed up to yield a RAP score per patient in the test set. The process was repeated 80 times, each time with a random division of development set patients into train and test sets. RAP scores were averaged per patient in the development set and linearly scaled. Model output is CB probability (a value between 0 and 1) that is translated to PROphet score. The model was then locked and tested on the independent validation set.

[0098] Figures 32A-32D: The PROphet predicts overall survival for patients receiving ICI- based therapy and outperforms PD-L1 based prediction. (32A) Kaplan-Meier plot for PD- Ll>50% versus PD-Ll<50%. (32B) Predicted CB probability based on PD-L1 prediction as a function of observed CB rate. Each dot represents a patient. The observed CB rate for each predicted CB probability datapoint refers to the proportion of observed CB patients within a patient group assigned CB probability ±0.04. The goodness of fit (R2) is indicated. (32C) Kaplan-Meier plot for patient stratification based on PROphet model. (32D) Predicted CB probability based on PROphet model as a function of observed CB rate. The observed CB rate for each predicted CB probability datapoint refers to the proportion of observed CB patients within a patient group assigned CB probability ±0.05.

[0099] Figures 33A-33B: PROphet is not predictive for chemotherapy patients. (33A) Kaplan-Meier plot for chemotherapy patients classified as PROphet positive or negative, with no significant difference between these two subgroups. (33B) Predicted CB probability based on PROphet model as a function of observed CB rate. Each dot represents a patient. The observed CB rate for each predicted CB probability datapoint refers to the proportion ofobserved CB patients within a patient group assigned CB probability ±0.05. The goodness of fit (R2) is indicated.

[0100] Figure 34: Flowchart of the patients participating in the PROphet + PD-L1 analysis.

[0101] Figures 35A-35H: The PROphet model predicts differential overall survival outcome between different subgroups when combined with PD-L1 expression level. (35A- 35C) Kaplan-Meier plots for PROphet-positive prediction with PD-Ll>50% patients (35A), PD-L1 1-49% (35B) and PD-L1<1% (35C). In 35A, PD-Ll>50% patients received either ICI-chemotherapy combination therapy or ICI monotherapy. In 35B. and 35C., PD-L1 1- 49% and PD-L1<1% patients that received ICI-chemotherapy combination were compared to patients receiving chemotherapy alone. (35D-35F) Kaplan-Meier plots for PROphet®- negative prediction with PD-Ll>50% patients (35D), PD-L1 1-49% (35E) and PD-L1<1% (35F). In 35D, PD-Ll>50% patients received either ICI-chemotherapy combination therapy or ICI monotherapy. In 35E. and 35F., PD-L1 1-49% and PD-L1<1% patients that received ICI-chemotherapy combination were compared to patients receiving chemotherapy alone. HR, hazard ratio. CI, confidence interval. (35G-35H) Kaplan Meier plots for PROphet positive (35G) or PROphet negative (35H) patients with PD-L1 expression level of 1-49%. The ICI-chemotherapy combination is compared either to ICI-monotherapy or to chemotherapy alone. HR, hazard ratio. CI, confidence interval.

[0102] Figures 36A-36F: The PROphet model predicts differential progression free survival outcome between different subgroups when combined with PD-L1 expression level. (36A-36C) Kaplan-Meier plots for PROphet-positive prediction with PD-Ll>50% patients (36A), PD-L1 1-49% (36B) and PD-L1<1% (36C). In 36A, PD-Ll>50% patients received either ICI-chemotherapy combination therapy or ICI monotherapy. In 36B. and 36C., PD- L1 1-49% and PD-L1<1% patients that received ICI-chemotherapy combination were compared to patients receiving chemotherapy alone. (36D-36F) Kaplan-Meier plots for PROphet®-negative prediction with PD-Ll>50% patients (36D), PD-L1 1-49% (36E) and PD-L1<1% (36F). In 36D, PD-Ll>50% patients received either ICI-chemotherapy combination therapy or ICI monotherapy. In 36E. and 36F., PD-L1 1-49% and PD-L1<1% patients that received ICI-chemotherapy combination were compared to patients receiving chemotherapy alone. HR, hazard ratio. CI, confidence interval.

[0103] Figures 37A-37C: Forest plot for multivariate analysis of the PROphet model. (37A) PD-Ll>50% patients. (37B) PD-L1 1-49% patients. (37C) PD-L1< 1% patients.

[0104] Figures 38A-38D: Comparison between PROphet-positive and -negative results. (38A-38B) Comparison for PD-Ll>50% patients. (38A), Patients receiving ICI monotherapy. (38B) Patients receiving ICI-chemotherapy combination. (38C) Comparison for PD-L1 1-49% patients receiving ICI-chemotherapy combination. (38D) Comparison for PD-L1<1% patients receiving ICI-chemotherapy combination.

[0105] Figures 39A-39H Applicability of the response prediction using the PROphet model in patients with melanoma, SCLC, and HPV-related malignancies. (39A-C) PROphet model prediction in melanoma cohort. (39A) Model ROC AUC for 1-year durable clinical benefit. (39B) Predicted response probability based on the PROphet model versus observed response probability. Each point indicates a specific patient. (39C) Kaplan Meier plots for PROphet positive and PROphet negative patients. (39D) Kaplan- Meier curves showing survival of PROphet positive and PROphet negative patients with SCLC. (39E-G) Kaplan-Meier curves showing survival of PROphet positive and PROphet negative patients with HPV-related malignancies including (39E) anogenital SCC, (39F) cervical cancer, and (39G) head and neck cancer. (39H) Kaplan-Meier curve for all HPV-related malignancies. Hazard ratio with 95% confidence intervals and p-values are indicated.

[0106] Figures 40A-40B: Applicability of the response prediction using the PROphet model for NSCLC patients with targetable mutations. (40A) Correlation between the PROphet score and overall survival duration of NSCLC patients with targeted mutations treated with PD-1 inhibitors. R2=0.41, p = 0.0073. (40B) Kaplan-Meier curves showing survival of PROphet positive and PROphet negative patients, HR=0.36, p=0.07.

[0107] Figures 41A-41I: A new PROphet response predictor based on number of RAPs only. (41A) Model ROC AUC is AUC=0.70 for 1-year durable clinical benefit. (41B) The agreement between the predicted and observed response probability. Each point indicates a specific patient. The black diagonal line indicates the line y = x, and the red dashed line is the regression result. The goodness of fit of the regression, the slope, the intercept and the p-value of the regression are indicated in the plot. (41C-D) Kaplan Meier plots for PROphet positive (41C) and PROphet negative (41D) patients with high PD-L1 expression (PD- Ll>50%). (41E-F) Kaplan Meier plots for PROphet positive (41E) and PROphet negative (41F) patients with low PD-L1 expression (PD-L1 1-49%). (41G-H) Kaplan Meier plots for PROphet positive (41G) and PROphet negative (41H) patients with negative PD-L1 expression (PD-L1<1%). (411) Model ROC AUC is AUC=0.65 for 85 chemotherapy-treated patients at 1-year DCB.

[0108] Figures 42A-E: Applicability of new PROphet response predictor in patients with renal cell carcinoma. (42A) Cohort description. Overall, data from 298 patients were explored. Following several filtration steps, 201 patients remained in the analysis. Basic clinical characteristics of the 201 patients participating in the analysis are described in the charts. (42B) Treatment combinations included Tyrosine Kinase Inhibitors (TKI), Immune Checkpoint Inhibitors (ICI), or a combination of both TKI and ICI. The breakdown of patients treated with each of these different regimens is provided. (42C) Kaplan-Meier survival analysis and Cox proportional hazards models comparing the overall survival (OS) and progression-free survival (PFS) in the PROphet-positive and PROphet-negative groups. (42D) Kaplan-Meier survival analysis and Cox proportional hazards models comparing the OS and PFS between the PROphet-positive and PROphet-negative groups in the different types of treatment (ICI only, TKI only or a combination of ICI and TKI or bevacizumab). (42E) Multivariate Cox proportional hazard regression analysis.DETAILED DESCRIPTION OF THE INVENTION

[0109] The present invention, in some embodiments, provides methods of predicting response of a subject comprising a tumor with high, low or negative levels of PD-L1 to immunotherapy. Here, we developed a novel and inherently robust machine learning (ML)- based model that analyzes proteomic profiles in pre-treatment blood plasma to predict benefit from ICI therapy in cancer patients. By integrating predictions from a large collection of proteomic biomarkers, the model accurately predicts clinical benefit at three time points along the treatment course and stratifies patients according to survival outcomes, or PFS, outperforming PD-Ll-based prediction. Furthermore, the model shows potential for further optimizing treatment selection when used together with PD-L1 classification. Overall, the model provides clinically valuable information to support treatment decisions in cancer.

[0110] The invention is based, at least in part on the discovery of a novel tool for supporting treatment decision for cancer patients. The RAP (PROphet) model provides two main clinical utilities. First, it successfully predicts therapeutic benefit at 12 months, displaying superior predictive capabilities over PD-L1 based models. Second, when used in combination with PD-L1 testing, the model helps in determining whether a patient should receive ICI alone, an Id-chemotherapy combination or an alternative therapy. Specifically, subjects with high PD-L1 levels and a high total response score are predicted to respond toICI as a monotherapy and need not be exposed to the adverse side effects resultant from chemotherapy. Subjects with high PD-L1 but a low total response score are advised to proceed with combination ICI-chemotherapy. In patients with low PD-L1 but a high total response score, treatment with combination ICI-chemotherapy is predicted to be effective, but patients with low PD-L1 and low total response score would be advised to consider alternative therapies.

[0111] By a first aspect, there is provided a method of predicting response of a subject suffering from a cancer to a therapy, the method comprising a. receiving expression levels for a plurality of factors i. in a population of subjects known to respond to the therapy (responders); ii. in a population of subjects known to not respond to the therapy (non-responders); and iii. in the subject; b. calculate for factors of the plurality of factors a resistance score, wherein the calculating comprises applying a machine learning algorithm and wherein the machine learning algorithm outputs the resistance score; and c. combine the calculated resistance scores to produce a total resistance score; wherein a subject with a total resistance score beyond a predetermined threshold is predicted to not respond to the therapy and a subject with a total resistance score within the predetermined threshold is predicted to respond to the therapy; and wherein the plurality of factors is selected from the factors provided in Table 4 and Table 6; thereby predicting the response of a subject to a therapy.

[0112] By another aspect, there is provided a method of predicting response of a subject suffering from a cancer to a therapy, the method comprising a. receiving expression levels for a plurality of factorsi. in a population of subjects known to respond to the therapy (responders); ii. in a population of subjects known to not respond to the therapy (non-responders); and iii. in the subject; b. calculate for factors of the plurality of factors a resistance score, wherein the calculating comprises applying a machine learning algorithm and wherein the machine learning algorithm outputs the resistance score; and c. determine the number of factors with a resistance score above a predetermined threshold to produce a total number of resistance- associated factors in the subject and convert that total number to a total resistance score; wherein a subject with a total resistance score beyond a predetermined threshold is predicted to not respond to the therapy and a subject with a total resistance score within the predetermined threshold is predicted to respond to the therapy; and wherein the plurality of factors is selected from the factors provided in Tables 4 and 6; thereby predicting the response of a subject to a therapy.

[0113] By another aspect, there is provided a method of predicting response of a subject suffering from a PD-L1 high cancer to a monotherapy comprising immunotherapy, the method comprising a. receiving expression levels for a plurality of factors i. in a population of subjects known to respond to the therapy (responders); ii. in a population of subjects known to not respond to the therapy (non-responders); and iii. in the subject;b. calculate for factors of the plurality of factors a resistance score, wherein the calculating comprises applying a machine learning algorithm and wherein the machine learning algorithm outputs the resistance score; and c. combine the calculated resistance scores to produce a total resistance score; wherein a subject with a total resistance score beyond a predetermined threshold is predicted to not respond to the monotherapy and a subject with a total resistance score within the predetermined threshold is predicted to respond to the monotherapy; and wherein the plurality of factors is selected from the factors provided in Tables 4 and 6; thereby predicting the response of a subject to a monotherapy.

[0114] By another aspect, there is provided a method of predicting response of a subject suffering from a PD-L1 low or negative cancer to a combination therapy comprising immunotherapy and chemotherapy, the method comprising a. receiving expression levels for a plurality of factors i. in a population of subjects known to respond to the therapy (responders); ii. in a population of subjects known to not respond to the therapy (non-responders); and iii. in the subject; b. calculate for factors of the plurality of factors a resistance score, wherein the calculating comprises applying a machine learning algorithm and wherein the machine learning algorithm outputs the resistance score; and c. combine the calculated resistance scores to produce a total resistance score; wherein a subject with a total resistance score beyond a predetermined threshold is predicted to not respond to the combination therapy and a subject with a total resistance score within the predetermined threshold is predicted to respond to the combination therapy;and wherein the plurality of factors is selected from the factors provided in Tables 4 and 6; thereby predicting the response of a subject to the combination therapy.

[0115] By another aspect, there is provided a method of predicting response of a subject to a therapy, the method comprising: a. receiving expression levels for a plurality of factors i. in a population of subjects known to respond to the therapy (responders); ii. in a population of subjects known to not respond to the therapy (non-responders); and iii. in the subject; b. calculate for at least one factor of the plurality of factors a resistance score; and c. classify a factor with a resistance score beyond a threshold as a resistance- associated factor; wherein a subject with a number of resistance-associated factors beyond a predetermined number is predicted to be resistant to the therapy, and wherein the plurality of factors is selected from the factors provided in Tables 4 and 6; thereby predicting the response of a subject to a therapy.

[0116] By another aspect, there is provided a method of predicting response of a subject to a therapy, the method comprising: a. receiving expression levels for a plurality of factors i. in a population of subjects known to respond to the therapy (responders); ii. in a population of subjects known to not respond to the therapy (non-responders); and iii. in the subject;b. calculate for at least one factor of the plurality of factors a resistance score; c. classify a factor with a resistance score beyond a threshold as a resistance- associated factor; d. sum the number resistance-associated factors; and e. apply a trained machine learning algorithm to the number of resistance- associated factors, wherein the trained machine learning algorithm outputs a total resistance score and a total resistance score beyond a predetermined threshold indicates the subject is resistant to the therapy; wherein the plurality of factors is selected from the factors provided in Tables 4 and 6; thereby predicting the response of a subject to a therapy.

[0117] By another aspect, there is provided a method of predicting response of a subject to a therapy, the method comprising: a. receiving expression levels for a plurality of factors i. in a population of subjects known to respond to the therapy (responders); ii. in a population of subjects known to not respond to the therapy (non-responders); and iii. in the subject; b. calculate for at least one factor of the plurality of factors a resistance score; c. classify a factor with a resistance score beyond a threshold as a resistance- associated factor; d. sum the number resistance-associated factors; and e. apply a trained machine learning algorithm to the number of resistance- associated factors and at least one clinical parameter, wherein the trained machine learning algorithm outputs a total resistance score and a total resistance score beyond a predetermined threshold indicates the subject is resistant to the therapy;wherein the plurality of factors is selected from the factors provided in Tables 4 and 6; thereby predicting the response of a subject to a therapy.

[0118] By another aspect, there is provided a method of predicting response of a subject to a therapy, the method comprising: a. receiving expression levels for a plurality of factors i. in a population of subjects known to respond to the therapy (responders); ii. in a population of subjects known to not respond to the therapy (non-responders); and iii. in the subject; b. calculate for factors of the plurality of factors a resistance score, wherein the resistance score is based on the similarity of the factor expression level in the subject to the factor expression level in the responders and the similarity to the factor expression level in the subject to the factor expression level in the non-responders and wherein the calculating comprises applying a trained machine learning algorithm that outputs the resistance score; and c. sum the calculated resistance scores to produce a total resistance score; wherein a subject with a total resistance score beyond a predetermined threshold is predicted to be resistant to said therapy; wherein the plurality of factors is selected from the factors provided in Tables 4 and 6; thereby predicting the response of a subject to a therapy.

[0119] By another aspect, there is provided a method comprising: at a training stage, training a machine learning algorithm on a training set comprising:(i) a number of resistance-associated factors expressed in samples from subjects suffering from a disease and known to be responsive to a therapy and from subjects suffering from the disease and known to be non-responsive to the therapy; and(ii) labels associated with the responsiveness of the subjects; to produce a trained machine learning algorithm.

[0120] By another aspect, there is provided a method comprising: at a training stage, training a machine learning algorithm on a training set comprising:(i) a number of resistance-associated factors expressed in samples from subjects suffering from a disease and known to be responsive to a therapy and from subjects suffering from the disease and known to be non-responsive to the therapy;(ii) at least one clinical parameter of the subjects; and(iii) labels associated with the responsiveness of the subjects; to produce a trained machine learning algorithm.

[0121] By another aspect, there is provided a method comprising: at a training stage, training a machine learning algorithm on a training set comprising:(i) factor expression levels of resistance-associated factors in samples from subjects suffering from a disease and known to be responsive to a therapy and from subjects suffering from the disease and known to be non-responsive to the therapy; and(ii) labels associated with the responsiveness of the subjects; to produce a trained machine learning algorithm.

[0122] By another aspect, there is provided a method comprising: at a training stage, training a machine learning algorithm on a training set comprising:(i) factor expression levels of resistance-associated factors in samples from subjects suffering from a disease and known to be responsive to a therapy and from subjects suffering from the disease and known to be non-responsive to the therapy;(ii) at least one clinical parameter of the subjects; and(iii) labels associated with the responsiveness of the subjects; to produce a trained machine learning algorithm.

[0123] In some embodiments, the method is a diagnostic method. In some embodiments, the method is an in vitro method. In some embodiments, the method is an ex vivo method. In some embodiments, the method is a computer implemented method. In some embodiments, the method is a statistical method. In some embodiments, the method is a method that cannot be performed in a human mind. In some embodiments, the method is a computerized method. In some embodiments, the processor is a computer processor. In some embodiments, the processor is a computer.

[0124] In some embodiments, the method is for predicting response to therapy. In some embodiments, the method is for determining response to therapy. In some embodiments, the method is for determining response score. In some embodiments, the method is for determining response probability. In some embodiments, a response probability is a response score. In some embodiments, the method is for determining clinical benefit probability. In some embodiments, the method is for determining overall survival. In some embodiments, the method is for determining progression free survival (PFS). In some embodiments, the method is for determining overall survival (OS). In some embodiments, the method is for determining survival probability. In some embodiments, determining is predicting. According to some embodiments, resistance score is determined. According to other embodiments, prediction of resistance probability is determined. According to some other embodiments, resistance probability below 20% indicates the subject is responsive to therapy. According to some other embodiments, resistance probability below 50% indicates the subject is responsive to therapy. According to some embodiments, response score is determined. According to other embodiments, prediction of response probability is determined. According to some other embodiments, response probability beyond 80% indicates the subject is responsive to therapy. According to some other embodiments, response probability beyond 50% indicates the subject is responsive to therapy. In some embodiments, beyond is above. In some embodiments, beyond is below. It will be understood by a skilled artisan that a scale can be designed to be measured in either direction and so above / below depends on the construction of the scale.

[0125] In some embodiments, the method is for determining if a subject is a responder to the therapy. In some embodiments, the method is for determining if a subject is a nonresponder to the therapy. In some embodiments, the method is for predicting a subject’s response to therapy. In some embodiments, the method is for monitoring response to the therapy. In some embodiments, the method is for determining if the therapy should continue,be adjusted (e.g., by further treating the subject with an additional therapy including but not limited to an agent determined by the RAP analysis provided hereinbelow) or changed. In some embodiments, the method is for determining a subject as being a responder to the therapy, or a non-responder to the therapy. In some embodiments, the method is for determining a subject as being a responder to the therapy, a non-responder to the therapy, or as having a stable diseased state. In some embodiments, the method is for predicting if a subject will respond to the therapy, or not respond to the therapy. In some embodiments, a responder is a responder to a monotherapy (mono-responder). In some embodiments, a responder is a responder to combination therapy (combo-responder). In some embodiments, a non-responder is a non-responder to monotherapy (mono-non-responder). In some embodiments, a non-responder is a non-responder to combination therapy (combo-non- responder). In some embodiments, the method is for determining if the subject will benefit or not benefit from the treatment. In some embodiments benefit is clinical benefit. In some embodiments, the method is a prognostic method.

[0126] In some embodiments, non-response comprises progressive disease. In some embodiments, non-response comprises cancer progression. In some embodiments, nonresponse comprises stable disease. In some embodiments, non-response comprises a worsening of symptoms of the disease. In some embodiments, non-response is not the development of side effects. In some embodiments, non-response comprises growth, metastasis and / or continued proliferation of a cancer. In some embodiments, non-response comprises no clinical benefit (NCB). In some embodiments, non-response is non-survival. In some embodiments, non-response is non-survival and / or cancer progression. In some embodiments, response is stable disease. In some embodiments, response comprises remission. In some embodiments, remission is minimal remission. In some embodiments, remission is partial remission. In some embodiments, remission is complete remission. In some embodiments, response is survival. In some embodiments, response is progression free survival. In some embodiments, response is long progression free survival. In some embodiments, response is measured using the overall response rate (ORR). A trained physician will be familiar with methods of determining response and specifically the ORR. In some embodiments, response is measured using Response Evaluation Criteria In Solid Tumors (RECIST). In some embodiments, response comprises survival. In some embodiments, survival is overall survival. In some embodiments, survival is progression free survival. In some embodiments, survival is overall survival. In some embodiments, responsecomprises a clinical benefit (CB). In some embodiments, response comprises a durable clinical benefit (DCB). In some embodiments, CB is DCB. In some embodiments, CB is PFS. In some embodiments, CB is PFS at 12 months after the commencement of treatment. In some embodiments, CB is PFS at 7 months after the commencement of treatment. In some embodiments, the population of subject known to respond and known not to respond are determined based on PFS and the predicted response comprises OS. In some embodiments, PFS is PFS at 12 months. In some embodiments, PFS is PFS at 7 months. In some embodiments, PFS is PFS at 6 months. In some embodiments, PFS is PFS at 3 months. In some embodiments, OS is OS at 12 months. In some embodiments, OS is OS at 7 months. In some embodiments, OS is OS at 6 months. In some embodiments, OS is OS at 3 months. In some embodiments, no clinical benefit or non-clinical benefit is the absent of a clinical benefit described herein.

[0127] In some embodiments, the subject is a mammal. In some embodiments, the subject is a human. In some embodiments, the subject suffers from a disease. In some embodiments, the disease is treatable by the therapy. In some embodiments, the disease is cancer. In some embodiments, the disease is treatable by an immune checkpoint inhibitor (ICI). In some embodiments, the cancer is a PD-L1 positive cancer. In some embodiments, the cancer is a PD-L1 high cancer. In some embodiments, the cancer is a PD-L1 low cancer. In some embodiments, the cancer is a PD-L1 negative cancer. In some embodiments, the cancer is a PD-L1 low or negative cancer. In some embodiments, the cancer is solid cancer. In some embodiments, the cancer is a tumor. In some embodiments, the cancer is selected from hepato-biliary cancer, cervical cancer, urogenital cancer (e.g., urothelial cancer), anogenital, testicular cancer, prostate cancer, thyroid cancer, ovarian cancer, nervous system cancer, ocular cancer, lung cancer, soft tissue cancer, bone cancer, pancreatic cancer, bladder cancer, skin cancer, intestinal cancer, hepatic cancer, rectal cancer, colorectal cancer, esophageal cancer, gastric cancer, gastroesophageal cancer, breast cancer (e.g., triple negative breast cancer), renal cancer (e.g., renal carcinoma), skin cancer, head and neck cancer, leukemia and lymphoma. In some embodiments, the cancer is selected from skin cancer, and lung cancer. In some embodiments, the cancer is skin cancer. In some embodiments, the cancer is lung cancer. In some embodiments, the skin cancer is melanoma. In some embodiments, the lung cancer is small cell lung cancer. In some embodiments, the lung cancer is non-small cell lung cancer. In some embodiments, the melanoma is non-resectable melanoma. In some embodiments, the melanoma is metastatic melanoma. In some embodiments, the cancer isan HPV (Human Papilloma Virus) positive cancer. In some embodiments, the cancer is an HPV-related cancer. In some embodiments, the cancer is anogenital cancer. In some embodiments, the anogenital cancer is anogenital squamous-cell carcinoma (SCC). In some embodiments, anogenital cancer comprises anal, cervical, penile, vaginal, and vulvar cancer. In some embodiments, the cancer is cervical cancer. In some embodiments, the cervical cancer is small-cell cervical cancer. In some embodiments, the cancer is a head and neck cancer. In some embodiments, the head and neck cancer is head and neck SCC (HNSCC). In some embodiments, the cancer is selected from lung cancer, skin cancer, anogenital cancer, cervical cancer and head and neck cancer. In some embodiments, the cancer is renal cancer. In some embodiments, the renal cancer is renal cell carcinoma.

[0128] In some embodiments, the cancer is resistant to a therapy. In some embodiments, the therapy is a non-immuno therapy. In some embodiments, the therapy is another therapy. In some embodiments, the therapy is targeted therapy. In some embodiments, the therapy is not anti-PD-l / Ll immunotherapy. In some embodiments, the cancer is resistant to a targeted therapy. In some embodiments, the targeted therapy is a tyrosine kinase inhibitor (TKI). In some embodiments, the targeted therapy is Bevacizumab. In some embodiments, the targeted therapy is TKI, Bevacizumab or both. In some embodiments, the subject has been previously treated with a TKI. In some embodiments, the subject was treated with and found resistant to a TKI. In some embodiments, the method is a method of determining if a subject resistant to a targeted therapy will respond to a PD-1 / L1 immunotherapy. In some embodiments, the subject comprises a TKI resistant cancer. In some embodiments, cancer is a TKI resistant NSCLC. In some embodiments, the cancer comprises a mutation of a tyrosine kinase receptor gene. In some embodiments, the tyrosine kinase receptor gene is selected from epidermal growth factor receptor (EGFR), Anaplastic lymphoma kinase (ALK) and Proto-oncogene tyrosine-protein kinase ROS (ROS1). In some embodiments, the therapy is chemotherapy. In some embodiments, the therapy is immunotherapy. In some embodiments, the therapy is TKI therapy. In some embodiments, the combination therapy is immunotherapy and chemotherapy. In some embodiments, the combination therapy is immunotherapy and TKI therapy. In some embodiments, the therapy is immunotherapy, TKI therapy and Bevacizumab.

[0129] In some embodiments, the subject is naive to therapy before the first determining. In some embodiments, the subject has not received the therapy before the first determining. In some embodiments, the subject has received the therapy previously. In some embodiments,the subject has previously been treated by a therapy other than the therapy. In some embodiments, the subject is simultaneously treated by a therapy other than the therapy. In some embodiments, the other therapy is a TGFB-trap fusion protein. In some embodiments, the other therapy is tyrosine kinase inhibitor. In some embodiments, the subject is naive to any therapy. In some embodiments, the subject is naive to immunotherapy. In some embodiments, the therapy is the first line of treatment. In some embodiments, the therapy is an advanced line of treatment.

[0130] In some embodiments, the therapy is an anticancer therapy. In some embodiments, the anticancer therapy is radiation. In some embodiments, the anticancer therapy is chemotherapy. In some embodiments, the therapy is immunotherapy. In some embodiments, the anticancer therapy is immunotherapy. In some embodiments, the anticancer therapy is targeted therapy. In some embodiments, the anticancer therapy is selected from radiation, chemotherapy, immunotherapy, targeted therapy, hormonal therapy, anti- angiogenic therapy and photodynamic therapy, thermo therapy, surgery, and a combination thereof. In some embodiments, the immunotherapy is selected from immune checkpoint inhibition, immune checkpoint modulation, immune checkpoint blockade, adoptive-cell transfer therapy, oncolytic virus therapy, vaccine therapy, immune system modulation and therapy using monoclonal antibodies. In some embodiments, an immunotherapy is selected from immune checkpoint inhibitors, immune checkpoint modulators, immune checkpoint blockers, adoptive-cell transfer therapy, oncolytic virus therapy, treatment vaccines, immune system modulators and monoclonal antibodies. In some embodiments, the immunotherapy is an immune checkpoint inhibitor. In some embodiments, the immunotherapy is immune checkpoint blockade. In some embodiments, the targeted therapy is tyrosine kinase inhibitors. In some embodiments, the targeted therapy is a TGFB-trap fusion protein.

[0131] In some embodiments, an immunotherapy is administered in combination with one or more conventional cancer therapy including chemotherapy, targeted therapy, steroids, and radiotherapy. Combinations of ICI and chemotherapy / radiotherapy / targeted therapy have been studied in multiple clinical trials. It will be understood by a skilled artisan that the predictive proteins disclosed herein are predictive in immunotherapy as a monotherapy, as well as part of a combination therapy. In some embodiments, the therapy is a monotherapy. In some embodiments, the monotherapy comprises an immunotherapy. In some embodiments, the monotherapy comprises a chemotherapy. In some embodiments, the monotherapy comprises a targeted therapy. In some embodiments, the monotherapy consistsof immunotherapy. In some embodiments, the monotherapy consists of chemotherapy. In some embodiments, the monotherapy consists of targeted therapy. In some embodiments, the monotherapy does not comprise chemotherapy. In some embodiments, the monotherapy is an anti-PD-l / PD-Ll immunotherapy. In some embodiments, the therapy is a combination therapy. In some embodiments, the combination therapy comprises an immunotherapy and another therapy. In some embodiments, the combination therapy comprises an immunotherapy and a chemotherapy. In some embodiments, the combination therapy comprises an immunotherapy and a targeted therapy. In some embodiments, the targeted therapy is a tyrosine kinase inhibitor. In some embodiments, the targeted therapy is an antitransforming growth factor beta (TGFB) agent. In some embodiments, the TGFB agent is a TGFB-trap fusion protein. TGFB-trap fusion proteins are well-known in the art and are disclosed for example in Knudson et al., “M7824, a novel bifunctional anti-PD-Ll / TGFp Trap fusion protein, promotes anti-tumor efficacy as monotherapy and in combination with vaccine”, Oncoimmunology. 2018 Feb 14;7(5):el426519 and Morris et al., “Bintrafusp alfa, an anti-PD-Ll:TGF-P trap fusion protein, in patients with ctDNA-positive, liver-limited metastatic colorectal cancer”, Cancer Res Commun. 2022 Sep;2(9):979-986, the contents of which are hereby incorporated by reference in their entirety. In some embodiments, the combination therapy further comprises radiation. In some embodiments, the combination therapy further comprises a non-anti-PD-l / PD-Ll immunotherapy. In some embodiments, the anti-PD-l / PD-Ll immunotherapy is selected from Pembrolizumab, Nivolumab, Durvalumab and Atezolizumab. In some embodiments, the anti-PD-l / PD-Ll immunotherapy is selected from Pembrolizumab, Nivolumab, Durvalumab, Atezolizumab, and Cemiplimab. In some embodiments, the immunotherapy comprises Pembrolizumab. In some embodiments, the immunotherapy comprises Nivolumab. In some embodiments, the immunotherapy comprises Durvalumab. In some embodiments, the immunotherapy comprises Atezolizumab. In some embodiments, the chemotherapy is selected from Carboplatin, Paclitaxel, Nab-Paclitaxel, Pemetrexed, Vinorelbine, and Cisplatin. In some embodiments, the chemotherapy is selected from Carboplatin, Paclitaxel, Nab-Paclitaxel, Pemetrexed, Vinorelbine, Cisplatin, dacarbazine, temozolomide, albumin-bound paclitaxel, and vinblastine. In some embodiments, the chemotherapy is Carboplatin. In some embodiments, the chemotherapy is Paclitaxel. In some embodiments, the chemotherapy is Nab-Paclitaxel. In some embodiments, the chemotherapy is Pemetrexed. In some embodiments, the chemotherapy is Vinorelbine. In some embodiments, the chemotherapy is Cisplatin. In some embodiments, the combination therapy comprises Carboplatin,Durvalumab, and Paclitaxel. In some embodiments, the combination therapy comprises Atezolizumab, Bevacizumab, Carboplatin, and Paclitaxel. In some embodiments, the combination therapy comprises Carboplatin, Nab-Paclitaxel, and Pembrolizumab. In some embodiments, the combination therapy comprises Carboplatin, Nivolumab, and Paclitaxel. In some embodiments, the combination therapy comprises Carboplatin, Paclitaxel, Pembrolizumab. In some embodiments, the combination therapy comprises Carboplatin, Nivolumab, Pemetrexed. In some embodiments, the combination therapy comprises Carboplatin, Paclitaxel, Pembrolizumab, and radiation. In some embodiments, the combination therapy comprises Carboplatin, and Pembrolizumab. In some embodiments, the combination therapy comprises Carboplatin, Pembrolizumab, and Pemetrexed. In some embodiments, the combination therapy comprises Carboplatin, Pembrolizumab, and Vinorelbine. In some embodiments, the combination therapy comprises Cisplatin, Pembrolizumab, and Pemetrexed. In some embodiments, the combination therapy comprises an anti-CTLA-4 antibody. In some embodiments, the CTLA-4 antibody is ipilimumab. In some embodiments, the CTLA-4 antibody is Tremelimumab. In some embodiments, the combination therapy comprises an anti-LAG3 antibody. In some embodiments, the LAG3 antibody is relatlimab. In some embodiments, the TKI is selected from Osimertinib, Erlotinib, Afatinib, Gefitinib, Dacomitinib, dacomitinib, Amivantamab-vmjw, Mobocertinib, Sotorasib, Adagrasib, Alectinib, Brigatinib, Lorlatinib, Ceritinib, Crizotinib, entrectinib, Dabrafenib, ceritinib, trametinib, Vemurafenib, Tepotinib, Capmatinib, Selpercatinib, Pralsetinib, Fam-trastuzumab, deruxtecan-nxki, Ado-trastuzumab, emtansine, Cabozantinib, Ado-trastuzumab emtansine, Larotrectinib, alectinib, Cetuximab, cobimetinib, Encorafenib, binimetinib, Lenvatinib, imatinib, dasatinib, nilotinib, sunitinib, pazopanib, axitinib, sorafenib, Belzutifan, tivozinib, Everolimus, temsirolimus, and ripretinib.

[0132] The NCCN guidelines for 2023 provide the following lists of treatment which may be used alone or in combination to treat NSCLC, Melanoma, RCC or SCLC as follows: NSCLC-ICI: Atezolizumab, pembrolizumab, Durvalumab, nivolumab, ipilimumab, Cemiplimab, Cemiplimab-rwlc, and Tremelimumab. TKIs: Osimertinib, Erlotinib, Afatinib, Gefitinib, Dacomitinib, dacomitinib, Amivantamab-vmjw, Mobocertinib, Sotorasib, Adagrasib, Alectinib, Brigatinib, Lorlatinib, Ceritinib, Crizotinib, entrectinib, Dabrafenib, ceritinib, trametinib, Vemurafenib, Tepotinib, Capmatinib, Selpercatinib, Pralsetinib, Famtrastuzumab, deruxtecan-nxki, Ado-trastuzumab, emtansine, Cabozantinib, Ado-trastuzumab emtansine, Larotrectinib, alectinib, and Cetuximab. Anti-VEGF: Ramucirumab, and bevacizumab. Chemotherapy: Carboplatin, paclitaxel, pemetrexed, gemcitabine, Cisplatin, docetaxel, vinorelbine, etoposide, and albumin-bound paclitaxel. Melanoma- ICI: Nivolumab, Pembrolizumab, Ipilimumab, and relatlimab. Targeted therapy: Dabrafenib, trametinib, Vemurafenib, cobimetinib, Encorafenib, binimetinib, and lenvatinib. KIT inhibitors: imatinib, dasatinib, nilotinib, and ripretinib. ROS1 fusions drugs: Crizotinib, and entrectinib. NTRK fusions drugs: Larotrectinib, and entrectinib. NRAS drugs: Binimetinib. Chemotherapy: dacarbazine, temozolomide, albumin-bound paclitaxel, carboplatin, paclitaxel, cisplatin, vinblastine, and dacarbazine. RCC-ICI: pembrolizumab, nivolumab, ipilimumab and avelumab. TKIs: Axitinib, Cabozantinib, Lenvatinib, Pazopanib, Sunitinib, Everolimus, Tivozanib, Erlotinib, and Belzutifan. SCLC- Chemotherapy: Cisplatin, etoposide, Carboplatin, irinotecan, Topotecan, Lurbinectedin, Cyclophosphamide, doxorubicin, vincristine, Docetaxel, Gemcitabine, Temozolomide, Vinorelbine, Bendamustine, platinum, and paclitaxel. ICI: atezolizumab, durvalumab, nivolumab, pembrolizumab, and ipilimumab.

[0133] In some embodiments, the immunotherapy is a plurality of immunotherapies. In some embodiments, the immunotherapy is immune checkpoint blockade. In some embodiments, the immunotherapy is immune checkpoint protein inhibition. In some embodiments, the immunotherapy is immune checkpoint protein modulation. In some embodiments, the immunotherapy comprises immune checkpoint inhibition. In some embodiments, the immunotherapy comprises immune checkpoint modulation. In some embodiments, immune checkpoint blockade and / or immune checkpoint inhibition comprises administering to the subject an immune checkpoint inhibitor. In some embodiments, inhibition comprises administering an immune checkpoint inhibitor. In some embodiments, the inhibitor is a blocking antibody. In some embodiments, the immunotherapy comprises immune checkpoint blockade. In some embodiments, modulation comprises administering an immune checkpoint modulator. In some embodiments, immune checkpoint modulation comprises administering to the subject an immune checkpoint modulator.

[0134] As used herein, the term “an immune checkpoint inhibitor (ICI)” refers to a single ICI, a combination of ICIs and a combination of an ICI with another cancer therapy. The ICI may be a monoclonal antibody, a dual-specific antibody, a humanized antibody, a fully human antibody, a fusion protein, or a combination thereof directed to blocking, inhibition or modulation of immune checkpoint proteins. In some embodiments, an immune checkpointinhibitor is an immune checkpoint modulator. In some embodiments, an immune checkpoint inhibitor is an immune checkpoint blocker. In some embodiments, the immune checkpoint protein is selected from PD-1 (Programmed Death-1); PD-L1; PD-L2; CTLA-4 (Cytotoxic T-Lymphocyte-Associated protein 4); A2AR (Adenosine A2A receptor), also known as AD0RA2A; B7-H3, also called CD276; B7-H4, also called VTCN1; B7-H5; BTLA (B and T Lymphocyte Attenuator), also called CD272; IDO (Indoleamine 2,3 -dioxygenase); KIR (Killer-cell Immunoglobulin-like Receptor); LAG-3 (Lymphocyte Activation Gene-3); TDO (Tryptophan 2,3 -dioxygenase); TIM-3 (T-cell Immunoglobulin domain and Mucin domain 3); VISTA (V-domain Ig suppressor of T cell activation); N0X2 (nicotinamide adenine dinucleotide phosphate NADPH oxidase isoform 2); SIGLEC7 (Sialic acid -binding immunoglobulin-type lectin 7), also called CD328; SIGLEC9 (Sialic acid-binding immunoglobulin-type lectin 9), also called CD329; 0X40 (Tumor necrosis factor receptor superfamily, member 4) also called CD134; and TIGIT. In some embodiments, the immune checkpoint protein is selected from PD-1, PD-L1 and PD-L2. In some embodiments, the immune checkpoint protein is selected from PD-1 and PD-L1. In some embodiments, the immune checkpoint protein is CTLA-4. In some embodiments, the immune checkpoint protein is PD-1. In some embodiments, immune checkpoint blockade comprises an anti-PD- 1 / PD-L1 / PD-L2 immunotherapy. In some embodiments, immune checkpoint blockade comprises an anti-PD-1 immunotherapy. In some embodiments, immune checkpoint blockade comprises an anti-PD-1 and / or anti-PD-Ll immunotherapy. In some embodiments, immune checkpoint blockade comprises an anti-CTLA-4 immunotherapy. In some embodiments, immune checkpoint blockade comprises an anti-PD-1 and / or anti-PD-Ll immunotherapy and an anti-CTLA-4 immunotherapy. In some embodiments, the immunotherapy is anti-PD-l / PD-Ll immunotherapy. In some embodiments, the immunotherapy is anti-PD-l / PD-Ll axis immunotherapy. In some embodiments, immune checkpoint blockade comprises an anti-LAG-3. In some embodiments, immune checkpoint blockade comprises an anti-PD-1 and / or anti-PD-Ll immunotherapy and an anti-LAG-3 immunotherapy .

[0135] In some embodiments, the resistance-associated factor is determined by a method comprising: a. receiving expression levels for a plurality of factors i. in a population of subjects known to respond to the therapy (responders);ii. in a population of subjects known to not respond to the therapy (non-responders); and iii. in the subject; b. calculate for at least one factor of the plurality of factors a resistance score; and c. classify a factor with a resistance score beyond a threshold as a resistance- associated factor.

[0136] In some embodiments, resistance-associated factors are in each subject. In some embodiments, resistance-associated factors are in the responders. In some embodiments, resistance-associated factors are in the non-responders. In some embodiments, the resistance-associated factors are labeled with the labels. In some embodiments, the expression levels of the resistance-associated factors are labeled with the labels. In some embodiments, the resistance-associated factors are resistance-associated proteins.

[0137] In some embodiments, the immunotherapy is a blocking antibody. In some embodiments, the immunotherapy is administration of a blocking antibody to the subject.

[0138] In some embodiments, the ICI is a monoclonal antibody (mAb) against PD-1 or PD- Ll. In some embodiments, the ICI is a mAb that neutralizes / blocks / inhibits / modulates the PD-1 pathway. In some embodiments, the ICI is a mAb against PD-1. In some embodiments, the anti-PD-1 mAb is Pembrolizumab (Keytruda; formerly called lambrolizumab). In some embodiments, the anti-PD-1 mAb is Nivolumab (Opdivo). In some embodiments, the anti- PD-1 mAb is Pidilizumab (CT0011). In some embodiments, the anti-PD-1 mAb is Cemiplimab (Libtayo, REGN2810). In some embodiments, the anti-PD-1 mAb is Toripalimab (Loqtorzi). In some embodiments, the anti-PD-1 mAb is Retifanlimab (Zynyz). In some embodiments, the anti-PD-1 mAb is Dostarlimab (Jemperli). In some embodiments, the anti-PD-1 mAb is any one of AMP-224, MED 10680, or PDR001. In some embodiments, the ICI is a mAb against PD-L1. In some embodiments, the anti-PD-Ll mAb is selected from Atezolizumab (Tecentriq), Avelumab (Bavencio), and Durvalumab (Imfinzi). In some embodiments, the anti-PD-Ll mAb is Atezolizumab. In some embodiments, the anti-PD-Ll mAb is Durvalumab. In some embodiments, the ICI is a mAb against CTLA-4. In some embodiments, the anti-CTLA-4 mAb is ipilimumab. In some embodiments, the anti-CTLA- 4 mAb is tremelimumab (Imjuno).In some embodiments, the ICI is a mAb against LAG-3. In some embodiments, the anti-LAG-3 mAb is Relatlimab.

[0139] As used herein, the term “factor” refers to any measurable biological molecule produced by the subject. In some embodiments, the factor is a protein. In some embodiments, the factor is an RNA. In some embodiments, the factor is a gene. In some embodiments, the factor is a secreted factor. In some embodiments, the secreted factors are selected from cytokines, chemokines, growth factors, soluble receptors and enzymes. In some embodiments, the factor is a soluble factor. In some embodiment, the factor is cellular factor. In some embodiments, the factor is membranal factor. In some embodiments, the factor is a cell adhesion molecule. In some embodiments, the factor is a factor found in blood. In some embodiments, the factor is a host-generated factor. In some embodiments, the factor is a resistance factor.

[0140] In some embodiments, the expression is protein expression. In some embodiments, the expression is secreted protein expression. In some embodiments, protein expression is soluble protein expression. In some embodiment, the expression is cellular protein expression. In some embodiments, the expression is membranal protein expression. In some embodiments, the expression is mRNA expression. In some embodiments, the expression is protein expression or mRNA expression. In some embodiments, expression level is concentration. In some embodiments, concentration is concentration level. It will be understood by a skilled artisan that when the presence of factor is measured in a liquid sample the expression can be provided as a concentration such as mg / ml or in arbitrary units according to the method of determining the factor’s expression. Arbitrary units can be selected from relative fluorescence unit (RFU) and Normalized Protein expression (NPX), or any other arbitrary units used as measurement of expression. The terms “expression” and “expression levels” are used herein interchangeably and refer to the amount of a gene product present in the sample. In some embodiments, gene product includes polynucleotide, e.g., tumor DNA, circulating tumor DNA, or circulating DNA. In some embodiments, the DNA is cell-free DNA. In some embodiments, determining comprises quantification of expression levels. In some embodiments, determining comprises normalization of expression levels. Determining of the expression level of the factor can be performed by any method known in the art. Methods of determining protein expression include, for example, antibody arrays, immunoblotting, immunohistochemistry, flow cytometry (FACS), ELISA, proximity extension assay (PEA), aptamer-based assays, proteomics arrays, proteome sequencing, flow cytometry (CyTOF), multiplex assays, mass spectrometry and chromatography. In some embodiments, determining protein expression levels comprises ELISA. In someembodiments, determining protein expression levels comprises protein array hybridization. In some embodiments, determining protein expression levels comprises mass-spectrometry quantification. In some embodiments, determining protein expression levels comprises PEA. In some embodiments, determining protein expression levels comprises aptamers. Methods of determining mRNA expression include, for example, RT-PCR, quantitative PCR, realtime PCR, microarrays, northern blotting, in situ hybridization, next generation sequencing, and massively parallel sequencing.

[0141] In some embodiments, the receiving factor expression levels is providing factor expression levels. In some embodiments, the receiving factor expression levels is determining factor expression levels. In some embodiments, determining is measuring. In some embodiments, the measuring is in a sample. In some embodiments, the expression levels were detected in a sample. In some embodiments, the sample is a biological sample. In some embodiments, the sample is provided by the subjects. In some embodiments, the sample is provided by the subject. In some embodiments, the sample is provided by a responder. In some embodiments, the sample is provided by a non-responder. In some embodiments, each subject of the population of responders provided a sample. In some embodiments, each subject of the population of non-responders provided a sample. In some embodiments, the sample is provided by a subject before receiving the therapy. In some embodiments, the factor expression level is from a time point before administration of the therapy. In some embodiments, the therapy is a monotherapy. In some embodiments, the therapy is an anti-PD-l / PD-Ll immunotherapy. In some embodiments, the therapy is a combination therapy. In some embodiments, the therapy is an anti-PD-l / PD-Ll immunotherapy and chemotherapy. In some embodiments, the sample is provided by a subject after receiving the therapy.In some embodiments, the determining is directly in the sample. In some embodiments, the determining is in the unprocessed sample. In some embodiments, the determining is in a processed sample. In some embodiments, the method further comprises processing the sample. In some embodiments, processing comprises isolating proteins from the sample. In some embodiments, processing comprises isolating nucleic acids from the sample. In some embodiments, the nucleic acid is RNA. In some embodiments, the RNA is mRNA. In some embodiments, the processing comprises lysing cells in the sample. In some embodiments, the nucleic acid is cell free DNA. In some embodiments, the nucleic acid is tumor cell DNA.

[0142] As used herein, the terms “peptide”, "polypeptide" and "protein" are used interchangeably to refer to a polymer of amino acid residues. In another embodiment, the terms "peptide", "polypeptide" and "protein" as used herein encompass native peptides, peptidomimetics (typically including non-peptide bonds or other synthetic modifications) and the peptide analogues peptoids and semipeptoids or any combination thereof. In another embodiment, the peptides polypeptides and proteins described have modifications rendering them more stable while in the body or more capable of penetrating into cells. In one embodiment, the terms “peptide”, "polypeptide" and "protein" apply to naturally occurring amino acid polymers. In another embodiment, the terms “peptide”, "polypeptide" and "protein" apply to amino acid polymers in which one or more amino acid residue is an artificial chemical analogue of a corresponding naturally occurring amino acid.

[0143] In some embodiments, the sample is a biological sample. In some embodiments, the sample is tissue. In some embodiments, the tissue sample is tumor sample. In some embodiments, the sample is a fluid. In some embodiments, the fluid is a biological fluid. In some embodiments, the sample is from the subject. In some embodiments, the sample is not a tumor sample. In some embodiments, the sample is a tumor sample. In some embodiments, the sample is not a hematopoietic cancer and the sample is a blood sample. In some embodiments, the sample is a sample that does not comprise cancer cells. In some embodiments, a blood sample comprises a peripheral blood sample, serum sample and a plasma sample. In some embodiments, the sample is a plasma sample. In some embodiments, the sample is a serum sample. In some embodiments, processing comprises isolating plasma. In some embodiments, processing comprises isolating serum. In some embodiments, the biological fluid is selected from, blood, plasma, serum, lymph, cerebral spinal fluid, urine, feces, semen, tumor fluid and gastric fluid. In some embodiments, the sample obtained from the subject and the responders are the same type of sample. In some embodiments, the sample obtained from the subject and the responders are different types of samples. In some embodiments, the sample obtained from the subject and the non-responders are the same type of sample. In some embodiments, the sample obtained from the subject and the non- responders are different types of samples. In some embodiments, the sample obtained from the non-responders and the responders are the same type of sample. In some embodiments, the sample obtained from the non-responders and the responders are different types of samples. In some embodiments, the sample obtained from the subject, the non-responders and the responders are the same type of sample. In some embodiments, the sample obtainedfrom the subject, the non-responders and the responders are blood samples. In some embodiments, the sample obtained from the subject, the non-responders and the responders are plasma samples. In some embodiments, the sample obtained from the subject, the non- responders and the responders are serum samples. In some embodiments, the sample obtained from the subject, the non-responders and the responders are different types of samples.

[0144] In some embodiments, a factor is a factor of the plurality of factors. In some embodiments, expression levels of a plurality of factors are received. In some embodiments, expression levels of at least 2, 3, 4, 5, 6 ,7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 30, 35, 40, 45, 50, 100, 150, 200, 250, 300, 400, 500, 600, 700, 800, 900, 1000, 1100, 1200, 1300, 1400, 1500, 1600, 1700, 1800, 1900, 2000, 3000, 4000, 5000, 6000, 7000, 8000, 9000, 10000, 12000, 15000, 20000, 25000, 30000, 35000, or 40000 factors is received. Each possibility represents a separate embodiment of the invention. In some embodiments, expression levels of at least 50 factors are received. In some embodiments, expression levels of at least 100 factors are received. In some embodiments, expression levels of at least 200 factors are received. In some embodiments, expression levels of at least 300 factors are received. In some embodiments, expression levels of at least 350 factors are received. In some embodiments, expression levels of at least 375 factors are received. In some embodiments, expression levels of at least 380 factors are received. In some embodiments, expression levels of at least 385 factors are received. In some embodiments, expression levels of at least 388 factors are received. In some embodiments, a plurality is at least 2, 3, 4, 5, 6 ,7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 30, 35, 40, 45, 50, 100, 150, 200, 250, 300, 400, 500, 600, 700, 800, 900, 1000, 1100, 1200, 1300, 1400, 1500, 1600, 1700, 1800, 1900, 2000, 3000, 4000, 5000, 6000, 7000, 8000, 9000, 10000, 12000, 15000, 20000, 25000, 30000, 35000, or 40000. Each possibility represents a separate embodiment of the invention. In some embodiments, a plurality is at least 50 factors. In some embodiments, a plurality is at least 100 factors. In some embodiments, a plurality is at least 200 factors. In some embodiments, a plurality is at least 221 factors. In some embodiments, a plurality is at least 300 factors. In some embodiments, a plurality is at least 350 factors. In some embodiments, a plurality is at least 375 factors. In some embodiments, a plurality is at least 380 factors. In some embodiments, a plurality is at least 385 factors. In some embodiments, a plurality is at least 388 factors. In some embodiments, expression levels of at least 50 factors are received. In some embodiments, expression levels of at least100 factors are received. In some embodiments, expression levels of at least 200 factors are received. In some embodiments, expression levels of at least 221 factors are received. In some embodiments, expression levels of at least 300 factors are received. In some embodiments, expression levels of at least 350 factors are received. In some embodiments, expression levels of at least 375 factors are received. In some embodiments, expression levels of at least 380 factors are received. In some embodiments, expression levels of at least 385 factors are received. In some embodiments, expression levels of at least 388 factors are received. In some embodiments, expression levels of at least 400 factors are received. In some embodiments, expression levels of at least 1000 factor are received. In some embodiments, expression levels of at least 5000 factors are received. In some embodiments, expression levels of at least 6000 factors are received. In some embodiments, expression levels of at least 7000 factors are received. In some embodiments, expression levels of at least 8000 factors are received.

[0145] In some embodiments, the factor is selected from a factor provided in Table 4. In some embodiments, the plurality of factors is selected from the factors provided in Table 4. In some embodiments, the plurality of factors comprises at least two factors selected from those provided in Table 4. In some embodiments, the plurality of factors consists of factors selected from Table 4. In some embodiments, the factor is selected from the group consisting of: KCNAB2, IL12B, IL23A, MCL1, KIR2DS2, AGA, RPN1, LAT, MFAP2, PUF60, MPZ, ACE, RNF122, TXNDC5, CDH15, FGFBP3, COL11A2, INPP5E, ADH7, MVK, RNF146, SOCS3, RBFOX2, ARFGAP1, SRSF6, RBM23, DDR1, APOF, TRA2B, MCTS1, TBCA, RGS7, PTPN9, CSNK1G2, ILF3, TPPP2, ARHGEF2, SRSF7, EWSR1, FSTL1, SPP1, FLRT2, FLRT3, VTN, ATP1B 1, WFIKKN2, NRAC, PKD2, HSPA9, EMC4, ASAP2, NAP1L2, HTR7, DCUN1D3, RBL2, MAD1L1, GRB14, RBBP5, NAB2, CSF1R, CCN4, GPD1, KLK3, CXCL13, GZMA, C9, IL12B, RAP1GAP, IGFBP1, DHX58, COPS2, IL1RAP, CCL25, HPX, ADM, CD93, ISG15, MYL6B, HSPA1A, MBD1, TRAPPC3, AKT2, CRLF1, FTL, RBBP4, BMPER, SERPINB5, PMP2, OTC, OTOR, AOC1, FGFBP1, ATRN, NAGLU, SAA1, SAA4, CLSTN1, GSS, DLD, EPHB4, PRSS27, MUC16, CFHR2, HTRA1, KRT19, RBP4, SMOC2, BTD, TXLNA, MZB 1, FADD, GSN, CDH17, LECT2, ADAMTSL1, RNASET2, SEMA4A, DDOST, BDH2, SNRPB2, G0LM1, RAB3A, CD46, SEPTIN6, WWOX, WDR5, HPCAL1, ALDH5A1, VAT1, SARS1, AFM, CDA, ITLN1, LRIG1, GREM1, PTGR2, UBE2L6, CLTA, GSR, PDCD6, SNCG, CRH, RGS21, UBE2R2, BASP1, GBP5, LMNB2, POP7, RAET1L, SEMA5B, CNTN3, UBL3, MMACHC, GTF2B,GCHFR, LRATD2, SGK1, TSEN15, SAR1B, CDK5RAP3, HAUS1, NKIRAS1, PH0SPH02, PCDH17, TRIM5, ALDH7A1, TXNL4A, CEP20, PDE1B, ITGA4, ITGB 1, LRFN3, ADGRB 1, SGSH, MGAT5, B3GAT1, MGAT5, FBLN7, APBB1IP, P0N2, PPP2R5D, RBF0X1, TIMP1, GEMIN7, CSNK1A1L, PHF11, BTN2A2, SKP2, SPATA46, LIN7A, BORCS5, ARRDC5, PCYT1A, PHYH, ANKRD63, VCX, NTAN1, STARD7, APOL2, FLT4, RCSD1, INIP, VMAC, XPNPEP3, IFNE, NEEFA, KDM8, NCBP1, USF2, LRRC75A, APCS, PLCD1, ESPN, RFX5, RPS6KB2, N0M02, TCEAL2, CES3, DYRK1A, CYP2C19, CFI, IGFBP3, IL6, LEP, CRTC3, VEGFA, IL1RAP, HGF, PLA2G2A, CCL25, SERPINA7, POR, CCN3, HPX, IGFBP1, MMP3, FGA, FGB, FGG, BCAM, SPINT1, HAT1, GHR, CFP, CNTN1, SERPINF2, IL 19, MB, C9, IGHM, LBP, NAAA, HAPLN1, IDS, NIDI, ACAN, TGFBI, DLL4, FCGR3B, ACY1, IBSP, SERPINA4, POSTN, SELE, B2M, HAMP, SERPINA1, AHSG, CKB, CKM, PROC, ANGPTL4, MBD4, PSMD7, IGHE, CXCL10, KLKB1, CFH, PFDN5, RBM39, DCTPP1, PRSS22, KYNU, IL6, SERPINA6, ITIH4, SFN, CCL7, LYZ, MMP13, STC1, CAPG, PI3, GPC5, HRG, SCGB2A1, SIRT2, TNFAIP6, CD300C, GPNMB, KRT18, TNFSF14, LEPR, PRKCG, FGL1, PGLYRP2, NPFF, MFAP4, TMX3, PRKCSH, DEFB I 12, SEMA4D, ACP6, AFP, NGF, FTH1, FTL, DMKN, EPHA10, CHRDL2, TP53, AOC1, IFNA8, CSH1, CSH2, TNC, PLTP, CCN1, CLSTN3, OIT3, GGT2, FMOD, C5orf38, VWA1, INHBC, ADGRF5, C1QL2, PCYOX1, AOC2, CFHR4, LRRC15, POSTN, UBE2J1, GFRAL, IGF2, LILRB5, LILRA6, APOA2, VWA2, DEPPI, C1QTNF3, SERPINA9, CFHR5, DLG3, GLTPD2, HBQ1, ENTPD1, AGGF1, NRG2, SPON2, FAM241B, JAML, BCHE, GPNMB, APOD, DLL1, PEAR1, RSPO4, LEP, ARL8B, PCDH10, MFAP3L, CD14, COL15A1, HAVCR1, ARHGEF10, MAN1A2, CRYZL1, TFPI2, PLXDC1, ACP2, BTD, MFAP2, ITIH2, EFCAB14, PLA1A, GZMK, YBX1, IDO1, NQO1, SPOCK3, and NXTL The amino acid sequences of these factors can be found in the Uniprot database, for example, and each factor’s Uniprot accession number is provided in Table 4. Further, methods, reagents, and assays for measuring expression levels of these factors are well known in the art and are commercially available.

[0146] In some embodiments, the plurality of factors comprises Voltage-gated potassium channel subunit beta-2 (KCNAB2). The human KCNAB2 gene can be found at Entrez gene #8514. The human KCNAB2 protein can be found at Uniprot ID Q133O3. In some embodiments, the plurality of factors comprises Interleukin 12 subunit beta (IL12B). The human IL12B gene can be found at Entrez gene #3593. The human IL12B protein can befound at Uniprot ID P29460. In some embodiments, the plurality of factors comprises Interleukin 23 subunit alpha (IL23A). The human IL23A gene can be found at Entrez gene #51561. The human IL23A protein can be found at Uniprot ID Q9NPF7. In some embodiments, the plurality of factors comprises Induced myeloid leukemia cell differentiation protein Mcl-1 (MCL1). The human MCLlgene can be found at Entrez gene #4170. The human MCLlprotein can be found at Uniprot ID Q07820. In some embodiments, the plurality of factors comprises Killer cell immunoglobulin-like receptor, two domains, short cytoplasmic tail, 1 (KIR2DS2). The human KIR2DS2gene can be found at Entrez gene #3806. The human KIR2DS2protein can be found at Uniprot ID Q14954. In some embodiments, the plurality of factors comprises N(4)-(beta-N-acetylglucosaminyl)-L- asparaginase (AGA). The human AGA gene can be found at Entrez gene #175. The human AGA protein can be found at Uniprot ID P20933. In some embodiments, the plurality of factors comprises Dolichyl-diphosphooligo saccharide — protein glycosyltransferase subunit 1 (RPN1). The human RPNlgene can be found at Entrez gene #6184. The human RPNlprotein can be found at Uniprot ID P04843. In some embodiments, the plurality of factors comprises Linker for activation of T cells (LAT). The human LAT gene can be found at Entrez gene #27040. The human LAT protein can be found at Uniprot ID 043561. In some embodiments, the plurality of factors comprises Microfibrillar-associated protein 2 (MFAP2). The human MFAP2 gene can be found at Entrez gene #4237. The human MFAP2protein can be found at Uniprot ID P55001. In some embodiments, the plurality of factors comprises Poly(U)-binding-splicing factor PUF60 (PUF60). The human PUF60 gene can be found at Entrez gene #22827. The human PUF60 protein can be found at Uniprot ID Q9UHX1. In some embodiments, the plurality of factors comprises Myelin protein zero (MPZ). The human MPZ gene can be found at Entrez gene #4359. The human MPZ protein can be found at Uniprot ID P25189. In some embodiments, the plurality of factors comprises Angiotensin-converting enzyme (ACE). The human ACE gene can be found at Entrez gene #1636. The human ACE protein can be found at Uniprot ID P12821. In some embodiments, the plurality of factors comprises RING finger protein 122 (RNF122). The human RNF122 gene can be found at Entrez gene #79845. The human RNF122 protein can be found at Uniprot ID Q9H9V4. In some embodiments, the plurality of factors comprises Thioredoxin domain-containing protein 5 (TXNDC5). The human TXNDC5gene can be found at Entrez gene #81567. The human TXNDC5protein can be found at Uniprot ID Q8NBS9. In some embodiments, the plurality of factors comprises Cadherin-15 (CDH15). The human CDH15 gene can be found at Entrez gene #1013. The human CDH15 protein can be found at UniprotID P55291. In some embodiments, the plurality of factors comprises Fibroblast growth factor-binding protein 1 (FGFBP3). The human FGFBP3 gene can be found at Entrez gene #9982. The human FGFBP3 protein can be found at Uniprot ID Q14512. In some embodiments, the plurality of factors comprises Collagen alpha-2(XI) chain (C0L11A2). The human C0L11A2 gene can be found at Entrez gene #1302. The human COL11A2 protein can be found at Uniprot ID P13942. In some embodiments, the plurality of factors comprises 72 kDa inositol polyphosphate 5-phosphatase (INPP5E). The human INPP5E gene can be found at Entrez gene #56623. The human INPP5E protein can be found at Uniprot ID Q2YD81. In some embodiments, the plurality of factors comprises Alcohol dehydrogenase class 4 mu / sigma chain (ADH7). The human ADH7gene can be found at Entrez gene #131. The human ADH7protein can be found at Uniprot ID P40394. In some embodiments, the plurality of factors comprises Mevalonate kinase (MVK). The human MVK gene can be found at Entrez gene #4598. The human MVK protein can be found at Uniprot ID Q03426. In some embodiments, the plurality of factors comprises RING finger protein 146 (RNF146). The human RNF146 gene can be found at Entrez gene #81847. The human RNF146 protein can be found at Uniprot ID Q9NTX7. In some embodiments, the plurality of factors comprises Suppressor of cytokine signaling 3 (SOCS3). The human SOCS3 gene can be found at Entrez gene #9021. The human SOCS3 protein can be found at Uniprot ID 014543 or Q6FI39. In some embodiments, the plurality of factors comprises RNA binding motif protein 9 (RBFOX2). The human RBFOX2 gene can be found at Entrez gene #23543. The human RBFOX2 protein can be found at Uniprot ID 043251. In some embodiments, the plurality of factors comprises ADP-ribosylation factor GTPase-activating protein 1 (ARFGAP1). The human ARFGAP1 gene can be found at Entrez gene #55738. The human ARFGAP1 protein can be found at Uniprot ID Q8N6T3. In some embodiments, the plurality of factors comprises Splicing factor, arginine / serine-rich 6 (SRSF6). The human SRSF6 gene can be found at Entrez gene #6431. The human SRSF6 protein can be found at Uniprot ID Q13247. In some embodiments, the plurality of factors comprises Probable RNA-binding protein 23 (RBM23). The human RBM23 gene can be found at Entrez gene #55147. The human RBM23 protein can be found at Uniprot ID Q88U06. In some embodiments, the plurality of factors comprises Discoidin domain receptor family, member 1 (DDR1). The human DDR1 gene can be found at Entrez gene #780. The human DDR1 protein can be found at Uniprot ID Q0345. In some embodiments, the plurality of factors comprises Apolipoprotein F (APOF). The human APOF gene can be found at Entrez gene #319. The human APOF protein can be found at Uniprot ID Q13790. In some embodiments,the plurality of factors comprises Transformer-2 protein homolog beta (TRA2B). The human TRA2B gene can be found at Entrez gene #6434. The human TRA2B protein can be found at Uniprot ID P62995. In some embodiments, the plurality of factors comprises MCTS1, reinitiation and release factor (MCTS 1). The human MCTS 1 gene can be found at Entrez gene #28985. The human MCTS1 protein can be found at Uniprot ID Q9ULC4. In some embodiments, the plurality of factors comprises Tubulin- specific chaperone A (TBCA). The human TBCA gene can be found at Entrez gene #6902. The human TBCA protein can be found at Uniprot ID 075347. In some embodiments, the plurality of factors comprises Regulator of G-protein signaling 7 (RGS7). The human RGS7 gene can be found at Entrez gene #6000. The human RGS7 protein can be found at Uniprot ID P49802. In some embodiments, the plurality of factors comprises Tyro sine-protein phosphatase non-receptor type (PTPN9). The human PTPN9 gene can be found at Entrez gene #5780. The human PTPN9 protein can be found at Uniprot ID P43378. In some embodiments, the plurality of factors comprises casein kinase 1 gamma 2 (CSNK1G2). The human CSNK1G2 gene can be found at Entrez gene #1455. The human CSNK1G2 protein can be found at Uniprot ID P78368. In some embodiments, the plurality of factors comprises Interleukin enhancerbinding factor 3 (ILF3). The human ILF3 gene can be found at Entrez gene #3609. The human ILF3 protein can be found at Uniprot ID Q 12906. In some embodiments, the plurality of factors comprises Tubulin polymerization-promoting protein (TPPP2). The human TPPP2 gene can be found at Entrez gene #11076. The human TPPP2 protein can be found at Uniprot ID 094811. In some embodiments, the plurality of factors comprises Rho guanine nucleotide exchange factor 2 (ARHGEF2). The human ARHGEF2 gene can be found at Entrez gene #9181. The human ARHGEF2 protein can be found at Uniprot ID Q92974. In some embodiments, the plurality of factors comprises Serine / arginine-rich splicing factor 7 (SRSF7). The human SRSF7 gene can be found at Entrez gene #6432. The human SRSF7 protein can be found at Uniprot ID Q16629. In some embodiments, the plurality of factors comprises RNA-binding protein EWS (EWSR1). The human EWSR1 gene can be found at Entrez gene #2130. The human EWSR1 protein can be found at Uniprot ID Q01844. In some embodiments, the plurality of factors comprises Follistatin-related protein 1 (FSTL1). The human FSTL1 gene can be found at Entrez gene #11167. The human FSTL1 protein can be found at Uniprot ID Q 12841. In some embodiments, the plurality of factors comprises Osteopontin also known as secreted phosphoprotein 1 (SPP1). The human SPP1 gene can be found at Entrez gene #6696. The human SPP1 protein can be found at Uniprot ID P10451. In some embodiments, the plurality of factors comprises Fibronectin leucine richtransmembrane protein 2 (FLRT2). The human FLRT2 gene can be found at Entrez gene #23768. The human FLRT2 protein can be found at Uniprot ID 043155. In some embodiments, the plurality of factors comprises Fibronectin leucine rich transmembrane protein 3 (FLRT3). The human FLRT3 gene can be found at Entrez gene #23767. The human FLRT3 protein can be found at Uniprot ID Q9NZU0. In some embodiments, the plurality of factors comprises Vitronectin (VTN). The human VTN gene can be found at Entrez gene #7448. The human VTN protein can be found at Uniprot ID P04004. In some embodiments, the plurality of factors comprises Sodium / potassium-transporting ATPase subunit beta-1 (ATP1B1). The human ATP1B 1 gene can be found at Entrez gene #481. The human ATP1B 1 protein can be found at Uniprot ID P05026. In some embodiments, the plurality of factors comprises WAP, follistatin / kazal, immunoglobulin, kunitz and netrin domain containing 2 (WFIKKN2). The human WFIKKN2 gene can be found at Entrez gene #134857. The human WFIKKN2 protein can be found at Uniprot ID Q8TEU8. In some embodiments, the plurality of factors comprises Nutritionally-regulated adipose and cardiac enriched protein homolog (NRAC). The human NRAC gene can be found at Entrez gene #400258. The human NRAC protein can be found at Uniprot ID Q8N912. In some embodiments, the plurality of factors comprises Polycystin-2 (PKD2). The human PKD2 gene can be found at Entrez gene #5311. The human PKD2 protein can be found at Uniprot ID Q13563. In some embodiments, the plurality of factors comprises Mitochondrial 70kDa heat shock protein (HSPA9). The human HSPA9 gene can be found at Entrez gene #3313. The human HSPA9 protein can be found at Uniprot ID P38646. In some embodiments, the plurality of factors comprises ER membrane protein complex subunit 4 (EMC4). The human EMC4 gene can be found at Entrez gene #51234. The human EMC4 protein can be found at Uniprot ID Q5J8M3. In some embodiments, the plurality of factors comprises Arf-GAP with SH3 domain, ANK repeat and PH domain-containing protein 2 (ASAP2). The human ASAP2 gene can be found at Entrez gene #8853. The human ASAP2 protein can be found at Uniprot ID 043150. In some embodiments, the plurality of factors comprises Nucleosome assembly protein 1-like 2 (NAP1L2). The human NAP1L2 gene can be found at Entrez gene #4674. The human NAP1L2 protein can be found at Uniprot ID Q9ULW6. In some embodiments, the plurality of factors comprises 5-HT7 receptor (HTR7). The human HTR7 gene can be found at Entrez gene #3363. The human HTR7 protein can be found at Uniprot ID P34969. In some embodiments, the plurality of factors comprises DCNl-like protein 3 (DCUN 1D3). The human DCUN 1D3 gene can be found at Entrez gene #123879. The human DCUN 1D3 protein can be found at Uniprot ID Q8IWE4. In some embodiments, the pluralityof factors comprises Retinoblastoma-like protein 2 (RBL2). The human RBL2 gene can be found at Entrez gene #5934. The human RBL2 protein can be found at Uniprot ID Q08999. In some embodiments, the plurality of factors comprises Mitotic spindle assembly checkpoint protein MAD1 (MAD IL 1). The human MAD IL 1 gene can be found at Entrez gene #8379. The human MAD1L1 protein can be found at Uniprot ID Q9Y6D9. In some embodiments, the plurality of factors comprises Growth factor receptor-bound protein 14 (GRB14). The human GRB14 gene can be found at Entrez gene #2888. The human GRB14 protein can be found at Uniprot ID Q14449. In some embodiments, the plurality of factors comprises Retinoblastoma-binding protein 5 (RBBP5). The human RBBP5 gene can be found at Entrez gene #5929. The human RBBP5 protein can be found at Uniprot ID Q15291. In some embodiments, the plurality of factors comprises NGFI-A-binding protein 2 (NAB2). The human NAB2 gene can be found at Entrez gene #4665. The human NAB2 protein can be found at Uniprot ID Q15742. In some embodiments, the plurality of factors comprises Colony stimulating factor 1 receptor (CSF1R). The human CSF1R gene can be found at Entrez gene #1436. The human CSF1R protein can be found at Uniprot ID P07333. In some embodiments, the plurality of factors comprises CCN family member 4 (CCN4). The human CCN4 gene can be found at Entrez gene #8840. The human CCN4 protein can be found at Uniprot ID 095388. In some embodiments, the plurality of factors comprises Glycerol-3- phosphate dehydrogenase 1-like protein (GPD1). The human GPD1 gene can be found at Entrez gene #2819. The human GPD1 protein can be found at Uniprot ID Q8N335. In some embodiments, the plurality of factors comprises Pro state- specific antigen also known as kallikrein-3 (KLK3). The human KLK3 gene can be found at Entrez gene #354. The human KLK3 protein can be found at Uniprot ID P07288. In some embodiments, the plurality of factors comprises Chemokine (C-X-C motif) ligand 13 (CXCL13). The human CXCL13 gene can be found at Entrez gene #10563. The human CXCL13 protein can be found at Uniprot ID 043927. In some embodiments, the plurality of factors comprises Granzyme A (GZMA). The human GZMA gene can be found at Entrez gene #3001. The human GZMA protein can be found at Uniprot ID P12544. In some embodiments, the plurality of factors comprises Complement component C9 (C9). The human C9 gene can be found at Entrez gene #735. The human C9 protein can be found at Uniprot ID P06683. In some embodiments, the plurality of factors comprises Rapl GTPase-activating protein 1 (RAP 1 GAP). The human RAP 1 GAP gene can be found at Entrez gene #5909. The human RAP1GAP protein can be found at Uniprot ID P47736. In some embodiments, the plurality of factors comprises Insulin-like growth factor-binding protein 1 (IGFBP1). The humanIGFBP1 gene can be found at Entrez gene #3484. The human IGFBP1 protein can be found at Uniprot ID PO8833. In some embodiments, the plurality of factors comprises Probable ATP-dependent RNA helicase DHX58 (DHX58). The human DHX58 gene can be found at Entrez gene #79132. The human DHX58 protein can be found at Uniprot ID Q96C10. In some embodiments, the plurality of factors comprises COP9 signalosome complex subunit 2 (COPS2). The human COPS2 gene can be found at Entrez gene #9318. The human COPS2 protein can be found at Uniprot ID P61201. In some embodiments, the plurality of factors comprises Interleukin- 1 receptor accessory protein (IL1RAP). The human IL1RAP gene can be found at Entrez gene #3556. The human IL 1 RAP protein can be found at Uniprot ID Q9NPH3. In some embodiments, the plurality of factors comprises Chemokine (C-C motif) ligand 25 (CCL25). The human CCL25 gene can be found at Entrez gene #6370. The human CCL25 protein can be found at Uniprot ID Q68A93. In some embodiments, the plurality of factors comprises Hemopexin (HPX). The human HPX gene can be found at Entrez gene #3263. The human HPX protein can be found at Uniprot ID P02790. In some embodiments, the plurality of factors comprises Adrenomedullin (ADM). The human ADM gene can be found at Entrez gene #133. The human ADM protein can be found at Uniprot ID P35318. In some embodiments, the plurality of factors comprises Cluster of Differentiation 93 (CD93). The human CD93 gene can be found at Entrez gene #22918. The human CD93 protein can be found at Uniprot ID Q9NPY3. In some embodiments, the plurality of factors comprises Interferon- stimulated gene 15 (ISG15). The human ISG15 gene can be found at Entrez gene #9636. The human ISG15 pro tein c an be found at Uniprot ID P05161. In some embodiments, the plurality of factors comprises Myosin light chain 6B (MYL6B). The human MYL6B gene can be found at Entrez gene #140465. The human MYL6B protein can be found at Uniprot ID P14649. In some embodiments, the plurality of factors comprises Heat shock 70 kDa protein 1 (HSPA1A). The human HSPA1A gene can be found at Entrez gene #3303. The human HSPA1A protein can be found at Uniprot ID P0DMV8 or P0DMV9. In some embodiments, the plurality of factors comprises Methyl-CpG-binding domain protein 1 (MBD1). The human MBD1 gene can be found at Entrez gene #4152. The human MBD1 protein can be found at Uniprot ID Q9UIS9. In some embodiments, the plurality of factors comprises Trafficking protein particle complex subunit 3 (TRAPPC3). The human TRAPPC3 gene can be found at Entrez gene #27095. The human TRAPPC3 protein can be found at Uniprot ID 043617. In some embodiments, the plurality of factors comprises RAC- beta serine / threonine-protein kinase (AKT2). The human AKT2 gene can be found at Entrez gene #208. The human AKT2 protein can be found at Uniprot ID P31751. In someembodiments, the plurality of factors comprises Cytokine receptor-like factor 1 (CRLF1). The human CRLF1 gene can be found at Entrez gene #9244. The human CRLF1 protein can be found at Uniprot ID 075462. In some embodiments, the plurality of factors comprises Ferritin light chain (FTL). The human FTL gene can be found at Entrez gene #2512. The human FTL protein can be found at Uniprot ID P02792. In some embodiments, the plurality of factors comprises Histone-binding protein RBBP4 (RBBP4). The human RBBP4 gene can be found at Entrez gene #5928. The human RBBP4 protein can be found at Uniprot ID Q09028. In some embodiments, the plurality of factors comprises BMP binding endothelial regulator (BMPER). The human BMPER gene can be found at Entrez gene #168667. The human BMPER protein can be found at Uniprot ID Q8N8U9. In some embodiments, the plurality of factors comprises Maspin (SERPINB5). The human SERPINB5 gene can be found at Entrez gene #5268. The human SERPINB5 protein can be found at Uniprot ID P36952. In some embodiments, the plurality of factors comprises Myelin P2 protein (PMP2). The human PMP2 gene can be found at Entrez gene #5375. The human PMP2 protein can be found at Uniprot ID P02689. In some embodiments, the plurality of factors comprises Ornithine transcarbamylase (OTC). The human OTC gene can be found at Entrez gene #5009. The human OTC protein can be found at Uniprot ID P00480. In some embodiments, the plurality of factors comprises Otoraplin (OTOR). The human OTOR gene can be found at Entrez gene #56914. The human OTOR protein can be found at Uniprot ID Q9NRC9. In some embodiments, the plurality of factors comprises Diamine oxidase [copper-containing] (AOC1). The human AOC1 gene can be found at Entrez gene #26. The human AOC1 protein can be found at Uniprot ID Q8JZQ5. In some embodiments, the plurality of factors comprises Fibroblast growth factor-binding protein 1 (FGFBP1). The human FGFBP1 gene can be found at Entrez gene #9982. The human FGFBP1 protein can be found at Uniprot ID Q14512. In some embodiments, the plurality of factors comprises Attractin (ATRN). The human ATRN gene can be found at Entrez gene #8455. The human ATRN protein can be found at Uniprot ID 075882. In some embodiments, the plurality of factors comprises N- acetylglucosaminidase, alpha (NAGLU). The human NAGLU gene can be found at Entrez gene #4669. The human NAGLU protein can be found at Uniprot ID P54802. In some embodiments, the plurality of factors comprises Serum amyloid Al (SAA1). The human SAA1 gene can be found at Entrez gene #6288. The human SAA1 protein can be found at Uniprot ID P0DJI8. In some embodiments, the plurality of factors comprises Serum amyloid A4 (SAA4). The human SAA4 gene can be found at Entrez gene #6291. The human SAA4 protein can be found at Uniprot ID P35542. In some embodiments, the plurality of factorscomprises Calsyntenin- 1 (CLSTN1). The human CLSTN 1 gene can be found at Entrez gene #22883. The human CLSTN1 protein can be found at Uniprot ID 094985. In some embodiments, the plurality of factors comprises Glutathione synthetase (GSS). The human GSS gene can be found at Entrez gene #2937. The human GSS protein can be found at Uniprot ID C9K4X8. In some embodiments, the plurality of factors comprises Dihydrolipoamide dehydrogenase (DLD). The human DLD gene can be found at Entrez gene #1738. The human DLD protein can be found at Uniprot ID P09622. In some embodiments, the plurality of factors comprises Ephrin type-B receptor 4 (EPHB4). The human EPHB4 gene can be found at Entrez gene #2050. The human EPHB4 protein can be found at Uniprot ID P54760. In some embodiments, the plurality of factors comprisesSerine protease 27 (PRSS27). The human PRSS27 gene can be found at Entrez gene #83886. The human PRSS27 protein can be found at Uniprot ID Q9BQR3. In some embodiments, the plurality of factors comprises Mucin- 16 (MUC16). The human MUC16 gene can be found at Entrez gene #94025. The human MUC16 protein can be found at Uniprot ID Q8WXI7. In some embodiments, the plurality of factors comprises Complement factor H-related protein 2 (CFHR2). The human CFHR2 gene can be found at Entrez gene #3080. The human CFHR2 protein can be found at Uniprot ID P36980. In some embodiments, the plurality of factors comprises Serine protease HTRA1 (HTRA1). The human HTRA1 gene can be found at Entrez gene #5654. The human HTRA1 protein can be found at Uniprot ID Q92743. In some embodiments, the plurality of factors comprises Keratin, type I cytoskeletal 19 (KRT19). The human KRT19 gene can be found at Entrez gene #3880. The human KRT19 protein can be found at Uniprot ID P08727. In some embodiments, the plurality of factors comprises Retinol binding protein 4 (RBP4). The human RBP4 gene can be found at Entrez gene #5950. The human RBP4 protein can be found at Uniprot ID P02753. In some embodiments, the plurality of factors comprises SPARC-related modular calcium-binding protein 2 (SMOC2). The human SMOC2 gene can be found at Entrez gene #64094. The human SMOC2 protein can be found at Uniprot ID Q9H3U7. In some embodiments, the plurality of factors comprises Biotinidase (BTD). The human BTD gene can be found at Entrez gene #686. The human BTD protein can be found at Uniprot ID P43251. In some embodiments, the plurality of factors comprises Alpha-taxili (TXLNA). The human TXLNA gene can be found at Entrez gene #200081. The human TXLNA protein can be found at Uniprot ID P40222. In some embodiments, the plurality of factors comprises Marginal zone B and B 1 cell-specific protein (MZB 1). The human MZB 1 gene can be found at Entrez gene #51237. The human MZB 1 protein can be found at UniprotID Q8WU39. In some embodiments, the plurality of factors comprises FAS -associated death domain protein (FADD). The human FADD gene can be found at Entrez gene #8772. The human FADD protein can be found at Uniprot ID Q13158. In some embodiments, the plurality of factors comprises Gelsolin (GSN). The human GSN gene can be found at Entrez gene #2934. The human GSN protein can be found at Uniprot ID P06396. In some embodiments, the plurality of factors comprises Cadherin-17 (CDH17). The human CDH17 gene can be found at Entrez gene #1015. The human CDH17 protein can be found at Uniprot ID Q12864. In some embodiments, the plurality of factors comprises Leukocyte cell-derived chemotaxin-2 (LECT2). The human LECT2 gene can be found at Entrez gene #3950. The human LECT2 protein can be found at Uniprot ID 014960. In some embodiments, the plurality of factors comprises ADAMTS-like protein 1 (ADAMTSL1). The human ADAMTSL1 gene can be found at Entrez gene #92949. The human ADAMTSL1 protein can be found at Uniprot ID Q8N6G6. In some embodiments, the plurality of factors comprises Ribonuclease T2 (RNASET2). The human RNASET2 gene can be found at Entrez gene #8635. The human RNASET2 protein can be found at Uniprot ID 000584. In some embodiments, the plurality of factors comprises Semaphorin-4A (SEMA4A). The human SEMA4A gene can be found at Entrez gene #64218. The human SEMA4A protein can be found at Uniprot ID Q9H3S1. In some embodiments, the plurality of factors comprises Dolichyl-diphosphooligo saccharide — protein glycosyltransferase 48 kDa subunit (DDOST). The human DDOST gene can be found at Entrez gene #1650. The human DDOST protein can be found at Uniprot ID P39656. In some embodiments, the plurality of factors comprises Dehydrogenase / reductase SDR family member 6 (BDH2 / DHRS6). The human BDH2 gene can be found at Entrez gene # 56898. The human BDH2 protein can be found at Uniprot ID Q9BUT1. In some embodiments, the plurality of factors comprises U2 small nuclear ribonucleoprotein B (SNRPB2). The human SNRPB2 gene can be found at Entrez gene #6629. The human SNRPB2 protein can be found at Uniprot ID P08579. In some embodiments, the plurality of factors comprises Golgi membrane protein 1 (G0LM1). The human G0LM1 gene can be found at Entrez gene #51280. The human G0LM1 protein can be found at Uniprot ID Q8NBJ4. In some embodiments, the plurality of factors comprises Ras-related protein Rab-3A (RAB3A). The human RAB3A gene can be found at Entrez gene #5864. The human RAB3A protein can be found at Uniprot ID P20336. In some embodiments, the plurality of factors comprises CD46 complement regulatory protein (CD46). The human CD46 gene can be found at Entrez gene #4179. The human CD46 protein can be found at Uniprot ID P15529. In some embodiments, the plurality of factorscomprises Septin-6 (SEPTIN6). The human SEPTIN6 gene can be found at Entrez gene # 23157. The human SEPTIN6 protein can be found at Uniprot ID Q3SZN0. In some embodiments, the plurality of factors comprises WW domain-containing oxidoreductase (WWOX). The human WWOX gene can be found at Entrez gene #51741. The human WWOX protein can be found at Uniprot ID Q9NZC7. In some embodiments, the plurality of factors comprises WD repeat-containing protein 5 (WDR5). The human WDR5 gene can be found at Entrez gene #11091. The human WDR5 protein can be found at Uniprot ID P61964. In some embodiments, the plurality of factors comprises Hippocalcin-like protein 1 (HPCAL1). The human HPCAL1 gene can be found at Entrez gene #3241. The human HPCAL1 protein can be found at Uniprot ID P37235. In some embodiments, the plurality of factors comprises Aldehyde dehydrogenase 5 family, member Al (ALDH5A1). The human ALDH5A1 gene can be found at Entrez gene #7915. The human ALDH5A1 protein can be found at Uniprot ID P51649. In some embodiments, the plurality of factors comprises Synaptic vesicle membrane protein VAT-1 homolog (VAT1). The human VAT1 gene can be found at Entrez gene #10493. The human VAT1 protein can be found at Uniprot ID Q99536. In some embodiments, the plurality of factors comprises Cytoplasmic seryl-tRNA synthetase (SARS1). The human SARS1 gene can be found at Entrez gene #6301. The human SARS1 protein can be found at Uniprot ID P49591. In some embodiments, the plurality of factors comprises Afamin (AFM). The human AFM gene can be found at Entrez gene #173. The human AFM protein can be found at Uniprot ID P43652. In some embodiments, the plurality of factors comprises Cytidine deaminase (CD A). The human CD A gene can be found at Entrez gene #978. The human CD A protein can be found at Uniprot ID P32320. In some embodiments, the plurality of factors comprises Intelectin-1 (ITLN1). The human ITLN1 gene can be found at Entrez gene #55600. The human ITLN1 protein can be found at Uniprot ID Q8WWA0. In some embodiments, the plurality of factors comprises Leucine-rich repeats and immunoglobulin-like domains protein 1 (LRIG1). The human LRIG1 gene can be found at Entrez gene #26018. The human LRIG1 protein can be found at Uniprot ID Q96JA1. In some embodiments, the plurality of factors comprises Gremlin (GREM1). The human GREM1 gene can be found at Entrez gene #26585. The human GREM1 protein can be found at Uniprot ID 060565. In some embodiments, the plurality of factors comprises Prostaglandin reductase 2 (PTGR2). The human PTGR2 gene can be found at Entrez gene # 145482. The human PTGR2 protein can be found at Uniprot ID Q8N8N7. In some embodiments, the plurality of factors comprises Gelsolin (GSN). The human GSN gene can be found at Entrez gene #2934. The human GSN protein can be foundat Uniprot ID P06396. In some embodiments, the plurality of factors comprises Ubiquitin / ISG15-conjugating enzyme E2 L6 (UBE2L6). The human UBE2L6 gene can be found at Entrez gene #9246. The human UBE2L6 protein can be found at Uniprot ID 014933. In some embodiments, the plurality of factors comprises Clathrin light chain A (CLTA). The human CLTA gene can be found at Entrez gene # 1211. The human CLTA protein can be found at Uniprot ID P09496. In some embodiments, the plurality of factors comprises Glutathione reductase (GSR). The human GSR gene can be found at Entrez gene #2936. The human GSR protein can be found at Uniprot ID P00390. In some embodiments, the plurality of factors comprises Programmed cell death protein 6 (PDCD6). The human PDCD6 gene can be found at Entrez gene #10016. The human PDCD6 protein can be found at Uniprot ID 075340. In some embodiments, the plurality of factors comprises Gamma- synuclein (SNCG). The human SNCG gene can be found at Entrez gene #6623. The human SNCG protein can be found at Uniprot ID 076070. In some embodiments, the plurality of factors comprises Corticotropin-releasing hormone (CRH). The human CRH gene can be found at Entrez gene #1392. The human CRH protein can be found at Uniprot ID P06850. In some embodiments, the plurality of factors comprises Regulator of G-protein signaling 21 (RGS21). The human RGS21 gene can be found at Entrez gene # 431704. The human RGS21 protein can be found at Uniprot ID Q2M5E4. In some embodiments, the plurality of factors comprises Ubiquitin-conjugating enzyme E2 R2 (UBE2R2). The human UBE2R2 gene can be found at Entrez gene #54926. The human UBE2R2 protein can be found at Uniprot ID Q712K3. In some embodiments, the plurality of factors comprises Brain acid soluble protein 1 (BASP1). The human BASP1 gene can be found at Entrez gene #10409. The human BASP1 protein can be found at Uniprot ID P80723. In some embodiments, the plurality of factors comprises Guanylate binding protein 5 (GBP5). The human GBP5 gene can be found at Entrez gene #115362. The human GBP5 protein can be found at Uniprot ID Q96PP8. In some embodiments, the plurality of factors comprises Lamin B2 (LMNB2). The human LMNB2 gene can be found at Entrez gene #84823. The human LMNB2 protein can be found at Uniprot ID Q03252. In some embodiments, the plurality of factors comprises Ribonuclease P protein subunit p20 (POP7). The human POP7 gene can be found at Entrez gene #10248. The human POP7 protein can be found at Uniprot ID 075817. In some embodiments, the plurality of factors comprises Retinoic acid early transcript IL (RAET1L / ULBP6). The human RAET1L gene can be found at Entrez gene # 154064. The human RAET1L protein can be found at Uniprot ID Q5VY80. In some embodiments, the plurality of factors comprises Semaphorin-5B (SEMA5B). The human SEMA5B gene canbe found at Entrez gene # 54437. The human SEMA5B protein can be found at Uniprot ID Q9P283. In some embodiments, the plurality of factors comprises Contactin-3 (CNTN3). The human CNTN3 gene can be found at Entrez gene #5067. The human CNTN3 protein can be found at Uniprot ID Q9P232. In some embodiments, the plurality of factors comprises Ubiquitin-like protein 3 (UBL3). The human UBL3 gene can be found at Entrez gene #5412. The human UBL3 protein can be found at Uniprot ID 095164. In some embodiments, the plurality of factors comprises Methylmalonic aciduria and homocystinuria type C protein (MMACHC). The human MMACHC gene can be found at Entrez gene #25974. The human MMACHC protein can be found at Uniprot ID Q9Y4U 1. In some embodiments, the plurality of factors comprises Transcription factor II B (GTF2B). The human GTF2B gene can be found at Entrez gene #2959. The human GTF2B protein can be found at Uniprot ID Q00403. In some embodiments, the plurality of factors comprises GTP cyclohydrolase 1 feedback regulatory protein (GCHFR). The human GCHFR gene can be found at Entrez gene #2644. The human GCHFR protein can be found at Uniprot ID P30047. In some embodiments, the plurality of factors comprises Protein LRATD2 (LRATD2). The human LRATD2 gene can be found at Entrez gene # 157638. The human LRATD2 protein can be found at Uniprot ID Q96KN1. In some embodiments, the plurality of factors comprises Serine / threonine -protein kinase Sgkl (SGK1). The human SGK1 gene can be found at Entrez gene #6446. The human SGK1 protein can be found at Uniprot ID 000141. In some embodiments, the plurality of factors comprises tRNA-splicing endonuclease subunit Senl5 (TSEN15). The human TSEN15 gene can be found at Entrez gene #116461. The human TSEN15 protein can be found at Uniprot ID Q8WW01. In some embodiments, the plurality of factors comprises GTP-binding protein SARlb (SAR1B). The human SAR1B gene can be found at Entrez gene # 51128. The human SAR1B protein can be found at Uniprot ID Q9Y6B6. In some embodiments, the plurality of factors comprises CDK5 regulatory subunit- associated protein 3 (CDK5RAP3). The human CDK5RAP3 gene can be found at Entrez gene #80279. The human CDK5RAP3 protein can be found at Uniprot ID Q96JB5. In some embodiments, the plurality of factors comprises HAUS augmin-like complex subunit 1 (HAUS1). The human HAUS 1 gene can be found at Entrez gene #115106. The human HAUS 1 protein can be found at Uniprot ID Q96CS2. In some embodiments, the plurality of factors comprises NF-kappa- B inhibitor alpha (NKIRAS1). The human NKIRAS1 gene can be found at Entrez gene # 28512. The human NKIRAS1 protein can be found at Uniprot ID P25963. In some embodiments, the plurality of factors comprises Pyridoxal phosphate phosphatase PHOSPHO2 (PHOSPHO2). The human PHOSPHO2 gene can be found at Entrez gene #493911. The human PH0SPH02 protein can be found at Uniprot ID Q8TCD6. In some embodiments, the plurality of factors comprises Protocadherin-17 (PCDH17). The human PCDH17 gene can be found at Entrez gene #27253. The human PCDH17 protein can be found at Uniprot ID 014917. In some embodiments, the plurality of factors comprises Tripartite motif-containing protein 5 (TRIM5). The human TRIM5 gene can be found at Entrez gene #85363. The human TRIM5 protein can be found at Uniprot ID Q9C035. In some embodiments, the plurality of factors comprises Aldehyde dehydrogenase 7 family, member Al (ALDH7A1). The human ALDH7A1 gene can be found at Entrez gene #501. The human ALDH7A1 protein can be found at Uniprot ID P49419. In some embodiments, the plurality of factors comprises Thioredoxin-like protein 4 A (TXNL4A). The human TXNL4A gene can be found at Entrez gene #10907. The human TXNL4A protein can be found at Uniprot ID P83876. In some embodiments, the plurality of factors comprises Centrosomal protein 20 (CEP20). The human CEP20 gene can be found at Entrez gene # 123811. The human CEP20 protein can be found at Uniprot ID Q96NB 1. In some embodiments, the plurality of factors comprises Calcium / calmodulin-dependent 3',5'-cyclic nucleotide phosphodiesterase IB (PDE1B). The human PDE1B gene can be found at Entrez gene #5153. The human PDE1B protein can be found at Uniprot ID Q01064. In some embodiments, the plurality of factors comprises Integrin alpha-4 (ITGA4). The human ITGA4 gene can be found at Entrez gene #3676. The human ITGA4 protein can be found at Uniprot ID P13612. In some embodiments, the plurality of factors comprises Integrin beta- 1 (ITGB 1). The human ITGB 1 gene can be found at Entrez gene #3688. The human ITGB 1 protein can be found at Uniprot ID P05556. In some embodiments, the plurality of factors comprises Leucine-rich repeat neuronal protein 3 (LRFN3). The human LRFN3 gene can be found at Entrez gene #54674. The human LRFN3 protein can be found at Uniprot ID Q9H3W5. In some embodiments, the plurality of factors comprises Adhesion G protein- coupled receptor B 1 (ADGRB 1). The human ADGRB 1 gene can be found at Entrez gene # 575. The human ADGRB 1 protein can be found at Uniprot ID 014514. In some embodiments, the plurality of factors comprises N- sulphoglucosamine sulphohydrolase (SGSH). The human SGSH gene can be found at Entrez gene #6448. The human SGSH protein can be found at Uniprot ID P51688. In some embodiments, the plurality of factors comprises Alpha- 1,6-mannosylglycoprotein 6-beta-N-acetylglucosaminyltransferase A (MGAT5). The human MGAT5 gene can be found at Entrez gene #4249. The human MGAT5 protein can be found at Uniprot ID Q09328. In some embodiments, the plurality of factors comprises 3-beta-glucuronosyltransferase 1 (B3GAT1). The human B3GAT1 genecan be found at Entrez gene #27087. The human B3GAT1 protei In some embodiments, the plurality of factors comprises Alpha- 1,6-mannosylglycoprotein 6-beta-N- acetylglucosaminyltransferase A (MGAT5). The human MGAT5 gene can be found at Entrez gene #4249. The human MGAT5 protein can be found at Uniprot ID P06396. n can be found at Uniprot ID Q09328. In some embodiments, the plurality of factors comprises Fibulin-7 (FBLN7). The human FBLN7 gene can be found at Entrez gene # 129804. The human FBLN7 protein can be found at Uniprot ID Q501P1. In some embodiments, the plurality of factors comprises Amyloid beta A4 precursor protein-binding family B member 1 -interacting protein ( APB B IIP). The human APB B IIP gene can be found at Entrez gene #54518. The human APBB1IP protein can be found at Uniprot ID Q7Z5R6. In some embodiments, the plurality of factors comprises Serum paraoxonase / arylesterase 2 (PON2). The human PON2 gene can be found at Entrez gene #5445. The human PON2 protein can be found at Uniprot ID Q15165. In some embodiments, the plurality of factors comprises Serine / threonine-protein phosphatase 2A 56 kDa regulatory subunit delta isoform (PPP2R5D). The human PPP2R5D gene can be found at Entrez gene #5528. The human PPP2R5D protein can be found at Uniprot ID Q14738. In some embodiments, the plurality of factors comprises Fox-1 homolog A (RBFOX1). The human RBFOX1 gene can be found at Entrez gene #54715. The human RBFOX1 protein can be found at Uniprot ID Q9NWB1. In some embodiments, the plurality of factors comprises TIMP metallopeptidase inhibitor 1 (TIMP1). The human TIMP1 gene can be found at Entrez gene #7076. The human TIMP1 protein can be found at Uniprot ID PO1O33 or Q6FGX5. In some embodiments, the plurality of factors comprises Gem-associated protein 7 (GEMIN7). The human GEMIN7 gene can be found at Entrez gene #79760. The human GEMIN7 protein can be found at Uniprot ID Q9H840. In some embodiments, the plurality of factors comprises Casein kinase I isoform alpha (CSNK1A1L). The human CSNK1A1L gene can be found at Entrez gene #1452. The human CSNK1A1L protein can be found at Uniprot ID P48729. In some embodiments, the plurality of factors comprises PHD finger protein 11 (PHF11). The human PHF11 gene can be found at Entrez gene # 51131. The human PHF11 protein can be found at Uniprot ID Q9UIL8. In some embodiments, the plurality of factors comprises Butyrophilin subfamily 2 member A2 (BTN2A2). The human BTN2A2 gene can be found at Entrez gene #10385. The human BTN2A2 protein can be found at Uniprot ID Q8WVV5. In some embodiments, the plurality of factors comprises S-phase kinase-associated protein 2 (SKP2). The human SKP2 gene can be found at Entrez gene #6502. The human SKP2 protein can be found at Uniprot ID Q13309. In some embodiments, the plurality of factors comprises Spermatogenesis-associated protein 16 (SPATA46). The human SPATA46 gene can be found at Entrez gene # 284680. The human SPATA46 protein can be found at Uniprot ID Q5T0L3. In some embodiments, the plurality of factors comprises Lin-7 homolog A (LIN7A). The human LIN7A gene can be found at Entrez gene #8825. The human LIN7A protein can be found at Uniprot ID 014910. In some embodiments, the plurality of factors comprises BLOC-1- related complex subunit 5 (BORCS5). The human BORCS5 gene can be found at Entrez gene # 118426. The human BORCS5 protein can be found at Uniprot ID Q969J3. In some embodiments, the plurality of factors comprises Arrestin domain-containing protein 5 (ARRDC5). The human ARRDC5 gene can be found at Entrez gene # 645432. The human ARRDC5 protein can be found at Uniprot ID A6NEK1. In some embodiments, the plurality of factors comprises Choline-phosphate cytidylyltransferase A (PCYT1A). The human PCYT1A gene can be found at Entrez gene #5130. The human PCYT1A protein can be found at Uniprot ID P49585. In some embodiments, the plurality of factors comprises Phytanoyl-CoA dioxygenase, peroxisomal (PHYH). The human PHYH gene can be found at Entrez gene # 5264. The human PHYH protein can be found at Uniprot ID 014832. In some embodiments, the plurality of factors comprises Ankyrin repeat domaincontaining protein 63 (ANKRD63). The human ANKRD63 gene can be found at Entrez gene #100131244. The human ANKRD63 protein can be found at Uniprot ID C9JTQ0. In some embodiments, the plurality of factors comprises Variable charge X-linked protein 1 (VCX). The human VCX gene can be found at Entrez gene #26609. The human VCX protein can be found at Uniprot ID Q9H320. In some embodiments, the plurality of factors comprises Protein N-terminal asparagine amidohydrolase (NTAN1). The human NTAN1 gene can be found at Entrez gene #123803. The human NTAN1 protein can be found at Uniprot ID Q96AB6. In some embodiments, the plurality of factors comprises StAR-related lipid transfer domain protein 7 (STARD7). The human STARD7 gene can be found at Entrez gene #56910. The human STARD7 protein can be found at Uniprot ID Q9NQZ5. In some embodiments, the plurality of factors comprises Apolipoprotein L2 (APOL2). The human APOL2 gene can be found at Entrez gene #23780. The human APOL2 protein can be found at Uniprot ID Q9BQE5. In some embodiments, the plurality of factors comprises Fms- related tyrosine kinase 4 (FLT4). The human FLT4 gene can be found at Entrez gene #2324. The human FLT4 protein can be found at Uniprot ID P35916. In some embodiments, the plurality of factors comprises CapZ-interacting protein (RCSD1). The human RCSD1 gene can be found at Entrez gene # 92241. The human RCSD1 protein can be found at Uniprot ID Q6JBY9. In some embodiments, the plurality of factors comprises INTS3 and NABPinteracting protein (INIP). The human INIP gene can be found at Entrez gene # 58493. The human INIP protein can be found at Uniprot ID Q9NRY2. In some embodiments, the plurality of factors comprises Vimentin-type intermediate filament-associated coiled-coil protein (VMAC). The human VMAC gene can be found at Entrez gene # 400673. The human VMAC protein can be found at Uniprot ID Q2NL98. In some embodiments, the plurality of factors comprises Xaa-Pro aminopeptidase 3 (XPNPEP3). The human XPNPEP3 gene can be found at Entrez gene # 63929. The human XPNPEP3 protein can be found at Uniprot ID Q9NQH7. In some embodiments, the plurality of factors comprises Interferon epsilon (IFNE). The human IFNE gene can be found at Entrez gene # 338376. The human IFNE protein can be found at Uniprot ID Q80ZF2. In some embodiments, the plurality of factors comprises Negative elongation factor A (NELFA). The human NELFA gene can be found at Entrez gene # 7469. The human NELFA protein can be found at Uniprot ID Q8BG3O. In some embodiments, the plurality of factors comprises Lysine demethylase 8 (KDM8). The human KDM8 gene can be found at Entrez gene #79831. The human KDM8 protein can be found at Uniprot ID Q8N371. In some embodiments, the plurality of factors comprises Nuclear cap-binding protein complex (NCBP1). The human NCBP1 gene can be found at Entrez gene #4686. The human NCBP1 protein can be found at Uniprot ID Q56A27. In some embodiments, the plurality of factors comprises Upstream stimulatory factor 2 (USF2). The human USF2 gene can be found at Entrez gene #7392. The human USF2 protein can be found at Uniprot ID Q15853. In some embodiments, the plurality of factors comprises Leucine-rich repeat-containing protein 75 A (LRRC75A). The human LRRC75A gene can be found at Entrez gene # 388341. The human LRRC75A protein can be found at Uniprot ID Q8NAA5. In some embodiments, the plurality of factors comprises Amyloid P component, serum (APCS). The human APCS gene can be found at Entrez gene #325. The human APCS protein can be found at Uniprot ID P02743. In some embodiments, the plurality of factors comprises l-Phosphatidylinositol-4,5-bisphosphate phosphodiesterase delta-1 (PLCD1). The human PLCD1 gene can be found at Entrez gene #5333. The human PLCD1 protein can be found at Uniprot ID P51178. In some embodiments, the plurality of factors comprises Espin (ESPN). The human ESPN gene can be found at Entrez gene #83715. The human ESPN protein can be found at Uniprot ID Q5JYL1. In some embodiments, the plurality of factors comprises DNA-binding protein RFX5 (RFX5). The human RFX5 gene can be found at Entrez gene #5993. The human RFX5 protein can be found at Uniprot ID P48382. In some embodiments, the plurality of factors comprises Ribosomal protein S6 kinase beta-2 (RPS6KB2). The human RPS6KB2 gene canbe found at Entrez gene # 6199. The human RPS6KB2 protein can be found at Uniprot ID Q9UBS0. In some embodiments, the plurality of factors comprises BOS complex subunit N0M02 (N0M02). The human N0M02 gene can be found at Entrez gene # 283820. The human N0M02 protein can be found at Uniprot ID Q5JPE7. In some embodiments, the plurality of factors comprises Transcription elongation factor A protein-like 2 (TCEAL2). The human TCEAL2 gene can be found at Entrez gene # 140597. The human TCEAL2 protein can be found at Uniprot ID Q9H3H9. In some embodiments, the plurality of factors comprises Carboxylesterase 3 (CES3). The human CES3 gene can be found at Entrez gene #23491. The human CES3 protein can be found at Uniprot ID Q6UWW8. In some embodiments, the plurality of factors comprises Dual specificity tyrosine-phosphorylation- regulated kinase 1A (DYRK1A). The human DYRK1A gene can be found at Entrez gene #1859. The human DYRK1A protein can be found at Uniprot ID Q13627. In some embodiments, the plurality of factors comprises Cytochrome P450 2C19 (CYP2C19). The human CYP2C19 gene can be found at Entrez gene #1557. The human CYP2C19 protein can be found at Uniprot ID P33261. In some embodiments, the plurality of factors comprises Complement factor I (CFI). The human CFI gene can be found at Entrez gene #3426. The human CFI protein can be found at Uniprot ID P05156 or Q8WW88. In some embodiments, the plurality of factors comprises Insulin-like growth factor-binding protein 3 (IGFBP3). The human IGFBP3 gene can be found at Entrez gene #3486. The human IGFBP3 protein can be found at Uniprot ID P17936. In some embodiments, the plurality of factors comprises Interleukin 6 (IL6). The human IL6 gene can be found at Entrez gene #3569. The human IL6 protein can be found at Uniprot ID P05231. In some embodiments, the plurality of factors comprises Leptin (LEP). The human LEP gene can be found at Entrez gene #3952. The human LEP protein can be found at Uniprot ID P41159. In some embodiments, the plurality of factors comprises CREB -regulated transcription coactivator 3 (CRTC3). The human CRTC3 gene can be found at Entrez gene #64784. The human CRTC3 protein can be found at Uniprot ID Q6UUV7. In some embodiments, the plurality of factors comprises Vascular endothelial growth factor A (VEGFA). The human VEGFA gene can be found at Entrez gene #7422. The human VEGFA protein can be found at Uniprot ID Pl 5692. In some embodiments, the plurality of factors comprises Interleukin- 1 receptor accessory protein (IL 1 RAP). The human IL 1 RAP gene can be found at Entrez gene #3556. The human IL1RAP protein can be found at Uniprot ID Q9NPH3. In some embodiments, the plurality of factors comprises Hepatocyte growth factor (HGF). The human HGF gene can be found at Entrez gene #3082. The human HGF protein can be found at Uniprot ID P14210. In someembodiments, the plurality of factors comprises Phospholipase A2, membrane associated (PLA2G2A). The human PLA2G2A gene can be found at Entrez gene #5320. The human PLA2G2A protein can be found at Uniprot ID P14555. In some embodiments, the plurality of factors comprises Chemokine (C-C motif) ligand 25 (CCL25). The human CCL25 gene can be found at Entrez gene #6370. The human CCL25 protein can be found at Uniprot ID 015444. In some embodiments, the plurality of factors comprises serpin family A member 7 (SERPINA7). The human SERPINA7 gene can be found at Entrez gene #6906. The human SERPINA7 protein can be found at Uniprot ID P05543. In some embodiments, the plurality of factors comprises Cytochrome P450 reductase (POR). The human POR gene can be found at Entrez gene #5447. The human POR protein can be found at Uniprot ID P16435. In some embodiments, the plurality of factors comprises CCN family member 3 (CCN3). The human CCN3 gene can be found at Entrez gene # 4856. The human CCN3 protein can be found at Uniprot ID P48745. In some embodiments, the plurality of factors comprises Hemopexin (HPX). The human HPX gene can be found at Entrez gene #3263. The human HPX protein can be found at Uniprot ID P02790. In some embodiments, the plurality of factors comprises Insulin-like growth factor-binding protein 1 (IGFBP1). The human IGFBP1 gene can be found at Entrez gene #3484. The human IGFBP1 protein can be found at Uniprot ID PO8833. In some embodiments, the plurality of factors comprises matrix metalloproteinase-3 (MMP3). The human MMP3 gene can be found at Entrez gene #4314. The human MMP3 protein can be found at Uniprot ID P08254. In some embodiments, the plurality of factors comprises Fibrinogen alpha chain (FGA). The human FGA gene can be found at Entrez gene #2243. The human FGA protein can be found at Uniprot ID P02671. In some embodiments, the plurality of factors comprises Fibrinogen beta chain (FGB). The human FGB gene can be found at Entrez gene #2244. The human FGB protein can be found at Uniprot ID P02675. In some embodiments, the plurality of factors comprises Fibrinogen gamma chain (FGG). The human FGG gene can be found at Entrez gene #2266. The human FGG protein can be found at Uniprot ID P02679. In some embodiments, the plurality of factors comprises Basal cell adhesion molecule (BCAM). The human BCAM gene can be found at Entrez gene #4059. The human BCAM protein can be found at Uniprot ID P50895. In some embodiments, the plurality of factors comprises Kunitz-type protease inhibitor 1 (SPINT1). The human SPINT1 gene can be found at Entrez gene #6692. The human SPINT1 protein can be found at Uniprot ID 043278. In some embodiments, the plurality of factors comprises Histone acetyltransferase 1 (HAT1). The human HAT1 gene can be found at Entrez gene #8520. The human HAT1 protein can be found at Uniprot ID 014929. In someembodiments, the plurality of factors comprises Growth hormone receptor (GHR). The human GHR gene can be found at Entrez gene #2690. The human GHR protein can be found at Uniprot ID P10912. In some embodiments, the plurality of factors comprises Complement factor properdin (CFP). The human CFP gene can be found at Entrez gene #5199. The human CFP protein can be found at Uniprot ID P27918. In some embodiments, the plurality of factors comprises Contactin 1 (CNTN1). The human CNTN1 gene can be found at Entrez gene #1272. The human CNTN1 protein can be found at Uniprot ID Q12860. In some embodiments, the plurality of factors comprises Alpha 2-antiplasmin (SERPINF2). The human SERPINF2 gene can be found at Entrez gene #5345. The human SERPINF2 protein can be found at Uniprot ID P08697. In some embodiments, the plurality of factors comprises Interleukin- 19 (IL 19). The human IL 19 gene can be found at Entrez gene # 29949. The human IL19 protein can be found at Uniprot ID Q9UHD0. In some embodiments, the plurality of factors comprises Myoglobin (MB). The human MB gene can be found at Entrez gene #4151. The human MB protein can be found at Uniprot ID P02144. In some embodiments, the plurality of factors comprises Immunoglobulin heavy constant mu (IGHM). The human IGHM gene can be found at Entrez gene #3507. The human IGHM protein can be found at Uniprot ID P01871. In some embodiments, the plurality of factors comprises Lipopolysaccharide binding protein (LBP). The human LBP gene can be found at Entrez gene #3929. The human LBP protein can be found at Uniprot ID P18428. In some embodiments, the plurality of factors comprises N-acylethanolamine-hydrolyzing acid amidase (NAAA). The human NAAA gene can be found at Entrez gene # 27163. The human NAAA protein can be found at Uniprot ID Q02083. In some embodiments, the plurality of factors comprises Hyaluronan and proteoglycan link protein 1 (HAPLN1). The human HAPLN1 gene can be found at Entrez gene #1404. The human HAPLN1 protein can be found at Uniprot ID P10915. In some embodiments, the plurality of factors comprises Iduronate 2-sulfatase (IDS). The human IDS gene can be found at Entrez gene # 3423. The human IDS protein can be found at Uniprot ID P22304. In some embodiments, the plurality of factors comprises Nidogen-1 (NIDI). The human NIDI gene can be found at Entrez gene #4811. The human NIDI protein can be found at Uniprot ID P14543. In some embodiments, the plurality of factors comprises Aggrecan (ACAN). The human ACAN gene can be found at Entrez gene #176. The human ACAN protein can be found at Uniprot ID P16112. In some embodiments, the plurality of factors comprises Transforming growth factor, beta-induced, 68kDa (TGFBI). The human TGFBI gene can be found at Entrez gene #7045. The human TGFBI protein can be found at Uniprot ID Q15582. In some embodiments, the plurality offactors comprises Delta- like 4 (DLL4). The human DLL4 gene can be found at Entrez gene #54567. The human DLL4 protein can be found at Uniprot ID Q9NR61. In some embodiments, the plurality of factors comprises Fc fragment of IgG, low affinity Illb, receptor (FCGR3B). The human FCGR3B gene can be found at Entrez gene #2215. The human FCGR3B protein can be found at Uniprot ID 075015. In some embodiments, the plurality of factors comprises Aminoacylase-1 (ACY1). The human ACY1 gene can be found at Entrez gene #95. The human ACY 1 protein can be found at Uniprot ID Q03154. In some embodiments, the plurality of factors comprises Integrin binding sialoprotein (IBSP). The human IBSP gene can be found at Entrez gene #3381. The human IBSP protein can be found at Uniprot ID P21815. In some embodiments, the plurality of factors comprises Kallistatin (SERPINA4). The human SERPINA4 gene can be found at Entrez gene #5267. The human SERPINA4 protein can be found at Uniprot ID P29622. In some embodiments, the plurality of factors comprises Periostin (POSTN). The human POSTN gene can be found at Entrez gene #10631. The human POSTN protein can be found at Uniprot ID Q15063. In some embodiments, the plurality of factors comprises E-selectin (SELE). The human SELE gene can be found at Entrez gene #6401. The human SELE protein can be found at Uniprot ID P16581. In some embodiments, the plurality of factors comprises P2 microglobulin (B2M). The human B2M gene can be found at Entrez gene #567. The human B2M protein can be found at Uniprot ID P61769. In some embodiments, the plurality of factors comprises Hepcidin (HAMP). The human HAMP gene can be found at Entrez gene #57817. The human HAMP protein can be found at Uniprot ID P81172. In some embodiments, the plurality of factors comprises Alpha-1 antitrypsin (SERPINA1). The human SERPINA1 gene can be found at Entrez gene #5265. The human SERPINA1 protein can be found at Uniprot ID P01009. In some embodiments, the plurality of factors comprises alpha-2-HS -glycoprotein (AHSG). The human AHSG gene can be found at Entrez gene #197. The human AHSG protein can be found at Uniprot ID P02765. In some embodiments, the plurality of factors comprises Brain-type creatine kinase (CKB). The human CKB gene can be found at Entrez gene #1152. The human CKB protein can be found at Uniprot ID P12277. In some embodiments, the plurality of factors comprises Creatine kinase, muscle (CKM). The human CKM gene can be found at Entrez gene #1158. The human CKM protein can be found at Uniprot ID P06732. In some embodiments, the plurality of factors comprises Protein C (PROC). The human PROC gene can be found at Entrez gene #5624. The human PROC protein can be found at Uniprot ID P04070. In some embodiments, the plurality of factors comprises Angiopoietin-like 4 (ANGPTL4). The human ANGPTL4 gene can be found atEntrez gene #51129. The human ANGPTL4 protein can be found at Uniprot ID Q9BY76. In some embodiments, the plurality of factors comprises Methyl-CpG-binding domain protein 4 (MBD4). The human MBD4 gene can be found at Entrez gene #8930. The human MBD4 protein can be found at Uniprot ID 095243. In some embodiments, the plurality of factors comprises 26S proteasome non-ATPase regulatory subunit 7 (PSMD7). The human PSMD7 gene can be found at Entrez gene #5713. The human PSMD7 protein can be found at Uniprot ID P51665. In some embodiments, the plurality of factors comprises immunoglobulin heavy constant epsilon (IGHE). The human IGHE gene can be found at Entrez gene #3497. The human IGHE protein can be found at Uniprot ID P01854. In some embodiments, the plurality of factors comprises C-X-C motif chemokine ligand 10 (CXCL10). The human CXCL10 gene can be found at Entrez gene #3627. The human CXCL10 protein can be found at Uniprot ID P02778. In some embodiments, the plurality of factors comprises Plasma kallikrein (KLKB1). The human KLKB1 gene can be found at Entrez gene # 3818. The human KLKB1 protein can be found at Uniprot ID P03952. In some embodiments, the plurality of factors comprises Factor H (CFH). The human CFH gene can be found at Entrez gene #3075. The human CFH protein can be found at Uniprot ID P08603. In some embodiments, the plurality of factors comprises Prefoldin subunit 5 (PFDN5). The human PFDN5 gene can be found at Entrez gene #5204. The human PFDN5 protein can be found at Uniprot ID Q99471. In some embodiments, the plurality of factors comprises RNA- binding protein 39 (RBM39). The human RBM39 gene can be found at Entrez gene #9584 The human RBM39 protein can be found at Uniprot ID Q14498 or Q5QP23. In some embodiments, the plurality of factors comprises dCTP pyrophosphatase 1 (DCTPP1). The human DCTPP1 gene can be found at Entrez gene #79077. The human DCTPP1 protein can be found at Uniprot ID Q9H773. In some embodiments, the plurality of factors comprises Brain-specific serine protease 4 (PRSS22). The human PRSS22 gene can be found at Entrez gene #64063. The human PRSS22 protein can be found at Uniprot ID Q9GZN4. In some embodiments, the plurality of factors compriseskynureninase (KYNU). The human KYNU gene can be found at Entrez gene # 8942. The human KYNU protein can be found at Uniprot ID QI 6719. In some embodiments, the plurality of factors comprises Interleukin 6 (IL6). The human IL6 gene can be found at Entrez gene #3569. The human IL6 protein can be found at Uniprot ID P05231. In some embodiments, the plurality of factors comprises Transcortin (SERPINA6). The human SERPINA6 gene can be found at Entrez gene #866. The human SERPINA6 protein can be found at Uniprot ID P08185. In some embodiments, the plurality of factors comprises Inter-alpha-trypsin inhibitor heavy chain H4 (ITIH4). Thehuman ITIH4 gene can be found at Entrez gene # 3700. The human ITIH4 protein can be found at Uniprot ID Q 14624. In some embodiments, the plurality of factors comprises Stratifin (SFN). The human SFN gene can be found at Entrez gene # 2810. The human SFN protein can be found at Uniprot ID P31947. In some embodiments, the plurality of factors comprises Chemokine (C-C motif) ligand 7 (CCE7). The human CCE7 gene can be found at Entrez gene # 6354. The human CCL7 protein can be found at Uniprot ID P80098. In some embodiments, the plurality of factors comprises Lysozyme C (LYZ). The human LYZ gene can be found at Entrez gene # 4069. The human LYZ protein can be found at Uniprot ID P61626. In some embodiments, the plurality of factors comprises Collagenase 3 (MMP13). The human MMP13 gene can be found at Entrez gene # 4322. The human MMP13 protein can be found at Uniprot ID P45452. In some embodiments, the plurality of factors comprises Stanniocalcin-1 (STC1). The human STC1 gene can be found at Entrez gene # 6781. The human STC1 protein can be found at Uniprot ID P52823. In some embodiments, the plurality of factors comprises Macrophage-capping protein (CAPG). The human CAPG gene can be found at Entrez gene # 822. The human CAPG protein can be found at Uniprot ID P40121. In some embodiments, the plurality of factors comprises Peptidase inhibitor 3 (PI3). The human PI3 gene can be found at Entrez gene # 5266. The human PI3 protein can be found at Uniprot ID P19957. In some embodiments, the plurality of factors comprises Glypican-5 (GPC5). The human GPC5 gene can be found at Entrez gene #2262. The human GPC5 protein can be found at Uniprot ID P78333. In some embodiments, the plurality of factors comprises Histidine-rich glycoprotein (HRG). The human HRG gene can be found at Entrez gene # 3273. The human HRG protein can be found at Uniprot ID P04196. In some embodiments, the plurality of factors comprises Mammaglobin-B (SCGB2A1). The human SCGB2A1 gene can be found at Entrez gene # 4246. The human SCGB2A1 protein can be found at Uniprot ID 075556. In some embodiments, the plurality of factors comprises NAD-dependent deacetylase sirtuin 2 (SIRT2). The human SIRT2 gene can be found at Entrez gene # 22933. The human SIRT2 protein can be found at Uniprot ID Q8IXJ6. In some embodiments, the plurality of factors comprises Tumor necrosis factor-inducible gene 6 protein (TNFAIP6). The human TNFAIP6 gene can be found at Entrez gene #2934. The human TNFAIP6 protein can be found at Uniprot ID P06396. In some embodiments, the plurality of factors comprises CMRF35-like molecule 6 (CD300C). The human CD300C gene can be found at Entrez gene # 10871. The human CD300C protein can be found at Uniprot ID Q08708. In some embodiments, the plurality of factors comprises Transmembrane glycoprotein NMB(GPNMB). The human GPNMB gene can be found at Entrez gene # 10457. The human GPNMB protein can be found at Uniprot ID QI 4956. In some embodiments, the plurality of factors comprises Keratin 18 (KRT18). The human KRT18 gene can be found at Entrez gene # 3875. The human KRT18 protein can be found at Uniprot ID P05783. In some embodiments, the plurality of factors comprises tumor necrosis factor superfamily member 14 (TNFSF14). The human TNFSF14 gene can be found at Entrez gene # 8740. The human TNFSF14 protein can be found at Uniprot ID 043557. In some embodiments, the plurality of factors comprises Leptin receptor (LEPR). The human LEPR gene can be found at Entrez gene # 3953. The human LEPR protein can be found at Uniprot ID P48357. In some embodiments, the plurality of factors comprises Protein kinase C gamma type (PRKCG). The human PRKCG gene can be found at Entrez gene # 5582. The human PRKCG protein can be found at Uniprot ID P05129. In some embodiments, the plurality of factors comprises Fibrinogen-like protein 1 (FGL1). The human FGLlgene can be found at Entrez gene # 2267. The human FGL1 protein can be found at Uniprot ID QO883O. In some embodiments, the plurality of factors comprises Peptidoglycan recognition protein 2 (PGLYRP2). The human PGLYRP2 gene can be found at Entrez gene # 114770. The human PGLYRP2 protein can be found at Uniprot ID Q96PD5. In some embodiments, the plurality of factors comprises Neuropeptide FF (NPFF). The human NPFF gene can be found at Entrez gene # 8620. The human NPFF protein can be found at Uniprot ID 015130. In some embodiments, the plurality of factors comprises Microfibril-associated glycoprotein 4 (MFAP4). The human MFAP4 gene can be found at Entrez gene # 4239. The human MFAP4 protein can be found at Uniprot ID P55083. In some embodiments, the plurality of factors comprises Protein disulfide-isomerase (TMX3). The human TMX3 gene can be found at Entrez gene # 54495. The human TMX3 protein can be found at Uniprot ID Q96JJ7. In some embodiments, the plurality of factors comprises Glucosidase 2 subunit beta (PRKCSH). The human PRKCSH gene can be found at Entrez gene # 5589. The human PRKCSH protein can be found at Uniprot ID P14314. In some embodiments, the plurality of factors comprises Beta- def ensin 112 (DEFB 112). The human DEFB 112 gene can be found at Entrez gene # 245915. The human DEFB 112 protein can be found at Uniprot ID Q3OKQ8. In some embodiments, the plurality of factors comprises Semaphorin-4D (SEMA4D). The human SEMA4D gene can be found at Entrez gene # 10507. The human SEMA4D protein can be found at Uniprot ID Q92854. In some embodiments, the plurality of factors comprises Lysophosphatidic acid phosphatase type 6 (ACP6). The human ACP6 gene can be found at Entrez gene # 51205. The human ACP6 protein can be found at Uniprot ID Q9NPH0. In some embodiments, theplurality of factors comprises Alpha-fetoprotein (AFP). The human AFP gene can be found at Entrez gene # 174. The human AFP protein can be found at Uniprot ID P02771. In some embodiments, the plurality of factors comprises Nerve growth factor (NGF). The human NGF gene can be found at Entrez gene # 4803. The human NGF protein can be found at Uniprot ID PO1138. In some embodiments, the plurality of factors comprises Ferritin heavy chain (FTH1). The human FTH1 gene can be found at Entrez gene # 2495. The human FTH1 protein can be found at Uniprot ID P02794. In some embodiments, the plurality of factors comprises Ferritin light chain (FTL). The human FTL gene can be found at Entrez gene # 2512. The human FTL protein can be found at Uniprot ID P02792. In some embodiments, the plurality of factors comprises Dermokine (DMKN). The human DMKN gene can be found at Entrez gene # 93099. The human DMKN protein can be found at Uniprot ID Q6E0U4. In some embodiments, the plurality of factors comprises EPH receptor A10 (EPHA10). The human EPHA10 gene can be found at Entrez gene # 284656. The human EPHA10 protein can be found at Uniprot ID Q5JZY3. In some embodiments, the plurality of factors comprises Chordin-like protein 2 (CHRDL2). The human CHRDL2 gene can be found at Entrez gene # 25884. The human CHRDL2 protein can be found at Uniprot ID Q6WN34. In some embodiments, the plurality of factors comprises Tumor protein P53 (TP53). The human TP53 gene can be found at Entrez gene # 7157. The human TP53 protein can be found at Uniprot ID P04637. In some embodiments, the plurality of factors comprises Diamine oxidase [copper-containing] (AOC1). The human AOC1 gene can be found at Entrez gene # 26. The human AOC1 protein can be found at Uniprot ID P19801. In some embodiments, the plurality of factors comprises Interferon alpha-8 (IFNA8). The human IFNA8 gene can be found at Entrez gene # 3445. The human IFNA8 protein can be found at Uniprot ID P32881. In some embodiments, the plurality of factors comprises Chorionic somatomammotropin hormone 1 (CSH1). The human CSH1 gene can be found at Entrez gene # 1442. The human CSH1 protein can be found at Uniprot ID P0DML2. In some embodiments, the plurality of factors comprises Chorionic somatomammotropin hormone 2 (CSH2). The human CSH2 gene can be found at Entrez gene # 1443. The human CSH2 protein can be found at Uniprot ID P0DML3. In some embodiments, the plurality of factors comprises Tenascin C (TNC). The human TNC gene can be found at Entrez gene # 3371. The human TNC protein can be found at Uniprot ID P24821. In some embodiments, the plurality of factors comprises Phospholipid transfer protein (PLTP). The human PLTP gene can be found at Entrez gene # 5360. The human PLTP protein can be found at Uniprot ID P55058. In some embodiments, the plurality of factors comprises cellular communicationnetwork factor 1 (CCN1). The human CCN1 gene can be found at Entrez gene # 3491. The human CCN1 protein can be found at Uniprot ID 000622. In some embodiments, the plurality of factors comprises Calsyntenin-3 (CLSTN3). The human CLSTN3 gene can be found at Entrez gene # 9746. The human CLSTN3 protein can be found at Uniprot ID Q9BQT9. In some embodiments, the plurality of factors comprises Oncoprotein-induced transcript 3 protein (OIT3). The human OIT3 gene can be found at Entrez gene # 170392. The human OIT3 protein can be found at Uniprot ID Q8WWZ8. In some embodiments, the plurality of factors comprises Inactive glutathione hydrolase 2 (GGT2). The human GGT2 gene can be found at Entrez gene # 102724197. The human GGT2 protein can be found at Uniprot ID P36268. In some embodiments, the plurality of factors comprises Fibromodulin (FMOD). The human FMOD gene can be found at Entrez gene # 2331. The human FMOD protein can be found at Uniprot ID Q06828. In some embodiments, the plurality of factors comprises Putative uncharacterized protein IRX2-DT (C5orf38). The human C5orf38 gene can be found at Entrez gene # 153571. The human C5orf38 protein can be found at Uniprot ID Q86SI9. In some embodiments, the plurality of factors comprises von Willebrand factor A domain-containing protein 1 (VWA1). The human VWA1 gene can be found at Entrez gene # 64856. The human VWA1 protein can be found at Uniprot ID Q6PCB0. In some embodiments, the plurality of factors comprises Inhibin beta C chain (INHBC). The human INHBC gene can be found at Entrez gene # 3626. The human INHBC protein can be found at Uniprot ID P55103. In some embodiments, the plurality of factors comprises Adhesion G protein-coupled receptor F5 (ADGRF5). The human ADGRF5 gene can be found at Entrez gene # 221395. The human ADGRF5 protein can be found at Uniprot ID Q8IZF2. In some embodiments, the plurality of factors comprises Complement Clq-like protein 2 (C1QL2). The human C1QL2 gene can be found at Entrez gene # 165257. The human C1QL2 protein can be found at Uniprot ID Q7Z5L3. In some embodiments, the plurality of factors comprises Prenylcysteine oxidase 1 (PCYOX1). The human PCYOX1 gene can be found at Entrez gene # 51449. The human PCYOX1 protein can be found at Uniprot ID Q9UHG3. In some embodiments, the plurality of factors comprises Amine oxidase, copper containing 2 (AOC2). The human AOC2 gene can be found at Entrez gene # 314. The human AOC2 protein can be found at Uniprot ID 075106. In some embodiments, the plurality of factors comprises Complement factor H-related protein 4 (CFHR4). The human CFHR4 gene can be found at Entrez gene # 10877. The human CFHR4 protein can be found at Uniprot ID Q92496. In some embodiments, the plurality of factors comprises Leucine rich repeat containing 15 (LRRC15). The human LRRC15 gene can be found at Entrez gene # 131578.The human LRRC15 protein can be found at Uniprot ID Q8TF66. In some embodiments, the plurality of factors comprises Periostin (POSTN). The human POSTN gene can be found at Entrez gene # 10631. The human POSTN protein can be found at Uniprot ID Q15063. In some embodiments, the plurality of factors comprises Ubiquitin-conjugating enzyme E2 JI (UBE2J1). The human UBE2J1 gene can be found at Entrez gene # 51465. The human UBE2J1 protein can be found at Uniprot ID Q9Y385. In some embodiments, the plurality of factors comprises GDNF family receptor alpha like (GFRAL). The human GFRAL gene can be found at Entrez gene # 389400. The human GFRAL protein can be found at Uniprot ID Q6UXV0. In some embodiments, the plurality of factors comprises Insulin-like growth factor 2 (IGF2). The human IGF2 gene can be found at Entrez gene # 3481. The human IGF2 protein can be found at Uniprot ID P01344. In some embodiments, the plurality of factors comprises Leukocyte immunoglobulin-like receptor subfamily B member 5 (LILRB5). The human LILRB5 gene can be found at Entrez gene # 10990. The human LILRB5 protein can be found at Uniprot ID 075023. In some embodiments, the plurality of factors comprises Leukocyte immunoglobulin-like receptor subfamily A member 6 (LILRA6). The human LILRA6 gene can be found at Entrez gene # 79168. The human LILRA6 protein can be found at Uniprot ID Q6PI73. In some embodiments, the plurality of factors comprises Apolipoprotein A-II (APOA2). The human APOA2 gene can be found at Entrez gene # 336. The human APOA2 protein can be found at Uniprot ID P02652. In some embodiments, the plurality of factors comprises von Willebrand factor A domain-containing protein 2 (VWA2). The human VWA2 gene can be found at Entrez gene # 340706. The human VWA2 protein can be found at Uniprot ID Q5GFL6. In some embodiments, the plurality of factors comprises Protein DEPPI (DEPPI). The human DEPPI gene can be found at Entrez gene # 11067. The human DEPPI protein can be found at Uniprot ID Q9NTK1. In some embodiments, the plurality of factors comprises Complement Clq tumor necrosis factor- related protein 3 (C1QTNF3). The human C1QTNF3 gene can be found at Entrez gene # 114899. The human C1QTNF3 protein can be found at Uniprot ID Q9BXJ4. In some embodiments, the plurality of factors comprises Serpin A9 (SERPINA9). The human SERPINA9 gene can be found at Entrez gene # 327657. The human SERPINA9 protein can be found at Uniprot ID Q86WD7. In some embodiments, the plurality of factors comprises Complement factor H-related protein 5 (CFHR5). The human CFHR5 gene can be found at Entrez gene # 81494. The human CFHR5 protein can be found at Uniprot ID Q9BXR6. In some embodiments, the plurality of factors comprises Disks large homolog 3 (DLG3). The human DLG3 gene can be found at Entrez gene # 1741. The human DLG3 protein can befound at Uniprot ID Q92796. In some embodiments, the plurality of factors comprises Glycolipid transfer protein domain-containing protein 2 (GLTPD2). The human GLTPD2 gene can be found at Entrez gene # 388323. The human GLTPD2 protein can be found at Uniprot ID A6NH11. In some embodiments, the plurality of factors comprises Hemoglobin subunit theta-1 (HBQ1). The human HBQ1 gene can be found at Entrez gene # 3049. The human HBQ1 protein can be found at Uniprot ID P09105. In some embodiments, the plurality of factors comprises Ectonucleoside triphosphate diphosphohydrolase- 1 (ENTPD1). The human ENTPD1 gene can be found at Entrez gene # 953. The human ENTPD1 protein can be found at Uniprot ID P49961. In some embodiments, the plurality of factors comprises Angiogenic factor with G patch and FHA domains 1 (AGGF1). The human AGGF1 gene can be found at Entrez gene # 55109. The human AGGF1 protein can be found at Uniprot ID Q8N302. In some embodiments, the plurality of factors comprises Neuregulin 2 (NRG2). The human NRG2 gene can be found at Entrez gene # 9542. The human NRG2 protein can be found at Uniprot ID 014511. In some embodiments, the plurality of factors comprises Spondin 2 (SPON2). The human SPON2 gene can be found at Entrez gene # 10417. The human SPON2 protein can be found at Uniprot ID Q9BUD6. In some embodiments, the plurality of factors comprises Protein FAM241B (FAM241B). The human FAM241B gene can be found at Entrez gene # 219738. The human FAM241B protein can be found at Uniprot ID Q96D05. In some embodiments, the plurality of factors comprises Junctional adhesion molecule-like (J AML). The human JAML gene can be found at Entrez gene # 120425. The human JAML protein can be found at Uniprot ID Q86YT9. In some embodiments, the plurality of factors comprises Butyrylcholinesterase (BCHE). The human BCHE gene can be found at Entrez gene # 590. The human BCHE protein can be found at Uniprot ID P06276. In some embodiments, the plurality of factors comprises Transmembrane glycoprotein NMB (GPNMB). The human GPNMB gene can be found at Entrez gene # 10457. The human GPNMB protein can be found at Uniprot ID Q14956. In some embodiments, the plurality of factors comprises Apolipoprotein D (APOD). The human APOD gene can be found at Entrez gene # 347. The human APOD protein can be found at Uniprot ID P05090. In some embodiments, the plurality of factors comprises Deltalike protein 1 (DLL1). The human DLL1 gene can be found at Entrez gene # 28514. The human DLL1 protein can be found at Uniprot ID 000548. In some embodiments, the plurality of factors comprises Platelet endothelial aggregation receptor 1 (PEAR1). The human PEAR1 gene can be found at Entrez gene # 375033. The human PEAR1 protein can be found at Uniprot ID Q5VY43. In some embodiments, the plurality of factors comprisesR-spondin-4 (RSPO4). The human RSPO4 gene can be found at Entrez gene # 343637. The human RSPO4 protein can be found at Uniprot ID Q2I0M5. In some embodiments, the plurality of factors comprises Leptin (LEP). The human LEP gene can be found at Entrez gene # 3952. The human LEP protein can be found at Uniprot ID P41159. In some embodiments, the plurality of factors comprises ADP-ribosylation factor-like protein 8B (ARL8B). The human ARL8B gene can be found at Entrez gene # 55207. The human ARL8B protein can be found at Uniprot ID Q9NVJ2. In some embodiments, the plurality of factors comprises Protocadherin-10 (PCDH10). The human PCDH10 gene can be found at Entrez gene # 57575. The human PCDH10 protein can be found at Uniprot ID Q9P2E7. In some embodiments, the plurality of factors comprises Microfibrillar-associated protein 3- like (MFAP3L). The human MFAP3L gene can be found at Entrez gene # 9848. The human MFAP3L protein can be found at Uniprot ID 075121. In some embodiments, the plurality of factors comprises Monocyte differentiation antigen CD 14 (CD 14). The human CD 14 gene can be found at Entrez gene # 929. The human CD 14 protein can be found at Uniprot ID P08571. In some embodiments, the plurality of factors comprises Collagen alpha-l(XV) chain (COL15A1). The human COL15A1 gene can be found at Entrez gene # 1306. The human COL15A1 protein can be found at Uniprot ID P39059. In some embodiments, the plurality of factors comprises Hepatitis A virus cellular receptor 1 (HAVCR1). The human HAVCR1 gene can be found at Entrez gene # 26762. The human HAVCR1 protein can be found at Uniprot ID Q96D42. In some embodiments, the plurality of factors comprises Rho guanine nucleotide exchange factor 10 (ARHGEF10). The human ARHGEF10 gene can be found at Entrez gene # 9639. The human ARHGEF10 protein can be found at Uniprot ID 015013. In some embodiments, the plurality of factors comprises Mannosyl-oligosaccharide 1,2-alpha-mannosidase IB (MAN1A2). The human MAN1A2 gene can be found at Entrez gene # 10905. The human MAN1A2 protein can be found at Uniprot ID 060476. In some embodiments, the plurality of factors comprises Quinone oxidoreductase-like protein 1 (CRYZL1). The human CRYZL1 gene can be found at Entrez gene # 9946. The human CRYZL1 protein can be found at Uniprot ID 095825. In some embodiments, the plurality of factors comprises Tissue factor pathway inhibitor 2 (TFPI2). The human TFPI2 gene can be found at Entrez gene # 7980. The human TFPI2 protein can be found at Uniprot ID P48307. In some embodiments, the plurality of factors comprises Plexin domain-containing protein 1 (PLXDC1). The human PLXDC1 gene can be found at Entrez gene # 57125. The human PLXDC1 protein can be found at Uniprot ID Q8IUK5. In some embodiments, the plurality of factors comprises Lysosomal acid phosphatase (ACP2). The human ACP2 genecan be found at Entrez gene # 53. The human ACP2 protein can be found at Uniprot ID Pl 1117. In some embodiments, the plurality of factors comprises Biotinidase (BTD). The human BTD gene can be found at Entrez gene # 686. The human BTD protein can be found at Uniprot ID P43251. In some embodiments, the plurality of factors comprises Microfibrillar- associated protein 2 (MFAP2). The human MFAP2 gene can be found at Entrez gene # 4237. The human MFAP2 protein can be found at Uniprot ID P55001. In some embodiments, the plurality of factors comprises Inter-alpha-trypsin inhibitor heavy chain H2 (ITIH2). The human ITIH2 gene can be found at Entrez gene # 3698. The human ITIH2 protein can be found at Uniprot ID P19823. In some embodiments, the plurality of factors comprises EF-hand calcium-binding domain-containing protein 14 (EFCAB14). The human EFCAB 14 gene can be found at Entrez gene # 9813. The human EFCAB 14 protein can be found at Uniprot ID 075071. In some embodiments, the plurality of factors comprises Phospholipase Al member A (PLA1A). The human PLA1A gene can be found at Entrez gene # 51365. The human PLA1A protein can be found at Uniprot ID Q53H76. In some embodiments, the plurality of factors comprises Granzyme K (GZMK). The human GZMK gene can be found at Entrez gene # 3003. The human GZMK protein can be found at Uniprot ID P49863. In some embodiments, the plurality of factors comprises Y box binding protein 1 (YBX1). The human YBX1 gene can be found at Entrez gene # 4904. The human YBX1 protein can be found at Uniprot ID P67809. In some embodiments, the plurality of factors comprises Indoleamine-pyrrole 2,3-dioxygenase (IDO1). The human IDO1 gene can be found at Entrez gene # 3620. The human IDO1 protein can be found at Uniprot ID P14902. In some embodiments, the plurality of factors comprises NAD(P)H dehydrogenase [quinone] 1 (NQO1). The human NQO1 gene can be found at Entrez gene # 1728. The human NQO1 protein can be found at Uniprot ID P15559. In some embodiments, the plurality of factors comprises Testican-3 (SPOCK3). The human SPOCK3 gene can be found at Entrez gene # 50859. The human SPOCK3 protein can be found at Uniprot ID Q9BQ16. In some embodiments, the plurality of factors comprises NTF2-related export protein 1 (NXT1). The human NXT1 gene can be found at Entrez gene # 29107. The human NXT1 protein can be found at Uniprot ID Q9UKK6.

[0147] In some embodiments, the factor is selected from a factor provided in Table 6. In some embodiments, the plurality of factors is selected from the factors provided in Table 6. In some embodiments, the plurality of factors comprises at least two factors selected from those provided in Table 6. In some embodiments, the plurality of factors consists of factorsselected from Table 6. In some embodiments, the factor is selected from the group consisting of: CYCS, SERPINA4, G0LM1, TIMP1, FUT5, ANGPT2, ITM2B, INHBA, CD5, THBS2, SAA2, CETP, KLK3|SERPINA3, C5, SERPINA5, MLN, TMEM132A, HDGFL3, TAC1, CXCL12, S100A12, L0XL2, RBP4, NFKBID, HHIP, UMOD, CDON, CPLX3, AD AMTS 1, KHSRP, SERPING1, ECM1, CTSB, NDUFV2, TLR5, SET, DNAJB9, APOL1, FGFR4, TNFRSF11B, KLK14, PLG, DCBLD1, BGLAP, CARNMT1, MAPKAPK5, RCAN2, CST2, INHBB, APOA4, IL17RB, BOC, LRP1, SF3B4, PARVA, GALNT9, CFHR1, WFDC1, BCHE, NT5DC3, PAPPA, KLK14, ANGPTL3, TBK1, CXCL8, LRP8, CCNH, CFB, and PPP4R3A. The amino acid sequences of these factors can be found in the Uniprot database, for example, and each factor’s Uniprot accession number is provided herein. Further, methods, reagents, and assays for measuring expression levels of these factors are well known in the art and are commercially available.

[0148] In some embodiments, the plurality of factors comprises A disintegrin and metalloproteinase with thrombospondin motifs 1 (AD AMTS 1). The human AD AMTS 1 gene can be found at Entrez gene # 9510. The human AD AMTS 1 protein can be found at Uniprot ID Q9UHI8. In some embodiments, the plurality of factors comprises Angiopoietin- 2 (ANGPT2). The human ANGPT2 gene can be found at Entrez gene # 285. The human ANGPT2 protein can be found at Uniprot ID 015123. In some embodiments, the plurality of factors comprises Angiopoietin-related protein 3 (ANGPTL3). The human ANGPTL3 gene can be found at Entrez gene # 27329. The human ANGPTL3 protein can be found at Uniprot ID Q9Y5C1. In some embodiments, the plurality of factors comprises Apolipoprotein A-IV (APOA4). The human APOA4 gene can be found at Entrez gene # 337. The human APOA4 protein can be found at Uniprot ID P06727. In some embodiments, the plurality of factors comprises Apolipoprotein LI (APOL1). The human APOL1 gene can be found at Entrez gene # 8542. The human APOL1 protein can be found at Uniprot ID 014791. In some embodiments, the plurality of factors comprises Cholinesterase (BCHE). The human BCHE gene can be found at Entrez gene # 590. The human BCHE protein can be found at Uniprot ID P06276. In some embodiments, the plurality of factors comprises Osteocalcin (BGLAP). The human BGLAP gene can be found at Entrez gene # 632. The human BGLAP protein can be found at Uniprot ID P02818. In some embodiments, the plurality of factors comprises Brother of CDO (BOC). The human BOC gene can be found at Entrez gene # 91653. The human BOC protein can be found at Uniprot ID Q9BWV1. In some embodiments, the plurality of factors comprises Complement C5 (C5). The human C5gene can be found at Entrez gene # 727. The human C5 protein can be found at Uniprot ID P01031. In some embodiments, the plurality of factors comprises Carnosine N- methyltransferase (CARNMT1). The human CARNMT1 gene can be found at Entrez gene # 138199. The human CARNMT1 protein can be found at Uniprot ID Q8N4J0. In some embodiments, the plurality of factors comprises Cyclin-H (CCNH). The human CCNH gene can be found at Entrez gene # 902. The human CCNH protein can be found at Uniprot ID P51946. In some embodiments, the plurality of factors comprises CD5 (CD5). The human CD5 gene can be found at Entrez gene # 921. The human CD5 protein can be found at Uniprot ID P06127. In some embodiments, the plurality of factors comprises Cell adhesion molecule-related / down-regulated by oncogenes (CDON). The human CDON gene can be found at Entrez gene # 50937. The human CDON protein can be found at Uniprot ID Q4KMG0. In some embodiments, the plurality of factors comprises Cholesteryl ester transfer protein (CETP). The human CETP gene can be found at Entrez gene # 1071. The human CETP protein can be found at Uniprot ID Pl 1597. In some embodiments, the plurality of factors comprises Complement factor B (CFB). The human CFB gene can be found at Entrez gene # 629. The human CFB protein can be found at Uniprot ID P00751. In some embodiments, the plurality of factors comprises Complement factor H-related protein 1 (CFHR1). The human CFHR1 gene can be found at Entrez gene # 3078. The human CFHR1 protein can be found at Uniprot ID Q03591. In some embodiments, the plurality of factors comprises Complexin-3 (CPLX3). The human CPLX3 gene can be found at Entrez gene # 594855. The human CPLX3 protein can be found at Uniprot ID Q8WVH0. In some embodiments, the plurality of factors comprises Cystatin-SA (CST2). The human CST2 gene can be found at Entrez gene # 1470. The human CST2 protein can be found at Uniprot ID P09228. In some embodiments, the plurality of factors comprises Cathepsin B (CTSB). The human CTSB gene can be found at Entrez gene # 1508. The human CTSB protein can be found at Uniprot ID P07858. In some embodiments, the plurality of factors comprises Stromal cell-derived factor 1 (CXCL12). The human CXCL12 gene can be found at Entrez gene # 6387. The human CXCL12 protein can be found at Uniprot ID P48061. In some embodiments, the plurality of factors comprises Interleukin-8 (CXCL8). The human CXCL8 gene can be found at Entrez gene # 3576. The human CXCL8 protein can be found at Uniprot ID P10145. In some embodiments, the plurality of factors comprises Cytochrome c (CYCS). The human CYCS gene can be found at Entrez gene # 54205. The human CYCS protein can be found at Uniprot ID P99999. In some embodiments, the plurality of factors comprises Discoidin, CUB and LCCL domain-containing protein 1 (DCBLD1). The human DCBLD1gene can be found at Entrez gene # 285761. The human DCBLD1 protein can be found at Uniprot ID Q8N8Z6. In some embodiments, the plurality of factors comprises DnaJ homolog subfamily B member 9 (DNAJB9). The human DNAJB9 gene can be found at Entrez gene # 4189. The human DNAJB9 protein can be found at Uniprot ID Q9UBS3. In some embodiments, the plurality of factors comprises Extracellular matrix protein 1 (ECM1). The human ECM1 gene can be found at Entrez gene # 1893. The human ECM1 protein can be found at Uniprot ID Q16610. In some embodiments, the plurality of factors comprises Fibroblast growth factor receptor 4 (FGFR4). The human FGFR4 gene can be found at Entrez gene # 2264. The human FGFR4 protein can be found at Uniprot ID P22455. In some embodiments, the plurality of factors comprises Alpha-(l,3)-fucosyltransferase 5 (FUT5). The human FUT5 gene can be found at Entrez gene # 2527. The human FUT5 protein can be found at Uniprot ID QI 1128. In some embodiments, the plurality of factors comprises Polypeptide N-acetylgalactosaminyltransferase 9 (GALNT9). The human GALNT9 gene can be found at Entrez gene # 50614. The human GALNT9 protein can be found at Uniprot ID Q9HCQ5. In some embodiments, the plurality of factors comprises Golgi membrane protein 1 (G0LM1). The human G0LM1 gene can be found at Entrez gene # 51280. The human G0LM1 protein can be found at Uniprot ID Q8NBJ4. In some embodiments, the plurality of factors comprises Hepatoma-derived growth factor-related protein 3 (HDGFL3). The human HDGFL3 gene can be found at Entrez gene # 50810. The human HDGFL3 protein can be found at Uniprot ID Q9Y3E1. In some embodiments, the plurality of factors comprises Hedgehog -interacting protein (HHIP). The human HHIP gene can be found at Entrez gene # 64399. The human HHIP protein can be found at Uniprot ID Q96QV1. In some embodiments, the plurality of factors comprises Interleukin- 17 receptor B (IL17RB). The human IL17RB gene can be found at Entrez gene # 55540. The human IL17RB protein can be found at Uniprot ID Q9NRM6. In some embodiments, the plurality of factors comprises Inhibin beta A chain (INHBA). The human INHBA gene can be found at Entrez gene # 3624. The human INHBA protein can be found at Uniprot ID P08476. In some embodiments, the plurality of factors comprises Inhibin beta B chain (INHBB). The human INHBB gene can be found at Entrez gene # 3625. The human INHBB protein can be found at Uniprot ID P09529. In some embodiments, the plurality of factors comprises Integral membrane protein 2B (ITM2B). The human ITM2B gene can be found at Entrez gene # 9445. The human ITM2B protein can be found at Uniprot ID Q9Y287. In some embodiments, the plurality of factors comprises Far upstream element-binding protein 2 (KHSRP). The human KHSRP gene can be found at Entrez gene # 8570. The human KHSRPprotein can be found at Uniprot ID Q92945. In some embodiments, the plurality of factors comprises Kallikrein-14 (KLK14). The human KLK14 gene can be found at Entrez gene # 43847. The human KLK14 protein can be found at Uniprot ID Q9P0G3. In some embodiments, the plurality of factors comprises Prostate-specific antigen (KLK3). The human KLK3gene can be found at Entrez gene # 354. The human KLK3protein can be found at Uniprot ID P07288. In some embodiments, the plurality of factors comprises Alpha-1- antichymotrypsin (SERPINA3). The human SERPINA3 gene can be found at Entrez gene # 12. The human SERPINA3 protein can be found at Uniprot ID P01011. In some embodiments, the plurality of factors comprises Lysyl oxidase homolog 2 (L0XL2). The human L0XL2 gene can be found at Entrez gene # 4017. The human LOXL2 protein can be found at Uniprot ID Q9Y4K0. In some embodiments, the plurality of factors comprises Low- density lipoprotein receptor-related protein 1, soluble:Cytoplasmic domain (LRP1). The human LRP1 gene can be found at Entrez gene # 4035. The human LRP1 protein can be found at Uniprot ID Q07954. In some embodiments, the plurality of factors comprises Low- density lipoprotein receptor-related protein 8 (LRP8). The human LRP8 gene can be found at Entrez gene # 7804. The human LRP8 protein can be found at Uniprot ID Q14114. In some embodiments, the plurality of factors comprises MAP kinase-activated protein kinase 5 (MAPKAPK5). The human MAPKAPK5 gene can be found at Entrez gene # 8550. The human MAPKAPK5 protein can be found at Uniprot ID Q8IW41. In some embodiments, the plurality of factors comprises Promotilin (MLN). The human MLN gene can be found at Entrez gene # 4295. The human MLN protein can be found at Uniprot ID P12872. In some embodiments, the plurality of factors comprises NADH dehydrogenase [ubiquinone] flavoprotein 2, mitochondrial (NDUFV2). The human NDUFV2 gene can be found at Entrez gene # 4729. The human NDUFV2 protein can be found at Uniprot ID P19404. In some embodiments, the plurality of factors comprises NF-kappa-B inhibitor delta (NFKBID). The human NFKBID gene can be found at Entrez gene # 84807. The human NFKBID protein can be found at Uniprot ID Q8NI38. In some embodiments, the plurality of factors comprises 5'-nucleotidase domain-containing protein 3 (NT5DC3). The human NT5DC3 gene can be found at Entrez gene # 51559. The human NT5DC3 protein can be found at Uniprot ID Q86UY8. In some embodiments, the plurality of factors comprises Pappalysin-1 (PAPPA). The human PAPPA gene can be found at Entrez gene # 5069. The human PAPPA protein can be found at Uniprot ID Q13219. In some embodiments, the plurality of factors comprises Alpha-parvin (PARVA). The human PARVA gene can be found at Entrez gene # 55742. The human PARVA protein can be found at Uniprot ID P43251. In some embodiments, theplurality of factors comprises Biotinidase (PLG). The human PLG gene can be found at Entrez gene # 5340. The human PLG protein can be found at Uniprot ID Q9NVD7. In some embodiments, the plurality of factors comprises Serine / threonine-protein phosphatase 4 regulatory subunit 3A (PPP4R3A). The human PPP4R3A gene can be found at Entrez gene# 55671. The human PPP4R3A protein can be found at Uniprot ID Q6IN85. In some embodiments, the plurality of factors comprises Retinol-binding protein 4 (RBP4). The human RBP4 gene can be found at Entrez gene # 5950. The human RBP4 protein can be found at Uniprot ID P02753. In some embodiments, the plurality of factors comprises Calcipressin-2 (RCAN2). The human RCAN2 gene can be found at Entrez gene # 10231. The human RCAN2 protein can be found at Uniprot ID Q14206. In some embodiments, the plurality of factors comprises Protein S100-A12 (S100A12). The human S100A12 gene can be found at Entrez gene # 6283. The human S100A12 protein can be found at Uniprot ID P80511. In some embodiments, the plurality of factors comprises Serum amyloid A-2 protein (SAA2). The human SAA2 gene can be found at Entrez gene # 6289. The human SAA2 protein can be found at Uniprot ID P0DJI9. In some embodiments, the plurality of factors comprises Kallistatin (SERPINA4). The human SERPINA4 gene can be found at Entrez gene # 5267. The human SERPINA4 protein can be found at Uniprot ID P29622. In some embodiments, the plurality of factors comprises Plasma serine protease inhibitor (SERPINA5). The human SERPINA5 gene can be found at Entrez gene # 5104. The human SERPINA5 protein can be found at Uniprot ID P05154. In some embodiments, the plurality of factors comprises Plasma protease Cl inhibitor (SERPING1). The human SERPING1 gene can be found at Entrez gene # 710. The human SERPING1 protein can be found at Uniprot ID P05155. In some embodiments, the plurality of factors comprises Protein SET (SET). The human SET gene can be found at Entrez gene # 6418. The human SET protein can be found at Uniprot ID Q01105. In some embodiments, the plurality of factors comprises Splicing factor 3B subunit 4 (SF3B4). The human SF3B4 gene can be found at Entrez gene# 10262. The human SF3B4 protein can be found at Uniprot ID Q15427. In some embodiments, the plurality of factors comprises Protachykinin- 1 (TAC1). The human TAC1 gene can be found at Entrez gene # 6863. The human TAC1 protein can be found at Uniprot ID P20366. In some embodiments, the plurality of factors comprises Serine / threonine- protein kinase TBK1 (TBK1). The human TBK1 gene can be found at Entrez gene # 29110. The human TBK1 protein can be found at Uniprot ID Q9UHD2. In some embodiments, the plurality of factors comprises Thrombospondin-2 (THBS2). The human THBS2 gene can be found at Entrez gene # 7058. The human THBS2 protein can be found at Uniprot ID P35442.In some embodiments, the plurality of factors comprises Metalloproteinase inhibitor 1 (TIMP1). The human TIMP1 gene can be found at Entrez gene # 7076. The human TIMP1 protein can be found at Uniprot ID PO1O33. In some embodiments, the plurality of factors comprises Toll-like receptor 5 (TLR5). The human TLR5 gene can be found at Entrez gene # 7100. The human TLR5 protein can be found at Uniprot ID 060602. In some embodiments, the plurality of factors comprises Transmembrane protein 132A (TMEM132A). The human TMEM132A gene can be found at Entrez gene # 54972. The human TMEM132A protein can be found at Uniprot ID Q24JP5. In some embodiments, the plurality of factors comprises Tumor necrosis factor receptor superfamily member 11B (TNFRSF11B). The human TNFRSF11B gene can be found at Entrez gene # 4982. The human TNFRSF1 IB protein can be found at Uniprot ID 000300. In some embodiments, the plurality of factors comprises Uromodulin (UMOD). The human UMOD gene can be found at Entrez gene # 7369. The human UMOD protein can be found at Uniprot ID P07911. In some embodiments, the plurality of factors comprises WAP four-disulfide core domain protein 1 (WFDC1). The human WFDC1 gene can be found at Entrez gene # 58189. The human WFDC1 protein can be found at Uniprot ID Q9HC57.

[0149] In some embodiments, the factor is selected from the factors provided in Table 4 and Table 6. In some embodiments, the factor is selected from the factors provided in Tables 4 and 6. In some embodiments, the plurality of factors is selected from the factors provided in Table 4 and Table 6. In some embodiments, the plurality of factors is selected from the factors provided in Tables 4 and 6. In some embodiments, the plurality of factors comprises at least two factors selected from those provided in Table 4 and Table 6. In some embodiments, the plurality of factors comprises at least two factors selected from those provided in Tables 4 and 6. In some embodiments, the plurality of factors consists of factors selected from those provided in Table 4 and Table 6. In some embodiments, the plurality of factors consists of factors selected from those provided in Tables 4 and 6. In some embodiments, the factor is selected from the group consisting of: KCNAB2, IL12B, IL23A, MCL1, KIR2DS2, AGA, RPN1, LAT, MFAP2, PUF60, MPZ, ACE, RNF122, TXNDC5, CDH15, FGFBP3, COL11A2, INPP5E, ADH7, MVK, RNF146, SOCS3, RBFOX2, ARFGAP1, SRSF6, RBM23, DDR1, APOF, TRA2B, MCTS1, TBCA, RGS7, PTPN9, CSNK1G2, ILF3, TPPP2, ARHGEF2, SRSF7, EWSR1, FSTL1, SPP1, FLRT2, FLRT3, VTN, ATP1B1, WFIKKN2, NRAC, PKD2, HSPA9, EMC4, ASAP2, NAP1L2, HTR7, DCUN1D3, RBL2, MAD1L1, GRB 14, RBBP5, NAB2, CSF1R, CCN4, GPD1, KLK3,CXCL13, GZMA, C9, IL12B, RAP1GAP, IGFBP1, DHX58, C0PS2, IL1RAP, CCL25, HPX, ADM, CD93, ISG15, MYL6B, HSPA1A, MBD1, TRAPPC3, AKT2, CRLF1, FTL, RBBP4, BMPER, SERPINB5, PMP2, OTC, OTOR, AOC1, FGFBP1, ATRN, NAGLU, SAA1, SAA4, CLSTN1, GSS, DLD, EPHB4, PRSS27, MUC16, CFHR2, HTRA1, KRT19, RBP4, SMOC2, BTD, TXLNA, MZB1, FADD, GSN, CDH17, LECT2, ADAMTSL1, RNASET2, SEMA4A, DDOST, BDH2, SNRPB2, G0LM1, RAB3A, CD46, SEPTIN6, WWOX, WDR5, HPCAL1, ALDH5A1, VAT1, SARS1, AFM, CD A, ITLN1, LRIG1, GREM1, PTGR2, UBE2L6, CLTA, GSR, PDCD6, SNCG, CRH, RGS21, UBE2R2, BASP1, GBP5, LMNB2, POP7, RAET1L, SEMA5B, CNTN3, UBL3, MMACHC, GTF2B, GCHFR, LRATD2, SGK1, TSEN15, SAR1B, CDK5RAP3, HAUS1, NKIRAS1, PHOSPHO2, PCDH17, TRIM5, ALDH7A1, TXNL4A, CEP20, PDE1B, ITGA4, ITGB1, LRFN3, ADGRB1, SGSH, MGAT5, B3GAT1, MGAT5, FBLN7, APBB1IP, PON2, PPP2R5D, RBFOX1, TIMP1, GEMIN7, CSNK1A1L, PHF11, BTN2A2, SKP2, SPATA46, LIN7A, BORCS5, ARRDC5, PCYT1A, PHYH, ANKRD63, VCX, NTAN1, STARD7, APOL2, FLT4, RCSD1, INIP, VMAC, XPNPEP3, IFNE, NELFA, KDM8, NCBP1, USF2, LRRC75A, APCS, PLCD1, ESPN, RFX5, RPS6KB2, N0M02, TCEAL2, CES3, DYRK1A, CYP2C19, CFI, IGFBP3, IL6, LEP, CRTC3, VEGFA, IL1RAP, HGF, PLA2G2A, CCL25, SERPINA7, POR, CCN3, HPX, IGFBP1, MMP3, FGA, FGB, FGG, BCAM, SPINT1, HAT1, GHR, CFP, CNTN1, SERPINF2, IL19, MB, C9, IGHM, LBP, NAAA, HAPLN1, IDS, NIDI, ACAN, TGFBI, DLL4, FCGR3B, ACY1, IBSP, SERPINA4, POSTN, SELE, B2M, HAMP, SERPINA1, AHSG, CKB, CKM, PROC, ANGPTL4, MBD4, PSMD7, IGHE, CXCL10, KLKB1, CFH, PFDN5, RBM39, DCTPP1, PRSS22, KYNU, IL6, SERPINA6, ITIH4, SFN, CCL7, LYZ, MMP13, STC1, CAPG, PI3, GPC5, HRG, SCGB2A1, SIRT2, TNFAIP6, CD300C, GPNMB, KRT18, TNFSF14, LEPR, PRKCG, FGL1, PGLYRP2, NPFF, MFAP4, TMX3, PRKCSH, DEFBI 12, SEMA4D, ACP6, AFP, NGF, FTH1, FTL, DMKN, EPHA10, CHRDL2, TP53, AOC1, IFNA8, CSH1, CSH2, TNC, PLTP, CCN1, CLSTN3, OIT3, GGT2, FMOD, C5orf38, VWA1, INHBC, ADGRF5, C1QL2, PCYOX1, AOC2, CFHR4, LRRC15, POSTN, UBE2J1, GFRAL, IGF2, LILRB5, LILRA6, APOA2, VWA2, DEPPI, C1QTNF3, SERPINA9, CFHR5, DLG3, GLTPD2, HBQ1, ENTPD1, AGGF1, NRG2, SPON2, FAM241B, JAML, BCHE, GPNMB, APOD, DLL1, PEAR1, RSPO4, LEP, ARL8B, PCDH10, MFAP3L, CD14, COL15A1, HAVCR1, ARHGEF10, MAN1A2, CRYZL1, TFPI2, PLXDC1, ACP2, BTD, MFAP2, ITIH2, EFCAB14, PLA1A, GZMK, YBX1, IDO1, NQO1, SPOCK3, NXT1, CYCS, SERPINA4, GOLM1, FUT5, ANGPT2, ITM2B, INHBA, CD5, THBS2, SAA2, CETP,KLK3|SERPINA3, C5, SERPINA5, MLN, TMEM132A, HDGFL3, TAC1, CXCL12, S100A12, L0XL2, RBP4, NFKBID, HHIP, UMOD, CDON, CPLX3, AD AMTS 1, KHSRP, SERPING1, ECM1, CTSB, NDUFV2, TLR5, SET, DNAJB9, APOL1, FGFR4, TNFRSF11B, KLK14, PLG, DCBLD1, BGLAP, CARNMT1, MAPKAPK5, RCAN2, CST2, INHBB, APOA4, IL17RB, BOC, LRP1, SF3B4, PARVA, GALNT9, CFHR1, WFDC1, BCHE, NT5DC3, PAPPA, KLK14, ANGPTL3, TBK1, CXCL8, LRP8, CCNH, CFB, and PPP4R3A. The amino acid sequences of these factors can be found in the Uniprot database, for example, and each factor’s Uniprot accession number is provided herein. Further, methods, reagents, and assays for measuring expression levels of these factors are well known in the art and are commercially available.

[0150] In some embodiments, the factor is selected from a factor provided in Table 12. In some embodiments, the plurality of factors is selected from the factors provided in Table 12. In some embodiments, the plurality of factors comprises at least two factors selected from those provided in Table 12. In some embodiments, the plurality of factors consists of factors selected from Table 12. In some embodiments, the factor is selected from the group consisting of: FGL1, SFN, C9, EPHA10, PRSS22, ITIH4, PUF60, RBFOX2, SPINT1, HAVCR1, OIT3, C9, TIMP1, AFM, CLSTN3, SERPINA1, EPHB4, CYCS, APOF, SERPINA4, LBP, DDR1, G0LM1, RBL2, EWSR1, SPP1, CFI, G0LM1, POSTN, IL6, PLXDC1, POSTN, RAET1L, LMNB2, AHSG, FUT5, CCN1, HAMP, SPON2, ANGPT2, ITM2B, SRSF7, SERPINB5, SERPINA4, PCYOX1, SAA1, MAN1A2, SAA4, HTRA1, MUC16, RBFOX1, ENTPD1, CFHR5, APCS, ALDH7A1, RBM23, INHBA, PTPN9, LEP, CD5, CD14, THBS2, MFAP3L, SAA2, CETP, ADM, KLK3|SERPINA3, C5, SERPINA5, MLN, ACP2, TMEM132A, HAT1, RBBP4, KDM8, PCDH17, RNASET2, HDGFL3, STC1, SERPINF2, TSEN15, FTL|FTH1, TAC1, CXCL12, C1QTNF3, YBX1, BTD, PCDH10, FMOD, LIN7A, S100A12, LOXL2, FLRT2, TFPI2, CHRDL2, RBP4, NFKBID, HHIP, NAP1L2, PROC, WDR5, DCTPP1, UMOD, BTD, KYNU, CDON, CPLX3, IGFBP1, RPN1, LEP, KLK3, TXNL4A, RBM39, PON2, CD300C, AD AMTS 1, DLD, IGHE, MFAP2, VTN, CD46, FTL, KHSRP, GHR, LRATD2, ILF3, APOA2, SERPING1, POR, ECM1, RGS7, CTSB, TRA2B, CAPG, TCEAL2, NKIRAS1, NDUFV2, CRH, TP53, TLR5, IL23A|IL12B, SET, DNAJB9, PTGR2, PI3, NAAA, NQO1, APOL1, FGFR4, TNFRSF11B, KLK14, LILRB5, PLG, TNFAIP6, DCBLD1, BGLAP, CARNMT1, LRRC15, LRRC75A, MGAT5, OTOR, CXCL13, MAPKAPK5, ARHGEF2, PDCD6, AFP, VEGFA, KCNAB2, RCAN2, CST2, HSPA9, INHBB, APOA4, MYL6B, UBE2J1, VWA2,IL17RB, BOC, LRP1, ISG15, SF3B4, PARVA, GALNT9, LRIG1, CFHR1, FLRT3, WFDC1, IDO1, BCHE, PGLYRP2, BASP1, NT5DC3, MZB1, SRSF6, PAPPA, KLK14, ANGPTL3, RPS6KB2, ATP1B1, CEP20, LEPR, TBK1, S0CS3, JAML, CXCL8, CRTC3, LRP8, CCNH, FGFBP3, VCX, CFB, PPP4R3A, SNCG, and TRIM5. The amino acid sequences of these factors can be found in the Uniprot database, for example, and each factor’s Uniprot accession number is provided herein. Further, methods, reagents, and assays for measuring expression levels of these factors are well known in the art and are commercially available.

[0151] In some embodiments, the factor is selected from a factor provided in Table 7. In some embodiments, the plurality of factors is selected from the factors provided in Table 7. In some embodiments, the plurality of factors comprises at least two factors selected from those provided in Table 7. In some embodiments, the plurality of factors consists of factors selected from Table 7. In some embodiments, the factor is selected from the group consisting of: FGL1, SFN, C9, EPHA10, PRSS22, ITIH4, PUF60, RBFOX2, SPINT1, HAVCR1, OIT3, C9, TIMP1, AFM, CLSTN3, SERPINA1, EPHB4, CYCS, APOF, SERPINA4, LBP, DDR1, G0LM1, RBL2, EWSR1, SPP1, CFI, G0LM1, POSTN, IL6, PLXDC1, POSTN, RAET1L, LMNB2, AHSG, FUT5, CCN1, HAMP, SPON2, ANGPT2, ITM2B, SRSF7, SERPINB5, SERPINA4, PCYOX1, SAA1, MAN1A2, SAA4, HTRA1, MUC16, RBFOX1, ENTPD1, CFHR5, APCS, ALDH7A1, RBM23, INHBA, PTPN9, LEP, and CD5. The amino acid sequences of these factors can be found in the Uniprot database, for example, and each factor’s Uniprot accession number is provided herein. Further, methods, reagents, and assays for measuring expression levels of these factors are well known in the art and are commercially available.

[0152] In some embodiments, the factor is selected from a factor provided in Table 8. In some embodiments, the plurality of factors is selected from the factors provided in Table 8. In some embodiments, the plurality of factors comprises at least two factors selected from those provided in Table 8. In some embodiments, the plurality of factors consists of factors selected from Table 8. In some embodiments, the factor is selected from the group consisting of: SERPINA1, SPINT1, ITIH4, LBP, PRSS22, ACP2, SPP1, ARHGEF2, SNCG,FLRT2, HAVCR1, IL6, OTOR, BTD, DLD, C9, PUF60, GHR, RBFOX2, DDR1, SFN, LRRC75A, SPON2, EPHB4, LRIG1, AFM, FGL1, ALDH7A1, EPHA10, APOF, C9, EWSR1, BTD, NAP1L2, TFPI2, POR, FMOD, RAET1L, LMNB2, LIN7A, MFAP3L, SERPINF2, CXCL13, PLXDC1, HAT1, CHRDL2, LRRC15, PON2, CFHR5, BASP1,WDR5, YBX1, PGLYRP2, KDM8, DCTPP1, TP53, RBM39, 0IT3, ILF3, TIMP1, HTRA1, ENTPD1, SRSF7, ISG15, IGFBP1, CRTC3, G0LM1, CRH, CAPG, LEP, MUC16, RPN1, KYNU, CLSTN3, KLK3, STC1, SERPINB5, FTL, MFAP2, FLRT3, CD46, FTL|FTH1, PROC, FGFBP3, CEP20, ADM, LILRB5, LRATD2, TRIM5, RGS7, RBFOX1, VTN, SRSF6, RBM23, PCDH10, CFI, ATP1B1, CD300C, TCEAL2, IDO1, RBL2, LEP, TRA2B, PTPN9, CCN1, PDCD6, LEPR, SERPINA4, PCYOX1, APCS, HSPA9, UBE2J1, RNASET2, TNFAIP6, PCDH17, MGAT5, IGHE, SOCS3, APOA2, HAMP, PI3, MZB1, TXNL4A, IL23A|IL12B, SAA4, NAAA, NKIRAS1, AFP, MYL6B, MAN1A2, VCX, AHSG, POSTN, TSEN15, PTGR2, POSTN, NQO1, SAA1, JAML, RPS6KB2, KCNAB2, CD14, VEGFA, C1QTNF3, RBBP4, and VWA2. The amino acid sequences of these factors can be found in the Uniprot database, for example, and each factor’s Uniprot accession number is provided herein. Further, methods, reagents, and assays for measuring expression levels of these factors are well known in the art and are commercially available.

[0153] In some embodiments, the factor is selected from a factor provided in Table 9. In some embodiments, the plurality of factors is selected from the factors provided in Table 9. In some embodiments, the plurality of factors comprises at least two factors selected from those provided in Table 9. In some embodiments, the plurality of factors consists of factors selected from Table 9. In some embodiments, the factor is selected from the group consisting of: SERPINA1, SPINT1, ITIH4, LBP, PRSS22, ACP2, SPP1, ARHGEF2, SNCG, FLRT2, HAVCR1, IL6, OTOR, BTD, DLD, C9, PUF60, GHR, RBFOX2, DDR1, SFN, LRRC75A, SPON2, EPHB4, LRIG1, AFM, FGL1, ALDH7A1, EPHA10, APOF, C9, EWSR1, BTD, NAP1L2, TFPI2, POR, FMOD, RAET1L, LMNB2, LIN7A, MFAP3L, SERPINF2, CXCL13, PLXDC1, HAT1, CHRDL2, LRRC15, PON2, and CFHR5. The amino acid sequences of these factors can be found in the Uniprot database, for example, and each factor’s Uniprot accession number is provided herein. Further, methods, reagents, and assays for measuring expression levels of these factors are well known in the art and are commercially available.

[0154] In some embodiments, the factor is selected from a factor provided in Table 10. In some embodiments, the plurality of factors is selected from the factors provided in Table 10. In some embodiments, the plurality of factors comprises at least two factors selected from those provided in Table 10. In some embodiments, the plurality of factors consists of factors selected from Table 10. In some embodiments, the factor is selected from the group consisting of: ITIH4, PRSS22, C9, SFN, FGL1, EPHA10, PUF60, RBFOX2, SPINT1,HAVCR1, C9, 0IT3, AFM, TIMP1, SERPINA1, CLSTN3, EPHB4, APOF, LBP, DDR1, G0LM1, RBL2, EWSR1, SPP1, CFI, POSTN, IL6, PLXDC1, POSTN, RAET1L, LMNB2, AHSG, SPON2, CCN1, HAMP, SRSF7, SERPINB5, SERPINA4, PCYOX1, MAN1A2, SAA1, HTRA1, SAA4, MUC16, RBFOX1, ENTPD1, CFHR5, ALDH7A1, APCS, RBM23, PTPN9, and LEP. The amino acid sequences of these factors can be found in the Uniprot database, for example, and each factor’s Uniprot accession number is provided herein. Further, methods, reagents, and assays for measuring expression levels of these factors are well known in the art and are commercially available.

[0155] In some embodiments, the factor is selected from a factor provided in Table 11. In some embodiments, the plurality of factors is selected from the factors provided in Table 11. In some embodiments, the plurality of factors comprises at least two factors selected from those provided in Table 11. In some embodiments, the plurality of factors consists of factors selected from Table 11. In some embodiments, the factor is selected from the group consisting of: SERPINA1, SPINT1, ITIH4, LBP, PRSS22, SPP1, HAVCR1, IL6, C9, PUF60, RBFOX2, DDR1, SFN, SPON2, EPHB4, AFM, FGL1, EPHA10, ALDH7A1, C9, APOF, EWSR1, RAET1L, LMNB2, PLXDC1, and CFHR5. The amino acid sequences of these factors can be found in the Uniprot database, for example, and each factor’s Uniprot accession number is provided herein. Further, methods, reagents, and assays for measuring expression levels of these factors are well known in the art and are commercially available.

[0156] In some embodiments, the population of responders suffers from the disease. In some embodiments, the disease is a proliferative disease. In some embodiments, the disease is cancer. In some embodiments, the responders all have the same disease. In some embodiments, the population of non-responders suffers from the disease. In some embodiments, the non-responders all suffer from the same disease. In some embodiments, the population of responders and non-responders all suffer from the same disease. In some embodiments, the population of responders and the subject suffer from the same disease. In some embodiments, the population of non-responders and the subject suffer from the same disease. In some embodiments, the population of non-responders, the population of responders and the subject suffer from the same disease.

[0157] In some embodiments, the expression levels are from the subject before receiving the therapy. In some embodiments, the expression levels are determined for the subject before receiving the therapy. In some embodiments, the expression levels are from time TO. In some embodiments, the expression levels are baseline expression levels. In some embodiments,the sample is provided by the subject before receiving the therapy. In some embodiments, the expression levels are from the subject before receiving a first treatment of the therapy. In some embodiments, the expression levels are from the subject before receiving the first cycle of the therapy. In some embodiments, a treatment is a dose. In some embodiments, a treatment is a regimen. In some embodiments, a treatment is a combination of dose and regimen.

[0158] In some embodiments, before is at least 1 hour, 2 hours, 3 hours, 6 hours, 8 hours, 12 hours, 1 day, 2 days, 3 days, 5 days, 1 week, 2 weeks, 3 weeks, 4 weeks, 1 month, 2 months, 3 months, 4 months, 5 months, or 6 months before the therapy or before administration of the therapy. Each possibility represents a separate embodiment of the invention. In some embodiments, before is at least 1 hour before. In some embodiments, before is just before the therapy or before administration of the therapy. In some embodiments, before is at most 1 hour, 2 hours, 3 hours, 4 hours, 6 hours, 9 hours, 12 hours, 18 hours, 24 hours, 2 days, 3 days, 5 days, 1 week, 2 weeks, 3 weeks, 4 weeks, 1 month, 2 months, 3 months, 4 months, 5 months, or 6 months before the therapy or before administration of the therapy. Each possibility represents a separate embodiment of the invention. In some embodiments, before is at most 24 hours before the therapy or before administration of the therapy. In some embodiments, administration of the therapy is the first administration of the therapy. In some embodiments, administration of the therapy is any administration of the therapy.

[0159] In some embodiments, the expression levels are from the subject after receiving the therapy. In some embodiments, the expression levels are from time Tl. In some embodiments, the sample is provided by the subject after receiving the therapy. In some embodiments, the expression levels are from the subject after receiving a first treatment of the therapy. In some embodiments, the expression levels are from the subject after receiving any treatment with the therapy.

[0160] In some embodiments, after is at a time after initiation of the therapy, or after administration of the therapy, sufficient for altered expression of the at least one factor. In some embodiments, after is at a time after initiation of the therapy, or after administration of the first treatment of the therapy. In some embodiments, after is at least 1 day, 2 days, 3 days, 4 days, 5 days, 6 days, 1 week, 2 weeks, 3 weeks, 4 weeks, 6 weeks, 1 month, 2 months, 3 months, 4 months, 5 months, 6 months, or a year after. Each possibility represents a separate embodiment of the invention. In some embodiments, after is at least 24 hours after. In someembodiments, after is at least 2 weeks after. In some embodiments, after is at least 3 weeks after. In some embodiments, after is at least 6 weeks after. In some embodiments, after is at most 1 week, 2 weeks, 3 weeks, 4 weeks, 6 weeks, 1 month, 2 months, 3 months, 4 months, 5 months, 6 months or a year after initiation of the therapy, or after administration of the therapy. Each possibility represents a separate embodiment of the invention.

[0161] In some embodiments, the receiving expression levels comprises receiving factor expression levels for a group of factors larger than the plurality of factors. In some embodiments, the received expression levels for the larger group are received for responders and non-responders. In some embodiments, a subgroup of proteins is selected from the group. In some embodiments, a subgroup is a subset. In some embodiments, the subgroup is designated the plurality of factors. In some embodiments, the method comprises designating. In some embodiments, the receiving further comprises for each factor of the group applying a machine learning algorithm. In some embodiments, the algorithm classifies factors as from responders and non-responders. In some embodiments, the algorithm outputs if a subject that provided the sample that had the measured factor expression level is a responder or nonresponder. In some embodiments, the receiving further comprises selecting a subgroup of factors for which the algorithm most evenly divides the subjects into responders and non- responders. In some embodiments, the subjects are all the subjects in the populations of responders and non-responders. In some embodiments, the factors are processed with an algorithm that most evenly divides all subjects, responders and non-responders, into groups of responders and non-responders (even if designations are incorrect) are selected as the subgroup. In some embodiments, the algorithm is trained on the received factor expression levels in responders and non-responders. In some embodiments, the algorithm is trained on a training set. In some embodiments, training is on expression levels and tags indicating if an expression level was from a responder or non-responder. In some embodiments, training is on expression levels, clinical information and tags indicating if an expression level was from a responder or non-responder. In some embodiments, training is on the number of the resistance associated factors. In some embodiments, training is on the number of the resistance associated factors and tags indicating if a number of resistance associate factors was from a responder or non-responder.

[0162] In some embodiments, the receiving further comprises for each factor of the group determining the average difference between responder and non-responders. In some embodiments, the receiving further comprises for each factor of the group determining thestatistical significance between the levels in responders and non-responders. In some embodiments, the statistical significance is between the averages. In some embodiments, the statistical significance is the p-value. In some embodiments, the receiving further comprises selecting a subgroup of factors with the greatest statistical significance. In some embodiments, a statistical test is applied to determine significance. In some embodiments, the test is a Kolmogorov-Smirnov test. In some embodiments, the subgroup comprises a predetermined number of factors with the greatest significance. In some embodiments, the predetermined number is about 50 factors. In some embodiments, the predetermined number is at least 50 factors.

[0163] In some embodiments, the subgroup comprises the factors whose algorithm most evenly divides the subjects. In some embodiments, evenly divides is into responders and non-responders. In some embodiments, the subgroup is the top 50, 100, 150, 200, 250, 300, 350, 400, 450, 500, 600, 700, 750, 800, 900, 1000, 1100, 1200, 1300, 1400, 1500, 1600, 1700, 1800, 1900, 2000, 3000, 4000, or 5000. Each possibility represents a separate embodiment of the invention. In some embodiments, the subgroup is the top 50. In some embodiments, the subgroup is the top 100. In some embodiments, the subgroup is the top 200. In some embodiments, the subgroup is the top 500.

[0164] In some embodiments, the method further comprises performing a dimensionality reduction step. In some embodiments, the reduction is with respect to the plurality of factors. In some embodiments, the reduction is reducing the number of factors in the plurality. In some embodiments, the dimensionality reduction step identifies a subgroup or a subset of factors. In some embodiments, factors are principal factors. In some embodiments, the training set comprises only the expression levels of the subset / subgroup of factors. In some embodiments, the subgroup or subset of factors are the factors that most evenly balance the predicted number of responders and non-responders. In some embodiments, predicted is predicted by the machine learning algorithm. In some embodiments, the machine learning algorithm is the trained machine learning algorithm. In some embodiments, the machine learning algorithm is the machine learning algorithm during training.

[0165] In some embodiments, a preprocessing stage may take place to preprocess the received expression levels. In some embodiments, the preprocessing stage may comprise at least one of data cleaning and normalizing, feature selection, feature extraction, dimensionality reduction, and / or any other suitable preprocessing method or technique.Feature selection can be performed by statistical tests, such as the Kolmogorov Smirnov (KS) test, or any other test known in the art.

[0166] In some embodiments, factor selection and / or dimensionality reduction steps may be performed, to reduce the number of factors in each sample and / or to obtain a set of principal factors, e.g., those factors that may have significant predictive power. In some embodiments, factor selection is RAP selection. Accordingly, in some embodiments, a factor selection and / or dimensionality reduction step may result in a reduction of the number of factors in each sample and / or set of values. In some embodiments, dimensionality reduction selects principal factors, e.g., proteins, based on the level of response predictive power a factor generates with respect to the desired prediction. In specific embodiments, the dimensionality reduction involves regarding all or some factors as vector components and calculating their norm.

[0167] In some embodiments, any suitable factor selection and / or dimensionality reduction method or technique may be employed, such as, but not limited to:• ANQVA with So parameter: Analysis of variance with an additional parameter (So) that controls for the relative importance of features based on resulted test p-values and difference between the group means (see, e.g., Tusher, Tibshirani and Chu, PNAS 98, pp5116-21, 2001).• Scalable EMpirical Bayes Model Selection (SEMMS): An empirical Bayes feature selection method which applies a parsimonious mixture model to identify significant predictors (see, e.g., Bar, Booth, and Wells. A scalable empirical Bayes approach to variable selection in generalized linear models, 2019).• L2N: A method for differential expression analysis that uses a three-component mixture model. The model consists of two log-normal components (L2) for differentially expressed features, one component for under-expressed features and the other for overexpressed features, and a single normal component (N) for non-differentially expressed features (see, e.g., Bar and Schifano. Differential variation and expression analysis. Stat 8, e237, doi:10.1002 / sta4.237, 2019).• Genetic algorithms: A family of heuristic optimization algorithms that employ organic evolutionary techniques such as random mutations, recombination, and natural selection as methods for achieving optimal configurations (see, e.g., Popovic, Sifrim, Pavlopoulos, Moreau, and Bart De Moor. A Simple Genetic Algorithm for Biomarker Mining. 2012).• Naive classifier: The naive classifier evaluates a response score by reducing the dimension to a single score. This is performed by regarding all features (e.g., specific profiles such as protein expression levels) as component of a vector and calculating its norm. The dimension reduction reduces the possible risk of an over-fitting. In some embodiments, the vector components are normalized according to the typical component value among patients that belong to the same response group (e.g., responders), such that the normalized norm quantifies the amount of deviation from the typical respective class value. In additional embodiments, the naive classifier enables training using data of subjects that belong only to part of the response groups.

[0168] As used herein, the term “responder” or a subject “known to respond” are used interchangeably and refer to a subject that when administered a therapy displays an improvement in at least one criteria of the disease being treated by the therapy or does not show an increase in severity of the disease. In some embodiments, a responder is a subject that when administered a therapy displays an improvement in the disease that is being treated by the therapy. In some embodiments, a responder is a subject that when administered a therapy displays a clinical benefit. In some embodiments, a responder is a subject that when administered a therapy does not show an increase in severity of the disease. In some embodiments, an increase in severity is over time. In some embodiments, does not show an increase in severity is stable disease. In some embodiments, a responder is a subject that when administered a therapy show mixed response. In some embodiments, a responder is a subject that when administered a therapy show mixed response, wherein mixed response is improvement in at least one criteria of the disease but does not show an improvement in other criteria of the disease. In some embodiments, mixed response is shrinkage of some lesions in combination with growth of new or existing lesions. In some embodiments, a responder is a subject for which the therapy produces an anti-disease response. In some embodiments, for a subject with cancer, a responder is a subject in which the therapy produces an anticancer response. In some embodiments, a response is not a reduction in side effects. In some embodiments, a response is a reduction in side effects. In some embodiments, a response is a response against the disease itself. In some embodiments, an anticancer response is an antitumor response. In some embodiments, an antitumor response comprises tumor regression. In some embodiments, an antitumor response comprises tumor shrinkage. In some embodiments, an antitumor response comprises a lack of tumor growth. In some embodiments, an antitumor response comprises a lack of tumor metastasis. In someembodiments, an antitumor response comprises a lack of tumor hyperproliferation. In some embodiments, an improvement is in at least one symptom of the disease. In some embodiments, response is complete response. In some embodiments, response is minimal response. In some embodiments, response is partial response. In some embodiments, response comprises stable disease. In some embodiments, responder is a subject with a favorable response to the therapy. In some embodiments, non-responder is a subject with a non-favorable response to the therapy. In some embodiments, a non-favorable response is an increase in tumor burden. Increases in tumor burden can encompass any increase in tumor size or total cancer cell number such as increase in tumor size, increase in tumor spread, increase in metastasis, increase in tumor cell proliferation or any other increase. In some embodiments, response is response to a monotherapy. In some embodiments, response is response to a combination therapy.

[0169] As used herein, a “favorable response” of the cancer patient indicates “responsiveness” of the cancer patient to the treatment with the therapy, namely, the treatment of the responsive cancer patient with the therapy will lead to the desired clinical outcome such as tumor regression, tumor shrinkage or tumor necrosis; reduction in tumor burden; an anti-tumor response by the immune system; preventing or delaying tumor recurrence, tumor growth or tumor metastasis. In some embodiments, the subject is complete responder or treatment with the cancer therapy leads to stable disease. In some embodiments, a complete responder is a subject in which there is an absence of detectable cancer after treatment with the therapy. In this case, it is possible and advised to continue the treatment of the responsive cancer patient with the therapy or if the patient is cancer free to discontinue treatment. In some embodiments, the method further comprises continuing to administer the therapy to a subject that is not a non-responder. In some embodiments, the subject is non- responder, a minimal responder, partial responder or has a stable disease, and the method further comprises continuing to administer the therapy to a subject, as well as treating the subject with an additional therapy (e.g., determined using the resistance associated protein (RAP) analysis provided herein) to increase responsiveness. In some embodiments, a subject that is not a non-responder is a responder.

[0170] As used herein, the term “non-responder” and a subject “known to not respond” are used interchangeably and refer to a subject that when administered a therapy displays no improvement or stabilization in disease. In some embodiments, a non-responder displays a worsening of disease when administered a therapy. In some embodiments, a non-responderis a subject that when administered a therapy displays no clinical benefit. In some embodiments, non-responder is not a subject that experiences a side effect of the therapy. In some embodiments, a non-responder is a subject in which the disease progresses. In some embodiments, a non-responder is a subject in which the disease does not stabilize after therapy. In some embodiments, a non-responder is a subject in which the disease does not improve after therapy. In some embodiments, a non-responder is a subject that is not a responder as defined hereinabove. In some embodiments, a non-responder is a subject with a non-favorable response to the therapy. In some embodiments, a non-responder is a subject resistant to the therapy. In some embodiments, a non-responder is a subject refractory to the therapy. In some embodiments, non-response is non-response to a monotherapy. In some embodiments, non-response is non-response to a combination therapy.

[0171] As used herein a “non-favorable response” of the cancer patient indicates “nonresponsiveness” of the cancer patient to the treatment with the therapy and thus the treatment of the non-responsive cancer patient with the therapy will not lead to the desired clinical outcome, and potentially to a non-desired outcomes such as tumor expansion, recurrence, or metastases. In some embodiments, the method further comprises discontinuing administration of the therapy to a subject that is a non-responder. In some embodiments the method further comprises continuing to administer the therapy to a subject, in combination with an additional therapy. In some embodiments, the additional therapy increases responsiveness of a non-responsive patient.

[0172] In some embodiments, the method is for determining whether the response is considered a durable response (e.g., a progression-free survival of more than 6 months). In some embodiments, response is response for at least 3-months. In some embodiments, the response is response at a time from treatment. In some embodiments, from treatment is from the commencement of treatment. In some embodiments, response is response at 3-months. In some embodiments, response is response for at least 6-months. In some embodiments, response is response at 6-months. In some embodiments, response is response for at least 7- months. In some embodiments, response is response at 7-months. In some embodiments, response is response for at least 1-year. In some embodiments, response is response at 1- year. In some embodiments, response is response for at least 2-year. In some embodiments, response is response at 2-year. In some embodiments, response is response for at least 3- year. In some embodiments, response is response at 3 -year. In some embodiments, response is response for at least 4-year. In some embodiments, response is response at 4-year. In someembodiments, response is response for at least 5-year. In some embodiments, response is response at 5-year. It will be understood by a skilled artisan that response for at least a given amount of time comprises at least monitoring response at that time point and also potentially monitoring response up until that time point.

[0173] In some embodiments, the method further comprises administering the therapy to the subject predicted to respond to the therapy. In some embodiments, the method further comprises continuing to administering the therapy to the subject predicted to respond to the therapy. In some embodiments, the method further comprises not administering the therapy to the subject predicted to not respond to the therapy. In some embodiments, the method further comprises discontinuing the therapy to the subject predicted to not respond to the therapy. In some embodiments, the method further comprises administering an alternative therapy to the subject predicted to be a non-responder. In some embodiments, the alternative therapy is an additional therapy. In some embodiments, the additional therapy is chemotherapy. In some embodiments, the method further comprises administering the therapy or continuing to administer the therapy in combination with an agent or therapy that blocks or inhibits at least one of the resistance-associated factors in the subject predicted to be resistant to the therapy. In some embodiments, an agent or therapy that blocks or inhibits at least one of the resistance-associated factors is an additional therapy. In some embodiments, an agent or therapy that blocks or inhibits the signaling pathway of at least one of the resistance-associated factors is an additional therapy. In some embodiments, the combination therapy is administered to a subject predicted to be a non-responder.

[0174] In some embodiments, the method further comprises administering the monotherapy to a subject predicted to respond to the monotherapy. In some embodiments, the method further comprises administering the monotherapy to a subject with PD-L1 high cancer predicted to respond to the monotherapy. In some embodiments, the method further comprises administering a combination therapy to a subject predicted to not respond to the monotherapy. In some embodiments, the method further comprises administering the combination therapy to a subject with PD-L1 high cancer predicted to not respond to the monotherapy.

[0175] In some embodiments, the method further comprises administering the combination therapy to a subject predicted to respond to the combination therapy. In some embodiments, the method further comprises administering the combination therapy to a subject with PD- L1 low or negative cancer predicted to respond to the combination therapy. In someembodiments, the method further comprises administering an alternative therapy to a subject predicted to not respond to the combination therapy. In some embodiments, the method further comprises administering an alternative therapy to a subject with PD-L1 low or negative cancer predicted to not respond to the combination therapy. Examples of alternative therapies include, but are not limited to other ICI combination (e.g., with anti-CTLA-4) and non-chemotherapeutic treatments.

[0176] In some embodiments, the method further comprises administering to the subject (e.g., a non-responder) an agent that modulates the at least one factor. In some embodiments, modulates comprises inhibits, blocks and regulates. In some embodiments, modulates is inhibits. In some embodiments, the method further comprises administering to the subject (e.g., a non-responder) an agent that modulates a pathway that comprises the at least one factor. In some embodiments, modulating the at least one factor is modulating a pathway comprising the at least one factor. In some embodiments, modulating a pathway comprising modulating a driver protein / gene that controls the at least one factor. In some embodiments, modulating a pathway comprising modulating a driver protein / gene that controls the pathway. In some embodiments, modulating a pathway comprising the at least one factor is modulating a receptor of the factor (e.g., using a receptor agonist or antagonists), a ligand or the factor, a paralog of the factor, or a combination thereof. In some embodiments, the modulating is modulating a plurality of factors. In some embodiments, the modulating is modulating a plurality of factors in the signature. In some embodiments, the modulation is modulating each factor in the signature. In some embodiments, the modulation achieves better response to therapy. In some embodiments the factor is a resistance-associated factor.

[0177] In some embodiments, a resistance score is a RAP score. In some embodiments, a resistance score is a response score. In some embodiments, a resistance score is 1 -response score. In some embodiments, a resistance score is 10-response score. In some embodiments, response score is 1-resistance score. In some embodiments, response score is 10-resistance score. It will be understood by a skilled artisan that the response score and resistance score are inverses. Thus, if the scale of the scores is 0-1 then the conversion of one score to the other is 1-score. Whereas if the scale of the scores is 0-10 then the conversion of one score to the other is 10-score. The same can be used for any scale being used for the two scores. In some embodiments, resistance score is total resistance score. In some embodiments, response score is total response score. In some embodiments, a RAP score is a total RAP score. In some embodiments, the resistance score is based on similarity of the factorexpression level in the subject to the factor expression level in the non-responders. In some embodiments, the resistance score is based on similarity of the factor expression level in the subject to the factor expression level in the responders. In some embodiments, based on is calculated based on. In some embodiments, similarity is lack of similarity. In some embodiments, similarity to responders is lack of similarity to non-responders. In some embodiments, similarity to non-responders is lack of similarity to responders. In some embodiments, similarity is measured on a scale.

[0178] In some embodiments, the scale is from 0 to 1, wherein 1 is perfectly similar to non- responders and 0 is perfectly similar to responders. In some embodiments, the resistance score is from 0 to 1 , wherein 1 is perfectly similar to non-responders and 0 is perfectly similar to responders. In some embodiments, the resistance score is based on similarity of the factor expression level in the subject to the factor expression level in the non-responders and the factor expression level in the responders. In some embodiments, the response score is from 0 to 1, wherein 1 is perfectly similar to responders and 0 is perfectly similar to non- responders. In some embodiments, the response score is the PROphet score. In some embodiments, a prophet positive subject is a subject with a response score above a predetermined threshold. In some embodiments, a prophet negative subject is a subject with a response score below a predetermined threshold. In some embodiments, the response score is based on similarity of the factor expression level in the subject to the factor expression level in the non-responders and the factor expression level in the responders. In some embodiments, a response score from 0.5 to 1 indicates the subject is a responder. In some embodiments, a response score above 0.5 indicates the subject is a responder. In some embodiments, a response score from 0.5 to 0 indicates the subject is a non-responder. In some embodiments, a response score below 0.5 indicates the subject is a non-responder.

[0179] In some embodiments, the scale is from 0 to 10, wherein 10 is perfectly similar to responders and 0 is perfectly similar to non-responders. In some embodiments, the resistance score is from 0 to 10, wherein 10 is perfectly similar to non-responders and 0 is perfectly similar to responders. In some embodiments, the resistance score is based on similarity of the factor expression level in the subject to the factor expression level in the non-responders and the factor expression level in the responders. In some embodiments, the response score is from 0 to 10, wherein 10 is perfectly similar to responders and 0 is perfectly similar to non-responders. In some embodiments, the response score is the PROphet score. In some embodiments, the response score is the total response score. In some embodiments, a prophetpositive subject is a subject with a response score above a predetermined threshold. In some embodiments, a prophet negative subject is a subject with a response score below a predetermined threshold. In some embodiments, the response score is based on similarity of the factor expression level in the subject to the factor expression level in the non-responders and the factor expression level in the responders. In some embodiments, a response score from 5 to 10 indicates the subject is a responder. In some embodiments, a response score above 5 indicates the subject is a responder. In some embodiments, a response score from 5 to 0 indicates the subject is a non-responder. In some embodiments, a response score below 5 indicates the subject is a non-responder.

[0180] In some embodiments, the method comprises before step (b) selecting a subset of factors. In some embodiments, the subset is a subset of the plurality of factors. In some embodiments, before step (b) is before the calculating. In some embodiments, the subset is a subset of the plurality of factors. In some embodiments, the subset comprises the factors that best differentiate between the responders and non-responders. In some embodiments, the factors that best differentiate are the top percentage. In some embodiments, the top percentage is the top 1, 3, 5, 10, 15, 20, 25, 30, 35, 40, 45 or 50% of factors. Each possibility represents a separate embodiment of the invention. In some embodiments, the top percentage is the top 20%. In some embodiments, the top factors are the top 10, 20, 25, 30, 40, 50, 60, 70, 75, 80, 90 or 100 factors. Each possibility represents a separate embodiment of the invention. In some embodiments, the top factors are the top 50 factors. In some embodiments, selection comprises applying a Kolmogorov-Smirnov test. In some embodiments, the Kolmogorov-Smirnov test is applied to the received factor expression levels. In some embodiments, the Kolmogorov -Smirnov test determines how well a factor differentiates between responders and non-responders. In some embodiments, the Kolmogorov-Smirnov test outputs a measure of how well a factor differentiates and the best factors are the factors with the highest scores. In some embodiments, selection comprises applying an XGBoost algorithm. In some embodiments, the calculating is for the subset. In some embodiments, the calculating is for each factor of the subset.

[0181] In some embodiments, calculating comprises applying a machine learning algorithm. In some embodiments, calculating comprises applying a machine learning model. In some embodiments, the machine learning model is a machine learning algorithm. In some embodiments, the machine learning model implements a machine learning algorithm. In some embodiments, the algorithm is a classifier. In some embodiments, the algorithm is aregression model. In some embodiments, the algorithm is supervised. In some embodiments, the algorithm is unsupervised. In some embodiments, the machine learning algorithm is trained on the expression levels in responders. In some embodiments, the machine learning algorithm is trained on the expression levels in non-responders. In some embodiments, the machine learning algorithm is trained on the expression levels in responders and non- responders. In some embodiments, the machine learning algorithm is trained on a training set. In some embodiments, the machine learning algorithm is trained by a method of the invention. In some embodiments, a machine learning algorithm is applied to factors of the plurality of factors. In some embodiments, a machine learning algorithm is applied to each factor of the plurality of factors. In some embodiments, a machine learning algorithm is applied to the subset. In some embodiments, a machine learning algorithm is applied to the subset of factors. In some embodiments, a machine learning algorithm is applied to each factor of the subset of factors. In some embodiments, each factor is analyzed and calculated separately, and the machine learning algorithm does not use expression levels of more than one factor as the training set. In some embodiments, a trained machine learning algorithm is applied to individual protein expression levels from the subject. In some embodiments, a machine learning algorithm trained on expression levels of a specific factor in responders and non-responders is applied to the expression level of that specific factor in the subject. It will be understood by a skilled artisan, that for each of the factors of the plurality of factors, a different algorithm will be trained and then applied to each expression level of the subject. Thus, if three algorithms are separately trained on expression in responders and non- responders for Factor A, Factor B and Factor C, then the algorithm trained on Factor A expression levels will be applied to the subject’s expression level of Factor A, the algorithm trained on Factor B expression levels will be applied to the subject’s expression level of Factor B, and the algorithm trained on Factor C expression levels will be applied to the subject’s expression level of Factor C. In some embodiments, during a training phase, the machine learning model is trained on a training set comprising expression data for a single factor from responders and non-responders, using corresponding annotations of “responder” or “non-responder” to predict or classify factor expression data according to classes “responder” and “non-responder”. In some embodiments, during an inference stage, the machine learning model is applied to expression data of the single factor from a subject to predict classification of the factor as similar to a responder or non-responder. In some embodiments, the classification is a resistance score. In some embodiments, theclassification is a response score. In some embodiments, the classification is a measure of how similar the factor is to non-responders and dissimilar to responders.

[0182] In some embodiments, the trained machine learning algorithm is trained to predict responsiveness of subjects suffering from the disease to the therapy. In some embodiments, the trained machine learning algorithm is trained to output a resistance score. In some embodiments, the trained machine learning algorithm is trained to output a resistance probability. In some embodiments, the trained machine learning algorithm is trained to output clinical benefit probability. In some embodiments, the trained machine learning algorithm is trained to output an activity score. In some embodiments, the trained machine learning algorithm is trained to predict activity of a resistance-associated factor in a subject. In some embodiments, the trained machine learning algorithm is trained to predict if a factor is a resistance-associated factor in the subject. In some embodiments, the trained machine learning algorithm is trained to predict if a factor of the subject is a resistance-associated factor in the subject.

[0183] In some embodiments, the trained machine learning algorithm is trained to predict responsiveness of subjects suffering from the disease to the therapy. In some embodiments, the trained machine learning algorithm is trained to output a response score. In some embodiments, the trained machine learning algorithm is trained to output a response probability. In some embodiments, the trained machine learning algorithm is trained to output clinical benefit probability. In some embodiments, the trained machine learning algorithm is trained to output an activity score. In some embodiments, the trained machine learning algorithm is trained to predict activity of a response-associated factor in a subject. In some embodiments, the trained machine learning algorithm is trained to predict if a factor is a response-associated factor in the subject. In some embodiments, the trained machine learning algorithm is trained to predict if a factor of the subject is a response-associated factor in the subject.

[0184] In some embodiments, the training set comprises received factor expression levels. In some embodiments, the training set comprises received factor expression levels in both responders and non-responders. In some embodiments, the training set comprises received factor expression levels in both mono-responders and mono-non-responders. In some embodiments, the training set comprises received factor expression levels in both comboresponders and combo-non-responders. In some embodiments, the training set comprises received factor expression levels in mono-responders, mono-non-responders, combo-responders and combo-non-responders. In some embodiments, the training set comprises received factor expression levels for only one factor. In some embodiments, the training set comprises the number of resistance-associated factors or response-associated factors expressed in samples. In some embodiments, the sample are from subjects suffering from the disease. In some embodiments, the sample are from responders. In some embodiments, the sample are from non-responders. In some embodiments, the training set comprises at least one clinical parameter. In some embodiments, the clinical parameter is from subjects. In some embodiments, subjects are responders and non-responders. In some embodiments, the training set comprises labels. In some embodiments, the labels are associated with the responsiveness of the subjects. In some embodiments, the labels are responder or nonresponder. In some embodiments, the resistance-associated factors are labeled with the labels. In some embodiments, the expression levels of the resistance-associated factors are labeled with the labels. In some embodiments, the at least one clinical parameter is labeled with the label.

[0185] According to some embodiments, the training set further comprises at least one clinical parameter of each responder and non-responder and the machine learning algorithm is applied to individual received factor expression levels from the subject and the subject’s at least one clinical parameter. In some embodiments, the at least one clinical parameter is the sex of the subjects. In some embodiments, the training set further comprises the sex of the subjects. In some embodiments, the subjects are each subject. In some embodiments, sex is gender. In some embodiments, the at least one clinical parameter is sex. In some embodiments, sex is a subject’s sex. In some embodiments, sex is male or female. In some embodiments, sex is sex at birth. In some embodiments, the training set comprises the sex of each responder. In some embodiments, the training set comprise the sex of each non- responder. In some embodiments, the training set comprises the sex of each mono-responder. In some embodiments, the training set comprise the sex of each mono-non-responder. In some embodiments, the training set comprises the sex of each combo-responder. In some embodiments, the training set comprise the sex of each combo-non-responder. In some embodiments, the clinical parameter is age. In some embodiments, age is a subject’s age. In some embodiments, the clinical parameter is the line of treatment. In some embodiments, the line of treatment parameter is whether the therapy was a first line of treatment or an advanced treatment. In some embodiments, a line of treatment is first line treatment. In some embodiments, a line of treatment is a secondary treatment. In some embodiments, secondarytreatment is an advanced treatment. It will be understood by a skilled artisan that advanced treatment may be any line of treatment after the first, e.g., second line, third line, fourth line, fifth line, etc. In some embodiments, the clinical parameter is whether the treatment is a first line treatment or an advanced treatment. In some embodiments, the clinical parameter is PD- L1 status. In some embodiments, PD-L1 status is PD-L1 status of the cancer. Methods of measuring PD-L1 levels in cancer cells (e.g., a tumor) are well known in the art and any such method may be employed. In some embodiments, PD-L1 status comprises high PD-L1 or low PD-L1. In some embodiments, PD-L1 status comprises high PD-L1, low PD-L1 or no PD-L1. In some embodiments, PD-L1 status comprises high PD-L1, medium PD-L1 or low PD-L1. In some embodiments, PD-L1 levels are numeric values between 0 to 100. In some embodiments, PD-L1 levels are percentages between 0 to 100. In some embodiments, PD- L1 status comprises PD-L1 expression in less than 1% of cancer cells, in 1-49% of cancer cells, or in 50% or more of cancer cells. In some embodiments, PD-L1 expression in less than 1% of cancer cells is no PD-L1 expression. In some embodiments, PD-L1 low or negative cancer comprises fewer than 50% of cancer cells being positive for PD-L1 expression. In some embodiments, expression is surface expression. In some embodiments, PD-L1 negative cancer comprises fewer than 1% of cancer cells being positive for PD-L1 expression. In some embodiments, PD-L1 expression in less than 1% of cancer cells is low PD-L1 expression. In some embodiments, PD-L1 expression in 1-49% of cancer cells is low PD-L1 expression. In some embodiments, PD-L1 low cancer comprises fewer than 1-49% of cancer cells being positive for PD-L1 expression. In some embodiments, PD-L1 expression in 1-49% of cancer cells is medium PD-L1 expression. In some embodiments, PD-L1 expression in 50% or more of cancer cells is high PD-L1 expression. In some embodiments, a high PD-L1 cancer comprises expression in at least 50% of cells. In some embodiments, PD-L1 high cancer comprises at least 50% of cancer cells being positive for PD-L1 expression. In some embodiments, a low PD-L1 cancer comprises expression in 1- 49% of cells. In some embodiments, a no PD-L1 cancer comprises expression in 0% of cells. In some embodiments, a no PD-L1 cancer comprises expression in less than 1% of cells. In some embodiments, the PD-L1 low or negative cancer is PD-L1 low cancer. In some embodiments, the PD-L1 low or negative cancer is PD-L1 negative cancer. In some embodiments, a no PD-L1 cancer is a PD-L1 negative cancer.

[0186] In some embodiments, the clinical parameter is a known biomarker of the disease or mutations in known biomarkers of the disease. In some embodiments, the biomarker isselected from MYC, NOTCH, EGFR, HER2, BRAF, KRAS, MAP2K1, MET, NRAS, NTRK1, NTRK2, NTRK3, PIK3CA, RET, ROS1, TP53, ALK, CDKN2A, KIT, NF1, BFAST, FGFR, LDH, PTEN, RB I, PD-L1, MSI (Micro satelite Instability), TMB (Tumor Mutational Burden), or a combination thereof. In some embodiments, the clinical parameter is expression of the biomarker. In some embodiments, expression is percent expression. In some embodiments, expression is mutational status.

[0187] In some embodiments, the training set further comprises the sex, age and PD-L1 status of each responder and non-responder. In some embodiments, the training set further comprises the sex of each responder and non-responder. In some embodiments, the training set further comprises the age and PD-L1 status of each responder and non-responder. In some embodiments, the machine learning algorithm is applied to individual received factor expression levels from the subject and the subject’s sex. In some embodiments, the machine learning algorithm is applied to individual received factor expression levels from the subject and the subject’s sex, age andPD-Ll status. In some embodiments, the calculating comprises applying a machine learning algorithm trained on a training set comprising the received factor expression levels in responders and non-responders and at least one clinical parameter, to the expression levels from the subject and the subject’s at least one clinical parameter and wherein the machine learning algorithm outputs the resistance score. In some embodiments, the training comprises the received factor expression levels in responders and non- responders and clinical parameters of each responder and non-responder and the machine learning algorithm is applied to individual received factor expression levels from the subject and the subject’s clinical parameters and wherein the machine learning algorithm outputs response score. In some embodiments, the training comprises the received factor expression levels in responders and non-responders and a clinical parameter selected from sex, age and PD-L1 expression, or any combination thereof, of each responder and non-responder and the machine learning algorithm is applied to individual received factor expression levels from the subject and the subject’s clinical parameters and wherein the machine learning algorithm outputs response prediction. In some embodiments, the training set comprises the number of resistance associated factors in each responder and non-responder and at least one clinical parameter and the machine learning algorithm is applied to the number of resistance associated factors from the subject and the subject’s at least one clinical parameters and wherein the machine learning algorithm outputs a response prediction. In some embodiments, the training set comprises the number of resistance associated factors in eachresponder and non-responder and sex of each responder and non-responder and the machine learning algorithm is applied to the number of resistance associated factors from the subject and the subject’s sex and wherein the machine learning algorithm outputs a response prediction. In some embodiments, the training set comprises the number of resistance associated factors in each responder and non-responder, age and PD-L1 status of each responder and non-responder and the machine learning algorithm is applied to the number of resistance associated factors from the subject and the subject’s age and PD-L1 status and wherein the machine learning algorithm outputs a response prediction.

[0188] In some embodiments, the training set comprises the received factor expression levels in responder and non-responders. In some embodiments, the training set comprises the received factor expression levels in responder and non-responders and a clinical parameter. In some embodiments, the training set comprises the received factor expression levels in responder and non-responders and sex of each of the responders and non- responders. In some embodiments, the trained machine learning algorithm is applied to individual received factor expression levels from the subject. In some embodiments, the trained machine learning algorithm is applied to each received factor expression levels from the subject. In some embodiments, the trained machine learning algorithm is applied to individual received factor expression levels from the subject and a clinical parameter from the subject. In some embodiments, the trained machine learning algorithm is applied to individual received factor expression levels from the subject and the subject’s sex.

[0189] In some embodiments, the clinical parameter is the type of treatment. In some embodiments, the clinical parameter is expression of a target of the therapy. In some embodiments, the clinical parameter is expression of a protein within a process that is a target of the therapy. In some embodiments, the process is a process comprising the target of the therapy. In some embodiments, expression is expression in the subject. In some embodiments, expression is expression in a diseased tissue. In some embodiments, expression is expression in a diseased tissue sample. In some embodiments, expression is expression in the tumor. In some embodiments, expression is expression in a tumor sample. In some embodiments, a tumor sample is a biopsy. In some embodiments, expression is expression not in the tumor. In some embodiments, expression is expression not in a tumor sample. In some embodiments, expression is expression in a liquid biopsy. In some embodiments, expression is percent expression. In some embodiments, percent is percent of cells. In some embodiments, the therapy is anti-PD-1 therapy and the protein in the processis PD-L1. In some embodiments, the therapy is anti-PD-Ll therapy, and the target protein is PD-L1. In some embodiments, the clinical parameter is PD-L1 expression. In some embodiments the training set comprises at least one clinical parameter selected from line of treatment, PD-L1 expression, sex and age. In some embodiments the training set comprises protein expression levels and sex. In some embodiments the training set comprises number of RAPs, age and PD-L1 status.

[0190] Additionally clinical parameters may also be included. A skilled artisan will be able to select relevant clinical parameters for inclusion in the training set. Examples of additional clinical parameters include, but are not limited to, histological type of the sample (e.g., adenocarcinoma, squamous cell carcinoma, etc.), metastatic location, tumor location, cancer staging (such as tumor, nodes and metastases, TNM, staging for example), performance status (such as ECOG performance status), genetic mutations, epigenetic status, general medical history, vital signs, blood measurements, renal and liver function, weight, height, pulse, blood pressure and smoking history.

[0191] In some embodiments, at an inference stage the trained machine learning algorithm is applied. In some embodiments, the trained machine learning algorithm is applied to individual received factor expression levels. In some embodiments, the trained machine learning algorithm is applied to a plurality of received factor expression levels. In some embodiments, the trained machine learning algorithm is applied to individual received factor expression levels and the at least one clinical parameter. In some embodiments, the trained machine learning algorithm is applied to individual received factor expression levels from the subjects and the subject’s sex. In some embodiments, the trained machine learning algorithm is applied to the number of resistance-associated proteins. In some embodiments, the trained machine learning algorithm is applied to the number of resistance-associated factors. In some embodiments, the trained machine learning algorithm is applied to the number of resistance-associated factors and at least one clinical parameter.

[0192] In some embodiments, at the inference stage an input is received. In some embodiments, the input comprises the number of resistance-associated factors expressed in a sample. In some embodiments, the sample is from a subject. In some embodiments, the input comprises at least one clinical parameter. In some embodiments, the subject suffers from the disease. In some embodiments, the subject has unknown responsiveness to the therapy. In some embodiments, the parameter is of the subject with unknown responsiveness. In some embodiments, at the inference stage the trained machine learning algorithm isapplied. In some embodiments, applied is applied to the input. In some embodiments, the input is the received input. In some embodiments, the inference stage is to predict responsiveness. In some embodiments, responsiveness is responsiveness to the therapy of the subject with unknown responsiveness.

[0193] In some embodiments, the machine learning algorithm outputs the resistance score. In some embodiments, the outputted resistance score is scaled from 0 to 1. In some embodiments, 1 is perfectly similar to non-responders and 0 is perfectly similar to responders. In some embodiments, the machine learning algorithm calculates similarity to responders. In some embodiments, the machine learning algorithm calculates similarity to non-responders. In some embodiments, the machine learning algorithm outputs a numeric value of similarity to responders and non-responders. In some embodiments, a protein is considered to be a RAP if its resistance score is beyond a certain threshold. In some embodiments, the threshold for the resistance score is calculated on a scale of 0 to 1. In some embodiments, the threshold for the resistance score of a certain protein is between 0.2 and 0.95. In some embodiments, the threshold for the resistance score of a certain protein is about 0.2, 0.25, 0.3, 0.35, 0.4, 0.42, 0.45, 0.5, 0.55, 0.6, 0.65, 0.7, 0.75, 0.8, 0.85, 0.9, or 0.95. Each possibility represents a separate embodiment of the invention. In some embodiments, the threshold for the resistance score is 0.25. In some embodiments, the threshold for the resistance score is 0.42. In some embodiments, the threshold for the resistance score is 0.6. In some embodiments, the threshold for the resistance score when calculated by a machine learning algorithm is about 0.2, 0.25, 0.3, 0.35, 0.4, 0.42, 0.45, 0.5, 0.55, 0.6, 0.65, 0.7, 0.75, 0.8, 0.85, 0.9, or 0.95. Each possibility represents a separate embodiment of the invention. In some embodiments, the threshold for the resistance score when calculated with a machine learning algorithm is 0.25. In some embodiments, the threshold for the resistance score when calculated with a machine learning algorithm is 0.42. In some embodiments, the threshold for the resistance score when calculated with a machine learning algorithm is 0.6. In some embodiments, the threshold is determined according to the percent of responders in the training set. It will be understood that if, for example, 25% of the training set were responders then the threshold would be set at 0.25.

[0194] In some embodiments, response probability is determined by the calculation (1- resistance score). In some embodiments, 1 -resistance score is 1 -total resistance score. In some embodiments, the resistance score is the total resistance score. In some embodiments, response probability is a response score. In some embodiments, the machine learningalgorithm outputs the response score. In some embodiments, the outputted response score is scaled from 0 to 1. In some embodiments, 1 is perfectly similar to responders and 0 is perfectly similar to non-responders. In some embodiments, the machine learning algorithm calculates similarity to responders. In some embodiments, the machine learning algorithm calculates similarity to non-responders. In some embodiments, the machine learning algorithm outputs a numeric value of similarity to responders and non-responders. In some embodiments, a protein is considered to be a RAP if its response score is beyond a certain threshold. In some embodiments, a protein is considered to be an active RAP if its response score is beyond a certain threshold. In some embodiments, the threshold for the response score is calculated on a scale of 0 to 1. In some embodiments, the threshold for the response score of a certain protein is between 0.2 and 0.95. In some embodiments, the threshold for the response score of a certain protein is about 0.2, 0.25, 0.3, 0.35, 0.4, 0.42, 0.45, 0.5, 0.55, 0.6, 0.65, 0.7, 0.75, 0.8, 0.85, 0.9, or 0.95. Each possibility represents a separate embodiment of the invention. In some embodiments, the threshold for the response score is 0.25. In some embodiments, the threshold for the response score is 0.276. In some embodiments, the threshold for the response score is 0.42. In some embodiments, the threshold for the response score is 0.5. In some embodiments, the threshold for the response score is 0.6. In some embodiments, the threshold for the response score when calculated by a machine learning algorithm is about 0.2, 0.25, 0.3, 0.35, 0.4, 0.42, 0.45, 0.5, 0.55, 0.6, 0.65, 0.7, 0.75, 0.8, 0.85, 0.9, or 0.95. Each possibility represents a separate embodiment of the invention. In some embodiments, the threshold for the response score when calculated with a machine learning algorithm is 0.25. In some embodiments, the threshold for the response score is 0.276. In some embodiments, the threshold for the response score when calculated with a machine learning algorithm is 0.42. In some embodiments, the threshold for the response score when calculated with a machine learning algorithm is 0.5. In some embodiments, the threshold for the response score when calculated with a machine learning algorithm is 0.6. In some embodiments, the algorithm outputs response probability, and the response probability is calculated on a scale of 0 to 1. In some embodiments, the algorithm outputs response probability, and the response probability is calculated on a scale of 0 to 10. In some embodiments, the algorithm outputs response probability, and the response probability is calculated on a scale of 0% to 100%, wherein 100% is a perfect responder and 0% is perfect non-responder. In some embodiments, a response probability above 50% indicates a subject likely to respond. In some embodiments, a response probability below 50% indicates a subject unlikely to respond. In some embodiments, the threshold for the response score whencalculated with a machine learning algorithm is 0.25. In some embodiments, a protein with a response score above 0.25 is active in the subject. In some embodiments, the threshold for the response score when calculated with a machine learning algorithm is 0.5. In some embodiments, a protein with a response score above 0.5 is active in the subject. In some embodiments, the algorithm outputs clinical benefit probability. In some embodiments, the clinical benefit probability is calculated on a scale of 0 to 1. In some embodiments, a clinical benefit probability of 0 indicates a 0% likelihood of clinical benefit to the subject. In some embodiments, a clinical benefit probability of 1 indicates a 100% likelihood of clinical benefit to the subject. In some embodiments, the algorithm outputs clinical benefit probability, and the clinical benefit probability is calculated on a scale of 0 to 10. In some embodiments, a clinical benefit probability of 10 indicates a 100% likelihood of clinical benefit to the subject. In some embodiments, the algorithm outputs clinical benefit probability, and the clinical benefit probability is calculated on a scale of 0% to 100%. In some embodiments, a clinical benefit probability of 100% indicates a 100% likelihood of clinical benefit to the subject. In some embodiments, a clinical benefit probability of 0% indicates a 0% likelihood of clinical benefit to the subject. In some embodiments, greater than 50% likelihood of clinical benefit to the subject indicates the subject should continue or be administered the therapy. In some embodiments, the therapy is a monotherapy. In some embodiments, the therapy is a combination therapy. In some embodiments, the threshold for the clinical benefit probability is the median clinical benefit probability in the development set. In some embodiments, the threshold for the clinical benefit probability is the median clinical benefit probability in the development set, wherein a clinical benefit probability higher than the median clinical benefit probability is responder and a clinical benefit probability lower than the median clinical benefit probability is non-responder. According to some other embodiments, response probability or clinical benefit probability beyond 50% indicates the subject is responsive to therapy. According to some other embodiments, response probability or clinical benefit probability below 50% indicates the subject is non- responsive to therapy. In some embodiments, the response probability or the clinical benefit probability is from 0-10, and response probability or clinical benefit probability beyond 5 indicates the subject is responsive to therapy. In some embodiments, the response probability or the clinical benefit probability is from 0-10, and response probability or clinical benefit probability below 5 indicates the subject is non-responsive to therapy.

[0195] In some embodiments, the score is between zero and 1. In some embodiments, active is active in the cancer. In some embodiments, active is active in the subject. In some embodiments, active is active in promoting resistance. In some embodiments, beyond a threshold is below a threshold. In some embodiments, beyond a threshold is above a threshold. In some embodiments, the predetermined threshold is 0.5, 0.4, 0.3, 0.25, 0.2, 0.15, 0.1, 0.05, 0.01, 0.005, 0.001, 0.0005 or 0.0001. Each possibility represents a separate embodiment of the invention. In some embodiments, the threshold is 0.05. In some embodiments, the threshold is 5%. In some embodiments, the number of active RAPs is combined to give a total number of RAPs active in the subject. In some embodiments, the number of active RAPs is linearized to provide a total score between 0 and 1. In some embodiments, linearized is linearly scaled. In some embodiments, linearizing comprises a linear regression. In some embodiments, the number of active RAPs is converted to a total score between 0 and 1. In some embodiments, the number of active RAPs is linearized to provide a total score between 0 and 10. In some embodiments, linearized is linearly scaled. In some embodiments, linearizing comprises a linear regression. In some embodiments, the number of active RAPs is converted to a total score between 0 and 10.

[0196] In some embodiments, the predetermined threshold is determined by performing a cross-validation within the training set. In some embodiments, the predetermined threshold is the median score in the training set. In some embodiments, the predetermined threshold is the score that best distinguishes between responders and non-responders in the training set. In some embodiments, the training set is development set.

[0197] In some embodiments, the machine learning algorithm outputs the resistance score. In some embodiments, the resistance score is the RAP score. In some embodiments, the outputted resistance score is scaled from 0 to 1. In some embodiments, 1 is perfectly similar to non-responders and 0 is perfectly similar to responders. In some embodiments, for a response score 1 is perfectly similar to responders and 0 is perfectly similar to non- responders. In some embodiments, the machine learning algorithm calculates similarity to responders. In some embodiments, the machine learning algorithm calculates similarity to non-responders. In some embodiments, the machine learning algorithm outputs a numeric value of similarity to responders and non-responders. In some embodiments, a protein is considered to be a RAP if its resistance score is beyond a certain threshold. In some embodiments, the threshold for the resistance score is calculated on a scale of 0 to 1. In some embodiments, the threshold for the resistance score of a certain protein is between 0.2 and0.95. In some embodiments, the threshold for the resistance score of a certain protein is about 0.01, 0.05, 0.1, 0.15, 0.2, 0.25, 0.3, 0.35, 0.4, 0.42, 0.45, 0.5, 0.55, 0.6, 0.65, 0.7, 0.75, 0.8, 0.85, 0.9, or 0.95. Each possibility represents a separate embodiment of the invention. In some embodiments, the threshold for the resistance score is 0.25. In some embodiments, the threshold for the resistance score is 0.42. In some embodiments, the threshold for the resistance score is 0.6. In some embodiments, the threshold for the resistance score when calculated by a machine learning algorithm is about 0.2, 0.25, 0.3, 0.35, 0.4, 0.42, 0.45, 0.5, 0.55, 0.6, 0.65, 0.7, 0.75, 0.8, 0.85, 0.9, or 0.95. Each possibility represents a separate embodiment of the invention. In some embodiments, the threshold for the resistance score when calculated with a machine learning algorithm is 0.25. In some embodiments, the threshold for the resistance score when calculated with a machine learning algorithm is 0.42. In some embodiments, the threshold for the resistance score when calculated with a machine learning algorithm is 0.6.

[0198] In some embodiments, response probability is determined by the calculation (1- resistance score). In some embodiments, 1 -resistance score is 1 -total resistance score. In some embodiments, the resistance score is the total resistance score. In some embodiments, response probability is a response score. In some embodiments, the machine learning algorithm outputs the response score. In some embodiments, the outputted response score is scaled from 0 to 1. In some embodiments, 1 is perfectly similar to responders and 0 is perfectly similar to non-responders. In some embodiments, the machine learning algorithm calculates similarity to responders. In some embodiments, the machine learning algorithm calculates similarity to non-responders. In some embodiments, the machine learning algorithm outputs a numeric value of similarity to responders and non-responders. In some embodiments, a protein is considered to be a RAP if its response score is beyond a certain threshold. In some embodiments, beyond is above. In some embodiments, beyond is below. In some embodiments, the threshold for the response score is calculated on a scale of 0 to 1. In some embodiments, the threshold for the response score of a certain protein is between 0.2 and 0.95. In some embodiments, the threshold for the response score of a certain protein is about 0.2, 0.25, 0.3, 0.35, 0.4, 0.42, 0.45, 0.5, 0.55, 0.6, 0.65, 0.7, 0.75, 0.8, 0.85, 0.9, or 0.95. Each possibility represents a separate embodiment of the invention. In some embodiments, the threshold for the response score is 0.25. In some embodiments, the threshold for the response score is 0.42. In some embodiments, the threshold for the response score is 0.6. In some embodiments, the threshold for the response score when calculated bya machine learning algorithm is about 0.2, 0.25, 0.3, 0.35, 0.4, 0.42, 0.45, 0.5, 0.55, 0.6, 0.65, 0.7, 0.75, 0.8, 0.85, 0.9, or 0.95. Each possibility represents a separate embodiment of the invention. In some embodiments, the threshold for the response score when calculated with a machine learning algorithm is 0.25. In some embodiments, the threshold for the response score when calculated with a machine learning algorithm is 0.42. In some embodiments, the threshold for the response score when calculated with a machine learning algorithm is 0.6.

[0199] In some embodiments, the calculated resistance scores are combined to produce a total resistance score. In some embodiments, the calculated response scores are combined to produce a total response score. It will be understood by a skilled artisan that as the response and resistance scores are just 1 minus the other, they are always interchangeable. The conversion of resistance to response can be performed on the individual factor level or after the scores are combined and performed on the total level. In some embodiments, combine is sum. In some embodiments, the resistance scores are summed to produce a total resistance score. In some embodiments, combine is average. In some embodiments, the resistance scores are averaged to produce a total resistance score. In some embodiments, the scores are weighted when combined.

[0200] In some embodiments, the method comprises determining the number of factors of the plurality of factors that are active in the subject. In some embodiments, an active factor is a factor with a resistance score above a predetermined threshold. In some embodiments, the threshold is 0.25. In some embodiments, a factor with a resistance score above 0.25 is a factor active in the subject. In some embodiments, the threshold is 0.276. In some embodiments, a factor with a resistance score above 0.276 is a factor active in the subject. In some embodiments, only the active factors are combined. In some embodiments, the combining the calculated resistance scores is combining the active resistance scores. In some embodiments, combining comprises adding up the number of factors that are active in the subject. In some embodiments, the number of factors active in the subject is converted into a score from 0 to 1. In some embodiments, the number of factors active in the subject is converted into a score from 0 to 10. In some embodiments, converted comprises applying a linear regression model. In some embodiments, the number of active factors is linearized to provide a total score between 0 and 1. In some embodiments, the number of active factors is linearized to provide a total score between 0 and 10. In some embodiments, linearized is linearly scaled. In some embodiments, linearizing comprises a linear regression. In some embodiments, the threshold is 5.

[0201] In some embodiments, the machine learning model is a machine learning algorithm. In some embodiments, the algorithm is a supervised learning algorithm. In some embodiments, the algorithm is an unsupervised learning algorithm. In some embodiments, the algorithm is a reinforcement learning algorithm. In some embodiments, the machine learning model is a Convolutional Neural Network (CNN). In some embodiments, the at least one hardware processor trains a machine learning model. In some embodiments, the model is based, at least in part, on a training set. In some embodiments, the model is based on a training set. In some embodiments, the model is trained on a training set. In some embodiments, the at least one hardware processor applies the machine learning model to a factor expression level from a subject.

[0202] In some embodiments, the calculating comprises calculating a mean expression for each protein in responders. In some embodiments, the calculating comprises calculating a mean expression for each protein in non-responders. In some embodiments, the calculating comprises calculating a mean expression for each protein in responders and a mean expression for each protein in non-responders. In some embodiments, the calculating comprises calculating a distribution of the expression for each protein in responders and non- responders. In some embodiments, the calculating comprises calculating a standard deviation of expression for each protein in responders and non-responders. In some embodiments, in responders is in the responders population. In some embodiments, in non- responders is in the non-responders population. In some embodiments, the resistance score is based on the ratio of deviation of the factor expression in the subject from the calculated mean in responders to the deviation of the factor expression in the subject from the calculated mean in non-responders. Calculation of deviation is well known to one skilled in the art. It will be understood that the more dissimilar the expression in the subject is from a mean the larger the deviation will be. Thus, factors that are very dissimilar to the mean in responders will have a large numerator in the calculation of this ratio and factors that are lowly dissimilar to the mean in non-responders will have a small denominator. Thus, the more dissimilar to responder expression and the more similar to non-responder expression is expression of a factor in a subject the higher the resistance score will be. In some embodiments, a resistance score beyond a predetermined threshold indicates a factor is a resistance-associated factor. In some embodiments, a resistance-associated factor is a resistance-associated protein (RAP). In some embodiments, resistance-associated factor is a RAP if its expression in responders is statistically different from its expression in non-responders.

[0203] In some embodiments, the calculating further comprises calculating a distribution for each factor in responders. In some embodiments, the calculating further comprises calculating a distribution for each factor in non-responders. In some embodiments, the calculating further comprises calculating a distribution for each factor in responders and a distribution for each factor in non-responders. In some embodiments, the calculating further comprises calculating a standard deviation for each factor in responders. In some embodiments, the calculating further comprises calculating a standard deviation for each factor in non-responders. In some embodiments, the calculating further comprises calculating a standard deviation for each factor in responders and a standard deviation for each protein in non-responders. In some embodiments, the calculating further comprises calculating a standard deviation for each factor in a mix of responders and non-responders. In some embodiments, the deviation is measured as a multiple of the calculated standard deviation. It will be understood by a skilled artisan that by scaling the deviation to the standard deviation for a group of expression values the deviation can be given in more absolute terms allow for the comparison of factors and populations with very small and very large stand deviations (which may also have very low and very high expression levels).

[0204] In some embodiments, the resistance score is based on a Z-score for the expression level of each factor in the subject. In some embodiments, the resistance score is based on the Z-score relative to responders. In some embodiments, the resistance score is based on the Z- score relative to non-responders. In some embodiments, the resistance score is based on both the Z-score relative to responders and the Z-score relative to non-responders. In some embodiments, the resistance score is based on the ratio of the Z-score relative to responders to the Z-score relative to non-responders. It will be well known to a skilled artisan that a Z- score counts the distance of the individual level from the population mean in units of the population standard deviation. In some embodiments, the Z-score is calculated by Equation 1.some embodiments, ZR is the deviation of the factor expression in the subject from the calculated mean in responders. In some embodiments, ZNR is the deviation of the factor expression in the subject from the calculated mean in non-responders. In some embodiments, | | is the Z-score of the deviation. In some embodiments, | | is the standardizing of the deviation to a multiple of the standard deviation. In some embodiments, c is a constant. In some embodiments, constant is a regulation constant that prevents the score from divergencefor ZNR= 0. In some embodiments, the resistance score is calculated by Equation 2. In some embodiments, monotonoic is an ad-hoc function that prevents the resistance score from decreasing for extreme values within the non-responder distributions. In some embodiments, function is the function provided in Algorithm 1.

[0206] In some embodiments, a resistance score beyond a predetermined threshold indicates a factor is a RAP. In some embodiments, beyond is above. In some embodiments, the threshold is a predetermined threshold. In some embodiments, threshold is a threshold value. In some embodiments, the threshold for the resistance score is about 1.0, 1.1, 1.2, 1.3, 1.4,1.5, 1.6, 1.7, 1.8, 1.9, 2.0, 2.1, 2.2, 2.3, 2.4, 2.5, 2.6, 2.7, 2.8, 2.9, 3.0, 3.1, 3.2, 3.3, 3.4, 3.5,3.6, 3.7, 3.8, 3.9, 4.0, 4.1, 4.2, 4.3, 4.4, 4.5, 4.6, 4.7, 4.8, 5.0, 5.1, 5.2, 5.3, 5.4, 5.5, 5.6, 5.7, 5.8, 5.9, 6.0. 6.1, 6.2, 6.3, 6.4, 6.5, 6.6, 6.7, 6.8, 6.9, or 7.0. Each possibility represents a separate embodiment of the invention. In some embodiments, the threshold is about 0.05, 0.1, 0.15, 0.2, 0.25, 0.3, 0.35, 0.4, 0.45, 0.5, 0.55, 0.6, 0.67, 0.7, 0.75, 0.8, 0.85 or 0.9. Each possibility represents a separate embodiment of the invention. In some embodiments, the threshold for the resistance score is about 2.9. In some embodiments, the threshold for the resistance score is 2.9. In some embodiments, the threshold for the resistance score is about 3.0. In some embodiments, the threshold for the resistance score is 3.0. In some embodiments, the threshold for the resistance score is calculated on a scale of arbitrary units. In some embodiments, the threshold for the resistance score when calculated by a mathematical calculation is about 1.0, 1.1, 1.2, 1.3, 1.4, 1.5, 1.6, 1.7, 1.8, 1.9, 2.0, 2.1, 2.2,2.3, 2.4, 2.5, 2.6, 2.7, 2.8, 2.9, 3.0, 3.1, 3.2, 3.3, 3.4, 3.5, 3.6, 3.7, 3.8, 3.9, 4.0, 4.1, 4.2, 4.3,4.4, 4.5, 4.6, 4.7, 4.8, or 5.0. Each possibility represents a separate embodiment of the invention. In some embodiments, the threshold for the resistance score when calculated with a mathematical calculation is about 2.9. In some embodiments, the threshold for the resistance score when calculated with a mathematical calculation is 2.9. In some embodiments, the threshold for the resistance score when calculated with a mathematical calculation is about 3.0. In some embodiments, the threshold for the resistance score when calculated with a mathematical calculation is 3.0. In some embodiments, a mathematical calculation is a method that comprises calculating a mean expression for each protein.

[0207] In some embodiments, a subject with a number of resistance-associated factors (e.g., RAPs) above a predetermined number is predicted to be resistant to the therapy. In some embodiments, a subject with a number of resistance-associated factors above a predetermined number is predicted to not respond to the therapy. In some embodiments, asubject with a number of resistance-associated factors above a predetermined number is predicted to be a non-responder to the therapy. In some embodiments, a subject with a number of resistance-associated factors below a predetermined number is predicted to be suitable to the therapy. In some embodiments, a subject with a number of resistance- associated factors below a predetermined number is predicted to respond to the therapy. In some embodiments, a subject with a number of resistance-associated factors below a predetermined number is predicted to be a responder to the therapy. In some embodiments, a subject with a number of resistance-associated factors at or below a predetermined number is predicted to be suitable to the therapy. In some embodiments, a subject with a number of resistance-associated factors at or below a predetermined number is predicted to respond to the therapy. In some embodiments, a subject with a number of resistance-associated factors at or below a predetermined number is predicted to be a responder to the therapy.

[0208] In some embodiments, the predetermined number is a threshold number. In some embodiments, the predetermined number is 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19 or 20. Each possibility represents a separate embodiment of the invention. In some embodiments, the predetermined number is 3. In some embodiments, the predetermined number is 4. In some embodiments, the predetermined number is 7. In some embodiments, the predetermined number is 13.

[0209] In some embodiments, the number of resistance-associated proteins is determined. In some embodiments, the number of resistance-associated proteins is the total number of resistance-associated factors. In some embodiments, the number of resistance-associated proteins / factors is converted into a total resistance score. In some embodiments, the converting is with a machine learning algorithm. In some embodiments, the converting comprises applying a machine learning algorithm to the total number of resistance- associated protein. In some embodiments, the converting comprises transformation. In some embodiments, the transformation is by linear regression. In some embodiments, the total number is linearized onto a scale of total resistance scores.

[0210] In some embodiments, the method further comprises classifications of the resistance- associated factors into at least one pathway, process, or network. In some embodiments, the method further comprises performing analysis on resistance associated factors to determine at least one pathway, process, or network in which the resistance-associated factors are involved. In some embodiments, the pathway, process, or network causes nonresponsiveness to the therapy. In some embodiments, the analysis is selected from pathwayanalysis, process analysis and network analysis. In some embodiments, the method further comprises performing pathway analysis on RAPs. In some embodiments, the method further comprises performing process analysis on RAPs. In some embodiments, the method further comprises performing network analysis on RAPs. In some embodiments, at least one pathway, process or network comprises at least 2, 3, 4, 5, 6, 7, 8, 9 or 10 pathways, processes, or networks. Each possibility represents a separate embodiment of the invention. In some embodiments, at least one pathway, process or network is all the pathways, processes or networks known to include the resistance associated factors. In some embodiments, at least one pathway, process or network is all the pathways, processes or networks enriched with resistance associated factors. In some embodiments, enriched is the most enriched. In some embodiments, enriched comprises contains the most RAPs of any or the pathways, processes or networks.

[0211] In some embodiments, the method comprises selecting a pathway, process or network. In some embodiments, the selected pathway, process or network is hypothesized to affect non-response to the therapy. In some embodiments, the selected pathway, process or network is hypothesized to cause non-response to the therapy. In some embodiments, the selected pathway, process or network is known to be druggable. In some embodiments, known to be druggable comprises a known therapeutic agent that modulates the pathway, process or network. In some embodiments, the known therapeutic agent is in or has concluded clinical trials. In some embodiments, the known therapeutic agent is approved for human use. In some embodiments, approved for human use is approved for use in treating the disease in a human. In some embodiments, the disease is cancer. In some embodiments, the method further comprises administering to a subject that is a non-responder, or predicted to be a non-responder, an agent that modulates the at least one pathway, process, or network containing a resistance associated factor. In some embodiments, the agent inhibits a target in said pathway, process, or network. In some embodiments, the target is a gene. In some embodiments, the target is a protein. In some embodiments, the protein is a regulatory RNA. In some embodiments, the target is a response associated factor. In some embodiments, the target is not a response associated factor. In some embodiments, the agent activates a target in the pathway, process, or network. In some embodiments, the agent modulates the pathway, process or network. In some embodiments, the pathway’s activity induces nonresponse, and the agent inhibits the pathway. In some embodiments, the pathway’s activity reduces non-response, and the agent activates the pathway. It will be understood by a skilledI l lartisan that a response associated factor is identified by its expression in a subject being more similar to the expression in non-responders than responders. Thus, for example, if the factor is more highly expressed in non-responders and increases activity of the pathway / process / network then the agent would inhibit the pathway. If, for example, the factor is more highly expressed in non-responders, but decreases activity of the pathway / process / network then the agent would activate the pathway / process / network. Similarly, if the factor, for example, is more lowly expressed in non-responders and decreases activity of the pathway / process / network the agent would inhibit the pathway / process / network. And lastly, if, for example, the factor is more lowly expressed in non-responders but increases activity of the pathway / process / network the agent would activate the pathway / process / network. Essentially, the agent should induce the pathway / process / network to function more as it does in responders. In some embodiments, the agent targets a hub target in the pathway. In some embodiments, the agent targets a regulator target in the pathway. In some embodiments, the process activity induces nonresponse, and the agent inhibits the process. In some embodiments, the processes’ activity reduces non-response, and the agent activates the process. In some embodiments, the agent targets a hub target in the process. In some embodiments, the agent targets a regulator target in the process. In some embodiments, the network activity induces non-response, and the agent inhibits the network. In some embodiments, the network activity reduces nonresponse, and the agent activates the network. In some embodiments, the agent targets a hub factor in the network. In some embodiments, the agent targets a regulator factor in the network. In some embodiments, the regulator is a master regulator. The factors can be classified into pathways, protein interaction or signals using any analysis tool known in the art. Examples include, but are not limited to, GO analysis, Ingenuity analysis, Metacore analysis (Clarivate Analytics), reactome pathway analysis and functional analysis.

[0212] In some embodiments, clusters of RAPs are determined. In some embodiments, the plurality comprises at least one RAP from each cluster. In some embodiments, the first cluster comprises or consists of: SERPINA1, ITIH4, LBP, HSPA1A, GHR, CDH15, ITIH2, SPON2, CFI, AFM, ALDH7A1, BCHE, ATP1B1, APOF, C9, BTD, VWA1, MGAT5, RBL2, LEP, PLXDC1, VCX, RBP4, SERPINA4, CFHR5, PCYOX1, PGLYRP2, IL6, TP53, AHSG, PLA2G2A, HTRA1, B2M, MGAT5, SAA1, CNTN3, LYZ, CD14, SIRT2, LEP, C1QTNF3, and TXNL4A. In some embodiments, the second cluster comprises or consists of: SNRPB2, MMACHC, APBB1IP, PUF60, RBFOX1, PLTP, RBFOX2, NXT1,SFN, SRSF6, RBM23, NELFA, RFX5, EPHA10, EWSR1, LMNB2, TRA2B, WDR5,YBX1, RCSD1, DCTPP1, RBM39, OIT3, ILF3, SRSF7, CAPG, RBBP4, CLSTN3,PFDN5, and SERPINB5. In some embodiments, the third cluster comprises or consists of:C0L15A1, LRRC15, CD46, SM0C2, FLRT2, FMOD, CD93, ADAMTSL1, ITLN1, andEPHB4. In some embodiments, the fourth cluster comprises or consists of: TMX3, HSPA9, DDOST, RPN1 and TXNDC5. In some embodiments, the fifth cluster comprises or consists of: ACY1, OTC and CYP2C19. In some embodiments, the sixth cluster comprises or consists of: FTL, HAMP and FTH1. In some embodiments, the plurality comprises at least one factor from each of cluster 1, cluster 2, cluster 3, cluster 4, cluster 5 and cluster 6. In some embodiments, the plurality comprises at least one factor from cluster 1. In some embodiments, the plurality comprises at least one factor from cluster 2. In some embodiments, the plurality comprises at least one factor from cluster 3. In some embodiments, the plurality comprises at least one factor from cluster 4. In some embodiments, the plurality comprises at least one factor from cluster 5. In some embodiments, the plurality comprises at least one factor from cluster 6. By another aspect there is provided, a computer program product comprising a non-transitory computer- readable storage medium having program code embodied thereon, the program code executable by at least one hardware processor to perform a method of the invention.

[0213] The present invention may be a system, a method, and / or a computer program product. The computer program product may include a computer readable storage medium (or media) having computer readable program instructions thereon for causing a processor to carry out aspects of the present invention.

[0214] The computer readable storage medium can be a tangible device that can retain and store instructions for use by an instruction execution device. The computer readable storage medium may be, for example, but is not limited to, an electronic storage device, a magnetic storage device, an optical storage device, an electromagnetic storage device, a semiconductor storage device, or any suitable combination of the foregoing. A non- exhaustive list of more specific examples of the computer readable storage medium includes the following: a portable computer diskette, a hard disk, a random access memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or Flash memory), a static random access memory (SRAM), a portable compact disc read-only memory (CD-ROM), a digital versa...

Claims

CLAIMS:

1. A method of predicting response of a subject suffering from cancer to an anticancer therapy, the method comprising: a. receiving factor expression levels for a plurality of factors i. in a population of subjects suffering from cancer and known to respond to said anticancer therapy (responders); ii. in a population of subjects suffering from cancer and known to not respond to said anticancer therapy (non-responders); and iii. in said subject; b. calculate for factors of said plurality of factors a resistance score, wherein said calculating comprises applying a machine learning algorithm trained on a training set comprising said received factor expression levels in responders and non-responders to individual received factor expression levels from said subject and wherein said machine learning algorithm outputs said resistance score; and c. combine said calculated resistance scores to produce a total resistance score or determining the number of factors with a resistance score above a predetermined threshold to produce a total number of resistance-associated factors in the subject and convert said total number to a total resistance score; wherein a subject with a total resistance score beyond a predetermined threshold is predicted to not respond to said anticancer therapy and a subject with a total resistance score within said predetermined threshold is predicted to respond to said anticancer therapy; wherein said plurality of factors is selected from the factors provided in Tables 4 and 6; thereby predicting response of a subject to a monotherapy.

2. The method of claim 1, wherein said total resistance score is converted to a total response score and wherein a total response score above a predetermined threshold indicates the subject is responsive to said anticancer therapy and a total response score below a predetermined threshold indicates the subject is not responsive to said anticancer therapy.

3. The method of claim 1 or 2, wherein said anticancer therapy is a monotherapy comprising chemotherapy, targeted therapy or immunotherapy.

4. The method of claim 3, wherein said cancer is a PD-L1 high cancer and said monotherapy is an anti-PD-l / PD-Ll immunotherapy.

5. The method of claim 3, wherein said monotherapy is chemotherapy.

6. The method of claim 4 or 5, wherein said training set further comprises received factor expression levels in subjects suffering from cancer and known to respond to a combination therapy comprising an anti-PD-l / PD-Ll immunotherapy and chemotherapy (combo-responders) and received factor expression levels in subject suffering from cancer and known to not respond to said combination therapy (combo- non-responders).

7. The method of claim 1 or 2, wherein said cancer is a PD-L1 low or negative cancer and said anticancer therapy is a combination therapy comprising an anti-PD-l / PD-Ll immunotherapy and chemotherapy.

8. The method of any one of claims 1 to 7, wherein said plurality of factors comprises at least two factors selected from the factors provided in Table 4.

9. The method of claim 8, wherein said plurality of factors consists of factors selected from Table 4.

10. The method of any one of claims 1 to 9, wherein said plurality of factors comprises at least two factors selected from the factors provided in Table 5.

11. The method of claim 10, wherein said plurality of factors consists of factors selected from Table 5.

12. The method of any one of claims 1 to 9, wherein said plurality of factors comprises at least two factors selected from the factors provided in Table 8.

13. The method of claim 12, wherein said plurality of factors consists of factors selected from Table 8.

14. The method of any one of claims 1 to 13, wherein said plurality of factors comprises at least two factors selected from the factors provided in Table 9.

15. The method of claim 14, wherein said plurality of factors consists of factors selected from Table 9.

16. The method of any one of claims 1 to 9 and 12 to 15, wherein said plurality of factors comprises at least two factors selected from the factors provided in Table 10.

17. The method of claim 16, wherein said plurality of factors consists of factors selected from Table 10.

18. The method of any one of claims 1 to 9 and 12 to 17, wherein said plurality of factors comprises at least two factors selected from the factors provided in Table 11.

19. The method of claim 18, wherein said plurality of factors consists of factors selected from Table 11.

20. The method of any one of claims 1 to 7, wherein said plurality of factors comprises at least two factors selected from the factors provided in Table 12.

21. The method of claim 20, wherein said plurality of factors consists of factors selected from Table 12.

22. The method of any one of claims 1 to 7 and 20 to 21, wherein said plurality of factors comprises at least two factors selected from the factors provided in Table 6.

23. The method of claim 22, wherein said plurality of factors consists of factors selected from Table 6.

24. The method of any one of claims 1 to 7 and 20 to 21, wherein said plurality of factors comprises at least two factors selected from the factors provided in Table 7.

25. The method of claim 24, wherein said plurality of factors consists of factors selected from Table 7.

26. The method of any one of claims 1 to 25, wherein said responders and non-responders are determined based on progression free survival (PFS) at 1 year after initiation of said monotherapy or combination therapy.

27. The method of any one of claims 1 to 26, comprising before (b) selecting a subset of said plurality of factors, wherein said subset comprises factors that best differentiate between said responders and non-responders, and wherein said calculating is for each factor of said subset.

28. The method of claim 27, wherein said selecting comprises applying a statistical test to said received factor expression levels, optionally wherein said statistical test is a Kolmogorov-Smirnov test.

29. The method of claim 27 or 28, wherein said subset consists of at least 50 factors.

30. The method of any one of claims 1 to 29 wherein said factor expression level is from a time point before administration of the anticancer therapy to said subject.

31. The method of any one of claims 1 to 30, wherein said combining is averaging.

32. The method of any one of claims 1 to 30, wherein said combining comprises determining the total number of factors with a resistance score above a predetermined threshold and producing a total resistance score proportional to said total number.

33. The method of any one of claims 1 to 32, wherein said converting comprises transformation by linear regression.

34. The method of any one of claims 1 to 33, wherein said cancer is selected from hepatobiliary cancer, cervical cancer, urogenital cancer, anogenital cancer, prostate cancer, thyroid cancer, ovarian cancer, nervous system cancer, ocular cancer, lung cancer, soft tissue cancer, bone cancer, pancreatic cancer, bladder cancer, skin cancer, intestinal cancer, hepatic cancer, rectal cancer, colorectal cancer, esophageal cancer, gastric cancer, gastroesophageal cancer, breast cancer, renal cancer, skin cancer, head and neck cancer, leukemia and lymphoma.

35. The method of claim 34, wherein said cancer is selected from lung cancer, skin cancer, anogenital cancer, cervical cancer, renal cancer and head and neck cancer.

36. The method of claim 34 or 35, wherein said cancer is non-small cell lung cancer (NSCLC).

37. The method of any one of claims 1 to 36, wherein said cancer is a tyrosine kinase inhibitor resistant cancer.

38. The method of any one of claims 1 to 37, wherein said predetermined threshold is determined by performing a cross-validation within said training set or is the median score of said training set.

39. The method of any one of claims 1 to 38, wherein said plurality of factors is at least 200 factors.

40. The method of any one of claims 1 to 39, wherein said factors expression levels are factors expression levels in a biological sample provided by said subjects.

41. The method of claim 40, wherein said biological sample is selected from blood plasma, whole blood, blood serum or peripheral blood mononuclear cells.

42. The method of claim 41, wherein said biological sample is blood plasma or blood serum.

43. The method of any one of claims 1 to 42, further comprising administering said anticancer therapy to said subject predicted to respond to said anticancer therapy or administering an alternative therapy to said subject predicted to not respond to said anticancer therapy.

44. The method of any one of claims 3 to 6 and 8 to 42, further comprising administering said monotherapy to said subject predicted to respond to said monotherapy or administering a combined therapy comprising said anti-PD-l / PD-Ll immunotherapy and chemotherapy to said subject predicted to not respond to said monotherapy.

45. The method of claim 7 to 44, further comprising administering said combination therapy to said subject predicted to respond to said combination therapy or administering an alternative therapy to said subject predicted to not respond to said combination therapy.

46. The method of any one of claims 3 to 45, wherein said anti-PD-l / PD-Ll immunotherapy is selected from Pembrolizumab, Nivolumab, Durvalumab and Atezolizumab.

47. The method of any one of claims 5 to 46, wherein said chemotherapy is selected from Carboplatin, Paclitaxel, Nab-Paclitaxel, Pemetrexed, Vinorelbine, and Cisplatin.

48. The method of claim 47, wherein said combination therapy is selected from: a. Carboplatin, Durvalumab, and Paclitaxel; b. Atezolizumab, Bevacizumab, Carboplatin, and Paclitaxel; c. Carboplatin, Nab-Paclitaxel, and Pembrolizumab; d. Carboplatin, Nivolumab, and Paclitaxel; e. Carboplatin, Nivolumab, Pemetrexed; f. Carboplatin, Paclitaxel, Pembrolizumab; g. Carboplatin, Paclitaxel, Pembrolizumab, and radiation; h. Carboplatin, and Pembrolizumab; i. Carboplatin, Pembrolizumab, and Pemetrexed; j. Carboplatin, Pembrolizumab, and Vinorelbine; and k. Cisplatin, Pembrolizumab, and Pemetrexed.

49. The method of any one of claims 1 to 48, wherein predicting response comprises predicting overall survival.

50. The method of any one of claims 1 to 49, wherein predicting response comprises predicting progression free survival.

51. The method of claim 50, wherein progression free survival is at 1 year after initiation of said monotherapy or combination therapy.

52. The method of any one of claims 7 to 51, wherein the subject suffers from a negative PD-L1 cancer.

53. The method of any one of claims 1 to 52, wherein PD-L1 high cancer comprises at least 50% of cancer cells being positive for surface expression of PD-L1 and PD-L1low or negative cancer comprises fewer than 50% of cancer cells being positive for surface expression of PD-L1.

54. The method of any one of claims 7 to 53, wherein said PD-L1 low or negative cancer is PD-L1 negative cancer comprising less than 1% of cells being positive for surface expression of PD-L1.

55. The method of any one of claims 1 to 54, wherein said trained machine learning algorithm is trained by a method comprising: at a training stage, training a machine learning algorithm on a training set comprising:(i) factor expression levels of resistance-associated factors in samples from subjects suffering from cancer and known to be responsive to said anticancer therapy and factor expression levels of resistance-associated factors in samples from subjects suffering from said cancer and known to be non- responsive to said anticancer therapy; and(ii) labels associated with the responsiveness of said subjects suffering from said cancer; to produce a trained machine learning algorithm, wherein said trained machine learning algorithm is trained to output said resistance score and wherein said resistance-associated factors are selected from those provided in Tables 4 and 6.

56. The method of claim 55, wherein said expression levels of resistance-associated factors are labeled with said labels.

57. The method of claim 55 or 56, wherein said total resistance score predetermined threshold is 5 and a resistance score above 5 indicates the subject is resistant to the therapy or said total resistance score is converted to a total response score by the equation (10-total resistance score) and wherein a total response score above a predetermined threshold indicates the subject is responsive to therapy, optionally wherein said total response score predetermined threshold is 5.

Images