Novel robust biomarkers for human bladder cancer based on activation of intracellular signaling pathways

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.

Effects and Mechanisms of Kaempferol in the Management of Cancers through Modulation of Inflammation and Signal Transduction Pathways

Cancer is the principal cause of death and its incidence is increasing continuously worldwide. Various treatment approaches are in practice to treat cancer, but these treatment strategies may be associated with severe side effects and also produce drug resistance. However, natural compounds have established their role in cancer management with minimal side effects. In this vista, kaempferol, a natural polyphenol, mainly found in vegetables and fruits, has been revealed to have many health-promoting effects. Besides its health-promoting potential, its anti-cancer potential has also been described in in vivo as well as in in vitro studies. The anti-cancer potential of kaempferol has been proven through modulation of cell signaling pathways in addition to the induction of apoptosis and cell cycle arrest in cancer cells. It leads to the activation of tumor suppressor genes, inhibition of angiogenesis, PI3K/AKT pathways, STAT3, transcription factor AP-1, Nrf2 and other cell signaling molecules. Poor bioavailability of this compound is one of the major limitations for its proper and effective disease management actions. Recently, some novel nanoparticle-based formulations have been used to overcome these limitations. The aim of this review is to provide a clear picture regarding the mechanism of action of kaempferol in different cancers through the modulation of cell signaling molecules. Besides this, strategies to improve the efficacy and synergistic effects of this compound have also been described. However, more studies are needed based on clinical trials to fully explore the therapeutic role of this compound, especially in cancer treatment.

Personalized targeted therapy prescription in colorectal cancer using algorithmic analysis of RNA sequencing data

BACKGROUND

Overall survival of advanced colorectal cancer (CRC) patients remains poor, and gene expression analysis could potentially complement detection of clinically relevant mutations to personalize CRC treatments.

METHODS

We performed RNA sequencing of formalin-fixed, paraffin-embedded (FFPE) cancer tissue samples of 23 CRC patients and interpreted the data obtained using bioinformatic method Oncobox for expression-based rating of targeted therapeutics. Oncobox ranks cancer drugs according to the efficiency score calculated using target genes expression and molecular pathway activation data. The patients had primary and metastatic CRC with metastases in liver, peritoneum, brain, adrenal gland, lymph nodes and ovary. Two patients had mutations in NRAS, seven others had mutated KRAS gene. Patients were treated by aflibercept, bevacizumab, bortezomib, cabozantinib, cetuximab, crizotinib, denosumab, panitumumab and regorafenib as monotherapy or in combination with chemotherapy, and information on the success of totally 39 lines of therapy was collected.

RESULTS

Oncobox drug efficiency score was effective biomarker that could predict treatment outcomes in the experimental cohort (AUC 0.77 for all lines of therapy and 0.91 for the first line after tumor sampling). Separately for bevacizumab, it was effective in the experimental cohort (AUC 0.87) and in 3 independent literature CRC datasets, n = 107 (AUC 0.84-0.94). It also predicted progression-free survival in univariate (Hazard ratio 0.14) and multivariate (Hazard ratio 0.066) analyses. Difference in AUC scores evidences importance of using recent biosamples for the prediction quality.

CONCLUSION

Our results suggest that RNA sequencing analysis of tumor FFPE materials may be helpful for personalizing prescriptions of targeted therapeutics in CRC.

Exploring the recent trends in perturbing the cellular signaling pathways in cancer by natural products

Cancer is commonly thought to be the product of irregular cell division. According to the World Health Organization (WHO), cancer is the major cause of death globally. Nature offers an abundant supply of bioactive compounds with high therapeutic efficacy. Anticancer effects have been studied in a variety of phytochemicals found in nature. When Food and Drug Administration (FDA)-approved anticancer drugs are combined with natural compounds, the effectiveness improves. Several agents have already progressed to clinical trials based on these promising results of natural compounds against various cancer forms. Natural compounds prevent cancer cell proliferation, development, and metastasis by inducing cell cycle arrest, activating intrinsic and extrinsic apoptosis pathways, generating reactive oxygen species (ROS), and down-regulating activated signaling pathways. These natural chemicals are known to affect numerous important cellular signaling pathways, such as NF-B, MAPK, Wnt, Notch, Akt, p53, AR, ER, and many others, to cause cell death signals and induce apoptosis in pre-cancerous or cancer cells without harming normal cells. As a result, non-toxic “natural drugs” taken from nature’s bounty could be effective for the prevention of tumor progression and/or therapy of human malignancies, either alone or in combination with conventional treatments. Natural compounds have also been shown in preclinical studies to improve the sensitivity of resistant cancers to currently available chemotherapy agents. To summarize, preclinical and clinical findings against cancer indicate that natural-sourced compounds have promising anticancer efficacy. The vital purpose of these studies is to target cellular signaling pathways in cancer by natural compounds.

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.

