scholarly journals Modeling a gene regulatory network of EMT hybrid states for mouse embryonic skin cells

2019 ◽  
Author(s):  
Dan Ramirez ◽  
Vivek Kohar ◽  
Ataur Katebi ◽  
Mingyang Lu

AbstractEpithelial-mesenchymal transition (EMT) plays a crucial role in embryonic development and tumorigenesis. Although EMT has been extensively studied with both computational and experimental methods, the gene regulatory mechanisms governing the transition are not yet well understood. Recent investigations have begun to better characterize the complex phenotypic plasticity underlying EMT using a computational systems biology approach. Here, we analyzed recently published single-cell RNA sequencing data from E9.5 to E11.5 mouse embryonic skin cells and identified the gene expression patterns of both epithelial and mesenchymal phenotypes, as well as a clear hybrid state. By integrating the scRNA-seq data and gene regulatory interactions from the literature, we constructed a gene regulatory network model governing the decision-making of EMT in the context of the developing mouse embryo. We simulated the network using a recently developed mathematical modeling method, named RACIPE, and observed three distinct phenotypic states whose gene expression patterns can be associated with the epithelial, hybrid, and mesenchymal states in the scRNA-seq data. Additionally, the model is in agreement with published results on the composition of EMT phenotypes and regulatory networks. We identified Wnt signaling as a major pathway in inducing the EMT and its role in driving cellular state transitions during embryonic development. Our findings demonstrate a new method of identifying and incorporating tissue-specific regulatory interactions into gene regulatory network modeling.Author SummaryEpithelial-mesenchymal transition (EMT) is a cellular process wherein cells become disconnected from their surroundings and acquire the ability to migrate through the body. EMT has been observed in biological contexts including development, wound healing, and cancer, yet the regulatory mechanisms underlying it are not well understood. Of particular interest is a purported hybrid state, in which cells can retain some adhesion to their surroundings but also show mesenchymal traits. Here, we examine the prevalence and composition of the hybrid state in the context of the embryonic mouse, integrating gene regulatory interactions from published experimental results as well as from the specific single cell RNA sequencing dataset of interest. Using mathematical modeling, we simulated a regulatory network based on these sources and aligned the simulated phenotypes with those in the data. We identified a hybrid EMT phenotype and revealed the inducing effect of Wnt signaling on EMT in this context. Our regulatory network construction process can be applied beyond EMT to illuminate the behavior of any biological phenomenon occurring in a specific context, allowing better identification of therapeutic targets and further research directions.

2021 ◽  
Author(s):  
Sreemol Gokuladhas ◽  
William Schierding ◽  
Roan Eltigani Zaied ◽  
Tayaza Fadason ◽  
Murim Choi ◽  
...  

Background & Aims: Non-alcoholic fatty liver disease (NAFLD) is a multi-system metabolic disease that co-occurs with various hepatic and extra-hepatic diseases. The phenotypic manifestation of NAFLD is primarily observed in the liver. Therefore, identifying liver-specific gene regulatory interactions between variants associated with NAFLD and multimorbid conditions may help to improve our understanding of underlying shared aetiology. Methods: Here, we constructed a liver-specific gene regulatory network (LGRN) consisting of genome-wide spatially constrained expression quantitative trait loci (eQTLs) and their target genes. The LGRN was used to identify regulatory interactions involving NAFLD-associated genetic modifiers and their inter-relationships to other complex traits. Results and Conclusions: We demonstrate that MBOAT7 and IL32, which are associated with NAFLD progression, are regulated by spatially constrained eQTLs that are enriched for an association with liver enzyme levels. MBOAT7 transcript levels are also linked to eQTLs associated with cirrhosis, and other traits that commonly co-occur with NAFLD. In addition, genes that encode interacting partners of NAFLD-candidate genes within the liver-specific protein-protein interaction network were affected by eQTLs enriched for phenotypes relevant to NAFLD (e.g. IgG glycosylation patterns, OSA). Furthermore, we identified distinct gene regulatory networks formed by the NAFLD-associated eQTLs in normal versus diseased liver, consistent with the context-specificity of the eQTLs effects. Interestingly, genes targeted by NAFLD-associated eQTLs within the LGRN were also affected by eQTLs associated with NAFLD-related traits (e.g. obesity and body fat percentage). Overall, the genetic links identified between these traits expand our understanding of shared regulatory mechanisms underlying NAFLD multimorbidities.


eLife ◽  
2016 ◽  
Vol 5 ◽  
Author(s):  
Erik Clark ◽  
Michael Akam

The Drosophila embryo transiently exhibits a double-segment periodicity, defined by the expression of seven 'pair-rule' genes, each in a pattern of seven stripes. At gastrulation, interactions between the pair-rule genes lead to frequency doubling and the patterning of 14 parasegment boundaries. In contrast to earlier stages of Drosophila anteroposterior patterning, this transition is not well understood. By carefully analysing the spatiotemporal dynamics of pair-rule gene expression, we demonstrate that frequency-doubling is precipitated by multiple coordinated changes to the network of regulatory interactions between the pair-rule genes. We identify the broadly expressed but temporally patterned transcription factor, Odd-paired (Opa/Zic), as the cause of these changes, and show that the patterning of the even-numbered parasegment boundaries relies on Opa-dependent regulatory interactions. Our findings indicate that the pair-rule gene regulatory network has a temporally modulated topology, permitting the pair-rule genes to play stage-specific patterning roles.


