Guided evolution of in silico microbial populations in complex environments accelerates evolutionary rates through a step-wise adaptation

Abstract Background During their lifetime, microbes are exposed to environmental variations, each with its distinct spatio-temporal dynamics. Microbial communities display a remarkable degree of phenotypic plasticity, and highly-fit individuals emerge quite rapidly during microbial adaptation to novel environments. However, there exists a high variability when it comes to adaptation potential, and while adaptation occurs rapidly in certain environmental transitions, in others organisms struggle to adapt. Here, we investigate the hypothesis that the rate of evolution can both increase or decrease, depending on the similarity and complexity of the intermediate and final environments. Elucidating such dependencies paves the way towards controlling the rate and direction of evolution, which is of interest to industrial and medical applications. Results Our results show that the rate of evolution can be accelerated by evolving cell populations in sequential combinations of environments that are increasingly more complex. To quantify environmental complexity, we evaluate various information-theoretic metrics, and we provide evidence that multivariate mutual information between environmental signals in a given environment correlates well with the rate of evolution in that environment, as measured in our simulations. We find that strong positive and negative correlations between the intermediate and final environments lead to the increase of evolutionary rates, when the environmental complexity increases. Horizontal Gene Transfer is shown to further augment this acceleration, under certain conditions. Interestingly, our simulations show that weak environmental correlations lead to deceleration of evolution, regardless of environmental complexity. Further analysis of network evolution provides a mechanistic explanation of this phenomenon, as exposing cells to intermediate environments can trap the population to local neighborhoods of sub-optimal fitness

Horizontal gene transfer dynamics and distribution of fitness effects during microbial in silico evolution

<p>Abstract</p> <p>Background</p> <p>Horizontal gene transfer (HGT) is a process that facilitates the transfer of genetic material between organisms that are not directly related, and thus can affect both the rate of evolution and emergence of traits. Recent phylogenetic studies reveal HGT events are likely ubiquitous in the Tree of Life. However, our knowledge of HGT's role in evolution and biological organization is very limited, mainly due to the lack of ancestral evolutionary signatures and the difficulty to observe complex evolutionary dynamics in a laboratory setting. Here, we utilize a multi-scale microbial evolution model to comprehensively study the effect of HGT on the evolution of complex traits and organization of gene regulatory networks.</p> <p>Results</p> <p>Large-scale simulations reveal a distinct signature of the Distribution of Fitness Effect (DFE) for HGT events: during evolution, while mutation fitness effects become more negative and neutral, HGT events result in a balanced effect distribution. In either case, lethal events are significantly decreased during evolution (33.0% to 3.2%), a clear indication of mutational robustness. Interestingly, evolution was accelerated when populations were exposed to correlated environments of increasing complexity, especially in the presence of HGT, a phenomenon that warrants further investigation. High HGT rates were found to be disruptive, while the average transferred fragment size was linked to functional module size in the underlying biological network. Network analysis reveals that HGT results in larger regulatory networks, but with the same sparsity level as those evolved in its absence. Observed phenotypic variability and co-existing solutions were traced to individual gain/loss of function events, while subsequent re-wiring after fragment integration was necessary for complex traits to emerge.</p

PHOTOELECTRON SPECTRA FROM WAVE PACKET DYNAMICS

Author Institution: Department of Chemistry, University of Southern California, Los Angeles, CA 90007Franck-Condon factors for photoelectron spectra can be computed from: (1) the overlap between initial and target vibrational wave functions; and (2) Fourier transform of a wave packet time autocorrelation function. These techniques were implemented in the new spectra modeling software using harmonic well approximation and full quantum mechanical treatment. The photoelectron spectrum of N${}_3$ was modeled using \emph{ab-initio} potential energy surfaces of the cation electronic states. Anharmonic effects are discussed

A flood-based information flow analysis and network minimization method for gene regulatory networks

Abstract Background Biological networks tend to have high interconnectivity, complex topologies and multiple types of interactions. This renders difficult the identification of sub-networks that are involved in condition- specific responses. In addition, we generally lack scalable methods that can reveal the information flow in gene regulatory and biochemical pathways. Doing so will help us to identify key participants and paths under specific environmental and cellular context. Results This paper introduces the theory of network flooding, which aims to address the problem of network minimization and regulatory information flow in gene regulatory networks. Given a regulatory biological network, a set of source (input) nodes and optionally a set of sink (output) nodes, our task is to find (a) the minimal sub-network that encodes the regulatory program involving all input and output nodes and (b) the information flow from the source to the sink nodes of the network. Here, we describe a novel, scalable, network traversal algorithm and we assess its potential to achieve significant network size reduction in both synthetic and E. coli networks. Scalability and sensitivity analysis show that the proposed method scales well with the size of the network, and is robust to noise and missing data. Conclusions The method of network flooding proves to be a useful, practical approach towards information flow analysis in gene regulatory networks. Further extension of the proposed theory has the potential to lead in a unifying framework for the simultaneous network minimization and information flow analysis across various “omics” levels

