The Unfolded Protein Response Plays a Predominant Homeostatic Role in Response to Mitochondrial Stress in Pancreatic Stellate Cells

Activated pancreatic stellate cells (PaSC) are key participants in the stroma of pancreatic cancer, secreting extracellular matrix proteins and inflammatory mediators. Tumors are poorly vascularized, creating metabolic stress conditions in cancer and stromal cells that necessitate adaptive homeostatic cellular programs. Activation of autophagy and the endoplasmic reticulum unfolded protein response (UPR) have been described in hepatic stellate cells, but the role of these processes in PaSC responses to metabolic stress is unknown. We reported that the PI3K/mTOR pathway, which AMPK can regulate through multiple inputs, modulates PaSC activation and fibrogenic potential. Here, using primary and immortalized mouse PaSC, we assess the relative contributions of AMPK/mTOR signaling, autophagy and the UPR to cell fate responses during metabolic stress induced by mitochondrial dysfunction. The mitochondrial uncoupler rottlerin at low doses (0.5–2.5 μM) was added to cells cultured in 10% FBS complete media. Mitochondria rapidly depolarized, followed by altered mitochondrial dynamics and decreased cellular ATP levels. This mitochondrial dysfunction elicited rapid, sustained AMPK activation, mTOR pathway inhibition, and blockade of autophagic flux. Rottlerin treatment also induced rapid, sustained PERK/CHOP UPR signaling. Subsequently, high doses (>5 μM) induced loss of cell viability and cell death. Interestingly, AMPK knock-down using siRNA did not prevent rottlerin-induced mTOR inhibition, autophagy, or CHOP upregulation, suggesting that AMPK is dispensable for these responses. Moreover, CHOP genetic deletion, but not AMPK knock-down, prevented rottlerin-induced apoptosis and supported cell survival, suggesting that UPR signaling is a major modulator of cell fate in PaSC during metabolic stress. Further, short-term rottlerin treatment reduced both PaSC fibrogenic potential and IL-6 mRNA expression. In contrast, expression levels of the angiogenic factors HGF and VEGFα were unaffected, and the immune modulator IL-4 was markedly upregulated. These data imply that metabolic stress-induced PaSC reprogramming differentially modulates neighboring cells in the tumor microenvironment.     

A Factor Graph Nested Effects Model To Identify Networks from Genetic Perturbations

Complex phenotypes such as the transformation of a normal population of cells into cancerous tissue result from a series of molecular triggers gone awry. We describe a method that searches for a genetic network consistent with expression changes observed under the knock-down of a set of genes that share a common role in the cell, such as a disease phenotype. The method extends the Nested Effects Model of Markowetz et al. (2005) by using a probabilistic factor graph to search for a network representing interactions among these silenced genes. The method also expands the network by attaching new genes at specific downstream points, providing candidates for subsequent perturbations to further characterize the pathway. We investigated an extension provided by the factor graph approach in which the model distinguishes between inhibitory and stimulatory interactions. We found that the extension yielded significant improvements in recovering the structure of simulated and Saccharomyces cerevisae networks. We applied the approach to discover a signaling network among genes involved in a human colon cancer cell invasiveness pathway. The method predicts several genes with new roles in the invasiveness process. We knocked down two genes identified by our approach and found that both knock-downs produce loss of invasive potential in a colon cancer cell line. Nested effects models may be a powerful tool for inferring regulatory connections and genes that operate in normal and disease-related processes.     

Criticality Is an Emergent Property of Genetic Networks that Exhibit Evolvability

Accumulating experimental evidence suggests that the gene regulatory networks of living organisms operate in the critical phase, namely, at the transition between ordered and chaotic dynamics. Such critical dynamics of the network permits the coexistence of robustness and flexibility which are necessary to ensure homeostatic stability (of a given phenotype) while allowing for switching between multiple phenotypes (network states) as occurs in development and in response to environmental change. However, the mechanisms through which genetic networks evolve such critical behavior have remained elusive. Here we present an evolutionary model in which criticality naturally emerges from the need to balance between the two essential components of evolvability: phenotype conservation and phenotype innovation under mutations. We simulated the Darwinian evolution of random Boolean networks that mutate gene regulatory interactions and grow by gene duplication. The mutating networks were subjected to selection for networks that both (i) preserve all the already acquired phenotypes (dynamical attractor states) and (ii) generate new ones. Our results show that this interplay between extending the phenotypic landscape (innovation) while conserving the existing phenotypes (conservation) suffices to cause the evolution of all the networks in a population towards criticality. Furthermore, the networks produced by this evolutionary process exhibit structures with hubs (global regulators) similar to the observed topology of real gene regulatory networks. Thus, dynamical criticality and certain elementary topological properties of gene regulatory networks can emerge as a byproduct of the evolvability of the phenotypic landscape.     

Homeostatic Interplay between Bacterial Cell-Cell Signaling and Iron in Virulence

