Bayesian analysis of signaling networks governing embryonic stem cell fate decisions

Bridging systems biology and tissue engineering: Unleashing the full potential of complex 3D in vitro tissue models of disease

Rapid advances in tissue engineering have resulted in more complex and physiologically relevant 3D in vitro tissue models with applications in fundamental biology and therapeutic development. However, the complexity provided by these models is often not leveraged fully due to the reductionist methods used to analyze them. Computational and mathematical models developed in the field of systems biology can address this issue. Yet, traditional systems biology has been mostly applied to simpler in vitro models with little physiological relevance and limited cellular complexity. Therefore, integrating these two inherently interdisciplinary fields can result in new insights and move both disciplines forward. In this review, we provide a systematic overview of how systems biology has been integrated with 3D in vitro tissue models and discuss key application areas where the synergies between both fields have led to important advances with potential translational impact. We then outline key directions for future research and discuss a framework for further integration between fields.

Signalling pathway crosstalk stimulated by L-proline drives mouse embryonic stem cells to primitive-ectoderm-like cells

The amino acid L-proline exhibits growth factor-like properties during development - from improving blastocyst development to driving neurogenesis in vitro. Addition of 400 μM L-proline to self-renewal medium drives naïve mouse embryonic stem cells (ESCs) to early primitive ectoderm-like (EPL) cells - a transcriptionally distinct primed or partially primed pluripotent state. EPL cells retain expression of pluripotency genes, upregulate primitive ectoderm markers, undergo a morphological change and have increased cell number. These changes are facilitated by a complex signalling network hinging on the Mapk, Fgfr, Pi3k and mTor pathways. Here, we use a factorial experimental design coupled with statistical modelling to understand which signalling pathways are involved in the transition between ESCs and EPL cells, and how they underpin changes in morphology, cell number, apoptosis, proliferation and gene expression. This approach reveals pathways which work antagonistically or synergistically. Most properties were affected by more than one inhibitor, and each inhibitor blocked specific aspects of the naïve-to-primed transition. These mechanisms underpin progression of stem cells across the in vitro pluripotency continuum and serve as a model for pre-, peri- and post-implantation embryogenesis.

Discovering novel cancer bio-markers in acquired lapatinib resistance using Bayesian methods

Signalling transduction pathways (STPs) are commonly hijacked by many cancers for their growth and malignancy, but demystifying their underlying mechanisms is difficult. Here, we developed methodologies with a fully Bayesian approach in discovering novel driver bio-markers in aberrant STPs given high-throughput gene expression (GE) data. This project, namely 'PathTurbEr' (Pathway Perturbation Driver) uses the GE dataset derived from the lapatinib (an EGFR/HER dual inhibitor) sensitive and resistant samples from breast cancer cell lines (SKBR3). Differential expression analysis revealed 512 differentially expressed genes (DEGs) and their pathway enrichment revealed 13 highly perturbed singalling pathways in lapatinib resistance, including PI3K-AKT, Chemokine, Hippo and TGF-$\beta $ singalling pathways. Next, the aberration in TGF-$\beta $ STP was modelled as a causal Bayesian network (BN) using three MCMC sampling methods, i.e. Neighbourhood sampler (NS) and Hit-and-Run (HAR) sampler that potentially yield robust inference with lower chances of getting stuck at local optima and faster convergence compared to other state-of-art methods. Next, we examined the structural features of the optimal BN as a statistical process that generates the global structure using $p_1$-model, a special class of Exponential Random Graph Models (ERGMs), and MCMC methods for their hyper-parameter sampling. This step enabled key drivers identification that drive the aberration within the perturbed BN structure of STP, and yielded 34, 34 and 23 perturbation driver genes out of 80 constituent genes of three perturbed STP models of TGF-$\beta $ signalling inferred by NS, HAR and MH sampling methods, respectively. Functional-relevance and disease-relevance analyses suggested their significant associations with breast cancer progression/resistance.

Use of Computational Modeling to Study Joint Degeneration: A Review

Osteoarthritis (OA), a degenerative joint disease, is the most common chronic condition of the joints, which cannot be prevented effectively. Computational modeling of joint degradation allows to estimate the patient-specific progression of OA, which can aid clinicians to estimate the most suitable time window for surgical intervention in osteoarthritic patients. This paper gives an overview of the different approaches used to model different aspects of joint degeneration, thereby focusing mostly on the knee joint. The paper starts by discussing how OA affects the different components of the joint and how these are accounted for in the models. Subsequently, it discusses the different modeling approaches that can be used to answer questions related to OA etiology, progression and treatment. These models are ordered based on their underlying assumptions and technologies: musculoskeletal models, Finite Element models, (gene) regulatory models, multiscale models and data-driven models (artificial intelligence/machine learning). Finally, it is concluded that in the future, efforts should be made to integrate the different modeling techniques into a more robust computational framework that should not only be efficient to predict OA progression but also easily allow a patient’s individualized risk assessment as screening tool for use in clinical practice.

In silico disease model: from simple networks to complex diseases

In this chapter we consider in silico modeling of diseases starting from some simple to some complex (and mathematical) concepts. Examples and applications of in silico modeling for some important categories of diseases (such as for cancers, infectious diseases, and neuronal diseases) are also given.

Computational Modeling and Reverse Engineering to Reveal Dominant Regulatory Interactions Controlling Osteochondral Differentiation: Potential for Regenerative Medicine

The specialization of cartilage cells, or chondrogenic differentiation, is an intricate and meticulously regulated process that plays a vital role in both bone formation and cartilage regeneration. Understanding the molecular regulation of this process might help to identify key regulatory factors that can serve as potential therapeutic targets, or that might improve the development of qualitative and robust skeletal tissue engineering approaches. However, each gene involved in this process is influenced by a myriad of feedback mechanisms that keep its expression in a desirable range, making the prediction of what will happen if one of these genes defaults or is targeted with drugs, challenging. Computer modeling provides a tool to simulate this intricate interplay from a network perspective. This paper aims to give an overview of the current methodologies employed to analyze cell differentiation in the context of skeletal tissue engineering in general and osteochondral differentiation in particular. In network modeling, a network can either be derived from mechanisms and pathways that have been reported in the literature (knowledge-based approach) or it can be inferred directly from the data (data-driven approach). Combinatory approaches allow further optimization of the network. Once a network is established, several modeling technologies are available to interpret dynamically the relationships that have been put forward in the network graph (implication of the activation or inhibition of certain pathways on the evolution of the system over time) and to simulate the possible outcomes of the established network such as a given cell state. This review provides for each of the aforementioned steps (building, optimizing, and modeling the network) a brief theoretical perspective, followed by a concise overview of published works, focusing solely on applications related to cell fate decisions, cartilage differentiation and growth plate biology. Particular attention is paid to an in-house developed example of gene regulatory network modeling of growth plate chondrocyte differentiation as all the aforementioned steps can be illustrated. In summary, this paper discusses and explores a series of tools that form a first step toward a rigorous and systems-level modeling of osteochondral differentiation in the context of regenerative medicine.

