Cancer modeling and network biology: Accelerating toward personalized medicine

Innovations in Genomics and Big Data Analytics for Personalized Medicine and Health Care: A Review

Big data in health care is a fast-growing field and a new paradigm that is transforming case-based studies to large-scale, data-driven research. As big data is dependent on the advancement of new data standards, technology, and relevant research, the future development of big data applications holds foreseeable promise in the modern day health care revolution. Enormously large, rapidly growing collections of biomedical omics-data (genomics, proteomics, transcriptomics, metabolomics, glycomics, etc.) and clinical data create major challenges and opportunities for their analysis and interpretation and open new computational gateways to address these issues. The design of new robust algorithms that are most suitable to properly analyze this big data by taking into account individual variability in genes has enabled the creation of precision (personalized) medicine. We reviewed and highlighted the significance of big data analytics for personalized medicine and health care by focusing mostly on machine learning perspectives on personalized medicine, genomic data models with respect to personalized medicine, the application of data mining algorithms for personalized medicine as well as the challenges we are facing right now in big data analytics.

Molecular differential analysis of uterine leiomyomas and leiomyosarcomas through weighted gene network and pathway tracing approaches

Uterine smooth muscular neoplastic growths like benign leiomyomas (UL) and metastatic leiomyosarcomas (ULMS) share similar clinical symptoms, radiological and histological appearances making their clinical distinction a difficult task. Therefore, the objective of this study is to identify key genes and pathways involved in transformation of UL to ULMS through molecular differential analysis. Global gene expression profiles of 25 ULMS, 25 UL, and 29 myometrium (Myo) tissues generated on Affymetrix U133A 2.0 human genome microarrays were analyzed by deploying robust statistical, molecular interaction network, and pathway enrichment methods. The comparison of expression signals across Myo vs UL, Myo vs ULMS, and UL vs ULMS groups identified 249, 1037, and 716 significantly expressed genes, respectively (p ≤ 0.05). The analysis of 249 DEGs from Myo vs UL confirms multistage dysregulation of various key pathways in extracellular matrix, collagen, cell contact inhibition, and cytokine receptors transform normal myometrial cells to benign leiomyomas (p value ≤ 0.01). The 716 DEGs between UL vs ULMS were found to affect cell cycle, cell division related Rho GTPases and PI3K signaling pathways triggering uncontrolled growth and metastasis of tumor cells (p value ≤ 0.01). Integration of gene networking data, with additional parameters like estimation of mutation burden of tumors and cancer driver gene identification, has led to the finding of 4 hubs (JUN, VCAN, TOP2A, and COL1A1) and 8 bottleneck genes (PIK3R1, MYH11, KDR, ESR1, WT1, CCND1, EZH2, and CDKN2A), which showed a clear distinction in their distribution pattern among leiomyomas and leiomyosarcomas. This study provides vital clues for molecular distinction of UL and ULMS which could further assist in identification of specific diagnostic markers and therapeutic targets.Abbreviations UL: Uterine Leiomyomas; ULMS: Uterine Leiomyosarcoma; Myo: Myometrium; DEGs: Differential Expressed Genes; RMA: Robust Multiarray Average; DC: Degree of Centrality; BC: Betweenness of Centrality; CGC: Cancer Gene Census; FDR: False Discovery Rate; TCGA: Cancer Genome Atlas; BP: Biological Process; CC: Cellular Components; MF: Molecular Function; PPI: Protein-Protein Interaction.

In silico mouse study identifies tumour growth kinetics as biomarkers for the outcome of anti-angiogenic treatment

Angiogenesis is a crucial step in tumour progression, as this process allows tumours to recruit new blood vessels and obtain oxygen and nutrients to sustain growth. Therefore, inhibiting angiogenesis remains a viable strategy for cancer therapy. However, anti-angiogenic therapy has not proved to be effective in reducing tumour growth across a wide range of tumours, and no reliable predictive biomarkers have been found to determine the efficacy of anti-angiogenic treatment. Using our previously established computational model of tumour-bearing mice, we sought to determine whether tumour growth kinetic parameters could be used to predict the outcome of anti-angiogenic treatment. A model trained with datasets from six in vivo mice studies was used to generate a randomized in silico tumour-bearing mouse population. We analysed tumour growth in untreated mice (control) and mice treated with an anti-angiogenic agent and determined the Kaplan–Meier survival estimates based on simulated tumour volume data. We found that the ratio between two kinetic parameters, k₀ and k₁, which characterize the tumour's exponential and linear growth rates, as well as k₁ alone, can be used as prognostic biomarkers of the population survival outcome. Our work demonstrates a robust, quantitative approach for identifying tumour growth kinetic parameters as prognostic biomarkers and serves as a template that can be used to identify other biomarkers for anti-angiogenic treatment.

