Natural selection canalizes expression variation of environmentally induced plasticity-enabling genes

Adaptive divergence generates distinct plastic responses in two closely related Senecio species

The evolution of plastic responses to external cues allows species to maintain fitness in response to the environmental variations they regularly experience. However, it remains unclear how plasticity evolves during adaptation. To test whether distinct patterns of plasticity are associated with adaptive divergence, we quantified plasticity for two closely related but ecologically divergent Sicilian daisy species (Senecio, Asteraceae). We sampled 40 representative genotypes of each species from their native range on Mt. Etna and then reciprocally transplanted multiple clones of each genotype into four field sites along an elevational gradient that included the native elevational range of each species, and two intermediate elevations. At each elevation, we quantified survival and measured leaf traits that included investment (specific leaf area), morphology, chlorophyll fluorescence, pigment content, and gene expression. Traits and differentially expressed genes that changed with elevation in one species often showed little changes in the other species, or changed in the opposite direction. As evidence of adaptive divergence, both species performed better at their native site and better than the species from the other habitat. Adaptive divergence is, therefore, associated with the evolution of distinct plastic responses to environmental variation, despite these two species sharing a recent common ancestor.

Optimal composition of chloride cells for osmoregulation in a randomly fluctuating environment

Fish live in water with a different osmotic pressure from that in the body. Their gills have chloride cells that transport ions to maintain an appropriate level of osmotic pressure in the body. The direction of ion transport is different between seawater and freshwater. There are two types of chloride cells that specialize in unidirectional transport and generalist cells that can switch their function quickly in response to environmental salinity. In species that experience salinity changes throughout life (euryhaline species), individuals may replace some chloride cells with cells of different types upon a sudden change in environmental salinity. In this paper, we develop a dynamic optimization model for the chloride cell composition of an individual living in an environment with randomly fluctuating salinity. The optimal solution is to minimize the sum of the workload of chloride cells in coping with the difference in osmotic pressure, the maintenance cost, and the temporal cost due to environmental change. The optimal fraction of generalist chloride cells increases with the frequency of salinity changes and the time needed for new cells to be fully functional but decreases with excess maintenance cost.

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.

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.

Restoring ancestral phenotypes is a general pattern in gene expression evolution during adaptation to new environments in Tribolium castaneum

Plasticity and evolution are two processes allowing populations to respond to environmental changes, but how both are related and impact each other remains controversial. We studied plastic and evolutionary responses in gene expression of Tribolium castaneum after exposure of the beetles to new environments that differed from ancestral conditions in temperature, humidity or both. Using experimental evolution with 10 replicated lines per condition, we were able to demonstrate adaptation after 20 generations. We measured whole-transcriptome gene expression with RNA-sequencing to infer evolutionary and plastic changes. We found more evidence for changes in mean expression (shift in the intercept of reaction norms) in adapted lines than for changes in plasticity (shifts in slopes). Plasticity was mainly preserved in selected lines and was responsible for a large part of the phenotypic divergence in expression between ancestral and new conditions. However, we found that genes with the largest evolutionary changes in expression also evolved reduced plasticity and often showed expression levels closer to the ancestral stage. Results obtained in the three different conditions were similar, suggesting that restoration of ancestral expression levels during adaptation is a general evolutionary pattern. With a larger sample in the most stressful condition, we were able to detect a positive correlation between the proportion of genes with reversion of the ancestral plastic response and mean fitness per selection line.

Chronic stress, physiological adaptation and developmental programming of the neuroendocrine stress system

Chronic stress undermines physical and mental health, in part via dysregulation of the neuroendocrine stress system. Key to understand this dysregulation is recognizing that the problem is not stress per se, but rather its chronicity. The optimally functioning stress system is highly dynamic, and negative feedback regulation enforces transient responses to acute stressors. Chronic stress overrides this, and adaptation to the chronicity can result in persistent dysregulation by altering sensitivity thresholds critical for control of system dynamics. Such adaptation involves plasticity within the central nervous system (CNS) as well as epigenetic regulation. When it occurs during development, it can have persistent effects on neuroendocrine regulation. Understanding how chronic stress programs development of the neuroendocrine stress system requires elucidation of stress-responsive gene regulatory networks that control CNS plasticity and development.

