Computational modelling of plasticity-led evolution

Plasticity-led evolution is a form of evolution where a change in the environment induces novel traits via phenotypic plasticity, after which the novel traits are genetically accommodated over generations under the novel environment. This mode of evolution is expected to resolve the problem of gradualism (i.e., evolution by the slow accumulation of mutations that induce phenotypic variation) implied by the Modern Evolutionary Synthesis, in the face of a large environmental change. While experimental works are essential for validating that plasticity-led evolution indeed happened, we need computational models to gain insight into its underlying mechanisms and make qualitative predictions. Such computational models should include the developmental process and gene-environment interactions in addition to genetics and natural selection. We point out that gene regulatory network models can incorporate all the above notions. In this review, we highlight results from computational modelling of gene regulatory networks that consolidate the criteria of plasticity-led evolution. Since gene regulatory networks are mathematically equivalent to artificial recurrent neural networks, we also discuss their analogies and discrepancies, which may help further understand the mechanisms underlying plasticity-led evolution.

Early effects of gene duplication on the robustness and phenotypic variability of gene regulatory networks

Background

Research on gene duplication is abundant and comes from a wide range of approaches, from high-throughput analyses and experimental evolution to bioinformatics and theoretical models. Notwithstanding, a consensus is still lacking regarding evolutionary mechanisms involved in evolution through gene duplication as well as the conditions that affect them. We argue that a better understanding of evolution through gene duplication requires considering explicitly that genes do not act in isolation. It demands studying how the perturbation that gene duplication implies percolates through the web of gene interactions. Due to evolution’s contingent nature, the paths that lead to the final fate of duplicates must depend strongly on the early stages of gene duplication, before gene copies have accumulated distinctive changes.

Methods

Here we use a widely-known model of gene regulatory networks to study how gene duplication affects network behavior in early stages. Such networks comprise sets of genes that cross-regulate. They organize gene activity creating the gene expression patterns that give cells their phenotypic properties. We focus on how duplication affects two evolutionarily relevant properties of gene regulatory networks: mitigation of the effect of new mutations and access to new phenotypic variants through mutation.

Results

Among other observations, we find that those networks that are better at maintaining the original phenotype after duplication are usually also better at buffering the effect of single interaction mutations and that duplication tends to enhance further this ability. Moreover, the effect of mutations after duplication depends on both the kind of mutation and genes involved in it. We also found that those phenotypes that had easier access through mutation before duplication had higher chances of remaining accessible through new mutations after duplication.

Conclusion

Our results support that gene duplication often mitigates the impact of new mutations and that this effect is not merely due to changes in the number of genes. The work that we put forward helps to identify conditions under which gene duplication may enhance evolvability and robustness to mutations.

Evolution enhances mutational robustness and suppresses the emergence of a new phenotype: A new computational approach for studying evolution

The aim of this paper is two-fold. First, we propose a new computational method to investigate the particularities of evolution. Second, we apply this method to a model of gene regulatory networks (GRNs) and explore the evolution of mutational robustness and bistability. Living systems have developed their functions through evolutionary processes. To understand the particularities of this process theoretically, evolutionary simulation (ES) alone is insufficient because the outcomes of ES depend on evolutionary pathways. We need a reference system for comparison. An appropriate reference system for this purpose is an ensemble of the randomly sampled genotypes. However, generating high-fitness genotypes by simple random sampling is difficult because such genotypes are rare. In this study, we used the multicanonical Monte Carlo method developed in statistical physics to construct a reference ensemble of GRNs and compared it with the outcomes of ES. We obtained the following results. First, mutational robustness was significantly higher in ES than in the reference ensemble at the same fitness level. Second, the emergence of a new phenotype, bistability, was delayed in evolution. Third, the bistable group of GRNs contains many mutationally fragile GRNs compared with those in the non-bistable group. This suggests that the delayed emergence of bistability is a consequence of the mutation-selection mechanism.

Living systems are products of evolution, and their present forms reflect their evolutionary history. Thus, to investigate the particularity of the evolutionary process by computer simulations, an appropriate reference system is needed for comparison with the outcomes of evolutionary simulations. In this study, we considered a model of gene regulatory networks (GRNs). Our idea was to construct a reference ensemble comprising randomly generated GRNs. To produce GRNs with high fitness values, which are rare, we employed a “rare event sampling” method developed in statistical physics. In particular, we focused on the evolution of mutational robustness. Living systems do not lose viability readily, even when some genes are mutated. This trait, called mutational robustness, has developed throughout evolution, along with functionality. Using the abovementioned method, we found that mutational robustness resulting from evolution exceeded that of the reference set. Therefore, mutational robustness is enhanced by evolution. We also found that the emergence of a new phenotype was significantly delayed in evolution. Our results suggest that this delay is a consequence of the fact that mutationally robust GRNs are favored by evolution.

On the effects of the modularity of gene regulatory networks on phenotypic variability and its association with robustness