A flood-based information flow analysis and network minimization method for gene regulatory networks

Abstract Background Biological networks tend to have high interconnectivity, complex topologies and multiple types of interactions. This renders difficult the identification of sub-networks that are involved in condition- specific responses. In addition, we generally lack scalable methods that can reveal the information flow in gene regulatory and biochemical pathways. Doing so will help us to identify key participants and paths under specific environmental and cellular context. Results This paper introduces the theory of network flooding, which aims to address the problem of network minimization and regulatory information flow in gene regulatory networks. Given a regulatory biological network, a set of source (input) nodes and optionally a set of sink (output) nodes, our task is to find (a) the minimal sub-network that encodes the regulatory program involving all input and output nodes and (b) the information flow from the source to the sink nodes of the network. Here, we describe a novel, scalable, network traversal algorithm and we assess its potential to achieve significant network size reduction in both synthetic and E. coli networks. Scalability and sensitivity analysis show that the proposed method scales well with the size of the network, and is robust to noise and missing data. Conclusions The method of network flooding proves to be a useful, practical approach towards information flow analysis in gene regulatory networks. Further extension of the proposed theory has the potential to lead in a unifying framework for the simultaneous network minimization and information flow analysis across various “omics” levels

A flood-based information flow analysis and network minimization method for gene regulatory networks

BACKGROUND: Biological networks tend to have high interconnectivity, complex topologies and multiple types of interactions. This renders difficult the identification of sub-networks that are involved in condition- specific responses. In addition, we generally lack scalable methods that can reveal the information flow in gene regulatory and biochemical pathways. Doing so will help us to identify key participants and paths under specific environmental and cellular context. RESULTS: This paper introduces the theory of network flooding, which aims to address the problem of network minimization and regulatory information flow in gene regulatory networks. Given a regulatory biological network, a set of source (input) nodes and optionally a set of sink (output) nodes, our task is to find (a) the minimal sub-network that encodes the regulatory program involving all input and output nodes and (b) the information flow from the source to the sink nodes of the network. Here, we describe a novel, scalable, network traversal algorithm and we assess its potential to achieve significant network size reduction in both synthetic and E. coli networks. Scalability and sensitivity analysis show that the proposed method scales well with the size of the network, and is robust to noise and missing data. CONCLUSIONS: The method of network flooding proves to be a useful, practical approach towards information flow analysis in gene regulatory networks. Further extension of the proposed theory has the potential to lead in a unifying framework for the simultaneous network minimization and information flow analysis across various “omics” levels

Evolutionary potential, cross-stress behavior and the genetic basis of acquired stress resistance in Escherichia coli.

Bacterial populations have a remarkable capacity to cope with extreme environmental fluctuations in their natural environments. In certain cases, adaptation to one stressful environment provides a fitness advantage when cells are exposed to a second stressor, a phenomenon that has been coined as cross-stress protection. A tantalizing question in bacterial physiology is how the cross-stress behavior emerges during evolutionary adaptation and what the genetic basis of acquired stress resistance is. To address these questions, we evolved Escherichia coli cells over 500 generations in five environments that include four abiotic stressors. Through growth profiling and competition assays, we identified several cases of positive and negative cross-stress behavior that span all strain-stress combinations. Resequencing the genomes of the evolved strains resulted in the identification of several mutations and gene amplifications, whose fitness effect was further assessed by mutation reversal and competition assays. Transcriptional profiling of all strains under a specific stress, NaCl-induced osmotic stress, and integration with resequencing data further elucidated the regulatory responses and genes that are involved in this phenomenon. Our results suggest that cross-stress dependencies are ubiquitous, highly interconnected, and can emerge within short timeframes. The high adaptive potential that we observed argues that bacterial populations occupy a genotypic space that enables a high phenotypic plasticity during adaptation in fluctuating environments

Conical for Stepwise, Glancing for Concerted: The Role of the Excited- State Topology in the Three-Body Dissociation of sym-Triazine Conical for Stepwise, Glancing for Concerted: The Role of the Excited-State Topology in the Three-Body Dissociation of sym

The highly debated three-body dissociation of sym-triazine to three HCN products has been investigated by translational spectroscopy and high-level ab initio calculations. Dissociation was induced by charge exchange between the sym-triazine radical cation and cesium. Calculated state energies and electronic couplings suggest that sym-triazine is produced in the 3s Rydberg and π* r n manifolds. Analysis of the topology of these manifolds along with momentum correlation in the dissociation products suggest that the 3s Rydberg manifold characterized by a conical intersection of two potential energy surfaces leads to stepwise dissociation, while the π* r n manifold consisting of a four-fold glancing intersection leads to a symmetric concerted reaction