Analyzing of Molecular Networks for Human Diseases and Drug Discovery

Molecular networks represent the interactions and relations of genes/proteins, and also encode molecular mechanisms of biological processes, development and diseases. Among the molecular networks, protein-protein Interaction Networks (PINs) have become effective platforms for uncovering the molecular mechanisms of diseases and drug discovery. PINs have been constructed for various organisms and utilized to solve many biological problems. In human, most proteins present their complex functions by interactions with other proteins, and the sum of these interactions represents the human protein interactome. Especially in the research on human disease and drugs, as an emerging tool, the PIN provides a platform to systematically explore the molecular complexities of specific diseases and the references for drug design. In this review, we summarized the commonly used approaches to aid disease research and drug discovery with PINs, including the network topological analysis, identification of novel pathways, drug targets and sub-network biomarkers for diseases. With the development of bioinformatic techniques and biological networks, PINs will play an increasingly important role in human disease research and drug discovery.

IRWRLDA: improved random walk with restart for lncRNA-disease association prediction

In recent years, accumulating evidences have shown that the dysregulations of lncRNAs are associated with a wide range of human diseases. It is necessary and feasible to analyze known lncRNA-disease associations, predict potential lncRNA-disease associations, and provide the most possible lncRNA-disease pairs for experimental validation. Considering the limitations of traditional Random Walk with Restart (RWR), the model of Improved Random Walk with Restart for LncRNA-Disease Association prediction (IRWRLDA) was developed to predict novel lncRNA-disease associations by integrating known lncRNA-disease associations, disease semantic similarity, and various lncRNA similarity measures. The novelty of IRWRLDA lies in the incorporation of lncRNA expression similarity and disease semantic similarity to set the initial probability vector of the RWR. Therefore, IRWRLDA could be applied to diseases without any known related lncRNAs. IRWRLDA significantly improved previous classical models with reliable AUCs of 0.7242 and 0.7872 in two known lncRNA-disease association datasets downloaded from the lncRNADisease database, respectively. Further case studies of colon cancer and leukemia were implemented for IRWRLDA and 60% of lncRNAs in the top 10 prediction lists have been confirmed by recent experimental reports.

FMLNCSIM: fuzzy measure-based lncRNA functional similarity calculation model

Accumulating experimental studies have indicated the influence of lncRNAs on various critical biological processes as well as disease development and progression. Calculating lncRNA functional similarity is of high value in inferring lncRNA functions and identifying potential lncRNA-disease associations. However, little effort has been attempt to measure the functional similarity among lncRNAs on a large scale. In this study, we developed a Fuzzy Measure-based LNCRNA functional SIMilarity calculation model (FMLNCSIM) based on the assumption that functionally similar lncRNAs tend to be associated with similar diseases. The performance improvement of FMLNCSIM mainly comes from the combination of information content and the concept of fuzzy measure, which was applied to the directed acyclic graphs of disease MeSH descriptors. To evaluate the effectiveness of FMLNCSIM, we further combined it with the previously proposed model of Laplacian Regularized Least Squares for lncRNA-Disease Association (LRLSLDA). As a result, the integrated model, LRLSLDA-FMLNCSIM, achieve good performance in the frameworks of global LOOCV (AUCs of 0.8266 and 0.9338 based on LncRNADisease and MNDR database) and 5-fold cross validation (average AUCs of 0.7979 and 0.9237 based on LncRNADisease and MNDR database), which significantly improve the performance of previous classical models. It is anticipated that FMLNCSIM could be used for searching functionally similar lncRNAs and inferring lncRNA functions in the future researches.

Mechanistic modeling quantifies the influence of tumor growth kinetics on the response to anti-angiogenic treatment

