1
|
Rastegar-Sedehi HR, Radhakrishnan C, Intissar Nehme S, Birman L, Velasquez Lau PM, Byrnes T. Nonequilibrium time dynamics of genetic evolution. Phys Rev E 2018; 98:022403. [PMID: 30253598 DOI: 10.1103/physreve.98.022403] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/30/2017] [Indexed: 11/06/2022]
Abstract
Biological systems are typically highly open, nonequilibrium systems that are very challenging to understand from a statistical mechanics perspective. While statistical treatments of evolutionary biological systems have a long and rich history, examination of the time-dependent nonequilibrium dynamics has been less studied. In this paper we first derive a generalized master equation in the genotype space for diploid organisms incorporating the processes of selection, mutation, recombination, and reproduction. The master equation is defined in terms of continuous time and can handle an arbitrary number of gene loci and alleles and can be defined in terms of an absolute population or probabilities. We examine and analytically solve several prototypical cases which illustrate the interplay of the various processes and discuss the timescales of their evolution. The entropy production during the evolution towards steady state is calculated and we find that it agrees with predictions from nonequilibrium statistical mechanics where it is large when the population distribution evolves towards a more viable genotype. The stability of the nonequilibrium steady state is confirmed using the Glansdorff-Prigogine criterion.
Collapse
Affiliation(s)
- Hamid-Reza Rastegar-Sedehi
- State Key Laboratory of Precision Spectroscopy, School of Physical and Material Sciences, East China Normal University, Shanghai 200062, China and Department of Physics, Jahrom University, 74137-66171, Jahrom, Iran
| | - Chandrashekar Radhakrishnan
- New York University Shanghai, 1555 Century Ave, Pudong, Shanghai 200122, China and NYU-ECNU Institute of Physics at NYU Shanghai, 3663 Zhongshan Road North, Shanghai 200062, China
| | - Samer Intissar Nehme
- New York University Shanghai, 1555 Century Ave, Pudong, Shanghai 200122, China and Department of Chemical Engineering and Biotechnology, University of Cambridge, Pembroke St, Cambridge CB2 3RA, United Kingdom
| | - Liev Birman
- New York University Shanghai, 1555 Century Ave, Pudong, Shanghai 200122, China and Department of Physics, University of Wisconsin-Madison, 1150 University Avenue, Madison, Wisconsin 53706, USA
| | | | - Tim Byrnes
- New York University Shanghai, 1555 Century Ave, Pudong, Shanghai 200122, China; State Key Laboratory of Precision Spectroscopy, School of Physical and Material Sciences, East China Normal University, Shanghai 200062, China; NYU-ECNU Institute of Physics at NYU Shanghai, 3663 Zhongshan Road North, Shanghai 200062, China; National Institute of Informatics, 2-1-2 Hitotsubashi, Chiyoda-ku, Tokyo 101-8430, Japan; and Department of Physics, New York University, New York, NY 10003, USA
| |
Collapse
|
2
|
Mandal S, Roychowdhury T, Chirom K, Bhattacharya A, Brojen Singh RK. Complex multifractal nature in Mycobacterium tuberculosis genome. Sci Rep 2017; 7:46395. [PMID: 28440326 PMCID: PMC5404331 DOI: 10.1038/srep46395] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/25/2016] [Accepted: 03/15/2017] [Indexed: 11/08/2022] Open
Abstract
The mutifractal and long range correlation (C(r)) properties of strings, such as nucleotide sequence can be a useful parameter for identification of underlying patterns and variations. In this study C(r) and multifractal singularity function f(α) have been used to study variations in the genomes of a pathogenic bacteria Mycobacterium tuberculosis. Genomic sequences of M. tuberculosis isolates displayed significant variations in C(r) and f(α) reflecting inherent differences in sequences among isolates. M. tuberculosis isolates can be categorised into different subgroups based on sensitivity to drugs, these are DS (drug sensitive isolates), MDR (multi-drug resistant isolates) and XDR (extremely drug resistant isolates). C(r) follows significantly different scaling rules in different subgroups of isolates, but all the isolates follow one parameter scaling law. The richness in complexity of each subgroup can be quantified by the measures of multifractal parameters displaying a pattern in which XDR isolates have highest value and lowest for drug sensitive isolates. Therefore C(r) and multifractal functions can be useful parameters for analysis of genomic sequences.
Collapse
Affiliation(s)
- Saurav Mandal
- School of Computational and Integrative Sciences, Jawaharlal Nehru University, New Delhi-110067, India
| | | | - Keilash Chirom
- School of Computational and Integrative Sciences, Jawaharlal Nehru University, New Delhi-110067, India
| | - Alok Bhattacharya
- School of Computational and Integrative Sciences, Jawaharlal Nehru University, New Delhi-110067, India
- School of Life Sciences, Jawaharlal Nehru University, New Delhi-110067, India
| | - R. K. Brojen Singh
- School of Computational and Integrative Sciences, Jawaharlal Nehru University, New Delhi-110067, India
| |
Collapse
|
3
|
Li W, Fontanelli O, Miramontes P. Size distribution of function-based human gene sets and the split-merge model. ROYAL SOCIETY OPEN SCIENCE 2016; 3:160275. [PMID: 27853602 PMCID: PMC5108952 DOI: 10.1098/rsos.160275] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 04/22/2016] [Accepted: 07/01/2016] [Indexed: 06/06/2023]
Abstract
The sizes of paralogues-gene families produced by ancestral duplication-are known to follow a power-law distribution. We examine the size distribution of gene sets or gene families where genes are grouped by a similar function or share a common property. The size distribution of Human Gene Nomenclature Committee (HGNC) gene sets deviate from the power-law, and can be fitted much better by a beta rank function. We propose a simple mechanism to break a power-law size distribution by a combination of splitting and merging operations. The largest gene sets are split into two to account for the subfunctional categories, and a small proportion of other gene sets are merged into larger sets as new common themes might be realized. These operations are not uncommon for a curator of gene sets. A simulation shows that iteration of these operations changes the size distribution of Ensembl paralogues and could lead to a distribution fitted by a rank beta function. We further illustrate application of beta rank function by the example of distribution of transcription factors and drug target genes among HGNC gene families.
Collapse
Affiliation(s)
- Wentian Li
- The Robert S. Boas Center for Genomics and Human Genetics, The Feinstein Institute for Medical Research, Northwell Health, Manhasset, NY, USA
| | - Oscar Fontanelli
- Departamento de Matemáticas, Facultad de Ciencias, Universidad Nacional Autónoma de México, Circuito Exterior, Ciudad Universitaria, México 04510 DF, México
| | - Pedro Miramontes
- Departamento de Matemáticas, Facultad de Ciencias, Universidad Nacional Autónoma de México, Circuito Exterior, Ciudad Universitaria, México 04510 DF, México
- Bioinformatics Group and Interdisciplinary Center for Bioinformatics, University of Leipzig, Haertelstrasse 16–18, 04107 Leipzig, Germany
| |
Collapse
|
4
|
Saakian DB, Hu CK. Mathematical Models of Quasi-Species Theory and Exact Results for the Dynamics. Curr Top Microbiol Immunol 2016; 392:121-39. [PMID: 26342705 DOI: 10.1007/82_2015_471] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/11/2023]
Abstract
We formulate the Crow-Kimura, discrete-time Eigen model, and continuous-time Eigen model. These models are interrelated and we established an exact mapping between them. We consider the evolutionary dynamics for the single-peak fitness and symmetric smooth fitness. We applied the quantum mechanical methods to find the exact dynamics of the evolution model with a single-peak fitness. For the smooth symmetric fitness landscape, we map exactly the evolution equations into Hamilton-Jacobi equation (HJE). We apply the method to the Crow-Kimura (parallel) and Eigen models. We get simple formulas to calculate the dynamics of the maximum of distribution and the variance. We review the existing mathematical tools of quasi-species theory.
Collapse
Affiliation(s)
- David B Saakian
- A.I. Alikhanyan National Science Laboratory (Yerevan Physics Institute) Foundation, 2 Alikhanian Brothers St., Yerevan, 375036, Armenia.
| | - Chin-Kun Hu
- Institute of Physics, Academia Sinica, Nankang, Taipei, 11529, Taiwan
| |
Collapse
|
5
|
Galstyan V, Saakian DB. Dynamics of the chemical master equation, a strip of chains of equations in d-dimensional space. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2012; 86:011125. [PMID: 23005386 DOI: 10.1103/physreve.86.011125] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/03/2011] [Revised: 05/09/2012] [Indexed: 06/01/2023]
Abstract
We investigate the multichain version of the chemical master equation, when there are transitions between different states inside the long chains, as well as transitions between (a few) different chains. In the discrete version, such a model can describe the connected diffusion processes with jumps between different types. We apply the Hamilton-Jacobi equation to solve some aspects of the model. We derive exact (in the limit of infinite number of particles) results for the dynamic of the maximum of the distribution and the variance of distribution.
Collapse
Affiliation(s)
- Vahe Galstyan
- Quantum College, Bagratunyats 23/2, Yerevan 0046, Armenia
| | | |
Collapse
|
6
|
Woo HJ, Vijaya Satya R, Reifman J. Thermodynamic basis for the emergence of genomes during prebiotic evolution. PLoS Comput Biol 2012; 8:e1002534. [PMID: 22693440 PMCID: PMC3364946 DOI: 10.1371/journal.pcbi.1002534] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/02/2011] [Accepted: 04/08/2012] [Indexed: 01/02/2023] Open
Abstract
The RNA world hypothesis views modern organisms as descendants of RNA molecules. The earliest RNA molecules must have been random sequences, from which the first genomes that coded for polymerase ribozymes emerged. The quasispecies theory by Eigen predicts the existence of an error threshold limiting genomic stability during such transitions, but does not address the spontaneity of changes. Following a recent theoretical approach, we applied the quasispecies theory combined with kinetic/thermodynamic descriptions of RNA replication to analyze the collective behavior of RNA replicators based on known experimental kinetics data. We find that, with increasing fidelity (relative rate of base-extension for Watson-Crick versus mismatched base pairs), replications without enzymes, with ribozymes, and with protein-based polymerases are above, near, and below a critical point, respectively. The prebiotic evolution therefore must have crossed this critical region. Over large regions of the phase diagram, fitness increases with increasing fidelity, biasing random drifts in sequence space toward ‘crystallization.’ This region encloses the experimental nonenzymatic fidelity value, favoring evolutions toward polymerase sequences with ever higher fidelity, despite error rates above the error catastrophe threshold. Our work shows that experimentally characterized kinetics and thermodynamics of RNA replication allow us to determine the physicochemical conditions required for the spontaneous crystallization of biological information. Our findings also suggest that among many potential oligomers capable of templated replication, RNAs may have evolved to form prebiotic genomes due to the value of their nonenzymatic fidelity. A leading hypothesis for the origin of life describes a prebiotic world where RNA molecules started carrying genetic information for catalyzing their own replication. This origin of biological information is akin to the crystallization of ice from water, where ‘order’ emerges from ‘disorder.’ What does the science of such phase transformations tell us about the emergence of genomes? In this paper, we show that such thermodynamic considerations of RNA synthesis, when combined with kinetics and population dynamics, lead to the conclusion that the ‘crystallization’ of genomes from its basic elements would have been spontaneous for RNAs, but not necessarily for other potential building blocks of genomes in the prebiotic soup.
Collapse
Affiliation(s)
| | | | - Jaques Reifman
- DoD Biotechnology High Performance Computing Software Applications Institute, Telemedicine and Advanced Technology Research Center, United States Army Medical Research and Materiel Command, Fort Detrick, Maryland, United States of America
- * E-mail:
| |
Collapse
|
7
|
Ramos AF, Innocentini GCP, Hornos JEM. Exact time-dependent solutions for a self-regulating gene. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2011; 83:062902. [PMID: 21797428 DOI: 10.1103/physreve.83.062902] [Citation(s) in RCA: 27] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/15/2010] [Indexed: 05/31/2023]
Abstract
The exact time-dependent solution for the stochastic equations governing the behavior of a binary self-regulating gene is presented. Using the generating function technique to rephrase the master equations in terms of partial differential equations, we show that the model is totally integrable and the analytical solutions are the celebrated confluent Heun functions. Self-regulation plays a major role in the control of gene expression, and it is remarkable that such a microscopic model is completely integrable in terms of well-known complex functions.
Collapse
Affiliation(s)
- A F Ramos
- Escola de Artes, Ciências e Humanidades, Universidade de São Paulo, São Paulo, SP, Brazil.
| | | | | |
Collapse
|
8
|
Woo HJ, Wallqvist A. Nonequilibrium phase transitions associated with DNA replication. PHYSICAL REVIEW LETTERS 2011; 106:060601. [PMID: 21405451 DOI: 10.1103/physrevlett.106.060601] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/26/2010] [Indexed: 05/30/2023]
Abstract
Thermodynamics governing the synthesis of DNA and RNA strands under a template is considered analytically and applied to the population dynamics of competing replicators. We find a nonequilibrium phase transition for high values of polymerase fidelity in a single replicator, where the two phases correspond to stationary states with higher elongation velocity and lower error rate than the other. At the critical point, the susceptibility linking velocity to thermodynamic force diverges. The overall behavior closely resembles the liquid-vapor phase transition in equilibrium. For a population of self-replicating macromolecules, Eigen's error catastrophe transition precedes this thermodynamic phase transition during starvation. For a given thermodynamic force, the fitness of replicators increases with increasing polymerase fidelity above a threshold.
Collapse
Affiliation(s)
- Hyung-June Woo
- Biotechnology High Performance Computing Software Applications Institute, Telemedicine and Advanced Technology Research Center, US Army Medical Research and Materiel Command, Fort Dietrick, Maryland 21702, USA.
| | | |
Collapse
|