Shaping Cell Fate: Influence of Topographical Substratum Properties on Embryonic Stem Cells

Development of multicellular organisms is a highly orchestrated process, with cells responding to factors and features present in the extracellular milieu. Changes in the surrounding environment help decide the fate of cells at various stages of development. This review highlights recent research that details the effects of mechanical properties of the surrounding environment and extracellular matrix and the underlying molecular mechanisms that regulate the behavior of embryonic stem cells (ESCs). In this study, we review the role of mechanical properties during embryogenesis and discuss the effect of engineered microtopographies on ESC pluripotency.

A system of recurrent neural networks for modularising, parameterising and dynamic analysis of cell signalling networks

In this paper, we show how to extend our previously proposed novel continuous time Recurrent Neural Networks (RNN) approach that retains the advantage of continuous dynamics offered by Ordinary Differential Equations (ODE) while enabling parameter estimation through adaptation, to larger signalling networks using a modular approach. Specifically, the signalling network is decomposed into several sub-models based on important temporal events in the network. Each sub-model is represented by the proposed RNN and trained using data generated from the corresponding ODE model. Trained sub-models are assembled into a whole system RNN which is then subjected to systems dynamics and sensitivity analyses. The concept is illustrated by application to G1/S transition in cell cycle using Iwamoto et al. (2008) ODE model. We decomposed the G1/S network into 3 sub-models: (i) E2F transcription factor release; (ii) E2F and CycE positive feedback loop for elevating cyclin levels; and (iii) E2F and CycA negative feedback to degrade E2F. The trained sub-models accurately represented system dynamics and parameters were in good agreement with the ODE model. The whole system RNN however revealed couple of parameters contributing to compounding errors due to feedback and required refinement to sub-model 2. These related to the reversible reaction between CycE/CDK2 and p27, its inhibitor. The revised whole system RNN model very accurately matched dynamics of the ODE system. Local sensitivity analysis of the whole system model further revealed the most dominant influence of the above two parameters in perturbing G1/S transition, giving support to a recent hypothesis that the release of inhibitor p27 from Cyc/CDK complex triggers cell cycle stage transition. To make the model useful in a practical setting, we modified each RNN sub-model with a time relay switch to facilitate larger interval input data (≈20min) (original model used data for 30s or less) and retrained them that produced parameters and protein concentrations similar to the original RNN system. Results thus demonstrated the reliability of the proposed RNN method for modelling relatively large networks by modularisation for practical settings. Advantages of the method are its ability to represent accurate continuous system dynamics and ease of: parameter estimation through training with data from a practical setting, model analysis (40% faster than ODE), fine tuning parameters when more data are available, sub-model extension when new elements and/or interactions come to light and model expansion with addition of sub-models.

Non-obvious correlations to disease management unraveled by Bayesian artificial intelligence analyses of CMS data

OBJECTIVE

Given the availability of extensive digitized healthcare data from medical records, claims and prescription information, it is now possible to use hypothesis-free, data-driven approaches to mine medical databases for novel insight. The goal of this analysis was to demonstrate the use of artificial intelligence based methods such as Bayesian networks to open up opportunities for creation of new knowledge in management of chronic conditions.

MATERIALS AND METHODS

Hospital level Medicare claims data containing discharge numbers for most common diagnoses were analyzed in a hypothesis-free manner using Bayesian networks learning methodology.

RESULTS

While many interactions identified between discharge rates of diagnoses using this data set are supported by current medical knowledge, a novel interaction linking asthma and renal failure was discovered. This interaction is non-obvious and had not been looked at by the research and clinical communities in epidemiological or clinical data. A plausible pharmacological explanation of this link is proposed together with a verification of the risk significance by conventional statistical analysis.

CONCLUSION

Potential clinical and molecular pathways defining the relationship between commonly used asthma medications and renal disease are discussed. The study underscores the need for further epidemiological research to validate this novel hypothesis. Validation will lead to advancement in clinical treatment of asthma & bronchitis, thereby, improving patient outcomes and leading to long term cost savings. In summary, this study demonstrates that application of advanced artificial intelligence methods in healthcare has the potential to enhance the quality of care by discovering non-obvious, clinically relevant relationships and enabling timely care intervention.

Mathematical Models of Pluripotent Stem Cells: At the Dawn of Predictive Regenerative Medicine

Regenerative medicine, ranging from stem cell therapy to organ regeneration, is promising to revolutionize treatments of diseases and aging. These approaches require a perfect understanding of cell reprogramming and differentiation. Predictive modeling of cellular systems has the potential to provide insights about the dynamics of cellular processes, and guide their control. Moreover in many cases, it provides alternative to experimental tests, difficult to perform for practical or ethical reasons. The variety and accuracy of biological processes represented in mathematical models grew in-line with the discovery of underlying molecular mechanisms. High-throughput data generation led to the development of models based on data analysis, as an alternative to more established modeling based on prior mechanistic knowledge. In this chapter, we give an overview of existing mathematical models of pluripotency and cell fate, to illustrate the variety of methods and questions. We conclude that current approaches are yet to overcome a number of limitations: Most of the computational models have so far focused solely on understanding the regulation of pluripotency, and the differentiation of selected cell lineages. In addition, models generally interrogate only a few biological processes. However, a better understanding of the reprogramming process leading to ESCs and iPSCs is required to improve stem-cell therapies. One also needs to understand the links between signaling, metabolism, regulation of gene expression, and the epigenetics machinery.

Bayesian network analysis of multi-compartmentalized immune responses in a murine model of sepsis and direct lung injury

Background