Tumors exploit angiogenesis, the formation of new blood vessels from pre-existing vasculature, in order to obtain nutrients required for continued growth and proliferation. Targeting factors that regulate angiogenesis, including the potent promoter vascular endothelial growth factor (VEGF), is therefore an attractive strategy for inhibiting tumor growth. Computational modeling can be used to identify tumor-specific properties that influence the response to anti-angiogenic strategies. Here, we build on our previous systems biology model of VEGF transport and kinetics in tumor-bearing mice to include a tumor compartment whose volume depends on the “angiogenic signal” produced when VEGF binds to its receptors on tumor endothelial cells. We trained and validated the model using published in vivo measurements of xenograft tumor volume, producing a model that accurately predicts the tumor’s response to anti-angiogenic treatment. We applied the model to investigate how tumor growth kinetics influence the response to anti-angiogenic treatment targeting VEGF. Based on multivariate regression analysis, we found that certain intrinsic kinetic parameters that characterize the growth of tumors could successfully predict response to anti-VEGF treatment, the reduction in tumor volume. Lastly, we use the trained model to predict the response to anti-VEGF therapy for tumors expressing different levels of VEGF receptors. The model predicts that certain tumors are more sensitive to treatment than others, and the response to treatment shows a nonlinear dependence on the VEGF receptor expression. Overall, this model is a useful tool for predicting how tumors will respond to anti-VEGF treatment, and it complements pre-clinical in vivo mouse studies.

One hallmark of cancer is angiogenesis, the formation of new blood capillaries from pre-existing vessels. Angiogenesis promotes tumor growth by enabling the tumor to obtain oxygen and nutrients from the surrounding microenvironment. Cancer drugs that inhibit angiogenesis ("anti-angiogenic therapies") have focused on inhibiting proteins that promote the growth of new blood vessels. The response to anti-angiogenic therapy is highly variable, and some tumors do not respond at all. Therefore, identifying a biomarker that predicts how specific tumors will respond would be extremely valuable. This work uses a computational model of tumor-bearing mice to investigate the response to anti-angiogenic treatment that targets the potent promoter of angiogenesis, vascular endothelial growth factor (VEGF), and how the response is influenced by tumor growth kinetics. We show that certain properties of tumor growth can be used to predict how much the tumor volume will be reduced upon administration of an anti-VEGF drug. This work identifies tumor growth parameters that may be reliable biomarkers for predicting how tumors will respond to anti-VEGF therapy. Our computational model generates novel, testable hypotheses and nicely complements pre-clinical studies of anti-angiogenic therapeutics.

Individualized network-based drug repositioning infrastructure for precision oncology in the panomics era

Advances in next-generation sequencing technologies have generated the data supporting a large volume of somatic alterations in several national and international cancer genome projects, such as The Cancer Genome Atlas and the International Cancer Genome Consortium. These cancer genomics data have facilitated the revolution of a novel oncology drug discovery paradigm from candidate target or gene studies toward targeting clinically relevant driver mutations or molecular features for precision cancer therapy. This focuses on identifying the most appropriately targeted therapy to an individual patient harboring a particularly genetic profile or molecular feature. However, traditional experimental approaches that are used to develop new chemical entities for targeting the clinically relevant driver mutations are costly and high-risk. Drug repositioning, also known as drug repurposing, re-tasking or re-profiling, has been demonstrated as a promising strategy for drug discovery and development. Recently, computational techniques and methods have been proposed for oncology drug repositioning and identifying pharmacogenomics biomarkers, but overall progress remains to be seen. In this review, we focus on introducing new developments and advances of the individualized network-based drug repositioning approaches by targeting the clinically relevant driver events or molecular features derived from cancer panomics data for the development of precision oncology drug therapies (e.g. one-person trials) to fully realize the promise of precision medicine. We discuss several potential challenges (e.g. tumor heterogeneity and cancer subclones) for precision oncology. Finally, we highlight several new directions for the precision oncology drug discovery via biotherapies (e.g. gene therapy and immunotherapy) that target the 'undruggable' cancer genome in the functional genomics era.

Inferring clonal structure in HTLV-1-infected individuals: towards bridging the gap between analysis and visualization

Background

Human T cell leukemia virus type 1 (HTLV-1) causes adult T cell leukemia (ATL) in a proportion of infected individuals after a long latency period. Development of ATL is a multistep clonal process that can be investigated by monitoring the clonal expansion of HTLV-1-infected cells by isolation of provirus integration sites. The clonal composition (size, number, and combinations of clones) during the latency period in a given infected individual has not been clearly elucidated.

Methods

We used high-throughput sequencing technology coupled with a tag system for isolating integration sites and measuring clone sizes from 60 clinical samples. We assessed the role of clonality and clone size dynamics in ATL onset by modeling data from high-throughput monitoring of HTLV-1 integration sites using single- and multiple-time-point samples.

