Blood-based identification of non-responders to anti-TNF therapy in rheumatoid arthritis.

BACKGROUND: Faced with an increasing number of choices for biologic therapies, rheumatologists have a critical need for better tools to inform rheumatoid arthritis (RA) disease management. The ability to identify patients who are unlikely to respond to first-line biologic anti-TNF therapies prior to their treatment would allow these patients to seek alternative therapies, providing faster relief and avoiding complications of disease. METHODS: We identified a gene expression classifier to predict, pre-treatment, which RA patients are unlikely to respond to the anti-TNF infliximab. The classifier was trained and independently evaluated using four published whole blood gene expression data sets, in which RA patients (n = 116 = 44 + 15 + 30 + 27) were treated with infliximab, and their response assessed 14-16 months post treatment according to the European League Against Rheumatism (EULAR) response criteria. For each patient, prior knowledge was used to group gene expression measurements into disease-relevant biological signaling mechanisms that were used as the input features for regularized logistic regression. RESULTS: The classifier produced a substantial enrichment of non-responders (59 %, given by the cross validated test precision) compared to the full population (27 % non-responders), while identifying nearly a third of non-responders. Given this classifier performance, treatment of predicted non-responders with alternative biologics would decrease their chance of non-response by between a third and a half, substantially improving their odds of effective treatment and stemming further disease progression. The classifier consisted of 18 signaling mechanisms, which together indicated that higher inflammatory signaling mediated by TNF and other cytokines was present pre-treatment in the blood of patients who responded to infliximab treatment. In contrast, non-responders were classified by relatively higher levels of specific metabolic activities in the blood prior to treatment. CONCLUSIONS: We were able to successfully produce a classifier to identify a population of RA patients significantly enriched in anti-TNF non-responders across four different patient cohorts. Additional prospective studies are needed to validate and refine the classifier for clinical use.

A vascular biology network model focused on inflammatory processes to investigate atherogenesis and plaque instability.

BACKGROUND: Numerous inflammation-related pathways have been shown to play important roles in atherogenesis. Rapid and efficient assessment of the relative influence of each of those pathways is a challenge in the era of "omics" data generation. The aim of the present work was to develop a network model of inflammation-related molecular pathways underlying vascular disease to assess the degree of translatability of preclinical molecular data to the human clinical setting. METHODS: We constructed and evaluated the Vascular Inflammatory Processes Network (V-IPN), a model representing a collection of vascular processes modulated by inflammatory stimuli that lead to the development of atherosclerosis. RESULTS: Utilizing the V-IPN as a platform for biological discovery, we have identified key vascular processes and mechanisms captured by gene expression profiling data from four independent datasets from human endothelial cells (ECs) and human and murine intact vessels. Primary ECs in culture from multiple donors revealed a richer mapping of mechanisms identified by the V-IPN compared to an immortalized EC line. Furthermore, an evaluation of gene expression datasets from aortas of old ApoE-/- mice (78 weeks) and human coronary arteries with advanced atherosclerotic lesions identified significant commonalities in the two species, as well as several mechanisms specific to human arteries that are consistent with the development of unstable atherosclerotic plaques. CONCLUSIONS: We have generated a new biological network model of atherogenic processes that demonstrates the power of network analysis to advance integrative, systems biology-based knowledge of cross-species translatability, plaque development and potential mechanisms leading to plaque instability.

A modular cell-type focused inflammatory process network model for non-diseased pulmonary tissue.

Exposure to environmental stressors such as cigarette smoke (CS) elicits a variety of biological responses in humans, including the induction of inflammatory responses. These responses are especially pronounced in the lung, where pulmonary cells sit at the interface between the body's internal and external environments. We combined a literature survey with a computational analysis of multiple transcriptomic data sets to construct a computable causal network model (the Inflammatory Process Network (IPN)) of the main pulmonary inflammatory processes. The IPN model predicted decreased epithelial cell barrier defenses and increased mucus hypersecretion in human bronchial epithelial cells, and an attenuated pro-inflammatory (M1) profile in alveolar macrophages following exposure to CS, consistent with prior results. The IPN provides a comprehensive framework of experimentally supported pathways related to CS-induced pulmonary inflammation. The IPN is freely available to the scientific community as a resource with broad applicability to study the pathogenesis of pulmonary disease.

Construction of a computable network model for DNA damage, autophagy, cell death, and senescence.