Using proteomic and transcriptomic data to assess activation of intracellular molecular pathways

Analysis of molecular pathway activation is the recent instrument that helps to quantize activities of various intracellular signaling, structural, DNA synthesis and repair, and biochemical processes. This may have a deep impact in fundamental research, bioindustry, and medicine. Unlike gene ontology analyses and numerous qualitative methods that can establish whether a pathway is affected in principle, the quantitative approach has the advantage of exactly measuring the extent of a pathway up/downregulation. This results in emergence of a new generation of molecular biomarkers-pathway activation levels, which reflect concentration changes of all measurable pathway components. The input data can be the high-throughput proteomic or transcriptomic profiles, and the output numbers take both positive and negative values and positively reflect overall pathway activation. Due to their nature, the pathway activation levels are more robust biomarkers compared to the individual gene products/protein levels. Here, we review the current knowledge of the quantitative gene expression interrogation methods and their applications for the molecular pathway quantization. We consider enclosed bioinformatic algorithms and their applications for solving real-world problems. Besides a plethora of applications in basic life sciences, the quantitative pathway analysis can improve molecular design and clinical investigations in pharmaceutical industry, can help finding new active biotechnological components and can significantly contribute to the progressive evolution of personalized medicine. In addition to the theoretical principles and concepts, we also propose publicly available software for the use of large-scale protein/RNA expression data to assess the human pathway activation levels.

Algorithmic Annotation of Functional Roles for Components of 3,044 Human Molecular Pathways

Current methods of high-throughput molecular and genomic analyses enabled to reconstruct thousands of human molecular pathways. Knowledge of molecular pathways structure and architecture taken along with the gene expression data can help interrogating the pathway activation levels (PALs) using different bioinformatic algorithms. In turn, the pathway activation profiles can characterize molecular processes, which are differentially regulated and give numeric characteristics of the extent of their activation or inhibition. However, different pathway nodes may have different functions toward overall pathway regulation, and calculation of PAL requires knowledge of molecular function of every node in the pathway in terms of its activator or inhibitory role. Thus, high-throughput annotation of functional roles of pathway nodes is required for the comprehensive analysis of the pathway activation profiles. We proposed an algorithm that identifies functional roles of the pathway components and applied it to annotate 3,044 human molecular pathways extracted from the Biocarta, Reactome, KEGG, Qiagen Pathway Central, NCI, and HumanCYC databases and including 9,022 gene products. The resulting knowledgebase can be applied for the direct calculation of the PALs and establishing large scale profiles of the signaling, metabolic, and DNA repair pathway regulation using high throughput gene expression data. We also provide a bioinformatic tool for PAL data calculations using the current pathway knowledgebase.

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 profiles and diagnostic signatures linked with response to ramucirumab in gastric cancer

Gastric cancer (GC) is the fifth-ranked cancer type by associated mortality. The proportion of early diagnosis is low, and most patients are diagnosed at the advanced stages. First-line therapy standardly includes fluoropyrimidines and platinum compounds with trastuzumab for HER2-positive cases. For recurrent disease, there are several alternative options including ramucirumab, a monoclonal therapeutic antibody that inhibits VEGF-mediated tumor angiogenesis by binding with VEGFR2, alone or in combination with other cancer drugs. However, overall response rate following ramucirumab or its combinations is 30%-80% of the patients, suggesting that personalization of drug prescription is needed to increase efficacy of treatment. We report here original tumor RNA sequencing profiles for 15 advanced GC patients linked with data on clinical response to ramucirumab or its combinations. Three genes showed differential expression in the tumors for responders versus nonresponders: CHRM3, LRFN1, and TEX15 Of them, CHRM3 was up-regulated in the responders. Using the bioinformatic platform Oncobox we simulated ramucirumab efficiency and compared output model results with actual tumor response data. An agreement was observed between predicted and real clinical outcomes (AUC ≥ 0.7). These results suggest that RNA sequencing may be used to personalize the prescription of ramucirumab for GC and indicate potential molecular mechanisms underlying ramucirumab resistance. The RNA sequencing profiles obtained here are fully compatible with the previously published Oncobox Atlas of Normal Tissue Expression (ANTE) data.

Gene signatures associated with genomic aberrations predict prognosis in neuroblastoma

BACKGROUND

