New treatment alternatives for primary and metastatic colorectal cancer by an integrated transcriptome and network analyses

Metastatic colorectal cancer (CRC) is still in need of effective treatments. This study applies a holistic approach to propose new targets for treatment of primary and liver metastatic CRC and investigates their therapeutic potential in-vitro. An integrative analysis of primary and metastatic CRC samples was implemented for alternative target and treatment proposals. Integrated microarray samples were grouped based on a co-expression network analysis. Significant gene modules correlated with primary CRC and metastatic phenotypes were identified. Network clustering and pathway enrichments were applied to gene modules to prioritize potential targets, which were shortlisted by independent validation. Finally, drug-target interaction search led to three agents for primary and liver metastatic CRC phenotypes. Hesperadin and BAY-1217389 suppress colony formation over a 14-day period, with Hesperadin showing additional efficacy in reducing cell viability within 48 h. As both candidates target the G2/M phase proteins NEK2 or TTK, we confirmed their anti-proliferative properties by Ki-67 staining. Hesperadinin particular arrested the cell cycle at the G2/M phase. IL-29A treatment reduced migration and invasion capacities of TGF-β induced metastatic cell lines. In addition, this anti-metastatic treatment attenuated TGF-β dependent mesenchymal transition. Network analysis suggests IL-29A induces the JAK/STAT pathway in a preventive manner.

In Silico Approach to Molecular Profiling of the Transition from Ovarian Epithelial Cells to Low-Grade Serous Ovarian Tumors for Targeted Therapeutic Insights

This paper aims to elucidate the differentially coexpressed genes, their potential mechanisms, and possible drug targets in low-grade invasive serous ovarian carcinoma (LGSC) in terms of the biologic continuity of normal, borderline, and malignant LGSC. We performed a bioinformatics analysis, integrating datasets generated using the GPL570 platform from different studies from the GEO database to identify changes in this transition, gene expression, drug targets, and their relationships with tumor microenvironmental characteristics. In the transition from ovarian epithelial cells to the serous borderline, the FGFR3 gene in the "Estrogen Response Late" pathway, the ITGB2 gene in the "Cell Adhesion Molecule", the CD74 gene in the "Regulation of Cell Migration", and the IGF1 gene in the "Xenobiotic Metabolism" pathway were upregulated in the transition from borderline to LGSC. The ERBB4 gene in "Proteoglycan in Cancer", the AR gene in "Pathways in Cancer" and "Estrogen Response Early" pathways, were upregulated in the transition from ovarian epithelial cells to LGSC. In addition, SPP1 and ITGB2 genes were correlated with macrophage infiltration in the LGSC group. This research provides a valuable framework for the development of personalized therapeutic approaches in the context of LGSC, with the aim of improving patient outcomes and quality of life. Furthermore, the main goal of the current study is a preliminary study designed to generate in silico inferences, and it is also important to note that subsequent in vitro and in vivo studies will be necessary to confirm the results before considering these results as fully reliable.

Network neighborhood operates as a drug repositioning method for cancer treatment

Computational drug repositioning approaches are important, as they cost less compared to the traditional drug development processes. This study proposes a novel network-based drug repositioning approach, which computes similarities between disease-causing genes and drug-affected genes in a network topology to suggest candidate drugs with highest similarity scores. This new method aims to identify better treatment options by integrating systems biology approaches. It uses a protein-protein interaction network that is the main topology to compute a similarity score between candidate drugs and disease-causing genes. The disease-causing genes were mapped on this network structure. Transcriptome profiles of drug candidates were taken from the LINCS project and mapped individually on the network structure. The similarity of these two networks was calculated by different network neighborhood metrics, including Adamic-Adar, PageRank and neighborhood scoring. The proposed approach identifies the best candidates by choosing the drugs with significant similarity scores. The method was experimented on melanoma, colorectal, and prostate cancers. Several candidate drugs were predicted by applying AUC values of 0.6 or higher. Some of the predictions were approved by clinical phase trials or other in-vivo studies found in literature. The proposed drug repositioning approach would suggest better treatment options with integration of functional information between genes and transcriptome level effects of drug perturbations and diseases.

