Algorithmically Reconstructed Molecular Pathways as the New Generation of Prognostic Molecular Biomarkers in Human Solid Cancers

Individual gene expression and molecular pathway activation profiles were shown to be effective biomarkers in many cancers. Here, we used the human interactome model to algorithmically build 7470 molecular pathways centered around individual gene products. We assessed their associations with tumor type and survival in comparison with the previous generation of molecular pathway biomarkers (3022 "classical" pathways) and with the RNA transcripts or proteomic profiles of individual genes, for 8141 and 1117 samples, respectively. For all analytes in RNA and proteomic data, respectively, we found a total of 7441 and 7343 potential biomarker associations for gene-centric pathways, 3020 and 2950 for classical pathways, and 24,349 and 6742 for individual genes. Overall, the percentage of RNA biomarkers was statistically significantly higher for both types of pathways than for individual genes (p < 0.05). In turn, both types of pathways showed comparable performance. The percentage of cancer-type-specific biomarkers was comparable between proteomic and transcriptomic levels, but the proportion of survival biomarkers was dramatically lower for proteomic data. Thus, we conclude that pathway activation level is the advanced type of biomarker for RNA and proteomic data, and momentary algorithmic computer building of pathways is a new credible alternative to time-consuming hypothesis-driven manual pathway curation and reconstruction.

Human Blood Serum Can Diminish EGFR-Targeted Inhibition of Squamous Carcinoma Cell Growth through Reactivation of MAPK and EGFR Pathways

Regardless of the presence or absence of specific diagnostic mutations, many cancer patients fail to respond to EGFR-targeted therapeutics, and a personalized approach is needed to identify putative (non)responders. We found previously that human peripheral blood and EGF can modulate the activities of EGFR-specific drugs on inhibiting clonogenity in model EGFR-positive A431 squamous carcinoma cells. Here, we report that human serum can dramatically abolish the cell growth rate inhibition by EGFR-specific drugs cetuximab and erlotinib. We show that this phenomenon is linked with derepression of drug-induced G1S cell cycle transition arrest. Furthermore, A431 cell growth inhibition by cetuximab, erlotinib, and EGF correlates with a decreased activity of ERK1/2 proteins. In turn, the EGF- and human serum-mediated rescue of drug-treated A431 cells restores ERK1/2 activity in functional tests. RNA sequencing revealed 1271 and 1566 differentially expressed genes (DEGs) in the presence of cetuximab and erlotinib, respectively. Erlotinib- and cetuximab-specific DEGs significantly overlapped. Interestingly, the expression of 100% and 75% of these DEGs restores to the no-drug level when EGF or a mixed human serum sample, respectively, is added along with cetuximab. In the case of erlotinib, EGF and human serum restore the expression of 39% and 83% of DEGs, respectively. We further assessed differential molecular pathway activation levels and propose that EGF/human serum-mediated A431 resistance to EGFR drugs can be largely explained by reactivation of the MAPK signaling cascade.

Large-scale assessment of pros and cons of autopsy-derived or tumor-matched tissues as the norms for gene expression analysis in cancers

Normal tissues are essential for studying disease-specific differential gene expression. However, healthy human controls are typically available only in postmortal/autopsy settings. In cancer research, fragments of pathologically normal tissue adjacent to tumor site are frequently used as the controls. However, it is largely underexplored how cancers can systematically influence gene expression of the neighboring tissues. Here we performed a comprehensive pan-cancer comparison of molecular profiles of solid tumor-adjacent and autopsy-derived "healthy" normal tissues. We found a number of systemic molecular differences related to activation of the immune cells, intracellular transport and autophagy, cellular respiration, telomerase activation, p38 signaling, cytoskeleton remodeling, and reorganization of the extracellular matrix. The tumor-adjacent tissues were deficient in apoptotic signaling and negative regulation of cell growth including G2/M cell cycle transition checkpoint. We also detected an extensive rearrangement of the chemical perception network. Molecular targets of 32 and 37 cancer drugs were over- or underexpressed, respectively, in the tumor-adjacent norms. These processes may be driven by molecular events that are correlated between the paired cancer and adjacent normal tissues, that mostly relate to inflammation and regulation of intracellular molecular pathways such as the p38, MAPK, Notch, and IGF1 signaling. However, using a model of macaque postmortal tissues we showed that for the 30 min - 24-hour time frame at 4ºC, an RNA degradation pattern in lung biosamples resulted in an artifact "differential" expression profile for 1140 genes, although no differences could be detected in liver. Thus, such concerns should be addressed in practice.