Neuroblastoma (NB) is a heterogeneous disease with respect to genomic abnormalities and clinical behaviors. Despite recent advances in our understanding of the association between the genetic aberrations and clinical features, it remains one of the major challenges to predict prognosis and stratify patients for determining personalized therapy in this disease. The aim of this study was to develop an effective prognosis prediction model for NB patients.

METHODS

We integrated diverse computational analyses to define gene signatures that reflect MYCN activity and chromosomal aberrations including deletion of chromosome 1p (Chr1p_del) and chromosome 11q (Chr11q_del) as well as chromosome 11q whole loss (Chr11q_wls). We evaluated the prognostic and predictive values of these signatures in seven NB gene expression datasets (the number of samples ranges from 94 to 498, with a total of 2120) generated from both RNA sequencing and microarray platforms.

RESULTS

MYCN signature was a more effective prognostic marker than MYCN amplification status and MYCN expression. Similarly, the Chr1p_del score was more prognostic than Chr1p status. The activity scores of MYCN, Chr1p_del and Chr11q_del were associated with poor prognosis, while the Chr11q_wls score was linked to good outcome. We integrated the activity scores of MYCN, Chr1p_del, Chr11q_del, and Chr11q_wls and clinical variables into an integrative prognostic model, which displayed significant performance over the clinical variables or each genomic aberration alone.

CONCLUSIONS

Our integrative gene signature model shows a significantly improved forecast performance with prognostic and predictive information, and thereby can be served as a biomarker to stratify NB patients for prognosis evaluation and surveillance programs.

Modulation of Multiple Signaling Pathways of the Plant-Derived Natural Products in Cancer

Natural compounds are highly effective anticancer chemotherapeutic agents, and the targets of plant-derived anticancer agents have been widely reported. In this review, we focus on the main signaling pathways of apoptosis, proliferation, invasion, and metastasis that are regulated by polyphenols, alkaloids, saponins, and polysaccharides. Alkaloids primarily affect apoptosis-related pathways, while polysaccharides primarily target pathways related to proliferation, invasion, and metastasis. Other compounds, such as flavonoids and saponins, affect all of these aspects. The association between compound structures and signaling pathways may play a critical role in drug discovery.

Comparative Bioinformatics Characteristic of Bladder Cancer Stage 2 from Stage 4 Expression Profile: A Network-Based Study

BACKGROUND

Bladder cancer (BC) has remained as one of the most challenging issues in medicine. The aim of this study was to investigate the differential network analysis of stages 2 and 4 of BC to better understand the molecular pathology of these states.

MATERIALS AND METHODS

We chose gene expression data of GSE52519 from Gene Expression Omnibus (GEO) database analyzed by the GEO2R online tool. Cytoscape version 3.6.1 and its algorithms are the methods applied for the network construction and investigation of differentially expressed genes (DEG) in these states.

RESULT

Our result revealed that the analysis DEGs provides useful information about a common molecular feature of stages 2 and 4 of BC.

CONCLUSION

Consequently, the network finding revealed that more investigation about stage 2 is required to achieve an effective therapeutic protocol to block the transition from stage 2 to stage 4.

Chemo brain or tumor brain - that is the question: the presence of extracranial tumors profoundly affects molecular processes in the prefrontal cortex of TumorGraft mice

Cancer chemotherapy causes numerous persistent central nervous system complications. This condition is known as chemo brain. Cognitive impairments occur even before treatment, and hence are referred to as cancer associated cognitive changes, or tumor brain. There is much yet to be learned about the mechanisms of both chemo brain and tumor brain. The frequency and timing of chemo brain and tumor brain occurrence and persistence strongly suggest they may be epigenetic in nature and associated with altered gene expression. Here we used TumorGraft^TM models wherein part of a patient's tumor is removed and grafted into immune-deficient mice and conducted global gene expression and DNA methylation analysis. We show that malignant non-central nervous system tumor growth causes profound molecular alterations in the brain. Mice harbouring triple negative or progesterone positive breast cancer TumorGrafts exhibited altered gene expression, decreased levels of DNA methylation, increased levels of DNA hydroxymethylation, and oxidative stress in the prefrontal cortex. Interestingly, chemotherapy did not have any additional synergistic effects on the analyzed processes. The molecular changes observed in this study are known signs of neurodegeneration and brain aging. This study provides an important roadmap for future large-scale analysis of the molecular and cellular mechanisms of tumor brain.

Acquired resistance to tyrosine kinase inhibitors may be linked with the decreased sensitivity to X-ray irradiation