Results

From four size categories analyzed, we found that big clones (B; 513–2048 infected cells) and very big clones (VB; >2048 infected cells) had prognostic value. No sample harbored two or more VB clones or three or more B clones. We examined the role of clone size, clone combination, and the number of integration sites in the prognosis of infected individuals. We found a moderate reverse correlation between the total number of clones and the size of the largest clone. We devised a data-driven model that allows intuitive representation of clonal composition.

Conclusions

This integration site-based clonality tree model represents the complexity of clonality and provides a global view of clonality data that facilitates the analysis, interpretation, understanding, and visualization of the behavior of clones on inter- and intra-individual scales. It is fully data-driven, intuitively depicts the clonality patterns of HTLV-1-infected individuals and can assist in early risk assessment of ATL onset by reflecting the prognosis of infected individuals. This model should assist in assimilating information on clonal composition and understanding clonal expansion in HTLV-1-infected individuals.

Electronic supplementary material

The online version of this article (doi:10.1186/s40246-017-0112-8) contains supplementary material, which is available to authorized users.

A Generic Individual-Based Spatially Explicit Model as a Novel Tool for Investigating Insect-Plant Interactions: A Case Study of the Behavioural Ecology of Frugivorous Tephritidae

Computational modelling of mechanisms underlying processes in the real world can be of great value in understanding complex biological behaviours. Uptake in general biology and ecology has been rapid. However, it often requires specific data sets that are overly costly in time and resources to collect. The aim of the current study was to test whether a generic behavioural ecology model constructed using published data could give realistic outputs for individual species. An individual-based model was developed using the Pattern-Oriented Modelling (POM) strategy and protocol, based on behavioural rules associated with insect movement choices. Frugivorous Tephritidae (fruit flies) were chosen because of economic significance in global agriculture and the multiple published data sets available for a range of species. The Queensland fruit fly (Qfly), Bactrocera tryoni, was identified as a suitable individual species for testing. Plant canopies with modified architecture were used to run predictive simulations. A field study was then conducted to validate our model predictions on how plant architecture affects fruit flies’ behaviours. Characteristics of plant architecture such as different shapes, e.g., closed-canopy and vase-shaped, affected fly movement patterns and time spent on host fruit. The number of visits to host fruit also differed between the edge and centre in closed-canopy plants. Compared to plant architecture, host fruit has less contribution to effects on flies’ movement patterns. The results from this model, combined with our field study and published empirical data suggest that placing fly traps in the upper canopy at the edge should work best. Such a modelling approach allows rapid testing of ideas about organismal interactions with environmental substrates in silico rather than in vivo, to generate new perspectives. Using published data provides a saving in time and resources. Adjustments for specific questions can be achieved by refinement of parameters based on targeted experiments.

Analysis of the intricate relationship between chronic inflammation and cancer

Deregulated inflammatory response plays a pivotal role in the initiation, development and progression of tumours. Potential molecular mechanism(s) that drive the establishment of an inflammatory-tumour microenvironment is not entirely understood owing to the complex cross-talk between pro-inflammatory and tumorigenic mediators such as cytokines, chemokines, oncogenes, enzymes, transcription factors and immune cells. These molecular mediators are critical linchpins between inflammation and cancer, and their activation and/or deactivation are influenced by both extrinsic (i.e. environmental and lifestyle) and intrinsic (i.e. hereditary) factors. At present, the research pertaining to inflammation-associated cancers is accumulating at an exponential rate. Interest stems from hope that new therapeutic strategies against molecular mediators can be identified to assist in cancer treatment and patient management. The present review outlines the various molecular and cellular inflammatory mediators responsible for tumour initiation, progression and development, and discusses the critical role of chronic inflammation in tumorigenesis.

Deep sea as a source of novel-anticancer drugs: update on discovery and preclinical/clinical evaluation in a systems medicine perspective

The deep-sea habitat is a source of very potent marine-derived agents that may inhibit the growth of human cancer cells "in vitro" and "in vivo". Salinosporamide-A, Marizomib, by Salinispora species is a proteasome inhibitor with promising anticancer activity (Phase I/II trials). Different deep-sea-derived drugs are under preclinical evaluation. Cancer is a complex disease that may be represented by network medicine. A simple consequence is the change of the concept of target entity from a single protein to a whole molecular pathway and or cellular network. Deep-sea-derived drugs fit well to this new concept.