From Phage lambda to human cancer: endogenous molecular-cellular network hypothesis

Learning a Structural and Functional Representation for Gene Expressions: To Systematically Dissect Complex Cancer Phenotypes

Cancer is a heterogeneous disease, thus one of the central problems is how to dissect the resulting complex phenotypes in terms of their biological building blocks. Computationally, this is to represent and interpret high dimensional observations through a structural and conceptual abstraction into the most influential determinants underlying the problem. The working hypothesis of this report is to consider gene interaction to be largely responsible for the manifestation of complex cancer phenotypes, thus where the representation is to be conceptualized. Here, we report a representation learning strategy combined with regularizations, in which gene expressions are described in terms of a regularized product of meta-genes and their expression levels. The meta-genes are constrained by gene interactions thus representing their original topological contexts. The expression levels are supervised by their conditional dependencies among the observations thus providing a cluster-specific constraint. We obtain both of these structural constraints using a node-based graphical model. Our representation allows the selection of more influential modules, thus implicating their possible roles in neoplastic transformations. We validate our representation strategy by its robust recognitions of various cancer phenotypes comparing with various classical methods. The modules discovered are either shared or specify for different types or stages of human cancers, all of which are consistent with literature and biology.

A Graphical Model of Smoking-Induced Global Instability in Lung Cancer

Smoking is the major cause of lung cancer and the leading cause of cancer-related death in the world. The most current view about lung cancer is no longer limited to individual genes being mutated by any carcinogenic insults from smoking. Instead, tumorigenesis is a phenotype conferred by many systematic and global alterations, leading to extensive heterogeneity and variation for both the genotypes and phenotypes of individual cancer cells. Thus, strategically it is foremost important to develop a methodology to capture any consistent and global alterations presumably shared by most of the cancerous cells for a given population. This is particularly true that almost all of the data collected from solid cancers (including lung cancers) are usually distant apart over a large span of temporal or even spatial contexts. Here, we report a multiple non-Gaussian graphical model to reconstruct the gene interaction network using two previously published gene expression datasets. Our graphical model aims to selectively detect gross structural changes at the level of gene interaction networks. Our methodology is extensively validated, demonstrating good robustness, as well as the selectivity and specificity expected based on our biological insights. In summary, gene regulatory networks are still relatively stable during presumably the early stage of neoplastic transformation. But drastic structural differences can be found between lung cancer and its normal control, including the gain of functional modules for cellular proliferations such as EGFR and PDGFRA, as well as the lost of the important IL6 module, supporting their roles as potential drug targets. Interestingly, our method can also detect early modular changes, with the ALDH3A1 and its associated interactions being strongly implicated as a potential early marker, whose activations appear to alter LCN2 module as well as its interactions with the important TP53-MDM2 circuitry. Our strategy using the graphical model to reconstruct gene interaction work with biologically-inspired constraints exemplifies the importance and beauty of biology in developing any bio-computational approach.

Endogenous Molecular-Cellular Network Cancer Theory: A Systems Biology Approach

In light of ever apparent limitation of the current dominant cancer mutation theory, a quantitative hypothesis for cancer genesis and progression, endogenous molecular-cellular network hypothesis has been proposed from the systems biology perspective, now for more than 10 years. It was intended to include both the genetic and epigenetic causes to understand cancer. Its development enters the stage of meaningful interaction with experimental and clinical data and the limitation of the traditional cancer mutation theory becomes more evident. Under this endogenous network hypothesis, we established a core working network of hepatocellular carcinoma (HCC) according to the hypothesis and quantified the working network by a nonlinear dynamical system. We showed that the two stable states of the working network reproduce the main known features of normal liver and HCC at both the modular and molecular levels. Using endogenous network hypothesis and validated working network, we explored genetic mutation pattern in cancer and potential strategies to cure or relieve HCC from a totally new perspective. Patterns of genetic mutations have been traditionally analyzed by posteriori statistical association approaches in light of traditional cancer mutation theory. One may wonder the possibility of a priori determination of any mutation regularity. Here, we found that based on the endogenous network theory the features of genetic mutations in cancers may be predicted without any prior knowledge of mutation propensities. Normal hepatocyte and cancerous hepatocyte stable states, specified by distinct patterns of expressions or activities of proteins in the network, provide means to directly identify a set of most probable genetic mutations and their effects in HCC. As the key proteins and main interactions in the network are conserved through cell types in an organism, similar mutational features may also be found in other cancers. This analysis yielded straightforward and testable predictions on an accumulated and preferred mutation spectrum in normal tissue. The validation of predicted cancer state mutation patterns demonstrates the usefulness and potential of a causal dynamical framework to understand and predict genetic mutations in cancer. We also obtained the following implication related to HCC therapy, (1) specific positive feedback loops are responsible for the maintenance of normal liver and HCC; (2) inhibiting proliferation and inflammation-related positive feedback loops, and simultaneously inducing liver-specific positive feedback loop is predicated as the potential strategy to cure or relieve HCC; (3) the genesis and regression of HCC is asymmetric. In light of the characteristic property of the nonlinear dynamical system, we demonstrate that positive feedback loops must be existed as a simple and general molecular basis for the maintenance of phenotypes such as normal liver and HCC, and regulating the positive feedback loops directly or indirectly provides potential strategies to cure or relieve HCC.

