BACKGROUND
the immediate goal of this work was to construct a computable network model for cell proliferation in non-diseased lung.
lung epithelial cells are stimulated to proliferate upon injury as a mechanism for renewal  <cit> .
alterations in the control of cell proliferation play a pivotal role in lung diseases including cancer, copd, and pulmonary fibrosis.
cancer results from both gains of inappropriate growth signaling as well as the loss of mechanisms inhibiting proliferation  <cit> .
hyperplasia of mucus-producing goblet cells and airway smooth muscle contribute to copd pathology  <cit> .
pulmonary fibrosis is characterized by excessive proliferation of lung fibroblasts, resulting in impaired lung function  <cit> .
thus, increasing the molecular understanding of the regulation of cell proliferation in the lung will serve to aid in the treatment and prevention of several lung diseases.
comprehensive and detailed pathway or network models of the processes that contribute to lung disease pathology are needed to effectively interpret modern "omics" data and to qualitatively and quantitatively compare signaling across diverse data sets.
the ultimate goal of this work is to evaluate the biological impact of xenobiotics and environmental toxins on experimental systems such as lung cell cultures or whole rodent lung.
network models representing key biological processes as they occur in non-diseased cells are crucial for this effort.
tumor cell lines and other cell contexts representing advanced disease states have genetic changes and altered signaling networks that may not be present in normal, non-diseased cells.
thus, the network model described in this report is focused on biological signaling pathways expected to be functional and to regulate cell proliferation in non-diseased lung.
many different approaches can be taken to develop biological models.
biological pathways such as those captured by kegg   <cit>  are manually drawn pathway maps linking genes to pathways; kegg pathways have limited computational value for analysis of systems biology data sets beyond directly mapping observed changes to pathways and assessing over-representation.
dynamic biochemical models, such as those commonly encoded in sbml   <cit> , are useful for assessing the dynamic behavior of biochemical systems.
however, because dynamic biochemical models require a large number of parameters, they are generally limited to representation of simplified and well-constrained biological processes, and are thus not well suited to the comprehensive evaluation of complex systems consisting of multiple inter-related signaling processes.
reverse causal reasoning  is a systems biology methodology that evaluates the statistical merit that a biological entity is active in a given system, based on automated reasoning to extrapolate back from observed biological data to plausible explanations for its cause.
rcr requires an extensive knowledgebase of biological cause and effect relationships as a substrate.
rcr has been successfully applied to identify and evaluate molecular mechanisms involved in diverse biological processes, including hypoxia-induced hemangiosarcoma, sirtuin 1-induced keratinocyte differentiation, and tumor sensitivity to akt inhibition  <cit> .
these previously published applications of rcr to experimental data have involved the analysis of diseased states.
here, we apply rcr to evaluate the biological process of cell proliferation in normal, non-diseased pulmonary cells.
the lung-focused cell proliferation network described in this paper was constructed and evaluated by applying rcr to published gene expression profiling data sets associated with measured cell proliferation endpoints in lung and related cell types.
the cell proliferation network reported here provides a detailed description of molecular processes leading to cell proliferation in the lung based on causal relationships obtained from extensive evaluation of the literature.
this novel pathway model is comprehensive and integrates core cell cycle machinery with other signaling pathways which control cell proliferation in the lung, including egf signaling, circadian clock, and hedgehog.
this pathway model is computable, and can be used for the qualitative systems-level evaluation of the complex biological processes contributing to cell proliferation pathway signaling from experimental gene expression profiling data.
construction of additional pathway models for key lung disease processes such as inflammatory signaling and response to oxidative stress is planned in order to build a comprehensive network of pathway models of lung biology relevant to lung disease.
scoring algorithms are under development to enable application of this cell proliferation network and other pathway models to the quantitative evaluation of biological impact across data sets for different lung diseases, time points, or environmental perturbations.
RESULTS
cell proliferation network construction overview
the construction of the cell proliferation network was an iterative process, summarized in figure  <dig>  the selection of biological boundaries of the model was guided by literature investigation of signaling pathways relevant to cell proliferation in the lung.
causal relationships describing cell proliferation  were added to the network model from the selventa knowledgebase , with those relationships coming from lung or lung-relevant cell types prioritized .
to avoid unintentional circularity, we excluded the causal information from the specific evaluation data sets used in this study when building and evaluating the network.
these data sets were analyzed using reverse causal reasoning , a method for identifying predictions of the activity states of biological entities  that are statistically significant and consistent with the measurements taken for a given high-throughput data set .
the rcr prediction of literature model nodes in directions consistent with the observations of cell proliferation in the experiments used to generate the gene expression data verified that the model is competent to capture mechanisms regulating proliferation.
additionally, proliferation-relevant nodes predicted by rcr which were not already represented in the literature model were used to extend the model.
using this approach, we generated a more comprehensive network with nodes derived from existing literature, as well as nodes derived from cell proliferation data sets, to create an integrated cell proliferation network .
cell proliferation network content
the cell proliferation network represents a broad collection of biological mechanisms that regulate cell proliferation in the lung, and was built using a framework that is amenable to computational analyses.
the cell proliferation network  contains  <dig> nodes,  <dig> edges ), and was constructed using information from  <dig> unique pubmed-abstracted literature sources .
nodes in the network are biological entities, such as the mrna, protein, or enzymatic activity linked to a given gene; nodes may also be cellular processes such as "cell proliferation" or phases of the cell cycle.
this fine-grained representation of biological entities allows for highly accurate qualitative modeling of biological mechanisms.
an example can be seen from the sub-network detail in figure  <dig>  showing several representative network node types, including root protein nodes , modified protein nodes  and activity nodes ) and transcriptional activity of rb <dig> ).
figure  <dig> contains a key relating the prefixes  shown in the sub-network detail to their biological meaning/interpretation.
edges are relationships between nodes and may be either non-causal or causal.
non-causal edges connect different forms of a biological entity, such as an mrna or protein complex, to its base protein  without an implied causal relationship.
causal edges are cause-effect relationships between biological entities, for example the increased kinase activity of cdk <dig> causally increases phosphorylation of rb <dig> at serine  <dig>  each causal edge is supported by a text line of evidence from a specific source reference.
additional contextual details of the relationship, such as the species and tissue/cell type in which the relationship was experimentally identified, are associated with causal edges.
for this work, we used causal edges derived only from published experiments performed in human, mouse, and rat model systems, both in vitro and in vivo.
this lung-focused, fully referenced cell proliferation network provides the most comprehensive publicly available connectivity map of the molecular mechanisms regulating proliferative processes in the lung.
summary of relevant statistics describing the content of the cell proliferation network
network boundaries, assumptions, and structure
when constructing the model using content derived from the selventa knowledgebase, some initial boundary conditions and a priori assumptions relating to tissue context and biological content were established to constrain the substance of the model to its most salient details.
tissue context boundaries
our goal was to build a network model that captures the biological mechanisms controlling cell proliferation in non-diseased mammalian lung.
to maintain the focus of the network on these elements, we determined and applied a set of rules for selecting network content.
ideally, all causal relationships comprising the network would be supported by published data from experiments conducted in non-diseased human, mouse, or rat whole lung.
thus, causal relationships with literature support coming from whole lung or normal lung cell types  were prioritized.
however, in many cases, the results of the relevant detailed experiments have not been published.
thus, as a second priority, relationships derived from cell types that are found in the normal lung , but not explicitly from lung were used.
the network was focused on relationships derived from experiments done in human systems, though relationships from mouse and rat were also included.
canonical mechanisms, such as the regulation of e2f transcription factor family members by the retinoblastoma protein rb <dig>  were included in the network even if literature support explicitly demonstrating the presence of the mechanism in lung-related cells was not identified.
it was assumed that the individual relationships within canonical mechanisms  can occur in the lung.
however, if canonical relationships with specific lung contexts were found in the literature, they were used.
if needed for completing critical mechanisms within the network, relationships with other tissue contexts were used, provided they reflected proliferative processes that can occur in the normal lung.
causal relationships derived from embryonic tissue contexts were included, as the embryonic lung represents a model for non-diseased lung cell proliferation  <cit> .
as a general rule, the use of causal relationships with tissue contexts from immortalized cell lines was limited to providing the molecular details for mechanisms in the network when these specific relationships were not available from normal cells; immortalized cell lines are highly amenable to experimental manipulation and are thus a valuable system for identifying signaling pathway details that are most likely conserved in normal cells.
relationships with tissue contexts derived from tumors or other diseased tissues were used sparingly in order to focus the content of the network to the pathways involved in normal lung cell proliferation.
biological mechanism boundaries
the cell proliferation network represents the biological mechanisms leading to cell proliferation in a specific organ, the lung.
thus, biological boundaries were designed to focus the network on the cellular processes and signaling pathways with a described role in regulating lung cell proliferation, with a particular emphasis on the proximal connections to core cell cycle machinery.
following an exhaustive search of the literature, a set of pathways were selected for inclusion, while other pathways with less direct relevance for proliferation were excluded, creating the mechanistic biological boundaries of the network.
these biological mechanism boundaries were used to ensure that the cell proliferation network represented the most relevant proliferative mechanisms that occur in the non-diseased lung.
cell proliferation can be directly or indirectly influenced by a wide range of factors, including external biological stimuli  and internal metabolic alterations .
the broad range of factors that can influence cell proliferation, coupled with the observation that many proteins involved in regulating cell proliferation have varying degrees of biological promiscuity , necessitated some additional delineations framing the biological boundaries of the network.
therefore, in addition to defining the biological content included in the network, certain processes and pathways were explicitly excluded.
specifically, inflammatory cytokine signaling, the p53-dependent dna damage response, and pathways regulating the induction of/escape from apoptosis were not included in the network.
finally, components of the core replication, transcription, and translation machinery  were considered outside the boundaries of the network.
the cell proliferation network was constructed in a modular fashion using a "building block" framework in which a core cell cycle building block is connected to additional biological pathways that contribute to cell proliferation in the lung .
these supporting blocks are peripheral to, but connected to the core cell cycle machinery regulating proliferative processes in the lung.
briefly, the five building blocks are:
cell cycle
includes canonical elements of the core machinery regulating entry and exit from the mammalian cell cycle, including but not limited to cyclin, cdk, and e2f family members.
growth factors
includes common extracellular growth factors involved in regulating lung cell proliferation, namely egf, tgf-beta, vegf, and fgf family members.
the egf family members egf and tgf-alpha play critical roles in regulating the proliferation of airway epithelial cells through egf receptor activation  <cit> .
fgf <dig> and fgf <dig>  largely through activation of fgfr2-iiib signaling, stimulate lung epithelial cell proliferation as well as regulate branching morphogenesis in the developing lung  <cit> .
vegf, a key regulator of normal angiogenesis and involved in regulating proliferation of human fetal airway epithelial cells,  <cit>  was also included.
intra- and extracellular  signaling
this block contains diverse elements of the common intra- and extracellular pathways involved in mediating lung cell proliferation, including the hedgehog, wnt, and notch signaling pathways.
hedgehog signaling regulates cell proliferation and branching morphogenesis in the developing mammalian lung  <cit> .
similarly, notch signaling controls lung cell proliferation as well as differentiation  <cit> .
elements of the wnt signaling pathway are important for mediating the proliferative processes seen following lung injury  <cit> .
the remaining areas covered by this building block are calcium signaling, mapk, hox, jak/stat, mtor, prostaglandin e <dig> , clock, and nuclear receptor signaling as relevant to lung cell proliferation.
cell interaction
includes the signal transduction pathways leading to cell proliferation that originate from the interactions of common cell adhesion molecules  and extracellular matrix components .
epigenetics
includes the main known epigenetic modulators of lung cell proliferation including the histone deacetylase  family and dna methyltransferase  family member dnmt <dig>  for this block, connections from these epigenetic mediators to the core cell cycle components  were prioritized.
network verification and expansion
selection of published cell proliferation transcriptomic data sets for verification
in order to verify the content of the network, we used publicly available data from experiments in which cell proliferation was modulated in the lung or lung relevant cell types.
specifically, we analyzed transcriptomic data sets using reverse causal reasoning , which identifies upstream controllers  that can explain the significant mrna state changes in a given transcriptomic data set.
upon completing the literature model, a search was initiated for transcriptomic data sets to verify and expand the model using public data repositories such as geo  and arrayexpress.
the ideal data set would have been collected from either whole lung or a specific untransformed lung cell type, involves a simple perturbation affecting cell proliferation , have cell proliferation phenotypic endpoint data , and have raw data available with at least three biological replicates for each sample group to clearly identify statistically significant changes in gene expression.
although this ideal data set was not found, these criteria were used to identify four "next best" data sets for these purposes .
the eif4g <dig> data set  examines gene expression changes associated with decreased cell proliferation resulting from eif4g <dig> knockdown in human breast epithelial cells   <cit> .
the rhoa data set  examines gene expression changes associated with increased cell proliferation in nih3t <dig> mouse fibroblasts, caused by the introduction of the dominant activating rhoa q63l mutation  <cit> .
the ctnnb <dig> data set  examines gene expression changes resulting from expression of constitutively active ctnnb1-lef <dig> fusion protein in embryonic lung, which causes increased cell proliferation and altered cell differentiation  <cit> .
finally, the nr3c <dig> data set  examines gene expression changes resulting from glucocorticoid receptor  knockout in embryonic mouse lung, which leads to increased cell proliferation  <cit> .
the eif4g <dig> and rhoa experiments were not performed in lung-derived cells , however were used in the network construction process due to 1) the proximity of the perturbation used to modulate cell proliferation to the mechanisms which are known to occur in lung cells and 2) the knowledge that these cell types  can be found in the normal lung.
by this reasoning, even though the gene expression studies in the eif4g <dig> and rhoa data sets were not performed in lung cells directly, we expected to observe the shared or common mechanisms regulating proliferation in the cell types commonly found in lung tissue.
reverse causal reasoning on transcriptomic data sets identifies proliferative mechanisms and verifies the literature model
we performed rcr analysis on each of these four cell proliferation transcriptomic data sets and evaluated the resulting hypotheses.
foremost, we assessed whether nodes in the cell proliferation literature model were predicted as hypotheses in directions consistent with their biological roles .
this analysis served as a means to verify the content of the literature model, as hypothesis predictions for a literature node can be taken as evidence that the particular proliferation-relevant mechanism are operating in the context of known experimentally modulated cell proliferation.
figure  <dig> shows the genstruct® technology platform heatmap key for figure  <dig>  figure  <dig>  and  <dig>  figure  <dig> and  <dig> show the rcr-predicted hypotheses from the four verification data sets which were present in the literature model.
figure  <dig> shows the predictions for many nodes in the core cell cycle block, including increased e2f <dig>   <dig>  and  <dig> activities, consistent with their published role in regulating cell proliferation in lung relevant cell types  <cit> .
in addition, predictions for increased myc activity in the rhoa and ctnnb <dig> data sets are consistent with the reported role of myc in positively regulating cell proliferation in lung and lung relevant cell types  <cit> .
in addition to predictions for increased activity of positive cell proliferation mediators in data sets where cell proliferation was experimentally induced to increase, rcr also predicted decreased activities of negative regulators of proliferation.
specifically, decreases in the transcriptional activity of rb <dig> and e2f <dig>  both known negative regulators of cell cycle progression  <cit> , were predicted in multiple data sets.
likewise, decreases in the abundance of cdkn1a or cdkn2a, cell cycle checkpoint proteins with potent anti-proliferative effects, were also predicted in all three data sets where proliferation was observed increased   <cit> .
one interesting prediction was that of decreased hras mutated at g12v.
although hras activity would be expected to increase, the hras g12v mutation leads to oncogene-induced senescence  <cit> ; therefore, this hypothesis likely reflects a transcriptional signature of decreased senescence.
rcr-predicted hypotheses appearing within the cell cycle block of literature model nodes provided verification that the proximal mechanisms regulating cell proliferation were 1) correctly present in the literature model and 2) detectable using this computational approach.
however, equally important were the predictions for nodes in the peripheral building blocks, which 1) identify additional mechanistic detail for the proliferative pathways modulated and 2) can be used together with the hypothesis predictions in the core cell cycle block to assess the coverage of the literature model by all four data sets .
for the purposes of highlighting the peripheral mechanisms involved in lung cell proliferation, hypotheses within the growth factors building block were especially well represented, including predicted increases in pdgf, fgfs  <dig>   <dig> and  <dig>  hgf, and egf and its receptors .
in particular, hypotheses for decreased fgf <dig> and fgf <dig> ) were predicted in the eif4g <dig> data set, directionally consistent with the experimental observation of decreased proliferation observed in mcf10a epithelial cells.
both fgf <dig> and fgf <dig> are critical for promoting epithelial cell proliferation in the developing respiratory epithelium  <cit> .
several egf receptor complexes and their ligands, which also play central roles in regulating normal lung cell proliferation, were also predicted as hypotheses in this analysis  <cit> .
these hypotheses were especially noticeable in the rhoa data set, which used nih3t <dig> cells as an experimental model.
although nih3t <dig> cells normally express low levels of egf family receptors and are minimally responsive to egf, rhoa activation has been shown to decrease egfr endocytosis, which could lead to increased levels of egf family responsiveness in rhoa overexpressing cells  <cit> .
hypotheses from many of the other blocks of the cell proliferation literature model are also predicted in directions consistent with the observed direction of cell proliferation in the four data sets, with nodes from the cell interaction , mapk signaling , hedgehog , and wnt/beta-catenin  blocks being particularly well represented.
despite the large number of rcr-derived hypotheses corresponding to nodes in the cell proliferation network predicted in directions consistent with increased cell proliferation, some showed a different pattern.
figure  <dig> shows the rcr-derived hypotheses corresponding to nodes in the cell proliferation network that were predicted in a direction that is opposite to what we expected based on their literature-described roles in regulating lung cell proliferation.
many of these hypotheses are pleiotropic signaling molecules, which are involved in other processes in addition to proliferation, and may result from the perturbation of non-proliferative areas of biology in the data sets examined.
for example, the "response to hypoxia" and transcriptional activity of hif1a ) predictions may be more indicative of angiogenesis than proliferation.
additionally, some of these hypotheses may be predicted in unexpected directions due to feedback mechanisms or other forms of regulation.
finally, these predictions may also result from alternative activities of these signaling molecules that have not been described in the literature, such as the microrna mir <dig>  which is still in the early stages of research into its functions.
it is important to note that none of the hypotheses predicted in unexpected directions are nodes in the core cell cycle block, an observation that further verifies the cell proliferation literature model.
this analysis supported the model as an accurate and comprehensive representation of cell proliferation in the lung.
predictions for nodes in the core cell cycle and growth factor blocks are especially robust, consistent with the key role these elements play in cell proliferation.
the analysis also confirms the ability of rcr to predict proliferative mechanisms based on transcriptomic data from multiple, independent data sets.
therefore, the proliferation literature model  appears to be very well-suited for the evaluation of mechanisms guiding lung cell proliferation using gene expression microarray data sets.
expansion of the literature model using data set-derived nodes to create the integrated model
in addition to verifying the cell proliferation literature model, rcr on the four cell proliferation data sets was used to identify other mechanisms impacting cell proliferation in the lung.
the prediction of a hypothesis in a cell proliferation data set may suggest involvement in proliferation; however, they may also reflect other biological processes that are affected by the experimental perturbations in these data sets.
therefore, each of the hypotheses predicted by rcr in these four data sets that were not already included in the model was investigated to determine its role in lung proliferation.
hypotheses that were determined to play a role in lung proliferation based on surveys of the literature were then further examined to determine how they could best be integrated into the existing literature model.
these nodes  were then added to the model, creating a more robust and comprehensive network of lung proliferation.
the literature model supplemented with these data set-derived nodes is referred to in this paper as the integrated cell proliferation network, as it takes into account not only known proliferative mechanisms operating in the lung from the literature, but also additional mechanisms determined to play a role in lung cell proliferation identified by rcr on cell proliferation data sets.
for example, the transcriptional activity of zbtb <dig> , was predicted to be increased in the ctnnb <dig> data set .
miz- <dig> is ubiquitously expressed during embryonic development and has the ability to induce growth arrest  <cit> .
recently, it has been reported that the physical interaction of miz- <dig> with myc blocks the ability of miz- <dig> to induce growth arrest, partially through removing the ability of miz- <dig> to activate p15ink4b gene expression  <cit> .
while zbtb <dig> is known to influence the transcriptional activity of myc  <cit> , and cell proliferation in other cell types, it does not yet have a direct literature-described role in regulating normal lung cell proliferation.
the data set-derived nodes added to the literature model as a result of their prediction as hypotheses in the cell proliferation data sets are designated in figure  <dig> and  <dig> by the "d" in the 'origin' column.
the content of the knowledgebase  used in this study is constantly updated with the latest scientific information.
as such, the proliferation model itself is dynamic, and has the flexibility to represent a contemporary view of lung cell proliferation as scientific knowledge advances.
rcr prediction of a given node using gene expression data sets requires a minimum of four observed rna expression changes that are consistent with the predicted change in node activity in the knowledgebase.
thus, one reason that a network node may not be predicted as a hypothesis using rcr on the cell proliferation data sets is that the knowledgebase contains too few causal connections from the node to downstream rna expressions.
to address this, we took advantage of the dynamic property of the knowledgebase to perform targeted knowledge curation around specific nodes in order to increase the likelihood of detecting them as hypotheses using rcr.
the extent of these curation efforts was limited to a subset of nodes in the proliferation network, however the structural framework used in the construction of this network allows for additional knowledge to be incorporated in the future.
evaluation of the cell proliferation network
in order to evaluate the content of the cell proliferation network we assessed the coverage of network nodes predicted by rcr  as a percentage of total network nodes that were capable of being predicted.
in all,  <dig> of the  <dig> nodes in the cell proliferation network met the minimum criteria to be predicted changed by rcr  and are termed the "possible nodes".
of these  <dig> "possible nodes", rcr predicted changes in  <dig>  in at least one of the four cell proliferation data sets.
seventy one  were predicted based on the rhoa data set alone, while  <dig> ,  <dig>  and  <dig>  were predicted based on the ctnnb <dig>  nr3c <dig>  and eif4g <dig> data sets, respectively.
notably, many of the nodes for which a prediction was not possible exert their influences on proliferation via non-transcriptional events, such as phosphorylation, degradation, etc., or have limited published information regarding their influences on gene expression.
as such, these nodes would be more likely predicted to increase or decrease when using a combination of systems biology data types .
these results further verify the cell proliferation network, as well as the method of using rcr to predict proliferative mechanisms using systems biology data.
as noted in the "network verification and expansion" section, the ideal publicly available data set for verifying the network would have adhered to collection of quality control criteria including 1) non-diseased lung tissue focus, 2) simple perturbation of primarily cell proliferation , 3) relevant endpoint data, and 4) statistical soundness.
the data sets used for evaluating the model were chosen because they all met criteria 2- <dig> detailed above, and were also done in lung cell relevant contexts.
in fact, two of the data sets were derived from experiments done in embryonic lung, and two were done in cell types that best approximated the biology occurring in lung cells.
the network described here is the first step in the larger objective of creating an integrated network of lung biology.
the cell proliferation network portrays the signaling pathways involved in normal lung cell proliferation with expanded coverage relative to existing representations.
however, it relates only a subset of the processes involved in many complex lung diseases.
for example, lung cancer is a disease of uncontrolled cell proliferation, but also involves response to dna damage and apoptosis evasion components, among others  <cit> .
similarly, chronic lung diseases such as asthma and copd involve not only alterations in the proliferative aspects of cell populations, but also profound alterations in the inflammatory response  <cit> .
in this light, a truly systemic evaluation of diseases such as these will require networks that cover multiple biological processes in a lung focused and interconnected manner.
as such, the cell proliferation network is the first of several planned networks that will be built over the coming months to capture the known universe of biological processes relevant for lung disease in a comprehensive, centralized, and computable structure .
CONCLUSIONS
cell proliferation is a complex biological process with relevance to several common lung diseases.
modern systems biology data, such as transcriptomics, are useful in unraveling the detail embedded in complex processes like cell proliferation, but require the appropriate tools.
the publicly available lung focused cell proliferation network described here represents the most comprehensive and fully referenced mechanistic representation of the signaling pathways that regulate normal lung cell proliferation in existence, and is compatible with analysis using systems biology data.
the adaptable and computable structure of the network makes it a useful tool for a wide variety of research investigators across broad scientific disciplines.

 the immediate goal of this work was to construct a computable network model for cell proliferation in non-diseased lung. lung epithelial cells are stimulated to proliferate upon injury as a mechanism for renewal. alterations in the control of cell proliferation play a pivotal role in lung diseases including cancer, copd, and pulmonary fibrosis. cancer results from both gains of inappropriate growth signaling as well as the loss of mechanisms inhibiting proliferation. hyperplasia of mucus-producing goblet cells and airway smooth muscle contribute to copd pathology. pulmonary fibrosis is characterized by excessive proliferation of lung fibroblasts, resulting in impaired lung function. thus, increasing the molecular understanding of the regulation of cell proliferation in the lung will serve to aid in the treatment and prevention of several lung diseases. comprehensive and detailed pathway or network models of the processes that contribute to lung disease pathology are needed to effectively interpret modern "omics" data and to qualitatively and quantitatively compare signaling across diverse data sets. the ultimate goal of this work is to evaluate the biological impact of xenobiotics and environmental toxins on experimental systems such as lung cell cultures or whole rodent lung. network models representing key biological processes as they occur in non-diseased cells are crucial for this effort. tumor cell lines and other cell contexts representing advanced disease states have genetic changes and altered signaling networks that may not be present in normal, non-diseased cells. thus, the network model described in this report is focused on biological signaling pathways expected to be functional and to regulate cell proliferation in non-diseased lung. many different approaches can be taken to develop biological models. biological pathways such as those captured by kegg are manually drawn pathway maps linking genes to pathways; kegg pathways have limited computational value for analysis of systems biology data sets beyond directly mapping observed changes to pathways and assessing over-representation. dynamic biochemical models, such as those commonly encoded in sbml, are useful for assessing the dynamic behavior of biochemical systems. however, because dynamic biochemical models require a large number of parameters, they are generally limited to representation of simplified and well-constrained biological processes, and are thus not well suited to the comprehensive evaluation of complex systems consisting of multiple inter-related signaling processes. reverse causal reasoning is a systems biology methodology that evaluates the statistical merit that a biological entity is active in a given system, based on automated reasoning to extrapolate back from observed biological data to plausible explanations for its cause. rcr requires an extensive knowledgebase of biological cause and effect relationships as a substrate. rcr has been successfully applied to identify and evaluate molecular mechanisms involved in diverse biological processes, including hypoxia-induced hemangiosarcoma, sirtuin 1-induced keratinocyte differentiation, and tumor sensitivity to akt inhibition. these previously published applications of rcr to experimental data have involved the analysis of diseased states. here, we apply rcr to evaluate the biological process of cell proliferation in normal, non-diseased pulmonary cells. the lung-focused cell proliferation network described in this paper was constructed and evaluated by applying rcr to published gene expression profiling data sets associated with measured cell proliferation endpoints in lung and related cell types. the cell proliferation network reported here provides a detailed description of molecular processes leading to cell proliferation in the lung based on causal relationships obtained from extensive evaluation of the literature. this novel pathway model is comprehensive and integrates core cell cycle machinery with other signaling pathways which control cell proliferation in the lung, including egf signaling, circadian clock, and hedgehog. this pathway model is computable, and can be used for the qualitative systems-level evaluation of the complex biological processes contributing to cell proliferation pathway signaling from experimental gene expression profiling data. construction of additional pathway models for key lung disease processes such as inflammatory signaling and response to oxidative stress is planned in order to build a comprehensive network of pathway models of lung biology relevant to lung disease. scoring algorithms are under development to enable application of this cell proliferation network and other pathway models to the quantitative evaluation of biological impact across data sets for different lung diseases, time points, or environmental perturbations.  cell proliferation network construction overview the construction of the cell proliferation network was an iterative process, summarized in figure the selection of biological boundaries of the model was guided by literature investigation of signaling pathways relevant to cell proliferation in the lung. causal relationships describing cell proliferation were added to the network model from the selventa knowledgebase, with those relationships coming from lung or lung-relevant cell types prioritized. to avoid unintentional circularity, we excluded the causal information from the specific evaluation data sets used in this study when building and evaluating the network. these data sets were analyzed using reverse causal reasoning, a method for identifying predictions of the activity states of biological entities that are statistically significant and consistent with the measurements taken for a given high-throughput data set. the rcr prediction of literature model nodes in directions consistent with the observations of cell proliferation in the experiments used to generate the gene expression data verified that the model is competent to capture mechanisms regulating proliferation. additionally, proliferation-relevant nodes predicted by rcr which were not already represented in the literature model were used to extend the model. using this approach, we generated a more comprehensive network with nodes derived from existing literature, as well as nodes derived from cell proliferation data sets, to create an integrated cell proliferation network. cell proliferation network content the cell proliferation network represents a broad collection of biological mechanisms that regulate cell proliferation in the lung, and was built using a framework that is amenable to computational analyses. the cell proliferation network contains nodes, edges ), and was constructed using information from unique pubmed-abstracted literature sources. nodes in the network are biological entities, such as the mrna, protein, or enzymatic activity linked to a given gene; nodes may also be cellular processes such as "cell proliferation" or phases of the cell cycle. this fine-grained representation of biological entities allows for highly accurate qualitative modeling of biological mechanisms. an example can be seen from the sub-network detail in figure showing several representative network node types, including root protein nodes, modified protein nodes and activity nodes ) and transcriptional activity of rb ). figure contains a key relating the prefixes shown in the sub-network detail to their biological meaning/interpretation. edges are relationships between nodes and may be either non-causal or causal. non-causal edges connect different forms of a biological entity, such as an mrna or protein complex, to its base protein without an implied causal relationship. causal edges are cause-effect relationships between biological entities, for example the increased kinase activity of cdk causally increases phosphorylation of rb at serine each causal edge is supported by a text line of evidence from a specific source reference. additional contextual details of the relationship, such as the species and tissue/cell type in which the relationship was experimentally identified, are associated with causal edges. for this work, we used causal edges derived only from published experiments performed in human, mouse, and rat model systems, both in vitro and in vivo. this lung-focused, fully referenced cell proliferation network provides the most comprehensive publicly available connectivity map of the molecular mechanisms regulating proliferative processes in the lung. summary of relevant statistics describing the content of the cell proliferation network network boundaries, assumptions, and structure when constructing the model using content derived from the selventa knowledgebase, some initial boundary conditions and a priori assumptions relating to tissue context and biological content were established to constrain the substance of the model to its most salient details. tissue context boundaries our goal was to build a network model that captures the biological mechanisms controlling cell proliferation in non-diseased mammalian lung. to maintain the focus of the network on these elements, we determined and applied a set of rules for selecting network content. ideally, all causal relationships comprising the network would be supported by published data from experiments conducted in non-diseased human, mouse, or rat whole lung. thus, causal relationships with literature support coming from whole lung or normal lung cell types were prioritized. however, in many cases, the results of the relevant detailed experiments have not been published. thus, as a second priority, relationships derived from cell types that are found in the normal lung, but not explicitly from lung were used. the network was focused on relationships derived from experiments done in human systems, though relationships from mouse and rat were also included. canonical mechanisms, such as the regulation of e2f transcription factor family members by the retinoblastoma protein rb were included in the network even if literature support explicitly demonstrating the presence of the mechanism in lung-related cells was not identified. it was assumed that the individual relationships within canonical mechanisms can occur in the lung. however, if canonical relationships with specific lung contexts were found in the literature, they were used. if needed for completing critical mechanisms within the network, relationships with other tissue contexts were used, provided they reflected proliferative processes that can occur in the normal lung. causal relationships derived from embryonic tissue contexts were included, as the embryonic lung represents a model for non-diseased lung cell proliferation. as a general rule, the use of causal relationships with tissue contexts from immortalized cell lines was limited to providing the molecular details for mechanisms in the network when these specific relationships were not available from normal cells; immortalized cell lines are highly amenable to experimental manipulation and are thus a valuable system for identifying signaling pathway details that are most likely conserved in normal cells. relationships with tissue contexts derived from tumors or other diseased tissues were used sparingly in order to focus the content of the network to the pathways involved in normal lung cell proliferation. biological mechanism boundaries the cell proliferation network represents the biological mechanisms leading to cell proliferation in a specific organ, the lung. thus, biological boundaries were designed to focus the network on the cellular processes and signaling pathways with a described role in regulating lung cell proliferation, with a particular emphasis on the proximal connections to core cell cycle machinery. following an exhaustive search of the literature, a set of pathways were selected for inclusion, while other pathways with less direct relevance for proliferation were excluded, creating the mechanistic biological boundaries of the network. these biological mechanism boundaries were used to ensure that the cell proliferation network represented the most relevant proliferative mechanisms that occur in the non-diseased lung. cell proliferation can be directly or indirectly influenced by a wide range of factors, including external biological stimuli and internal metabolic alterations. the broad range of factors that can influence cell proliferation, coupled with the observation that many proteins involved in regulating cell proliferation have varying degrees of biological promiscuity, necessitated some additional delineations framing the biological boundaries of the network. therefore, in addition to defining the biological content included in the network, certain processes and pathways were explicitly excluded. specifically, inflammatory cytokine signaling, the p53-dependent dna damage response, and pathways regulating the induction of/escape from apoptosis were not included in the network. finally, components of the core replication, transcription, and translation machinery were considered outside the boundaries of the network. the cell proliferation network was constructed in a modular fashion using a "building block" framework in which a core cell cycle building block is connected to additional biological pathways that contribute to cell proliferation in the lung. these supporting blocks are peripheral to, but connected to the core cell cycle machinery regulating proliferative processes in the lung. briefly, the five building blocks are: cell cycle includes canonical elements of the core machinery regulating entry and exit from the mammalian cell cycle, including but not limited to cyclin, cdk, and e2f family members. growth factors includes common extracellular growth factors involved in regulating lung cell proliferation, namely egf, tgf-beta, vegf, and fgf family members. the egf family members egf and tgf-alpha play critical roles in regulating the proliferation of airway epithelial cells through egf receptor activation. fgf and fgf largely through activation of fgfr2-iiib signaling, stimulate lung epithelial cell proliferation as well as regulate branching morphogenesis in the developing lung. vegf, a key regulator of normal angiogenesis and involved in regulating proliferation of human fetal airway epithelial cells, was also included. intra- and extracellular signaling this block contains diverse elements of the common intra- and extracellular pathways involved in mediating lung cell proliferation, including the hedgehog, wnt, and notch signaling pathways. hedgehog signaling regulates cell proliferation and branching morphogenesis in the developing mammalian lung. similarly, notch signaling controls lung cell proliferation as well as differentiation. elements of the wnt signaling pathway are important for mediating the proliferative processes seen following lung injury. the remaining areas covered by this building block are calcium signaling, mapk, hox, jak/stat, mtor, prostaglandin e, clock, and nuclear receptor signaling as relevant to lung cell proliferation. cell interaction includes the signal transduction pathways leading to cell proliferation that originate from the interactions of common cell adhesion molecules and extracellular matrix components. epigenetics includes the main known epigenetic modulators of lung cell proliferation including the histone deacetylase family and dna methyltransferase family member dnmt for this block, connections from these epigenetic mediators to the core cell cycle components were prioritized. network verification and expansion selection of published cell proliferation transcriptomic data sets for verification in order to verify the content of the network, we used publicly available data from experiments in which cell proliferation was modulated in the lung or lung relevant cell types. specifically, we analyzed transcriptomic data sets using reverse causal reasoning, which identifies upstream controllers that can explain the significant mrna state changes in a given transcriptomic data set. upon completing the literature model, a search was initiated for transcriptomic data sets to verify and expand the model using public data repositories such as geo and arrayexpress. the ideal data set would have been collected from either whole lung or a specific untransformed lung cell type, involves a simple perturbation affecting cell proliferation, have cell proliferation phenotypic endpoint data, and have raw data available with at least three biological replicates for each sample group to clearly identify statistically significant changes in gene expression. although this ideal data set was not found, these criteria were used to identify four "next best" data sets for these purposes. the eif4g data set examines gene expression changes associated with decreased cell proliferation resulting from eif4g knockdown in human breast epithelial cells. the rhoa data set examines gene expression changes associated with increased cell proliferation in nih3t mouse fibroblasts, caused by the introduction of the dominant activating rhoa q63l mutation. the ctnnb data set examines gene expression changes resulting from expression of constitutively active ctnnb1-lef fusion protein in embryonic lung, which causes increased cell proliferation and altered cell differentiation. finally, the nr3c data set examines gene expression changes resulting from glucocorticoid receptor knockout in embryonic mouse lung, which leads to increased cell proliferation. the eif4g and rhoa experiments were not performed in lung-derived cells, however were used in the network construction process due to 1) the proximity of the perturbation used to modulate cell proliferation to the mechanisms which are known to occur in lung cells and 2) the knowledge that these cell types can be found in the normal lung. by this reasoning, even though the gene expression studies in the eif4g and rhoa data sets were not performed in lung cells directly, we expected to observe the shared or common mechanisms regulating proliferation in the cell types commonly found in lung tissue. reverse causal reasoning on transcriptomic data sets identifies proliferative mechanisms and verifies the literature model we performed rcr analysis on each of these four cell proliferation transcriptomic data sets and evaluated the resulting hypotheses. foremost, we assessed whether nodes in the cell proliferation literature model were predicted as hypotheses in directions consistent with their biological roles. this analysis served as a means to verify the content of the literature model, as hypothesis predictions for a literature node can be taken as evidence that the particular proliferation-relevant mechanism are operating in the context of known experimentally modulated cell proliferation. figure shows the genstruct® technology platform heatmap key for figure figure and figure and show the rcr-predicted hypotheses from the four verification data sets which were present in the literature model. figure shows the predictions for many nodes in the core cell cycle block, including increased e2f and activities, consistent with their published role in regulating cell proliferation in lung relevant cell types. in addition, predictions for increased myc activity in the rhoa and ctnnb data sets are consistent with the reported role of myc in positively regulating cell proliferation in lung and lung relevant cell types. in addition to predictions for increased activity of positive cell proliferation mediators in data sets where cell proliferation was experimentally induced to increase, rcr also predicted decreased activities of negative regulators of proliferation. specifically, decreases in the transcriptional activity of rb and e2f both known negative regulators of cell cycle progression, were predicted in multiple data sets. likewise, decreases in the abundance of cdkn1a or cdkn2a, cell cycle checkpoint proteins with potent anti-proliferative effects, were also predicted in all three data sets where proliferation was observed increased. one interesting prediction was that of decreased hras mutated at g12v. although hras activity would be expected to increase, the hras g12v mutation leads to oncogene-induced senescence; therefore, this hypothesis likely reflects a transcriptional signature of decreased senescence. rcr-predicted hypotheses appearing within the cell cycle block of literature model nodes provided verification that the proximal mechanisms regulating cell proliferation were 1) correctly present in the literature model and 2) detectable using this computational approach. however, equally important were the predictions for nodes in the peripheral building blocks, which 1) identify additional mechanistic detail for the proliferative pathways modulated and 2) can be used together with the hypothesis predictions in the core cell cycle block to assess the coverage of the literature model by all four data sets. for the purposes of highlighting the peripheral mechanisms involved in lung cell proliferation, hypotheses within the growth factors building block were especially well represented, including predicted increases in pdgf, fgfs and hgf, and egf and its receptors. in particular, hypotheses for decreased fgf and fgf ) were predicted in the eif4g data set, directionally consistent with the experimental observation of decreased proliferation observed in mcf10a epithelial cells. both fgf and fgf are critical for promoting epithelial cell proliferation in the developing respiratory epithelium. several egf receptor complexes and their ligands, which also play central roles in regulating normal lung cell proliferation, were also predicted as hypotheses in this analysis. these hypotheses were especially noticeable in the rhoa data set, which used nih3t cells as an experimental model. although nih3t cells normally express low levels of egf family receptors and are minimally responsive to egf, rhoa activation has been shown to decrease egfr endocytosis, which could lead to increased levels of egf family responsiveness in rhoa overexpressing cells. hypotheses from many of the other blocks of the cell proliferation literature model are also predicted in directions consistent with the observed direction of cell proliferation in the four data sets, with nodes from the cell interaction, mapk signaling, hedgehog, and wnt/beta-catenin blocks being particularly well represented. despite the large number of rcr-derived hypotheses corresponding to nodes in the cell proliferation network predicted in directions consistent with increased cell proliferation, some showed a different pattern. figure shows the rcr-derived hypotheses corresponding to nodes in the cell proliferation network that were predicted in a direction that is opposite to what we expected based on their literature-described roles in regulating lung cell proliferation. many of these hypotheses are pleiotropic signaling molecules, which are involved in other processes in addition to proliferation, and may result from the perturbation of non-proliferative areas of biology in the data sets examined. for example, the "response to hypoxia" and transcriptional activity of hif1a ) predictions may be more indicative of angiogenesis than proliferation. additionally, some of these hypotheses may be predicted in unexpected directions due to feedback mechanisms or other forms of regulation. finally, these predictions may also result from alternative activities of these signaling molecules that have not been described in the literature, such as the microrna mir which is still in the early stages of research into its functions. it is important to note that none of the hypotheses predicted in unexpected directions are nodes in the core cell cycle block, an observation that further verifies the cell proliferation literature model. this analysis supported the model as an accurate and comprehensive representation of cell proliferation in the lung. predictions for nodes in the core cell cycle and growth factor blocks are especially robust, consistent with the key role these elements play in cell proliferation. the analysis also confirms the ability of rcr to predict proliferative mechanisms based on transcriptomic data from multiple, independent data sets. therefore, the proliferation literature model appears to be very well-suited for the evaluation of mechanisms guiding lung cell proliferation using gene expression microarray data sets. expansion of the literature model using data set-derived nodes to create the integrated model in addition to verifying the cell proliferation literature model, rcr on the four cell proliferation data sets was used to identify other mechanisms impacting cell proliferation in the lung. the prediction of a hypothesis in a cell proliferation data set may suggest involvement in proliferation; however, they may also reflect other biological processes that are affected by the experimental perturbations in these data sets. therefore, each of the hypotheses predicted by rcr in these four data sets that were not already included in the model was investigated to determine its role in lung proliferation. hypotheses that were determined to play a role in lung proliferation based on surveys of the literature were then further examined to determine how they could best be integrated into the existing literature model. these nodes were then added to the model, creating a more robust and comprehensive network of lung proliferation. the literature model supplemented with these data set-derived nodes is referred to in this paper as the integrated cell proliferation network, as it takes into account not only known proliferative mechanisms operating in the lung from the literature, but also additional mechanisms determined to play a role in lung cell proliferation identified by rcr on cell proliferation data sets. for example, the transcriptional activity of zbtb, was predicted to be increased in the ctnnb data set. miz- is ubiquitously expressed during embryonic development and has the ability to induce growth arrest. recently, it has been reported that the physical interaction of miz- with myc blocks the ability of miz- to induce growth arrest, partially through removing the ability of miz- to activate p15ink4b gene expression. while zbtb is known to influence the transcriptional activity of myc, and cell proliferation in other cell types, it does not yet have a direct literature-described role in regulating normal lung cell proliferation. the data set-derived nodes added to the literature model as a result of their prediction as hypotheses in the cell proliferation data sets are designated in figure and by the "d" in the'origin' column. the content of the knowledgebase used in this study is constantly updated with the latest scientific information. as such, the proliferation model itself is dynamic, and has the flexibility to represent a contemporary view of lung cell proliferation as scientific knowledge advances. rcr prediction of a given node using gene expression data sets requires a minimum of four observed rna expression changes that are consistent with the predicted change in node activity in the knowledgebase. thus, one reason that a network node may not be predicted as a hypothesis using rcr on the cell proliferation data sets is that the knowledgebase contains too few causal connections from the node to downstream rna expressions. to address this, we took advantage of the dynamic property of the knowledgebase to perform targeted knowledge curation around specific nodes in order to increase the likelihood of detecting them as hypotheses using rcr. the extent of these curation efforts was limited to a subset of nodes in the proliferation network, however the structural framework used in the construction of this network allows for additional knowledge to be incorporated in the future. evaluation of the cell proliferation network in order to evaluate the content of the cell proliferation network we assessed the coverage of network nodes predicted by rcr as a percentage of total network nodes that were capable of being predicted. in all, of the nodes in the cell proliferation network met the minimum criteria to be predicted changed by rcr and are termed the "possible nodes". of these "possible nodes", rcr predicted changes in in at least one of the four cell proliferation data sets. seventy one were predicted based on the rhoa data set alone, while, and were predicted based on the ctnnb nr3c and eif4g data sets, respectively. notably, many of the nodes for which a prediction was not possible exert their influences on proliferation via non-transcriptional events, such as phosphorylation, degradation, etc., or have limited published information regarding their influences on gene expression. as such, these nodes would be more likely predicted to increase or decrease when using a combination of systems biology data types. these results further verify the cell proliferation network, as well as the method of using rcr to predict proliferative mechanisms using systems biology data. as noted in the "network verification and expansion" section, the ideal publicly available data set for verifying the network would have adhered to collection of quality control criteria including 1) non-diseased lung tissue focus, 2) simple perturbation of primarily cell proliferation, 3) relevant endpoint data, and 4) statistical soundness. the data sets used for evaluating the model were chosen because they all met criteria 2- detailed above, and were also done in lung cell relevant contexts. in fact, two of the data sets were derived from experiments done in embryonic lung, and two were done in cell types that best approximated the biology occurring in lung cells. the network described here is the first step in the larger objective of creating an integrated network of lung biology. the cell proliferation network portrays the signaling pathways involved in normal lung cell proliferation with expanded coverage relative to existing representations. however, it relates only a subset of the processes involved in many complex lung diseases. for example, lung cancer is a disease of uncontrolled cell proliferation, but also involves response to dna damage and apoptosis evasion components, among others. similarly, chronic lung diseases such as asthma and copd involve not only alterations in the proliferative aspects of cell populations, but also profound alterations in the inflammatory response. in this light, a truly systemic evaluation of diseases such as these will require networks that cover multiple biological processes in a lung focused and interconnected manner. as such, the cell proliferation network is the first of several planned networks that will be built over the coming months to capture the known universe of biological processes relevant for lung disease in a comprehensive, centralized, and computable structure.  cell proliferation is a complex biological process with relevance to several common lung diseases. modern systems biology data, such as transcriptomics, are useful in unraveling the detail embedded in complex processes like cell proliferation, but require the appropriate tools. the publicly available lung focused cell proliferation network described here represents the most comprehensive and fully referenced mechanistic representation of the signaling pathways that regulate normal lung cell proliferation in existence, and is compatible with analysis using systems biology data. the adaptable and computable structure of the network makes it a useful tool for a wide variety of research investigators across broad scientific disciplines. 