Acquired resistance to chemotherapy and radiation therapy is one of the major obstacles decreasing efficiency of treatment of the oncologic diseases. In this study, on the two cell lines (ovarian carcinoma SKOV-3 and neuroblastoma NGP-127), we modeled acquired resistance to five target anticancer drugs. The cells were grown on gradually increasing concentrations of the clinically relevant tyrosine kinase inhibitors (TKIs) Sorafenib, Pazopanib and Sunitinib, and rapalogs Everolimus and Temsirolimus, for 20 weeks. After 20 weeks of culturing, the half-inhibitory concentrations (IC₅₀) increased by 25 – 186% for the particular combinations of the drugs and cell types. We next subjected cells to 10 Gy irradiation, a dose frequently used in clinical radiation therapy. For the SKOV-3, but not NGP-127 cells, for the TKIs Sorafenib, Pazopanib and Sunitinib, we noticed statistically significant increase in capacity to repair radiation-induced DNA double strand breaks compared to naïve control cells not previously treated with TKIs. These peculiarities were linked with the increased activation of ATM DNA repair pathway in the TKI-treated SKOV-3, but not NGP-127 cells. Our results provide a new cell culture model for studying anti-cancer therapy efficiency and evidence that there may be a tissue-specific radioresistance emerging as a side effect of treatment with TKIs.

Towards natural mimetics of metformin and rapamycin

Aging is now at the forefront of major challenges faced globally, creating an immediate need for safe, widescale interventions to reduce the burden of chronic disease and extend human healthspan. Metformin and rapamycin are two FDA-approved mTOR inhibitors proposed for this purpose, exhibiting significant anti-cancer and anti-aging properties beyond their current clinical applications. However, each faces issues with approval for off-label, prophylactic use due to adverse effects. Here, we initiate an effort to identify nutraceuticals-safer, naturally-occurring compounds-that mimic the anti-aging effects of metformin and rapamycin without adverse effects. We applied several bioinformatic approaches and deep learning methods to the Library of Integrated Network-based Cellular Signatures (LINCS) dataset to map the gene- and pathway-level signatures of metformin and rapamycin and screen for matches among over 800 natural compounds. We then predicted the safety of each compound with an ensemble of deep neural network classifiers. The analysis revealed many novel candidate metformin and rapamycin mimetics, including allantoin and ginsenoside (metformin), epigallocatechin gallate and isoliquiritigenin (rapamycin), and withaferin A (both). Four relatively unexplored compounds also scored well with rapamycin. This work revealed promising candidates for future experimental validation while demonstrating the applications of powerful screening methods for this and similar endeavors.

Data aggregation at the level of molecular pathways improves stability of experimental transcriptomic and proteomic data

High throughput technologies opened a new era in biomedicine by enabling massive analysis of gene expression at both RNA and protein levels. Unfortunately, expression data obtained in different experiments are often poorly compatible, even for the same biologic samples. Here, using experimental and bioinformatic investigation of major experimental platforms, we show that aggregation of gene expression data at the level of molecular pathways helps to diminish cross- and intra-platform bias otherwise clearly seen at the level of individual genes. We created a mathematical model of cumulative suppression of data variation that predicts the ideal parameters and the optimal size of a molecular pathway. We compared the abilities to aggregate experimental molecular data for the 5 alternative methods, also evaluated by their capacity to retain meaningful features of biologic samples. The bioinformatic method OncoFinder showed optimal performance in both tests and should be very useful for future cross-platform data analyses.

The cornucopia of meaningful leads: Applying deep adversarial autoencoders for new molecule development in oncology

Recent advances in deep learning and specifically in generative adversarial networks have demonstrated surprising results in generating new images and videos upon request even using natural language as input. In this paper we present the first application of generative adversarial autoencoders (AAE) for generating novel molecular fingerprints with a defined set of parameters. We developed a 7-layer AAE architecture with the latent middle layer serving as a discriminator. As an input and output the AAE uses a vector of binary fingerprints and concentration of the molecule. In the latent layer we also introduced a neuron responsible for growth inhibition percentage, which when negative indicates the reduction in the number of tumor cells after the treatment. To train the AAE we used the NCI-60 cell line assay data for 6252 compounds profiled on MCF-7 cell line. The output of the AAE was used to screen 72 million compounds in PubChem and select candidate molecules with potential anti-cancer properties. This approach is a proof of concept of an artificially-intelligent drug discovery engine, where AAEs are used to generate new molecular fingerprints with the desired molecular properties.

Chemo brain: From discerning mechanisms to lifting the brain fog-An aging connection

