Gene expression plasticity followed by genetic change during colonization in a high-elevation environment

Phenotypic plasticity facilitates organismal invasion of novel environments, and the resultant phenotypic change may later be modified by genetic change, so called 'plasticity first.' Herein, we quantify gene expression plasticity and regulatory adaptation in a wild bird (Eurasian Tree Sparrow) from its original lowland (ancestral stage), experimentally implemented hypoxia acclimation (plastic stage), and colonized highland (colonized stage). Using a group of co-expressed genes from the cardiac and flight muscles, respectively, we demonstrate that gene expression plasticity to hypoxia tolerance is more often reversed than reinforced at the colonized stage. By correlating gene expression change with muscle phenotypes, we show that colonized tree sparrows reduce maladaptive plasticity that largely associated with decreased hypoxia tolerance. Conversely, adaptive plasticity that is congruent with increased hypoxia tolerance is often reinforced in the colonized tree sparrows. Genes displaying large levels of reinforcement or reversion plasticity (i.e. 200% of original level) show greater genetic divergence between ancestral and colonized populations. Overall, our work demonstrates that gene expression plasticity at the initial stage of high-elevation colonization can be reversed or reinforced through selection-driven adaptive modification.

Environmentally robust cis-regulatory changes underlie rapid climatic adaptation

Changes in gene expression are thought to play a major role in adaptive evolution. While it is known that gene expression is highly sensitive to the environment, very few studies have determined the influence of genetic and environmental effects on adaptive gene expression differences in natural populations. Here, we utilize allele-specific expression to characterize cis and trans gene regulatory divergence in temperate and tropical house mice in two metabolic tissues under two thermal conditions. First, we show that gene expression divergence is pervasive between populations and across thermal conditions, with roughly 5 to 10% of genes exhibiting genotype-by-environment interactions. Second, we found that most expression divergence was due to cis-regulatory changes that were stable across temperatures. In contrast, patterns of expression plasticity were largely attributable to trans-effects, which showed greater sensitivity to temperature. Nonetheless, we found a small subset of temperature-dependent cis-regulatory changes, thereby identifying loci underlying expression plasticity. Finally, we performed scans for selection in wild house mice to identify genomic signatures of rapid adaptation. Genomic outliers were enriched in genes with evidence for cis-regulatory divergence. Notably, these genes were associated with phenotypes that affected body weight and metabolism, suggesting that cis-regulatory changes are a possible mechanism for adaptive body size evolution between populations. Our results show that gene expression plasticity, largely controlled in trans, may facilitate the colonization of new environments, but that evolved changes in gene expression are largely controlled in cis, illustrating the genetic and nongenetic mechanisms underlying the establishment of populations in new environments.

Genome architecture and selective signals compensatorily shape plastic response to a new environment

Transcriptional plasticity interacts with natural selection in complex ways and is crucial for the survival of species under rapid climate change. How 3D genome architecture affects transcriptional plasticity and its interaction with genetic adaptation are unclear. We transplanted estuarine oysters to a new environment and found that genes located in active chromatin regions exhibited greater transcriptional plasticity, and changes in these regions were negatively correlated with selective signals. This indicates a trade-off between 3D active regions and selective signals in shaping plastic responses to a new environment. Specifically, a mutation, lincRNA, and changes in the accessibility of a distal enhancer potentially affect its interaction with the ManⅡa gene, which regulates the muscle function and survival of oysters. Our findings reveal that 3D genome architecture compensates for the role of genetic adaptation in environmental response to new environments and provide insights into synergetic genetic and epigenetic interactions critical for fitness-related trait and survival in a model marine species.

Gene expression studies of plastic and evolutionary responses to global warming

Phenotypic plasticity can be a rapid response for coping with global warming, yet may be insufficient to protect species from extinction. Evolutionary adaptation may reinforce adaptive or oppose maladaptive plastic responses. With advances in technology whole transcriptomes can provide us with an unprecedented overview of genes and functional processes underlying the interplay between plasticity and evolution. We advocate that insects provide ideal opportunities to study plasticity in non-adapted and thermally adapted populations to infer reaction norms across the whole transcriptome ('reactionomes'). This can advance our understanding of how the interplay between plasticity and evolution shapes responses to warming. So far, a limited number of studies suggest predominantly maladaptive plastic responses to novel environments that are reduced with time, but much more research is needed to infer general patterns.

