Theoretical and in silico Analyses Reveal MYC as a Dynamic Network Biomarker in Colon and Rectal Cancer

In this article, we make a theoretical and in silico study for uncovering and evaluating biomarkers in colon and rectal cancer (CRC) by the dynamic network biomarker (DNB) theory. We propose a strategy to employ the theoretical concept of UICC TNM classification in CRC. To reveal the critical transition of CRC, the DNB algorithm was implemented to analyze the genome-wide dynamic network through temporal gene expression data. The relationship between gene sets and clinical features was evaluated by weighted gene co-expression network analysis. The results show that MYC was significantly associated with tumor amplification, tumor immune cells, and survival times. The candidate tumor suppressor genes were ZBTB16, MAL, LIFR, and SLIT2. Protein-protein interaction (PPI) analysis shows that these candidate tumor suppressor genes were significant in immune cells. Data from the Human Protein Atlas showed that a high expression of these candidate tumor suppressor genes was associated with favorable prognosis in TNM stages I-IV. In conclusion, this work provides significant and novel information regarding the TNM stage, cause, and consequences of elevated MYC expression in CRC. MYC expression levels had significant negative correlations with tumor suppressor genes and immune cells.

Computer-aided biomarker discovery for precision medicine: data resources, models and applications

Biomarkers are a class of measurable and evaluable indicators with the potential to predict disease initiation and progression. In contrast to disease-associated factors, biomarkers hold the promise to capture the changeable signatures of biological states. With methodological advances, computer-aided biomarker discovery has now become a burgeoning paradigm in the field of biomedical science. In recent years, the 'big data' term has accumulated for the systematical investigation of complex biological phenomena and promoted the flourishing of computational methods for systems-level biomarker screening. Compared with routine wet-lab experiments, bioinformatics approaches are more efficient to decode disease pathogenesis under a holistic framework, which is propitious to identify biomarkers ranging from single molecules to molecular networks for disease diagnosis, prognosis and therapy. In this review, the concept and characteristics of typical biomarker types, e.g. single molecular biomarkers, module/network biomarkers, cross-level biomarkers, etc., are explicated on the guidance of systems biology. Then, publicly available data resources together with some well-constructed biomarker databases and knowledge bases are introduced. Biomarker identification models using mathematical, network and machine learning theories are sequentially discussed. Based on network substructural and functional evidences, a novel bioinformatics model is particularly highlighted for microRNA biomarker discovery. This article aims to give deep insights into the advantages and challenges of current computational approaches for biomarker detection, and to light up the future wisdom toward precision medicine and nation-wide healthcare.

Computational medicine: quantitative modeling of complex diseases

Biological complex systems are composed of numerous components that interact within and across different scales. The ever-increasing generation of high-throughput biomedical data has given us an opportunity to develop a quantitative model of nonlinear biological systems having implications in health and diseases. Multidimensional molecular data can be modeled using various statistical methods at different scales of biological organization, such as genome, transcriptome and proteome. I will discuss recent advances in the application of computational medicine in complex diseases such as network-based studies, genome-scale metabolic modeling, kinetic modeling and support vector machines with specific examples in the field of cancer, psychiatric disorders and type 2 diabetes. The recent advances in translating these computational models in diagnosis and identification of drug targets of complex diseases are discussed, as well as the challenges researchers and clinicians are facing in taking computational medicine from the bench to bedside.

Dysfunction of PLA2G6 and CYP2C44-associated network signals imminent carcinogenesis from chronic inflammation to hepatocellular carcinoma

Little is known about how chronic inflammation contributes to the progression of hepatocellular carcinoma (HCC), especially the initiation of cancer. To uncover the critical transition from chronic inflammation to HCC and the molecular mechanisms at a network level, we analyzed the time-series proteomic data of woodchuck hepatitis virus/c-myc mice and age-matched wt-C57BL/6 mice using our dynamical network biomarker (DNB) model. DNB analysis indicated that the 5th month after birth of transgenic mice was the critical period of cancer initiation, just before the critical transition, which is consistent with clinical symptoms. Meanwhile, the DNB-associated network showed a drastic inversion of protein expression and coexpression levels before and after the critical transition. Two members of DNB, PLA2G6 and CYP2C44, along with their associated differentially expressed proteins, were found to induce dysfunction of arachidonic acid metabolism, further activate inflammatory responses through inflammatory mediator regulation of transient receptor potential channels, and finally lead to impairments of liver detoxification and malignant transition to cancer. As a c-Myc target, PLA2G6 positively correlated with c-Myc in expression, showing a trend from decreasing to increasing during carcinogenesis, with the minimal point at the critical transition or tipping point. Such trend of homologous PLA2G6 and c-Myc was also observed during human hepatocarcinogenesis, with the minimal point at high-grade dysplastic nodules (a stage just before the carcinogenesis). Our study implies that PLA2G6 might function as an oncogene like famous c-Myc during hepatocarcinogenesis, while downregulation of PLA2G6 and c-Myc could be a warning signal indicating imminent carcinogenesis.

