How do phenotypic plasticity and adaptive evolution interact in a novel or changing environment? Does evolution by natural selection generally reinforce initially plastic phenotypic responses, or does it instead oppose them? And to what extent does evolution of a trait involve evolution of its plasticity? These questions have lied at the heart of research on phenotypic evolution in heterogeneous environments ever since it was realized that the environment is likely to affect the expression of many (perhaps most) characters of an individual. Importantly, this broad definition of phenotypic plasticity as change in the average phenotype of a given genotype in response to its environment of development (or expression) does not involve any statement about the adaptiveness of the plastic changes. Theory on the evolution of plasticity has devoted much effort to understanding how reaction norm should evolve under different regimes of environmental change in space and time, and depending on genetic constraints on reaction norm shapes. However on an empirically ground, the questions above have mostly been addressed for individual traits, often chosen a priori for their likeliness to exhibit adaptive plasticity, and we still lack more systematic answers. These can be provided by so-called ‘phenomic’ approaches, where a large number of traits are tracked without prior information on their biological or ecological function. A problem is that the number of phenotypic characters that can be measured in an organism is virtually infinite (and to some extent arbitrary), and that scaling issues makes it difficult to compare different sets of traits. Gene-expression levels offer a partial solution to this dilemma, as they can be considered as a very large number of traits (one per typed gene) that can be measured easily and uniformly (fold change in the number of reads in RNAseq). As for any traits, expression levels of different genes may be genetically correlated, to an extent that depends on their regulation mechanism: cis-regulatory sequences that only affect expression of neighboring genes are likely to cause independent gene expression, while more systematic modifiers of expression (e.g. trans-regulators such as transcription factors) may cause correlated genetic responses of the expression of many genes. Huang and Agrawal [1] have studied plasticity and evolution of gene expression level in young larvae of populations of Drosophila melanogaster that have evolved for about 130 generations under either a constant environment (salt or cadmium), or an environment that is heterogeneous in time or space (combining salt and cadmium). They report a wealth of results, of which we summarize the most striking here. First, among genes that (i) were initially highly plastic and (ii) evolved significant divergence in expression levels between constant environment treatments, the evolved divergence is predominantly in the opposite direction to the initial plastic response. This suggests that either plasticity was initially maladaptive, or the selective pressure changed during the evolutionary process (see below). This somewhat unexpected result strikingly mirrors that from a study published last year in Nature [2], where the same pattern was found for responses of guppies to the presence of predators. However, Huang and Agrawal [1] went beyond this study by deciphering the underlying mechanisms in several interesting ways. First, they showed that change in gene expression often occurred at genes close to SNPs with differentiated frequencies across treatments (but not at genes with differentiated SNPs in their coding sequences), suggesting that cis-regulatory sequences are involved. This is also suggested by the fact that changes in gene expression are mostly caused by the increased expression of only one allele at polymorphic loci, and is a first step towards investigating the genetic underpinnings of (co)variation in gene expression levels. Another interesting set of findings concerns evolution of plasticity in treatments with variable environments. To compare the gene-expression plasticity that evolved in these treatments to an expectation, the authors considered that the expression levels in populations maintained for a long time under constant salt or cadmium had reached an optimum. The differences between these expression levels were thus assumed to predict the level of plasticity that should evolve in a heterogeneous environment (with both cadmium and salt) under perfect environmental predictability. The authors showed that plasticity did evolve more in the expected direction in heterogeneous than in constant environments, resulting in better adapted final expression levels across environments. Taken collectively, these results provide an unprecedented set of patterns that are greatly informative on how plasticity and evolution interact in constant versus changing environments. But of course, interpretations in terms of adaptive versus maladaptive plasticity are more challenging, as the authors themselves admit. Even though environmentally determined gene expression is the basic mechanism underlying the phenotypic plasticity of most traits, it is extremely difficult to relate to more integrated phenotypes for which we can understand the selection pressures, especially in multicellular organisms. The authors have recently investigated evolutionary change of quantitative traits in these selected lines, so it might be possible to establish links between reaction norms for macroscopic traits to those for gene expression levels. Such an approach would also involve tracking gene expression throughout life, rather than only in young larvae as done here, thus putting phenotypic complexity back in the picture also for expression levels. Another difficulty is that a plastic response that was originally adaptive may be replaced by an opposite evolutionary response in the long run, without having to invoke initially maladaptive plasticity. For instance, the authors mention the possibility that a generic stress response is initially triggered by cadmium, but is eventually unnecessary and costly after evolution of genetic mechanisms for cadmium detoxification (a case of so-called genetic accommodation). In any case, this study by Huang and Agrawal [1], together with the one by Ghalambor et al. last year [2], reports novel and unexpected results, which are likely to stimulate researchers interested in plasticity and evolution in heterogeneous environments for the years to come.

References