Towards the development of a systems biology-based risk assessment approach for environmental toxicants, including tobacco products in a systems toxicology setting such as the "21st Century Toxicology", we are building a series of computable biological network models specific to non-diseased pulmonary and cardiovascular cells/tissues which capture the molecular events that can be activated following exposure to environmental toxicants. Here we extend on previous work and report on the construction and evaluation of a mechanistic network model focused on DNA damage response and the four main cellular fates induced by stress: autophagy, apoptosis, necroptosis, and senescence. In total, the network consists of 34 sub-models containing 1052 unique nodes and 1538 unique edges which are supported by 1231 PubMed-referenced literature citations. Causal node-edge relationships are described using the Biological Expression Language (BEL), which allows for the semantic representation of life science relationships in a computable format. The Network is provided in .XGMML format and can be viewed using freely available network visualization software, such as Cytoscape.

Comparative transcriptional network modeling of three PPAR-α/γ co-agonists reveals distinct metabolic gene signatures in primary human hepatocytes.

AIMS: To compare the molecular and biologic signatures of a balanced dual peroxisome proliferator-activated receptor (PPAR)-α/γ agonist, aleglitazar, with tesaglitazar (a dual PPAR-α/γ agonist) or a combination of pioglitazone (Pio; PPAR-γ agonist) and fenofibrate (Feno; PPAR-α agonist) in human hepatocytes. METHODS AND RESULTS: Gene expression microarray profiles were obtained from primary human hepatocytes treated with EC(50)-aligned low, medium and high concentrations of the three treatments. A systems biology approach, Causal Network Modeling, was used to model the data to infer upstream molecular mechanisms that may explain the observed changes in gene expression. Aleglitazar, tesaglitazar and Pio/Feno each induced unique transcriptional signatures, despite comparable core PPAR signaling. Although all treatments inferred qualitatively similar PPAR-α signaling, aleglitazar was inferred to have greater effects on high- and low-density lipoprotein cholesterol levels than tesaglitazar and Pio/Feno, due to a greater number of gene expression changes in pathways related to high-density and low-density lipoprotein metabolism. Distinct transcriptional and biologic signatures were also inferred for stress responses, which appeared to be less affected by aleglitazar than the comparators. In particular, Pio/Feno was inferred to increase NFE2L2 activity, a key component of the stress response pathway, while aleglitazar had no significant effect. All treatments were inferred to decrease proliferative signaling. CONCLUSIONS: Aleglitazar induces transcriptional signatures related to lipid parameters and stress responses that are unique from other dual PPAR-α/γ treatments. This may underlie observed favorable changes in lipid profiles in animal and clinical studies with aleglitazar and suggests a differentiated gene profile compared with other dual PPAR-α/γ agonist treatments.

Genes associated with intestinal permeability in ulcerative colitis: changes in expression following infliximab therapy.

BACKGROUND: Alterations in intestinal permeability have been implicated in ulcerative colitis (UC). Infliximab, a monoclonal anti-tumor necrosis factor alpha (TNFα) antibody, can induce clinical response in UC. Gene expression in colonic biopsies taken from responders and nonresponders to infliximab can provide insight into the mechanisms of the altered intestinal permeability at a molecular level. METHODS: Colonic biopsies (n = 18 anti-TNFα naïve UC patients; n = 8 normal controls; n = 80 Active Ulcerative Colitis Trial [ACT] 1 patients) were analyzed for mRNA expression using gene expression microarrays. Computational reverse causal reasoning was applied to build causal network models of UC and response and nonresponse of UC to treatment. Quantitative reverse-transcription polymerase chain reaction (qPCR) was used to confirm differentially expressed genes. RESULTS: Reverse causal reasoning on mRNA expression data from anti-TNFα-naïve UC and normal samples provided a mechanistic disease model of the biology of gene expression observed in UC. mRNA expression data from the ACT 1 study enabled construction of a mechanistic model describing the biology of nonresponders to infliximab, including evidence for increased intestinal permeability compared with normal and responder samples. Gene expression changes identified as central to intestinal permeability dysregulation were confirmed in normal, UC, and infliximab-treated patients by qPCR analysis. Gene expression returned toward normal levels in infliximab responders, but not in nonresponders. CONCLUSION: Gene expression analysis and causal network modeling in combination showed that aberrant mRNA expression of genes involved in intestinal epithelial permeability for infliximab responders was restored toward levels observed in normal samples. Infliximab nonresponders showed no equivalent restoration in the expression of these genes.

Early patient stratification and predictive biomarkers in drug discovery and development: a case study of ulcerative colitis anti-TNF therapy.