Physiological and behavioural strategies of aquatic animals living in fluctuating environments

Shallow or near-shore environments, such as ponds, estuaries and intertidal zones, are among the most physiologically challenging of all aquatic settings. Animals inhabiting these environments experience conditions that fluctuate markedly over relatively short temporal and spatial scales. Living in these habitats requires the ability to tolerate the physiological disturbances incurred by these environmental fluctuations. This tolerance is achieved through a suite of physiological and behavioural responses that allow animals to maintain homeostasis, including the ability to dynamically modulate their physiology through reversible phenotypic plasticity. However, maintaining the plasticity to adjust to some stresses in a dynamic environment may trade off with the capacity to deal with other stressors. This paper will explore studies on select fishes and invertebrates exposed to fluctuations in dissolved oxygen, salinity and pH. We assess the physiological mechanisms these species employ to achieve homeostasis, with a focus on the plasticity of their responses, and consider the resulting physiological trade-offs in function. Finally, we discuss additional factors that may influence organismal responses to fluctuating environments, such as the presence of multiple stressors, including parasites. We echo recent calls from experimental biologists to consider physiological responses to life in naturally fluctuating environments, not only because they are interesting in their own right but also because they can reveal mechanisms that may be crucial for living with increasing environmental instability as a consequence of climate change.

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.

Selection on adaptive and maladaptive gene expression plasticity during thermal adaptation to urban heat islands

Phenotypic plasticity enables a single genotype to produce multiple phenotypes in response to environmental variation. Plasticity may play a critical role in the colonization of novel environments, but its role in adaptive evolution is controversial. Here we suggest that rapid parallel regulatory adaptation of Anolis lizards to urban heat islands is due primarily to selection for reduced and/or reversed heat-induced plasticity that is maladaptive in urban thermal conditions. We identify evidence for polygenic selection across genes of the skeletal muscle transcriptome associated with heat tolerance. Forest lizards raised in common garden conditions exhibit heat-induced changes in expression of these genes that largely correlate with decreased heat tolerance, consistent with maladaptive regulatory response to high-temperature environments. In contrast, urban lizards display reduced gene expression plasticity after heat challenge in common garden and a significant increase in gene expression change that is congruent with greater heat tolerance, a putatively adaptive state in warmer urban environments. Genes displaying maladaptive heat-induced plasticity repeatedly show greater genetic divergence between urban and forest habitats than those displaying adaptive plasticity. These results highlight the role of selection against maladaptive regulatory plasticity during rapid adaptive modification of complex systems in the wild.

Axes of multivariate sexual signal divergence among incipient species: Concordance with selection, genetic variation and phenotypic plasticity

Sexual signalling traits are often observed to diverge rapidly among populations, thereby playing a potentially key early role in the evolution of reproductive isolation. While often assumed to reflect divergent sexual selection among populations, patterns of sexual trait diversification might sometimes be biased along axes of standing additive genetic variation and covariation among trait components. Additionally, theory predicts that environmentally induced phenotypic variation might facilitate rapid trait evolution, suggesting that patterns of divergence between populations should mirror phenotypic plasticity within populations. Here, we evaluate the concordance between observed axes of multivariate sexual trait divergence and predicted divergence based on (1) interpopulation variation in sexual selection, (2) additive genetic variances and (3) temperature-related phenotypic plasticity in male courtship song among geographically isolated populations of the Hawaiian swordtail cricket, Laupala cerasina, which exhibit sexual isolation due acoustic signalling traits. The major axis of multivariate divergence, d_max , accounted for 76% of variation among population male song trait means and was moderately correlated with interpopulation differences in directional sexual selection based on female preferences. However, the majority of additive genetic variance was largely oriented away from the direction of divergence, suggesting that standing genetic variation may not play a dominant role in the patterning of signal divergence. In contrast, the axis of phenotypic plasticity strongly mirrored patterns of interpopulation phenotypic divergence, which is consistent with a role for temperature-related plasticity in facilitating instead of inhibiting male song evolution and sexual isolation in these incipient species. We propose potential mechanisms by which sexual selection might interact with phenotypic plasticity to facilitate the rapid acoustic diversification observed in this species and clade.