2019 ◽  
Author(s):  
Daniel Morgan ◽  
Matthew Studham ◽  
Andreas Tjärnberg ◽  
Holger Weishaupt ◽  
Fredrik J. Swartling ◽  
...  

AbstractThe gene regulatory network (GRN) of human cells encodes mechanisms to ensure proper functioning. However, if this GRN is dysregulated, the cell may enter into a disease state such as cancer. Understanding the GRN as a system can therefore help identify novel mechanisms underlying disease, which can lead to new therapies. Reliable inference of GRNs is however still a major challenge in systems biology.To deduce regulatory interactions relevant to cancer, we applied a recent computational inference framework to data from perturbation experiments in squamous carcinoma cell line A431. GRNs were inferred using several methods, and the false discovery rate was controlled by the NestBoot framework. We developed a novel approach to assess the predictiveness of inferred GRNs against validation data, despite the lack of a gold standard. The best GRN was significantly more predictive than the null model, both in crossvalidated benchmarks and for an independent dataset of the same genes under a different perturbation design. It agrees with many known links, in addition to predicting a large number of novel interactions from which a subset was experimentally validated. The inferred GRN captures regulatory interactions central to cancer-relevant processes and thus provides mechanistic insights that are useful for future cancer research.Data available at GSE125958Inferred GRNs and inference statistics available at https://dcolin.shinyapps.io/CancerGRN/ Software available at https://bitbucket.org/sonnhammergrni/genespider/src/BFECV/Author SummaryCancer is the second most common cause of death globally, and although cancer treatments have improved in recent years, we need to understand how regulatory mechanisms are altered in cancer to combat the disease efficiently. By applying gene perturbations and inference of gene regulatory networks to 40 genes known or suspected to have a role in cancer due to interactions with the oncogene MYC, we deduce their underlying regulatory interactions. Using a recent computational framework for inference together with a novel method for cross validation, we infer a reliable regulatory model of this system in a completely data driven manner, not reliant on literature or priors. The novel interactions add to the understanding of the progressive oncogenic regulatory process and may provide new targets for therapy.


2021 ◽  
Author(s):  
Juan Camilo Arboleda Rivera ◽  
Gloria Machado Rodriguez ◽  
Boris Anghelo Rodriguez Rey ◽  
Jayson Gutierrez Betancur

Background: A central problem in developmental and synthetic biology is understanding the mechanisms by which cells in a tissue or a Petri dish process external cues and transform such information into a coherent response, e.g., a terminal differentiation state. It was long believed that this type of positional information could be entirely attributed to a gradient of concentration of a specific signaling molecule (i.e., a morphogen). However, advances in experimental methodologies and computer modeling have demonstrated the crucial role of the dynamics of a cell's gene regulatory network (GRN) in decoding the information carried by the morphogen, which is eventually translated into a spatial pattern. This morphogen interpretation mechanism has gained much attention in systems biology as a tractable system to investigate the emergent properties of complex genotype-phenotype maps. Methods: In this study, we apply a Markov chain Monte Carlo (MCMC)-like algorithm to probe the design space of three-node GRNs with the ability to generate a band-like expression pattern (target phenotype) in the middle of an arrangement of 30 cells, which resemble a simple (1-D) morphogenetic field in a developing embryo. Unlike most modeling studies published so far, here we explore the space of GRN topologies with nodes having the potential to perceive the same input signal differently. This allows for a lot more flexibility during the search space process, and thus enables us to identify a larger set of potentially interesting and realizable morphogen interpretation mechanisms. Results: Out of 2061 GRNs selected using the search space algorithm, we found 714 classes of network topologies that could correctly interpret the morphogen. Notably, the main network motif that generated the target phenotype in response to the input signal was the type 3 Incoherent Feed-Forward Loop (I3-FFL), which agrees with previous theoretical expectations and experimental observations. Particularly, compared to a previously reported pattern forming GRN topologies, we have uncovered a great variety of novel network designs, some of which might be worth inquiring through synthetic biology methodologies to test for the ability of network design with minimal regulatory complexity to interpret a developmental cue robustly.


Oncogene ◽  
2021 ◽  
Author(s):  
Patrick Frey ◽  
Antoine Devisme ◽  
Katja Rose ◽  
Monika Schrempp ◽  
Vivien Freihen ◽  
...  