Inflammatory disease processes involve complex and interrelated systems of mediators. Determining the causal relationships among these mediators becomes more complicated when two, concurrent inflammatory conditions occur. In those cases, the outcome may also be dependent upon the timing, severity and compartmentalization of the insults. Unfortunately, standard methods of experimentation and analysis of data sets may investigate a single scenario without uncovering many potential associations among mediators. However, Bayesian network analysis is able to model linear, nonlinear, combinatorial, and stochastic relationships among variables to explore complex inflammatory disease systems. In these studies, we modeled the development of acute lung injury from an indirect insult (sepsis induced by cecal ligation and puncture) complicated by a direct lung insult (aspiration). To replicate multiple clinical situations, the aspiration injury was delivered at different severities and at different time intervals relative to the septic insult. For each scenario, we measured numerous inflammatory cell types and cytokines in samples from the local compartments (peritoneal and bronchoalveolar lavage fluids) and the systemic compartment (plasma). We then analyzed these data by Bayesian networks and standard methods.

Results

Standard data analysis demonstrated that the lung injury was actually reduced when two insults were involved as compared to one lung injury alone. Bayesian network analysis determined that both the severity of lung insult and presence of sepsis influenced neutrophil recruitment and the amount of injury to the lung. However, the levels of chemoattractant cytokines responsible for neutrophil recruitment were more strongly linked to the timing and severity of the lung insult compared to the presence of sepsis. This suggests that something other than sepsis-driven exacerbation of chemokine levels was influencing the lung injury, contrary to previous theories.

Conclusions

To our knowledge, these studies are the first to use Bayesian networks together with experimental studies to examine the pathogenesis of sepsis-associated lung injury. Compared to standard statistical analysis and inference, these analyses elucidated more intricate relationships among the mediators, immune cells and insult-related variables (timing, compartmentalization and severity) that cause lung injury. Bayesian networks are an effective tool for evaluating complex models of inflammation.

Inferring Aberrant Signal Transduction Pathways in Ovarian Cancer from TCGA Data

This paper concerns a new method for identifying aberrant signal transduction pathways (STPs) in cancer using case/control gene expression-level datasets, and applying that method and an existing method to an ovarian carcinoma dataset. Both methods identify STPs that are plausibly linked to all cancers based on current knowledge. Thus, the paper is most appropriate for the cancer informatics community. Our hypothesis is that STPs that are altered in tumorous tissue can be identified by applying a new Bayesian network (BN)-based method (causal analysis of STP aberration (CASA)) and an existing method (signaling pathway impact analysis (SPIA)) to the cancer genome atlas (TCGA) gene expression-level datasets. To test this hypothesis, we analyzed 20 cancer-related STPs and 6 randomly chosen STPs using the 591 cases in the TCGA ovarian carcinoma dataset, and the 102 controls in all 5 TCGA cancer datasets. We identified all the genes related to each of the 26 pathways, and developed separate gene expression datasets for each pathway. The results of the two methods were highly correlated. Furthermore, many of the STPs that ranked highest according to both methods are plausibly linked to all cancers based on current knowledge. Finally, CASA ranked the cancer-related STPs over the randomly selected STPs at a significance level below 0.05 (P = 0.047), but SPIA did not (P = 0.083).

Modeling the altered expression levels of genes on signaling pathways in tumors as causal bayesian networks

This paper concerns a study indicating that the expression levels of genes in signaling pathways can be modeled using a causal Bayesian network (BN) that is altered in tumorous tissue. These results open up promising areas of future research that can help identify driver genes and therapeutic targets. So, it is most appropriate for the cancer informatics community. Our central hypothesis is that the expression levels of genes that code for proteins on a signal transduction network (STP) are causally related and that this causal structure is altered when the STP is involved in cancer. To test this hypothesis, we analyzed 5 STPs associated with breast cancer, 7 STPs associated with other cancers, and 10 randomly chosen pathways, using a breast cancer gene expression level dataset containing 529 cases and 61 controls. We identified all the genes related to each of the 22 pathways and developed separate gene expression datasets for each pathway. We obtained significant results indicating that the causal structure of the expression levels of genes coding for proteins on STPs, which are believed to be implicated in both breast cancer and in all cancers, is more altered in the cases relative to the controls than the causal structure of the randomly chosen pathways.

Cancer systems biology and modeling: microscopic scale and multiscale approaches

Cancer has become known as a complex and systematic disease on macroscopic, mesoscopic and microscopic scales. Systems biology employs state-of-the-art computational theories and high-throughput experimental data to model and simulate complex biological procedures such as cancer, which involves genetic and epigenetic, in addition to intracellular and extracellular complex interaction networks. In this paper, different systems biology modeling techniques such as systems of differential equations, stochastic methods, Boolean networks, Petri nets, cellular automata methods and agent-based systems are concisely discussed. We have compared the mentioned formalisms and tried to address the span of applicability they can bear on emerging cancer modeling and simulation approaches. Different scales of cancer modeling, namely, microscopic, mesoscopic and macroscopic scales are explained followed by an illustration of angiogenesis in microscopic scale of the cancer modeling. Then, the modeling of cancer cell proliferation and survival are examined on a microscopic scale and the modeling of multiscale tumor growth is explained along with its advantages.

Single timepoint models of dynamic systems

Many interesting studies aimed at elucidating the connectivity structure of biomolecular pathways make use of abundance measurements, and employ statistical and information theoretic approaches to assess connectivities. These studies often do not address the effects of the dynamics of the underlying biological system, yet dynamics give rise to impactful issues such as timepoint selection and its effect on structure recovery. In this work, we study conditions for reliable retrieval of the connectivity structure of a dynamic system, and the impact of dynamics on structure-learning efforts. We encounter an unexpected problem not previously described in elucidating connectivity structure from dynamic systems, show how this confounds structure learning of the system and discuss possible approaches to overcome the confounding effect. Finally, we test our hypotheses on an accurate dynamic model of the IGF signalling pathway. We use two structure-learning methods at four time points to contrast the performance and robustness of those methods in terms of recovering correct connectivity.

Systematically studying kinase inhibitor induced signaling network signatures by integrating both therapeutic and side effects

Substantial effort in recent years has been devoted to analyzing data based large-scale biological networks, which provide valuable insight into the topologies of complex biological networks but are rarely context specific and cannot be used to predict the responses of cell signaling proteins to specific ligands or compounds. In this work, we proposed a novel strategy to investigate kinase inhibitor induced pathway signatures by integrating multiplex data in Library of Integrated Network-based Cellular Signatures (LINCS), e.g. KINOMEscan data and cell proliferation/mitosis imaging data. Using this strategy, we first established a PC9 cell line specific pathway model to investigate the pathway signatures in PC9 cell line when perturbed by a small molecule kinase inhibitor GW843682. This specific pathway revealed the role of PI3K/AKT in modulating the cell proliferation process and the absence of two anti-proliferation links, which indicated a potential mechanism of abnormal expansion in PC9 cell number. Incorporating the pathway model for side effects on primary human hepatocytes, it was used to screen 27 kinase inhibitors in LINCS database and PF02341066, known as Crizotinib, was finally suggested with an optimal concentration 4.6 uM to suppress PC9 cancer cell expansion while avoiding severe damage to primary human hepatocytes. Drug combination analysis revealed that the synergistic effect region can be predicted straightforwardly based on a threshold which is an inherent property of each kinase inhibitor. Furthermore, this integration strategy can be easily extended to other specific cell lines to be a powerful tool for drug screen before clinical trials.