Inferring multilayer interactome networks shaping phenotypic plasticity and evolution

Phenotypic plasticity represents a capacity by which the organism changes its phenotypes in response to environmental stimuli. Despite its pivotal role in adaptive evolution, how phenotypic plasticity is genetically controlled remains elusive. Here, we develop a unified framework for coalescing all single nucleotide polymorphisms (SNPs) from a genome-wide association study (GWAS) into a quantitative graph. This framework integrates functional genetic mapping, evolutionary game theory, and predator-prey theory to decompose the net genetic effect of each SNP into its independent and dependent components. The independent effect arises from the intrinsic capacity of a SNP, only expressed when it is in isolation, whereas the dependent effect results from the extrinsic influence of other SNPs. The dependent effect is conceptually beyond the traditional definition of epistasis by not only characterizing the strength of epistasis but also capturing the bi-causality of epistasis and the sign of the causality. We implement functional clustering and variable selection to infer multilayer, sparse, and multiplex interactome networks from any dimension of genetic data. We design and conduct two GWAS experiments using Staphylococcus aureus, aimed to test the genetic mechanisms underlying the phenotypic plasticity of this species to vancomycin exposure and Escherichia coli coexistence. We reconstruct the two most comprehensive genetic networks for abiotic and biotic phenotypic plasticity. Pathway analysis shows that SNP-SNP epistasis for phenotypic plasticity can be annotated to protein-protein interactions through coding genes. Our model can unveil the regulatory mechanisms of significant loci and excavate missing heritability from some insignificant loci. Our multilayer genetic networks provide a systems tool for dissecting environment-induced evolution.

Noncoding variation and transcriptional plasticity promote thermal adaptation in oysters by altering energy metabolism

Genetic variation and phenotypic plasticity are both important to adaptive evolution. However, how they act together on particular traits remains poorly understood. Here, we integrated phenotypic, genomic, and transcriptomic data from two allopatric but closely related congeneric oyster species, Crassostrea angulata from southern/warm environments and Crassostrea gigas from northern/cold environments, to investigate the roles of genetic divergence and plasticity in thermal adaptation. Reciprocal transplantation experiments showed that both species had higher fitness in their native habitats than in nonnative environments, indicating strong adaptive divergence. The southern species evolved higher transcriptional plasticity, and the plasticity was adaptive, suggesting that increased plasticity is important for thermal adaptation to warm climates. Genome-wide comparisons between the two species revealed that genes under selection tended to respond to environmental changes and showed higher sequence divergence in noncoding regions. All genes under selection and related to energy metabolism exhibited habitat-specific expression with genes involved in ATP production and lipid catabolism highly expressed in warm/southern habitats, and genes involved in ATP consumption and lipid synthesis were highly expressed in cold/northern habitats. The gene for acyl-CoA desaturase, a key enzyme for lipid synthesis, showed strong selective sweep in the upstream noncoding region and lower transcription in the southern species. These results were further supported by the lower free fatty acid (FFA) but higher ATP content in southern species and habitat, pointing to significance of ATP/FFA trade-off. Our findings provide evidence that noncoding variation and transcriptional plasticity play important roles in shaping energy metabolism for thermal adaptation in oysters.

Cis-regulatory evolution spotlights species differences in the adaptive potential of gene expression plasticity