[1] Huang Y, Agrawal AF. 2016. Experimental Evolution of Gene Expression and Plasticity in Alternative Selective Regimes. PLoS Genetics 12:e1006336. doi: 10.1371/journal.pgen.1006336

[2] Ghalambor CK, Hoke KL, Ruell EW, Fischer EK, Reznick DN, Hughes KA. 2015. Non-adaptive plasticity potentiates rapid adaptive evolution of gene expression in nature. Nature 525: 372-375. doi: 10.1038/nature15256

Studying the effects of speciation, hybridisation, and evolutionary outcomes following reproduction from divergent populations is a major research area in evolutionary genetics [1]. There are two phenomena that have been the focus of contemporary research. First, a classic concept is the formation of ‘Bateson-Dobzhansky-Muller’ incompatibilities (BDMi) [2–4] that negatively affect hybrid fitness. Here, two diverging populations accumulate mutations over time that are unique to that subpopulation. If they subsequently meet, then these mutations might negatively interact, leading to a loss in fitness or even a complete lack of reproduction. BDMi formation can be complex, involving multiple genes and the fitness changes can depend on the direction of introgression [5]. Second, such secondary contact can instead lead to heterosis, where offspring are fitter than their parental progenitors [6].

Understanding which outcomes are likely to arise require one to know the potential fitness effects of mutations underlying reproductive isolation, to determine whether they are likely to reduce or enhance fitness when hybrids are formed. This is far from an easy task, as it requires one to track mutations at several loci, along with their effects, across a fitness landscape.

The work of De Sanctis et al. [7] neatly fills in this knowledge gap, by creating a general mathematical framework for describing the consequences of a cross from two divergent populations. The derivations are based on Fisher’s Geometric Model, which is widely used to quantify selection acting on a general fitness landscape that is affected by several biological traits [8,9], and has previously been used in theoretical studies of hybridisation [10–12]. By doing so, they are able to decompose how divergence at multiple loci affects offspring fitness through both additive and dominance effects.

A key result arising from their analyses is demonstrating how offspring fitness can be captured by two main functions. The first one is the ‘net effect of evolutionary change’ that, broadly defined, measures how phenotypically divergent two populations are. The second is the ‘total amount of evolutionary change’, which reflects how many mutations contribute to divergence and the effect sizes captured by each of them. The authors illustrate these measurements using simulations covering different scenarios, demonstrating how different parental states can lead to similar fitness outcomes. They also propose experimental methods to measure the underlying mutational effects.

This study neatly demonstrates how complex genetic phenomena underlying hybridisation can be captured using fairly simple mathematical formulae. This powerful approach will thus open the door for future research to investigate hybridisation in more detail, whether it is by expanding on these theoretical models or using the elegant outcomes to quantify fitness effects in experiments.

References

1. Coyne JA, Orr HA. Speciation. Sunderland, Mass: Sinauer Associates; 2004.
2. Bateson W, Seward A. Darwin and modern science. Heredity and variation in modern lights. 1909;85: 101. https://doi.org/10.1017/CBO9780511693953.007
3. Dobzhansky T. Genetics and the Origin of Species. Columbia university press; 1937.
4. Muller HJ. Isolating mechanisms, evolution and temperature. Biol Symp. 1942;6: 71-125.
5. Fraïsse C, Elderfield JAD, Welch JJ. The genetics of speciation: are complex incompatibilities easier to evolve? J Evol Biol. 2014;27: 688-699. https://doi.org/10.1111/jeb.12339
6. Birchler JA, Yao H, Chudalayandi S, Vaiman D, Veitia RA. Heterosis. The Plant Cell. 2010;22: 2105-2112. https://doi.org/10.1105/tpc.110.076133
7. De Sanctis B, Schneemann H, Welch JJ. How does the mode of evolutionary divergence affect reproductive isolation? bioRxiv. 2022. 2022.03.08.483443 version 4. https://doi.org/10.1101/2022.03.08.483443 
8. Fisher RA. The genetical theory of natural selection. Oxford: The Clarendon Press; 1930. https://doi.org/10.5962/bhl.title.27468 
9. Tenaillon O. The Utility of Fisher's Geometric Model in Evolutionary Genetics. Annu Rev Ecol Evol Syst. 2014;45: 179-201. https://doi.org/10.1146/annurev-ecolsys-120213-091846
10. Barton NH. The role of hybridization in evolution. Molecular Ecology. 2001;10: 551-568. https://doi.org/10.1046/j.1365-294x.2001.01216.x 
11. Chevin L-M, Decorzent G, Lenormand T. Niche Dimensionality and The Genetics of Ecological Speciation. Evolution. 2014;68: 1244-1256. https://doi.org/10.1111/evo.12346 
12. Fraïsse C, Gunnarsson PA, Roze D, Bierne N, Welch JJ. The genetics of speciation: Insights from Fisher's geometric model. Evolution. 2016;70: 1450-1464. https://doi.org/10.1111/evo.12968