Higher than … or lower than ….? Evidence for the validity of the extrapolation of laboratory toxicity test results to predict the effects of chemicals and ionising radiation in the field

Single species laboratory tests and associated species sensitivity distributions (SSDs) that utilise the resulting data can make a key contribution to efforts to prospective hazard assessments for pesticides, biocides, metals and ionising radiation for research and regulatory risk assessment. An assumption that underlies the single species based toxicity testing approach when combined in SSD models is that the assessments of sensitivities to chemical and ionising radiation measured across a range of species in the laboratory can inform on the likely effects on communities present in the field. Potential issues with the validity of this assumption were already recognised by Van Straalen and Denneman (1989) in their landmark paper on the SSD methodology. In this work, they identified eight major factors that could potentially compromise the extrapolation of laboratory toxicity data to the field. Factors covered a range of issues related to differences in chemistry (e.g. bioavailability, mixtures); environmental conditions (optimal, variable), ecological (compensatory, time-scale) and population genetic structure (adaptation, meta-population dynamics). This paper outlines the evidence pertaining to the influence of these different factors on toxicity in the laboratory as compared to the field focussing especially on terrestrial ecosystems. Through radiological and ecotoxicological research, evidence of the influence of each factor on the translation of observed toxicity from the laboratory to field is available in all cases. The importance of some factors, such as differences in chemical bioavailability between laboratory tests and the field and the ubiquity of exposure to mixtures is clearly established and has some relevance to radiological protection. However, other factors such as the differences in test conditions (optimal vs sub-optimal) and the development of tolerance may be relevant on a case by case basis. When SSDs generated from laboratory tests have been used to predict chemical and ionising radiation effects in the field, results have indicated that they may often seem to under-predict impacts, although this may also be due to other factors such as the effects of other non-chemical stressors also affecting communities at polluted sites. A better understanding of the main factors affecting this extrapolation can help to reduce uncertainty during risk assessment.

Dynamics of Cadmium Acclimation in Daphnia pulex: Linking Fitness Costs, Cross-Tolerance, and Hyper-Induction of Metallothionein

Acclimation increases tolerance to stress in individuals but is assumed to contribute fitness costs when the stressor is absent, though data supporting this widely held claim are sparse. Therefore, using clonal (i.e., genetically identical) cultures of Daphnia pulex, we isolated the contributions of acclimation to the regulation of the metal response gene, metallothionein 1 (MT1), and defined the reproductive benefits and costs of cadmium (Cd)-acclimation. Daphnia pulex were exposed for 50 parthenogenetic generations to environmentally realistic levels (1 μg Cd/L), and tolerance to Cd and other metals assessed during this period via standard toxicity tests. These tests revealed (1) increased tolerance to Cd compared to genetically identical nonacclimated cultures, (2) fitness costs in Cd-acclimated Daphnia when Cd was removed, and (3) cross-tolerance of Cd-acclimated Daphnia to zinc and silver, but not arsenic, thereby defining a functional role for metallothionein. Indeed, Cd-acclimated clones had significantly higher expression of MT1 mRNA than nonacclimated clones, when Cd exposed. Both the enhanced induction of MT1 and tolerant phenotype were rapidly lost when Cd was removed (1-2 generations), which is further evidence of acclimation costs. These findings provide evidence for the widely held view that acclimation is costly and are important for investigating evolutionary principles of genetic assimilation and the survival mechanisms of natural populations that face changing environments.

Current evidence for a role of epigenetic mechanisms in response to ionizing radiation in an ecotoxicological context

The issue of potential long-term or hereditary effects for both humans and wildlife exposed to low doses (or dose rates) of ionising radiation is a major concern. Chronic exposure to ionising radiation, defined as an exposure over a large fraction of the organism's lifespan or even over several generations, can possibly have consequences in the progeny. Recent work has begun to show that epigenetics plays an important role in adaptation of organisms challenged to environmental stimulae. Changes to so-called epigenetic marks such as histone modifications, DNA methylation and non-coding RNAs result in altered transcriptomes and proteomes, without directly changing the DNA sequence. Moreover, some of these environmentally-induced epigenetic changes tend to persist over generations, and thus, epigenetic modifications are regarded as the conduits for environmental influence on the genome. Here, we review the current knowledge of possible involvement of epigenetics in the cascade of responses resulting from environmental exposure to ionising radiation. In addition, from a comparison of lab and field obtained data, we investigate evidence on radiation-induced changes in the epigenome and in particular the total or locus specific levels of DNA methylation. The challenges for future research and possible use of changes as an early warning (biomarker) of radiosensitivity and individual exposure is discussed. Such a biomarker could be used to detect and better understand the mechanisms of toxic action and inter/intra-species susceptibility to radiation within an environmental risk assessment and management context.