Biological adaptations depend on natural selection sorting out those individuals that exhibit characters fit to their environment. Selection, in turn, depends on the phenotypic variation present in a population. Thus, evolutionary outcomes depend, to a certain extent, on the kind of variation that organisms can produce through random genetic perturbation, that is, their phenotypic variability. Moreover, the properties of developmental mechanisms that produce the organisms affect their phenotypic variability. Two of these properties are modularity and robustness. Modularity is the degree to which interactions occur mostly within groups of the system's elements and scarcely between elements in different groups. Robustness is the propensity of a system to endure perturbations while preserving its phenotype. In this paper, we used a model of gene regulatory networks (GRNs) to study the relationship between modularity and robustness in developmental processes and how modularity affects the variation that random genetic mutations produce in the expression patterns of GRNs. Our results show that modularity and robustness are correlated in multifunctional GRNs and that selection for one of these properties affects the other as well. We contend that these observations may help to understand why modularity and robustness are widespread in biological systems. Additionally, we found that modular networks tend to produce new expression patterns with subtle changes localized in the expression of a few groups of genes. This effect in the phenotypic variability of modular GRNs may bear important consequences for adaptive evolution: it may help to adjust the expression of one group of genes at a time, with few alterations on other previously evolved expression patterns.

Lanza V, Rodríguez-Beltrán J, Galán JC, San Millán A, Cantón R, Coque TM. Evolutionary Pathways and Trajectories in Antibiotic Resistance

Evolution is the hallmark of life. Descriptions of the evolution of microorganisms have provided a wealth of information, but knowledge regarding "what happened" has precluded a deeper understanding of "how" evolution has proceeded, as in the case of antimicrobial resistance. The difficulty in answering the "how" question lies in the multihierarchical dimensions of evolutionary processes, nested in complex networks, encompassing all units of selection, from genes to communities and ecosystems. At the simplest ontological level (as resistance genes), evolution proceeds by random (mutation and drift) and directional (natural selection) processes; however, sequential pathways of adaptive variation can occasionally be observed, and under fixed circumstances (particular fitness landscapes), evolution is predictable. At the highest level (such as that of plasmids, clones, species, microbiotas), the systems' degrees of freedom increase dramatically, related to the variable dispersal, fragmentation, relatedness, or coalescence of bacterial populations, depending on heterogeneous and changing niches and selective gradients in complex environments. Evolutionary trajectories of antibiotic resistance find their way in these changing landscapes subjected to random variations, becoming highly entropic and therefore unpredictable. However, experimental, phylogenetic, and ecogenetic analyses reveal preferential frequented paths (highways) where antibiotic resistance flows and propagates, allowing some understanding of evolutionary dynamics, modeling and designing interventions. Studies on antibiotic resistance have an applied aspect in improving individual health, One Health, and Global Health, as well as an academic value for understanding evolution. Most importantly, they have a heuristic significance as a model to reduce the negative influence of anthropogenic effects on the environment.

Using phenotypic plasticity to understand the structure and evolution of the genotype-phenotype map

Deciphering the genotype-phenotype map necessitates relating variation at the genetic level to variation at the phenotypic level. This endeavour is inherently limited by the availability of standing genetic variation, the rate of spontaneous mutation to novo genetic variants, and possible biases associated with induced mutagenesis. An interesting alternative is to instead rely on the environment as a source of variation. Many phenotypic traits change plastically in response to the environment, and these changes are generally underlain by changes in gene expression. Relating gene expression plasticity to the phenotypic plasticity of more integrated organismal traits thus provides useful information about which genes influence the development and expression of which traits, even in the absence of genetic variation. We here appraise the prospects and limits of such an environment-for-gene substitution for investigating the genotype-phenotype map. We review models of gene regulatory networks, and discuss the different ways in which they can incorporate the environment to mechanistically model phenotypic plasticity and its evolution. We suggest that substantial progress can be made in deciphering this genotype-environment-phenotype map, by connecting theory on gene regulatory network to empirical patterns of gene co-expression, and by more explicitly relating gene expression to the expression and development of phenotypes, both theoretically and empirically.

Genotype by environment interaction for gene expression in Drosophila melanogaster

The genetics of phenotypic responses to changing environments remains elusive. Using whole-genome quantitative gene expression as a model, here we study how the genetic architecture of regulatory variation in gene expression changed in a population of fully sequenced inbred Drosophila melanogaster strains when flies developed in different environments (25 °C and 18 °C). We find a substantial fraction of the transcriptome exhibited genotype by environment interaction, implicating environmentally plastic genetic architecture of gene expression. Genetic variance in expression increases at 18 °C relative to 25 °C for most genes that have a change in genetic variance. Although the majority of expression quantitative trait loci (eQTLs) for the gene expression traits in the two environments are shared and have similar effects, analysis of the environment-specific eQTLs reveals enrichment of binding sites for two transcription factors. Finally, although genotype by environment interaction in gene expression could potentially disrupt genetic networks, the co-expression networks are highly conserved across environments. Genes with higher network connectivity are under stronger stabilizing selection, suggesting that stabilizing selection on expression plays an important role in promoting network robustness.