Phenotypic plasticity is the variation in phenotype that a single genotype can produce in different environments and, as such, is an important component of individual fitness. However, whether the effect of new mutations, and hence evolution, depends on the direction of plasticity remains controversial. Here, we identify the cis-acting modifications that have reshaped gene expression in response to dehydration stress in three Arabidopsis species. Our study shows that the direction of effects of most cis-regulatory variants differentiating the response between A. thaliana and the sister species A. lyrata and A. halleri depends on the direction of pre-existing plasticity in gene expression. A comparison of the rate of cis-acting variant accumulation in each lineage indicates that the selective forces driving adaptive evolution in gene expression favors regulatory changes that magnify the stress response in A. lyrata. The evolutionary constraints measured on the amino-acid sequence of these genes support this interpretation. In contrast, regulatory changes that mitigate the plastic response to stress evolved more frequently in A. halleri. Our results demonstrate that pre-existing plasticity may be a stepping stone for adaptation, but its selective remodeling differs between lineages.

Predator- and competitor-induced responses in amphibian populations that evolved different levels of pesticide tolerance

Exposure to agrochemicals can drive rapid phenotypic and genetic changes in exposed populations. For instance, amphibian populations living far from agriculture (a proxy for agrochemical exposure) exhibit low pesticide tolerance, but they can be induced to possess high tolerance following a sublethal pesticide exposure. In contrast, amphibian populations close to agriculture exhibit high, constitutive tolerance to pesticides. A recent study has demonstrated that induced pesticide tolerance appears to have arisen from plastic responses to predator cues. As a result, we might expect that selection for constitutive pesticide tolerance in populations near agriculture (i.e., genetic assimilation) will lead to the evolution of constitutive responses to natural stressors. Using 15 wood frog (Rana sylvatica) populations from across an agricultural gradient, we conducted an outdoor mesocosm experiment to examine morphological (mass, body length, and tail depth) and behavioral responses (number of tadpoles observed and overall activity) of tadpoles exposed to three stressor environments (no-stressor, competitors, or predator cues). We discovered widespread differences in tadpole traits among populations and stressor environments, but no population-by-environment interaction. Subsequent linear models revealed that population distance to agriculture (DTA) was occasionally correlated with tadpole traits in a given environment and with magnitudes of plasticity, but none of the correlations were significant after Bonferroni adjustment. The magnitudes of predator and competitor plasticity were never correlated with the magnitude of pesticide-induced plasticity that we documented in a companion study. These results suggest that while predator-induced plasticity appears to have laid the foundation for the evolution of pesticide-induced plasticity and its subsequent genetic assimilation, inspection of population-level differences in plastic responses show that the evolution of pesticide-induced plasticity has not had a reciprocal effect on the evolved plastic responses to natural stressors.

Gene expression plasticity and desert adaptation in house mice

Understanding how organisms adapt to new environments is a key problem in evolution, yet it remains unclear whether phenotypic plasticity generally facilitates or hinders this process. Here we studied evolved and plastic responses to water-stress in lab-born descendants of wild house mice (Mus musculus domesticus) collected from desert and non-desert environments and measured gene expression and organismal phenotypes under control and water-stressed conditions. After many generations in the lab, desert mice consumed significantly less water than mice from other localities, indicating that this difference has a genetic basis. Under water-stress, desert mice maintained more weight than non-desert mice, and exhibited differences in blood chemistry related to osmoregulatory function. Gene expression in the kidney revealed evolved differences between mice from different environments as well as plastic responses between hydrated and dehydrated mice. Desert mice showed reduced expression plasticity under water-stress compared to non-desert mice. Importantly, non-desert mice under water-stress generally showed shifts toward desert-like expression, consistent with adaptive plasticity. Finally, we identify several co-expression modules linked to phenotypes of interest. These findings provide evidence for local adaptation after a recent invasion and suggest that adaptive plasticity may have facilitated colonization of the desert environment.

Nonparallel transcriptional divergence during parallel adaptation

How underlying mechanisms bias evolution toward predictable outcomes remains an area of active debate. In this study, we leveraged phenotypic plasticity and parallel adaptation across independent lineages of Trinidadian guppies (Poecilia reticulata) to assess the predictability of gene expression evolution during parallel adaptation. Trinidadian guppies have repeatedly and independently adapted to high- and low-predation environments in the wild. We combined this natural experiment with a laboratory breeding design to attribute transcriptional variation to the genetic influences of population of origin and developmental plasticity in response to rearing with or without predators. We observed substantial gene expression plasticity, as well as the evolution of expression plasticity itself, across populations. Genes exhibiting expression plasticity within populations were more likely to also differ in expression between populations, with the direction of population differences more likely to be opposite those of plasticity. While we found more overlap than expected by chance in genes differentially expressed between high- and low-predation populations from distinct evolutionary lineages, the majority of differentially expressed genes were not shared between lineages. Our data suggest alternative transcriptional configurations associated with shared phenotypes, highlighting a role for transcriptional flexibility in the parallel phenotypic evolution of a species known for rapid adaptation.