Profiling microRNA expression in Atlantic killifish (Fundulus heteroclitus) gill and responses to arsenic and hyperosmotic stress

The Atlantic killifish (Fundulus heteroclitus), native to estuarine areas of the Atlantic coast of the United States, has become a valuable ecotoxicological model as a result of its ability to acclimate to rapid environmental changes and adapt to polluted habitats. MicroRNAs (miRNAs) are highly conserved small RNAs that regulate gene expression and play critical roles in stress responses in a variety of organisms. Global miRNA expression in killifish and the potential roles miRNA have in environmental acclimation have yet to be characterized. Accordingly, we profiled miRNA expression in killifish gill for the first time and identified a small group of highly expressed, well-conserved miRNAs as well as 16 novel miRNAs not yet identified in other organisms. Killifish respond to large fluctuations in salinity with rapid changes in gene expression and protein trafficking to maintain osmotic balance, followed by a secondary phase of gene and protein expression changes that enable remodeling of the gills. Arsenic, a major environmental toxicant, was previously shown to inhibit gene expression responses in killifish gill, as well the ability of killifish to acclimate to a rapid increase in salinity. Thus, we examined the individual and combined effects of salinity and arsenic on miRNA expression in killifish gill. Using small RNA sequencing, we identified 270 miRNAs expressed in killifish, and found that miR-135b was differentially expressed in response to arsenic and at 24 h following transfer to salt water. Predicted targets of miR-135b are involved in ion transport, cell motility and migration, GTPase mediated signal transduction and organelle assembly. Consistent with previous studies of these two environmental stressors, we found a significant interaction (i.e., arsenic dependent salinity effect), whereby killifish exposed to arsenic exhibited an opposite response in miR-135b expression at 24 h post hyperosmotic challenge compared to controls. By examining mRNA expression of predicted miRNA targets during salinity acclimation and arsenic exposure, we found that miR-135b targets were significantly more likely to decrease during salinity acclimation than non-targets. Our identification of a significant interaction effect of arsenic and salinity on miR-135b expression supports the hypothesis that arsenic alters upstream regulators of stress response networks, which may adversely affect the killifish response to osmotic stress. The characterization of miRNAs in this ecotoxicological model will be a valuable resource for future studies investigating the role of miRNAs in response to environmental stress.

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.

Functional and evolutionary genomic inferences in Populus through genome and population sequencing of American and European aspen

We performed de novo, full-genome sequence analysis of two Populus species, North American quaking and Eurasian trembling aspen, that contain striking levels of genetic variation. Our results showed that positive and negative selection broadly affects patterns of genomic variation, but to varying degrees across coding and noncoding regions. The strength of selection and rates of sequence divergence were strongly related to differences in gene expression and coexpression network connectivity. These results highlight the importance of both positive and negative selection in shaping genome-wide levels of genetic variation in an obligately outcrossing, perennial plant. The resources we present establish aspens as a powerful study system enabling future studies for understanding the genomic determinants of adaptive evolution.

The Populus genus is one of the major plant model systems, but genomic resources have thus far primarily been available for poplar species, and primarily Populus trichocarpa (Torr. & Gray), which was the first tree with a whole-genome assembly. To further advance evolutionary and functional genomic analyses in Populus, we produced genome assemblies and population genetics resources of two aspen species, Populus tremula L. and Populus tremuloides Michx. The two aspen species have distributions spanning the Northern Hemisphere, where they are keystone species supporting a wide variety of dependent communities and produce a diverse array of secondary metabolites. Our analyses show that the two aspens share a similar genome structure and a highly conserved gene content with P. trichocarpa but display substantially higher levels of heterozygosity. Based on population resequencing data, we observed widespread positive and negative selection acting on both coding and noncoding regions. Furthermore, patterns of genetic diversity and molecular evolution in aspen are influenced by a number of features, such as expression level, coexpression network connectivity, and regulatory variation. To maximize the community utility of these resources, we have integrated all presented data within the PopGenIE web resource (PopGenIE.org).