Gene Expression-Based Signature Can Predict Sorafenib Response in Kidney Cancer

Sorafenib is a tyrosine kinase inhibitory drug with multiple molecular specificities that is approved for clinical use in second-line treatments of metastatic and advanced renal cell carcinomas (RCCs). However, only 10–40% of RCC patients respond on sorafenib-containing therapies, and personalization of its prescription may help in finding an adequate balance of clinical efficiency, cost-effectiveness, and side effects. We investigated whether expression levels of known molecular targets of sorafenib in RCC can serve as prognostic biomarker of treatment response. We used Illumina microarrays to profile RNA expression in pre-treatment formalin-fixed paraffin-embedded (FFPE) samples of 22 metastatic or advanced RCC cases with known responses on next-line sorafenib monotherapy. Among them, nine patients showed partial response (PR), three patients—stable disease (SD), and 10 patients—progressive disease (PD) according to Response Evaluation Criteria In Solid Tumors (RECIST) criteria. We then classified PR + SD patients as “responders” and PD patients as “poor responders”. We found that gene signature including eight sorafenib target genes was congruent with the drug response characteristics and enabled high-quality separation of the responders and poor responders [area under a receiver operating characteristic curve (AUC) 0.89]. We validated these findings on another set of 13 experimental annotated FFPE RCC samples (for 2 PR, 1 SD, and 10 PD patients) that were profiled by RNA sequencing and observed AUC 0.97 for 8-gene signature as the response classifier. We further validated these results in a series of qRT-PCR experiments on the third experimental set of 12 annotated RCC biosamples (for 4 PR, 3 SD, and 5 PD patients), where 8-gene signature showed AUC 0.83.

Machine Learning Applicability for Classification of PAD/VCD Chemotherapy Response Using 53 Multiple Myeloma RNA Sequencing Profiles

Multiple myeloma (MM) affects ~500,000 people and results in ~100,000 deaths annually, being currently considered treatable but incurable. There are several MM chemotherapy treatment regimens, among which eleven include bortezomib, a proteasome-targeted drug. MM patients respond differently to bortezomib, and new prognostic biomarkers are needed to personalize treatments. However, there is a shortage of clinically annotated MM molecular data that could be used to establish novel molecular diagnostics. We report new RNA sequencing profiles for 53 MM patients annotated with responses on two similar chemotherapy regimens: bortezomib, doxorubicin, dexamethasone (PAD), and bortezomib, cyclophosphamide, dexamethasone (VCD), or with responses to their combinations. Fourteen patients received both PAD and VCD; six received only PAD, and 33 received only VCD. We compared profiles for the good and poor responders and found five genes commonly regulated here and in the previous datasets for other bortezomib regimens (all upregulated in the good responders): FGFR3, MAF, IGHA2, IGHV1-69, and GRB14. Four of these genes are linked with known immunoglobulin locus rearrangements. We then used five machine learning (ML) methods to build a classifier distinguishing good and poor responders for two cohorts: PAD + VCD (53 patients), and separately VCD (47 patients). We showed that the application of FloWPS dynamic data trimming was beneficial for all ML methods tested in both cohorts, and also in the previous MM bortezomib datasets. However, the ML models build for the different datasets did not allow cross-transferring, which can be due to different treatment regimens, experimental profiling methods, and MM heterogeneity.

DNA repair pathway activation features in follicular and papillary thyroid tumors, interrogated using 95 experimental RNA sequencing profiles