A categorical network approach for discovering differentially expressed regulations in cancer

Background

The problem of efficient utilization of genome-wide expression profiles for identification and prediction of complex disease conditions is both important and challenging. Polygenic pathologies such as most types of cancer involve disregulation of many interacting genes which has prompted search for suitable statistical models for their representation. By accounting for changes in gene regulations between comparable conditions, graphical statistical models are expected to improve prediction precision.

Methods

In comparison problems with two or more experimental conditions, we represent the classes by categorical Bayesian networks that share one and the same graph structure but have class-specific probability parameters. The graph structure is learned by a score-based procedure that maximizes the difference between class probabilities using a suitable measure of divergence. The proposed framework includes an indirect model selection by adhering to a principle of optimal class separation and identifies interactions presenting significant difference between the compared conditions.

Results

We evaluate the performance of the new model against some benchmark algorithms such as support vector machine, penalized linear regression and linear Gaussian networks. The classifiers are compared by prediction accuracy across 15 different data sets from breast, lung, gastric and renal cancer studies. In addition to the demonstrated strong performance against the competitors, the proposed method is able to identify disease specific changes in gene regulations which are inaccessible by other approaches. The latter is illustrated by analyzing some gene interactions differentiating adenocarcinoma and squamous cell lung cancers.

Network-based approaches in drug discovery and early development

Identification of novel targets is a critical first step in the drug discovery and development process. Most diseases such as cancer, metabolic disorders, and neurological disorders are complex, and their pathogenesis involves multiple genetic and environmental factors. Finding a viable drug target-drug combination with high potential for yielding clinical success within the efficacy-toxicity spectrum is extremely challenging. Many examples are now available in which network-based approaches show potential for the identification of novel targets and for the repositioning of established targets. The objective of this article is to highlight network approaches for identifying novel targets with greater chances of gaining approved drugs with maximal efficacy and minimal side effects. Further enhancement of these approaches may emerge from effectively integrating computational systems biology with pharmacodynamic systems analysis. Coupling genomics, proteomics, and metabolomics databases with systems pharmacology modeling may aid in the development of disease-specific networks that can be further used to build confidence in target identification.

Modeling to optimize terminal stem cell differentiation

Embryonic stem cell (ESC), iPCs, and adult stem cells (ASCs) all are among the most promising potential treatments for heart failure, spinal cord injury, neurodegenerative diseases, and diabetes. However, considerable uncertainty in the production of ESC-derived terminally differentiated cell types has limited the efficiency of their development. To address this uncertainty, we and other investigators have begun to employ a comprehensive statistical model of ESC differentiation for determining the role of intracellular pathways (e.g., STAT3) in ESC differentiation and determination of germ layer fate. The approach discussed here applies the Baysian statistical model to cell/developmental biology combining traditional flow cytometry methodology and specific morphological observations with advanced statistical and probabilistic modeling and experimental design. The final result of this study is a unique tool and model that enhances the understanding of how and when specific cell fates are determined during differentiation. This model provides a guideline for increasing the production efficiency of therapeutically viable ESCs/iPSCs/ASC derived neurons or any other cell type and will eventually lead to advances in stem cell therapy.

Network modeling reveals steps in angiotensin peptide processing

New insights into the intrarenal renin-angiotensin (Ang) system have modified our traditional view of the system. However, many finer details of this network of peptides and associated peptidases remain unclear. We hypothesized that a computational systems biology approach, applied to peptidomic data, could help to unravel the network of enzymatic conversions. We built and refined a Bayesian network model and a dynamic systems model starting from a skeleton created with established elements of the renin-Ang system and further developed it with archived matrix-assisted laser desorption ionization-time of flight mass spectra from experiments conducted in mouse podocytes exposed to exogenous Ang substrates. The model-building process suggested previously unrecognized steps, 3 of which were confirmed in vitro, including the conversion of Ang(2-10) to Ang(2-7) by neprilysin, Ang(1-9) to Ang(2-9), and Ang(1-7) to Ang(2-7) by aminopeptidase A. These data suggest a wider role of neprilysin and aminopeptidase A in glomerular formation of bioactive Ang peptides and shunting their formation. Other steps were also suggested by the model, and supporting evidence for those steps was evaluated using model-comparison methods. Our results demonstrate that systems biology methods applied to peptidomic data are effective in identifying novel steps in the Ang peptide processing network, and these findings improve our understanding of the glomerular renin-Ang system.

Signaling network prediction by the Ontology Fingerprint enhanced Bayesian network

BACKGROUND

Despite large amounts of available genomic and proteomic data, predicting the structure and response of signaling networks is still a significant challenge. While statistical method such as Bayesian network has been explored to meet this challenge, employing existing biological knowledge for network prediction is difficult. The objective of this study is to develop a novel approach that integrates prior biological knowledge in the form of the Ontology Fingerprint to infer cell-type-specific signaling networks via data-driven Bayesian network learning; and to further use the trained model to predict cellular responses.

RESULTS

We applied our novel approach to address the Predictive Signaling Network Modeling challenge of the fourth (2009) Dialog for Reverse Engineering Assessment's and Methods (DREAM4) competition. The challenge results showed that our method accurately captured signal transduction of a network of protein kinases and phosphoproteins in that the predicted protein phosphorylation levels under all experimental conditions were highly correlated (R2 = 0.93) with the observed results. Based on the evaluation of the DREAM4 organizer, our team was ranked as one of the top five best performers in predicting network structure and protein phosphorylation activity under test conditions.

CONCLUSIONS

Bayesian network can be used to simulate the propagation of signals in cellular systems. Incorporating the Ontology Fingerprint as prior biological knowledge allows us to efficiently infer concise signaling network structure and to accurately predict cellular responses.

Identifying biological network structure, predicting network behavior, and classifying network state with High Dimensional Model Representation (HDMR)