In silico identification of novel biomarkers for key players in transition from normal colon tissue to adenomatous polyps

Adenomatous polyps of the colon are the most common neoplastic polyps. Although most of adenomatous polyps do not show malign transformation, majority of colorectal carcinomas originate from neoplastic polyps. Therefore, understanding of this transformation process would help in both preventive therapies and evaluation of malignancy risks. This study uncovers alterations in gene expressions as potential biomarkers that are revealed by integration of several network-based approaches. In silico analysis performed on a unified microarray cohort, which is covering 150 normal colon and adenomatous polyp samples. Significant gene modules were obtained by a weighted gene co-expression network analysis. Gene modules with similar profiles were mapped to a colon tissue specific functional interaction network. Several clustering algorithms run on the colon-specific network and the most significant sub-modules between the clusters were identified. The biomarkers were selected by filtering differentially expressed genes which also involve in significant biological processes and pathways. Biomarkers were also validated on two independent datasets based on their differential gene expressions. To the best of our knowledge, such a cascaded network analysis pipeline was implemented for the first time on a large collection of normal colon and polyp samples. We identified significant increases in TLR4 and MSX1 expressions as well as decrease in chemokine profiles with mostly pro-tumoral activities. These biomarkers might appear as both preventive targets and biomarkers for risk evaluation. As a result, this research proposes novel molecular markers that might be alternative to endoscopic approaches for diagnosis of adenomatous polyps.

Integrative Biological Network Analysis to Identify Shared Genes in Metabolic Disorders

Identification of common molecular mechanisms in interrelated diseases is essential for better prognoses and targeted therapies. However, complexity of metabolic pathways makes it difficult to discover common disease genes underlying metabolic disorders; and it requires more sophisticated bioinformatics models that combine different types of biological data and computational methods. Accordingly, we built an integrative network analysis model to identify shared disease genes in metabolic syndrome (MS), type 2 diabetes (T2D), and coronary artery disease (CAD). We constructed weighted gene co-expression networks by combining gene expression, protein-protein interaction, and gene ontology data from multiple sources. For 90 different configurations of disease networks, we detected the significant modules by using MCL, SPICi, and Linkcomm graph clustering algorithms. We also performed a comparative evaluation on disease modules to determine the best method providing the highest biological validity. By overlapping the disease modules, we identified 22 shared genes for MS-CAD and T2D-CAD. Moreover, 19 out of these genes were directly or indirectly associated with relevant diseases in the previous medical studies. This study does not only demonstrate the performance of different biological data sources and computational methods in disease-gene discovery, but also offers potential insights into common genetic mechanisms of the metabolic disorders.

Predicting drug synergy for precision medicine using network biology and machine learning

Identification of effective drug combinations for patients is an expensive and time-consuming procedure, especially for in vitro experiments. To accelerate the synergistic drug discovery process, we present a new classification model to identify more effective anti-cancer drug pairs using in silico network biology approach. Based on the hypotheses that the drug synergy comes from the collective effects on the biological network, therefore, we developed six network biology features, including overlap and distance of drug perturbation network, that were derived by using individual drug-perturbed transcriptome profiles and the relevant biological network analysis. Using publicly available drug synergy databases and three machine-learning (ML) methods, the model was trained to discriminate the positive (synergistic) and negative (nonsynergistic) drug combinations. The proposed models were evaluated on the test cases to predict the most promising network biology feature, which is the network degree activity, i.e. the synergistic effect between drug pairs is mainly accounted by the complementary signaling pathways or molecular networks from two drugs.

Ubiquitin-Proteasome Axis, Especially Ubiquitin-Specific Protease-17 (USP17) Gene Family, is a Potential Target for Epithelial-Mesenchymal Transition in High-Grade Serous Ovarian Cancer

OBJECTIVES

To investigate gene expression differences and related functions between primary tumor, malignant cells in ascites, and metastatic peritoneal implant in high-grade serous ovarian cancer.