Comparative Studies of Gene Expression Kinetics: Methodologies and Insights on Development and Evolution

Across the animal kingdom, embryos of closely related species show high morphological similarity despite genetic and environmental distances. Deciphering the molecular mechanisms that underlie morphological conservation and those that support embryonic adaptation are keys to understand developmental robustness and evolution. Comparative studies of developmental gene regulatory networks can track the genetic changes that lead to evolutionary novelties. However, these studies are limited to a relatively small set of genes and demand extensive experimental efforts. An alternative approach enabled by next-generation sequencing, is to compare the expression kinetic of large sets of genes between different species. The advantages of these comparisons are that they can be done relatively easily, for any species and they provide information of all expressed genes. The challenge in these experiments is to compare the kinetic profiles of thousands of genes between species that develop in different rates. Here we review recent comparative studies that tackled the challenges of accurate staging and large-scale analyses using different computational approaches. These studies reveal how correct temporal scaling exposes the striking conservation of developmental gene expression between morphologically similar species. Different clustering approaches are used to address various comparative questions and identify the conservation and divergence of large gene sets. We discuss the unexpected contribution of housekeeping genes to the interspecies correlations and how this contribution distorts the hourglass pattern generated by developmental genes. Overall, we demonstrate how comparative studies of gene expression kinetics can provide novel insights into the developmental constraints and plasticity that shape animal body plans.

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.

Arsenic Reduces Gene Expression Response to Changing Salinity in Killifish