This work presents an adapted Random Sampling - High Dimensional Model Representation (RS-HDMR) algorithm for synergistically addressing three key problems in network biology: (1) identifying the structure of biological networks from multivariate data, (2) predicting network response under previously unsampled conditions, and (3) inferring experimental perturbations based on the observed network state. RS-HDMR is a multivariate regression method that decomposes network interactions into a hierarchy of non-linear component functions. Sensitivity analysis based on these functions provides a clear physical and statistical interpretation of the underlying network structure. The advantages of RS-HDMR include efficient extraction of nonlinear and cooperative network relationships without resorting to discretization, prediction of network behavior without mechanistic modeling, robustness to data noise, and favorable scalability of the sampling requirement with respect to network size. As a proof-of-principle study, RS-HDMR was applied to experimental data measuring the single-cell response of a protein-protein signaling network to various experimental perturbations. A comparison to network structure identified in the literature and through other inference methods, including Bayesian and mutual-information based algorithms, suggests that RS-HDMR can successfully reveal a network structure with a low false positive rate while still capturing non-linear and cooperative interactions. RS-HDMR identified several higher-order network interactions that correspond to known feedback regulations among multiple network species and that were unidentified by other network inference methods. Furthermore, RS-HDMR has a better ability to predict network response under unsampled conditions in this application than the best statistical inference algorithm presented in the recent DREAM3 signaling-prediction competition. RS-HDMR can discern and predict differences in network state that arise from sources ranging from intrinsic cell-cell variability to altered experimental conditions, such as when drug perturbations are introduced. This ability ultimately allows RS-HDMR to accurately classify the experimental conditions of a given sample based on its observed network state.

Progress and prospects for stem cell engineering

Stem cells offer tremendous biomedical potential owing to their abilities to self-renew and differentiate into cell types of multiple adult tissues. Researchers and engineers have increasingly developed novel discovery technologies, theoretical approaches, and cell culture systems to investigate microenvironmental cues and cellular signaling events that control stem cell fate. Many of these technologies facilitate high-throughput investigation of microenvironmental signals and the intracellular signaling networks and machinery processing those signals into cell fate decisions. As our aggregate empirical knowledge of stem cell regulation grows, theoretical modeling with systems and computational biology methods has and will continue to be important for developing our ability to analyze and extract important conceptual features of stem cell regulation from complex data. Based on this body of knowledge, stem cell engineers will continue to develop technologies that predictably control stem cell fate with the ultimate goal of being able to accurately and economically scale up these systems for clinical-grade production of stem cell therapeutics.

Mass spectrometry-based proteomics: qualitative identification to activity-based protein profiling

Mass spectrometry has become the method of choice for proteome characterization, including multicomponent protein complexes (typically tens to hundreds of proteins) and total protein expression (up to tens of thousands of proteins), in biological samples. Qualitative sequence assignment based on MS/MS spectra is relatively well-defined, while statistical metrics for relative quantification have not completely stabilized. Nonetheless, proteomics studies have progressed to the point whereby various gene-, pathway-, or network-oriented computational frameworks may be used to place mass spectrometry data into biological context. Despite this progress, the dynamic range of protein expression remains a significant hurdle, and impedes comprehensive proteome analysis. Methods designed to enrich specific protein classes have emerged as an effective means to characterize enzymes or other catalytically active proteins that are otherwise difficult to detect in typical discovery mode proteomics experiments. Collectively, these approaches will facilitate identification of biomarkers and pathways relevant to diagnosis and treatment of human disease.

Can case study approaches speed implementation of the NRC report: "toxicity testing in the 21st century: a vision and a strategy?"

The 2007 report "Toxicity Testing in the 21st Century: A Vision and a Strategy" argued for a change in toxicity testing for environmental agents and discussed federal funding mechanisms that could be used to support this transformation within the USA. The new approach would test for in vitro perturbations of toxicity pathways using human cells with high-throughput testing platforms. The NRC report proposed a deliberate timeline, spanning about 20 years, to implement a wholesale replacement of current in-life toxicity test approaches focused on apical responses with in vitro assays. One approach to accelerating implementation is to focus on well-studied prototype compounds with known toxicity pathway targets. Through a series of carefully executed case studies with four or five pathway prototypes, the various steps required for implementation of an in vitro toxicity pathway approach to risk assessment could be developed and refined. In this article, we discuss alternative approaches for implementation and also outline advantages of a case study approach and the manner in which the case studies could be pursued using current methodologies. A case study approach would be complementary to recently proposed efforts to map the human toxome, while representing a significant extension toward more formal risk assessment compared to the profiling and prioritization approaches inherent in programs such as the EPA's ToxCast effort.

Targeted protein-omic methods are bridging the gap between proteomic and hypothesis-driven protein analysis approaches

While proteomic methods have illuminated many areas of biological protein space, many fundamental questions remain with regard to systems-level relationships between mRNAs, proteins and cell behaviors. While mass spectrometric methods offer a panoramic picture of the relative expression and modification of large numbers of proteins, they are neither optimal for the analysis of predefined targets across large numbers of samples nor for assessing differences in proteins between individual cells or cell compartments. Conversely, traditional antibody-based methods are effective at sensitively analyzing small numbers of proteins across small numbers of conditions, and can be used to analyze relative differences in protein abundance and modification between cells and cell compartments. However, traditional antibody-based approaches are not optimal for analyzing large numbers of protein abundances and modifications across many samples. In this article, we will review recent advances in methodologies and philosophies behind several microarray-based, intermediate-level, 'protein-omic' methods, including a focus on reverse-phase lysate arrays and micro-western arrays, which have been helpful for bridging gaps between large- and small-scale protein analysis approaches and have provided insight into the roles that protein systems play in several biological processes.

Querying quantitative logic models (Q2LM) to study intracellular signaling networks and cell-cytokine interactions

Mathematical models have substantially improved our ability to predict the response of a complex biological system to perturbation, but their use is typically limited by difficulties in specifying model topology and parameter values. Additionally, incorporating entities across different biological scales ranging from molecular to organismal in the same model is not trivial. Here, we present a framework called “querying quantitative logic models” (Q2LM) for building and asking questions of constrained fuzzy logic (cFL) models. cFL is a recently developed modeling formalism that uses logic gates to describe influences among entities, with transfer functions to describe quantitative dependencies. Q2LM does not rely on dedicated data to train the parameters of the transfer functions, and it permits straight-forward incorporation of entities at multiple biological scales. The Q2LM framework can be employed to ask questions such as: Which therapeutic perturbations accomplish a designated goal, and under what environmental conditions will these perturbations be effective? We demonstrate the utility of this framework for generating testable hypotheses in two examples: (i) a intracellular signaling network model; and (ii) a model for pharmacokinetics and pharmacodynamics of cell-cytokine interactions; in the latter, we validate hypotheses concerning molecular design of granulocyte colony stimulating factor.