Interactions Between Genetics, Lifestyle, and Environmental Factors for Healthcare

The occurrence and progression of diseases are strongly associated with a combination of genetic, lifestyle, and environmental factors. Understanding the interplay between genetic and nongenetic components provides deep insights into disease pathogenesis and promotes personalized strategies for people healthcare. Recently, the paradigm of systems medicine, which integrates biomedical data and knowledge at multidimensional levels, is considered to be an optimal way for disease management and clinical decision-making in the era of precision medicine. In this chapter, epigenetic-mediated genetics-lifestyle-environment interactions within specific diseases and different ethnic groups are systematically discussed, and data sources, computational models, and translational platforms for systems medicine research are sequentially presented. Moreover, feasible suggestions on precision healthcare and healthy longevity are kindly proposed based on the comprehensive review of current studies.

Identifying critical differentiation state of MCF-7 cells for breast cancer by dynamical network biomarkers

Identifying the pre-transition state just before a critical transition during a complex biological process is a challenging task, because the state of the system may show neither apparent changes nor clear phenomena before this critical transition during the biological process. By exploring rich correlation information provided by high-throughput data, the dynamical network biomarker (DNB) can identify the pre-transition state. In this work, we apply DNB to detect an early-warning signal of breast cancer on the basis of gene expression data of MCF-7 cell differentiation. We find a number of the related modules and pathways in the samples, which can be used not only as the biomarkers of cancer cells but also as the drug targets. Both functional and pathway enrichment analyses validate the results.

Unravelling personalized dysfunctional gene network of complex diseases based on differential network model

In the conventional analysis of complex diseases, the control and case samples are assumed to be of great purity. However, due to the heterogeneity of disease samples, many disease genes are even not always consistently up-/down-regulated, leading to be under-estimated. This problem will seriously influence effective personalized diagnosis or treatment. The expression variance and expression covariance can address such a problem in a network manner. But, these analyses always require multiple samples rather than one sample, which is generally not available in clinical practice for each individual. To extract the common and specific network characteristics for individual patients in this paper, a novel differential network model, e.g. personalized dysfunctional gene network, is proposed to integrate those genes with different features, such as genes with the differential gene expression (DEG), genes with the differential expression variance (DEVG) and gene-pairs with the differential expression covariance (DECG) simultaneously, to construct personalized dysfunctional networks. This model uses a new statistic-like measurement on differential information, i.e., a differential score (DEVC), to reconstruct the differential expression network between groups of normal and diseased samples; and further quantitatively evaluate different feature genes in the patient-specific network for each individual. This DEVC-based differential expression network (DEVC-net) has been applied to the study of complex diseases for prostate cancer and diabetes. (1) Characterizing the global expression change between normal and diseased samples, the differential gene networks of those diseases were found to have a new bi-coloured topological structure, where their non hub-centred sub-networks are mainly composed of genes/proteins controlling various biological processes. (2) The differential expression variance/covariance rather than differential expression is new informative sources, and can be used to identify genes or gene-pairs with discriminative power, which are ignored by traditional methods. (3) More importantly, DEVC-net is effective to measure the expression state or activity of different feature genes and their network or modules in one sample for an individual. All of these results support that DEVC-net indeed has a clear advantage to effectively extract discriminatively interpretable features of gene/protein network of one sample (i.e. personalized dysfunctional network) even when disease samples are heterogeneous, and thus can provide new features like gene-pairs, in addition to the conventional individual genes, to the analysis of the personalized diagnosis and prognosis, and a better understanding on the underlying biological mechanisms.

Identifying network biomarkers based on protein-protein interactions and expression data

Identifying effective biomarkers to battle complex diseases is an important but challenging task in biomedical research today. Molecular data of complex diseases is increasingly abundant due to the rapid advance of high throughput technologies. However, a great gap remains in identifying the massive molecular data to phenotypic changes, in particular, at a network level, i.e., a novel method for identifying network biomarkers is in pressing need to accurately classify and diagnose diseases from molecular data and shed light on the mechanisms of disease pathogenesis. Rather than seeking differential genes at an individual-molecule level, here we propose a novel method for identifying network biomarkers based on protein-protein interaction affinity (PPIA), which identify the differential interactions at a network level. Specifically, we firstly define PPIAs by estimating the concentrations of protein complexes based on the law of mass action upon gene expression data. Then we select a small and non-redundant group of protein-protein interactions and single proteins according to the PPIAs, that maximizes the discerning ability of cases from controls. This method is mathematically formulated as a linear programming, which can be efficiently solved and guarantees a globally optimal solution. Extensive results on experimental data in breast cancer demonstrate the effectiveness and efficiency of the proposed method for identifying network biomarkers, which not only can accurately distinguish the phenotypes but also provides significant biological insights at a network or pathway level. In addition, our method provides a new way to integrate static protein-protein interaction information with dynamical gene expression data.