Pathogenic bacteria use interconnected multi-layered regulatory networks, such as quorum sensing (QS) networks to sense and respond to environmental cues and external and internal bacterial cell signals, and thereby adapt to and exploit target hosts. Despite the many advances that have been made in understanding QS regulation, little is known regarding how these inputs are integrated and processed in the context of multi-layered QS regulatory networks. Here we report the examination of the Pseudomonas aeruginosa QS 4-hydroxy-2-alkylquinolines (HAQs) MvfR regulatory network and determination of its interaction with the QS acyl-homoserine-lactone (AHL) RhlR network. The aim of this work was to elucidate paradigmatically the complex relationships between multi-layered regulatory QS circuitries, their signaling molecules, and the environmental cues to which they respond. Our findings revealed positive and negative homeostatic regulatory loops that fine-tune the MvfR regulon via a multi-layered dependent homeostatic regulation of the cell-cell signaling molecules PQS and HHQ, and interplay between these molecules and iron. We discovered that the MvfR regulon component PqsE is a key mediator in orchestrating this homeostatic regulation, and in establishing a connection to the QS rhlR system in cooperation with RhlR. Our results show that P. aeruginosa modulates the intensity of its virulence response, at least in part, through this multi-layered interplay. Our findings underscore the importance of the homeostatic interplay that balances competition within and between QS systems via cell-cell signaling molecules and environmental cues in the control of virulence gene expression. Elucidation of the fine-tuning of this complex relationship offers novel insights into the regulation of these systems and may inform strategies designed to limit infections caused by P. aeruginosa and related human pathogens.     

Identification of the angiogenesis signaling domain in pleiotrophin defines a mechanism of the angiogenic switch

Neoplasms progress through genetic and epigenetic mutations that deregulate pathways in the malignant cell that stimulate more aggressive growth of the malignant cell itself and/or remodel the tumor microenvironment to support the developing tumor mass. The appearance of new blood vessels in malignant tumors is known as the "angiogenic switch." The angiogenic switch triggers a stage of rapid tumor growth supported by extensive tumor angiogenesis and a more aggressive tumor phenotype and its onset is a poor prognostic indicator for host survival. Identification of the factors that stimulate the angiogenic switch thus is of high importance. Pleiotrophin (PTN the protein, Ptn the gene) is an angiogenic factor and the Ptn gene has been found to be constitutively expressed in many human tumors of different cell types. These studies use a nude mouse model to test if Ptn constitutively expressed in premalignant cells is sufficient to trigger an angiogenic switch in vivo. We introduced an ectopic Ptn gene into "premalignant" SW-13 cells and analyzed the phenotype of SW-13 Ptn cell tumor implants in the flanks of nude mice. SW-13 Ptn cell subcutaneous tumor implants grew very rapidly and had a striking increase in the density of new blood vessels compared to the SW-13 cell tumor implants, suggesting that constitutive PTN signaling in the premalignant SW-13 cell implants in the nude mouse recapitulates fully the angiogenic switch. It was found also that ectopic expression of the C-terminal domain of PTN in SW-13 cell implants was equally effective in initiating an angiogenic switch as the full-length PTN whereas implants of SW-13 cells in nude mice that express the N-terminal domain of PTN grew rapidly but failed to develop tumor angiogenesis. The data suggest the possibility that mutations that activate Ptn in premalignant cells are sufficient to stimulate an angiogenic switch in vivo and, since these mutations are frequently found in human malignancies, that constitutive PTN signaling may be an important contributor to progression of human tumors. The data also suggest that the C-terminal and the N-terminal domains of PTN equally initiate switches in premalignant cells to cells of a more aggressive tumor phenotype but the separate domains of PTN signal different mechanisms and perhaps signal through activation of a separate receptor-like protein.     

Modeling the proteome of a Marek's disease transformed cell line: a natural animal model for CD30 overexpressing lymphomas

Marek's disease (MD) in the chicken, caused by the highly infectious MD alpha-herpesvirus (MDV), is both commercially important and a unique, naturally occurring model for human T-cell lymphomas overexpressing the Hodgkin's disease antigen, CD30. Here, we used proteomics as a basis for modeling the molecular functions and biological processes involved in MDV-induced lymphomagenesis. Proteins were extracted from an MDV-transformed cell line and were then identified using 2-D LC-ESI-MS/MS. From the resulting 3870 cellular and 21 MDV proteins we confirm the existence of 3150 "predicted" and 12 "hypothetical" chicken proteins. The UA-01 proteome is proliferative, differentiated, angiogenic, pro-metastatic and pro-immune-escape but anti-programmed cell death, -anergy, -quiescence and -senescence and is consistent with a cancer phenotype. In particular, the pro-metastatic integrin signaling pathway and the ERK/MAPK signaling pathways were the two predominant signaling pathways represented. The cytokines, cytokine receptors, and their related proteins suggest that UA-01 has a regulatory T-cell phenotype.     