Developmental plasticity shapes social traits and selection in a facultatively eusocial bee

Developmental processes are an important source of phenotypic variation, but the extent to which this variation contributes to evolutionary change is unknown. We used integrative genomic analyses to explore the relationship between developmental and social plasticity in a bee species that can adopt either a social or solitary lifestyle. We find genes regulating this social flexibility also regulate development, and positive selection on these genes is influenced by their function during development. This suggests that developmental plasticity may influence the evolution of sociality. Our additional finding of genetic variants linked to differences in social behavior sheds light on how phenotypic variation derived from development may become encoded into the genome, and thus contribute to evolutionary change.

Developmental plasticity generates phenotypic variation, but how it contributes to evolutionary change is unclear. Phenotypes of individuals in caste-based (eusocial) societies are particularly sensitive to developmental processes, and the evolutionary origins of eusociality may be rooted in developmental plasticity of ancestral forms. We used an integrative genomics approach to evaluate the relationships among developmental plasticity, molecular evolution, and social behavior in a bee species (Megalopta genalis) that expresses flexible sociality, and thus provides a window into the factors that may have been important at the evolutionary origins of eusociality. We find that differences in social behavior are derived from genes that also regulate sex differentiation and metamorphosis. Positive selection on social traits is influenced by the function of these genes in development. We further identify evidence that social polyphenisms may become encoded in the genome via genetic changes in regulatory regions, specifically in transcription factor binding sites. Taken together, our results provide evidence that developmental plasticity provides the substrate for evolutionary novelty and shapes the selective landscape for molecular evolution in a major evolutionary innovation: Eusociality.

Additive and mostly adaptive plastic responses of gene expression to multiple stress in Tribolium castaneum

Gene expression is known to be highly responsive to the environment and important for adjustment of metabolism but there is also growing evidence that differences in gene regulation contribute to species divergence and differences among locally adapted populations. However, most studies so far investigated populations when divergence had already occurred. Selection acting on expression levels at the onset of adaptation to an environmental change has not been characterized. Understanding the mechanisms is further complicated by the fact that environmental change is often multivariate, meaning that organisms are exposed to multiple stressors simultaneously with potentially interactive effects. Here we use a novel approach by combining fitness and whole-transcriptome data in a large-scale experiment to investigate responses to drought, heat and their combination in Tribolium castaneum. We found that fitness was reduced by both stressors and their combined effect was almost additive. Expression data showed that stressor responses were acting independently and did not interfere physiologically. Since we measured expression and fitness within the same individuals, we were able to estimate selection on gene expression levels. We found that variation in fitness can be attributed to gene expression variation and that selection pressures were environment dependent and opposite between control and stress conditions. We could further show that plastic responses of expression were largely adaptive, i.e. in the direction that should increase fitness.

Developmental integration and evolution of labile plasticity in a complex quantitative character in a multiperiodic environment

Labile plasticity in a complex quantitative character is modeled, with multiple components contributing to net plasticity in the character. Each component has a specific development rate, norm of reaction, and cost of plasticity. For example, thermal adaptation in mammals includes seasonal fat deposition and fur growth, short-term shivering and sweating or panting, and movement between warm and cold sites. Norms of reaction do not reveal patterns of developmental integration, which must be investigated by studies of developmental dynamics in a changing environment. In a periodic environment, a labile character with a single component of plasticity is constrained by filtering environmental frequencies above the development rate and by the cost of plasticity. With multiple components of plasticity, some patterns of integration can alleviate these constraints to greatly improve fidelity of the mean phenotype tracking multiperiodic cycles in the optimum phenotype. This occurs by environmental signal amplification or inhibition through developmental integration among components and by an augmented development rate of net plasticity in the character that reduces environmental frequency filtering. When development of a component with high cost of plasticity is regulated partly by the norm of reaction of another component, evolution can diminish the reaction norm slope of the costly component without curtailing its development, thereby reducing the loss of fitness from its cost of plasticity. Apparent maladaptation in a component of plasticity may be an integral part of an adaptive pattern of developmental integration by mutual inhibition between components and compensatory evolution of a negative component reaction norm slope.

