1
|
Lin S, Lim B. Multifaceted effects on even-skipped transcriptional dynamics upon Krüppel dosage changes. Development 2024; 151:dev202132. [PMID: 38345298 PMCID: PMC10948998 DOI: 10.1242/dev.202132] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/27/2023] [Accepted: 02/08/2024] [Indexed: 03/05/2024]
Abstract
Although fluctuations in transcription factor (TF) dosage are often well tolerated, TF dosage modulation can change the target gene expression dynamics and result in significant non-lethal developmental phenotypes. Using MS2/MCP-mediated quantitative live imaging in early Drosophila embryos, we analyzed how changing levels of the gap gene Krüppel (Kr) affects transcriptional dynamics of the pair-rule gene even-skipped (eve). Halving the Kr dosage leads to a transient posterior expansion of the eve stripe 2 and an anterior shift of stripe 5. Surprisingly, the most significant changes are observed in eve stripes 3 and 4, the enhancers of which do not contain Kr-binding sites. In Kr heterozygous embryos, both stripes 3 and 4 display narrower widths, anteriorly shifted boundaries and reduced mRNA production levels. We show that Kr dosage indirectly affects stripe 3 and 4 dynamics by modulating other gap gene dynamics. We quantitatively correlate moderate body segment phenotypes of Kr heterozygotes with spatiotemporal changes in eve expression. Our results indicate that nonlinear relationships between TF dosage and phenotypes underlie direct TF-DNA and indirect TF-TF interactions.
Collapse
Affiliation(s)
- Shufan Lin
- Department of Bioengineering, University of Pennsylvania, Philadelphia, PA 19104, USA
| | - Bomyi Lim
- Department of Chemical and Biomolecular Engineering, University of Pennsylvania, Philadelphia, PA 19104, USA
| |
Collapse
|
2
|
Patra P, B R D, Kundu P, Das M, Ghosh A. Recent advances in machine learning applications in metabolic engineering. Biotechnol Adv 2023; 62:108069. [PMID: 36442697 DOI: 10.1016/j.biotechadv.2022.108069] [Citation(s) in RCA: 9] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/17/2022] [Revised: 10/18/2022] [Accepted: 11/22/2022] [Indexed: 11/27/2022]
Abstract
Metabolic engineering encompasses several widely-used strategies, which currently hold a high seat in the field of biotechnology when its potential is manifesting through a plethora of research and commercial products with a strong societal impact. The genomic revolution that occurred almost three decades ago has initiated the generation of large omics-datasets which has helped in gaining a better understanding of cellular behavior. The itinerary of metabolic engineering that has occurred based on these large datasets has allowed researchers to gain detailed insights and a reasonable understanding of the intricacies of biosystems. However, the existing trail-and-error approaches for metabolic engineering are laborious and time-intensive when it comes to the production of target compounds with high yields through genetic manipulations in host organisms. Machine learning (ML) coupled with the available metabolic engineering test instances and omics data brings a comprehensive and multidisciplinary approach that enables scientists to evaluate various parameters for effective strain design. This vast amount of biological data should be standardized through knowledge engineering to train different ML models for providing accurate predictions in gene circuits designing, modification of proteins, optimization of bioprocess parameters for scaling up, and screening of hyper-producing robust cell factories. This review briefs on the premise of ML, followed by mentioning various ML methods and algorithms alongside the numerous omics datasets available to train ML models for predicting metabolic outcomes with high-accuracy. The combinative interplay between the ML algorithms and biological datasets through knowledge engineering have guided the recent advancements in applications such as CRISPR/Cas systems, gene circuits, protein engineering, metabolic pathway reconstruction, and bioprocess engineering. Finally, this review addresses the probable challenges of applying ML in metabolic engineering which will guide the researchers toward novel techniques to overcome the limitations.
Collapse
Affiliation(s)
- Pradipta Patra
- School School of Energy Science and Engineering, Indian Institute of Technology Kharagpur, West Bengal 721302, India
| | - Disha B R
- B.M.S College of Engineering, Basavanagudi, Bengaluru, Karnataka 560019, India
| | - Pritam Kundu
- School School of Energy Science and Engineering, Indian Institute of Technology Kharagpur, West Bengal 721302, India
| | - Manali Das
- School of Bioscience, Indian Institute of Technology Kharagpur, West Bengal 721302, India
| | - Amit Ghosh
- School School of Energy Science and Engineering, Indian Institute of Technology Kharagpur, West Bengal 721302, India; P.K. Sinha Centre for Bioenergy and Renewables, Indian Institute of Technology Kharagpur, West Bengal 721302, India.
| |
Collapse
|
3
|
Duk MA, Gursky VV, Samsonova MG, Surkova SY. Application of Domain- and Genotype-Specific Models to Infer Post-Transcriptional Regulation of Segmentation Gene Expression in Drosophila. Life (Basel) 2021; 11:life11111232. [PMID: 34833107 PMCID: PMC8618293 DOI: 10.3390/life11111232] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/24/2021] [Revised: 11/05/2021] [Accepted: 11/10/2021] [Indexed: 11/16/2022] Open
Abstract
Unlike transcriptional regulation, the post-transcriptional mechanisms underlying zygotic segmentation gene expression in early Drosophila embryo have been insufficiently investigated. Condition-specific post-transcriptional regulation plays an important role in the development of many organisms. Our recent study revealed the domain- and genotype-specific differences between mRNA and the protein expression of Drosophila hb, gt, and eve genes in cleavage cycle 14A. Here, we use this dataset and the dynamic mathematical model to recapitulate protein expression from the corresponding mRNA patterns. The condition-specific nonuniformity in parameter values is further interpreted in terms of possible post-transcriptional modifications. For hb expression in wild-type embryos, our results predict the position-specific differences in protein production. The protein synthesis rate parameter is significantly higher in hb anterior domain compared to the posterior domain. The parameter sets describing Gt protein dynamics in wild-type embryos and Kr mutants are genotype-specific. The spatial discrepancy between gt mRNA and protein posterior expression in Kr mutants is well reproduced by the whole axis model, thus rejecting the involvement of post-transcriptional mechanisms. Our models fail to describe the full dynamics of eve expression, presumably due to its complex shape and the variable time delays between mRNA and protein patterns, which likely require a more complex model. Overall, our modeling approach enables the prediction of regulatory scenarios underlying the condition-specific differences between mRNA and protein expression in early embryo.
Collapse
Affiliation(s)
- Maria A. Duk
- Mathematical Biology and Bioinformatics Laboratory, Peter the Great Saint Petersburg Polytechnic University, 195251 St. Petersburg, Russia; (M.A.D.); (M.G.S.)
- Theoretical Department, Ioffe Institute, 194021 St. Petersburg, Russia;
| | - Vitaly V. Gursky
- Theoretical Department, Ioffe Institute, 194021 St. Petersburg, Russia;
| | - Maria G. Samsonova
- Mathematical Biology and Bioinformatics Laboratory, Peter the Great Saint Petersburg Polytechnic University, 195251 St. Petersburg, Russia; (M.A.D.); (M.G.S.)
| | - Svetlana Yu. Surkova
- Mathematical Biology and Bioinformatics Laboratory, Peter the Great Saint Petersburg Polytechnic University, 195251 St. Petersburg, Russia; (M.A.D.); (M.G.S.)
- Correspondence:
| |
Collapse
|
4
|
Jaeger J, Monk N. Dynamical modules in metabolism, cell and developmental biology. Interface Focus 2021; 11:20210011. [PMID: 34055307 PMCID: PMC8086940 DOI: 10.1098/rsfs.2021.0011] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 03/10/2021] [Indexed: 02/06/2023] Open
Abstract
Modularity is an essential feature of any adaptive complex system. Phenotypic traits are modules in the sense that they have a distinguishable structure or function, which can vary (quasi-)independently from its context. Since all phenotypic traits are the product of some underlying regulatory dynamics, the generative processes that constitute the genotype-phenotype map must also be functionally modular. Traditionally, modular processes have been identified as structural modules in regulatory networks. However, structure only constrains, but does not determine, the dynamics of a process. Here, we propose an alternative approach that decomposes the behaviour of a complex regulatory system into elementary activity-functions. Modular activities can occur in networks that show no structural modularity, making dynamical modularity more widely applicable than structural decomposition. Furthermore, the behaviour of a regulatory system closely mirrors its functional contribution to the outcome of a process, which makes dynamical modularity particularly suited for functional decomposition. We illustrate our approach with numerous examples from the study of metabolism, cellular processes, as well as development and pattern formation. We argue that dynamical modules provide a shared conceptual foundation for developmental and evolutionary biology, and serve as the foundation for a new account of process homology, which is presented in a separate contribution by DiFrisco and Jaeger to this focus issue.
Collapse
Affiliation(s)
- Johannes Jaeger
- Complexity Science Hub (CSH) Vienna, Josefstädter Strasse 39, 1080 Vienna, Austria
| | - Nick Monk
- School of Mathematics and Statistics, University of Sheffield, Hicks Building, Sheffield S3 7RH, UK
| |
Collapse
|
5
|
Volk MJ, Lourentzou I, Mishra S, Vo LT, Zhai C, Zhao H. Biosystems Design by Machine Learning. ACS Synth Biol 2020; 9:1514-1533. [PMID: 32485108 DOI: 10.1021/acssynbio.0c00129] [Citation(s) in RCA: 52] [Impact Index Per Article: 13.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
Biosystems such as enzymes, pathways, and whole cells have been increasingly explored for biotechnological applications. However, the intricate connectivity and resulting complexity of biosystems poses a major hurdle in designing biosystems with desirable features. As -omics and other high throughput technologies have been rapidly developed, the promise of applying machine learning (ML) techniques in biosystems design has started to become a reality. ML models enable the identification of patterns within complicated biological data across multiple scales of analysis and can augment biosystems design applications by predicting new candidates for optimized performance. ML is being used at every stage of biosystems design to help find nonobvious engineering solutions with fewer design iterations. In this review, we first describe commonly used models and modeling paradigms within ML. We then discuss some applications of these models that have already shown success in biotechnological applications. Moreover, we discuss successful applications at all scales of biosystems design, including nucleic acids, genetic circuits, proteins, pathways, genomes, and bioprocesses. Finally, we discuss some limitations of these methods and potential solutions as well as prospects of the combination of ML and biosystems design.
Collapse
|
6
|
Wu S, Cui T, Zhang X, Tian T. A non-linear reverse-engineering method for inferring genetic regulatory networks. PeerJ 2020; 8:e9065. [PMID: 32391205 PMCID: PMC7195839 DOI: 10.7717/peerj.9065] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/06/2020] [Accepted: 04/05/2020] [Indexed: 12/19/2022] Open
Abstract
Hematopoiesis is a highly complex developmental process that produces various types of blood cells. This process is regulated by different genetic networks that control the proliferation, differentiation, and maturation of hematopoietic stem cells (HSCs). Although substantial progress has been made for understanding hematopoiesis, the detailed regulatory mechanisms for the fate determination of HSCs are still unraveled. In this study, we propose a novel approach to infer the detailed regulatory mechanisms. This work is designed to develop a mathematical framework that is able to realize nonlinear gene expression dynamics accurately. In particular, we intended to investigate the effect of possible protein heterodimers and/or synergistic effect in genetic regulation. This approach includes the Extended Forward Search Algorithm to infer network structure (top-down approach) and a non-linear mathematical model to infer dynamical property (bottom-up approach). Based on the published experimental data, we study two regulatory networks of 11 genes for regulating the erythrocyte differentiation pathway and the neutrophil differentiation pathway. The proposed algorithm is first applied to predict the network topologies among 11 genes and 55 non-linear terms which may be for heterodimers and/or synergistic effect. Then, the unknown model parameters are estimated by fitting simulations to the expression data of two different differentiation pathways. In addition, the edge deletion test is conducted to remove possible insignificant regulations from the inferred networks. Furthermore, the robustness property of the mathematical model is employed as an additional criterion to choose better network reconstruction results. Our simulation results successfully realized experimental data for two different differentiation pathways, which suggests that the proposed approach is an effective method to infer the topological structure and dynamic property of genetic regulations.
Collapse
Affiliation(s)
- Siyuan Wu
- School of Mathematics, Monash University, Clayton, VIC, Australia
| | - Tiangang Cui
- School of Mathematics, Monash University, Clayton, VIC, Australia
| | - Xinan Zhang
- School of Mathematics and Statistics, Central China Normal University, Wuhan, PR China
| | - Tianhai Tian
- School of Mathematics, Monash University, Clayton, VIC, Australia
| |
Collapse
|
7
|
Roy J, Cheung E, Bhatti J, Muneem A, Lobo D. Curation and annotation of planarian gene expression patterns with segmented reference morphologies. Bioinformatics 2020; 36:2881-2887. [DOI: 10.1093/bioinformatics/btaa023] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/18/2019] [Revised: 12/07/2019] [Accepted: 01/14/2020] [Indexed: 12/30/2022] Open
Abstract
Abstract
Motivation
Morphological and genetic spatial data from functional experiments based on genetic, surgical and pharmacological perturbations are being produced at an extraordinary pace in developmental and regenerative biology. However, our ability to extract knowledge from these large datasets are hindered due to the lack of formalization methods and tools able to unambiguously describe, centralize and interpret them. Formalizing spatial phenotypes and gene expression patterns is especially challenging in organisms with highly variable morphologies such as planarian worms, which due to their extraordinary regenerative capability can experimentally result in phenotypes with almost any combination of body regions or parts.
Results
Here, we present a computational methodology and mathematical formalism to encode and curate the morphological outcomes and gene expression patterns in planaria. Worm morphologies are encoded with mathematical graphs based on anatomical ontology terms to automatically generate reference morphologies. Gene expression patterns are registered to these standard reference morphologies, which can then be annotated automatically with anatomical ontology terms by analyzing the spatial expression patterns and their textual descriptions. This methodology enables the curation and annotation of complex experimental morphologies together with their gene expression patterns in a centralized standardized dataset, paving the way for the extraction of knowledge and reverse-engineering of the much sought-after mechanistic models in planaria and other regenerative organisms.
Availability and implementation
We implemented this methodology in a user-friendly graphical software tool, PlanGexQ, freely available together with the data in the manuscript at https://lobolab.umbc.edu/plangexq.
Supplementary information
Supplementary data are available at Bioinformatics online.
Collapse
Affiliation(s)
- Joy Roy
- Department of Biological Sciences, University of Maryland, Baltimore County, Baltimore, MD 21250, USA
| | - Eric Cheung
- Department of Biological Sciences, University of Maryland, Baltimore County, Baltimore, MD 21250, USA
| | - Junaid Bhatti
- Department of Biological Sciences, University of Maryland, Baltimore County, Baltimore, MD 21250, USA
| | - Abraar Muneem
- Department of Biological Sciences, University of Maryland, Baltimore County, Baltimore, MD 21250, USA
| | - Daniel Lobo
- Department of Biological Sciences, University of Maryland, Baltimore County, Baltimore, MD 21250, USA
| |
Collapse
|
8
|
Jaeger J, Verd B. Dynamic positional information: Patterning mechanism versus precision in gradient-driven systems. Curr Top Dev Biol 2019; 137:219-246. [PMID: 32143744 DOI: 10.1016/bs.ctdb.2019.11.017] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Abstract
There is much talk about information in biology. In developmental biology, this takes the form of "positional information," especially in the context of morphogen-based pattern formation. Unfortunately, the concept of "information" is rarely defined in any precise manner. Here, we provide two alternative interpretations of "positional information," and examine the complementary meanings and uses of each concept. Positional information defined as Shannon information helps us understand decoding and error propagation in patterning systems. General relativistic positional information, in contrast, provides a metric to assess the output of pattern-forming mechanisms. Both interpretations provide powerful conceptual tools that do not compete, but are best used in combination to gain a proper mechanistic understanding of robust patterning.
Collapse
Affiliation(s)
- Johannes Jaeger
- Complexity Science Hub (CSH), Vienna, Austria; Department of Molecular Evolution & Development, University of Vienna, Vienna, Austria.
| | - Berta Verd
- Department of Genetics, University of Cambridge, Cambridge, United Kingdom
| |
Collapse
|
9
|
Classification-Based Inference of Dynamical Models of Gene Regulatory Networks. G3-GENES GENOMES GENETICS 2019; 9:4183-4195. [PMID: 31624138 PMCID: PMC6893186 DOI: 10.1534/g3.119.400603] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 12/05/2022]
Abstract
Cell-fate decisions during development are controlled by densely interconnected gene regulatory networks (GRNs) consisting of many genes. Inferring and predictively modeling these GRNs is crucial for understanding development and other physiological processes. Gene circuits, coupled differential equations that represent gene product synthesis with a switch-like function, provide a biologically realistic framework for modeling the time evolution of gene expression. However, their use has been limited to smaller networks due to the computational expense of inferring model parameters from gene expression data using global non-linear optimization. Here we show that the switch-like nature of gene regulation can be exploited to break the gene circuit inference problem into two simpler optimization problems that are amenable to computationally efficient supervised learning techniques. We present FIGR (Fast Inference of Gene Regulation), a novel classification-based inference approach to determining gene circuit parameters. We demonstrate FIGR’s effectiveness on synthetic data generated from random gene circuits of up to 50 genes as well as experimental data from the gap gene system of Drosophila melanogaster, a benchmark for inferring dynamical GRN models. FIGR is faster than global non-linear optimization by a factor of 600 and its computational complexity scales much better with GRN size. On a practical level, FIGR can accurately infer the biologically realistic gap gene network in under a minute on desktop-class hardware instead of requiring hours of parallel computing. We anticipate that FIGR would enable the inference of much larger biologically realistic GRNs than was possible before.
Collapse
|
10
|
Verd B, Monk NAM, Jaeger J. Modularity, criticality, and evolvability of a developmental gene regulatory network. eLife 2019; 8:e42832. [PMID: 31169494 PMCID: PMC6645726 DOI: 10.7554/elife.42832] [Citation(s) in RCA: 47] [Impact Index Per Article: 9.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/14/2019] [Accepted: 06/05/2019] [Indexed: 01/16/2023] Open
Abstract
The existence of discrete phenotypic traits suggests that the complex regulatory processes which produce them are functionally modular. These processes are usually represented by networks. Only modular networks can be partitioned into intelligible subcircuits able to evolve relatively independently. Traditionally, functional modularity is approximated by detection of modularity in network structure. However, the correlation between structure and function is loose. Many regulatory networks exhibit modular behaviour without structural modularity. Here we partition an experimentally tractable regulatory network-the gap gene system of dipteran insects-using an alternative approach. We show that this system, although not structurally modular, is composed of dynamical modules driving different aspects of whole-network behaviour. All these subcircuits share the same regulatory structure, but differ in components and sensitivity to regulatory interactions. Some subcircuits are in a state of criticality, while others are not, which explains the observed differential evolvability of the various expression features in the system.
Collapse
Affiliation(s)
- Berta Verd
- EMBL/CRG Systems Biology Research Unit, Centre for Genomic Regulation (CRG)The Barcelona Institute of Science and TechnologyBarcelonaSpain
- Universitat Pompeu Fabra (UPF)BarcelonaSpain
- Konrad Lorenz Institute for Evolution and Cognition Research (KLI)KlosterneuburgAustria
- Department of GeneticsUniversity of CambridgeCambridgeUnited Kingdom
| | - Nicholas AM Monk
- School of Mathematics and StatisticsUniversity of SheffieldSheffieldUnited States
| | - Johannes Jaeger
- EMBL/CRG Systems Biology Research Unit, Centre for Genomic Regulation (CRG)The Barcelona Institute of Science and TechnologyBarcelonaSpain
- Universitat Pompeu Fabra (UPF)BarcelonaSpain
- Konrad Lorenz Institute for Evolution and Cognition Research (KLI)KlosterneuburgAustria
- School of Mathematics and StatisticsUniversity of SheffieldSheffieldUnited States
- Wissenschaftskolleg zu BerlinBerlinGermany
- Center for Systems Biology Dresden (CSBD)DresdenGermany
- Complexity Science Hub (CSH)ViennaAustria
- Centre de Recherches Interdisciplinaires (CRI)ParisFrance
| |
Collapse
|
11
|
Hiscock TW. Adapting machine-learning algorithms to design gene circuits. BMC Bioinformatics 2019; 20:214. [PMID: 31029103 PMCID: PMC6487017 DOI: 10.1186/s12859-019-2788-3] [Citation(s) in RCA: 16] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/18/2019] [Accepted: 04/02/2019] [Indexed: 02/07/2023] Open
Abstract
BACKGROUND Gene circuits are important in many aspects of biology, and perform a wide variety of different functions. For example, some circuits oscillate (e.g. the cell cycle), some are bistable (e.g. as cells differentiate), some respond sharply to environmental signals (e.g. ultrasensitivity), and some pattern multicellular tissues (e.g. Turing's model). Often, one starts from a given circuit, and using simulations, asks what functions it can perform. Here we want to do the opposite: starting from a prescribed function, can we find a circuit that executes this function? Whilst simple in principle, this task is challenging from a computational perspective, since gene circuit models are complex systems with many parameters. In this work, we adapted machine-learning algorithms to significantly accelerate gene circuit discovery. RESULTS We use gradient-descent optimization algorithms from machine learning to rapidly screen and design gene circuits. With this approach, we found that we could rapidly design circuits capable of executing a range of different functions, including those that: (1) recapitulate important in vivo phenomena, such as oscillators, and (2) perform complex tasks for synthetic biology, such as counting noisy biological events. CONCLUSIONS Our computational pipeline will facilitate the systematic study of natural circuits in a range of contexts, and allow the automatic design of circuits for synthetic biology. Our method can be readily applied to biological networks of any type and size, and is provided as an open-source and easy-to-use python module, GeneNet.
Collapse
Affiliation(s)
- Tom W Hiscock
- Cancer Research UK, Cambridge Institute, Li Ka Shing Centre, Robinson Way, Cambridge, CB2 0RE, UK.
- Wellcome Trust/Cancer Research UK Gurdon Institute, University of Cambridge, Cambridge, UK.
| |
Collapse
|
12
|
Liang Y, Kelemen A. Dynamic modeling and network approaches for omics time course data: overview of computational approaches and applications. Brief Bioinform 2019; 19:1051-1068. [PMID: 28430854 DOI: 10.1093/bib/bbx036] [Citation(s) in RCA: 20] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/21/2016] [Indexed: 12/23/2022] Open
Abstract
Inferring networks and dynamics of genes, proteins, cells and other biological entities from high-throughput biological omics data is a central and challenging issue in computational and systems biology. This is essential for understanding the complexity of human health, disease susceptibility and pathogenesis for Predictive, Preventive, Personalized and Participatory (P4) system and precision medicine. The delineation of the possible interactions of all genes/proteins in a genome/proteome is a task for which conventional experimental techniques are ill suited. Urgently needed are rapid and inexpensive computational and statistical methods that can identify interacting candidate disease genes or drug targets out of thousands that can be further investigated or validated by experimentations. Moreover, identifying biological dynamic systems, and simultaneously estimating the important kinetic structural and functional parameters, which may not be experimentally accessible could be important directions for drug-disease-gene network studies. In this article, we present an overview and comparison of recent developments of dynamic modeling and network approaches for time-course omics data, and their applications to various biological systems, health conditions and disease statuses. Moreover, various data reduction and analytical schemes ranging from mathematical to computational to statistical methods are compared including their merits, drawbacks and limitations. The most recent software, associated web resources and other potentials for the compared methods are also presented and discussed in detail.
Collapse
Affiliation(s)
- Yulan Liang
- Department of Family and Community Health, University of Maryland, Baltimore, MD, USA
| | - Arpad Kelemen
- Department of Family and Community Health, University of Maryland, Baltimore, MD, USA
| |
Collapse
|
13
|
Surkova S, Golubkova E, Mamon L, Samsonova M. Dynamic maternal gradients and morphogenetic networks in Drosophila early embryo. Biosystems 2018; 173:207-213. [DOI: 10.1016/j.biosystems.2018.10.009] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/14/2018] [Revised: 10/06/2018] [Accepted: 10/08/2018] [Indexed: 11/30/2022]
|
14
|
|
15
|
A damped oscillator imposes temporal order on posterior gap gene expression in Drosophila. PLoS Biol 2018; 16:e2003174. [PMID: 29451884 PMCID: PMC5832388 DOI: 10.1371/journal.pbio.2003174] [Citation(s) in RCA: 38] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/31/2017] [Revised: 03/01/2018] [Accepted: 01/31/2018] [Indexed: 12/21/2022] Open
Abstract
Insects determine their body segments in two different ways. Short-germband insects, such as the flour beetle Tribolium castaneum, use a molecular clock to establish segments sequentially. In contrast, long-germband insects, such as the vinegar fly Drosophila melanogaster, determine all segments simultaneously through a hierarchical cascade of gene regulation. Gap genes constitute the first layer of the Drosophila segmentation gene hierarchy, downstream of maternal gradients such as that of Caudal (Cad). We use data-driven mathematical modelling and phase space analysis to show that shifting gap domains in the posterior half of the Drosophila embryo are an emergent property of a robust damped oscillator mechanism, suggesting that the regulatory dynamics underlying long- and short-germband segmentation are much more similar than previously thought. In Tribolium, Cad has been proposed to modulate the frequency of the segmentation oscillator. Surprisingly, our simulations and experiments show that the shift rate of posterior gap domains is independent of maternal Cad levels in Drosophila. Our results suggest a novel evolutionary scenario for the short- to long-germband transition and help explain why this transition occurred convergently multiple times during the radiation of the holometabolan insects. Different insect species exhibit one of two distinct modes of determining their body segments (known as segmentation) during development: they either use a molecular oscillator to position segments sequentially, or they generate segments simultaneously through a hierarchical gene-regulatory cascade. The sequential mode is ancestral, while the simultaneous mode has been derived from it independently several times during evolution. In this paper, we present evidence suggesting that simultaneous segmentation also involves an oscillator in the posterior end of the embryo of the vinegar fly, Drosophila melanogaster. This surprising result indicates that both modes of segment determination are much more similar than previously thought. Such similarity provides an important step towards our understanding of the frequent evolutionary transitions observed between sequential and simultaneous segmentation.
Collapse
|
16
|
Erdogdu U, Polat F, Alhajj R. Employing decomposable partially observable Markov decision processes to control gene regulatory networks. Artif Intell Med 2017; 83:14-34. [DOI: 10.1016/j.artmed.2017.06.007] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/06/2017] [Revised: 04/28/2017] [Accepted: 06/14/2017] [Indexed: 11/30/2022]
|
17
|
|
18
|
Verd B, Crombach A, Jaeger J. Dynamic Maternal Gradients Control Timing and Shift-Rates for Drosophila Gap Gene Expression. PLoS Comput Biol 2017; 13:e1005285. [PMID: 28158178 PMCID: PMC5291410 DOI: 10.1371/journal.pcbi.1005285] [Citation(s) in RCA: 45] [Impact Index Per Article: 6.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/22/2016] [Accepted: 12/06/2016] [Indexed: 11/18/2022] Open
Abstract
Pattern formation during development is a highly dynamic process. In spite of this, few experimental and modelling approaches take into account the explicit time-dependence of the rules governing regulatory systems. We address this problem by studying dynamic morphogen interpretation by the gap gene network in Drosophila melanogaster. Gap genes are involved in segment determination during early embryogenesis. They are activated by maternal morphogen gradients encoded by bicoid (bcd) and caudal (cad). These gradients decay at the same time-scale as the establishment of the antero-posterior gap gene pattern. We use a reverse-engineering approach, based on data-driven regulatory models called gene circuits, to isolate and characterise the explicitly time-dependent effects of changing morphogen concentrations on gap gene regulation. To achieve this, we simulate the system in the presence and absence of dynamic gradient decay. Comparison between these simulations reveals that maternal morphogen decay controls the timing and limits the rate of gap gene expression. In the anterior of the embyro, it affects peak expression and leads to the establishment of smooth spatial boundaries between gap domains. In the posterior of the embryo, it causes a progressive slow-down in the rate of gap domain shifts, which is necessary to correctly position domain boundaries and to stabilise the spatial gap gene expression pattern. We use a newly developed method for the analysis of transient dynamics in non-autonomous (time-variable) systems to understand the regulatory causes of these effects. By providing a rigorous mechanistic explanation for the role of maternal gradient decay in gap gene regulation, our study demonstrates that such analyses are feasible and reveal important aspects of dynamic gene regulation which would have been missed by a traditional steady-state approach. More generally, it highlights the importance of transient dynamics for understanding complex regulatory processes in development. Animal development is a highly dynamic process. Biochemical or environmental signals can cause the rules that shape it to change over time. We know little about the effects of such changes. For the sake of simplicity, we usually leave them out of our models and experimental assays. Here, we do exactly the opposite. We characterise precisely those aspects of pattern formation caused by changing signalling inputs to a gene regulatory network, the gap gene system of Drosophila melanogaster. Gap genes are involved in determining the body segments of flies and other insects during early development. Gradients of maternal morphogens activate the expression of the gap genes. These gradients are highly dynamic themselves, as they decay while being read out. We show that this decay controls the peak concentration of gap gene products, produces smooth boundaries of gene expression, and slows down the observed positional shifts of gap domains in the posterior of the embryo, thereby stabilising the spatial pattern. Our analysis demonstrates that the dynamics of gene regulation not only affect the timing, but also the positioning of gene expression. This suggests that we must pay closer attention to transient dynamic aspects of development than is currently the case.
Collapse
Affiliation(s)
- Berta Verd
- EMBL/CRG Systems Biology Research Unit, Centre for Genomic Regulation (CRG), The Barcelona Institute of Science and Technology, Barcelona, Spain
- Universitat Pompeu Fabra (UPF), Barcelona, Spain
- KLI Klosterneuburg, Klosterneuburg, Austria
- * E-mail: (BV); (JJ)
| | - Anton Crombach
- EMBL/CRG Systems Biology Research Unit, Centre for Genomic Regulation (CRG), The Barcelona Institute of Science and Technology, Barcelona, Spain
- Universitat Pompeu Fabra (UPF), Barcelona, Spain
- Center for Interdisciplinary Research in Biology (CIRB), Collège de France, CNRS, INSERM, PSL Research University, Paris, France
| | - Johannes Jaeger
- EMBL/CRG Systems Biology Research Unit, Centre for Genomic Regulation (CRG), The Barcelona Institute of Science and Technology, Barcelona, Spain
- Universitat Pompeu Fabra (UPF), Barcelona, Spain
- KLI Klosterneuburg, Klosterneuburg, Austria
- Wissenschaftskolleg zu Berlin, Berlin, Germany
- * E-mail: (BV); (JJ)
| |
Collapse
|
19
|
Discovering novel phenotypes with automatically inferred dynamic models: a partial melanocyte conversion in Xenopus. Sci Rep 2017; 7:41339. [PMID: 28128301 PMCID: PMC5269672 DOI: 10.1038/srep41339] [Citation(s) in RCA: 23] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/28/2016] [Accepted: 12/16/2016] [Indexed: 12/20/2022] Open
Abstract
Progress in regenerative medicine requires reverse-engineering cellular control networks to infer perturbations with desired systems-level outcomes. Such dynamic models allow phenotypic predictions for novel perturbations to be rapidly assessed in silico. Here, we analyzed a Xenopus model of conversion of melanocytes to a metastatic-like phenotype only previously observed in an all-or-none manner. Prior in vivo genetic and pharmacological experiments showed that individual animals either fully convert or remain normal, at some characteristic frequency after a given perturbation. We developed a Machine Learning method which inferred a model explaining this complex, stochastic all-or-none dataset. We then used this model to ask how a new phenotype could be generated: animals in which only some of the melanocytes converted. Systematically performing in silico perturbations, the model predicted that a combination of altanserin (5HTR2 inhibitor), reserpine (VMAT inhibitor), and VP16-XlCreb1 (constitutively active CREB) would break the all-or-none concordance. Remarkably, applying the predicted combination of three reagents in vivo revealed precisely the expected novel outcome, resulting in partial conversion of melanocytes within individuals. This work demonstrates the capability of automated analysis of dynamic models of signaling networks to discover novel phenotypes and predictively identify specific manipulations that can reach them.
Collapse
|
20
|
Lobo D, Levin M. Computing a Worm: Reverse-Engineering Planarian Regeneration. EMERGENCE, COMPLEXITY AND COMPUTATION 2017. [DOI: 10.1007/978-3-319-33921-4_24] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
|
21
|
Rothschild JB, Tsimiklis P, Siggia ED, François P. Predicting Ancestral Segmentation Phenotypes from Drosophila to Anopheles Using In Silico Evolution. PLoS Genet 2016; 12:e1006052. [PMID: 27227405 PMCID: PMC4882032 DOI: 10.1371/journal.pgen.1006052] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/24/2016] [Accepted: 04/23/2016] [Indexed: 12/23/2022] Open
Abstract
Molecular evolution is an established technique for inferring gene homology but regulatory DNA turns over so rapidly that inference of ancestral networks is often impossible. In silico evolution is used to compute the most parsimonious path in regulatory space for anterior-posterior patterning linking two Dipterian species. The expression pattern of gap genes has evolved between Drosophila (fly) and Anopheles (mosquito), yet one of their targets, eve, has remained invariant. Our model predicts that stripe 5 in fly disappears and a new posterior stripe is created in mosquito, thus eve stripe modules 3+7 and 4+6 in fly are homologous to 3+6 and 4+5 in mosquito. We can place Clogmia on this evolutionary pathway and it shares the mosquito homologies. To account for the evolution of the other pair-rule genes in the posterior we have to assume that the ancestral Dipterian utilized a dynamic method to phase those genes in relation to eve. The last common ancestor of the fruit fly (Drosophila) and mosquito (Anopheles) lived more than 200 Million years ago. Can we use available data on insects alive today to infer what their ancestor looked like? In this manuscript, we focus on early embryonic development, when stripes of genetic expression appear and define the location of insect segments (“segmentation”). We use an evolutionary algorithm to reconstruct and predict dynamics of genes controlling stripes in the last common ancestor of fly and mosquito. We predict a new and different combinatorial logic of stripe formation in mosquito compared to fly, which is fully consistent with development of intermediate species such as moth-fly (Clogmia). Our simulations further suggest that the dynamics of gene expression in this last common ancestor were similar to other insects, such as wasps (Nasonia). Our method illustrates how computational methods inspired by machine learning and non-linear physics can be used to infer gene dynamics in species that disappeared millions of years ago.
Collapse
Affiliation(s)
- Jeremy B. Rothschild
- Physics Department, McGill University, Ernest Rutherford Physics Building, Montreal, Quebec, Canada
| | - Panagiotis Tsimiklis
- Physics Department, McGill University, Ernest Rutherford Physics Building, Montreal, Quebec, Canada
| | - Eric D. Siggia
- Center for Studies in Physics and Biology, The Rockefeller University, New York, New York, United States of America
| | - Paul François
- Physics Department, McGill University, Ernest Rutherford Physics Building, Montreal, Quebec, Canada
- * E-mail:
| |
Collapse
|
22
|
Lobo D, Morokuma J, Levin M. Computational discovery andin vivovalidation ofhnf4as a regulatory gene in planarian regeneration. Bioinformatics 2016; 32:2681-5. [DOI: 10.1093/bioinformatics/btw299] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/08/2016] [Accepted: 05/04/2016] [Indexed: 11/14/2022] Open
|
23
|
Buetti-Dinh A, Pivkin IV, Friedman R. S100A4 and its role in metastasis – computational integration of data on biological networks. MOLECULAR BIOSYSTEMS 2016; 11:2238-46. [PMID: 26118552 DOI: 10.1039/c5mb00110b] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
Abstract
Characterising signal transduction networks is fundamental to our understanding of biology. However, redundancy and different types of feedback mechanisms make it difficult to understand how variations of the network components contribute to a biological process. In silico modelling of signalling interactions therefore becomes increasingly useful for the development of successful therapeutic approaches. Unfortunately, quantitative information cannot be obtained for all of the proteins or complexes that comprise the network, which limits the usability of computational models. We developed a flexible computational framework for the analysis of biological signalling networks. We demonstrate our approach by studying the mechanism of metastasis promotion by the S100A4 protein, and suggest therapeutic strategies. The advantage of the proposed method is that only limited information (interaction type between species) is required to set up a steady-state network model. This permits a straightforward integration of experimental information where the lack of details are compensated by efficient sampling of the parameter space. We investigated regulatory properties of the S100A4 network and the role of different key components. The results show that S100A4 enhances the activity of matrix metalloproteinases (MMPs), causing higher cell dissociation. Moreover, it leads to an increased stability of the pathological state. Thus, avoiding metastasis in S100A4-expressing tumours requires multiple target inhibition. Moreover, the analysis could explain the previous failure of MMP inhibitors in clinical trials. Finally, our method is applicable to a wide range of biological questions that can be represented as directional networks.
Collapse
Affiliation(s)
- Antoine Buetti-Dinh
- Department of Chemistry and Biomedical Sciences, Linnæus University, Kalmar, Sweden.
| | | | | |
Collapse
|
24
|
Lou Z, Reinitz J. Parallel Simulated Annealing Using an Adaptive Resampling Interval. PARALLEL COMPUTING 2016; 53:23-31. [PMID: 26941469 PMCID: PMC4770898 DOI: 10.1016/j.parco.2016.02.001] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/05/2023]
Abstract
This paper presents a parallel simulated annealing algorithm that is able to achieve 90% parallel efficiency in iteration on up to 192 processors and up to 40% parallel efficiency in time when applied to a 5000-dimension Rastrigin function. Our algorithm breaks scalability barriers in the method of Chu et al. (1999) by abandoning adaptive cooling based on variance. The resulting gains in parallel efficiency are much larger than the loss of serial efficiency from lack of adaptive cooling. Our algorithm resamples the states across processors periodically. The resampling interval is tuned according to the success rate for each specific number of processors. We further present an adaptive method to determine the resampling interval based on the adoption rate. This adaptive method is able to achieve nearly identical parallel efficiency but higher success rates compared to the fixed interval one using the best interval found.
Collapse
Affiliation(s)
- Zhihao Lou
- Department of Computer Science, the University of Chicago,
Chicago, Illinois, USA
| | - John Reinitz
- Department of Statistics, Department of Ecology and
Evolution, Department of Molecular Genetics and Cell Biology, Institute of Genomics
and Systems Biology, the University of Chicago, Chicago, Illinois, USA
| |
Collapse
|
25
|
Durant F, Lobo D, Hammelman J, Levin M. Physiological controls of large-scale patterning in planarian regeneration: a molecular and computational perspective on growth and form. REGENERATION (OXFORD, ENGLAND) 2016; 3:78-102. [PMID: 27499881 PMCID: PMC4895326 DOI: 10.1002/reg2.54] [Citation(s) in RCA: 39] [Impact Index Per Article: 4.9] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 12/15/2015] [Revised: 02/18/2016] [Accepted: 02/22/2016] [Indexed: 12/12/2022]
Abstract
Planaria are complex metazoans that repair damage to their bodies and cease remodeling when a correct anatomy has been achieved. This model system offers a unique opportunity to understand how large-scale anatomical homeostasis emerges from the activities of individual cells. Much progress has been made on the molecular genetics of stem cell activity in planaria. However, recent data also indicate that the global pattern is regulated by physiological circuits composed of ionic and neurotransmitter signaling. Here, we overview the multi-scale problem of understanding pattern regulation in planaria, with specific focus on bioelectric signaling via ion channels and gap junctions (electrical synapses), and computational efforts to extract explanatory models from functional and molecular data on regeneration. We present a perspective that interprets results in this fascinating field using concepts from dynamical systems theory and computational neuroscience. Serving as a tractable nexus between genetic, physiological, and computational approaches to pattern regulation, planarian pattern homeostasis harbors many deep insights for regenerative medicine, evolutionary biology, and engineering.
Collapse
Affiliation(s)
- Fallon Durant
- Department of Biology, Allen Discovery Center at Tufts University, Tufts Center for Regenerative and Developmental BiologyTufts UniversityMA02155USA
| | - Daniel Lobo
- Department of Biological SciencesUniversity of MarylandBaltimore County, 1000 Hilltop CircleBaltimoreMD21250USA
| | - Jennifer Hammelman
- Department of Biology, Allen Discovery Center at Tufts University, Tufts Center for Regenerative and Developmental BiologyTufts UniversityMA02155USA
| | - Michael Levin
- Department of Biology, Allen Discovery Center at Tufts University, Tufts Center for Regenerative and Developmental BiologyTufts UniversityMA02155USA
| |
Collapse
|
26
|
Lobo D, Hammelman J, Levin M. MoCha: Molecular Characterization of Unknown Pathways. J Comput Biol 2016; 23:291-7. [PMID: 26950055 DOI: 10.1089/cmb.2015.0211] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022] Open
Abstract
Automated methods for the reverse-engineering of complex regulatory networks are paving the way for the inference of mechanistic comprehensive models directly from experimental data. These novel methods can infer not only the relations and parameters of the known molecules defined in their input datasets, but also unknown components and pathways identified as necessary by the automated algorithms. Identifying the molecular nature of these unknown components is a crucial step for making testable predictions and experimentally validating the models, yet no specific and efficient tools exist to aid in this process. To this end, we present here MoCha (Molecular Characterization), a tool optimized for the search of unknown proteins and their pathways from a given set of known interacting proteins. MoCha uses the comprehensive dataset of protein-protein interactions provided by the STRING database, which currently includes more than a billion interactions from over 2,000 organisms. MoCha is highly optimized, performing typical searches within seconds. We demonstrate the use of MoCha with the characterization of unknown components from reverse-engineered models from the literature. MoCha is useful for working on network models by hand or as a downstream step of a model inference engine workflow and represents a valuable and efficient tool for the characterization of unknown pathways using known data from thousands of organisms. MoCha and its source code are freely available online under the GPLv3 license.
Collapse
Affiliation(s)
- Daniel Lobo
- 1 Department of Biological Sciences, University of Maryland , Baltimore County, Baltimore, Maryland
| | - Jennifer Hammelman
- 2 Center for Regenerative and Developmental Biology, and Department of Biology, Tufts University , Medford, Massachusetts
| | - Michael Levin
- 2 Center for Regenerative and Developmental Biology, and Department of Biology, Tufts University , Medford, Massachusetts
| |
Collapse
|
27
|
Hoermann A, Cicin-Sain D, Jaeger J. A quantitative validated model reveals two phases of transcriptional regulation for the gap gene giant in Drosophila. Dev Biol 2016; 411:325-338. [DOI: 10.1016/j.ydbio.2016.01.005] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/24/2015] [Revised: 12/22/2015] [Accepted: 01/08/2016] [Indexed: 01/05/2023]
|
28
|
Crombach A, Wotton KR, Jiménez-Guri E, Jaeger J. Gap Gene Regulatory Dynamics Evolve along a Genotype Network. Mol Biol Evol 2016; 33:1293-307. [PMID: 26796549 PMCID: PMC4839219 DOI: 10.1093/molbev/msw013] [Citation(s) in RCA: 51] [Impact Index Per Article: 6.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022] Open
Abstract
Developmental gene networks implement the dynamic regulatory mechanisms that pattern and shape the organism. Over evolutionary time, the wiring of these networks changes, yet the patterning outcome is often preserved, a phenomenon known as “system drift.” System drift is illustrated by the gap gene network—involved in segmental patterning—in dipteran insects. In the classic model organism Drosophila melanogaster and the nonmodel scuttle fly Megaselia abdita, early activation and placement of gap gene expression domains show significant quantitative differences, yet the final patterning output of the system is essentially identical in both species. In this detailed modeling analysis of system drift, we use gene circuits which are fit to quantitative gap gene expression data in M. abdita and compare them with an equivalent set of models from D. melanogaster. The results of this comparative analysis show precisely how compensatory regulatory mechanisms achieve equivalent final patterns in both species. We discuss the larger implications of the work in terms of “genotype networks” and the ways in which the structure of regulatory networks can influence patterns of evolutionary change (evolvability).
Collapse
Affiliation(s)
- Anton Crombach
- EMBL/CRG Systems Biology Research Unit, Centre for Genomic Regulation (CRG), The Barcelona Institute of Science and Technology, Dr. Aiguader 88, Barcelona 08003, Spain Universitat Pompeu Fabra (UPF), Barcelona, Spain
| | - Karl R Wotton
- EMBL/CRG Systems Biology Research Unit, Centre for Genomic Regulation (CRG), The Barcelona Institute of Science and Technology, Dr. Aiguader 88, Barcelona 08003, Spain Universitat Pompeu Fabra (UPF), Barcelona, Spain
| | - Eva Jiménez-Guri
- EMBL/CRG Systems Biology Research Unit, Centre for Genomic Regulation (CRG), The Barcelona Institute of Science and Technology, Dr. Aiguader 88, Barcelona 08003, Spain Universitat Pompeu Fabra (UPF), Barcelona, Spain
| | - Johannes Jaeger
- EMBL/CRG Systems Biology Research Unit, Centre for Genomic Regulation (CRG), The Barcelona Institute of Science and Technology, Dr. Aiguader 88, Barcelona 08003, Spain Universitat Pompeu Fabra (UPF), Barcelona, Spain
| |
Collapse
|
29
|
Temporal and spatial dynamics of scaling-specific features of a gene regulatory network in Drosophila. Nat Commun 2015; 6:10031. [PMID: 26644070 PMCID: PMC4686680 DOI: 10.1038/ncomms10031] [Citation(s) in RCA: 38] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/13/2015] [Accepted: 10/28/2015] [Indexed: 01/19/2023] Open
Abstract
A widely appreciated aspect of developmental robustness is pattern formation in proportion to size. But how such scaling features emerge dynamically remains poorly understood. Here we generate a data set of the expression profiles of six gap genes in Drosophila melanogaster embryos that differ significantly in size. Expression patterns exhibit size-dependent dynamics both spatially and temporally. We uncover a dynamic emergence of under-scaling in the posterior, accompanied by reduced expression levels of gap genes near the middle of large embryos. Simulation results show that a size-dependent Bicoid gradient input can lead to reduced Krüppel expression that can have long-range and dynamic effects on gap gene expression in the posterior. Thus, for emergence of scaled patterns, the entire embryo may be viewed as a single unified dynamic system where maternally derived size-dependent information interpreted locally can be propagated in space and time as governed by the dynamics of a gene regulatory network. How pattern formation is regulated relative to the size of an organism is unclear. Here, Wu et al. take data from gap gene expression in flies of different sizes together with simulations, identifying how scaling emerges dynamically and that local patterning influences global gene regulatory networks.
Collapse
|
30
|
Lobikin M, Lobo D, Blackiston DJ, Martyniuk CJ, Tkachenko E, Levin M. Serotonergic regulation of melanocyte conversion: A bioelectrically regulated network for stochastic all-or-none hyperpigmentation. Sci Signal 2015; 8:ra99. [PMID: 26443706 DOI: 10.1126/scisignal.aac6609] [Citation(s) in RCA: 45] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
Experimentally induced depolarization of resting membrane potential in "instructor cells" in Xenopus laevis embryos causes hyperpigmentation in an all-or-none fashion in some tadpoles due to excess proliferation and migration of melanocytes. We showed that this stochastic process involved serotonin signaling, adenosine 3',5'-monophosphate (cAMP), and the transcription factors cAMP response element-binding protein (CREB), Sox10, and Slug. Transcriptional microarray analysis of embryos taken at stage 15 (early neurula) and stage 45 (free-swimming tadpole) revealed changes in the abundance of 45 and 517 transcripts, respectively, between control embryos and embryos exposed to the instructor cell-depolarizing agent ivermectin. Bioinformatic analysis revealed that the human homologs of some of the differentially regulated genes were associated with cancer, consistent with the induced arborization and invasive behavior of converted melanocytes. We identified a physiological circuit that uses serotonergic signaling between instructor cells, melanotrope cells of the pituitary, and melanocytes to control the proliferation, cell shape, and migration properties of the pigment cell pool. To understand the stochasticity and properties of this multiscale signaling system, we applied a computational machine-learning method that iteratively explored network models to reverse-engineer a stochastic dynamic model that recapitulated the frequency of the all-or-none hyperpigmentation phenotype produced in response to various pharmacological and molecular genetic manipulations. This computational approach may provide insight into stochastic cellular decision-making that occurs during normal development and pathological conditions, such as cancer.
Collapse
Affiliation(s)
- Maria Lobikin
- Biology Department and Center for Regenerative and Developmental Biology, Tufts University, Medford, MA 02155, USA
| | - Daniel Lobo
- Department of Biological Sciences, University of Maryland, Baltimore County, Baltimore, MD 21250, USA
| | - Douglas J Blackiston
- Biology Department and Center for Regenerative and Developmental Biology, Tufts University, Medford, MA 02155, USA
| | - Christopher J Martyniuk
- Center for Environmental and Human Toxicology and Department of Physiological Sciences, UF Genetics Institute, University of Florida, Gainesville, FL 32611, USA
| | - Elizabeth Tkachenko
- Biology Department and Center for Regenerative and Developmental Biology, Tufts University, Medford, MA 02155, USA
| | - Michael Levin
- Biology Department and Center for Regenerative and Developmental Biology, Tufts University, Medford, MA 02155, USA.
| |
Collapse
|
31
|
Fusco G. For a new dialogue between theoretical and empirical studies in evo-devo. Front Ecol Evol 2015. [DOI: 10.3389/fevo.2015.00097] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/13/2022] Open
|
32
|
Uzkudun M, Marcon L, Sharpe J. Data-driven modelling of a gene regulatory network for cell fate decisions in the growing limb bud. Mol Syst Biol 2015; 11:815. [PMID: 26174932 PMCID: PMC4547844 DOI: 10.15252/msb.20145882] [Citation(s) in RCA: 32] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/03/2022] Open
Abstract
Parameter optimization coupled with model selection is a convenient approach to infer gene regulatory networks from experimental gene expression data, but so far it has been limited to single cells or static tissues where growth is not significant. Here, we present a computational study in which we determine an optimal gene regulatory network from the spatiotemporal dynamics of gene expression patterns in a complex 2D growing tissue (non-isotropic and heterogeneous growth rates). We use this method to predict the regulatory mechanisms that underlie proximodistal (PD) patterning of the developing limb bud. First, we map the expression patterns of the PD markers Meis1, Hoxa11 and Hoxa13 into a dynamic description of the tissue movements that drive limb morphogenesis. Secondly, we use reverse-engineering to test how different gene regulatory networks can interpret the opposing gradients of fibroblast growth factors (FGF) and retinoic acid (RA) to pattern the PD markers. Finally, we validate and extend the best model against various previously published manipulative experiments, including exogenous application of RA, surgical removal of the FGF source and genetic ectopic expression of Meis1. Our approach identifies the most parsimonious gene regulatory network that can correctly pattern the PD markers downstream of FGF and RA. This network reveals a new model of PD regulation which we call the “crossover model”, because the proximal morphogen (RA) controls the distal boundary of Hoxa11, while conversely the distal morphogens (FGFs) control the proximal boundary.
Collapse
Affiliation(s)
- Manu Uzkudun
- EMBL-CRG Systems Biology Program Centre for Genomic Regulation (CRG) Universitat Pompeu Fabra (UPF), Barcelona Spain
| | - Luciano Marcon
- EMBL-CRG Systems Biology Program Centre for Genomic Regulation (CRG) Universitat Pompeu Fabra (UPF), Barcelona Spain
| | - James Sharpe
- EMBL-CRG Systems Biology Program Centre for Genomic Regulation (CRG) Universitat Pompeu Fabra (UPF), Barcelona Spain Institucio Catalana de Recerca i Estudis Avancats (ICREA), Barcelona, Spain
| |
Collapse
|
33
|
Inferring regulatory networks from experimental morphological phenotypes: a computational method reverse-engineers planarian regeneration. PLoS Comput Biol 2015; 11:e1004295. [PMID: 26042810 PMCID: PMC4456145 DOI: 10.1371/journal.pcbi.1004295] [Citation(s) in RCA: 57] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/24/2014] [Accepted: 04/21/2015] [Indexed: 01/18/2023] Open
Abstract
Transformative applications in biomedicine require the discovery of complex regulatory networks that explain the development and regeneration of anatomical structures, and reveal what external signals will trigger desired changes of large-scale pattern. Despite recent advances in bioinformatics, extracting mechanistic pathway models from experimental morphological data is a key open challenge that has resisted automation. The fundamental difficulty of manually predicting emergent behavior of even simple networks has limited the models invented by human scientists to pathway diagrams that show necessary subunit interactions but do not reveal the dynamics that are sufficient for complex, self-regulating pattern to emerge. To finally bridge the gap between high-resolution genetic data and the ability to understand and control patterning, it is critical to develop computational tools to efficiently extract regulatory pathways from the resultant experimental shape phenotypes. For example, planarian regeneration has been studied for over a century, but despite increasing insight into the pathways that control its stem cells, no constructive, mechanistic model has yet been found by human scientists that explains more than one or two key features of its remarkable ability to regenerate its correct anatomical pattern after drastic perturbations. We present a method to infer the molecular products, topology, and spatial and temporal non-linear dynamics of regulatory networks recapitulating in silico the rich dataset of morphological phenotypes resulting from genetic, surgical, and pharmacological experiments. We demonstrated our approach by inferring complete regulatory networks explaining the outcomes of the main functional regeneration experiments in the planarian literature; By analyzing all the datasets together, our system inferred the first systems-biology comprehensive dynamical model explaining patterning in planarian regeneration. This method provides an automated, highly generalizable framework for identifying the underlying control mechanisms responsible for the dynamic regulation of growth and form. Developmental and regenerative biology experiments are producing a huge number of morphological phenotypes from functional perturbation experiments. However, existing pathway models do not generally explain the dynamic regulation of anatomical shape due to the difficulty of inferring and testing non-linear regulatory networks responsible for appropriate form, shape, and pattern. We present a method that automates the discovery and testing of regulatory networks explaining morphological outcomes directly from the resultant phenotypes, producing network models as testable hypotheses explaining regeneration data. Our system integrates a formalization of the published results in planarian regeneration, an in silico simulator in which the patterning properties of regulatory networks can be quantitatively tested in a regeneration assay, and a machine learning module that evolves networks whose behavior in this assay optimally matches the database of planarian results. We applied our method to explain the key experiments in planarian regeneration, and discovered the first comprehensive model of anterior-posterior patterning in planaria under surgical, pharmacological, and genetic manipulations. Beyond the planarian data, our approach is readily generalizable to facilitate the discovery of testable regulatory networks in developmental biology and biomedicine, and represents the first developmental model discovered de novo from morphological outcomes by an automated system.
Collapse
|
34
|
Wotton KR, Jiménez-Guri E, Jaeger J. Maternal co-ordinate gene regulation and axis polarity in the scuttle fly Megaselia abdita. PLoS Genet 2015; 11:e1005042. [PMID: 25757102 PMCID: PMC4355411 DOI: 10.1371/journal.pgen.1005042] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/20/2014] [Accepted: 01/30/2015] [Indexed: 02/01/2023] Open
Abstract
Axis specification and segment determination in dipteran insects are an excellent model system for comparative analyses of gene network evolution. Antero-posterior polarity of the embryo is established through systems of maternal morphogen gradients. In Drosophila melanogaster, the anterior system acts through opposing gradients of Bicoid (Bcd) and Caudal (Cad), while the posterior system involves Nanos (Nos) and Hunchback (Hb) protein. These systems act redundantly. Both Bcd and Hb need to be eliminated to cause a complete loss of polarity resulting in mirror-duplicated abdomens, so-called bicaudal phenotypes. In contrast, knock-down of bcd alone is sufficient to induce double abdomens in non-drosophilid cyclorrhaphan dipterans such as the hoverfly Episyrphus balteatus or the scuttle fly Megaselia abdita. We investigate conserved and divergent aspects of axis specification in the cyclorrhaphan lineage through a detailed study of the establishment and regulatory effect of maternal gradients in M. abdita. Our results show that the function of the anterior maternal system is highly conserved in this species, despite the loss of maternal cad expression. In contrast, hb does not activate gap genes in this species. The absence of this activatory role provides a precise genetic explanation for the loss of polarity upon bcd knock-down in M. abdita, and suggests a general scenario in which the posterior maternal system is increasingly replaced by the anterior one during the evolution of the cyclorrhaphan dipteran lineage. The basic head-to-tail polarity of an animal is established very early in development. In dipteran insects (flies, midges, and mosquitoes), polarity is established with the help of so-called morphogen gradients. Morphogens are regulatory proteins that are distributed as a concentration gradient, often involving diffusion from a localised source. This graded distribution then leads to the concentration-dependent activation of different target genes along the embryo’s axis. We examine this process, which differs to a surprising extent between dipteran species, in the scuttle fly Megaselia abdita, and compare our results to the model organism Drosophila melanogaster. In this way, we not only gain insights into how the mechanisms that establish polarity function differently in different species, but also how the system has evolved since these two flies shared a common ancestor. Specifically, we pin down the main difference between Drosophila and Megaselia in the altered function of the maternal Hunchback morphogen gradient, which activates target genes in the former, but not the latter species, where it has been completely replaced by the Bicoid morphogen during evolution.
Collapse
Affiliation(s)
- Karl R. Wotton
- EMBL/CRG Research Unit in Systems Biology, Centre for Genomic Regulation (CRG), Barcelona, Spain
- Universitat Pompeu Fabra (UPF), Barcelona, Spain
- * E-mail: (KW); (JJ)
| | - Eva Jiménez-Guri
- EMBL/CRG Research Unit in Systems Biology, Centre for Genomic Regulation (CRG), Barcelona, Spain
- Universitat Pompeu Fabra (UPF), Barcelona, Spain
| | - Johannes Jaeger
- EMBL/CRG Research Unit in Systems Biology, Centre for Genomic Regulation (CRG), Barcelona, Spain
- Universitat Pompeu Fabra (UPF), Barcelona, Spain
- * E-mail: (KW); (JJ)
| |
Collapse
|
35
|
Wotton KR, Jiménez-Guri E, Crombach A, Cicin-Sain D, Jaeger J. High-resolution gene expression data from blastoderm embryos of the scuttle fly Megaselia abdita. Sci Data 2015; 2:150005. [PMID: 25977812 PMCID: PMC4423355 DOI: 10.1038/sdata.2015.5] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/19/2014] [Accepted: 02/03/2015] [Indexed: 12/21/2022] Open
Abstract
Gap genes are involved in segment determination during early development in dipteran insects (flies, midges, and mosquitoes). We carried out a systematic quantitative comparative analysis of the gap gene network across different dipteran species. Our work provides mechanistic insights into the evolution of this pattern-forming network. As a central component of our project, we created a high-resolution quantitative spatio-temporal data set of gap and maternal co-ordinate gene expression in the blastoderm embryo of the non-drosophilid scuttle fly, Megaselia abdita. Our data include expression patterns in both wild-type and RNAi-treated embryos. The data-covering 10 genes, 10 time points, and over 1,000 individual embryos-consist of original embryo images, quantified expression profiles, extracted positions of expression boundaries, and integrated expression patterns, plus metadata and intermediate processing steps. These data provide a valuable resource for researchers interested in the comparative study of gene regulatory networks and pattern formation, an essential step towards a more quantitative and mechanistic understanding of developmental evolution.
Collapse
Affiliation(s)
- Karl R Wotton
- EMBL/CRG Research Unit in Systems Biology, Centre for Genomic Regulation (CRG), 08003 Barcelona, Spain
- Universitat Pompeu Fabra (UPF), 08002 Barcelona, Spain
| | - Eva Jiménez-Guri
- EMBL/CRG Research Unit in Systems Biology, Centre for Genomic Regulation (CRG), 08003 Barcelona, Spain
- Universitat Pompeu Fabra (UPF), 08002 Barcelona, Spain
| | - Anton Crombach
- EMBL/CRG Research Unit in Systems Biology, Centre for Genomic Regulation (CRG), 08003 Barcelona, Spain
- Universitat Pompeu Fabra (UPF), 08002 Barcelona, Spain
| | - Damjan Cicin-Sain
- EMBL/CRG Research Unit in Systems Biology, Centre for Genomic Regulation (CRG), 08003 Barcelona, Spain
- Universitat Pompeu Fabra (UPF), 08002 Barcelona, Spain
| | - Johannes Jaeger
- EMBL/CRG Research Unit in Systems Biology, Centre for Genomic Regulation (CRG), 08003 Barcelona, Spain
- Universitat Pompeu Fabra (UPF), 08002 Barcelona, Spain
| |
Collapse
|
36
|
BioPreDyn-bench: a suite of benchmark problems for dynamic modelling in systems biology. BMC SYSTEMS BIOLOGY 2015; 9:8. [PMID: 25880925 PMCID: PMC4342829 DOI: 10.1186/s12918-015-0144-4] [Citation(s) in RCA: 52] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 07/21/2014] [Accepted: 01/15/2015] [Indexed: 11/21/2022]
Abstract
Background Dynamic modelling is one of the cornerstones of systems biology. Many research efforts are currently being invested in the development and exploitation of large-scale kinetic models. The associated problems of parameter estimation (model calibration) and optimal experimental design are particularly challenging. The community has already developed many methods and software packages which aim to facilitate these tasks. However, there is a lack of suitable benchmark problems which allow a fair and systematic evaluation and comparison of these contributions. Results Here we present BioPreDyn-bench, a set of challenging parameter estimation problems which aspire to serve as reference test cases in this area. This set comprises six problems including medium and large-scale kinetic models of the bacterium E. coli, baker’s yeast S. cerevisiae, the vinegar fly D. melanogaster, Chinese Hamster Ovary cells, and a generic signal transduction network. The level of description includes metabolism, transcription, signal transduction, and development. For each problem we provide (i) a basic description and formulation, (ii) implementations ready-to-run in several formats, (iii) computational results obtained with specific solvers, (iv) a basic analysis and interpretation. Conclusions This suite of benchmark problems can be readily used to evaluate and compare parameter estimation methods. Further, it can also be used to build test problems for sensitivity and identifiability analysis, model reduction and optimal experimental design methods. The suite, including codes and documentation, can be freely downloaded from the BioPreDyn-bench website, https://sites.google.com/site/biopredynbenchmarks/. Electronic supplementary material The online version of this article (doi:10.1186/s12918-015-0144-4) contains supplementary material, which is available to authorized users.
Collapse
|
37
|
Wotton KR, Jiménez-Guri E, Crombach A, Janssens H, Alcaine-Colet A, Lemke S, Schmidt-Ott U, Jaeger J. Quantitative system drift compensates for altered maternal inputs to the gap gene network of the scuttle fly Megaselia abdita. eLife 2015; 4:e04785. [PMID: 25560971 PMCID: PMC4337606 DOI: 10.7554/elife.04785] [Citation(s) in RCA: 40] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/16/2014] [Accepted: 01/02/2015] [Indexed: 12/20/2022] Open
Abstract
The segmentation gene network in insects can produce equivalent phenotypic outputs despite differences in upstream regulatory inputs between species. We investigate the mechanistic basis of this phenomenon through a systems-level analysis of the gap gene network in the scuttle fly Megaselia abdita (Phoridae). It combines quantification of gene expression at high spatio-temporal resolution with systematic knock-downs by RNA interference (RNAi). Initiation and dynamics of gap gene expression differ markedly between M. abdita and Drosophila melanogaster, while the output of the system converges to equivalent patterns at the end of the blastoderm stage. Although the qualitative structure of the gap gene network is conserved, there are differences in the strength of regulatory interactions between species. We term such network rewiring 'quantitative system drift'. It provides a mechanistic explanation for the developmental hourglass model in the dipteran lineage. Quantitative system drift is likely to be a widespread mechanism for developmental evolution.
Collapse
Affiliation(s)
- Karl R Wotton
- European Molecular Biology Laboratory, CRG Systems Biology Research Unit, Centre for Genomic Regulation, Barcelona, Spain
- Universitat Pompeu Fabra, Barcelona, Spain
| | - Eva Jiménez-Guri
- European Molecular Biology Laboratory, CRG Systems Biology Research Unit, Centre for Genomic Regulation, Barcelona, Spain
- Universitat Pompeu Fabra, Barcelona, Spain
| | - Anton Crombach
- European Molecular Biology Laboratory, CRG Systems Biology Research Unit, Centre for Genomic Regulation, Barcelona, Spain
- Universitat Pompeu Fabra, Barcelona, Spain
| | - Hilde Janssens
- European Molecular Biology Laboratory, CRG Systems Biology Research Unit, Centre for Genomic Regulation, Barcelona, Spain
- Universitat Pompeu Fabra, Barcelona, Spain
| | - Anna Alcaine-Colet
- European Molecular Biology Laboratory, CRG Systems Biology Research Unit, Centre for Genomic Regulation, Barcelona, Spain
- Universitat Pompeu Fabra, Barcelona, Spain
- Universitat de Barcelona, Barcelona, Spain
| | - Steffen Lemke
- Department of Organismal Biology and Anatomy, University of Chicago, Chicago, United States
| | - Urs Schmidt-Ott
- Department of Organismal Biology and Anatomy, University of Chicago, Chicago, United States
| | - Johannes Jaeger
- European Molecular Biology Laboratory, CRG Systems Biology Research Unit, Centre for Genomic Regulation, Barcelona, Spain
- Universitat Pompeu Fabra, Barcelona, Spain
| |
Collapse
|
38
|
Davila-Velderrain J, Martinez-Garcia JC, Alvarez-Buylla ER. Descriptive vs. mechanistic network models in plant development in the post-genomic era. Methods Mol Biol 2015; 1284:455-79. [PMID: 25757787 DOI: 10.1007/978-1-4939-2444-8_23] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/26/2023]
Abstract
Network modeling is now a widespread practice in systems biology, as well as in integrative genomics, and it constitutes a rich and diverse scientific research field. A conceptually clear understanding of the reasoning behind the main existing modeling approaches, and their associated technical terminologies, is required to avoid confusions and accelerate the transition towards an undeniable necessary more quantitative, multidisciplinary approach to biology. Herein, we focus on two main network-based modeling approaches that are commonly used depending on the information available and the intended goals: inference-based methods and system dynamics approaches. As far as data-based network inference methods are concerned, they enable the discovery of potential functional influences among molecular components. On the other hand, experimentally grounded network dynamical models have been shown to be perfectly suited for the mechanistic study of developmental processes. How do these two perspectives relate to each other? In this chapter, we describe and compare both approaches and then apply them to a given specific developmental module. Along with the step-by-step practical implementation of each approach, we also focus on discussing their respective goals, utility, assumptions, and associated limitations. We use the gene regulatory network (GRN) involved in Arabidopsis thaliana Root Stem Cell Niche patterning as our illustrative example. We show that descriptive models based on functional genomics data can provide important background information consistent with experimentally supported functional relationships integrated in mechanistic GRN models. The rationale of analysis and modeling can be applied to any other well-characterized functional developmental module in multicellular organisms, like plants and animals.
Collapse
Affiliation(s)
- J Davila-Velderrain
- Laboratorio de Genética Molecular, Desarrollo y Evolución de Plantas, Instituto de Ecología, Universidad Nacional Autónoma de México, Ciudad Universitaria, Av. Universidad 3000, México D.F., 04510, Mexico
| | | | | |
Collapse
|
39
|
Cicin-Sain D, Pulido AH, Crombach A, Wotton KR, Jiménez-Guri E, Taly JF, Roma G, Jaeger J. SuperFly: a comparative database for quantified spatio-temporal gene expression patterns in early dipteran embryos. Nucleic Acids Res 2014; 43:D751-5. [PMID: 25404137 PMCID: PMC4383950 DOI: 10.1093/nar/gku1142] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/25/2022] Open
Abstract
We present SuperFly (http://superfly.crg.eu), a relational database for quantified spatio-temporal expression data of segmentation genes during early development in different species of dipteran insects (flies, midges and mosquitoes). SuperFly has a special focus on emerging non-drosophilid model systems. The database currently includes data of high spatio-temporal resolution for three species: the vinegar fly Drosophila melanogaster, the scuttle fly Megaselia abdita and the moth midge Clogmia albipunctata. At this point, SuperFly covers up to 9 genes and 16 time points per species, with a total of 1823 individual embryos. It provides an intuitive web interface, enabling the user to query and access original embryo images, quantified expression profiles, extracted positions of expression boundaries and integrated datasets, plus metadata and intermediate processing steps. SuperFly is a valuable new resource for the quantitative comparative study of gene expression patterns across dipteran species. Moreover, it provides an interesting test set for systems biologists interested in fitting mathematical gene network models to data. Both of these aspects are essential ingredients for progress toward a more quantitative and mechanistic understanding of developmental evolution.
Collapse
Affiliation(s)
- Damjan Cicin-Sain
- EMBL/CRG Research Unit in Systems Biology, Centre for Genomic Regulation (CRG), 08003 Barcelona, Spain Bioinformatics Core Facility, Centre for Genomic Regulation (CRG), 08003 Barcelona, Spain
| | - Antonio Hermoso Pulido
- Universitat Pompeu Fabra (UPF), 08002 Barcelona, Spain Bioinformatics Core Facility, Centre for Genomic Regulation (CRG), 08003 Barcelona, Spain
| | - Anton Crombach
- EMBL/CRG Research Unit in Systems Biology, Centre for Genomic Regulation (CRG), 08003 Barcelona, Spain Bioinformatics Core Facility, Centre for Genomic Regulation (CRG), 08003 Barcelona, Spain
| | - Karl R Wotton
- EMBL/CRG Research Unit in Systems Biology, Centre for Genomic Regulation (CRG), 08003 Barcelona, Spain Bioinformatics Core Facility, Centre for Genomic Regulation (CRG), 08003 Barcelona, Spain
| | - Eva Jiménez-Guri
- EMBL/CRG Research Unit in Systems Biology, Centre for Genomic Regulation (CRG), 08003 Barcelona, Spain Bioinformatics Core Facility, Centre for Genomic Regulation (CRG), 08003 Barcelona, Spain
| | - Jean-François Taly
- Universitat Pompeu Fabra (UPF), 08002 Barcelona, Spain Bioinformatics Core Facility, Centre for Genomic Regulation (CRG), 08003 Barcelona, Spain
| | - Guglielmo Roma
- Universitat Pompeu Fabra (UPF), 08002 Barcelona, Spain Bioinformatics Core Facility, Centre for Genomic Regulation (CRG), 08003 Barcelona, Spain
| | - Johannes Jaeger
- EMBL/CRG Research Unit in Systems Biology, Centre for Genomic Regulation (CRG), 08003 Barcelona, Spain Bioinformatics Core Facility, Centre for Genomic Regulation (CRG), 08003 Barcelona, Spain
| |
Collapse
|
40
|
Wotton KR. Heterochronic shifts in germband movements contribute to the rapid embryonic development of the coffin fly Megaselia scalaris. ARTHROPOD STRUCTURE & DEVELOPMENT 2014; 43:589-594. [PMID: 25193503 DOI: 10.1016/j.asd.2014.08.001] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/16/2014] [Revised: 07/31/2014] [Accepted: 08/18/2014] [Indexed: 06/03/2023]
Abstract
The coffin fly, Megaselia scalaris, is a species of medical and forensic importance and is increasingly being used for the study of genetics. Postmortem interval can be estimated based on the life stage of M. scalaris recovered from corpses, therefore many studies have addressed the duration of each life stage. These studies demonstrate that embryogenesis completes significantly faster in M. scalaris than in the congener Megaselia abdita and faster even than the 24 h needed for Drosophila melanogaster embryogenesis. However, until now it has been unclear if this increased speed is achieved by reducing developmental time across all embryonic stages or by the acceleration of individual stages and processes. Here I use time-lapse imaging to create a staging scheme for M. scalaris embryogenesis. Comparison of stages between D. melanogaster and both Megaselia species reveals that heterochronic shifts, simultaneous morphogenetic movements and compression of individual stages all contribute to the rapid development of M. scalaris.
Collapse
Affiliation(s)
- Karl R Wotton
- EMBL/CRG Systems Biology Research Unit, Centre for Genomic Regulation (CRG), Dr. Aiguader 88, 08003 Barcelona, Spain; Universitat Pompeu Fabra (UPF), Barcelona, Spain.
| |
Collapse
|
41
|
He FQ, Wang W, Zheng P, Sudhakar P, Sun J, Zeng AP. Essential O2-responsive genes of Pseudomonas aeruginosa and their network revealed by integrating dynamic data from inverted conditions. Integr Biol (Camb) 2014; 6:215-23. [PMID: 24413814 DOI: 10.1039/c3ib40180d] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Abstract
Identification of the gene network through which Pseudomonas aeruginosa PAO1 (PA) adapts to altered oxygen-availability environments is essential for a better understanding of stress responses and pathogenicity of PA. We performed high-time-resolution (HTR) transcriptome analyses of PA in a continuous cultivation system during the transition from high oxygen tension to low oxygen tension (HLOT) and the reversed transition from low to high oxygen tension (LHOT). From those genes responsive to both transient conditions, we identified 85 essential oxygen-availability responsive genes (EORGs), including the expected ones (arcDABC) encoding enzymes for arginine fermentation. We then constructed the regulatory network for the EORGs of PA by integrating information from binding motif searching, literature and HTR data. Notably, our results show that only the sub-networks controlled by the well-known oxygen-responsive transcription factors show a very high consistency between the inferred network and literature knowledge, e.g. 87.5% and 83.3% of the obtained sub-network controlled by the anaerobic regulator (ANR) and a quorum sensing regulator RhIR, respectively. These results not only reveal stringent EORGs of PA and their transcription regulatory network, but also highlight that achieving a high accuracy of the inferred regulatory network might be feasible only for the apparently affected regulators under the given conditions but not for all the expressed regulators on a genome scale.
Collapse
Affiliation(s)
- Feng Q He
- Helmholtz Centre for Infection Research, D-38124, Braunschweig, Germany
| | | | | | | | | | | |
Collapse
|
42
|
Botman D, Röttinger E, Martindale MQ, de Jong J, Kaandorp JA. A computational approach towards a gene regulatory network for the developing Nematostella vectensis gut. PLoS One 2014; 9:e103341. [PMID: 25076223 PMCID: PMC4116165 DOI: 10.1371/journal.pone.0103341] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/24/2014] [Accepted: 06/26/2014] [Indexed: 11/19/2022] Open
Abstract
Background The starlet sea anemone Nematostella vectensis is a diploblastic cnidarian that expresses a set of conserved genes for gut formation during its early development. During the last decade, the spatial distribution of many of these genes has been visualized with RNA hybridization or protein immunolocalization techniques. However, due to N. vectensis' curved and changing morphology, quantification of these spatial data is problematic. A method is developed for two-dimensional gene expression quantification, which enables a numerical analysis and dynamic modeling of these spatial patterns. Methods/Result In this work, first standardized gene expression profiles are generated from publicly available N. vectensis embryo images that display mRNA and/or protein distributions. Then, genes expressed during gut formation are clustered based on their expression profiles, and further grouped based on temporal appearance of their gene products in embryonic development. Representative expression profiles are manually selected from these clusters, and used as input for a simulation-based optimization scheme. This scheme iteratively fits simulated profiles to the selected profiles, leading to an optimized estimation of the model parameters. Finally, a preliminary gene regulatory network is derived from the optimized model parameters. Outlook While the focus of this study is N. vectensis, the approach outlined here is suitable for inferring gene regulatory networks in the embryonic development of any animal, thus allowing to comparatively study gene regulation of gut formation in silico across various species.
Collapse
Affiliation(s)
- Daniel Botman
- Computational Science, University of Amsterdam, Amsterdam, The Netherlands
| | - Eric Röttinger
- Université Nice Sophia Antipolis, Institute for Research on Cancer and Aging, Nice (IRCAN), UMR 7284, Nice, France
- Centre National de la Recherche Scientifique (CNRS), Institute for Research on Cancer and Aging, Nice (IRCAN), UMR 7284, Nice, France
- Institut National de la Santé et de la Recherche Médicale (INSERM), Institute for Research on Cancer and Aging, Nice (IRCAN), U1081, Nice, France
| | - Mark Q. Martindale
- Whitney Lab for Marine Bioscience, University of Florida, St. Augustine, Florida, United States of America
| | - Johann de Jong
- Computational Cancer Biology Group, Division of Molecular Carcinogenesis, The Netherlands Cancer Institute, Amsterdam, the Netherlands
| | - Jaap A. Kaandorp
- Computational Science, University of Amsterdam, Amsterdam, The Netherlands
- * E-mail:
| |
Collapse
|
43
|
Marchand G, Huynh-Thu VA, Kane NC, Arribat S, Varès D, Rengel D, Balzergue S, Rieseberg LH, Vincourt P, Geurts P, Vignes M, Langlade NB. Bridging physiological and evolutionary time-scales in a gene regulatory network. THE NEW PHYTOLOGIST 2014; 203:685-696. [PMID: 24786523 DOI: 10.1111/nph.12818] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/15/2013] [Accepted: 03/17/2014] [Indexed: 06/03/2023]
Abstract
Gene regulatory networks (GRNs) govern phenotypic adaptations and reflect the trade-offs between physiological responses and evolutionary adaptation that act at different time-scales. To identify patterns of molecular function and genetic diversity in GRNs, we studied the drought response of the common sunflower, Helianthus annuus, and how the underlying GRN is related to its evolution. We examined the responses of 32,423 expressed sequences to drought and to abscisic acid (ABA) and selected 145 co-expressed transcripts. We characterized their regulatory relationships in nine kinetic studies based on different hormones. From this, we inferred a GRN by meta-analyses of a Gaussian graphical model and a random forest algorithm and studied the genetic differentiation among populations (FST ) at nodes. We identified two main hubs in the network that transport nitrate in guard cells. This suggests that nitrate transport is a critical aspect of the sunflower physiological response to drought. We observed that differentiation of the network genes in elite sunflower cultivars is correlated with their position and connectivity. This systems biology approach combined molecular data at different time-scales and identified important physiological processes. At the evolutionary level, we propose that network topology could influence responses to human selection and possibly adaptation to dry environments.
Collapse
Affiliation(s)
- Gwenaëlle Marchand
- INRA, Laboratoire des Interactions Plantes-Microorganismes (LIPM), UMR441, F-31326, Castanet-Tolosan, France
- CNRS, Laboratoire des Interactions Plantes-Microorganismes (LIPM), UMR2594, F-31326, Castanet-Tolosan, France
| | - Vân Anh Huynh-Thu
- Department of Electrical Engineering and Computer Science and GIGA-R, Systems and Modeling, University of Liège, Liège, Belgium
| | - Nolan C Kane
- Department of Ecology and Evolutionary Biology, University of Colorado at Boulder, Boulder, CO, 80309, USA
| | - Sandrine Arribat
- INRA, Unité de Recherche en Génomique Végétale (URGV), UMR1165 - Université d'Evry Val d'Essonne - ERL CNRS 8196, CP 5708, F-91057, Evry Cedex, France
| | - Didier Varès
- INRA, Laboratoire des Interactions Plantes-Microorganismes (LIPM), UMR441, F-31326, Castanet-Tolosan, France
- CNRS, Laboratoire des Interactions Plantes-Microorganismes (LIPM), UMR2594, F-31326, Castanet-Tolosan, France
| | - David Rengel
- INRA, Laboratoire des Interactions Plantes-Microorganismes (LIPM), UMR441, F-31326, Castanet-Tolosan, France
- CNRS, Laboratoire des Interactions Plantes-Microorganismes (LIPM), UMR2594, F-31326, Castanet-Tolosan, France
| | - Sandrine Balzergue
- INRA, Unité de Recherche en Génomique Végétale (URGV), UMR1165 - Université d'Evry Val d'Essonne - ERL CNRS 8196, CP 5708, F-91057, Evry Cedex, France
| | - Loren H Rieseberg
- Department of Botany, University of British Columbia, Vancouver, BC, V6T 1Z4, Canada
- Department of Biology, Indiana University, Bloomington, IN, 47405, USA
| | - Patrick Vincourt
- INRA, Laboratoire des Interactions Plantes-Microorganismes (LIPM), UMR441, F-31326, Castanet-Tolosan, France
- CNRS, Laboratoire des Interactions Plantes-Microorganismes (LIPM), UMR2594, F-31326, Castanet-Tolosan, France
| | - Pierre Geurts
- Department of Electrical Engineering and Computer Science and GIGA-R, Systems and Modeling, University of Liège, Liège, Belgium
| | - Matthieu Vignes
- INRA, Mathématiques et Informatique Appliquées (MIA), UPR875, F-31326, Castanet-Tolosan, France
| | - Nicolas B Langlade
- INRA, Laboratoire des Interactions Plantes-Microorganismes (LIPM), UMR441, F-31326, Castanet-Tolosan, France
- CNRS, Laboratoire des Interactions Plantes-Microorganismes (LIPM), UMR2594, F-31326, Castanet-Tolosan, France
| |
Collapse
|
44
|
Castro-González C, Luengo-Oroz MA, Duloquin L, Savy T, Rizzi B, Desnoulez S, Doursat R, Kergosien YL, Ledesma-Carbayo MJ, Bourgine P, Peyriéras N, Santos A. A digital framework to build, visualize and analyze a gene expression atlas with cellular resolution in zebrafish early embryogenesis. PLoS Comput Biol 2014; 10:e1003670. [PMID: 24945246 PMCID: PMC4063669 DOI: 10.1371/journal.pcbi.1003670] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/30/2013] [Accepted: 04/28/2014] [Indexed: 01/30/2023] Open
Abstract
A gene expression atlas is an essential resource to quantify and understand the multiscale processes of embryogenesis in time and space. The automated reconstruction of a prototypic 4D atlas for vertebrate early embryos, using multicolor fluorescence in situ hybridization with nuclear counterstain, requires dedicated computational strategies. To this goal, we designed an original methodological framework implemented in a software tool called Match-IT. With only minimal human supervision, our system is able to gather gene expression patterns observed in different analyzed embryos with phenotypic variability and map them onto a series of common 3D templates over time, creating a 4D atlas. This framework was used to construct an atlas composed of 6 gene expression templates from a cohort of zebrafish early embryos spanning 6 developmental stages from 4 to 6.3 hpf (hours post fertilization). They included 53 specimens, 181,415 detected cell nuclei and the segmentation of 98 gene expression patterns observed in 3D for 9 different genes. In addition, an interactive visualization software, Atlas-IT, was developed to inspect, supervise and analyze the atlas. Match-IT and Atlas-IT, including user manuals, representative datasets and video tutorials, are publicly and freely available online. We also propose computational methods and tools for the quantitative assessment of the gene expression templates at the cellular scale, with the identification, visualization and analysis of coexpression patterns, synexpression groups and their dynamics through developmental stages. We propose a workflow to map the expression domains of multiple genes onto a series of 3D templates, or “atlas”, during early embryogenesis. It was applied to the zebrafish at different stages between 4 and 6.3 hpf, generating 6 templates. Our system overcomes the lack of significant morphological landmarks in early development by relying on the expression of a reference gene (goosecoid, gsc) and nuclear staining to guide the registration of the analyzed genes. The proposed method also successfully maps gene domains from partially imaged embryos, thus allowing greater microscope magnification and cellular resolution. By using the workflow to construct a spatiotemporal database of zebrafish, we opened the way to a systematic analysis of vertebrate embryogenesis. The atlas database, together with the mapping software (Match-IT), a custom-made visualization platform (Atlas-IT), and step-by-step user guides are available from the Supplementary Material. We expect that this will encourage other laboratories to generate, map, visualize and analyze new gene expression datasets.
Collapse
Affiliation(s)
- Carlos Castro-González
- Biomedical Image Technologies, ETSIT, Universidad Politécnica de Madrid, CEIMoncloa, Madrid, Spain
- Research Center in Bioengineering, Biomaterials and Nanomedicine (CIBER-BBN), Madrid, Spain
- Madrid-MIT M+Visión Consortium, Massachusetts Institute of Technology, Cambridge, Massachusetts, United States of America
| | - Miguel A. Luengo-Oroz
- Biomedical Image Technologies, ETSIT, Universidad Politécnica de Madrid, CEIMoncloa, Madrid, Spain
- Research Center in Bioengineering, Biomaterials and Nanomedicine (CIBER-BBN), Madrid, Spain
| | - Louise Duloquin
- MDAM UPR3294, Institut de Neurobiologie Alfred Fessard, CNRS, Gif-sur-Yvette, France
- Institut des Systèmes Complexes, Paris, France
- BioEmergences-IBiSA, Institut de Neurobiologie Alfred Fessard, CNRS, Gif-sur-Yvette, France
| | - Thierry Savy
- MDAM UPR3294, Institut de Neurobiologie Alfred Fessard, CNRS, Gif-sur-Yvette, France
- Institut des Systèmes Complexes, Paris, France
- BioEmergences-IBiSA, Institut de Neurobiologie Alfred Fessard, CNRS, Gif-sur-Yvette, France
| | - Barbara Rizzi
- MDAM UPR3294, Institut de Neurobiologie Alfred Fessard, CNRS, Gif-sur-Yvette, France
- Institut des Systèmes Complexes, Paris, France
- BioEmergences-IBiSA, Institut de Neurobiologie Alfred Fessard, CNRS, Gif-sur-Yvette, France
| | - Sophie Desnoulez
- MDAM UPR3294, Institut de Neurobiologie Alfred Fessard, CNRS, Gif-sur-Yvette, France
- BioEmergences-IBiSA, Institut de Neurobiologie Alfred Fessard, CNRS, Gif-sur-Yvette, France
| | - René Doursat
- Institut des Systèmes Complexes, Paris, France
- BioEmergences-IBiSA, Institut de Neurobiologie Alfred Fessard, CNRS, Gif-sur-Yvette, France
- School of Biomedical Engineering, Drexel University, Philadelphia, Pennsylvania, United States of America
| | - Yannick L. Kergosien
- MDAM UPR3294, Institut de Neurobiologie Alfred Fessard, CNRS, Gif-sur-Yvette, France
- BioEmergences-IBiSA, Institut de Neurobiologie Alfred Fessard, CNRS, Gif-sur-Yvette, France
- LIMICS-INSERM UMR 1142, UFR SMBH, Université Paris 13, Bobigny, France
| | - María J. Ledesma-Carbayo
- Biomedical Image Technologies, ETSIT, Universidad Politécnica de Madrid, CEIMoncloa, Madrid, Spain
- Research Center in Bioengineering, Biomaterials and Nanomedicine (CIBER-BBN), Madrid, Spain
| | - Paul Bourgine
- MDAM UPR3294, Institut de Neurobiologie Alfred Fessard, CNRS, Gif-sur-Yvette, France
- Institut des Systèmes Complexes, Paris, France
- BioEmergences-IBiSA, Institut de Neurobiologie Alfred Fessard, CNRS, Gif-sur-Yvette, France
| | - Nadine Peyriéras
- MDAM UPR3294, Institut de Neurobiologie Alfred Fessard, CNRS, Gif-sur-Yvette, France
- Institut des Systèmes Complexes, Paris, France
- BioEmergences-IBiSA, Institut de Neurobiologie Alfred Fessard, CNRS, Gif-sur-Yvette, France
- * E-mail: (NP); (AS)
| | - Andrés Santos
- Biomedical Image Technologies, ETSIT, Universidad Politécnica de Madrid, CEIMoncloa, Madrid, Spain
- Research Center in Bioengineering, Biomaterials and Nanomedicine (CIBER-BBN), Madrid, Spain
- * E-mail: (NP); (AS)
| |
Collapse
|
45
|
Crombach A, García-Solache MA, Jaeger J. Evolution of early development in dipterans: reverse-engineering the gap gene network in the moth midge Clogmia albipunctata (Psychodidae). Biosystems 2014; 123:74-85. [PMID: 24911671 DOI: 10.1016/j.biosystems.2014.06.003] [Citation(s) in RCA: 34] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/24/2014] [Revised: 06/04/2014] [Accepted: 06/04/2014] [Indexed: 11/18/2022]
Abstract
Understanding the developmental and evolutionary dynamics of regulatory networks is essential if we are to explain the non-random distribution of phenotypes among the diversity of organismic forms. Here, we present a comparative analysis of one of the best understood developmental gene regulatory networks today: the gap gene network involved in early patterning of insect embryos. We use gene circuit models, which are fitted to quantitative spatio-temporal gene expression data for the four trunk gap genes hunchback (hb), Krüppel (Kr), giant (gt), and knirps (kni)/knirps-like (knl) in the moth midge Clogmia albipunctata, and compare them to equivalent reverse-engineered circuits from our reference species, the vinegar fly Drosophila melanogaster. In contrast to the single network structure we find for D. melanogaster, our models predict four alternative networks for C. albipunctata. These networks share a core structure, which includes the central regulatory feedback between hb and knl. Other interactions are only partially determined, as they differ between our four network structures. Nevertheless, our models make testable predictions and enable us to gain specific insights into gap gene regulation in C. albipunctata. They suggest a less central role for Kr in C. albipunctata than in D. melanogaster, and show that the mechanisms causing an anterior shift of gap domains over time are largely conserved between the two species, although shift dynamics differ. The set of C. albipunctata gene circuit models presented here will be used as the starting point for data-constrained in silico evolutionary simulations to study patterning transitions in the early development of dipteran species.
Collapse
Affiliation(s)
- Anton Crombach
- EMBL/CRG Research Unit in Systems Biology, Centre for Genomic Regulation (CRG), Dr. Aiguader 88, 08003 Barcelona, Spain; Universitat Pompeu Fabra (UPF), Barcelona, Spain.
| | - Mónica A García-Solache
- Laboratory for Development and Evolution, University Museum of Zoology and Department of Zoology, Downing Street, Cambridge CB2 3EJ, UK
| | - Johannes Jaeger
- EMBL/CRG Research Unit in Systems Biology, Centre for Genomic Regulation (CRG), Dr. Aiguader 88, 08003 Barcelona, Spain; Universitat Pompeu Fabra (UPF), Barcelona, Spain
| |
Collapse
|
46
|
Jug F, Pietzsch T, Preibisch S, Tomancak P. Bioimage Informatics in the context of Drosophila research. Methods 2014; 68:60-73. [PMID: 24732429 DOI: 10.1016/j.ymeth.2014.04.004] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/25/2014] [Revised: 04/02/2014] [Accepted: 04/04/2014] [Indexed: 01/05/2023] Open
Abstract
Modern biological research relies heavily on microscopic imaging. The advanced genetic toolkit of Drosophila makes it possible to label molecular and cellular components with unprecedented level of specificity necessitating the application of the most sophisticated imaging technologies. Imaging in Drosophila spans all scales from single molecules to the entire populations of adult organisms, from electron microscopy to live imaging of developmental processes. As the imaging approaches become more complex and ambitious, there is an increasing need for quantitative, computer-mediated image processing and analysis to make sense of the imagery. Bioimage Informatics is an emerging research field that covers all aspects of biological image analysis from data handling, through processing, to quantitative measurements, analysis and data presentation. Some of the most advanced, large scale projects, combining cutting edge imaging with complex bioimage informatics pipelines, are realized in the Drosophila research community. In this review, we discuss the current research in biological image analysis specifically relevant to the type of systems level image datasets that are uniquely available for the Drosophila model system. We focus on how state-of-the-art computer vision algorithms are impacting the ability of Drosophila researchers to analyze biological systems in space and time. We pay particular attention to how these algorithmic advances from computer science are made usable to practicing biologists through open source platforms and how biologists can themselves participate in their further development.
Collapse
Affiliation(s)
- Florian Jug
- Max Planck Institute of Molecular Cell Biology and Genetics, 01307 Dresden, Germany
| | - Tobias Pietzsch
- Max Planck Institute of Molecular Cell Biology and Genetics, 01307 Dresden, Germany
| | - Stephan Preibisch
- Janelia Farm Research Campus, Howard Hughes Medical Institute, Ashburn, VA 20147, USA; Department of Anatomy and Structural Biology, Gruss Lipper Biophotonics Center, Albert Einstein College of Medicine, Bronx, NY 10461, USA
| | - Pavel Tomancak
- Max Planck Institute of Molecular Cell Biology and Genetics, 01307 Dresden, Germany.
| |
Collapse
|
47
|
Mbodj A, Junion G, Brun C, Furlong EEM, Thieffry D. Logical modelling of Drosophila signalling pathways. MOLECULAR BIOSYSTEMS 2014; 9:2248-58. [PMID: 23868318 DOI: 10.1039/c3mb70187e] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
Abstract
A limited number of signalling pathways are involved in the specification of cell fate during the development of all animals. Several of these pathways were originally identified in Drosophila. To clarify their roles, and possible cross-talk, we have built a logical model for the nine key signalling pathways recurrently used in metazoan development. In each case, we considered the associated ligands, receptors, signal transducers, modulators, and transcription factors reported in the literature. Implemented using the logical modelling software GINsim, the resulting models qualitatively recapitulate the main characteristics of each pathway, in wild type as well as in various mutant situations (e.g. loss-of-function or gain-of-function). These models constitute pluggable modules that can be used to assemble comprehensive models of complex developmental processes. Moreover, these models of Drosophila pathways could serve as scaffolds for more complicated models of orthologous mammalian pathways. Comprehensive model annotations and GINsim files are provided for each of the nine considered pathways.
Collapse
Affiliation(s)
- Abibatou Mbodj
- Technological Advances for Genomics and Clinics (TAGC), INSERM UMR_S 1090, Aix-Marseille Université, Marseille, France.
| | | | | | | | | |
Collapse
|
48
|
Villaverde AF, Banga JR. Reverse engineering and identification in systems biology: strategies, perspectives and challenges. J R Soc Interface 2014; 11:20130505. [PMID: 24307566 PMCID: PMC3869153 DOI: 10.1098/rsif.2013.0505] [Citation(s) in RCA: 163] [Impact Index Per Article: 16.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/07/2013] [Accepted: 11/12/2013] [Indexed: 12/17/2022] Open
Abstract
The interplay of mathematical modelling with experiments is one of the central elements in systems biology. The aim of reverse engineering is to infer, analyse and understand, through this interplay, the functional and regulatory mechanisms of biological systems. Reverse engineering is not exclusive of systems biology and has been studied in different areas, such as inverse problem theory, machine learning, nonlinear physics, (bio)chemical kinetics, control theory and optimization, among others. However, it seems that many of these areas have been relatively closed to outsiders. In this contribution, we aim to compare and highlight the different perspectives and contributions from these fields, with emphasis on two key questions: (i) why are reverse engineering problems so hard to solve, and (ii) what methods are available for the particular problems arising from systems biology?
Collapse
Affiliation(s)
| | - Julio R. Banga
- BioProcess Engineering Group, IIM-CSIC, Spanish National Research Council, Vigo 36208, Spain
| |
Collapse
|
49
|
Zheng Z, Christley S, Chiu WT, Blitz IL, Xie X, Cho KWY, Nie Q. Inference of the Xenopus tropicalis embryonic regulatory network and spatial gene expression patterns. BMC SYSTEMS BIOLOGY 2014; 8:3. [PMID: 24397936 PMCID: PMC3896677 DOI: 10.1186/1752-0509-8-3] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 04/16/2013] [Accepted: 12/19/2013] [Indexed: 11/10/2022]
Abstract
BACKGROUND During embryogenesis, signaling molecules produced by one cell population direct gene regulatory changes in neighboring cells and influence their developmental fates and spatial organization. One of the earliest events in the development of the vertebrate embryo is the establishment of three germ layers, consisting of the ectoderm, mesoderm and endoderm. Attempts to measure gene expression in vivo in different germ layers and cell types are typically complicated by the heterogeneity of cell types within biological samples (i.e., embryos), as the responses of individual cell types are intermingled into an aggregate observation of heterogeneous cell types. Here, we propose a novel method to elucidate gene regulatory circuits from these aggregate measurements in embryos of the frog Xenopus tropicalis using gene network inference algorithms and then test the ability of the inferred networks to predict spatial gene expression patterns. RESULTS We use two inference models with different underlying assumptions that incorporate existing network information, an ODE model for steady-state data and a Markov model for time series data, and contrast the performance of the two models. We apply our method to both control and knockdown embryos at multiple time points to reconstruct the core mesoderm and endoderm regulatory circuits. Those inferred networks are then used in combination with known dorsal-ventral spatial expression patterns of a subset of genes to predict spatial expression patterns for other genes. Both models are able to predict spatial expression patterns for some of the core mesoderm and endoderm genes, but interestingly of different gene subsets, suggesting that neither model is sufficient to recapitulate all of the spatial patterns, yet they are complementary for the patterns that they do capture. CONCLUSION The presented methodology of gene network inference combined with spatial pattern prediction provides an additional layer of validation to elucidate the regulatory circuits controlling the spatial-temporal dynamics in embryonic development.
Collapse
Affiliation(s)
| | | | | | | | | | | | - Qing Nie
- Department of Mathematics, University of California, Irvine, CA 92697, USA.
| |
Collapse
|
50
|
Janssens H, Siggens K, Cicin-Sain D, Jiménez-Guri E, Musy M, Akam M, Jaeger J. A quantitative atlas of Even-skipped and Hunchback expression in Clogmia albipunctata (Diptera: Psychodidae) blastoderm embryos. EvoDevo 2014; 5:1. [PMID: 24393251 PMCID: PMC3897886 DOI: 10.1186/2041-9139-5-1] [Citation(s) in RCA: 36] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/17/2013] [Accepted: 11/22/2013] [Indexed: 11/13/2022] Open
Abstract
Background Comparative studies of developmental processes are one of the main approaches to evolutionary developmental biology (evo-devo). Over recent years, there has been a shift of focus from the comparative study of particular regulatory genes to the level of whole gene networks. Reverse-engineering methods can be used to computationally reconstitute and analyze the function and dynamics of such networks. These methods require quantitative spatio-temporal expression data for model fitting. Obtaining such data in non-model organisms remains a major technical challenge, impeding the wider application of data-driven mathematical modeling to evo-devo. Results We have raised antibodies against four segmentation gene products in the moth midge Clogmia albipunctata, a non-drosophilid dipteran species. We have used these antibodies to create a quantitative atlas of protein expression patterns for the gap gene hunchback (hb), and the pair-rule gene even-skipped (eve). Our data reveal differences in the dynamics of Hb boundary positioning and Eve stripe formation between C. albipunctata and Drosophila melanogaster. Despite these differences, the overall relative spatial arrangement of Hb and Eve domains is remarkably conserved between these two distantly related dipteran species. Conclusions We provide a proof of principle that it is possible to acquire quantitative gene expression data at high accuracy and spatio-temporal resolution in non-model organisms. Our quantitative data extend earlier qualitative studies of segmentation gene expression in C. albipunctata, and provide a starting point for comparative reverse-engineering studies of the evolutionary and developmental dynamics of the segmentation gene system.
Collapse
Affiliation(s)
- Hilde Janssens
- EMBL/CRG Research Unit in Systems Biology, Centre de Regulació Genòmica (CRG), and Universitat Pompeu Fabra (UPF), Dr. Aiguader 88, 08003 Barcelona, Spain
| | - Ken Siggens
- Department of Zoology, Downing Street, Cambridge CB2 3EJ UK
| | - Damjan Cicin-Sain
- EMBL/CRG Research Unit in Systems Biology, Centre de Regulació Genòmica (CRG), and Universitat Pompeu Fabra (UPF), Dr. Aiguader 88, 08003 Barcelona, Spain
| | - Eva Jiménez-Guri
- EMBL/CRG Research Unit in Systems Biology, Centre de Regulació Genòmica (CRG), and Universitat Pompeu Fabra (UPF), Dr. Aiguader 88, 08003 Barcelona, Spain
| | - Marco Musy
- EMBL/CRG Research Unit in Systems Biology, Centre de Regulació Genòmica (CRG), and Universitat Pompeu Fabra (UPF), Dr. Aiguader 88, 08003 Barcelona, Spain
| | - Michael Akam
- Department of Zoology, Downing Street, Cambridge CB2 3EJ UK
| | - Johannes Jaeger
- EMBL/CRG Research Unit in Systems Biology, Centre de Regulació Genòmica (CRG), and Universitat Pompeu Fabra (UPF), Dr. Aiguader 88, 08003 Barcelona, Spain
| |
Collapse
|