Understanding multiple ecological responses to anthropogenic disturbance: rivers and potential flow regime change

Human-induced alteration of the natural flow regime is a major threat to freshwater ecosystems and biodiversity. The effects of hydrological alteration on the structural and functional attributes of riverine communities are expected to be multiple and complex, and they may not be described easily by a single model. Based on existing knowledge of key hydrological and ecological attributes, we explored potential effects of a flow-regulation scenario on macroinvertebrate assemblage composition and diversity in two river systems in Australia's relatively undeveloped wet-dry tropics. We used a single Bayesian belief network (BBN) to model potential changes in multiple assemblage attributes within each river type during dry and wet seasons given two flow scenarios: the current, near-natural flow condition, and flow regulation. We then used multidimensional scaling (MDS) ordination to visually summarize and compare the most probable attributes of assemblages and their environment under the different scenarios. The flow-regulation scenario provided less certainty in the ecological responses of one river type during the dry season, which reduced the ability to make predictions from the BBN outputs directly. However, visualizing the BBN results in an ordination highlighted similarities and differences between the scenarios that may have been otherwise difficult to ascertain. In particular, the MDS showed that flow regulation would reduce the seasonal differentiation in hydrology and assemblage characteristics that is expected under the current low level of development. Our approach may have wider application in understanding ecosystem responses to different river management practices and should be transferred easily to other ecosystems or biotic assemblages to provide researchers, managers, and decision makers an enhanced understanding of ecological responses to potential anthropogenic disturbance.

Framework for network modularization and Bayesian network analysis to investigate the perturbed metabolic network

BACKGROUND

Genome-scale metabolic network models have contributed to elucidating biological phenomena, and predicting gene targets to engineer for biotechnological applications. With their increasing importance, their precise network characterization has also been crucial for better understanding of the cellular physiology.

RESULTS

We herein introduce a framework for network modularization and Bayesian network analysis (FMB) to investigate organism's metabolism under perturbation. FMB reveals direction of influences among metabolic modules, in which reactions with similar or positively correlated flux variation patterns are clustered, in response to specific perturbation using metabolic flux data. With metabolic flux data calculated by constraints-based flux analysis under both control and perturbation conditions, FMB, in essence, reveals the effects of specific perturbations on the biological system through network modularization and Bayesian network analysis at metabolic modular level. As a demonstration, this framework was applied to the genetically perturbed Escherichia coli metabolism, which is a lpdA gene knockout mutant, using its genome-scale metabolic network model.

CONCLUSIONS

After all, it provides alternative scenarios of metabolic flux distributions in response to the perturbation, which are complementary to the data obtained from conventionally available genome-wide high-throughput techniques or metabolic flux analysis.

Application of Bayesian network structure learning to identify causal variant SNPs from resequencing data

Using single-nucleotide polymorphism (SNP) genotypes from the 1000 Genomes Project pilot3 data provided for Genetic Analysis Workshop 17 (GAW17), we applied Bayesian network structure learning (BNSL) to identify potential causal SNPs associated with the Affected phenotype. We focus on the setting in which target genes that harbor causal variants have already been chosen for resequencing; the goal was to detect true causal SNPs from among the measured variants in these genes. Examining all available SNPs in the known causal genes, BNSL produced a Bayesian network from which subsets of SNPs connected to the Affected outcome were identified and measured for statistical significance using the hypergeometric distribution. The exploratory phase of analysis for pooled replicates sometimes identified a set of involved SNPs that contained more true causal SNPs than expected by chance in the Asian population. Analyses of single replicates gave inconsistent results. No nominally significant results were found in analyses of African or European populations. Overall, the method was not able to identify sets of involved SNPs that included a higher proportion of true causal SNPs than expected by chance alone. We conclude that this method, as currently applied, is not effective for identifying causal SNPs that follow the simulation model for the GAW17 data set, which includes many rare causal SNPs.

RegNetB: predicting relevant regulator-gene relationships in localized prostate tumor samples

Background

A central question in cancer biology is what changes cause a healthy cell to form a tumor. Gene expression data could provide insight into this question, but it is difficult to distinguish between a gene that causes a change in gene expression from a gene that is affected by this change. Furthermore, the proteins that regulate gene expression are often themselves not regulated at the transcriptional level. Here we propose a Bayesian modeling framework we term RegNetB that uses mechanistic information about the gene regulatory network to distinguish between factors that cause a change in expression and genes that are affected by the change. We test this framework using human gene expression data describing localized prostate cancer progression.

Results

The top regulatory relationships identified by RegNetB include the regulation of RLN1, RLN2, by PAX4, the regulation of ACPP (PAP) by JUN, BACH1 and BACH2, and the co-regulation of PGC and GDF15 by MAZ and TAF8. These target genes are known to participate in tumor progression, but the suggested regulatory roles of PAX4, BACH1, BACH2, MAZ and TAF8 in the process is new.

Conclusion

Integrating gene expression data and regulatory topologies can aid in identifying potentially causal mechanisms for observed changes in gene expression.

Partially observed bipartite network analysis to identify predictive connections in transcriptional regulatory networks

Background

Messenger RNA expression is regulated by a complex interplay of different regulatory proteins. Unfortunately, directly measuring the individual activity of these regulatory proteins is difficult, leaving us with only the resulting gene expression pattern as a marker for the underlying regulatory network or regulator-gene associations. Furthermore, traditional methods to predict these regulator-gene associations do not define the relative importance of each association, leading to a large number of connections in the global regulatory network that, although true, are not useful.

Results

Here we present a Bayesian method that identifies which known transcriptional relationships in a regulatory network are consistent with a given body of static gene expression data by eliminating the non-relevant ones. The Partially Observed Bipartite Network (POBN) approach developed here is tested using E. coli expression data and a transcriptional regulatory network derived from RegulonDB. When the regulatory network for E. coli was integrated with 266 E. coli gene chip observations, POBN identified 93 out of 570 connections that were either inconsistent or not adequately supported by the expression data.

Conclusion

POBN provides a systematic way to integrate known transcriptional networks with observed gene expression data to better identify which transcriptional pathways are likely responsible for the observed gene expression pattern.

Causal modeling using network ensemble simulations of genetic and gene expression data predicts genes involved in rheumatoid arthritis