to the best of our knowledge, this lung-focused cell proliferation network provides the most comprehensive connectivity map in existence of the molecular mechanisms regulating cell proliferation in the lung. the network is based on fully referenced causal relationships obtained from extensive evaluation of the literature. the computable structure of the network enables its application to the qualitative and quantitative evaluation of cell proliferation using systems biology data sets. the network was verified using four published gene expression profiling data sets associated with measured cell proliferation endpoints in lung and lung-related cell types. lung cell proliferation is a key biological process to capture in such a network model, given the pivotal role that proliferation plays in lung diseases including cancer, chronic obstructive pulmonary disease, and fibrosis. the network encompasses diverse biological areas that lead to the regulation of normal lung cell proliferation, and contains a total of nodes and edges. ideally, these networks will be specifically designed to capture the normal, non-diseased biology of the tissue or cell types under investigation, and can be used with experimentally generated systems biology data to assess the biological impact of perturbations like xenobiotics and other cellular stresses. to further a systems-level assessment of the biological impact of perturbations on non-diseased mammalian lung cells, we constructed a lung-focused network for cell proliferation. critical to advancing the systems-level evaluation of complex biological processes is the development of comprehensive networks and computational methods to apply to the analysis of systems biology data. predicted changes in the activity of core machinery involved in cell cycle regulation are statistically supported across multiple data sets, underscoring the general applicability of this approach for a network-wide biological impact assessment using systems biology data.
 critical to advancing the systems-level evaluation of complex biological processes is the development of comprehensive networks and computational methods to apply to the analysis of systems biology data. ideally, these networks will be specifically designed to capture the normal, non-diseased biology of the tissue or cell types under investigation, and can be used with experimentally generated systems biology data to assess the biological impact of perturbations like xenobiotics and other cellular stresses. lung cell proliferation is a key biological process to capture in such a network model, given the pivotal role that proliferation plays in lung diseases including cancer, chronic obstructive pulmonary disease, and fibrosis. unfortunately, no such network has been available prior to this work.  to further a systems-level assessment of the biological impact of perturbations on non-diseased mammalian lung cells, we constructed a lung-focused network for cell proliferation. the network encompasses diverse biological areas that lead to the regulation of normal lung cell proliferation, and contains a total of nodes and edges. the network was verified using four published gene expression profiling data sets associated with measured cell proliferation endpoints in lung and lung-related cell types. predicted changes in the activity of core machinery involved in cell cycle regulation are statistically supported across multiple data sets, underscoring the general applicability of this approach for a network-wide biological impact assessment using systems biology data.  to the best of our knowledge, this lung-focused cell proliferation network provides the most comprehensive connectivity map in existence of the molecular mechanisms regulating cell proliferation in the lung. the network is based on fully referenced causal relationships obtained from extensive evaluation of the literature. the computable structure of the network enables its application to the qualitative and quantitative evaluation of cell proliferation using systems biology data sets. the network is available for public use. 