DNA repair can prevent mutations and cancer development, but it can also restore damaged tumor cells after chemo and radiation therapy. We performed RNA sequencing on 95 human pathological thyroid biosamples including 17 follicular adenomas, 23 follicular cancers, 3 medullar cancers, 51 papillary cancers and 1 poorly differentiated cancer. The gene expression profiles are annotated here with the clinical and histological diagnoses and, for papillary cancers, with BRAF gene V600E mutation status. DNA repair molecular pathway analysis showed strongly upregulated pathway activation levels for most of the differential pathways in the papillary cancer and moderately upregulated pattern in the follicular cancer, when compared to the follicular adenomas. This was observed for the BRCA1, ATM, p53, excision repair, and mismatch repair pathways. This finding was validated using independent thyroid tumor expression dataset PRJEB11591. We also analyzed gene expression patterns linked with the radioiodine resistant thyroid tumors (n = 13) and identified 871 differential genes that according to Gene Ontology analysis formed two functional groups: (i) response to topologically incorrect protein and (ii) aldo-keto reductase (NADP) activity. We also found RNA sequencing reads for two hybrid transcripts: one in-frame fusion for well-known NCOA4-RET translocation, and another frameshift fusion of ALK oncogene with a new partner ARHGAP12. The latter could probably support increased expression of truncated ALK downstream from 4th exon out of 28. Both fusions were found in papillary thyroid cancers of follicular histologic subtype with node metastases, one of them (NCOA4-RET) for the radioactive iodine resistant tumor. The differences in DNA repair activation patterns may help to improve therapy of different thyroid cancer types under investigation and the data communicated may serve for finding additional markers of radioiodine resistance.

Clinical use of RNA sequencing and oncobox analytics to predict personalized targeted therapeutic efficacy

e13676

Background: Analysis of mutation profiles in cancer patients does not provide clinical benefits in 80-90% of cases in the US (Marquart et al., 2018). Gene expression analysis potentially complements standard detection of clinically relevant mutations. Methods: 239 adult late-stage cancer patients. RNA gene expression sequencing completed on solid tumor samples using FFPE blocks. Patient mRNA profiles were analyzed using Oncobox bioinformatics, prioritizing target drugs according to their personalized predicted efficacy. Summary reports were provided to oncologists and resulting treatment selection and outcomes were assessed. Results: As of February 2020, feedback was received from participating doctors for 224 patients; 34 patients died before therapy prescription, 52 patients received treatment other than targeted therapy (chemo, surgery, radiation, or palliative care), 75 patients received at least one targeted therapy (single or combination therapy) predicted to be effective based on Oncobox analysis (“RNAseq cohort”). 63 patients received chemo or other drug therapy predicted to be potentially ineffective from Oncobox analysis (“other cohort”). Therapeutic response was obtained on 46 patients with biopsies collected no longer than 6 months prior to analysis who had no further surgery (30 in the RNAseq cohort and 16 in the other cohort). 63% of the RNAseq cohort obtained either partial response or stable disease using Oncobox guided therapies, compared to 44% of the other cohort (19% increase of disease control). The RNAseq cohort had higher mean prior therapies (1.3) compared to the other cohort (0.8) indicating more advanced disease. The similarly designed WINTHER trial reported ~8% increase of disease control using gene expression-guided vs mutation-guided therapeutics in a cohort of advanced cancer patients averaging three prior therapies (Rodon et al., 2019). Conclusions: Collectively these data suggest that gene expression profiling provides a more clinically relevant therapeutic match, and better response rates, than mutation guided therapeutic treatments. This potentially results in improved clinical outcomes for cancer patients. Clinical trial information: NCT03724097.

RNA Sequencing in Comparison to Immunohistochemistry for Measuring Cancer Biomarkers in Breast Cancer and Lung Cancer Specimens