Mounting evidence indicates that cancer treatments cause numerous deleterious effects, including central nervous system (CNS) toxicity. Chemotherapy-caused CNS side effects encompass changes in cognitive function, memory, and attention, to name a few. Although chemotherapy treatment-induced side effects occur in 16-75% of all patients, the mechanisms of these effects are not well understood. We have recently proposed a new epigenetic theory of chemo brain and, in a pioneer study, determined that cytotoxic chemotherapy agents induce oxidative DNA damage and affect molecular and epigenetic processes in the brain, and may be associated with brain aging processes. In this paper, we discuss the implications of chemo brain epigenetic effects and future perspectives, as well as outline potential links with brain aging and future translational research opportunities.

An analysis of gene expression data involving examination of signaling pathways activation reveals new insights into the mechanism of action of minoxidil topical foam in men with androgenetic alopecia

Androgenetic alopecia is the most common form of hair loss. Minoxidil has been approved for the treatment of hair loss, however its mechanism of action is still not fully clarified. In this study, we aimed to elucidate the effects of 5% minoxidil topical foam on gene expression and activation of signaling pathways in vertex and frontal scalp of men with androgenetic alopecia. We identified regional variations in gene expression and perturbed signaling pathways using in silico Pathway Activation Network Decomposition Analysis (iPANDA) before and after treatment with minoxidil. Vertex and frontal scalp of patients showed a generally similar response to minoxidil. Both scalp regions showed upregulation of genes that encode keratin associated proteins and downregulation of ILK, Akt, and MAPK signaling pathways after minoxidil treatment. Our results provide new insights into the mechanism of action of minoxidil topical foam in men with androgenetic alopecia.

In silico analysis of pathways activation landscape in oral squamous cell carcinoma and oral leukoplakia

A subset of patients with oral squamous cell carcinoma (OSCC), the most common subtype of head and neck squamous cell carcinoma (HNSCC), harbor dysplastic lesions (often visually identified as leukoplakia) prior to cancer diagnosis. Although evidence suggest that leukoplakia represents an initial step in the progression to cancer, signaling networks driving this progression are poorly understood. Here, we applied in silico Pathway Activation Network Decomposition Analysis (iPANDA), a new bioinformatics software suite for qualitative analysis of intracellular signaling pathway activation using transcriptomic data, to assess a network of molecular signaling in OSCC and pre-neoplastic oral lesions. In tumor samples, our analysis detected major conserved mitogenic and survival signaling pathways strongly associated with HNSCC, suggesting that some of the pathways identified by our algorithm, but not yet validated as HNSCC related, may be attractive targets for future research. While pathways activation landscape in the majority of leukoplakias was different from that seen in OSCC, a subset of pre-neoplastic lesions has demonstrated some degree of similarity to the signaling profile seen in tumors, including dysregulation of the cancer-driving pathways related to survival and apoptosis. These results suggest that dysregulation of these signaling networks may be the driving force behind the early stages of OSCC tumorigenesis. While future studies with larger leukoplakia data sets are warranted to further estimate the values of this approach for capturing signaling features that characterize relevant lesions that actually progress to cancers, our platform proposes a promising new approach for detecting cancer-promoting pathways and tailoring the right therapy to prevent tumorigenesis.

Bioinformatics Meets Biomedicine: OncoFinder, a Quantitative Approach for Interrogating Molecular Pathways Using Gene Expression Data

We propose a biomathematical approach termed OncoFinder (OF) that enables performing both quantitative and qualitative analyses of the intracellular molecular pathway activation. OF utilizes an algorithm that distinguishes the activator/repressor role of every gene product in a pathway. This method is applicable for the analysis of any physiological, stress, malignancy, and other conditions at the molecular level. OF showed a strong potential to neutralize background-caused differences between experimental gene expression data obtained using NGS, microarray and modern proteomics techniques. Importantly, in most cases, pathway activation signatures were better markers of cancer progression compared to the individual gene products. OF also enables correlating pathway activation with the success of anticancer therapy for individual patients. We further expanded this approach to analyze impact of micro RNAs (miRs) on the regulation of cellular interactome. Many alternative sources provide information about miRs and their targets. However, instruments elucidating higher level impact of the established total miR profiles are still largely missing. A variant of OncoFinder termed MiRImpact enables linking miR expression data with its estimated outcome on the regulation of molecular processes, such as signaling, metabolic, cytoskeleton, and DNA repair pathways. MiRImpact was used to establish cancer-specific and cytomegaloviral infection-linked interactomic signatures for hundreds of molecular pathways. Interestingly, the impact of miRs appeared orthogonal to pathway regulation at the mRNA level, which stresses the importance of combining all available levels of gene regulation to build a more objective molecular model of cell.

Mathematical Justification of Expression-Based Pathway Activation Scoring (PAS)