The current drug discovery paradigm is long, costly, and prone to failure. For projects in early development, lack of efficacy in Phase II is a major contributor to the overall failure rate. Efficacy failures often occur from one of two major reasons: either the investigational agent did not achieve the required pharmacology or the mechanism targeted by the investigational agent did not significantly contribute to the disease in the tested patient population. The latter scenario can arise due to insufficient study power stemming from patient heterogeneity. If the subset of disease patients driven by the mechanism that is likely to respond to the drug can be identified and selected before enrollment begins, efficacy and response rates should improve. This will not only augment drug approval percentages, but will also minimize the number of patients at risk of side effects in the face of a suboptimal response to treatment. Here we describe a systems biology approach using molecular profiling data from patients at baseline for the development of predictive biomarker content to identify potential responders to a molecular targeted therapy before the drug is tested in humans. A case study is presented where a classifier to predict response to a TNF targeted therapy for ulcerative colitis is developed a priori and verified against a test set of patients where clinical outcomes are known. This approach will promote the tandem development of drugs with predictive response, patient selection biomarkers.

A computable cellular stress network model for non-diseased pulmonary and cardiovascular tissue.

BACKGROUND: Humans and other organisms are equipped with a set of responses that can prevent damage from exposure to a multitude of endogenous and environmental stressors. If these stress responses are overwhelmed, this can result in pathogenesis of diseases, which is reflected by an increased development of, e.g., pulmonary and cardiac diseases in humans exposed to chronic levels of environmental stress, including inhaled cigarette smoke (CS). Systems biology data sets (e.g., transcriptomics, phosphoproteomics, metabolomics) could enable comprehensive investigation of the biological impact of these stressors. However, detailed mechanistic networks are needed to determine which specific pathways are activated in response to different stressors and to drive the qualitative and eventually quantitative assessment of these data. A current limiting step in this process is the availability of detailed mechanistic networks that can be used as an analytical substrate. RESULTS: We have built a detailed network model that captures the biology underlying the physiological cellular response to endogenous and exogenous stressors in non-diseased mammalian pulmonary and cardiovascular cells. The contents of the network model reflect several diverse areas of signaling, including oxidative stress, hypoxia, shear stress, endoplasmic reticulum stress, and xenobiotic stress, that are elicited in response to common pulmonary and cardiovascular stressors. We then tested the ability of the network model to identify the mechanisms that are activated in response to CS, a broad inducer of cellular stress. Using transcriptomic data from the lungs of mice exposed to CS, the network model identified a robust increase in the oxidative stress response, largely mediated by the anti-oxidant NRF2 pathways, consistent with previous reports on the impact of CS exposure in the mammalian lung. CONCLUSIONS: The results presented here describe the construction of a cellular stress network model and its application towards the analysis of environmental stress using transcriptomic data. The proof-of-principle analysis described here, coupled with the future development of additional network models covering distinct areas of biology, will help to further clarify the integrated biological responses elicited by complex environmental stressors such as CS, in pulmonary and cardiovascular cells.

Construction of a computable cell proliferation network focused on non-diseased lung cells.

BACKGROUND: 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 (transcriptomics, proteomics/phosphoproteomics, metabolomics, etc.). 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 (COPD), and fibrosis. Unfortunately, no such network has been available prior to this work. RESULTS: 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 (Cell Cycle, Growth Factors, Cell Interaction, Intra- and Extracellular Signaling, and Epigenetics), and contains a total of 848 nodes (biological entities) and 1597 edges (relationships between biological entities). 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 (RB1, CDKN1A, and MYC/MYCN) 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. CONCLUSIONS: 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.

Centromere glue provides spindle cue.

During cell division, accurate distribution of the genome by the mitotic spindle requires that sister chromatids remain tethered until the right moment. Studies of an uncharacterized vertebrate protein, Sgo (Salic et al., 2004 [this issue of Cell]), reveal dual roles as a chromosome cohesion factor and a regulator of spindle microtubule dynamics.

The condensin complex is required for proper spindle assembly and chromosome segregation in Xenopus egg extracts.

Chromosome condensation is required for the physical resolution and segregation of sister chromatids during cell division, but the precise role of higher order chromatin structure in mitotic chromosome functions is unclear. Here, we address the role of the major condensation machinery, the condensin complex, in spindle assembly and function in Xenopus laevis egg extracts. Immunodepletion of condensin inhibited microtubule growth and organization around chromosomes, reducing the percentage of sperm nuclei capable of forming spindles, and causing dramatic defects in anaphase chromosome segregation. Although the motor CENP-E was recruited to kinetochores pulled poleward during anaphase, the disorganized chromosome mass was not resolved. Inhibition of condensin function during anaphase also inhibited chromosome segregation, indicating its continuous requirement. Spindle assembly around DNA-coated beads in the absence of kinetochores was also impaired upon condensin inhibition. These results support an important role for condensin in establishing chromosomal architecture necessary for proper spindle assembly and chromosome segregation.