Phenotypic plasticity and robustness: evolutionary stability theory, gene expression dynamics model, and laboratory experiments

Flux sampling is a powerful tool to study metabolism under changing environmental conditions

The development of high-throughput 'omic techniques has sparked a rising interest in genome-scale metabolic models, with applications ranging from disease diagnostics to crop adaptation. Efficient and accurate methods are required to analyze large metabolic networks. Flux sampling can be used to explore the feasible flux solutions in metabolic networks by generating probability distributions of steady-state reaction fluxes. Unlike other methods, flux sampling can be used without assuming a particular cellular objective. We have undertaken a rigorous comparison of several sampling algorithms and concluded that the coordinate hit-and-run with rounding (CHRR) algorithm is the most efficient based on both run-time and multiple convergence diagnostics. We demonstrate the power of CHRR by using it to study the metabolic changes that underlie photosynthetic acclimation to cold of Arabidopsis thaliana plant leaves. In combination with experimental measurements, we show how the regulated interplay between diurnal starch and organic acid accumulation defines the plant acclimation process. We confirm fumarate accumulation as a requirement for cold acclimation and further predict γ-aminobutyric acid to have a key role in metabolic signaling under cold conditions. These results demonstrate how flux sampling can be used to analyze the feasible flux solutions across changing environmental conditions, whereas eliminating the need to make assumptions which introduce observer bias.

Integrating molecular networks with genetic variant interpretation for precision medicine

More reliable and cheaper sequencing technologies have revealed the vast mutational landscapes characteristic of many phenotypes. The analysis of such genetic variants has led to successful identification of altered proteins underlying many Mendelian disorders. Nevertheless the simple one‐variant one‐phenotype model valid for many monogenic diseases does not capture the complexity of polygenic traits and disorders. Although experimental and computational approaches have improved detection of functionally deleterious variants and important interactions between gene products, the development of comprehensive models relating genotype and phenotypes remains a challenge in the field of genomic medicine. In this context, a new view of the pathologic state as significant perturbation of the network of interactions between biomolecules is crucial for the identification of biochemical pathways associated with complex phenotypes. Seminal studies in systems biology combined the analysis of genetic variation with protein–protein interaction networks to demonstrate that even as biological systems evolve to be robust to genetic variation, their topologies create disease vulnerabilities. More recent analyses model the impact of genetic variants as changes to the “wiring” of the interactome to better capture heterogeneity in genotype–phenotype relationships. These studies lay the foundation for using networks to predict variant effects at scale using machine‐learning or algorithmic approaches. A wealth of databases and resources for the annotation of genotype–phenotype relationships have been developed to support developments in this area. This overview describes how study of the molecular interactome has generated insights linking the organization of biological systems to disease mechanism, and how this information can enable precision medicine.

This article is categorized under:

Translational, Genomic, and Systems Medicine > Translational Medicine

Biological Mechanisms > Cell Signaling

Models of Systems Properties and Processes > Mechanistic Models

Analytical and Computational Methods > Computational Methods

Chaperones, Canalization, and Evolution of Animal Forms

Over half a century ago, British developmental biologist Conrad Hal Waddington proposed the idea of canalization, that is, homeostasis in development. Since the breakthrough that was made by Rutherford and Lindquist (1998), who proposed a role of Hsp90 in developmental buffering, chaperones have gained much attention in the study of canalization. However, recent studies have revealed that a number of other molecules are also potentially involved in canalization. Here, I introduce the emerging role of DnaJ chaperones in canalization. I also discuss how the expression levels of such buffering molecules can be altered, thereby altering organismal development. Since developmental robustness is maternally inherited in various organisms, I propose that dynamic bet hedging, an increase in within-clutch variation in offspring phenotypes that is caused by unpredictable environmental challenges to the mothers, plays a key role in altering the expression levels of buffering molecules. Investigating dynamic bet hedging at the molecular level and how it impacts upon morphological phenotypes will help our understanding of the molecular mechanisms of canalization and evolutionary processes.

Macroscopic Theory for Evolving Biological Systems Akin to Thermodynamics

We present a macroscopic theory to characterize the plasticity, robustness, and evolvability of biological responses and their fluctuations. First, linear approximation in intracellular reaction dynamics is used to demonstrate proportional changes in the expression of all cellular components in response to a given environmental stress, with the proportion coefficient determined by the change in growth rate as a consequence of the steady growth of cells. We further demonstrate that this relationship is supported through adaptation experiments of bacteria, perhaps too well as this proportionality is held even across cultures of different types of conditions. On the basis of simulations of cell models, we further show that this global proportionality is a consequence of evolution in which expression changes in response to environmental or genetic perturbations are constrained along a unique one-dimensional curve, which is a result of evolutionary robustness. It then follows that the expression changes induced by environmental changes are proportionally reduced across different components of a cell by evolution, which is akin to the Le Chatelier thermodynamics principle. Finally, with the aid of a fluctuation-response relationship, this proportionality is shown to hold between fluctuations caused by genetic changes and those caused by noise. Overall, these results and support from the theoretical and experimental literature suggest a formulation of cellular systems akin to thermodynamics, in which a macroscopic potential is given by the growth rate (or fitness) represented as a function of environmental and evolutionary changes.

Global relationships in fluctuation and response in adaptive evolution

Cells change their internal state to adapt to environmental changes, and evolve in response to the new conditions. The phenotype changes first via adaptation in response to environmental changes, and then through mutational changes in the genomic sequence, followed by selection in evolution. Here, we analysed simulated adaptive evolution using a simple cell model consisting of thousands of intracellular components, and found that the changes in their concentrations by adaptation are proportional to those by evolution across all the components, where the proportion coefficient between the two agreed well with the change in the growth rate of a cell. Furthermore, we demonstrate that the phenotypic variance in concentrations of cellular components due to (non-genetic) noise and to genomic alternations is proportional across all components. This implies that the specific phenotypes that are highly evolvable were already given by non-genetic fluctuations. These global relationships in cellular states were also supported by phenomenological theory based on steady reproduction and transcriptome analysis of laboratory evolution in Escherichia coli. These findings demonstrate that a possible evolutionary change in phenotypic state is highly restricted. Our results provide a basis for the development of a quantitative theory of plasticity and robustness in phenotypic evolution.

Challenges for complex microbial ecosystems: combination of experimental approaches with mathematical modeling

In the past couple of decades, molecular ecological techniques have been developed to elucidate microbial diversity and distribution in microbial ecosystems. Currently, modern techniques, represented by meta-omics and single cell observations, are revealing the incredible complexity of microbial ecosystems and the large degree of phenotypic variation. These studies propound that microbiological techniques are insufficient to untangle the complex microbial network. This minireview introduces the application of advanced mathematical approaches in combination with microbiological experiments to microbial ecological studies. These combinational approaches have successfully elucidated novel microbial behaviors that had not been recognized previously. Furthermore, the theoretical perspective also provides an understanding of the plasticity, robustness and stability of complex microbial ecosystems in nature.