Potential landscape of high dimensional nonlinear stochastic dynamics with large noise

Quantifying stochastic processes is essential to understand many natural phenomena, particularly in biology, including the cell-fate decision in developmental processes as well as the genesis and progression of cancers. While various attempts have been made to construct potential landscape in high dimensional systems and to estimate transition rates, they are practically limited to the cases where either noise is small or detailed balance condition holds. A general and practical approach to investigate real-world nonequilibrium systems, which are typically high-dimensional and subject to large multiplicative noise and the breakdown of detailed balance, remains elusive. Here, we formulate a computational framework that can directly compute the relative probabilities between locally stable states of such systems based on a least action method, without the necessity of simulating the steady-state distribution. The method can be applied to systems with arbitrary noise intensities through A-type stochastic integration, which preserves the dynamical structure of the deterministic counterpart dynamics. We demonstrate our approach in a numerically accurate manner through solvable examples. We further apply the method to investigate the role of noise on tumor heterogeneity in a 38-dimensional network model for prostate cancer, and provide a new strategy on controlling cell populations by manipulating noise strength.

Evolutionary Origins of Cancer Driver Genes and Implications for Cancer Prognosis

The cancer atavistic theory suggests that carcinogenesis is a reverse evolution process. It is thus of great interest to explore the evolutionary origins of cancer driver genes and the relevant mechanisms underlying the carcinogenesis. Moreover, the evolutionary features of cancer driver genes could be helpful in selecting cancer biomarkers from high-throughput data. In this study, through analyzing the cancer endogenous molecular networks, we revealed that the subnetwork originating from eukaryota could control the unlimited proliferation of cancer cells, and the subnetwork originating from eumetazoa could recapitulate the other hallmarks of cancer. In addition, investigations based on multiple datasets revealed that cancer driver genes were enriched in genes originating from eukaryota, opisthokonta, and eumetazoa. These results have important implications for enhancing the robustness of cancer prognosis models through selecting the gene signatures by the gene age information.

Decoding early myelopoiesis from dynamics of core endogenous network

A decade ago mainstream molecular biologists regarded it impossible or biologically ill-motivated to understand the dynamics of complex biological phenomena, such as cancer genesis and progression, from a network perspective. Indeed, there are numerical difficulties even for those who were determined to explore along this direction. Undeterred, seven years ago a group of Chinese scientists started a program aiming to obtain quantitative connections between tumors and network dynamics. Many interesting results have been obtained. In this paper we wish to test such idea from a different angle: the connection between a normal biological process and the network dynamics. We have taken early myelopoiesis as our biological model. A standard roadmap for the cell-fate diversification during hematopoiesis has already been well established experimentally, yet little was known for its underpinning dynamical mechanisms. Compounding this difficulty there were additional experimental challenges, such as the seemingly conflicting hematopoietic roadmaps and the cell-fate inter-conversion events. With early myeloid cell-fate determination in mind, we constructed a core molecular endogenous network from well-documented gene regulation and signal transduction knowledge. Turning the network into a set of dynamical equations, we found computationally several structurally robust states. Those states nicely correspond to known cell phenotypes. We also found the states connecting those stable states. They reveal the developmental routes-how one stable state would most likely turn into another stable state. Such interconnected network among stable states enabled a natural organization of cell-fates into a multi-stable state landscape. Accordingly, both the myeloid cell phenotypes and the standard roadmap were explained mechanistically in a straightforward manner. Furthermore, recent challenging observations were also explained naturally. Moreover, the landscape visually enables a prediction of a pool of additional cell states and developmental routes, including the non-sequential and cross-branch transitions, which are testable by future experiments. In summary, the endogenous network dynamics provide an integrated quantitative framework to understand the heterogeneity and lineage commitment in myeloid progenitors.

Cancer as robust intrinsic state shaped by evolution: a key issues review