RNA sequencing is considered the gold standard for high-throughput profiling of gene expression at the transcriptional level. Its increasing importance in cancer research and molecular diagnostics is reflected in the growing number of its mentions in scientific literature and clinical trial reports. However, the use of different reagents and protocols for RNA sequencing often produces incompatible results. Recently, we published the Oncobox Atlas of RNA sequencing profiles for normal human tissues obtained from healthy donors killed in road accidents. This is a database of molecular profiles obtained using uniform protocol and reagents settings that can be broadly used in biomedicine for data normalization in pathology, including cancer. Here, we publish new original 39 breast cancer (BC) and 19 lung cancer (LC) RNA sequencing profiles obtained for formalin-fixed paraffin-embedded (FFPE) tissue samples, fully compatible with the Oncobox Atlas. We performed the first correlation study of RNA sequencing and immunohistochemistry-measured expression profiles for the clinically actionable biomarker genes in FFPE cancer tissue samples. We demonstrated high (Spearman's rho 0.65-0.798) and statistically significant (p < 0.00004) correlations between the RNA sequencing (Oncobox protocol) and immunohistochemical measurements for HER2/ERBB2, ER/ESR1 and PGR genes in BC, and for PDL1 gene in LC; AUC: 0.963 for HER2, 0.921 for ESR1, 0.912 for PGR, and 0.922 for PDL1. To our knowledge, this is the first validation that total RNA sequencing of archived FFPE materials provides a reliable estimation of marker protein levels. These results show that in the future, RNA sequencing can complement immunohistochemistry for reliable measurements of the expression biomarkers in FFPE cancer samples.

RNA sequencing for research and diagnostics in clinical oncology

Molecular diagnostics is becoming one of the major drivers of personalized oncology. With hundreds of different approved anticancer drugs and regimens of their administration, selecting the proper treatment for a patient is at least nontrivial task. This is especially sound for the cases of recurrent and metastatic cancers where the standard lines of therapy failed. Recent trials demonstrated that mutation assays have a strong limitation in personalized selection of therapeutics, consequently, most of the drugs cannot be ranked and only a small percentage of patients can benefit from the screening. Other approaches are, therefore, needed to address a problem of finding proper targeted therapies. The analysis of RNA expression (transcriptomic) profiles presents a reasonable solution because transcriptomics stands a few steps closer to tumor phenotype than the genome analysis. Several recent studies pioneered using transcriptomics for practical oncology and showed truly encouraging clinical results. The possibility of directly measuring of expression levels of molecular drugs' targets and profiling activation of the relevant molecular pathways enables personalized prioritizing for all types of molecular-targeted therapies. RNA sequencing is the most robust tool for the high throughput quantitative transcriptomics. Its use, potentials, and limitations for the clinical oncology will be reviewed here along with the technical aspects such as optimal types of biosamples, RNA sequencing profile normalization, quality controls and several levels of data analysis.

Oncobox, gene expression-based second opinion system for predicting response to treatment in advanced solid tumors

e13143

Background: Anticancer Targeted Drugs (ATDs) specifically target one or a few types of tumor-related molecules in a cell. More than two hundred of ATDs were approved worldwide. They have different mechanisms of action and are effective for different cohorts of patients. However, many individual cases remain poorly responsive and it is of great importance to identify predictive markers of ATD efficacy. Our aim was to develop a platform enabling smart selection of the most efficient ATD therapies. Methods: We generated a second-opinion platform for clinical oncologists termed Oncobox. It is based on the analysis of gene expression profile of a cancer sample in comparison with the corresponding normal tissue biosamples in order to personalize selection of targeted drugs for individual cancer patients. Based on RNA-seq gene expression data, pathway activation levels are calculated and along with the concentrations of molecular target genes products used as predictors of tumor response to ATDs. Results: The Oncobox method was tested first on the retrospective samples of advanced tumors: gastric, renal cell, ovarian and colorectal cancers, thus showing ROC AUC values > 0.7. In currently ongoing prospective trial for advanced solid tumor patients, the Oncobox tests were completed for 239 patients, primary feedback information received for 144 patients. 25 patients (17%) died before prescription of the therapy, 19% received palliative care treatment, 39% received Oncobox-recommended therapies and 25% received other therapies (February 2019). Tumor responses were estimated for 30 patients (breast, colorectal, ovarian and other cancers) receiving therapies, with disease control rate of 71% for Oncobox-recommended ATDs (Bevacizumab, Crizotinib, Trastuzumab and others) and 44% for other therapies. Our results also suggest that cancer metastases and primary tumors may have different gene expression, molecular pathway activation and drug scoring patterns, thus pointing to the importance of testing multiple tumor sites. Conclusions: Our study suggests that profound analysis of gene expression profiles of tumor tissues may be useful for ATDs treatment of advanced and metastatic cancers. Clinical trial information: NCT03724097.

RNA sequencing analysis for profiling activation of cancer-associated molecular pathways

e13032

Background: Intracellular molecular pathways (IMPs) control all major events in the living cell. They are considered hotspots in contemporary oncology because knowledge of IMPs activation is essential for understanding mechanisms of molecular pathogenesis in oncology. Profiling IMPs requires RNA-seq data for tumors and for a collection of reference normal tissues. However, there is a shortage now in such profiles for normal tissues from healthy human donors, uniformly profiled in a single series of experiments. Access to the largest dataset of normal profiles GTEx is only partly available through the dbGaP. In TCGA database, norms are adjacent to surgically removed tumors and may be affected by tumor-linked growth factors, inflammation and altered vascularization. ENCODE datasets were for the autopsies of normal tissues, but they can’t form statistically significant reference groups. Methods: Tissue samples representing 20 organs were taken from post-mortal human healthy donors killed in road accidents no later than 36 hours after death, blood samples were taken from healthy volunteers. Gene expression was profiled in RNA-seq experiments using the same reagents, equipment and protocols. Bioinformatic algorithms for IMP analysis were developed and validated using experimental and public gene expression datasets. Results: From original sequencing data we constructed the biggest fully open reference expression database of normal human tissues including 465 profiles termed Oncobox Atlas of Normal Tissue Expression (ANTE, original data: GSE120795). We next developed a method termed Oncobox for interrogating activation of IMPs in human cancers. It includes modules of expression data harmonization and comparison and an algorithm for automatic annotation of molecular pathways. The Oncobox system enables accurate scoring of thousands molecular pathways using RNA-seq data. Oncobox pathway analysis is also applicable for quantitative proteomics and microRNA data in oncology. Conclusions: The Oncobox system can be used for a plethora of applications in cancer research including finding differentially regulated genes and IMPs, and for discovery of new pathway-related diagnostic and prognostic biomarkers.

Atlas of RNA sequencing profiles for normal human tissues

Comprehensive analysis of molecular pathology requires a collection of reference samples representing normal tissues from healthy donors. For the available limited collections of normal tissues from postmortal donors, there is a problem of data incompatibility, as different datasets generated using different experimental platforms often cannot be merged in a single panel. Here, we constructed and deposited the gene expression database of normal human tissues based on uniformly screened original sequencing data. In total, 142 solid tissue samples representing 20 organs were taken from post-mortal human healthy donors of different age killed in road accidents no later than 36 hours after death. Blood samples were taken from 17 healthy volunteers. We then compared them with the 758 transcriptomic profiles taken from the other databases. We found that overall 463 biosamples showed tissue-specific rather than platform- or database-specific clustering and could be aggregated in a single database termed Oncobox Atlas of Normal Tissue Expression (ANTE). Our data will be useful to all those working with the analysis of human gene expression.

Gene expression and molecular pathway activation signatures of MYCN-amplified neuroblastomas

Neuroblastoma is a pediatric cancer arising from sympathetic nervous system. Remarkable heterogeneity in outcomes is one of its widely known features. One of the traits strongly associated with the unfavorable subtype is the amplification of oncogene MYCN. Here, we performed cross-platform biomarker detection by comparing gene expression and pathway activation patterns from the two literature reports and from our experimental dataset, combining profiles for the 761 neuroblastoma patients with known MYCN amplification status. We identified 109 / 25 gene expression / pathway activation biomarkers strongly linked with the MYCN amplification. The marker genes/pathways are involved in the processes of purine nucleotide biosynthesis, ATP-binding, tetrahydrofolate metabolism, building mitochondrial matrix, biosynthesis of amino acids, tRNA aminoacylation and NADP-linked oxidation-reduction processes, as well as in the tyrosine phosphatase activity, p53 signaling, cell cycle progression and the G1/S and G2/M checkpoints. To connect molecular functions of the genes involved in MYCN-amplified phenotype, we built a new molecular pathway using known intracellular protein interaction networks. The activation of this pathway was highly selective in discriminating MYCN-amplified neuroblastomas in all three datasets. Our data also suggest that the phosphoinositide 3-kinase (PI3K) inhibitors may provide new opportunities for the treatment of the MYCN-amplified neuroblastoma subtype.

Activation of intracellular signaling pathways as a new type of biomarkers for selection of target anticancer drugs

e23142

Background: Anticancer target drugs (ATDs) specifically bind and inhibit molecular targets that play important roles in tumorigenesis. More than 150 different ATDs have been approved for clinical use worldwide, and the clinicians are faced with the problem of choosing the best therapeutic solution for each patient. The problem of efficient ATD selection remains largely unsolved and personalized approaches are needed to select the best ATD candidates for individual patients. Methods: We propose a new approach termed OncoFinder. It is based on digesting gene expression profiles for the analysis of activation of intracellular signalling pathways as a marker for the selection of target therapies. The original bioinformatic algorithms were integrated with the databases featuring molecular drug targets, compositions of signalling pathways, including the functional role of each gene product, for more than 1700 pathways (Buzdin, Front.Genet 2014; Ozerov, Nature Communications 2016). Results: We showed that pathway activation strengths are more stable and reliable biomarkers of cancer than the expressions of individual genes. OncoFinder allows to detect changes at the level of pathway activation and to predict the effectiveness of drugs based on the knowledge of their molecular targets. We applied it to find new biomarkers of clinical response to the ATD cetuximab; for modelling the combined chemotherapy of acute myeloid leukemia and combined anti-VEGF/BRAF therapy of melanoma. For two unrelated datasets obtained for colon cancer patients before treatment with the ATD bevacizumab, we were able to distinguish between those who responded to treatment and not (p < 0.01). We next assayed biopsies for kidney cancer patients with known responses to the ATD sorafenib. The responders and non-responders showed a significant difference (p = 0.02). Finally, the OncoFinder platform was prospectively used for decision making support to patients with advanced metastatic solid tumors (n = 23). The efficiency of the ATD treatment was 61% (complete + partial response, RECIST). Conclusions: OncoFinder method may be effective for predicting response to ATD based on high throughput gene expression profiles.

MiRImpact, a new bioinformatic method using complete microRNA expression profiles to assess their overall influence on the activity of intracellular molecular pathways

MicroRNAs (miRs) are short noncoding RNA molecules that regulate expression of target mRNAs. Many published sources provide information about miRs and their targets. However, bioinformatic tools elucidating higher level impact of the established total miR profiles, are still largely missing. Recently, we developed a method termed OncoFinder enabling quantification of the activities of intracellular molecular pathways basing on gene expression data. Here we propose a new technique, MiRImpact, which enables to link miR expression data with its estimated outcome on the regulation of molecular pathways, like signaling, metabolic, cytoskeleton rearrangement, and DNA repair pathways. MiRImpact uses OncoFinder rationale for pathway activity calculations, with the major distinctions that (i) it deals with the concentrations of miRs - known regulators of gene products participating in molecular pathways, and (ii) miRs are considered as negative regulators of target molecules, if other is not specified. MiRImpact operates with 2 types of databases: for molecular targets of miRs and for gene products participating in molecular pathways. We applied MiRImpact to compare regulation of human bladder cancer-specific signaling pathways at the levels of mRNA and miR expression. We took 2 most complete alternative databases of experimentally validated miR targets – miRTarBase and DianaTarBase, and an OncoFinder database featuring 2725 gene products and 271 signaling pathways. We showed that the impact of miRs is orthogonal to pathway regulation at the mRNA level, which stresses the importance of studying posttranscriptional regulation of gene expression. We also report characteristic set of miR and mRNA regulation features linked with bladder cancer.

Molecular pathway activation features linked with transition from normal skin to primary and metastatic melanomas in human

Melanoma is the most aggressive and dangerous type of skin cancer, but its molecular mechanisms remain largely unclear. For transcriptomic data of 478 primary and metastatic melanoma, nevi and normal skin samples, we performed high-throughput analysis of intracellular molecular networks including 592 signaling and metabolic pathways. We showed that at the molecular pathway level, the formation of nevi largely resembles transition from normal skin to primary melanoma. Using a combination of bioinformatic machine learning algorithms, we identified 44 characteristic signaling and metabolic pathways connected with the formation of nevi, development of primary melanoma, and its metastases. We created a model describing formation and progression of melanoma at the level of molecular pathway activation. We discovered six novel associations between activation of metabolic molecular pathways and progression of melanoma: for allopregnanolone biosynthesis, L-carnitine biosynthesis, zymosterol biosynthesis (inhibited in melanoma), fructose 2, 6-bisphosphate synthesis and dephosphorylation, resolvin D biosynthesis (activated in melanoma), D-myo-inositol hexakisphosphate biosynthesis (activated in primary, inhibited in metastatic melanoma). Finally, we discovered fourteen tightly coordinated functional clusters of molecular pathways. This study helps to decode molecular mechanisms underlying the development of melanoma.

A method for predicting target drug efficiency in cancer based on the analysis of signaling pathway activation

A new generation of anticancer therapeutics called target drugs has quickly developed in the 21st century. These drugs are tailored to inhibit cancer cell growth, proliferation, and viability by specific interactions with one or a few target proteins. However, despite formally known molecular targets for every "target" drug, patient response to treatment remains largely individual and unpredictable. Choosing the most effective personalized treatment remains a major challenge in oncology and is still largely trial and error. Here we present a novel approach for predicting target drug efficacy based on the gene expression signature of the individual tumor sample(s). The enclosed bioinformatic algorithm detects activation of intracellular regulatory pathways in the tumor in comparison to the corresponding normal tissues. According to the nature of the molecular targets of a drug, it predicts whether the drug can prevent cancer growth and survival in each individual case by blocking the abnormally activated tumor-promoting pathways or by reinforcing internal tumor suppressor cascades. To validate the method, we compared the distribution of predicted drug efficacy scores for five drugs (Sorafenib, Bevacizumab, Cetuximab, Sorafenib, Imatinib, Sunitinib) and seven cancer types (Clear Cell Renal Cell Carcinoma, Colon cancer, Lung adenocarcinoma, non-Hodgkin Lymphoma, Thyroid cancer and Sarcoma) with the available clinical trials data for the respective cancer types and drugs. The percent of responders to a drug treatment correlated significantly (Pearson's correlation 0.77 p = 0.023) with the percent of tumors showing high drug scores calculated with the current algorithm.

nMETR: technique for facile recovery of hypomethylation genomic tags

Genome-wide methylation studies frequently lack adequate controls to estimate proportions of background reads in the resulting datasets. To generate appropriate control pools, we developed technique termed nMETR (non-methylated tag recovery) based on digestion of genomic DNA with methylation-sensitive restriction enzyme, ligation of adapter oligonucleotide and PCR amplification of non-methylated sites associated with genomic repetitive elements. The protocol takes only two working days to generate amplicons for deep sequencing. We applied nMETR for human DNA using BspFNI enzyme and retrotransposon Alu-specific primers. 454-sequencing enabled identification of 1113 nMETR tag sites, of them ~65% were parts of CpG islands. Representation of reads inversely correlated with methylation levels, thus confirming nMETR fidelity. We created software that eliminates background reads and enables to map and annotate individual tags on human genome. nMETR tags may serve as the controls for large-scale epigenetic studies and for identifying unmethylated transposable elements located close to genomic CpG islands.