METHODS

Biopsies from primary tumor, peritoneal implant, and ascites were collected from 10 patients operated primarily for high-grade, advanced-staged serous ovarian cancer. Total RNA isolation was performed from collected tissue biopsy and fluid samples, and RNA expression profile was measured. Messenger RNA expression profiles of 3 different groups were compared. Functional analyses of candidate genes were carried out by gene ontology and pathway analysis.

RESULTS

There were significant differences in the expression of 5 genes between primary tumor and peritoneal implant, 979 genes between primary tumor and malignant cells in ascites, and 649 genes between peritoneal implant and malignant cells in ascites. Three commonly enriched gene ontology functions between "primary tumor and malignant cells in the ascites" and "peritoneal implant and malignant cells in the ascites" were protein deubiquitination, ubiquitin-dependent protein catabolism, and apoptotic processes. All genes related to these functions belonged to USP17 gene family.

CONCLUSION

Gene expression difference between primary tumor and the peritoneal implant is not as much as the difference between primary tumor and free cells in the ascites. These results show that malignant cells in the ascites return into its genetic origin after they invade on the peritoneum. Significantly increased expression of DUB-enzyme genes, SNAR gene family, and ribosomal pathway genes in epithelial-mesenchymal transition suggests that this regulation is ubiquitin-proteasome dependent. Especially, this is the first study that offers USP17 as a potential target for epithelial-mesenchymal transition.

Integration of RNA-Seq and RPPA data for survival time prediction in cancer patients

Integration of several types of patient data in a computational framework can accelerate the identification of more reliable biomarkers, especially for prognostic purposes. This study aims to identify biomarkers that can successfully predict the potential survival time of a cancer patient by integrating the transcriptomic (RNA-Seq), proteomic (RPPA), and protein-protein interaction (PPI) data. The proposed method -RPBioNet- employs a random walk-based algorithm that works on a PPI network to identify a limited number of protein biomarkers. Later, the method uses gene expression measurements of the selected biomarkers to train a classifier for the survival time prediction of patients. RPBioNet was applied to classify kidney renal clear cell carcinoma (KIRC), glioblastoma multiforme (GBM), and lung squamous cell carcinoma (LUSC) patients based on their survival time classes (long- or short-term). The RPBioNet method correctly identified the survival time classes of patients with between 66% and 78% average accuracy for three data sets. RPBioNet operates with only 20 to 50 biomarkers and can achieve on average 6% higher accuracy compared to the closest alternative method, which uses only RNA-Seq data in the biomarker selection. Further analysis of the most predictive biomarkers highlighted genes that are common for both cancer types, as they may be driver proteins responsible for cancer progression. The novelty of this study is the integration of a PPI network with mRNA and protein expression data to identify more accurate prognostic biomarkers that can be used for clinical purposes in the future.

Systematic drug perturbations on cancer cells reveal diverse exit paths from proliferative state

During a cell state transition, cells travel along trajectories in a gene expression state space. This dynamical systems framework complements the traditional concept of molecular pathways that drive cell phenotype switching. To expose the structure that hinders cancer cells from exiting robust proliferative state, we assessed the perturbation capacity of a drug library and identified 16 non-cytotoxic compounds that stimulate MCF7 breast cancer cells to exit from proliferative state to differentiated state. The transcriptome trajectories triggered by these drugs diverged, then converged. Chemical structures and drug targets of these compounds overlapped minimally. However, a network analysis of targeted pathways identified a core signaling pathway - indicating common stress-response and down-regulation of STAT1 before differentiation. This multi-trajectory analysis explores the cells' state transition with a multitude of perturbations in combination with traditional pathway analysis, leading to an encompassing picture of the dynamics of a therapeutically desired cell-state switching.

A signal transduction score flow algorithm for cyclic cellular pathway analysis, which combines transcriptome and ChIP-seq data