Cancer is a complex disease: its pathology cannot be properly understood in terms of independent players-genes, proteins, molecular pathways, or their simple combinations. This is similar to many-body physics of a condensed phase that many important properties are not determined by a single atom or molecule. The rapidly accumulating large 'omics' data also require a new mechanistic and global underpinning to organize for rationalizing cancer complexity. A unifying and quantitative theory was proposed by some of the present authors that cancer is a robust state formed by the endogenous molecular-cellular network, which is evolutionarily built for the developmental processes and physiological functions. Cancer state is not optimized for the whole organism. The discovery of crucial players in cancer, together with their developmental and physiological roles, in turn, suggests the existence of a hierarchical structure within molecular biology systems. Such a structure enables a decision network to be constructed from experimental knowledge. By examining the nonlinear stochastic dynamics of the network, robust states corresponding to normal physiological and abnormal pathological phenotypes, including cancer, emerge naturally. The nonlinear dynamical model of the network leads to a more encompassing understanding than the prevailing linear-additive thinking in cancer research. So far, this theory has been applied to prostate, hepatocellular, gastric cancers and acute promyelocytic leukemia with initial success. It may offer an example of carrying physics inquiring spirit beyond its traditional domain: while quantitative approaches can address individual cases, however there must be general rules/laws to be discovered in biology and medicine.

Learning Latent Variable Gaussian Graphical Model for Biomolecular Network with Low Sample Complexity

Learning a Gaussian graphical model with latent variables is ill posed when there is insufficient sample complexity, thus having to be appropriately regularized. A common choice is convex ℓ₁ plus nuclear norm to regularize the searching process. However, the best estimator performance is not always achieved with these additive convex regularizations, especially when the sample complexity is low. In this paper, we consider a concave additive regularization which does not require the strong irrepresentable condition. We use concave regularization to correct the intrinsic estimation biases from Lasso and nuclear penalty as well. We establish the proximity operators for our concave regularizations, respectively, which induces sparsity and low rankness. In addition, we extend our method to also allow the decomposition of fused structure-sparsity plus low rankness, providing a powerful tool for models with temporal information. Specifically, we develop a nontrivial modified alternating direction method of multipliers with at least local convergence. Finally, we use both synthetic and real data to validate the excellence of our method. In the application of reconstructing two-stage cancer networks, “the Warburg effect” can be revealed directly.

Endogenous network states predict gain or loss of functions for genetic mutations in hepatocellular carcinoma

Cancers have been typically characterized by genetic mutations. Patterns of such mutations have traditionally been analysed by posteriori statistical association approaches. One may ponder the possibility of a priori determination of any mutation regularity. Here by exploring biological processes implied in a mechanistic theory recently developed (the endogenous molecular-cellular network theory), we found that the features of genetic mutations in cancers may be predicted without any prior knowledge of mutation propensities. With hepatocellular carcinoma (HCC) as an example, we found that the normal hepatocyte and cancerous hepatocyte can be represented by robust stable states of one single endogenous network. These stable states, specified by distinct patterns of expressions or activities of proteins in the network, provide means to directly identify a set of most probable genetic mutations and their effects in HCC. As the key proteins and main interactions in the network are conserved through cell types in an organism, similar mutational features may also be found in other cancers. This analysis yielded straightforward and testable predictions on accumulated and preferred mutation spectra in normal tissue. The validation of predicted cancer state mutation patterns demonstrates the usefulness and potential of a causal dynamical framework to understand and predict genetic mutations in cancer.

The phenotypic equilibrium of cancer cells: From average-level stability to path-wise convergence

The phenotypic equilibrium, i.e. heterogeneous population of cancer cells tending to a fixed equilibrium of phenotypic proportions, has received much attention in cancer biology very recently. In the previous literature, some theoretical models were used to predict the experimental phenomena of the phenotypic equilibrium, which were often explained by different concepts of stabilities of the models. Here we present a stochastic multi-phenotype branching model by integrating conventional cellular hierarchy with phenotypic plasticity mechanisms of cancer cells. Based on our model, it is shown that: (i) our model can serve as a framework to unify the previous models for the phenotypic equilibrium, and then harmonizes the different kinds of average-level stabilities proposed in these models; and (ii) path-wise convergence of our model provides a deeper understanding to the phenotypic equilibrium from stochastic point of view. That is, the emergence of the phenotypic equilibrium is rooted in the stochastic nature of (almost) every sample path, the average-level stability just follows from it by averaging stochastic samples.

Biological Sources of Intrinsic and Extrinsic Noise in cI Expression of Lysogenic Phage Lambda