Remodeling Ancestral Phenotypic Plasticity in Local Adaptation: A New Framework to Explore the Role of Genetic Compensation in the Evolution of Homeostasis

Phenotypic plasticity is not universally adaptive. In certain cases, plasticity can result in phenotypic shifts that reduce fitness relative to the un-induced state. A common cause of such maladaptive plasticity is the co-option of ancestral developmental and physiological response systems to meet novel challenges. Because these systems evolved to meet specific challenges in an ancestral environment (e.g., localized and transient hypoxia), their co-option to meet a similar, but novel, stressor (e.g., reductions in ambient pO2 at high elevation) can lead to misdirected responses that reduce fitness. In such cases, natural selection should act to remodel phenotypic plasticity to suppress the expression of these maladaptive responses. Because these maladaptive responses reduce the fitness of colonizers in new environments, this remodeling of ancestral plasticity may be among the earliest steps in adaptive walks toward new local optima. Genetic compensation has been proposed as a general form of adaptive evolution that leads to the suppression of maladaptive plasticity to restore the ancestral trait value in the face of novel stimuli. Given their central role in the regulation of basic physiological functions, we argue that genetic compensation may often be achieved by modifications of homeostatic regulatory systems. We further suggest that genetic compensation to modify homeostatic systems can be achieved by two alternative strategies that differ in their mechanistic underpinnings; to our knowledge, these strategies have not been formally recognized by previous workers. We then consider how the mechanistic details of these alternative strategies may constrain their evolution. These considerations lead us to argue that genetic compensation is most likely to evolve by compensatory physiological changes that safeguard internal homeostatic conditions to prevent the expression of maladaptive portions of conserved reaction norms, rather than direct evolution of plasticity itself. Finally, we outline a simple experimental framework to test this hypothesis. Our goal is to stimulate research aimed at providing a deeper mechanistic understanding of whether and how phenotypic plasticity can be remodeled following environmental shifts that render ancestral responses maladaptive, an issue with increasing importance in our current era of rapid environmental change.

Maladaptive phenotypic plasticity in cardiac muscle growth is suppressed in high-altitude deer mice

How often phenotypic plasticity acts to promote or inhibit adaptive evolution is an ongoing debate among biologists. Recent work suggests that adaptive phenotypic plasticity promotes evolutionary divergence, though several studies have also suggested that maladaptive plasticity can potentiate adaptation. The role of phenotypic plasticity, adaptive, or maladaptive, in evolutionary divergence remains controversial. We examined the role of plasticity in evolutionary divergence between two species of Peromyscus mice that differ in native elevations. We used cardiac mass as a model phenotype, since ancestral hypoxia-induced responses of the heart may be both adaptive and maladaptive at high-altitude. While left ventricle growth should enhance oxygen delivery to tissues, hypertrophy of the right ventricle can lead to heart failure and death. We compared left- and right-ventricle plasticity in response to hypoxia between captive-bred P. leucopus (representing the ancestral lowland condition) and P. maniculatus from high-altitude. We found that maladaptive ancestral plasticity in right ventricle hypertrophy is reduced in high-altitude deer mice. Analysis of the heart transcriptome suggests that changes in expression of inflammatory signaling genes, particularly interferon regulatory factors, contribute to the suppression of right ventricle hypertrophy. We found weak evidence that adaptive plasticity of left ventricle mass contributes to evolution. Our results suggest that selection to suppress ancestral maladaptive plasticity plays a role in adaptation.

Modularity Facilitates Flexible Tuning of Plastic and Evolutionary Gene Expression Responses during Early Divergence