Although modeling of activation kinetics for various cell signaling pathways has reached a high grade of sophistication and thoroughness, most such kinetic models still remain of rather limited practical value for biomedicine. Nevertheless, recent advancements have been made in application of signaling pathway science for real needs of prescription of the most effective drugs for individual patients. The methods for such prescription evaluate the degree of pathological changes in the signaling machinery based on two types of data: first, on the results of high-throughput gene expression profiling, and second, on the molecular pathway graphs that reflect interactions between the pathway members. For example, our algorithm OncoFinder evaluates the activation of molecular pathways on the basis of gene/protein expression data in the objects of the interest.Yet, the question of assessment of the relative importance for each gene product in a molecular pathway remains unclear unless one call for the methods of parameter sensitivity /stiffness analysis in the interactomic kinetic models of signaling pathway activation in terms of total concentrations of each gene product.Here we show two principal points: 1. First, the importance coefficients for each gene in pathways that were obtained using the extremely time- and labor-consuming stiffness analysis of full-scaled kinetic models generally differ from much easier-to-calculate expression-based pathway activation score (PAS) not more than by 30%, so the concept of PAS is kinetically justified. 2. Second, the use of pathway-based approach instead of distinct gene analysis, due to the law of large numbers, allows restoring the correlation between the similar samples that were examined using different transcriptome investigation techniques.

Early stage of cytomegalovirus infection suppresses host microRNA expression regulation in human fibroblasts

Responses to human cytomegalovirus (HCMV) infection are largely individual and cell type specific. We investigated molecular profiles in 2 primary cell cultures of human fibroblasts, which are highly or marginally sensitive to HCMV infection, respectively. We screened expression of genes and microRNAs (miRs) at the early (3 hours) stage of infection. To assess molecular pathway activation profiles, we applied bioinformatic algorithms OncoFinder and MiRImpact. In both cell types, pathway regulation properties at mRNA and miR levels were markedly different. Surprisingly, in the infected highly sensitive cells, we observed a "freeze" of miR expression profiles compared to uninfected controls. Our results evidence that in the sensitive cells, HCMV blocks intracellular regulation of microRNA expression already at the earliest stage of infection. These data suggest somewhat new functions for HCMV products and demonstrate dependence of miR expression arrest on the host-encoded factors.

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.

Towards understanding the breast cancer epigenome: a comparison of genome-wide DNA methylation and gene expression data

Until recently, an elevated disease risk has been ascribed to a genetic predisposition, however, exciting progress over the past years has discovered alternate elements of inheritance that involve epigenetic regulation. Epigenetic changes are heritably stable alterations that include DNA methylation, histone modifications and RNA-mediated silencing. Aberrant DNA methylation is a common molecular basis for a number of important human diseases, including breast cancer. Changes in DNA methylation profoundly affect global gene expression patterns. What is emerging is a more dynamic and complex association between DNA methylation and gene expression than previously believed. Although many tools have already been developed for analyzing genome-wide gene expression data, tools for analyzing genome-wide DNA methylation have not yet reached the same level of refinement.

Here we provide an in-depth analysis of DNA methylation in parallel with gene expression data characteristics and describe the particularities of low-level and high-level analyses of DNA methylation data. Low-level analysis refers to pre-processing of methylation data (i.e. normalization, transformation and filtering), whereas high-level analysis is focused on illustrating the application of the widely used class comparison, class prediction and class discovery methods to DNA methylation data. Furthermore, we investigate the influence of DNA methylation on gene expression by measuring the correlation between the degree of CpG methylation and the level of expression and to explore the pattern of methylation as a function of the promoter region.

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.

Combinatorial high-throughput experimental and bioinformatic approach identifies molecular pathways linked with the sensitivity to anticancer target drugs

Effective choice of anticancer drugs is important problem of modern medicine. We developed a method termed OncoFinder for the analysis of new type of biomarkers reflecting activation of intracellular signaling and metabolic molecular pathways. These biomarkers may be linked with the sensitivity to anticancer drugs. In this study, we compared the experimental data obtained in our laboratory and in the Genomics of Drug Sensitivity in Cancer (GDS) project for testing response to anticancer drugs and transcriptomes of various human cell lines. The microarray-based profiling of transcriptomes was performed for the cell lines before the addition of drugs to the medium, and experimental growth inhibition curves were built for each drug, featuring characteristic IC₅₀ values. We assayed here four target drugs - Pazopanib, Sorafenib, Sunitinib and Temsirolimus, and 238 different cell lines, of which 11 were profiled in our laboratory and 227 - in GDS project. Using the OncoFinder-processed transcriptomic data on ∼600 molecular pathways, we identified pathways showing significant correlation between pathway activation strength (PAS) and IC₅₀ values for these drugs. Correlations reflect relationships between response to drug and pathway activation features. We intersected the results and found molecular pathways significantly correlated in both our assay and GDS project. For most of these pathways, we generated molecular models of their interaction with known molecular target(s) of the respective drugs. For the first time, our study uncovered mechanisms underlying cancer cell response to drugs at the high-throughput molecular interactomic level.