Tumor necrosis factor α (TNF-α) is a key regulator of inflammation and rheumatoid arthritis (RA). TNF-α blocker therapies can be very effective for a substantial number of patients, but fail to work in one third of patients who show no or minimal response. It is therefore necessary to discover new molecular intervention points involved in TNF-α blocker treatment of rheumatoid arthritis patients. We describe a data analysis strategy for predicting gene expression measures that are critical for rheumatoid arthritis using a combination of comprehensive genotyping, whole blood gene expression profiles and the component clinical measures of the arthritis Disease Activity Score 28 (DAS28) score. Two separate network ensembles, each comprised of 1024 networks, were built from molecular measures from subjects before and 14 weeks after treatment with TNF-α blocker. The network ensemble built from pre-treated data captures TNF-α dependent mechanistic information, while the ensemble built from data collected under TNF-α blocker treatment captures TNF-α independent mechanisms. In silico simulations of targeted, personalized perturbations of gene expression measures from both network ensembles identify transcripts in three broad categories. Firstly, 22 transcripts are identified to have new roles in modulating the DAS28 score; secondly, there are 6 transcripts that could be alternative targets to TNF-α blocker therapies, including CD86 - a component of the signaling axis targeted by Abatacept (CTLA4-Ig), and finally, 59 transcripts that are predicted to modulate the count of tender or swollen joints but not sufficiently enough to have a significant impact on DAS28.

The collection and analysis of clinical data has played a key role in providing insights into the diagnosis, prognosis and treatment of disease. However, it is imperative that molecular and genetic data also be collected and integrated into the creation of network models, which capture underlying mechanisms of disease and can be interrogated to elucidate previously unknown biology. Bringing data from the clinic to the bench completes the cycle of translational research, which we demonstrate with this work. We built disease models from genetics, whole blood gene expression profiles and the component clinical measures of rheumatoid arthritis using a data-driven approach that leverages supercomputing. Genetic factors can be utilized as a source of perturbation to the system such that causal connections between genetics, molecular entities and clinical outcomes can be inferred. The existing TNF-α blocker treatments for rheumatoid arthritis are only effective for approximately 2/3 of the affected population. We identified novel therapeutic intervention points that may lead to the development of alternatives to TNF-α blocker treatments. We believe this approach will provide improved drug discovery programs, new insights into disease progression, increased drug efficacy and novel biomarkers for chronic and complex diseases.

Toward a complete in silico, multi-layered embryonic stem cell regulatory network

Recent efforts in systematically profiling embryonic stem (ES) cells have yielded a wealth of high-throughput data. Complementarily, emerging databases and computational tools facilitate ES cell studies and further pave the way toward the in silico reconstruction of regulatory networks encompassing multiple molecular layers. Here, we briefly survey databases, algorithms, and software tools used to organize and analyze high-throughput experimental data collected to study mammalian cellular systems with a focus on ES cells. The vision of using heterogeneous data to reconstruct a complete multi-layered ES cell regulatory network is discussed. This review also provides an accompanying manually extracted dataset of different types of regulatory interactions from low-throughput experimental ES cell studies available at http://amp.pharm.mssm.edu/iscmid/literature.

Network inference from time-dependent Omics data

We provide a commented overview of the available databases for the systematic collection of pathway information and biological models essential for the interpretation of Omics data. Then, we present both the state of the art and the future challenges of network inference, a research area dealing with the deduction of reaction mechanisms from experimental Omics data. This approach represents one of the most challenging instances for making use of the huge amount of information gathered in the Omics era.

Bayesian network expansion identifies new ROS and biofilm regulators

Signaling and regulatory pathways that guide gene expression have only been partially defined for most organisms. However, given the increasing number of microarray measurements, it may be possible to reconstruct such pathways and uncover missing connections directly from experimental data. Using a compendium of microarray gene expression data obtained from Escherichia coli, we constructed a series of Bayesian network models for the reactive oxygen species (ROS) pathway as defined by EcoCyc. A consensus Bayesian network model was generated using those networks sharing the top recovered score. This microarray-based network only partially agreed with the known ROS pathway curated from the literature and databases. A top network was then expanded to predict genes that could enhance the Bayesian network model using an algorithm we termed ‘BN+1’. This expansion procedure predicted many stress-related genes (e.g., dusB and uspE), and their possible interactions with other ROS pathway genes. A term enrichment method discovered that biofilm-associated microarray data usually contained high expression levels of both uspE and gadX. The predicted involvement of gene uspE in the ROS pathway and interactions between uspE and gadX were confirmed experimentally using E. coli reporter strains. Genes gadX and uspE showed a feedback relationship in regulating each other's expression. Both genes were verified to regulate biofilm formation through gene knockout experiments. These data suggest that the BN+1 expansion method can faithfully uncover hidden or unknown genes for a selected pathway with significant biological roles. The presently reported BN+1 expansion method is a generalized approach applicable to the characterization and expansion of other biological pathways and living systems.

Cancer systems biology: a network modeling perspective

Cancer is now appreciated as not only a highly heterogenous pathology with respect to cell type and tissue origin but also as a disease involving dysregulation of multiple pathways governing fundamental cell processes such as death, proliferation, differentiation and migration. Thus, the activities of molecular networks that execute metabolic or cytoskeletal processes, or regulate these by signal transduction, are altered in a complex manner by diverse genetic mutations in concert with the environmental context. A major challenge therefore is how to develop actionable understanding of this multivariate dysregulation, with respect both to how it arises from diverse genetic mutations and to how it may be ameliorated by prospective treatments. While high-throughput experimental platform technologies ranging from genomic sequencing to transcriptomic, proteomic and metabolomic profiling are now commonly used for molecular-level characterization of tumor cells and surrounding tissues, the resulting data sets defy straightforward intuitive interpretation with respect to potential therapeutic targets or the effects of perturbation. In this review article, we will discuss how significant advances can be obtained by applying computational modeling approaches to elucidate the pathways most critically involved in tumor formation and progression, impact of particular mutations on pathway operation, consequences of altered cell behavior in tissue environments and effects of molecular therapeutics.

Using mechanistic Bayesian networks to identify downstream targets of the sonic hedgehog pathway

BACKGROUND

The topology of a biological pathway provides clues as to how a pathway operates, but rationally using this topology information with observed gene expression data remains a challenge.

RESULTS