Gene expression changes have been recognized as important drivers of adaptation to changing environmental conditions. Little is known about the relative roles of plastic and evolutionary responses in complex gene expression networks during the early stages of divergence. Large gene expression data sets coupled with in silico methods for identifying coexpressed modules now enable systems genetics approaches also in nonmodel species for better understanding of gene expression responses during early divergence. Here, we combined gene coexpression analyses with population genetics to separate plastic and population (evolutionary) effects in expression networks using small salmonid populations as a model system. We show that plastic and population effects were highly variable among the six identified modules and that the plastic effects explained larger proportion of the total eigengene expression than population effects. A more detailed analysis of the population effects using a Q_ST - F_ST comparison across 16,622 annotated transcripts revealed that gene expression followed neutral expectations within modules and at the global level. Furthermore, two modules showed enrichment for genes coding for early developmental traits that have been previously identified as important phenotypic traits in thermal responses in the same model system indicating that coexpression analysis can capture expression patterns underlying ecologically important traits. We suggest that module-specific responses may facilitate the flexible tuning of expression levels to local thermal conditions. Overall, our study indicates that plasticity and neutral evolution are the main drivers of gene expression variance in the early stages of thermal adaptation in this system.

Transcriptomics and Fitness Data Reveal Adaptive Plasticity of Thermal Tolerance in Oysters Inhabiting Different Tidal Zones

Fine-scale adaptive evolution is always constrained by strong gene flow at vertical level in marine organisms. Rapid environmental fluctuations and phenotypic plasticity through optimization of fitness-related traits in organisms play important roles in shaping intraspecific divergence. The coastal systems experience strong variations in multiple abiotic environmental factors, especially the temperature. We used a typical intertidal species, Pacific oyster (Crassostrea gigas), to investigate the interaction between plasticity and adaptive evolution. We collected intertidal and subtidal oysters from two ecological niches and carried out common garden experiments for one generation. We identified fine-scale vertical adaptive divergence between intertidal and subtidal F1 progeny at both sites, based on different hierarchical phenotypes, including morphological, physiological, and molecular traits. We further quantified the global plasticity to thermal stress through transcriptomic analysis. The intertidal oysters exhibited slow growth rate. However, they showed high survival and metabolic rates under heat stress, indicating vertically fine-scale phenotypic adaptive mechanisms and evolutionary trade-offs between growth and thermal tolerance. Transcriptomic analysis confirmed that the intertidal oysters have evolved high plasticity. The genes were classified into three types: evolutionarily divergent, concordantly plastic, and adaptive plastic genes. The evolved divergence between intertidal and subtidal oysters for these gene sets showed a significant positive correlation with plastic changes of subtidal populations in response to high temperature. Furthermore, the intertidal oysters exhibited delayed large-scale increase in expressional plasticity than that in subtidal counterparts. The same direction between plasticity and selection suggests that the oysters have evolved adaptive plasticity. This implies that adaptive plasticity facilitates the oyster to adapt to severe intertidal zones. The oysters exposed to strong environmental variability are thermal tolerant and have high adaptive potential to face the current global warming. Our findings will not only provide new insights into the significant role of plasticity in adaptive evolution that can be extended to other marine invertebrates, but also provide basic information for oyster resources conservation and reef reestablishment.

Integrating Phylogenetic and Network Approaches to Study Gene Family Evolution: The Case of the AGAMOUS Family of Floral Genes

The study of gene family evolution has benefited from the use of phylogenetic tools, which can greatly inform studies of both relationships within gene families and functional divergence. Here, we propose the use of a network-based approach that in combination with phylogenetic methods can provide additional support for models of gene family evolution. We dissect the contributions of each method to the improved understanding of relationships and functions within the well-characterized family of AGAMOUS floral development genes. The results obtained with the two methods largely agreed with one another. In particular, we show how network approaches can provide improved interpretations of branches with low support in a conventional gene tree. The network approach used here may also better reflect known and suspected patterns of functional divergence relative to phylogenetic methods. Overall, we believe that the combined use of phylogenetic and network tools provide a more robust assessment of gene family evolution.