ELK1 is up-regulated by androgen in bladder cancer cells and promotes tumor progression

Little is known about biological significance of ELK1, a transcriptional factor that activates downstream targets including c-fos proto-oncogene, in bladder cancer. Recent preclinical evidence also suggests the involvement of androgen receptor (AR) signaling in bladder cancer progression. In this study, we aim to investigate the functions of ELK1 in bladder cancer growth and their regulation by AR signals. Immunohistochemistry in bladder tumor specimens showed that the levels of phospho-ELK1 (p-ELK1) expression were significantly elevated in urothelial neoplasms, compared with non-neoplastic urothelium tissues, and were also correlated with AR positivity. Patients with p-ELK1-positive non-muscle-invasive and muscle-invasive tumors had significantly higher risks for tumor recurrence and progression, respectively. In AR-positive bladder cancer cell lines, dihydrotestosterone treatment increased ELK1 expression (mRNA, protein) and its nuclear translocation, ELK1 transcriptional activity, and c-fos expression, which was restored by an anti-androgen hydroxyflutamide. ELK1 silencing via short hairpin RNA (shRNA) resulted in decreases in cell viability/colony formation, and cell migration/invasion as well as an increase in apoptosis. Importantly, ELK1 appears to require activated AR to regulate bladder cancer cell proliferation, but not cell migration. Androgen also failed to significantly induce AR transactivation in ELK1-knockdown cells. In accordance with our in vitro findings, ELK1-shRNA expression considerably retarded tumor formation as well as its growth in xenograft-bearing male mice. Our results suggest that ELK1 plays an important role in bladder tumorigenesis and cancer progression, which is further induced by AR activation. Accordingly, ELK1 inhibition, together with AR inactivation, has the potential of being a therapeutic approach for bladder cancer.

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.

Quantifying signaling pathway activation to monitor the quality of induced pluripotent stem cells

Many attempts have been made to evaluate the safety and potency of human induced pluripotent stem cells (iPSCs) for clinical applications using transcriptome data, but results so far have been ambiguous or even contradictory. Here, we characterized stem cells at the pathway level, rather than at the gene level as has been the focus of previous work. We meta-analyzed publically-available gene expression data sets and evaluated signaling and metabolic pathway activation profiles for 20 human embryonic stem cell (ESC) lines, 12 human iPSC lines, five embryonic body lines, and six fibroblast cell lines. We demonstrated the close resemblance of iPSCs with ESCs at the pathway level, and provided examples of how pathway activity can be applied to identify iPSC line abnormalities or to predict in vitro differentiation potential. Our results indicate that pathway activation profiling is a promising strategy for evaluating the safety and potency of iPSC lines in translational medicine applications.

Common pathway signature in lung and liver fibrosis

Fibrosis, a progressive accumulation of extracellular matrix components, encompasses a wide spectrum of distinct organs, and accounts for an increasing burden of morbidity and mortality worldwide. Despite the tremendous clinical impact, the mechanisms governing the fibrotic process are not yet understood, and to date, no clinically reliable therapies for fibrosis have been discovered. Here we applied Regeneration Intelligence, a new bioinformatics software suite for qualitative analysis of intracellular signaling pathway activation using transcriptomic data, to assess a network of molecular signaling in lung and liver fibrosis. In both tissues, our analysis detected major conserved signaling pathways strongly associated with fibrosis, suggesting that some of the pathways identified by our algorithm but not yet wet-lab validated as fibrogenesis related, may be attractive targets for future research. While the majority of significantly disrupted pathways were specific to histologically distinct organs, several pathways have been concurrently activated or downregulated among the hepatic and pulmonary fibrosis samples, providing new evidence of evolutionary conserved pathways that may be relevant as possible therapeutic targets. While future confirmatory studies are warranted to validate these observations, our platform proposes a promising new approach for detecting fibrosis-promoting pathways and tailoring the right therapy to prevent fibrogenesis.

Applications of Deep Learning in Biomedicine