We introduce a new general-purpose analytic method called Mechanistic Bayesian Networks (MBNs) that allows for the integration of gene expression data and known constraints within a signal or regulatory pathway to predict new downstream pathway targets. The MBN framework is implemented in an open-source Bayesian network learning package, the Python Environment for Bayesian Learning (PEBL). We demonstrate how MBNs can be used by modeling the early steps of the sonic hedgehog pathway using gene expression data from different developmental stages and genetic backgrounds in mouse. Using the MBN approach we are able to automatically identify many of the known downstream targets of the hedgehog pathway such as Gas1 and Gli1, along with a short list of likely targets such as Mig12.

CONCLUSIONS

The MBN approach shown here can easily be extended to other pathways and data types to yield a more mechanistic framework for learning genetic regulatory models.

Kinetic analysis of neurotrophin-3-mediated differentiation of embryonic stem cells into neurons

The goal of this study was to develop a kinetic analysis that could predict the behavior of embryonic stem cell-derived neural progenitor cells (ESNPCs) in response to treatment with neurotrophin-3 (NT-3). Previous studies have shown that NT-3 activates the mitogen-activated protein (MAP) kinase cascade in embryonic stem cells and promotes differentiation of ESNPCs into neurons. MAP kinase activation after NT-3 stimulation was confirmed experimentally, and a kinetic analysis was developed using rate constants obtained from the literature. Concentrations of select signaling components were estimated for ESNPCs using real-time reverse transcription polymerase chain reaction by comparing mRNA levels to those of cell types with known protein concentrations. This assumption was validated using Western blots, and incorporated into the analysis. This analysis was used to predict the minimum NT-3 concentration necessary to promote neuronal differentiation of ESNPCs based on the activation of MAP kinase. These predictions were then tested experimentally to confirm the validity of the analysis. Finally, expression of the transcription factor mammalian achate schute homolog 1 and beta-tubulin III (an early neuronal marker) was examined in response to the different NT-3 doses to confirm the link between MAP kinase activation and neuronal differentiation. Overall, this study provides insight into the kinetics of the intracellular processes that promote ESNPC differentiation to neurons.

Translating biomaterial properties to intracellular signaling

Bioactive materials present important micro-environmental cues that induce specific intracellular signaling responses which ultimately determine cell behavior. For example, vascular endothelial cells on a normal vessel wall resist inflammation and thrombosis, but the same cells seeded on an artificial vascular graft or stent do not. What makes these cells behave so differently when they are adhered to different materials? Intracellular signaling from integrins and other cell-surface receptors is an important part of the answer, but these signaling responses constitute a highly-branched, interconnected network of molecules. In order to perform rational design of biomaterials, one must understand how altering the properties of the material (micro-environment) causes changes in cell behavior, and this in turn requires understanding the complex signaling response. Systems biology and mathematical modeling aid analysis of the connectivity of this network. This review summarizes applicable systems biology and mathematical modeling techniques including ordinary differential equations-based models, principal component analysis, and Bayesian networks. Next covered is biomaterials research which studies the intracellular signaling responses generated by variation of biomaterial properties. Finally, the review details ways in which modeling has been or could be applied to better understand the link between biomaterial properties and intracellular signaling.

Regulated fluctuations in nanog expression mediate cell fate decisions in embryonic stem cells

The notion that the differentiated state of a cell population is determined simply by expression of specific marker genes is changing. In this work, the authors reveal that a pluripotent cell population comprises cells with temporal fluctuations in the expression of Nanog.

There is evidence that pluripotency of mouse embryonic stem (ES) cells is associated with the activity of a network of transcription factors with Sox2, Oct4, and Nanog at the core. Using fluorescent reporters for the expression of Nanog, we observed that a population of ES cells is best described by a dynamic distribution of Nanog expression characterized by two peaks defined by high (HN) and low (LN) Nanog expression. Typically, the LN state is 5%–20% of the total population, depending on the culture conditions. Modelling of the activity of Nanog reveals that a simple network of Oct4/Sox2 and Nanog activity can account for the observed distribution and its properties as long as the transcriptional activity is tuned by transcriptional noise. The model also predicts that the LN state is unstable, something that is born out experimentally. While in this state, cells can differentiate. We suggest that transcriptional fluctuations in Nanog expression are an essential element of the pluripotent state and that the function of Sox2, Oct4, and Nanog is to act as a network that promotes and maintains transcriptional noise to interfere with the differentiation signals.

Embryonic stem (ES) cells are a pluripotent cell population derived from early mammalian embryos. An intrinsic feature of ES cells is their phenotypic heterogeneity: they display promiscuous activation of lineage-specific genes and exhibit a fluctuating flow of differentiating cells. A gene regulatory network (GRN) centred around the transcription factors Sox2, Oct4, and Nanog is essential for the establishment and the maintenance of the pluripotent state. Previous studies had suggested that ES cells can reversibly change their state of Nanog expression without losing pluripotency. Here, we extend these studies by quantifying and monitoring the expression of Nanog in a Nanog-GFP reporter cell line. We show that Nanog levels undergo slow, random fluctuations in ES cells, giving rise to heterogeneous cell populations. We identify two states, one stable, characterized by high levels of expression (HN), and another with low levels of Nanog expression (LN), which is highly unstable. While in the LN state, cells are more likely to differentiate depending on the culture medium. Mathematical modelling shows that a simple excitable system driven by transcriptional noise can account for the observed distributions and behaviours in gene expression. Our study suggests that rather than a discrete state dependent on the fixed expression of a small set of genes, pluripotency is best represented by a state of dynamic heterogeneity of a population driven by transcriptional noise, and that the function of the gene regulatory network centred around Nanog might be to generate dynamic heterogeneities at the population level.

Learning signaling network structures with sparsely distributed data

Flow cytometric measurement of signaling protein abundances has proved particularly useful for elucidation of signaling pathway structure. The single cell nature of the data ensures a very large dataset size, providing a statistically robust dataset for structure learning. Moreover, the approach is easily scaled to many conditions in high throughput. However, the technology suffers from a dimensionality constraint: at the cutting edge, only about 12 protein species can be measured per cell, far from sufficient for most signaling pathways. Because the structure learning algorithm (in practice) requires that all variables be measured together simultaneously, this restricts structure learning to the number of variables that constitute the flow cytometer's upper dimensionality limit. To address this problem, we present here an algorithm that enables structure learning for sparsely distributed data, allowing structure learning beyond the measurement technology's upper dimensionality limit for simultaneously measurable variables. The algorithm assesses pairwise (or n-wise) dependencies, constructs "Markov neighborhoods" for each variable based on these dependencies, measures each variable in the context of its neighborhood, and performs structure learning using a constrained search.