Toxicogenomic approaches can detect and classify adverse interactions between environmental toxicants and other environmental stressors but require more complex experimental designs and analytical approaches. Here we use novel toxicogenomic techniques to analyze the effect of arsenic exposure in wild killifish populations acclimating to changing salinity. Fish from three populations were acclimated to full strength seawater and transferred to fresh water for 1 or 24 h. Linear models of gene expression in gill tissue identified 31 genes that responded to osmotic shock at 1 h and 178 genes that responded at 24 h. Arsenic exposure (100 μg/L) diminished the responses (reaction norms) of these genes by 22% at 1 h ( p = 1.0 × 10^-6) and by 10% at 24 h ( p = 3.0 × 10^-10). Arsenic also significantly reduced gene coregulation in gene regulatory networks ( p = 0.002, paired Levene's test), and interactions between arsenic and salinity acclimation were uniformly antagonistic at the biological pathway level ( p < 0.05, binomial test). Arsenic's systematic interference with gene expression reaction norms was validated in a mouse multistressor experiment, demonstrating the ability of these toxicogenomic approaches to identify biologically relevant adverse interactions between environmental toxicants and other environmental stressors.

Nitrate: An Environmental Endocrine Disruptor? A Review of Evidence and Research Needs

Nitrate is heavily used as an agricultural fertilizer and is today a ubiquitous environmental pollutant. Environmental endocrine effects caused by nitrate have received increasing attention over the last 15 years. Nitrate is hypothesized to interfere with thyroid and steroid hormone homeostasis and developmental and reproductive end points. The current review focuses on aquatic ecotoxicology with emphasis on field and laboratory controlled in vitro and in vivo studies. Furthermore, nitrate is just one of several forms of nitrogen that is present in the environment and many of these are quickly interconvertible. Therefore, the focus is additionally confined to the oxidized nitrogen species (nitrate, nitrite and nitric oxide). We reviewed 26 environmental toxicology studies and our main findings are (1) nitrate has endocrine disrupting properties and hypotheses for mechanisms exist, which warrants for further investigations; (2) there are issues determining actual nitrate-speciation and abundance is not quantified in a number of studies, making links to speciation-specific effects difficult; and (3) more advanced analytical chemistry methodologies are needed both for exposure assessment and in the determination of endocrine biomarkers.

The roles of plasticity and evolutionary change in shaping gene expression variation in natural populations of extremophile fish

The notorious plasticity of gene expression responses and the complexity of environmental gradients complicate the identification of adaptive differences in gene regulation among populations. We combined transcriptome analyses in nature with common-garden and exposure experiments to establish cause-effect relationships between the presence of a physiochemical stressor and expression differences, as well as to test how evolutionary change and plasticity interact to shape gene expression variation in natural systems. We studied two evolutionarily independent population pairs of an extremophile fish (Poecilia mexicana) living in toxic, hydrogen sulphide (H₂ S)-rich springs and adjacent nontoxic habitats and assessed genomewide expression patterns of wild-caught and common-garden-raised individuals exposed to different concentrations of H₂ S. We found that 7.7% of genes that were differentially expressed between sulphidic and nonsulphidic ecotypes remained differentially expressed in the laboratory, indicating that sources of selection other than H₂ S-or plastic responses to other environmental factors-contribute substantially to gene expression patterns observed in the wild. Concordantly differentially expressed genes in the wild and the laboratory were primarily associated with H₂ S detoxification, sulphur processing and metabolic physiology. While shared, ancestral plasticity played a minor role in shaping gene expression variation observed in nature, we documented evidence for evolved population differences in the constitutive expression as well as the H₂ S inducibility of candidate genes. Mechanisms underlying gene expression variation also varied substantially across the two ecotype pairs. These results provide a springboard for studying evolutionary modifications of gene regulatory mechanisms that underlie expression variation in locally adapted populations.

Parallel embryonic transcriptional programs evolve under distinct constraints and may enable morphological conservation amidst adaptation

Embryonic development evolves by balancing stringent morphological constraints with genetic and environmental variation. The design principle that allows developmental transcriptional programs to conserve embryonic morphology while adapting to environmental changes is still not fully understood. To address this fundamental challenge, we compare developmental transcriptomes of two sea urchin species, Paracentrotus lividus and Strongylocentrotus purpuratus, that shared a common ancestor about 40 million years ago and are geographically distant yet show similar morphology. We find that both developmental and housekeeping genes show highly dynamic and strongly conserved temporal expression patterns. The expression of other gene sets, including homeostasis and response genes, show divergent expression which could result from either evolutionary drift or adaptation to local environmental conditions. The interspecies correlations of developmental gene expressions are highest between morphologically similar developmental time points whereas the interspecies correlations of housekeeping gene expression are high between all the late zygotic time points. Relatedly, the position of the phylotypic stage varies between these two groups of genes: developmental gene expression shows highest conservation at mid-developmental stage, in agreement with the hourglass model while the conservation of housekeeping genes keeps increasing with developmental time. When all genes are combined, the relationship between conservation of gene expression and morphological similarity is partially masked by housekeeping genes and genes with diverged expression. Our study illustrates various transcriptional programs that coexist in the developing embryo and evolve under different constraints. Apparently, morphological constraints underlie the conservation of developmental gene expression while embryonic fitness requires the conservation of housekeeping gene expression and the species-specific adjustments of homeostasis gene expression. The distinct evolutionary forces acting on these transcriptional programs enable the conservation of similar body plans while allowing adaption.

The interplay between phenotypic and ontogenetic plasticities can be assessed using reaction-diffusion models : The case of Pseudoplatystoma fishes

Every morphological, behavioral, or even developmental character expression of living beings is coded in its genotype and is expressed in its phenotype. Nevertheless, the interplay between phenotypic and ontogenetic plasticities, that is, the capability to manifest trait variations, is a current field of research that needs morphometric, numerical, or even mathematical modeling investigations. In the present work, we are searching for a phenotypic index able to identify the underlying correlation among phenotypic, ontogenetic, and geographic distribution of the evolutionary development of species of the same genus. By studying the case of Pseudoplatystoma fishes, we use their skin patterns as an auxiliary trait that can be reproduced by means of a reaction diffusion (RD) model. From this model, we infer the phenotypic index in terms of one of the parameters appearing in the mathematical equations. To achieve this objective, we perform extensive numerical simulations and analysis of the model equations and link the parameter variations with different environmental and physicochemical conditions in which the individuals develop, and which may be regulated by the ontogenetic plasticity of the species. Our numerical study indicates that the patterns predicted by a set of reaction diffusion equations are not uniquely determined by the value of the parameters of the equation, but also depend on how the process is initiated and on the spatial distribution of values of these parameters. These factors are therefore significant, since they show that an individual's growth dynamics and apparent secondary transport processes, like advection, can be determinant for the alignment of motifs in a skin pattern. Our results allow us to discern the correlation between phenotypic, ontogenetic, and geographic distribution of the different species of Pseudoplatystoma fishes, thus indicating that RD models represent a useful taxonomic tool able to quantify evolutionary indexes.

The landscape of extreme genomic variation in the highly adaptable Atlantic killifish

Understanding and predicting the fate of populations in changing environments require knowledge about the mechanisms that support phenotypic plasticity and the adaptive value and evolutionary fate of genetic variation within populations. Atlantic killifish (Fundulus heteroclitus) exhibit extensive phenotypic plasticity that supports large population sizes in highly fluctuating estuarine environments. Populations have also evolved diverse local adaptations. To yield insights into the genomic variation that supports their adaptability, we sequenced a reference genome and 48 additional whole genomes from a wild population. Evolution of genes associated with cell cycle regulation and apoptosis is accelerated along the killifish lineage, which is likely tied to adaptations for life in highly variable estuarine environments. Genome-wide standing genetic variation, including nucleotide diversity and copy number variation, is extremely high. The highest diversity genes are those associated with immune function and olfaction, whereas genes under greatest evolutionary constraint are those associated with neurological, developmental, and cytoskeletal functions. Reduced genetic variation is detected for tight junction proteins, which in killifish regulate paracellular permeability that supports their extreme physiological flexibility. Low-diversity genes engage in more regulatory interactions than high-diversity genes, consistent with the influence of pleiotropic constraint on molecular evolution. High genetic variation is crucial for continued persistence of species given the pace of contemporary environmental change. Killifish populations harbor among the highest levels of nucleotide diversity yet reported for a vertebrate species, and thus may serve as a useful model system for studying evolutionary potential in variable and changing environments.

Canalization of gene expression is a major signature of regulatory cold adaptation in temperate Drosophila melanogaster

BACKGROUND

Transcriptome analysis may provide means to investigate the underlying genetic causes of shared and divergent phenotypes in different populations and help to identify potential targets of adaptive evolution. Applying RNA sequencing to whole male Drosophila melanogaster from the ancestral tropical African environment and a very recently colonized cold-temperate European environment at both standard laboratory conditions and following a cold shock, we seek to uncover the transcriptional basis of cold adaptation.

RESULTS

In both the ancestral and the derived populations, the predominant characteristic of the cold shock response is the swift and massive upregulation of heat shock proteins and other chaperones. Although we find ~25 % of the genome to be differentially expressed following a cold shock, only relatively few genes (n = 16) are up- or down-regulated in a population-specific way. Intriguingly, 14 of these 16 genes show a greater degree of differential expression in the African population. Likewise, there is an excess of genes with particularly strong cold-induced changes in expression in Africa on a genome-wide scale.

CONCLUSIONS

The analysis of the transcriptional cold shock response most prominently reveals an upregulation of components of a general stress response, which is conserved over many taxa and triggered by a plethora of stressors. Despite the overall response being fairly similar in both populations, there is a definite excess of genes with a strong cold-induced fold-change in Africa. This is consistent with a detrimental deregulation or an overshooting stress response. Thus, the canalization of European gene expression might be responsible for the increased cold tolerance of European flies.

Environmental variability and the evolution of the glucocorticoid receptor (Nr3c1) in African starlings

One of the primary ways that organisms cope with environmental change is through regulation of the hypothalamo-pituitary-adrenal (HPA) axis, the neuroendocrine system that controls reactions to stress. Variation in genes regulating the HPA axis - particularly the glucocorticoid receptor - may facilitate adaptation to changing climatic conditions by altering expression. Here we examine signatures of selection on the glucocorticoid receptor gene (Nr3c1) in African starlings that inhabit a range of environments, including those with variable climatic conditions. To investigate potential adaptive mechanisms underlying the vertebrate stress response, we sequence the Nr3c1 gene in 27 species of African starlings. Although we find some evidence of positive selection, substitution rate is negatively correlated with variance in precipitation. This suggests climatic cycling in sub-Saharan Africa may have resulted in lower substitution rates to maintain a successful coping strategy. When environmental conditions fluctuate rapidly, variation in the strength of purifying selection can explain evolutionary rate variation.

Plasticity of the Berry Ripening Program in a White Grape Variety

Grapevine (Vitis vinifera L.) is considered one of the most environmentally sensitive crops and is characterized by broad phenotypic plasticity, offering important advantages such as the large range of different wines that can be produced from the same cultivar, and the adaptation of existing cultivars to diverse growing regions. The uniqueness of berry quality traits reflects complex interactions between the grapevine plant and the combination of natural factors and human cultural practices which leads to the expression of wine typicity. Despite the scientific and commercial importance of genotype interactions with growing conditions, few studies have characterized the genes and metabolites directly involved in this phenomenon. Here, we used two large-scale analytical approaches to explore the metabolomic and transcriptomic basis of the broad phenotypic plasticity of Garganega, a white berry variety grown at four sites characterized by different pedoclimatic conditions (altitudes, soil texture, and composition). These conditions determine berry ripening dynamics in terms of sugar accumulation and the abundance of phenolic compounds. Multivariate analysis unraveled a highly plastic metabolomic response to different environments, especially the accumulation of hydroxycinnamic and hydroxybenzoic acids and flavonols. Principal component analysis (PCA) revealed that the four sites strongly affected the berry transcriptome allowing the identification of environmentally-modulated genes and the plasticity of commonly-modulated transcripts at different sites. Many genes that control transcription, translation, transport, and carbohydrate metabolism showed different expression depending on the environmental conditions, indicating a key role in the observed transcriptomic plasticity of Garganega berries. Interestingly, genes representing the phenylpropanoid/flavonoid pathway showed plastic responses to the environment mirroring the accumulation of the corresponding metabolites. The comparison of Garganega and Corvina berries showed that the metabolism of phenolic compounds is more plastic in ripening Garganega berries under different pedoclimatic conditions.

Functional genomics to assess biological responses to marine pollution at physiological and evolutionary timescales: toward a vision of predictive ecotoxicology

Marine pollution is ubiquitous, and is one of the key factors influencing contemporary marine biodiversity worldwide. To protect marine biodiversity, how do we surveil, document and predict the short- and long-term impacts of pollutants on at-risk species? Modern genomics tools offer high-throughput, information-rich and increasingly cost-effective approaches for characterizing biological responses to environmental stress, and are important tools within an increasing sophisticated kit for surveiling and assessing impacts of pollutants on marine species. Through the lens of recent research in marine killifish, we illustrate how genomics tools may be useful for screening chemicals and pollutants for biological activity and to reveal specific mechanisms of action. The high dimensionality of transcriptomic responses enables their usage as highly specific fingerprints of exposure, and these fingerprints can be used to diagnose environmental problems. We also emphasize that molecular pathways recruited to respond at physiological timescales are the same pathways that may be targets for natural selection during chronic exposure to pollutants. Gene complement and sequence variation in those pathways can be related to variation in sensitivity to environmental pollutants within and among species. Furthermore, allelic variation associated with evolved tolerance in those pathways could be tracked to estimate the pace of environmental health decline and recovery. We finish by integrating these paradigms into a vision of how genomics approaches could anchor a modernized framework for advancing the predictive capacity of environmental and ecotoxicological science.

Genetic approaches in comparative and evolutionary physiology

Whole animal physiological performance is highly polygenic and highly plastic, and the same is generally true for the many subordinate traits that underlie performance capacities. Quantitative genetics, therefore, provides an appropriate framework for the analysis of physiological phenotypes and can be used to infer the microevolutionary processes that have shaped patterns of trait variation within and among species. In cases where specific genes are known to contribute to variation in physiological traits, analyses of intraspecific polymorphism and interspecific divergence can reveal molecular mechanisms of functional evolution and can provide insights into the possible adaptive significance of observed sequence changes. In this review, we explain how the tools and theory of quantitative genetics, population genetics, and molecular evolution can inform our understanding of mechanism and process in physiological evolution. For example, lab-based studies of polygenic inheritance can be integrated with field-based studies of trait variation and survivorship to measure selection in the wild, thereby providing direct insights into the adaptive significance of physiological variation. Analyses of quantitative genetic variation in selection experiments can be used to probe interrelationships among traits and the genetic basis of physiological trade-offs and constraints. We review approaches for characterizing the genetic architecture of physiological traits, including linkage mapping and association mapping, and systems approaches for dissecting intermediary steps in the chain of causation between genotype and phenotype. We also discuss the promise and limitations of population genomic approaches for inferring adaptation at specific loci. We end by highlighting the role of organismal physiology in the functional synthesis of evolutionary biology.

Considerations for the use of transcriptomics in identifying the ‘genes that matter’ for environmental adaptation

Transcriptomics has emerged as a powerful approach for exploring physiological responses to the environment. However, like any other experimental approach, transcriptomics has its limitations. Transcriptomics has been criticized as an inappropriate method to identify genes with large impacts on adaptive responses to the environment because: (1) genes with large impacts on fitness are rare; (2) a large change in gene expression does not necessarily equate to a large effect on fitness; and (3) protein activity is most relevant to fitness, and mRNA abundance is an unreliable indicator of protein activity. In this review, these criticisms are re-evaluated in the context of recent systems-level experiments that provide new insight into the relationship between gene expression and fitness during environmental stress. In general, these criticisms remain valid today, and indicate that exclusively using transcriptomics to screen for genes that underlie environmental adaptation will overlook constitutively expressed regulatory genes that play major roles in setting tolerance limits. Standard practices in transcriptomic data analysis pipelines may also be limiting insight by prioritizing highly differentially expressed and conserved genes over those genes that undergo moderate fold-changes and cannot be annotated. While these data certainly do not undermine the continued and widespread use of transcriptomics within environmental physiology, they do highlight the types of research questions for which transcriptomics is best suited and the need for more gene functional analyses. Such information is pertinent at a time when transcriptomics has become increasingly tractable and many researchers may be contemplating integrating transcriptomics into their research programs.

Extending the integrated phenotype: covariance and correlation in plasticity of behavioural traits

In the field of behavioural ecology there has been a longstanding interest in the evolution of phenotypic plasticity, as plasticity in behavioural traits such as foraging, mating, and reproduction governs the capacity of organisms to cope with environmental variability. In this paper we highlight the need for an integrated perspective to phenotypic plasticity of traits, taking into account covariation among plastic responses of traits. We discuss new perspectives on the importance of integrated plasticity of traits for adaptive behavioural strategies. We review empirical evidence for correlated plasticity across behavioural traits in insects, for example, through genetic correlation, a shared pool of resources or dependency on a common developmental path. Taking on an integrated plasticity perspective, we suggest an alternative explanation for the apparent lack of costs of plasticity, and offer a better understanding of the relative benefits of plasticity or canalization of traits.

A novel variant of aquaporin 3 is expressed in killifish (Fundulus heteroclitus) intestine

Killifish (Fundulus heteroclitus) are euryhaline teleosts that are widely used in environmental and toxicological studies, and they are tolerant to arsenic, in part due to very low assimilation of arsenic from the environment. The mechanism of arsenic uptake by the intestine, a major route of arsenic uptake in humans is unknown. Thus, the goal of this study was to determine if aquaglyceroporins (AQPs), which transport water and other small molecules including arsenite across cell membranes, are expressed in the killifish intestine, and whether AQP expression is affected by osmotic stress. Through RT-PCR and sequence analysis of PCR amplicons, we demonstrated that the intestine expresses kfAQP3a and kfAQP3b, two previously identified variants, and also identified a novel variant of killifish AQP3 (kfAQP3c) in the intestine. The variants likely represent alternate splice forms. A BLAST search of the F. heteroclitus reference genome revealed that the AQP3 gene resides on a single locus, while an alignment of the AQP3 sequence among 384 individuals from eight population ranging from Rhode Island to North Carolina revealed that its coding sequence was remarkably conserved with no fixed polymorphism residing in the region that distinguishes these variants. We further demonstrate that the novel variant transports arsenite into HEK293T cells. Whereas kfAQP3a, which does not transport arsenite, was expressed in both freshwater (FW) and saltwater (SW) acclimated fish, kfAQP3b, an arsenic transporter, was expressed only in FW acclimated fish, and kfAQP3c was expressed only in SW acclimated fish. Thus, we have identified a novel, putative splice variant of kfAQP3, kfAQP3c, which transports arsenic and is expressed only in SW acclimated fish.