AbstractTransforming growth factor beta (TGFβ) superfamily signaling is a prime inducer of epithelial-mesenchymal transitions (EMT) that foster cancer cell invasion and metastasis, a major cause of cancer-related deaths. Yet, TGFβ signaling is frequently inactivated in human tumor entities including colorectal cancer (CRC) and pancreatic adenocarcinoma (PAAD) with a high proportion of mutations incapacitating SMAD4, which codes for a transcription factor (TF) central to canonical TGFβ and bone morphogenetic protein (BMP) signaling. Beyond its role in initiating EMT, SMAD4 was reported to crucially contribute to subsequent gene regulatory events during EMT execution. It is therefore widely assumed that SMAD4-mutant (SMAD4mut) cancer cells are unable to undergo EMT. Here, we scrutinized this notion and probed for potential SMAD4-independent EMT execution using SMAD4mut CRC cell lines. We show that SMAD4mut cells exhibit morphological changes, become invasive, and regulate EMT marker genes upon induction of the EMT-TF SNAIL1. Furthermore, SNAIL1-induced EMT in SMAD4mut cells was found to be entirely independent of TGFβ/BMP receptor activity. Global assessment of the SNAIL1-dependent transcriptome confirmed the manifestation of an EMT gene regulatory program in SMAD4mut cells highly related to established EMT signatures. Finally, analyses of human tumor transcriptomes showed that SMAD4 mutations are not underrepresented in mesenchymal tumor samples and that expression patterns of EMT-associated genes are similar in SMAD4mut and SMAD4 wild-type (SMAD4wt) cases. Altogether, our findings suggest that alternative TFs take over the gene regulatory functions of SMAD4 downstream of EMT-TFs, arguing for considerable plasticity of gene regulatory networks operating in EMT execution. Further, they establish that EMT is not categorically precluded in SMAD4mut tumors, which is relevant for their diagnostic and therapeutic evaluation.


2021 ◽  
Vol 12 ◽  
Author(s):  
Manasa Kalya ◽  
Alexander Kel ◽  
Darius Wlochowitz ◽  
Edgar Wingender ◽  
Tim Beißbarth

Only 2% of glioblastoma multiforme (GBM) patients respond to standard therapy and survive beyond 36 months (long-term survivors, LTS), while the majority survive less than 12 months (short-term survivors, STS). To understand the mechanism leading to poor survival, we analyzed publicly available datasets of 113 STS and 58 LTS. This analysis revealed 198 differentially expressed genes (DEGs) that characterize aggressive tumor growth and may be responsible for the poor prognosis. These genes belong largely to the Gene Ontology (GO) categories “epithelial-to-mesenchymal transition” and “response to hypoxia.” In this article, we applied an upstream analysis approach that involves state-of-the-art promoter analysis and network analysis of the dysregulated genes potentially responsible for short survival in GBM. Binding sites for transcription factors (TFs) associated with GBM pathology like NANOG, NF-κB, REST, FRA-1, PPARG, and seven others were found enriched in the promoters of the dysregulated genes. We reconstructed the gene regulatory network with several positive feedback loops controlled by five master regulators [insulin-like growth factor binding protein 2 (IGFBP2), vascular endothelial growth factor A (VEGFA), VEGF165, platelet-derived growth factor A (PDGFA), adipocyte enhancer-binding protein (AEBP1), and oncostatin M (OSMR)], which can be proposed as biomarkers and as therapeutic targets for enhancing GBM prognosis. A critical analysis of this gene regulatory network gives insights into the mechanism of gene regulation by IGFBP2 via several TFs including the key molecule of GBM tumor invasiveness and progression, FRA-1. All the observations were validated in independent cohorts, and their impact on overall survival has been investigated.


2019 ◽  
Vol 20 (1) ◽  
Author(s):  
Patrick M. Staunton ◽  
Aleksandra A. Miranda-CasoLuengo ◽  
Brendan J. Loftus ◽  
Isobel Claire Gormley

Abstract Background Although many of the genic features in Mycobacterium abscessus have been fully validated, a comprehensive understanding of the regulatory elements remains lacking. Moreover, there is little understanding of how the organism regulates its transcriptomic profile, enabling cells to survive in hostile environments. Here, to computationally infer the gene regulatory network for Mycobacterium abscessus we propose a novel statistical computational modelling approach: BayesIan gene regulatory Networks inferreD via gene coExpression and compaRative genomics (BINDER). In tandem with derived experimental coexpression data, the property of genomic conservation is exploited to probabilistically infer a gene regulatory network in Mycobacterium abscessus.Inference on regulatory interactions is conducted by combining ‘primary’ and ‘auxiliary’ data strata. The data forming the primary and auxiliary strata are derived from RNA-seq experiments and sequence information in the primary organism Mycobacterium abscessus as well as ChIP-seq data extracted from a related proxy organism Mycobacterium tuberculosis. The primary and auxiliary data are combined in a hierarchical Bayesian framework, informing the apposite bivariate likelihood function and prior distributions respectively. The inferred relationships provide insight to regulon groupings in Mycobacterium abscessus. Results We implement BINDER on data relating to a collection of 167,280 regulator-target pairs resulting in the identification of 54 regulator-target pairs, across 5 transcription factors, for which there is strong probability of regulatory interaction. Conclusions The inferred regulatory interactions provide insight to, and a valuable resource for further studies of, transcriptional control in Mycobacterium abscessus, and in the family of Mycobacteriaceae more generally. Further, the developed BINDER framework has broad applicability, useable in settings where computational inference of a gene regulatory network requires integration of data sources derived from both the primary organism of interest and from related proxy organisms.


Sign in / Sign up

Export Citation Format

Share Document