Determination of cell signalling behaviour is crucial for understanding the physiological response to a specific stimulus or drug treatment. Current approaches for large-scale data analysis do not effectively incorporate critical topological information provided by the signalling network. We herein describe a novel model- and data-driven hybrid approach, or signal transduction score flow algorithm, which allows quantitative visualization of cyclic cell signalling pathways that lead to ultimate cell responses such as survival, migration or death. This score flow algorithm translates signalling pathways as a directed graph and maps experimental data, including negative and positive feedbacks, onto gene nodes as scores, which then computationally traverse the signalling pathway until a pre-defined biological target response is attained. Initially, experimental data-driven enrichment scores of the genes were computed in a pathway, then a heuristic approach was applied using the gene score partition as a solution for protein node stoichiometry during dynamic scoring of the pathway of interest. Incorporation of a score partition during the signal flow and cyclic feedback loops in the signalling pathway significantly improves the usefulness of this model, as compared to other approaches. Evaluation of the score flow algorithm using both transcriptome and ChIP-seq data-generated signalling pathways showed good correlation with expected cellular behaviour on both KEGG and manually generated pathways. Implementation of the algorithm as a Cytoscape plug-in allows interactive visualization and analysis of KEGG pathways as well as user-generated and curated Cytoscape pathways. Moreover, the algorithm accurately predicts gene-level and global impacts of single or multiple in silico gene knockouts.

Network information improves cancer outcome prediction

Disease progression in cancer can vary substantially between patients. Yet, patients often receive the same treatment. Recently, there has been much work on predicting disease progression and patient outcome variables from gene expression in order to personalize treatment options. Despite first diagnostic kits in the market, there are open problems such as the choice of random gene signatures or noisy expression data. One approach to deal with these two problems employs protein-protein interaction networks and ranks genes using the random surfer model of Google's PageRank algorithm. In this work, we created a benchmark dataset collection comprising 25 cancer outcome prediction datasets from literature and systematically evaluated the use of networks and a PageRank derivative, NetRank, for signature identification. We show that the NetRank performs significantly better than classical methods such as fold change or t-test. Despite an order of magnitude difference in network size, a regulatory and protein-protein interaction network perform equally well. Experimental evaluation on cancer outcome prediction in all of the 25 underlying datasets suggests that the network-based methodology identifies highly overlapping signatures over all cancer types, in contrast to classical methods that fail to identify highly common gene sets across the same cancer types. Integration of network information into gene expression analysis allows the identification of more reliable and accurate biomarkers and provides a deeper understanding of processes occurring in cancer development and progression.

Organellar genome analysis of rye (Secale cereale) representing diverse geographic regions

Rye (Secale cereale) is an important diploid (2n = 14, RR) crop species of the Triticeae and a better understanding of its organellar genome variation can aid in its improvement. Previous genetic analyses of rye focused on the nuclear genome. In the present study, the objective was to investigate the organellar genome diversity and relationships of 96 accessions representing diverse geographic regions using chloroplast (cp) and mitochondrial (mt) DNA PCR-RFLPs. Seven cpDNA and 4 mtDNA coding and noncoding regions were amplified using universal cpDNA and mtDNA primer pairs. Each amplified fragment was digested with 13 different restriction enzymes. mtDNA analysis indicated that the number of polymorphic loci (20) was low and genetic differentiation (GST) was 0.60, excluding the outgroups (hexaploid wheat, Triticum aestivum, 2n = 6x = 42, AABBDD; triticale, xTriticosecale Wittmack, 2n = 6x = 42, AABBRR). cpDNA analysis revealed a low level of polymorphism (40%) among the accessions, and GST was 0.39. Of the 96 genotypes studied, 70 could not be differentiated using cpDNA PCR-RFLPs even though they are from different geographic regions. This is most likely due to germplasm exchange, indicating that genotypes might have a common genetic background. Two cpDNA and 3 mtDNA fragments were significantly correlated to the site of germplasm collection. However, there was no clear trend. These results indicate that the level of organellar polymorphism is low among the cultivated rye genotypes. The cpDNA and mtDNA PCR-RFLP markers used in the present study could be used as molecular markers in rye genetics and breeding programs.