1
|
Li S, Liu Q, Wang E, Wang J. Global quantitative understanding of non-equilibrium cell fate decision-making in response to pheromone. iScience 2023; 26:107885. [PMID: 37766979 PMCID: PMC10520453 DOI: 10.1016/j.isci.2023.107885] [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: 05/23/2023] [Revised: 08/09/2023] [Accepted: 09/07/2023] [Indexed: 09/29/2023] Open
Abstract
Cell-cycle arrest and polarized growth are commonly used to characterize the response of yeast to pheromone. However, the quantitative decision-making processes underlying time-dependent changes in cell fate remain unclear. In this study, we conducted single-cell level experiments to observe multidimensional responses, uncovering diverse fates of yeast cells. Multiple states are revealed, along with the kinetic switching rates and pathways among them, giving rise to a quantitative landscape of mating response. To quantify the experimentally observed cell fates, we developed a theoretical framework based on non-equilibrium landscape and flux theory. Additionally, we performed stochastic simulations of biochemical reactions to elucidate signal transduction and cell growth. Notably, our experimental findings have provided the first global quantitative evidence of the real-time synchronization between intracellular signaling, physiological growth, and morphological functions. These results validate the proposed underlying mechanism governing the emergence of multiple cell fate states. This study introduces an emerging mechanistic approach to understand non-equilibrium cell fate decision-making in response to pheromone.
Collapse
Affiliation(s)
- Sheng Li
- College of Chemistry, Jilin University, Changchun, Jilin 130012, China
- State Key Laboratory of Electroanalytical Chemistry, Changchun Institute of Applied Chemistry, Chinese Academy of Sciences, Changchun, Jilin 130022, China
| | - Qiong Liu
- State Key Laboratory of Electroanalytical Chemistry, Changchun Institute of Applied Chemistry, Chinese Academy of Sciences, Changchun, Jilin 130022, China
| | - Erkang Wang
- College of Chemistry, Jilin University, Changchun, Jilin 130012, China
- State Key Laboratory of Electroanalytical Chemistry, Changchun Institute of Applied Chemistry, Chinese Academy of Sciences, Changchun, Jilin 130022, China
| | - Jin Wang
- Department of Chemistry and of Physics and Astronomy, State University of New York at Stony Brook, Stony Brook, NY 11794-3400, USA
| |
Collapse
|
2
|
Son D, Kim J. Estimation of Ordinary Differential Equation Models for Gene Regulatory Networks Through Data Cloning. J Comput Biol 2023; 30:609-618. [PMID: 36898058 DOI: 10.1089/cmb.2022.0201] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/12/2023] Open
Abstract
Ordinary differential equations (ODEs) are widely used for elucidating dynamic processes in various fields. One of the applications of ODEs is to describe dynamics of gene regulatory networks (GRNs), which is a critical step in understanding disease mechanisms. However, estimation of ODE models for GRNs is challenging because of inflexibility of the model and noisy data with complex error structures such as heteroscedasticity, correlations between genes, and time dependency. In addition, either a likelihood or Bayesian approach is commonly used for estimation of ODE models, but both approaches have benefits and drawbacks in their own right. Data cloning is a maximum likelihood (ML) estimation method through the Bayesian framework. Since it works in the Bayesian framework, it is free from local optimum problems that are common drawbacks of ML methods. Also, its inference is invariant for the selection of prior distributions, which is a major issue in Bayesian methods. This study proposes an estimation method of ODE models for GRNs through data cloning. The proposed method is demonstrated through simulation and it is applied to real gene expression time-course data.
Collapse
Affiliation(s)
- Donghui Son
- Department of Statistics, Sungkyunkwan University, Seoul, South Korea
| | - Jaejik Kim
- Department of Statistics, Sungkyunkwan University, Seoul, South Korea
| |
Collapse
|
3
|
Bajiya N, Dhall A, Aggarwal S, Raghava GPS. Advances in the field of phage-based therapy with special emphasis on computational resources. Brief Bioinform 2023; 24:6961791. [PMID: 36575815 DOI: 10.1093/bib/bbac574] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/28/2022] [Revised: 11/07/2022] [Accepted: 11/25/2022] [Indexed: 12/29/2022] Open
Abstract
In the current era, one of the major challenges is to manage the treatment of drug/antibiotic-resistant strains of bacteria. Phage therapy, a century-old technique, may serve as an alternative to antibiotics in treating bacterial infections caused by drug-resistant strains of bacteria. In this review, a systematic attempt has been made to summarize phage-based therapy in depth. This review has been divided into the following two sections: general information and computer-aided phage therapy (CAPT). In the case of general information, we cover the history of phage therapy, the mechanism of action, the status of phage-based products (approved and clinical trials) and the challenges. This review emphasizes CAPT, where we have covered primary phage-associated resources, phage prediction methods and pipelines. This review covers a wide range of databases and resources, including viral genomes and proteins, phage receptors, host genomes of phages, phage-host interactions and lytic proteins. In the post-genomic era, identifying the most suitable phage for lysing a drug-resistant strain of bacterium is crucial for developing alternate treatments for drug-resistant bacteria and this remains a challenging problem. Thus, we compile all phage-associated prediction methods that include the prediction of phages for a bacterial strain, the host for a phage and the identification of interacting phage-host pairs. Most of these methods have been developed using machine learning and deep learning techniques. This review also discussed recent advances in the field of CAPT, where we briefly describe computational tools available for predicting phage virions, the life cycle of phages and prophage identification. Finally, we describe phage-based therapy's advantages, challenges and opportunities.
Collapse
Affiliation(s)
- Nisha Bajiya
- Department of Computational Biology, Indraprastha Institute of Information Technology, Okhla Phase 3, New Delhi, 110020, India
| | - Anjali Dhall
- Department of Computational Biology, Indraprastha Institute of Information Technology, Okhla Phase 3, New Delhi, 110020, India
| | - Suchet Aggarwal
- Department of Computer Science and Engineering, Indraprastha Institute of Information Technology, Okhla Phase 3, New Delhi, 110020, India
| | - Gajendra P S Raghava
- Department of Computational Biology, Indraprastha Institute of Information Technology, Okhla Phase 3, New Delhi, 110020, India
| |
Collapse
|
4
|
Cockerell A, Wright L, Dattani A, Guo G, Smith A, Tsaneva-Atanasova K, Richards DM. Biophysical models of early mammalian embryogenesis. Stem Cell Reports 2023; 18:26-46. [PMID: 36630902 PMCID: PMC9860129 DOI: 10.1016/j.stemcr.2022.11.021] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/15/2022] [Revised: 11/02/2022] [Accepted: 11/24/2022] [Indexed: 01/12/2023] Open
Abstract
Embryo development is a critical and fascinating stage in the life cycle of many organisms. Despite decades of research, the earliest stages of mammalian embryogenesis are still poorly understood, caused by a scarcity of high-resolution spatial and temporal data, the use of only a few model organisms, and a paucity of truly multidisciplinary approaches that combine biological research with biophysical modeling and computational simulation. Here, we explain the theoretical frameworks and biophysical processes that are best suited to modeling the early mammalian embryo, review a comprehensive list of previous models, and discuss the most promising avenues for future work.
Collapse
Affiliation(s)
- Alaina Cockerell
- Living Systems Institute, University of Exeter, Stocker Road, Exeter EX4 4QD, UK
| | - Liam Wright
- Department of Mathematics, University of Exeter, North Park Road, Exeter EX4 4QF, UK
| | - Anish Dattani
- Living Systems Institute, University of Exeter, Stocker Road, Exeter EX4 4QD, UK
| | - Ge Guo
- Living Systems Institute, University of Exeter, Stocker Road, Exeter EX4 4QD, UK
| | - Austin Smith
- Living Systems Institute, University of Exeter, Stocker Road, Exeter EX4 4QD, UK
| | - Krasimira Tsaneva-Atanasova
- Living Systems Institute, University of Exeter, Stocker Road, Exeter EX4 4QD, UK; Department of Mathematics, University of Exeter, North Park Road, Exeter EX4 4QF, UK; EPSRC Hub for Quantitative Modelling in Healthcare, University of Exeter, Exeter EX4 4QJ, UK; Department of Bioinformatics and Mathematical Modelling, Institute of Biophysics and Biomedical Engineering, Bulgarian Academy of Sciences, 105 Acad. G. Bonchev Street, 1113 Sofia, Bulgaria
| | - David M Richards
- Living Systems Institute, University of Exeter, Stocker Road, Exeter EX4 4QD, UK; Department of Physics and Astronomy, University of Exeter, North Park Road, Exeter EX4 4QL, UK.
| |
Collapse
|
5
|
Ling MY, Chiu LJ, Hsieh CC, Shu CC. Dimerization induces bimodality in protein number distributions. Biosystems 2023; 223:104812. [PMID: 36427705 DOI: 10.1016/j.biosystems.2022.104812] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/27/2022] [Revised: 11/10/2022] [Accepted: 11/10/2022] [Indexed: 11/26/2022]
Abstract
We examined gene expression with DNA switching between two states, active and inactive. Subpopulations emerge from mechanisms that do not arise from trivial transcriptional heterogeneity. Although the RNA demonstrates a unimodal distribution, dimerization intriguingly causes protein bimodality. No control loop or deterministic bistability are present. In such a situation, increasing the degradation rate of the protein does not lead to bimodality. The bimodality is achieved through the interplay between the protein monomer and the formation of protein dimer. We applied Stochastic Simulation Algorithm (SSA) and found that cells spontaneously change states at the protein level. While sweeping parameters, decreasing the rate constant of dimerization severely impairs the bimodality. We also examined the influence of DNA switching. To have bimodality, the system requires a proper ratio of DNA in the active state to the inactive state. In addition to bimodality of the monomer, tetramerization also causes the bimodality of the dimer.
Collapse
Affiliation(s)
- Ming-Yang Ling
- Department of Chemical Engineering and Biotechnology, National Taipei University of Technology, Taipei City, Taiwan
| | - Lin-Jie Chiu
- Department of Chemical Engineering and Biotechnology, National Taipei University of Technology, Taipei City, Taiwan
| | - Ching-Chu Hsieh
- Department of Chemical Engineering and Biotechnology, National Taipei University of Technology, Taipei City, Taiwan
| | - Che-Chi Shu
- Department of Chemical Engineering and Biotechnology, National Taipei University of Technology, Taipei City, Taiwan.
| |
Collapse
|
6
|
Jia C, Grima R. Coupling gene expression dynamics to cell size dynamics and cell cycle events: Exact and approximate solutions of the extended telegraph model. iScience 2022; 26:105746. [PMID: 36619980 PMCID: PMC9813732 DOI: 10.1016/j.isci.2022.105746] [Citation(s) in RCA: 9] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/17/2022] [Revised: 11/02/2022] [Accepted: 12/02/2022] [Indexed: 12/12/2022] Open
Abstract
The standard model describing the fluctuations of mRNA numbers in single cells is the telegraph model which includes synthesis and degradation of mRNA, and switching of the gene between active and inactive states. While commonly used, this model does not describe how fluctuations are influenced by the cell cycle phase, cellular growth and division, and other crucial aspects of cellular biology. Here, we derive the analytical time-dependent solution of an extended telegraph model that explicitly considers the doubling of gene copy numbers upon DNA replication, dependence of the mRNA synthesis rate on cellular volume, gene dosage compensation, partitioning of molecules during cell division, cell-cycle duration variability, and cell-size control strategies. Based on the time-dependent solution, we obtain the analytical distributions of transcript numbers for lineage and population measurements in steady-state growth and also find a linear relation between the Fano factor of mRNA fluctuations and cell volume fluctuations. We show that generally the lineage and population distributions in steady-state growth cannot be accurately approximated by the steady-state solution of extrinsic noise models, i.e. a telegraph model with parameters drawn from probability distributions. This is because the mRNA lifetime is often not small enough compared to the cell cycle duration to erase the memory of division and replication. Accurate approximations are possible when this memory is weak, e.g. for genes with bursty expression and for which there is sufficient gene dosage compensation when replication occurs.
Collapse
Affiliation(s)
- Chen Jia
- Applied and Computational Mathematics Division, Beijing Computational Science Research Center, Beijing 100193, China
| | - Ramon Grima
- School of Biological Sciences, University of Edinburgh, Edinburgh EH9 3JH, UK,Corresponding author
| |
Collapse
|
7
|
Zhu H, Shen W, Luo C, Liu F. An integrated microfluidic device for multiplexed imaging of spatial gene expression patterns of Drosophila embryos. LAB ON A CHIP 2022; 22:4081-4092. [PMID: 36165088 DOI: 10.1039/d2lc00514j] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/16/2023]
Abstract
To reveal the underlying mechanism of the biological function of multicellular systems, it is important to obtain comprehensive spatial gene expression profiles. Among the emerging single-cell spatial-omics techniques, immunofluorescence (IF)-based iterative multiplexed imaging is a promising approach. However, the conventional method is usually costly, time-consuming, labor-intensive, and has low throughput. Moreover, it has yet to be demonstrated in intact multicellular organisms. Here, we developed an integrated microfluidic system to overcome these challenges for quantitatively measuring multiple protein profiles sequentially in situ in the same Drosophila embryo. We designed an array of hydrodynamic trapping sites to automatically capture over ten Drosophila embryos with orientation selectivity at more than 90% trapping rates. We also optimized the geometry of confinement and the on-chip IF protocol to achieve the same high signal-to-noise ratio as the off-chip traditional IF experiments. Moreover, we developed an efficient de-staining protocol by combining on-chip antibody stripping and fluorophore bleaching. Using the same secondary antibody to sequentially stain different genes, we confirmed that the de-stained genes have no detectable interference with the subsequently stained genes, and the gene expression profiles are preserved after multiple cycles of staining and de-staining processes. This preliminary test shows that our newly developed integrated microfluidic system can be a powerful tool for multiplexed imaging of Drosophila embryos. Our work opens a new avenue to design microfluidic chips for multicellular organisms and single-cell spatial-omics techniques.
Collapse
Affiliation(s)
- Hongcun Zhu
- The State Key Laboratory for Artificial Microstructures and Mesoscopic Physics, School of Physics, Peking University, Beijing, 100871, China.
| | - Wenting Shen
- Center for Quantitative Biology, Peking University, Beijing, 100871, China.
| | - Chunxiong Luo
- The State Key Laboratory for Artificial Microstructures and Mesoscopic Physics, School of Physics, Peking University, Beijing, 100871, China.
- Center for Quantitative Biology, Peking University, Beijing, 100871, China.
- Wenzhou Institute University of Chinese Academy of Sciences, Wenzhou, Zhejiang, China
| | - Feng Liu
- Center for Quantitative Biology, Peking University, Beijing, 100871, China.
- Key Laboratory of Hebei Province for Molecular Biophysics, Institute of Biophysics, School of Health Science & Biomedical Engineering, Hebei University of Technology, Tianjin, 300130, China
| |
Collapse
|
8
|
Concentration fluctuations in growing and dividing cells: Insights into the emergence of concentration homeostasis. PLoS Comput Biol 2022; 18:e1010574. [PMID: 36194626 PMCID: PMC9565450 DOI: 10.1371/journal.pcbi.1010574] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/07/2022] [Revised: 10/14/2022] [Accepted: 09/14/2022] [Indexed: 11/19/2022] Open
Abstract
Intracellular reaction rates depend on concentrations and hence their levels are often regulated. However classical models of stochastic gene expression lack a cell size description and cannot be used to predict noise in concentrations. Here, we construct a model of gene product dynamics that includes a description of cell growth, cell division, size-dependent gene expression, gene dosage compensation, and size control mechanisms that can vary with the cell cycle phase. We obtain expressions for the approximate distributions and power spectra of concentration fluctuations which lead to insight into the emergence of concentration homeostasis. We find that (i) the conditions necessary to suppress cell division-induced concentration oscillations are difficult to achieve; (ii) mRNA concentration and number distributions can have different number of modes; (iii) two-layer size control strategies such as sizer-timer or adder-timer are ideal because they maintain constant mean concentrations whilst minimising concentration noise; (iv) accurate concentration homeostasis requires a fine tuning of dosage compensation, replication timing, and size-dependent gene expression; (v) deviations from perfect concentration homeostasis show up as deviations of the concentration distribution from a gamma distribution. Some of these predictions are confirmed using data for E. coli, fission yeast, and budding yeast.
Collapse
|
9
|
Network Approaches for Charting the Transcriptomic and Epigenetic Landscape of the Developmental Origins of Health and Disease. Genes (Basel) 2022; 13:genes13050764. [PMID: 35627149 PMCID: PMC9141211 DOI: 10.3390/genes13050764] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/04/2022] [Revised: 04/04/2022] [Accepted: 04/13/2022] [Indexed: 02/04/2023] Open
Abstract
The early developmental phase is of critical importance for human health and disease later in life. To decipher the molecular mechanisms at play, current biomedical research is increasingly relying on large quantities of diverse omics data. The integration and interpretation of the different datasets pose a critical challenge towards the holistic understanding of the complex biological processes that are involved in early development. In this review, we outline the major transcriptomic and epigenetic processes and the respective datasets that are most relevant for studying the periconceptional period. We cover both basic data processing and analysis steps, as well as more advanced data integration methods. A particular focus is given to network-based methods. Finally, we review the medical applications of such integrative analyses.
Collapse
|
10
|
Pilati S, Malacarne G, Navarro-Payá D, Tomè G, Riscica L, Cavecchia V, Matus JT, Moser C, Blanzieri E. Vitis OneGenE: A Causality-Based Approach to Generate Gene Networks in Vitis vinifera Sheds Light on the Laccase and Dirigent Gene Families. Biomolecules 2021; 11:1744. [PMID: 34944388 PMCID: PMC8698957 DOI: 10.3390/biom11121744] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/09/2021] [Revised: 11/11/2021] [Accepted: 11/16/2021] [Indexed: 12/24/2022] Open
Abstract
The abundance of transcriptomic data and the development of causal inference methods have paved the way for gene network analyses in grapevine. Vitis OneGenE is a transcriptomic data mining tool that finds direct correlations between genes, thus producing association networks. As a proof of concept, the stilbene synthase gene regulatory network obtained with OneGenE has been compared with published co-expression analysis and experimental data, including cistrome data for MYB stilbenoid regulators. As a case study, the two secondary metabolism pathways of stilbenoids and lignin synthesis were explored. Several isoforms of laccase, peroxidase, and dirigent protein genes, putatively involved in the final oxidative oligomerization steps, were identified as specifically belonging to either one of these pathways. Manual curation of the predicted sequences exploiting the last available genome assembly, and the integration of phylogenetic and OneGenE analyses, identified a group of laccases exclusively present in grapevine and related to stilbenoids. Here we show how network analysis by OneGenE can accelerate knowledge discovery by suggesting new candidates for functional characterization and application in breeding programs.
Collapse
Affiliation(s)
- Stefania Pilati
- Research and Innovation Centre, Department of Genomics and Biology of Fruit Crops, Fondazione Edmund Mach, 38098 San Michele all’Adige, Italy; (G.M.); (C.M.)
| | - Giulia Malacarne
- Research and Innovation Centre, Department of Genomics and Biology of Fruit Crops, Fondazione Edmund Mach, 38098 San Michele all’Adige, Italy; (G.M.); (C.M.)
| | - David Navarro-Payá
- Institute for Integrative Systems Biology (I2SysBio), Universitat de València-CSIC, 46908 Paterna, Valencia, Spain; (D.N.-P.); (J.T.M.)
| | - Gabriele Tomè
- Centre for Integrative Biology (CIBIO), University of Trento, 38123 Trento, Italy;
| | - Laura Riscica
- Department of Information Engineering and Computer Science, University of Trento, 38123 Trento, Italy; (L.R.); (E.B.)
| | - Valter Cavecchia
- CNR-Institute of Materials for Electronics and Magnetism, 38123 Trento, Italy;
| | - José Tomás Matus
- Institute for Integrative Systems Biology (I2SysBio), Universitat de València-CSIC, 46908 Paterna, Valencia, Spain; (D.N.-P.); (J.T.M.)
| | - Claudio Moser
- Research and Innovation Centre, Department of Genomics and Biology of Fruit Crops, Fondazione Edmund Mach, 38098 San Michele all’Adige, Italy; (G.M.); (C.M.)
| | - Enrico Blanzieri
- Department of Information Engineering and Computer Science, University of Trento, 38123 Trento, Italy; (L.R.); (E.B.)
- CNR-Institute of Materials for Electronics and Magnetism, 38123 Trento, Italy;
| |
Collapse
|
11
|
Qin BW, Zhao L, Lin W. A frequency-amplitude coordinator and its optimal energy consumption for biological oscillators. Nat Commun 2021; 12:5894. [PMID: 34625549 PMCID: PMC8501100 DOI: 10.1038/s41467-021-26182-2] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/01/2021] [Accepted: 09/22/2021] [Indexed: 02/08/2023] Open
Abstract
Biorhythm including neuron firing and protein-mRNA interaction are fundamental activities with diffusive effect. Their well-balanced spatiotemporal dynamics are beneficial for healthy sustainability. Therefore, calibrating both anomalous frequency and amplitude of biorhythm prevents physiological dysfunctions or diseases. However, many works were devoted to modulate frequency exclusively whereas amplitude is usually ignored, although both quantities are equally significant for coordinating biological functions and outputs. Especially, a feasible method coordinating the two quantities concurrently and precisely is still lacking. Here, for the first time, we propose a universal approach to design a frequency-amplitude coordinator rigorously via dynamical systems tools. We consider both spatial and temporal information. With a single well-designed coordinator, they can be calibrated to desired levels simultaneously and precisely. The practical usefulness and efficacy of our method are demonstrated in representative neuronal and gene regulatory models. We further reveal its fundamental mechanism and optimal energy consumption providing inspiration for biorhythm regulation in future.
Collapse
Affiliation(s)
- Bo-Wei Qin
- School of Mathematical Sciences, Fudan University, 200433, Shanghai, China.
- State Key Laboratory of Medical Neurobiology and MOE Frontiers Center for Brain Science, Institutes of Brain Science, Fudan University, 200032, Shanghai, China.
| | - Lei Zhao
- School of Mathematical Sciences, Fudan University, 200433, Shanghai, China
- The GLOBE Institute, University of Copenhagen, Copenhagen, Denmark
| | - Wei Lin
- School of Mathematical Sciences, Fudan University, 200433, Shanghai, China.
- State Key Laboratory of Medical Neurobiology and MOE Frontiers Center for Brain Science, Institutes of Brain Science, Fudan University, 200032, Shanghai, China.
- Shanghai Center for Mathematical Sciences, 200438, Shanghai, China.
- Center for Computational Systems Biology of ISTBI, LCNBI, and Research Institute of Intelligent Complex Systems, Fudan University, 200433, Shanghai, China.
| |
Collapse
|
12
|
Mathematical Modeling of ceRNA-Based Interactions. Methods Mol Biol 2021. [PMID: 34165711 DOI: 10.1007/978-1-0716-1503-4_7] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/31/2023]
Abstract
Competing endogenous RNA (ceRNA) molecules have emerged as key players in regulating gene expression, increasing the complexity of the range of possible dynamics within a cell. The actions of competing RNA typically are sponging behaviors, in a manner that fine-tunes gene expression, but there are particular network structures that may show destabilization due to ceRNA interactions. In this chapter, we discuss how these interactions can be modeled and probed from a mathematical, first-principles perspective.
Collapse
|
13
|
Uriu K, Tei H. Complementary phase responses via functional differentiation of dual negative feedback loops. PLoS Comput Biol 2021; 17:e1008774. [PMID: 33684114 PMCID: PMC7971863 DOI: 10.1371/journal.pcbi.1008774] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/18/2020] [Revised: 03/18/2021] [Accepted: 02/05/2021] [Indexed: 11/18/2022] Open
Abstract
Multiple feedback loops are often found in gene regulations for various cellular functions. In mammalian circadian clocks, oscillations of Period1 (Per1) and Period2 (Per2) expression are caused by interacting negative feedback loops (NFLs) whose protein products with similar molecular functions repress each other. However, Per1 expression peaks earlier than Per2 in the pacemaker tissue, raising the question of whether the peak time difference reflects their different dynamical functions. Here, we address this question by analyzing phase responses of the circadian clock caused by light-induced transcription of both Per1 and Per2 mRNAs. Through mathematical analyses of dual NFLs, we show that phase advance is mainly driven by light inputs to the repressor with an earlier expression peak as Per1, whereas phase delay is driven by the other repressor with a later peak as Per2. Due to the complementary contributions to phase responses, the ratio of light-induced transcription rates between Per1 and Per2 determines the magnitude and direction of phase shifts at each time of day. Specifically, stronger Per1 light induction than Per2 results in a phase response curve (PRC) with a larger phase advance zone than delay zone as observed in rats and hamsters, whereas stronger Per2 induction causes a larger delay zone as observed in mice. Furthermore, the ratio of light-induced transcription rates required for entrainment is determined by the relation between the circadian and light-dark periods. Namely, if the autonomous period of a circadian clock is longer than the light-dark period, a larger light-induced transcription rate of Per1 than Per2 is required for entrainment, and vice versa. In short, the time difference between Per1 and Per2 expression peaks can differentiate their dynamical functions. The resultant complementary contributions to phase responses can determine entrainability of the circadian clock to the light-dark cycle.
Collapse
Affiliation(s)
- Koichiro Uriu
- Graduate School of Natural Science and Technology, Kanazawa University, Kanazawa, Japan
- * E-mail:
| | - Hajime Tei
- Graduate School of Natural Science and Technology, Kanazawa University, Kanazawa, Japan
| |
Collapse
|
14
|
Hsu IS, Moses AM. Stochastic models for single-cell data: Current challenges and the way forward. FEBS J 2021; 289:647-658. [PMID: 33570798 DOI: 10.1111/febs.15760] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/30/2020] [Revised: 12/22/2020] [Accepted: 02/10/2021] [Indexed: 11/28/2022]
Abstract
Although the quantity and quality of single-cell data have progressed rapidly, making quantitative predictions with single-cell stochastic models remains challenging. The stochastic nature of cellular processes leads to at least three challenges in building models with single-cell data: (a) because variability in single-cell data can be attributed to multiple different sources, it is difficult to rule out conflicting mechanistic models that explain the same data equally well; (b) the distinction between interesting biological variability and experimental variability is sometimes ambiguous; (c) the nonstandard distributions of single-cell data can lead to violations of the assumption of symmetric errors in least-squares fitting. In this review, we first discuss recent studies that overcome some of the challenges or set up a promising direction and then introduce some powerful statistical approaches utilized in these studies. We conclude that applying and developing statistical approaches could lead to further progress in building stochastic models for single-cell data.
Collapse
Affiliation(s)
- Ian S Hsu
- Department of Cell & Systems Biology, University of Toronto, ON, Canada
| | - Alan M Moses
- Department of Cell & Systems Biology, University of Toronto, ON, Canada
| |
Collapse
|
15
|
Khazaaleh M, Samarasinghe S, Kulasiri D. A new hierarchical approach to multi-level model abstraction for simplifying ODE models of biological networks and a case study: The G1/S Checkpoint/DNA damage signalling pathways of mammalian cell cycle. Biosystems 2021; 203:104374. [PMID: 33556446 DOI: 10.1016/j.biosystems.2021.104374] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/11/2020] [Revised: 01/19/2021] [Accepted: 01/25/2021] [Indexed: 11/15/2022]
Abstract
Model reduction is an important topic in studies of biological systems. By reducing the complexity of large models through multi-level models while keeping the essence (biological meaning) of the model, model reduction can help answer many important questions about these systems. In this paper, we present a new reduction method based on hierarchical representation and a lumping approach. We used G1/S checkpoint pathway represented in Ordinary Differential Equations (ODE) in Iwamoto et al. (2011) as a case study to present this reduction method. The approach consists of two parts; the first part represents the biological network as a hierarchy (multiple levels) based on protein binding relations, which allowed us to model the biological network at different levels of abstraction. The second part applies different levels (level 1, 2 and 3) of lumping the species together depending on the level of the hierarchy, resulting in a reduced and transformed model for each level. The model at each level is a representation of the whole system and can address questions pertinent to that level. We develop and simulate reduced models for levels-1, 2 and 3 of lumping for the G1/S checkpoint pathway and evaluate the biological plausibility of the proposed method by comparing the results with the original ODE model of Iwamoto et al. (2011). The results for continuous dynamics of the G1/S checkpoint pathway with or without DNA-damage for reduced models of level- 1, 2 and 3 of lumping are in very good agreement and consistent with the original model results and with biological findings. Therefore, the reduced models (level-1, 2 and 3) can be used to study cell cycle progression in G1 and the effects of DNA damage on it. It is suitable for reducing complex ODE biological network models while retaining accurate continuous dynamics of the system. The 3 levels of the reduced models respectively achieved 20%, 26% and 31% reduction of the base model. Moreover, the reduced model is more efficient to run (30%, 44% and 52% time reduction for the three levels) and generate solutions than the original ODE model. Simplification of complex mathematical models is possible and the proposed reduction method has the potential to make an impact across many fields of biomedical research.
Collapse
Affiliation(s)
- Mutaz Khazaaleh
- Complex Systems, Big Data and Informatics Initiative (CSBII), New Zealand
| | - Sandhya Samarasinghe
- Complex Systems, Big Data and Informatics Initiative (CSBII), New Zealand; Centre for Advanced Computational Solutions, Lincoln University, Christchurch, New Zealand.
| | - Don Kulasiri
- Complex Systems, Big Data and Informatics Initiative (CSBII), New Zealand; Centre for Advanced Computational Solutions, Lincoln University, Christchurch, New Zealand
| |
Collapse
|
16
|
Stochastic models coupling gene expression and partitioning in cell division in Escherichia coli. Biosystems 2020; 193-194:104154. [PMID: 32353481 DOI: 10.1016/j.biosystems.2020.104154] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/10/2019] [Revised: 04/03/2020] [Accepted: 04/16/2020] [Indexed: 12/18/2022]
Abstract
Regulation of future RNA and protein numbers is a key process by which cells continuously best fit the environment. In bacteria, RNA and proteins exist in small numbers and their regulatory processes are stochastic. Consequently, there is cell-to-cell variability in these numbers, even between sister cells. Traditionally, the two most studied sources of this variability are gene expression and RNA and protein degradation, with evidence suggesting that the latter is subject to little regulation, when compared to the former. However, time-lapse microscopy and single molecule fluorescent tagging have produced evidence that cell division can also be a significant source of variability due to asymmetries in the partitioning of RNA and proteins. Relevantly, the impact of this noise differs from noise in production and degradation since, unlike these, it is not continuous. Rather, it occurs at specific time points, at which moment it can introduce major fluctuations. Several models have now been proposed that integrate noise from cell division, in addition to noise in gene expression, to mimic the dynamics of RNA and protein numbers of cell lineages. This is expected to be particularly relevant in genetic circuits, where significant fluctuations in one component protein, at specific time moments, are expected to perturb near-equilibrium states of the circuits, which can have long-lasting consequences. Here we review stochastic models coupling these processes in Escherichia coli, from single genes to small circuits.
Collapse
|
17
|
Ahmed SS, Roy S, Kalita J. Assessing the Effectiveness of Causality Inference Methods for Gene Regulatory Networks. IEEE/ACM TRANSACTIONS ON COMPUTATIONAL BIOLOGY AND BIOINFORMATICS 2020; 17:56-70. [PMID: 29994618 DOI: 10.1109/tcbb.2018.2853728] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/08/2023]
Abstract
Causality inference is the use of computational techniques to predict possible causal relationships for a set of variables, thereby forming a directed network. Causality inference in Gene Regulatory Networks (GRNs) is an important, yet challenging task due to the limits of available data and lack of efficiency in existing causality inference techniques. A number of techniques have been proposed and applied to infer causal relationships in various domains, although they are not specific to regulatory network inference. In this paper, we assess the effectiveness of methods for inferring causal GRNs. We introduce seven different inference methods and apply them to infer directed edges in GRNs. We use time-series expression data from the DREAM challenges to assess the methods in terms of quality of inference and rank them based on performance. The best method is applied to Breast Cancer data to infer a causal network. Experimental results show that Causation Entropy is best, however, highly time-consuming and not feasible to use in a relatively large network. We infer Breast Cancer GRN with the second-best method, Transfer Entropy. The topological analysis of the network reveals that top out-degree genes such as SLC39A5 which are considered central genes, play important role in cancer progression.
Collapse
|
18
|
Wang H, Li C, Zhang J, Wang J, Ma Y, Lian Y. A new LSTM-based gene expression prediction model: L-GEPM. J Bioinform Comput Biol 2019; 17:1950022. [PMID: 31617459 DOI: 10.1142/s0219720019500227] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/03/2023]
Abstract
Molecular biology combined with in silico machine learning and deep learning has facilitated the broad application of gene expression profiles for gene function prediction, optimal crop breeding, disease-related gene discovery, and drug screening. Although the acquisition cost of genome-wide expression profiles has been steadily declining, the requirement generates a compendium of expression profiles using thousands of samples remains high. The Library of Integrated Network-Based Cellular Signatures (LINCS) program used approximately 1000 landmark genes to predict the expression of the remaining target genes by linear regression; however, this approach ignored the nonlinear features influencing gene expression relationships, limiting the accuracy of the experimental results. We herein propose a gene expression prediction model, L-GEPM, based on long short-term memory (LSTM) neural networks, which captures the nonlinear features affecting gene expression and uses learned features to predict the target genes. By comparing and analyzing experimental errors and fitting the effects of different prediction models, the LSTM neural network-based model, L-GEPM, can achieve low error and a superior fitting effect.
Collapse
Affiliation(s)
- Huiqing Wang
- College of Information and Computer, Taiyuan University of Technology, P. R. China
| | - Chun Li
- College of Information and Computer, Taiyuan University of Technology, P. R. China
| | - Jianhui Zhang
- College of Information and Computer, Taiyuan University of Technology, P. R. China
| | - Jingjing Wang
- College of Information and Computer, Taiyuan University of Technology, P. R. China
| | - Yue Ma
- College of Information and Computer, Taiyuan University of Technology, P. R. China
| | - Yuanyuan Lian
- College of Information and Computer, Taiyuan University of Technology, P. R. China
| |
Collapse
|
19
|
Ahsen ME, Chun Y, Grishin A, Grishina G, Stolovitzky G, Pandey G, Bunyavanich S. NeTFactor, a framework for identifying transcriptional regulators of gene expression-based biomarkers. Sci Rep 2019; 9:12970. [PMID: 31506535 PMCID: PMC6737052 DOI: 10.1038/s41598-019-49498-y] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/31/2018] [Accepted: 08/27/2019] [Indexed: 12/21/2022] Open
Abstract
Biological and regulatory mechanisms underlying many multi-gene expression-based disease biomarkers are often not readily evident. We describe an innovative framework, NeTFactor, that combines network analyses with gene expression data to identify transcription factors (TFs) that significantly and maximally regulate such a biomarker. NeTFactor uses a computationally-inferred context-specific gene regulatory network and applies topological, statistical, and optimization methods to identify regulator TFs. Application of NeTFactor to a multi-gene expression-based asthma biomarker identified ETS translocation variant 4 (ETV4) and peroxisome proliferator-activated receptor gamma (PPARG) as the biomarker's most significant TF regulators. siRNA-based knock down of these TFs in an airway epithelial cell line model demonstrated significant reduction of cytokine expression relevant to asthma, validating NeTFactor's top-scoring findings. While PPARG has been associated with airway inflammation, ETV4 has not yet been implicated in asthma, thus indicating the possibility of novel, disease-relevant discovery by NeTFactor. We also show that NeTFactor's results are robust when the gene regulatory network and biomarker are derived from independent data. Additionally, our application of NeTFactor to a different disease biomarker identified TF regulators of interest. These results illustrate that the application of NeTFactor to multi-gene expression-based biomarkers could yield valuable insights into regulatory mechanisms and biological processes underlying disease.
Collapse
Affiliation(s)
- Mehmet Eren Ahsen
- Icahn Institute for Genomics and Multiscale Biology and Department of Genetics and Genomic Sciences, Icahn School of Medicine at Mount Sinai, New York, NY, USA
| | - Yoojin Chun
- Icahn Institute for Genomics and Multiscale Biology and Department of Genetics and Genomic Sciences, Icahn School of Medicine at Mount Sinai, New York, NY, USA
| | - Alexander Grishin
- Division of Allergy & Immunology, Department of Pediatrics, Icahn School of Medicine at Mount Sinai, New York, NY, USA
| | - Galina Grishina
- Division of Allergy & Immunology, Department of Pediatrics, Icahn School of Medicine at Mount Sinai, New York, NY, USA
| | - Gustavo Stolovitzky
- Icahn Institute for Genomics and Multiscale Biology and Department of Genetics and Genomic Sciences, Icahn School of Medicine at Mount Sinai, New York, NY, USA
- IBM T.J. Watson Research Center, Yorktown Heights, New York, NY, USA
| | - Gaurav Pandey
- Icahn Institute for Genomics and Multiscale Biology and Department of Genetics and Genomic Sciences, Icahn School of Medicine at Mount Sinai, New York, NY, USA.
| | - Supinda Bunyavanich
- Icahn Institute for Genomics and Multiscale Biology and Department of Genetics and Genomic Sciences, Icahn School of Medicine at Mount Sinai, New York, NY, USA.
- Division of Allergy & Immunology, Department of Pediatrics, Icahn School of Medicine at Mount Sinai, New York, NY, USA.
| |
Collapse
|
20
|
Yuan Y, Liu J, Zhao P, Xing F, Huo H, Fang T. Structural Insights Into the Dynamic Evolution of Neuronal Networks as Synaptic Density Decreases. Front Neurosci 2019; 13:892. [PMID: 31507365 PMCID: PMC6714520 DOI: 10.3389/fnins.2019.00892] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/08/2019] [Accepted: 08/08/2019] [Indexed: 11/13/2022] Open
Abstract
The human brain is thought to be an extremely complex but efficient computing engine, processing vast amounts of information from a changing world. The decline in the synaptic density of neuronal networks is one of the most important characteristics of brain development, which is closely related to synaptic pruning, synaptic growth, synaptic plasticity, and energy metabolism. However, because of technical limitations in observing large-scale neuronal networks dynamically connected through synapses, how neuronal networks are organized and evolve as their synaptic density declines remains unclear. Here, by establishing a biologically reasonable neuronal network model, we show that despite a decline in the synaptic density, the connectivity, and efficiency of neuronal networks can be improved. Importantly, by analyzing the degree distribution, we also find that both the scale-free characteristic of neuronal networks and the emergence of hub neurons rely on the spatial distance between neurons. These findings may promote our understanding of neuronal networks in the brain and have guiding significance for the design of neuronal network models.
Collapse
Affiliation(s)
- Ye Yuan
- Department of Automation, Shanghai Jiao Tong University, Shanghai, China.,Key Laboratory of System Control and Information Processing, Ministry of Education, Shanghai, China
| | - Jian Liu
- Department of Automation, Shanghai Jiao Tong University, Shanghai, China.,Key Laboratory of System Control and Information Processing, Ministry of Education, Shanghai, China
| | - Peng Zhao
- Department of Automation, Shanghai Jiao Tong University, Shanghai, China.,Key Laboratory of System Control and Information Processing, Ministry of Education, Shanghai, China
| | - Fu Xing
- Department of Automation, Shanghai Jiao Tong University, Shanghai, China.,Key Laboratory of System Control and Information Processing, Ministry of Education, Shanghai, China
| | - Hong Huo
- Department of Automation, Shanghai Jiao Tong University, Shanghai, China.,Key Laboratory of System Control and Information Processing, Ministry of Education, Shanghai, China
| | - Tao Fang
- Department of Automation, Shanghai Jiao Tong University, Shanghai, China.,Key Laboratory of System Control and Information Processing, Ministry of Education, Shanghai, China
| |
Collapse
|
21
|
Ainuddin U, Khurram M, Hasan SMR. Cloning the λ Switch: Digital and Markov Representations. IEEE Trans Nanobioscience 2019; 18:428-436. [PMID: 30946673 DOI: 10.1109/tnb.2019.2908669] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/08/2022]
Abstract
The lysis-lysogeny switch in E. coli due to infection from lambda phage has been extensively studied and explained by scientists of molecular biology. The bacterium either survives with the viral strand of deoxyribonucleic acid (DNA) or dies producing hundreds of viruses for propagation of infection. Many proteins transcribed after infection by λ phage take part in determining the fate of the bacterium, but two proteins that play a key role in this regard are the cI and cro dimers, which are transcribed off the viral DNA. This paper presents a novel modeling mechanism for the lysis-lysogeny switch, by transferring the interactions of the main proteins, the lambda right operator and promoter regions and the ribonucleic acid (RNA) polymerase, to a finite state machine (FSM), to determine cell fate. The FSM, and thus derived is implemented in field-programmable gate array (FPGA), and simulations have been run in random conditions. A Markov model has been created for the same mechanism. Steady state analysis has been conducted for the transition matrix of the Markov model, and the results have been generated to show the steady state probability of lysis with various model values. In this paper, it is hoped to lay down guidelines to convert biological processes into computing machines.
Collapse
|
22
|
Roy A, Klumpp S. Simulating Genetic Circuits in Bacterial Populations with Growth Heterogeneity. Biophys J 2019; 114:484-492. [PMID: 29401445 DOI: 10.1016/j.bpj.2017.11.3745] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/07/2017] [Revised: 10/31/2017] [Accepted: 11/17/2017] [Indexed: 01/21/2023] Open
Abstract
We computationally study genetic circuits in bacterial populations with heterogeneities in the growth rate. To that end, we present a stochastic simulation method for gene circuits in populations of cells and propose an efficient implementation that we call the "Next Family Method". Within this approach, we implement different population setups, specifically Chemostat-type growth and growth in an ideal Mother Machine and show that the population structure and its statistics are different for the different setups whenever there is growth heterogeneity. Such dependence on the population setup is demonstrated, in the case of bistable systems with different growth rates in the stable states, to have distinctive signatures on quantities including the distributions of protein concentration and growth rates, and hysteresis curves. Applying this method to a bistable antibiotic resistance circuit, we find that as a result of the different statistics in different population setups, the estimated minimal inhibitory concentration of the antibiotic becomes dependent on the population setup in which it is measured.
Collapse
Affiliation(s)
- Anjan Roy
- Max Planck Institute of Colloids and Interfaces, Potsdam, Germany; Institute for Nonlinear Dynamics, University of Göttingen, Göttingen, Germany
| | - Stefan Klumpp
- Max Planck Institute of Colloids and Interfaces, Potsdam, Germany; Institute for Nonlinear Dynamics, University of Göttingen, Göttingen, Germany.
| |
Collapse
|
23
|
Cussat-Blanc S, Harrington K, Banzhaf W. Artificial Gene Regulatory Networks-A Review. ARTIFICIAL LIFE 2019; 24:296-328. [PMID: 30681915 DOI: 10.1162/artl_a_00267] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/09/2023]
Abstract
In nature, gene regulatory networks are a key mediator between the information stored in the DNA of living organisms (their genotype) and the structural and behavioral expression this finds in their bodies, surviving in the world (their phenotype). They integrate environmental signals, steer development, buffer stochasticity, and allow evolution to proceed. In engineering, modeling and implementations of artificial gene regulatory networks have been an expanding field of research and development over the past few decades. This review discusses the concept of gene regulation, describes the current state of the art in gene regulatory networks, including modeling and simulation, and reviews their use in artificial evolutionary settings. We provide evidence for the benefits of this concept in natural and the engineering domains.
Collapse
Affiliation(s)
| | - Kyle Harrington
- University of Idaho, Computational and Physical Systems Group, Virtual Technology and Design.
| | - Wolfgang Banzhaf
- Michigan State University, BEACON Center for the Study of Evolution in Action, Department of Computer Science and Engineering.
| |
Collapse
|
24
|
Turner A, Tyrrell A, Trefzer M, Lones M. Evolutionary acquisition of complex traits in artificial epigenetic networks. Biosystems 2018; 176:17-26. [PMID: 30557598 DOI: 10.1016/j.biosystems.2018.12.001] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/16/2018] [Revised: 11/30/2018] [Accepted: 12/03/2018] [Indexed: 11/29/2022]
Abstract
How complex traits arise within organisms over evolutionary time is an important question that has relevance both to the understanding of biological systems and to the design of bio-inspired computing systems. This paper investigates the process of acquiring complex traits within epiNet, a recurrent connectionist architecture capable of adapting its topology during execution. Inspired by the biological processes of gene regulation and epigenetics, epiNet captures biological organisms' ability to alter their regulatory topologies according to environmental stimulus. By applying epiNet to a series of computational tasks, each requiring a range of complex behaviours to solve, and capturing the evolutionary process in detail, we can show not only how the physical structure of epiNet changed when acquiring complex traits, but also how these changes in physical structure affected its dynamic behaviour. This is facilitated by using a lightweight optimisation method which makes minor iterative changes to the network structure so that when complex traits emerge for the first time, a direct lineage can be observed detailing exactly how they evolved. From this we can build an understanding of how complex traits evolve and which regulatory environments best allow for the emergence of these complex traits, pointing us towards computational models that allow more swift and robust acquisition of complex traits when optimised in an evolutionary computing setting.
Collapse
Affiliation(s)
| | - Andy Tyrrell
- Department of Electronic Engineering, University of York, UK
| | - Martin Trefzer
- Department of Electronic Engineering, University of York, UK
| | - Michael Lones
- School of Mathematical and Computer Sciences, Heriot-Watt University, UK
| |
Collapse
|
25
|
Vardi N, Chaturvedi S, Weinberger LS. Feedback-mediated signal conversion promotes viral fitness. Proc Natl Acad Sci U S A 2018; 115:E8803-E8810. [PMID: 30150412 PMCID: PMC6140503 DOI: 10.1073/pnas.1802905115] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
A fundamental signal-processing problem is how biological systems maintain phenotypic states (i.e., canalization) long after degradation of initial catalyst signals. For example, to efficiently replicate, herpesviruses (e.g., human cytomegalovirus, HCMV) rapidly counteract cell-mediated silencing using transactivators packaged in the tegument of the infecting virion particle. However, the activity of these tegument transactivators is inherently transient-they undergo immediate proteolysis but delayed synthesis-and how transient activation sustains lytic viral gene expression despite cell-mediated silencing is unclear. By constructing a two-color, conditional-feedback HCMV mutant, we find that positive feedback in HCMV's immediate-early 1 (IE1) protein is of sufficient strength to sustain HCMV lytic expression. Single-cell time-lapse imaging and mathematical modeling show that IE1 positive feedback converts transient transactivation signals from tegument pp71 proteins into sustained lytic expression, which is obligate for efficient viral replication, whereas attenuating feedback decreases fitness by promoting a reversible silenced state. Together, these results identify a regulatory mechanism enabling herpesviruses to sustain expression despite transient activation signals-akin to early electronic transistors-and expose a potential target for therapeutic intervention.
Collapse
Affiliation(s)
- Noam Vardi
- Gladstone-University of California, San Francisco (UCSF) Center for Cell Circuitry, Gladstone Institutes, San Francisco, CA 94158
| | - Sonali Chaturvedi
- Gladstone-University of California, San Francisco (UCSF) Center for Cell Circuitry, Gladstone Institutes, San Francisco, CA 94158
| | - Leor S Weinberger
- Gladstone-University of California, San Francisco (UCSF) Center for Cell Circuitry, Gladstone Institutes, San Francisco, CA 94158;
- Department of Biochemistry and Biophysics, University of California, San Francisco, CA 94158
- Department of Pharmaceutical Chemistry, University of California, San Francisco, CA 94158
| |
Collapse
|
26
|
Vlysidis M, Kaznessis YN. On Differences between Deterministic and Stochastic Models of Chemical Reactions: Schlögl Solved with ZI-Closure. ENTROPY 2018; 20:e20090678. [PMID: 33265767 PMCID: PMC7513203 DOI: 10.3390/e20090678] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 07/03/2018] [Revised: 08/29/2018] [Accepted: 09/04/2018] [Indexed: 11/16/2022]
Abstract
Deterministic and stochastic models of chemical reaction kinetics can give starkly different results when the deterministic model exhibits more than one stable solution. For example, in the stochastic Schlögl model, the bimodal stationary probability distribution collapses to a unimodal distribution when the system size increases, even for kinetic constant values that result in two distinct stable solutions in the deterministic Schlögl model. Using zero-information (ZI) closure scheme, an algorithm for solving chemical master equations, we compute stationary probability distributions for varying system sizes of the Schlögl model. With ZI-closure, system sizes can be studied that have been previously unattainable by stochastic simulation algorithms. We observe and quantify paradoxical discrepancies between stochastic and deterministic models and explain this behavior by postulating that the entropy of non-equilibrium steady states (NESS) is maximum.
Collapse
|
27
|
Khaluf Y, Ferrante E, Simoens P, Huepe C. Scale invariance in natural and artificial collective systems: a review. J R Soc Interface 2018; 14:rsif.2017.0662. [PMID: 29093130 DOI: 10.1098/rsif.2017.0662] [Citation(s) in RCA: 33] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/08/2017] [Accepted: 10/09/2017] [Indexed: 01/10/2023] Open
Abstract
Self-organized collective coordinated behaviour is an impressive phenomenon, observed in a variety of natural and artificial systems, in which coherent global structures or dynamics emerge from local interactions between individual parts. If the degree of collective integration of a system does not depend on size, its level of robustness and adaptivity is typically increased and we refer to it as scale-invariant. In this review, we first identify three main types of self-organized scale-invariant systems: scale-invariant spatial structures, scale-invariant topologies and scale-invariant dynamics. We then provide examples of scale invariance from different domains in science, describe their origins and main features and discuss potential challenges and approaches for designing and engineering artificial systems with scale-invariant properties.
Collapse
Affiliation(s)
- Yara Khaluf
- Ghent University-imec, IDLab-INTEC, Technologiepark 15, 9052 Gent, Belgium
| | - Eliseo Ferrante
- KU Leuven, Laboratory of Socioecology and Social Evolution, Naamsestraat 59, 3000 Leuven, Belgium
| | - Pieter Simoens
- Ghent University-imec, IDLab-INTEC, Technologiepark 15, 9052 Gent, Belgium
| | - Cristián Huepe
- CHuepe Labs, 814 W 19th Street 1F, Chicago, IL 60608, USA.,Northwestern Institute on Complex Systems & ESAM, Northwestern University, Evanston, IL 60208, USA
| |
Collapse
|
28
|
Hong CF, Chen YC, Chen WC, Tu KC, Tsai MH, Chan YK, Yu SS. Construction of diagnosis system and gene regulatory networks based on microarray analysis. J Biomed Inform 2018; 81:61-73. [PMID: 29550394 DOI: 10.1016/j.jbi.2018.03.008] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/02/2017] [Revised: 01/30/2018] [Accepted: 03/12/2018] [Indexed: 01/02/2023]
Abstract
A microarray analysis generally contains expression data of thousands of genes, but most of them are irrelevant to the disease of interest, making analyzing the genes concerning specific diseases complicated. Therefore, filtering out a few essential genes as well as their regulatory networks is critical, and a disease can be easily diagnosed just depending on the expression profiles of a few critical genes. In this study, a target gene screening (TGS) system, which is a microarray-based information system that integrates F-statistics, pattern recognition matching, a two-layer K-means classifier, a Parameter Detection Genetic Algorithm (PDGA), a genetic-based gene selector (GBG selector) and the association rule, was developed to screen out a small subset of genes that can discriminate malignant stages of cancers. During the first stage, F-statistic, pattern recognition matching, and a two-layer K-means classifier were applied in the system to filter out the 20 critical genes most relevant to ovarian cancer from 9600 genes, and the PDGA was used to decide the fittest values of the parameters for these critical genes. Among the 20 critical genes, 15 are associated with cancer progression. In the second stage, we further employed a GBG selector and the association rule to screen out seven target gene sets, each with only four to six genes, and each of which can precisely identify the malignancy stage of ovarian cancer based on their expression profiles. We further deduced the gene regulatory networks of the 20 critical genes by applying the Pearson correlation coefficient to evaluate the correlationship between the expression of each gene at the same stages and at different stages. Correlationships between gene pairs were calculated, and then, three regulatory networks were deduced. Their correlationships were further confirmed by the Ingenuity pathway analysis. The prognostic significances of the genes identified via regulatory networks were examined using online tools, and most represented biomarker candidates. In summary, our proposed system provides a new strategy to identify critical genes or biomarkers, as well as their regulatory networks, from microarray data.
Collapse
Affiliation(s)
- Chun-Fu Hong
- Department of Long-Term Care, National Quemoy University, Kinmen County 892, Taiwan, ROC
| | - Ying-Chen Chen
- Institute of Genomics and Bioinformatics, National Chung Hsing University, Taichung City 402, Taiwan, ROC
| | - Wei-Chun Chen
- Department of Management Information System, National Chung Hsing University, Taichung City 402, Taiwan, ROC
| | - Keng-Chang Tu
- Deparment of Computer Science and Engineering, National Chung Hsing University, Taichung City 402, Taiwan, ROC
| | - Meng-Hsiun Tsai
- Department of Management Information System, National Chung Hsing University, Taichung City 402, Taiwan, ROC.
| | - Yung-Kuan Chan
- Department of Management Information System, National Chung Hsing University, Taichung City 402, Taiwan, ROC.
| | - Shyr Shen Yu
- Deparment of Computer Science and Engineering, National Chung Hsing University, Taichung City 402, Taiwan, ROC
| |
Collapse
|
29
|
Xiong W, Wang C, Zhang X, Yang Q, Shao R, Lai J, Du C. Highly interwoven communities of a gene regulatory network unveil topologically important genes for maize seed development. THE PLANT JOURNAL : FOR CELL AND MOLECULAR BIOLOGY 2017; 92:1143-1156. [PMID: 29072883 DOI: 10.1111/tpj.13750] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/02/2017] [Revised: 10/10/2017] [Accepted: 10/17/2017] [Indexed: 06/07/2023]
Abstract
The complex interactions between transcription factors (TFs) and their target genes in a spatially and temporally specific manner are crucial to all cellular processes. Reconstruction of gene regulatory networks (GRNs) from gene expression profiles can help to decipher TF-gene regulations in a variety of contexts; however, the inevitable prediction errors of GRNs hinder optimal data mining of RNA-Seq transcriptome profiles. Here we perform an integrative study of Zea mays (maize) seed development in order to identify key genes in a complex developmental process. First, we reverse engineered a GRN from 78 maize seed transcriptome profiles. Then, we studied collective gene interaction patterns and uncovered highly interwoven network communities as the building blocks of the GRN. One community, composed of mostly unknown genes interacting with opaque2, brittle endosperm1 and shrunken2, contributes to seed phenotypes. Another community, composed mostly of genes expressed in the basal endosperm transfer layer, is responsible for nutrient transport. We further integrated our inferred GRN with gene expression patterns in different seed compartments and at various developmental stages and pathways. The integration facilitated a biological interpretation of the GRN. Our yeast one-hybrid assays verified six out of eight TF-promoter bindings in the reconstructed GRN. This study identified topologically important genes in interwoven network communities that may be crucial to maize seed development.
Collapse
Affiliation(s)
- Wenwei Xiong
- College of Agronomy, Henan Agricultural University, Zhengzhou, 450002, China
- Department of Biology, Montclair State University, Montclair, NJ, 07043, USA
| | - Chunlei Wang
- National Maize Improvement Center, China Agricultural University, Beijing, 100083, China
| | - Xiangbo Zhang
- National Maize Improvement Center, China Agricultural University, Beijing, 100083, China
| | - Qinghua Yang
- College of Agronomy, Henan Agricultural University, Zhengzhou, 450002, China
| | - Ruixin Shao
- College of Agronomy, Henan Agricultural University, Zhengzhou, 450002, China
| | - Jinsheng Lai
- National Maize Improvement Center, China Agricultural University, Beijing, 100083, China
| | - Chunguang Du
- College of Agronomy, Henan Agricultural University, Zhengzhou, 450002, China
- Department of Biology, Montclair State University, Montclair, NJ, 07043, USA
| |
Collapse
|
30
|
Cencetti G, Bagnoli F, Battistelli G, Chisci L, Fanelli D. Control of multidimensional systems on complex network. PLoS One 2017; 12:e0184431. [PMID: 28892493 PMCID: PMC5593194 DOI: 10.1371/journal.pone.0184431] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/12/2017] [Accepted: 08/20/2017] [Indexed: 11/19/2022] Open
Abstract
Multidimensional systems coupled via complex networks are widespread in nature and thus frequently invoked for a large plethora of interesting applications. From ecology to physics, individual entities in mutual interactions are grouped in families, homogeneous in kind. These latter interact selectively, through a sequence of self-consistently regulated steps, whose deeply rooted architecture is stored in the assigned matrix of connections. The asymptotic equilibrium eventually attained by the system, and its associated stability, can be assessed by employing standard nonlinear dynamics tools. For many practical applications, it is however important to externally drive the system towards a desired equilibrium, which is resilient, hence stable, to external perturbations. To this end we here consider a system made up of N interacting populations which evolve according to general rate equations, bearing attributes of universality. One species is added to the pool of interacting families and used as a dynamical controller to induce novel stable equilibria. Use can be made of the root locus method to shape the needed control, in terms of intrinsic reactivity and adopted protocol of injection. The proposed method is tested on both synthetic and real data, thus enabling to demonstrate its robustness and versatility.
Collapse
Affiliation(s)
- Giulia Cencetti
- Dipartimento di Ingegneria dell’Informazione, Università degli Studi di Firenze, Via S. Marta 3, Florence, Italy
- Dipartimento di Fisica e Astronomia and CSDC, Università degli Studi di Firenze, via G. Sansone 1, Sesto Fiorentino, Italy
- INFN Sezione di Firenze, via G. Sansone 1, Sesto Fiorentino, Italy
| | - Franco Bagnoli
- Dipartimento di Fisica e Astronomia and CSDC, Università degli Studi di Firenze, via G. Sansone 1, Sesto Fiorentino, Italy
- INFN Sezione di Firenze, via G. Sansone 1, Sesto Fiorentino, Italy
| | - Giorgio Battistelli
- Dipartimento di Ingegneria dell’Informazione, Università degli Studi di Firenze, Via S. Marta 3, Florence, Italy
| | - Luigi Chisci
- Dipartimento di Ingegneria dell’Informazione, Università degli Studi di Firenze, Via S. Marta 3, Florence, Italy
| | - Duccio Fanelli
- Dipartimento di Fisica e Astronomia and CSDC, Università degli Studi di Firenze, via G. Sansone 1, Sesto Fiorentino, Italy
- INFN Sezione di Firenze, via G. Sansone 1, Sesto Fiorentino, Italy
| |
Collapse
|
31
|
Kordmahalleh MM, Sefidmazgi MG, Harrison SH, Homaifar A. Identifying time-delayed gene regulatory networks via an evolvable hierarchical recurrent neural network. BioData Min 2017; 10:29. [PMID: 28785315 PMCID: PMC5543747 DOI: 10.1186/s13040-017-0146-4] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/07/2017] [Accepted: 07/14/2017] [Indexed: 01/25/2023] Open
Abstract
BACKGROUND The modeling of genetic interactions within a cell is crucial for a basic understanding of physiology and for applied areas such as drug design. Interactions in gene regulatory networks (GRNs) include effects of transcription factors, repressors, small metabolites, and microRNA species. In addition, the effects of regulatory interactions are not always simultaneous, but can occur after a finite time delay, or as a combined outcome of simultaneous and time delayed interactions. Powerful biotechnologies have been rapidly and successfully measuring levels of genetic expression to illuminate different states of biological systems. This has led to an ensuing challenge to improve the identification of specific regulatory mechanisms through regulatory network reconstructions. Solutions to this challenge will ultimately help to spur forward efforts based on the usage of regulatory network reconstructions in systems biology applications. METHODS We have developed a hierarchical recurrent neural network (HRNN) that identifies time-delayed gene interactions using time-course data. A customized genetic algorithm (GA) was used to optimize hierarchical connectivity of regulatory genes and a target gene. The proposed design provides a non-fully connected network with the flexibility of using recurrent connections inside the network. These features and the non-linearity of the HRNN facilitate the process of identifying temporal patterns of a GRN. RESULTS Our HRNN method was implemented with the Python language. It was first evaluated on simulated data representing linear and nonlinear time-delayed gene-gene interaction models across a range of network sizes and variances of noise. We then further demonstrated the capability of our method in reconstructing GRNs of the Saccharomyces cerevisiae synthetic network for in vivo benchmarking of reverse-engineering and modeling approaches (IRMA). We compared the performance of our method to TD-ARACNE, HCC-CLINDE, TSNI and ebdbNet across different network sizes and levels of stochastic noise. We found our HRNN method to be superior in terms of accuracy for nonlinear data sets with higher amounts of noise. CONCLUSIONS The proposed method identifies time-delayed gene-gene interactions of GRNs. The topology-based advancement of our HRNN worked as expected by more effectively modeling nonlinear data sets. As a non-fully connected network, an added benefit to HRNN was how it helped to find the few genes which regulated the target gene over different time delays.
Collapse
Affiliation(s)
- Mina Moradi Kordmahalleh
- Department of Electrical and Computer Engineering, North Carolina A&T State University, 1601 E. Market Street, Greensboro, 27411 NC USA
| | - Mohammad Gorji Sefidmazgi
- Department of Electrical and Computer Engineering, North Carolina A&T State University, 1601 E. Market Street, Greensboro, 27411 NC USA
| | - Scott H Harrison
- Department of Biology, North Carolina A&T State University, 1601 E. Market Street, Greensboro, 27411 NC USA
| | - Abdollah Homaifar
- Department of Electrical and Computer Engineering, North Carolina A&T State University, 1601 E. Market Street, Greensboro, 27411 NC USA
| |
Collapse
|
32
|
Zankoc C, Fanelli D, Ginelli F, Livi R. Intertangled stochastic motifs in networks of excitatory-inhibitory units. Phys Rev E 2017; 96:022308. [PMID: 28950520 DOI: 10.1103/physreve.96.022308] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/12/2017] [Indexed: 06/07/2023]
Abstract
A stochastic model of excitatory and inhibitory interactions which bears universality traits is introduced and studied. The endogenous component of noise, stemming from finite size corrections, drives robust internode correlations that persist at large distances. Antiphase synchrony at small frequencies is resolved on adjacent nodes and found to promote the spontaneous generation of long-ranged stochastic patterns that invade the network as a whole. These patterns are lacking under the idealized deterministic scenario, and could provide hints on how living systems implement and handle a large gallery of delicate computational tasks.
Collapse
Affiliation(s)
- Clement Zankoc
- Dipartimento di Fisica e Astronomia and CSDC, Università degli Studi di Firenze, Via G. Sansone 1, I-50019 Sesto Fiorentino, Italia
- INFN Sezione di Firenze, Via G. Sansone 1, I-50019 Sesto Fiorentino, Italia
| | - Duccio Fanelli
- Dipartimento di Fisica e Astronomia and CSDC, Università degli Studi di Firenze, Via G. Sansone 1, I-50019 Sesto Fiorentino, Italia
- INFN Sezione di Firenze, Via G. Sansone 1, I-50019 Sesto Fiorentino, Italia
| | - Francesco Ginelli
- SUPA, Institute for Complex Systems and Mathematical Biology, Kings College, University of Aberdeen, Aberdeen AB24 3UE, United Kingdom
| | - Roberto Livi
- Dipartimento di Fisica e Astronomia and CSDC, Università degli Studi di Firenze, Via G. Sansone 1, I-50019 Sesto Fiorentino, Italia
- INFN Sezione di Firenze, Via G. Sansone 1, I-50019 Sesto Fiorentino, Italia
| |
Collapse
|
33
|
Inhibitors Alter the Stochasticity of Regulatory Proteins to Force Cells to Switch to the Other State in the Bistable System. Sci Rep 2017; 7:4413. [PMID: 28667253 PMCID: PMC5493615 DOI: 10.1038/s41598-017-04596-7] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/28/2017] [Accepted: 05/17/2017] [Indexed: 12/19/2022] Open
Abstract
The cellular behaviors under the control of genetic circuits are subject to stochastic fluctuations, or noise. The stochasticity in gene regulation, far from a nuisance, has been gradually appreciated for its unusual function in cellular activities. In this work, with Chemical Master Equation (CME), we discovered that the addition of inhibitors altered the stochasticity of regulatory proteins. For a bistable system of a mutually inhibitory network, such a change of noise led to the migration of cells in the bimodal distribution. We proposed that the consumption of regulatory protein caused by the addition of inhibitor is not the only reason for pushing cells to the specific state; the change of the intracellular stochasticity is also the main cause for the redistribution. For the level of the inhibitor capable of driving 99% of cells, if there is no consumption of regulatory protein, 88% of cells were guided to the specific state. It implied that cells were pushed, by the inhibitor, to the specific state due to the change of stochasticity.
Collapse
|
34
|
Baran NM, McGrath PT, Streelman JT. Applying gene regulatory network logic to the evolution of social behavior. Proc Natl Acad Sci U S A 2017; 114:5886-5893. [PMID: 28584121 PMCID: PMC5468628 DOI: 10.1073/pnas.1610621114] [Citation(s) in RCA: 22] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/12/2023] Open
Abstract
Animal behavior is ultimately the product of gene regulatory networks (GRNs) for brain development and neural networks for brain function. The GRN approach has advanced the fields of genomics and development, and we identify organizational similarities between networks of genes that build the brain and networks of neurons that encode brain function. In this perspective, we engage the analogy between developmental networks and neural networks, exploring the advantages of using GRN logic to study behavior. Applying the GRN approach to the brain and behavior provides a quantitative and manipulative framework for discovery. We illustrate features of this framework using the example of social behavior and the neural circuitry of aggression.
Collapse
Affiliation(s)
- Nicole M Baran
- School of Biological Sciences, Georgia Institute of Technology, Atlanta, GA 30332
| | - Patrick T McGrath
- School of Biological Sciences, Georgia Institute of Technology, Atlanta, GA 30332
- The Petit Institute for Bioengineering and Bioscience, Georgia Institute of Technology, Atlanta, GA 30332
| | - J Todd Streelman
- School of Biological Sciences, Georgia Institute of Technology, Atlanta, GA 30332;
- The Petit Institute for Bioengineering and Bioscience, Georgia Institute of Technology, Atlanta, GA 30332
| |
Collapse
|
35
|
Nguyen Ba AN, Strome B, Osman S, Legere EA, Zarin T, Moses AM. Parallel reorganization of protein function in the spindle checkpoint pathway through evolutionary paths in the fitness landscape that appear neutral in laboratory experiments. PLoS Genet 2017; 13:e1006735. [PMID: 28410373 PMCID: PMC5409178 DOI: 10.1371/journal.pgen.1006735] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/07/2016] [Revised: 04/28/2017] [Accepted: 04/05/2017] [Indexed: 11/22/2022] Open
Abstract
Regulatory networks often increase in complexity during evolution through gene duplication and divergence of component proteins. Two models that explain this increase in complexity are: 1) adaptive changes after gene duplication, such as resolution of adaptive conflicts, and 2) non-adaptive processes such as duplication, degeneration and complementation. Both of these models predict complementary changes in the retained duplicates, but they can be distinguished by direct fitness measurements in organisms with short generation times. Previously, it has been observed that repeated duplication of an essential protein in the spindle checkpoint pathway has occurred multiple times over the eukaryotic tree of life, leading to convergent protein domain organization in its duplicates. Here, we replace the paralog pair in S. cerevisiae with a single-copy protein from a species that did not undergo gene duplication. Surprisingly, using quantitative fitness measurements in laboratory conditions stressful for the spindle-checkpoint pathway, we find no evidence that reorganization of protein function after gene duplication is beneficial. We then reconstruct several evolutionary intermediates from the inferred ancestral network to the extant one, and find that, at the resolution of our assay, there exist stepwise mutational paths from the single protein to the divergent pair of extant proteins with no apparent fitness defects. Parallel evolution has been taken as strong evidence for natural selection, but our results suggest that even in these cases, reorganization of protein function after gene duplication may be explained by neutral processes.
Collapse
Affiliation(s)
- Alex N. Nguyen Ba
- Department of Cell & Systems Biology, University of Toronto, Toronto, Ontario, Canada
- Center for Analysis of Evolution and Function, University of Toronto, Toronto, Ontario, Canada
| | - Bob Strome
- Department of Cell & Systems Biology, University of Toronto, Toronto, Ontario, Canada
| | - Selma Osman
- Department of Cell & Systems Biology, University of Toronto, Toronto, Ontario, Canada
| | - Elizabeth-Ann Legere
- Department of Ecology and Evolutionary Biology, University of Toronto, Toronto, Ontario, Canada
| | - Taraneh Zarin
- Department of Cell & Systems Biology, University of Toronto, Toronto, Ontario, Canada
| | - Alan M. Moses
- Department of Cell & Systems Biology, University of Toronto, Toronto, Ontario, Canada
- Center for Analysis of Evolution and Function, University of Toronto, Toronto, Ontario, Canada
- Department of Ecology and Evolutionary Biology, University of Toronto, Toronto, Ontario, Canada
| |
Collapse
|
36
|
Pazhamala LT, Purohit S, Saxena RK, Garg V, Krishnamurthy L, Verdier J, Varshney RK. Gene expression atlas of pigeonpea and its application to gain insights into genes associated with pollen fertility implicated in seed formation. JOURNAL OF EXPERIMENTAL BOTANY 2017; 68:2037-2054. [PMID: 28338822 PMCID: PMC5429002 DOI: 10.1093/jxb/erx010] [Citation(s) in RCA: 33] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/20/2023]
Abstract
Pigeonpea (Cajanus cajan) is an important grain legume of the semi-arid tropics, mainly used for its protein rich seeds. To link the genome sequence information with agronomic traits resulting from specific developmental processes, a Cajanus cajan gene expression atlas (CcGEA) was developed using the Asha genotype. Thirty tissues/organs representing developmental stages from germination to senescence were used to generate 590.84 million paired-end RNA-Seq data. The CcGEA revealed a compendium of 28 793 genes with differential, specific, spatio-temporal and constitutive expression during various stages of development in different tissues. As an example to demonstrate the application of the CcGEA, a network of 28 flower-related genes analysed for cis-regulatory elements and splicing variants has been identified. In addition, expression analysis of these candidate genes in male sterile and male fertile genotypes suggested their critical role in normal pollen development leading to seed formation. Gene network analysis also identified two regulatory genes, a pollen-specific SF3 and a sucrose-proton symporter, that could have implications for improvement of agronomic traits such as seed production and yield. In conclusion, the CcGEA provides a valuable resource for pigeonpea to identify candidate genes involved in specific developmental processes and to understand the well-orchestrated growth and developmental process in this resilient crop.
Collapse
Affiliation(s)
- Lekha T Pazhamala
- International Crops Research Institute for the Semi-Arid Tropics (ICRISAT), Hyderabad 502 324, India
| | - Shilp Purohit
- International Crops Research Institute for the Semi-Arid Tropics (ICRISAT), Hyderabad 502 324, India
| | - Rachit K Saxena
- International Crops Research Institute for the Semi-Arid Tropics (ICRISAT), Hyderabad 502 324, India
| | - Vanika Garg
- International Crops Research Institute for the Semi-Arid Tropics (ICRISAT), Hyderabad 502 324, India
| | - L Krishnamurthy
- International Crops Research Institute for the Semi-Arid Tropics (ICRISAT), Hyderabad 502 324, India
| | - Jerome Verdier
- INRA - Research Institute in Horticulture and Seeds (IRHS), 49071 Beaucouze, France
| | - Rajeev K Varshney
- International Crops Research Institute for the Semi-Arid Tropics (ICRISAT), Hyderabad 502 324, India
- School of Plant Biology and Institute of Agriculture, University of Western Australia, 35 Stirling Highway, Crawley, WA, 6009, Australia
| |
Collapse
|
37
|
Mohammadi P, Beerenwinkel N, Benenson Y. Automated Design of Synthetic Cell Classifier Circuits Using a Two-Step Optimization Strategy. Cell Syst 2017; 4:207-218.e14. [DOI: 10.1016/j.cels.2017.01.003] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/14/2016] [Revised: 10/25/2016] [Accepted: 01/06/2017] [Indexed: 10/20/2022]
|
38
|
Prabahar A, Natarajan J. MicroRNA mediated network motifs in autoimmune diseases and its crosstalk between genes, functions and pathways. J Immunol Methods 2017; 440:19-26. [DOI: 10.1016/j.jim.2016.10.002] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/16/2016] [Accepted: 10/05/2016] [Indexed: 12/27/2022]
|
39
|
Ait-Oudhia S, Ovacik MA, Mager DE. Systems pharmacology and enhanced pharmacodynamic models for understanding antibody-based drug action and toxicity. MAbs 2017; 9:15-28. [PMID: 27661132 PMCID: PMC5240652 DOI: 10.1080/19420862.2016.1238995] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/06/2016] [Revised: 09/02/2016] [Accepted: 09/14/2016] [Indexed: 10/21/2022] Open
Abstract
Pharmacokinetic (PK) and pharmacodynamic (PD) models seek to describe the temporal pattern of drug exposures and their associated pharmacological effects produced at micro- and macro-scales of organization. Antibody-based drugs have been developed for a large variety of diseases, with effects exhibited through a comprehensive range of mechanisms of action. Mechanism-based PK/PD and systems pharmacology models can play a major role in elucidating and integrating complex antibody pharmacological properties, such as nonlinear disposition and dynamical intracellular signaling pathways triggered by ligation to their cognate targets. Such complexities can be addressed through the use of robust computational modeling techniques that have proven powerful tools for pragmatic characterization of experimental data and for theoretical exploration of antibody efficacy and adverse effects. The primary objectives of such multi-scale mathematical models are to generate and test competing hypotheses and to predict clinical outcomes. In this review, relevant systems pharmacology and enhanced PD (ePD) models that are used as predictive tools for antibody-based drug action are reported. Their common conceptual features are highlighted, along with approaches used for modeling preclinical and clinically available data. Key examples illustrate how systems pharmacology and ePD models codify the interplay among complex biology, drug concentrations, and pharmacological effects. New hybrid modeling concepts that bridge cutting-edge systems pharmacology models with established PK/ePD models will be needed to anticipate antibody effects on disease in subpopulations and individual patients.
Collapse
Affiliation(s)
- Sihem Ait-Oudhia
- Center for Pharmacometrics and Systems Pharmacology, Department of Pharmaceutics, College of Pharmacy, University of Florida, Orlando, FL, USA
| | - Meric Ayse Ovacik
- Department of Pharmaceutical Sciences, University at Buffalo, State University of New York, Buffalo, NY, USA
| | - Donald E. Mager
- Department of Pharmaceutical Sciences, University at Buffalo, State University of New York, Buffalo, NY, USA
| |
Collapse
|
40
|
Ghafari M, Mashaghi A. On the role of topology in regulating transcriptional cascades. Phys Chem Chem Phys 2017; 19:25168-25179. [DOI: 10.1039/c7cp02671d] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Topology of interactions in a transcriptional cascade determines the behavior of its signal-response profile and the activation states of genes.
Collapse
Affiliation(s)
- Mahan Ghafari
- Leiden Academic Centre for Drug Research
- Faculty of Mathematics and Natural Sciences
- Leiden University
- Leiden
- The Netherlands
| | - Alireza Mashaghi
- Leiden Academic Centre for Drug Research
- Faculty of Mathematics and Natural Sciences
- Leiden University
- Leiden
- The Netherlands
| |
Collapse
|
41
|
Shu CC, Yeh CC, Jhang WS, Lo SC. Driving Cells to the Desired State in a Bimodal Distribution through Manipulation of Internal Noise with Biologically Practicable Approaches. PLoS One 2016; 11:e0167563. [PMID: 27911933 PMCID: PMC5135133 DOI: 10.1371/journal.pone.0167563] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/25/2016] [Accepted: 11/16/2016] [Indexed: 12/19/2022] Open
Abstract
The stochastic nature of gene regulatory networks described by Chemical Master Equation (CME) leads to the distribution of proteins. A deterministic bistability is usually reflected as a bimodal distribution in stochastic simulations. Within a certain range of the parameter space, a bistable system exhibits two stable steady states, one at the low end and the other at the high end. Consequently, it appears to have a bimodal distribution with one sub-population (mode) around the low end and the other around the high end. In most cases, only one mode is favorable, and guiding cells to the desired state is valuable. Traditionally, the population was redistributed simply by adjusting the concentration of the inducer or the stimulator. However, this method has limitations; for example, the addition of stimulator cannot drive cells to the desired state in a common bistable system studied in this work. In fact, it pushes cells only to the undesired state. In addition, it causes a position shift of the modes, and this shift could be as large as the value of the mode itself. Such a side effect might damage coordination, and this problem can be avoided by applying a new method presented in this work. We illustrated how to manipulate the intensity of internal noise by using biologically practicable methods and utilized it to prompt the population to the desired mode. As we kept the deterministic behavior untouched, the aforementioned drawback was overcome. Remarkably, more than 96% of cells has been driven to the desired state. This method is genetically applicable to biological systems exhibiting a bimodal distribution resulting from bistability. Moreover, the reaction network studied in this work can easily be extended and applied to many other systems.
Collapse
Affiliation(s)
- Che-Chi Shu
- Department of Chemical Engineering and Biotechnology, National Taipei University of Technology, Taipei City, Taiwan
- * E-mail:
| | - Chen-Chao Yeh
- Department of Chemical Engineering and Biotechnology, National Taipei University of Technology, Taipei City, Taiwan
| | - Wun-Sin Jhang
- Department of Chemical Engineering and Biotechnology, National Taipei University of Technology, Taipei City, Taiwan
| | - Shih-Chiang Lo
- Department of Chemical Engineering and Biotechnology, National Taipei University of Technology, Taipei City, Taiwan
| |
Collapse
|
42
|
de Anda-Jáuregui G, Velázquez-Caldelas TE, Espinal-Enríquez J, Hernández-Lemus E. Transcriptional Network Architecture of Breast Cancer Molecular Subtypes. Front Physiol 2016; 7:568. [PMID: 27920729 PMCID: PMC5118907 DOI: 10.3389/fphys.2016.00568] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/20/2016] [Accepted: 11/08/2016] [Indexed: 12/22/2022] Open
Abstract
Breast cancer heterogeneity is evident at the clinical, histological and molecular level. High throughput technologies allowed the identification of intrinsic subtypes that capture transcriptional differences among tumors. A remaining question is whether said differences are associated to a particular transcriptional program which involves different connections between the same molecules. In other words, whether particular transcriptional network architectures can be linked to specific phenotypes. In this work we infer, construct and analyze transcriptional networks from whole-genome gene expression microarrays, by using an information theory approach. We use 493 samples of primary breast cancer tissue classified in four molecular subtypes: Luminal A, Luminal B, Basal and HER2-enriched. For comparison, a network for non-tumoral mammary tissue (61 samples) is also inferred and analyzed. Transcriptional networks present particular architectures in each breast cancer subtype as well as in the non-tumor breast tissue. We find substantial differences between the non-tumor network and those networks inferred from cancer samples, in both structure and gene composition. More importantly, we find specific network architectural features associated to each breast cancer subtype. Based on breast cancer networks' centrality, we identify genes previously associated to the disease, either, generally (i.e., CNR2) or to a particular subtype (such as LCK). Similarly, we identify LUZP4, a gene barely explored in breast cancer, playing a role in transcriptional networks with subtype-specific relevance. With this approach we observe architectural differences between cancer and non-cancer at network level, as well as differences between cancer subtype networks which might be associated with breast cancer heterogeneity. The centrality measures of these networks allow us to identify genes with potential biomedical implications to breast cancer.
Collapse
Affiliation(s)
| | | | - Jesús Espinal-Enríquez
- Computational Genomics, National Institute of Genomic MedicineMexico City, Mexico
- Complejidad en Biología de Sistemas, Centro de Ciencias de la Complejidad, Universidad Nacional Autónoma de MéxicoMexico City, Mexico
| | - Enrique Hernández-Lemus
- Computational Genomics, National Institute of Genomic MedicineMexico City, Mexico
- Complejidad en Biología de Sistemas, Centro de Ciencias de la Complejidad, Universidad Nacional Autónoma de MéxicoMexico City, Mexico
| |
Collapse
|
43
|
CLIP-GENE: a web service of the condition specific context-laid integrative analysis for gene prioritization in mouse TF knockout experiments. Biol Direct 2016; 11:57. [PMID: 27776539 PMCID: PMC5078909 DOI: 10.1186/s13062-016-0158-x] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/02/2016] [Accepted: 10/10/2016] [Indexed: 02/06/2023] Open
Abstract
MOTIVATION Transcriptome data from the gene knockout experiment in mouse is widely used to investigate functions of genes and relationship to phenotypes. When a gene is knocked out, it is important to identify which genes are affected by the knockout gene. Existing methods, including differentially expressed gene (DEG) methods, can be used for the analysis. However, existing methods require cutoff values to select candidate genes, which can produce either too many false positives or false negatives. This hurdle can be addressed either by improving the accuracy of gene selection or by providing a method to rank candidate genes effectively, or both. Prioritization of candidate genes should consider the goals or context of the knockout experiment. As of now, there are no tools designed for both selecting and prioritizing genes from the mouse knockout data. Hence, the necessity of a new tool arises. RESULTS In this study, we present CLIP-GENE, a web service that selects gene markers by utilizing differentially expressed genes, mouse transcription factor (TF) network, and single nucleotide variant information. Then, protein-protein interaction network and literature information are utilized to find genes that are relevant to the phenotypic differences. One of the novel features is to allow researchers to specify their contexts or hypotheses in a set of keywords to rank genes according to the contexts that the user specify. We believe that CLIP-GENE will be useful in characterizing functions of TFs in mouse experiments. AVAILABILITY http://epigenomics.snu.ac.kr/CLIP-GENE REVIEWERS: This article was reviewed by Dr. Lee and Dr. Pongor.
Collapse
|
44
|
Peng X, Zhang S, Lei X. Multi-target trapping in constrained environments using gene regulatory network-based pattern formation. INT J ADV ROBOT SYST 2016. [DOI: 10.1177/1729881416670152] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2022] Open
Abstract
Inspired by the morphogenesis of biological organisms, gene regulatory network-based methods have been used in complex pattern formation of swarm robotic systems. In this article, obstacle information was embedded into the gene regulatory network model to make the robots trap targets with an expected pattern while avoiding obstacles in a distributed manner. Based on the modified gene regulatory network model, an implicit function method was adopted to represent the expected pattern which is easily adjusted by adding extra feature points. Considering environmental constraints (e.g. tunnels or gaps in which robots must adjust their pattern to conduct trapping task), a pattern adaptation strategy was proposed for the pattern modeler to adaptively adjust the expected pattern. Also to trap multiple targets, a splitting pattern adaptation strategy was proposed for diffusively moving targets so that the robots can trap each target separately with split sub-patterns. The proposed model and strategies were verified through a set of simulation with complex environmental constraints and non-consensus movements of targets.
Collapse
Affiliation(s)
- Xingguang Peng
- School of Marine Science and Technology, Northwestern Polytechnical University, Xi’an, China
| | - Shuai Zhang
- School of Marine Science and Technology, Northwestern Polytechnical University, Xi’an, China
| | - Xiaokang Lei
- School of Information and Control Engineering, Xi’an University of Architecture and Technology, Xi’an, China
| |
Collapse
|
45
|
Martin O, Krzywicki A, Zagorski M. Drivers of structural features in gene regulatory networks: From biophysical constraints to biological function. Phys Life Rev 2016; 17:124-58. [DOI: 10.1016/j.plrev.2016.06.002] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/21/2015] [Revised: 03/25/2016] [Accepted: 04/20/2016] [Indexed: 12/23/2022]
|
46
|
Wang Y, Qian X. Stochastic block coordinate Frank-Wolfe algorithm for large-scale biological network alignment. EURASIP JOURNAL ON BIOINFORMATICS & SYSTEMS BIOLOGY 2016; 2016:9. [PMID: 27110234 PMCID: PMC4826425 DOI: 10.1186/s13637-016-0041-1] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 08/24/2015] [Accepted: 03/17/2016] [Indexed: 11/10/2022]
Abstract
With increasingly “big” data available in biomedical research, deriving accurate and reproducible biology knowledge from such big data imposes enormous computational challenges. In this paper, motivated by recently developed stochastic block coordinate algorithms, we propose a highly scalable randomized block coordinate Frank-Wolfe algorithm for convex optimization with general compact convex constraints, which has diverse applications in analyzing biomedical data for better understanding cellular and disease mechanisms. We focus on implementing the derived stochastic block coordinate algorithm to align protein-protein interaction networks for identifying conserved functional pathways based on the IsoRank framework. Our derived stochastic block coordinate Frank-Wolfe (SBCFW) algorithm has the convergence guarantee and naturally leads to the decreased computational cost (time and space) for each iteration. Our experiments for querying conserved functional protein complexes in yeast networks confirm the effectiveness of this technique for analyzing large-scale biological networks.
Collapse
Affiliation(s)
- Yijie Wang
- Department of Electrical and Computer Engineering, Texas A&M University, College Station, USA
| | - Xiaoning Qian
- Department of Electrical and Computer Engineering, Texas A&M University, College Station, USA
| |
Collapse
|
47
|
Stable Gene Regulatory Network Modeling From Steady-State Data. Bioengineering (Basel) 2016; 3:bioengineering3020012. [PMID: 28952574 PMCID: PMC5597136 DOI: 10.3390/bioengineering3020012] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/19/2015] [Revised: 03/09/2016] [Accepted: 04/06/2016] [Indexed: 12/19/2022] Open
Abstract
Gene regulatory networks represent an abstract mapping of gene regulations in living cells. They aim to capture dependencies among molecular entities such as transcription factors, proteins and metabolites. In most applications, the regulatory network structure is unknown, and has to be reverse engineered from experimental data consisting of expression levels of the genes usually measured as messenger RNA concentrations in microarray experiments. Steady-state gene expression data are obtained from measurements of the variations in expression activity following the application of small perturbations to equilibrium states in genetic perturbation experiments. In this paper, the least absolute shrinkage and selection operator-vector autoregressive (LASSO-VAR) originally proposed for the analysis of economic time series data is adapted to include a stability constraint for the recovery of a sparse and stable regulatory network that describes data obtained from noisy perturbation experiments. The approach is applied to real experimental data obtained for the SOS pathway in Escherichia coli and the cell cycle pathway for yeast Saccharomyces cerevisiae. Significant features of this method are the ability to recover networks without inputting prior knowledge of the network topology, and the ability to be efficiently applied to large scale networks due to the convex nature of the method.
Collapse
|
48
|
Chang ETY, Lin YT, Galla T, Clayton RH, Eatock J. A Stochastic Individual-Based Model of the Progression of Atrial Fibrillation in Individuals and Populations. PLoS One 2016; 11:e0152349. [PMID: 27070920 PMCID: PMC4829251 DOI: 10.1371/journal.pone.0152349] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/03/2015] [Accepted: 03/11/2016] [Indexed: 12/19/2022] Open
Abstract
Models that represent the mechanisms that initiate and sustain atrial fibrillation (AF) in the heart are computationally expensive to simulate and therefore only capture short time scales of a few heart beats. It is therefore difficult to embed biophysical mechanisms into both policy-level disease models, which consider populations of patients over multiple decades, and guidelines that recommend treatment strategies for patients. The aim of this study is to link these modelling paradigms using a stylised population-level model that both represents AF progression over a long time-scale and retains a description of biophysical mechanisms. We develop a non-Markovian binary switching model incorporating three different aspects of AF progression: genetic disposition, disease/age related remodelling, and AF-related remodelling. This approach allows us to simulate individual AF episodes as well as the natural progression of AF in patients over a period of decades. Model parameters are derived, where possible, from the literature, and the model development has highlighted a need for quantitative data that describe the progression of AF in population of patients. The model produces time series data of AF episodes over the lifetimes of simulated patients. These are analysed to quantitatively describe progression of AF in terms of several underlying parameters. Overall, the model has potential to link mechanisms of AF to progression, and to be used as a tool to study clinical markers of AF or as training data for AF classification algorithms.
Collapse
Affiliation(s)
- Eugene T. Y. Chang
- Insigneo Institute for in-silico Medicine and Department of Computer Science, The University of Sheffield, Sheffield S1 4DP, United Kingdom
| | - Yen Ting Lin
- School of Physics and Astronomy, The University of Manchester, Manchester M13 9PL, United Kingdom
| | - Tobias Galla
- School of Physics and Astronomy, The University of Manchester, Manchester M13 9PL, United Kingdom
| | - Richard H. Clayton
- Insigneo Institute for in-silico Medicine and Department of Computer Science, The University of Sheffield, Sheffield S1 4DP, United Kingdom
| | - Julie Eatock
- Department of Computer Science, Brunel University London, Uxbridge UB8 3PH, Middlesex, United Kingdom
- * E-mail:
| |
Collapse
|
49
|
Chen Y, Li Y, Narayan R, Subramanian A, Xie X. Gene expression inference with deep learning. Bioinformatics 2016; 32:1832-9. [PMID: 26873929 DOI: 10.1093/bioinformatics/btw074] [Citation(s) in RCA: 192] [Impact Index Per Article: 24.0] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/05/2015] [Accepted: 02/03/2016] [Indexed: 01/14/2023] Open
Abstract
MOTIVATION Large-scale gene expression profiling has been widely used to characterize cellular states in response to various disease conditions, genetic perturbations, etc. Although the cost of whole-genome expression profiles has been dropping steadily, generating a compendium of expression profiling over thousands of samples is still very expensive. Recognizing that gene expressions are often highly correlated, researchers from the NIH LINCS program have developed a cost-effective strategy of profiling only ∼1000 carefully selected landmark genes and relying on computational methods to infer the expression of remaining target genes. However, the computational approach adopted by the LINCS program is currently based on linear regression (LR), limiting its accuracy since it does not capture complex nonlinear relationship between expressions of genes. RESULTS We present a deep learning method (abbreviated as D-GEX) to infer the expression of target genes from the expression of landmark genes. We used the microarray-based Gene Expression Omnibus dataset, consisting of 111K expression profiles, to train our model and compare its performance to those from other methods. In terms of mean absolute error averaged across all genes, deep learning significantly outperforms LR with 15.33% relative improvement. A gene-wise comparative analysis shows that deep learning achieves lower error than LR in 99.97% of the target genes. We also tested the performance of our learned model on an independent RNA-Seq-based GTEx dataset, which consists of 2921 expression profiles. Deep learning still outperforms LR with 6.57% relative improvement, and achieves lower error in 81.31% of the target genes. AVAILABILITY AND IMPLEMENTATION D-GEX is available at https://github.com/uci-cbcl/D-GEX CONTACT: xhx@ics.uci.edu SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics online.
Collapse
Affiliation(s)
- Yifei Chen
- Department of Computer Science, University of California, Irvine, CA 92697, USA Baidu Research-Big Data Lab, Beijing, 100085, China
| | - Yi Li
- Department of Computer Science, University of California, Irvine, CA 92697, USA
| | - Rajiv Narayan
- Broad Institute of MIT And Harvard, Cambridge, MA 02142, USA
| | | | - Xiaohui Xie
- Department of Computer Science, University of California, Irvine, CA 92697, USA Center for Complex Biological Systems, University of California, Irvine, CA 92697, USA
| |
Collapse
|
50
|
Dixon J, Lindemann A, McCoy JH. Transient amplification limits noise suppression in biochemical networks. Phys Rev E 2016; 93:012415. [PMID: 26871109 DOI: 10.1103/physreve.93.012415] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/01/2015] [Indexed: 06/05/2023]
Abstract
Cell physiology is orchestrated, on a molecular level, through complex networks of biochemical reactions. The propagation of random fluctuations through these networks can significantly impact cell behavior, raising challenging questions about how network design shapes the cell's ability to suppress or exploit these fluctuations. Here, drawing on insights from statistical physics, fluid dynamics, and systems biology, we explore how transient amplification phenomena arising from network connectivity naturally limit a biochemical system's ability to suppress small fluctuations around steady-state behaviors. We find that even a simple system consisting of two variables linked by a single interaction is capable of amplifying small fluctuations orders of magnitude beyond the levels predicted by linear stability theory. We also find that adding additional interactions can promote further amplification, even when these interactions implement classic design strategies known to suppress fluctuations. These results establish that transient amplification is an essential factor determining baseline noise levels in stable intracellular networks. Significantly, our analysis is not bound to specific systems or interaction mechanisms: we find that noise amplification is an emergent phenomenon found near steady states in any network containing sufficiently strong interactions, regardless of its form or function.
Collapse
Affiliation(s)
- John Dixon
- Department of Physics and Astronomy, Colby College, Waterville, Maine 04901, USA
| | - Anika Lindemann
- Department of Physics and Astronomy, Colby College, Waterville, Maine 04901, USA
| | - Jonathan H McCoy
- Department of Physics and Astronomy, Colby College, Waterville, Maine 04901, USA
| |
Collapse
|