Glucocorticoid signaling in pancreatic islets modulates gene regulatory programs and genetic risk of type 2 diabetes

Glucocorticoids are key regulators of glucose homeostasis and pancreatic islet function, but the gene regulatory programs driving responses to glucocorticoid signaling in islets and the contribution of these programs to diabetes risk are unknown. In this study we used ATAC-seq and RNA-seq to map chromatin accessibility and gene expression from eleven primary human islet samples cultured in vitro with the glucocorticoid dexamethasone at multiple doses and durations. We identified thousands of accessible chromatin sites and genes with significant changes in activity in response to glucocorticoids. Chromatin sites up-regulated in glucocorticoid signaling were prominently enriched for glucocorticoid receptor binding sites and up-regulated genes were enriched for ion transport and lipid metabolism, whereas down-regulated chromatin sites and genes were enriched for inflammatory, stress response and proliferative processes. Genetic variants associated with glucose levels and T2D risk were enriched in glucocorticoid-responsive chromatin sites, including fine-mapped variants at 51 known signals. Among fine-mapped variants in glucocorticoid-responsive chromatin, a likely casual variant at the 2p21 locus had glucocorticoid-dependent allelic effects on beta cell enhancer activity and affected SIX2 and SIX3 expression. Our results provide a comprehensive map of islet regulatory programs in response to glucocorticoids through which we uncover a role for islet glucocorticoid signaling in mediating genetic risk of T2D.     

Prediction of Oncogenic Interactions and Cancer-Related Signaling Networks Based on Network Topology

Cancer has been increasingly recognized as a systems biology disease since many investigators have demonstrated that this malignant phenotype emerges from abnormal protein-protein, regulatory and metabolic interactions induced by simultaneous structural and regulatory changes in multiple genes and pathways. Therefore, the identification of oncogenic interactions and cancer-related signaling networks is crucial for better understanding cancer. As experimental techniques for determining such interactions and signaling networks are labor-intensive and time-consuming, the development of a computational approach capable to accomplish this task would be of great value. For this purpose, we present here a novel computational approach based on network topology and machine learning capable to predict oncogenic interactions and extract relevant cancer-related signaling subnetworks from an integrated network of human genes interactions (INHGI). This approach, called graph2sig, is twofold: first, it assigns oncogenic scores to all interactions in the INHGI and then these oncogenic scores are used as edge weights to extract oncogenic signaling subnetworks from INHGI. Regarding the prediction of oncogenic interactions, we showed that graph2sig is able to recover 89% of known oncogenic interactions with a precision of 77%. Moreover, the interactions that received high oncogenic scores are enriched in genes for which mutations have been causally implicated in cancer. We also demonstrated that graph2sig is potentially useful in extracting oncogenic signaling subnetworks: more than 80% of constructed subnetworks contain more than 50% of original interactions in their corresponding oncogenic linear pathways present in the KEGG PATHWAY database. In addition, the potential oncogenic signaling subnetworks discovered by graph2sig are supported by experimental evidence. Taken together, these results suggest that graph2sig can be a useful tool for investigators involved in cancer research interested in detecting signaling networks most prone to contribute with the emergence of malignant phenotype.     

Enteric commensal bacteria induce extracellular signal-regulated kinase pathway signaling via formyl peptide receptor-dependent redox modulation of dual specific phosphatase 3

The normal microbial occupants of the mammalian intestine are crucial for maintaining gut homeostasis, yet the mechanisms by which intestinal cells perceive and respond to the microbiota are largely unknown. Intestinal epithelial contact with commensal bacteria and/or their products has been shown to activate noninflammatory signaling pathways, such as extracellular signal-related kinase (ERK), thus influencing homeostatic processes. We previously demonstrated that commensal bacteria stimulate ERK pathway activity via interaction with formyl peptide receptors (FPRs). In the current study, we expand on these findings and show that commensal bacteria initiate ERK signaling through rapid FPR-dependent reactive oxygen species (ROS) generation and subsequent modulation of MAP kinase phosphatase redox status. ROS generation induced by the commensal bacteria Lactobacillus rhamnosus GG and the FPR peptide ligand, N-formyl-Met-Leu-Phe, was abolished in the presence of selective inhibitors for G protein-coupled signaling and FPR ligand interaction. In addition, pretreatment of cells with inhibitors of ROS generation attenuated commensal bacteria-induced ERK signaling, indicating that ROS generation is required for ERK pathway activation. Bacterial colonization also led to oxidative inactivation of the redox-sensitive and ERK-specific phosphatase, DUSP3/VHR, and consequent stimulation of ERK pathway signaling. Together, these data demonstrate that commensal bacteria and their products activate ROS signaling in an FPR-dependent manner and define a mechanism by which cellular ROS influences the ERK pathway through a redox-sensitive regulatory circuit.