Increases in throughput and installed base of biomedical research equipment led to a massive accumulation of -omics data known to be highly variable, high-dimensional, and sourced from multiple often incompatible data platforms. While this data may be useful for biomarker identification and drug discovery, the bulk of it remains underutilized. Deep neural networks (DNNs) are efficient algorithms based on the use of compositional layers of neurons, with advantages well matched to the challenges -omics data presents. While achieving state-of-the-art results and even surpassing human accuracy in many challenging tasks, the adoption of deep learning in biomedicine has been comparatively slow. Here, we discuss key features of deep learning that may give this approach an edge over other machine learning methods. We then consider limitations and review a number of applications of deep learning in biomedical studies demonstrating proof of concept and practical utility.

Pathway activation profiling reveals new insights into age-related macular degeneration and provides avenues for therapeutic interventions

Age-related macular degeneration (AMD) is a major cause of blindness in older people and is caused by loss of the central region of the retinal pigment epithelium (RPE). Conventional methods of gene expression analysis have yielded important insights into AMD pathogenesis, but the precise molecular pathway alterations are still poorly understood. Therefore we developed a new software program, “AMD Medicine”, and discovered differential pathway activation profiles in samples of human RPE/choroid from AMD patients and controls. We identified 29 pathways in RPE-choroid AMD phenotypes: 27 pathways were activated in AMD compared to controls, and 2 pathways were activated in controls compared to AMD. In AMD, we identified a graded activation of pathways related to wound response, complement cascade, and cell survival. Also, there was downregulation of two pathways responsible for apoptosis. Furthermore, significant activation of pro-mitotic pathways is consistent with dedifferentiation and cell proliferation events, which occur early in the pathogenesis of AMD. Significantly, we discovered new global pathway activation signatures of AMD involved in the cell-based inflammatory response: IL-2, STAT3, and ERK. The ultimate aim of our research is to achieve a better understanding of signaling pathways involved in AMD pathology, which will eventually lead to better treatments.

Pathway activation strength is a novel independent prognostic biomarker for cetuximab sensitivity in colorectal cancer patients

Cetuximab, a monoclonal antibody against epidermal growth factor receptor (EGFR), was shown to be active in colorectal cancer. Although some patients who harbor K-ras wild-type tumors benefit from cetuximab treatment, 40 to 60% of patients with wild-type K-ras tumors do not respond to cetuximab. Currently, there is no universal marker or method of clinical utility that could guide the treatment of cetuximab in colorectal cancer. Here, we demonstrate a method to predict response to cetuximab in patients with colorectal cancer using OncoFinder pathway activation strength (PAS), based on the transcriptomic data of the tumors. We first evaluated our OncoFinder pathway activation strength model in a set of transcriptomic data obtained from patient-derived xenograft (PDx) models established from colorectal cancer biopsies. Then, the approach and models were validated using a clinical trial data set. PAS could efficiently predict patients’ response to cetuximab, and thus holds promise as a selection criterion for cetuximab treatment in metastatic colorectal cancer.

Altered radiation responses of breast cancer cells resistant to hormonal therapy

Endocrine therapy agents (the selective estrogen receptor (ER) modulators such as tamoxifen or the selective ER down-regulators such as ICI 182,780) are key treatment regimens for hormone receptor-positive breast cancers. While these drugs are very effective in controlling ER-positive breast cancer, many tumors that initially respond well to treatment often acquire drug resistance, which is a major clinical problem. In clinical practice, hormonal therapy agents are commonly used in combination or sequence with radiation therapy. Tamoxifen treatment and radiotherapy improve both local tumor control and patient survival. However, tamoxifen treatment may render cancer cells less responsive to radiation therapy. Only a handful of data exist on the effects of radiation on cells resistant to hormonal therapy agents. These scarce data show that cells that were resistant to tamoxifen were also resistant to radiation. Yet, the existence and mechanisms of cross-resistance to endocrine therapy and radiation therapy need to be established. Here, we for the first time examined and compared radiation responses of MCF-7 breast adenocarcinoma cells (MCF-7/S0.5) and two antiestrogen resistant cell lines derived from MCF-7/S0.5: the tamoxifen resistant MCF-7/TAMR-1 and ICI 182,780 resistant MCF-7/182R-6 cell lines. Specifically, we analyzed the radiation-induced changes in the expression of genes involved in DNA damage, apoptosis, and cell cycle regulation. We found that the tamoxifen-resistant cell line in contrast to the parental and ICI 182,780-resistant cell lines displayed a significantly less radiation-induced decrease in the expression of genes involved in DNA repair. Furthermore, we show that MCF-7/TAMR-1 and MCF-7/182R-6 cells were less susceptible to radiation-induced apoptosis as compared to the parental line. These data indicate that tamoxifen-resistant breast cancer cells have a reduced sensitivity to radiation treatment. The current study may therefore serve as a roadmap to the future analysis of the mechanisms of cross-resistance between hormonal therapy and radiation.