Genetically identical cells exposed to homogeneous environment can show remarkable phenotypic difference. To predict how phenotype is shaped, understanding of how each factor contributes is required. During gene expression processes, noise could arise either intrinsically in biochemical processes of gene expression or extrinsically from other cellular processes such as cell growth. In this work, important noise sources in gene expression of phage λ lysogen are quantified using models described by stochastic differential equations (SDEs). Results show that DNA looping has sophisticated impacts on gene expression noise: When DNA looping provides autorepression, like in wild type, it reduces noise in the system; When the autorepression is defected as it is in certain mutants, DNA looping increases expression noise. We also study how each gene operator affects the expression noise by changing the binding affinity between the gene and the transcription factor systematically. We find that the system shows extraordinarily large noise when the binding affinity is in certain range, which changes the system from monostable to bistable. In addition, we find that cell growth causes non-negligible noise, which increases with gene expression level. Quantification of noise and identification of new noise sources will provide deeper understanding on how stochasticity impacts phenotype.

Endogenous molecular-cellular hierarchical modeling of prostate carcinogenesis uncovers robust structure

We explored endogenous molecular-cellular network hypothesis for prostate cancer by constructing relevant endogenous interaction network model and analyzing its dynamical properties. Molecular regulations involved in cell proliferation, apoptosis, differentiation and metabolism are included in a hierarchical mathematical modeling scheme. This dynamical network organizes into multiple robust functional states, including physiological and pathological ones. Some states have characteristics of cancer: elevated metabolic and immune activities, high concentration of growth factors and different proliferative, apoptotic and adhesive behaviors. The molecular profile of calculated cancer state agrees with existing experiments. The modeling results have additional predictions which may be validated by further experiment: 1) Prostate supports both stem cell like and liver style proliferation; 2) While prostate supports multiple cell types, including basal, luminal and endocrine cell type differentiated from its stem cell, luminal cell is most likely to be transformed malignantly into androgen independent type cancer; 3) Retinoic acid pathway and C/EBPα are possible therapeutic targets.

Modular construction of mammalian gene circuits using TALE transcriptional repressors

An important goal of synthetic biology is the rational design and predictable implementation of synthetic gene circuits using standardized and interchangeable parts. However, engineering of complex circuits in mammalian cells is currently limited by the availability of well-characterized and orthogonal transcriptional repressors. Here, we introduce a library of 26 reversible transcription activator-like effector repressors (TALERs) that bind newly designed hybrid promoters and exert transcriptional repression through steric hindrance of key transcriptional initiation elements. We demonstrate that using the input-output transfer curves of our TALERs enables accurate prediction of the behavior of modularly assembled TALER cascade and switch circuits. We also show that TALER switches employing feedback regulation exhibit improved accuracy for microRNA-based HeLa cancer cell classification versus HEK293 cells. Our TALER library is a valuable toolkit for modular engineering of synthetic circuits, enabling programmable manipulation of mammalian cells and helping elucidate design principles of coupled transcriptional and microRNA-mediated post-transcriptional regulation.

Key network approach reveals new insight into Alzheimer's disease

Alzheimer's disease (AD) is a severe neurodegenerative disorder without curative treatment. Extensive data on pathological molecular processes have been accumulated over the last years. These data combined allows a systems biology approach to identify key regulatory elements of AD and to establish a model descriptive of the disease process which can be used for the development of therapeutic agents. In this study, the authors propose a closed network that uses a set of nodes (amyloid beta, tau, beta-secretase, glutamate, cyclin-dependent kinase 5, phosphoinositide 3-kinase and hypoxia-induced factor 1 alpha) as key elements of importance to the pathogenesis of AD. The proposed network, in total 39 nodes, is able to become a novel tool capable of providing new insights into AD, such as feedback loops. Further, it highlights interconnections between pathways and identifies their combination for therapy of AD.

Quantitative implementation of the endogenous molecular-cellular network hypothesis in hepatocellular carcinoma

A quantitative hypothesis for cancer genesis and progression-the endogenous molecular-cellular network hypothesis, intended to include both genetic and epigenetic causes of cancer-has been proposed recently. Using this hypothesis, here we address the molecular basis for maintaining normal liver and hepatocellular carcinoma (HCC), and the potential strategy to cure or relieve HCC. First, we elaborate the basic assumptions of the hypothesis and establish a core working network of HCC according to the hypothesis. Second, we quantify the working network by a nonlinear dynamical system. We show that the working network reproduces the main known features of normal liver and HCC at both the modular and molecular levels. Lastly, the validated working network reveals that (i) specific positive feedback loops are responsible for the maintenance of normal liver and HCC; (ii) inhibiting proliferation and inflammation-related positive feedback loops and simultaneously inducing a liver-specific positive feedback loop is predicated as a potential strategy to cure or relieve HCC; and (iii) the genesis and regression of HCC are asymmetric. In light of the characteristic properties of the nonlinear dynamical system, we demonstrate that positive feedback loops must exist as a simple and general molecular basis for the maintenance of heritable phenotypes, such as normal liver and HCC, and regulating the positive feedback loops directly or indirectly provides potential strategies to cure or relieve HCC.