Towards a gene regulatory network perspective on phenotypic plasticity, genetic accommodation and genetic assimilation

Is the genetic architecture of behavior exceptionally complex?

Are traits with high levels of plasticity more complex in their genetic architecture, as they can be modulated by numerous different environmental inputs? Many authors have assumed that behavioral traits, in part because they are highly plastic, have an exceptionally complex genetic basis. We quantitatively summarized data from 31 genome-wide association studies (GWAS) and 87 traits in Drosophila melanogaster and found no evidence that behavioral traits have fundamental differences in the number of single-nucleotide polymorphisms or the significance or effect size of those associations, compared with nonbehavioral (morphological or physiological) traits. We suggest the assertion that behavioral traits are inherently more complex on a genetic basis compared with other types of traits should not be assumed true, and merits further investigation.

Hybridization and gene expression: Beyond differentially expressed genes

Gene expression has a key role in reproductive isolation, and studies of hybrid gene expression have identified mechanisms causing hybrid sterility. Here, we review the evidence for altered gene expression following hybridization and outline the mechanisms shown to contribute to altered gene expression in hybrids. Transgressive gene expression, transcending that of both parental species, is pervasive in early generation sterile hybrids, but also frequently observed in viable, fertile hybrids. We highlight studies showing that hybridization can result in transgressive gene expression, also in established hybrid lineages or species. Such extreme patterns of gene expression in stabilized hybrid taxa suggest that altered hybrid gene expression may result in hybridization-derived evolutionary novelty. We also conclude that while patterns of misexpression in hybrids are well documented, the understanding of the mechanisms causing misexpression is lagging. We argue that jointly assessing differences in cell composition and cell-specific changes in gene expression in hybrids, in addition to assessing changes in chromatin and methylation, will significantly advance our understanding of the basis of altered gene expression. Moreover, uncovering to what extent evolution of gene expression results in altered expression for individual genes, or entire networks of genes, will advance our understanding of how selection moulds gene expression. Finally, we argue that jointly studying the dual roles of altered hybrid gene expression, serving both as a mechanism for reproductive isolation and as a substrate for hybrid ecological adaptation, will lead to significant advances in our understanding of the evolution of gene expression.

One genome, multiple phenotypes: decoding the evolution and mechanisms of environmentally induced developmental plasticity in insects

Plasticity in developmental processes gives rise to remarkable environmentally induced phenotypes. Some of the most striking and well-studied examples of developmental plasticity are seen in insects. For example, beetle horn size responds to nutritional state, butterfly eyespots are enlarged in response to temperature and humidity, and environmental cues also give rise to the queen and worker castes of eusocial insects. These phenotypes arise from essentially identical genomes in response to an environmental cue during development. Developmental plasticity is taxonomically widespread, affects individual fitness, and may act as a rapid-response mechanism allowing individuals to adapt to changing environments. Despite the importance and prevalence of developmental plasticity, there remains scant mechanistic understanding of how it works or evolves. In this review, we use key examples to discuss what is known about developmental plasticity in insects and identify fundamental gaps in the current knowledge. We highlight the importance of working towards a fully integrated understanding of developmental plasticity in a diverse range of species. Furthermore, we advocate for the use of comparative studies in an evo-devo framework to address how developmental plasticity works and how it evolves.

Cis- and Trans-variations of Stearoyl-CoA Desaturase Provide New Insights into the Mechanisms of Diverged Pattern of Phenotypic Plasticity for Temperature Adaptation in Two Congeneric Oyster Species

The evolution of phenotypic plasticity plays an essential role in adaptive responses to climate change; however, its regulatory mechanisms in marine organisms which exhibit high phenotypic plasticity still remain poorly understood. The temperature-responsive trait oleic acid content and its major gene stearoyl-CoA desaturase (Scd) expression have diverged in two allopatric congeneric oyster species, cold-adapted Crassostrea gigas and warm-adapted Crassostrea angulata. In this study, genetic and molecular methods were used to characterize fatty acid desaturation and membrane fluidity regulated by oyster Scd. Sixteen causative single-nucleotide polymorphisms (SNPs) were identified in the promoter/cis-region of the Scd between wild C. gigas and C. angulata. Further functional experiments showed that an SNP (g.-333C [C. gigas allele] >T [C. angulata allele]) may influence Scd transcription by creating/disrupting the binding motif of the positive trans-factor Y-box factor in C. gigas/C. angulata, which mediates the higher/lower constitutive expression of Scd in C. gigas/C. angulata. Additionally, the positive trans-factor sterol-regulatory element-binding proteins (Srebp) were identified to specifically bind to the promoter of Scd in both species, and were downregulated during cold stress in C. gigas compared to upregulated in C. angulata. This partly explains the relatively lower environmental sensitivity (plasticity) of Scd in C. gigas. This study serves as an experimental case to reveal that both cis- and trans-variations shape the diverged pattern of phenotypic plasticity, which provides new insights into the formation of adaptive traits and the prediction of the adaptive potential of marine organisms to future climate change.

Where are the seeds? Lack of floral morphs prevent seed production by the tristylous Pontederia cordata in South Africa

Reproduction is a crucial part of the successful establishment and spread of an invasive species. Invasive plants often produce seeds prolifically to spread into new ranges, yet the invasive macrophyte, Pontederia cordata L., does not appear to produce seeds in South Africa, limiting its invasive potential. Here, we aimed to determine what limits seed production of the tristylous P. cordata in South Africa, where it is widespread with impacts on the ecology of wetlands it invades, South Africa. We measured floral traits and pollen grain size from populations throughout the invasive range in South Africa to determine the relative proportion of tristylous morphs. We speculated that the absence of specialist native pollinators in the invasive range may be responsible for the absence of sexual reproduction and thus conducted a pollination study to determine whether flowers were visited. Thereafter, we hand pollinated 8865 flowers to conclude whether P. cordata exhibited an incompatibility system, which prevented seed production. The floral traits and pollen grain measurements were similar to those reported for short-morphed flowers from the native range. The pollination study confirmed the absence of specialist insect visitors, while the hand-pollination experiments resulted in no seed production. Only short-morphed plants are present in South Africa, and the illegitimate pollination of short-morphed plants prevents seed production. Vegetative spread through rhizome production is thus responsible for the invasion of P. cordata throughout South Africa. These findings suggest that control programs should target the plants' rhizomes to prevent and reduce spread. More importantly, preventing the introduction of medium- and long-morphed plants into South Africa is crucial to preclude P. cordata from producing seeds and enhancing invasion.

Shining light on the transcriptome: Molecular regulatory networks leading to a fast-growth phenotype by continuous light in an environmentally sensitive teleost (Atherinopsidae)

Photoperiod can profoundly affect the physiology of teleost fish, including accelerated growth here defined as "fast growth phenotypes". However, molecular regulatory networks (MRNs) and biological processes being affected by continuous illumination and which allow some teleost species evident plasticity to thrive under this condition are not yet clear. Therefore, to provide a broad perspective of such mechanisms, Chirostoma estor fish were raised and sampled for growth under a simulated control (LD) 12 h Light: 12 h Dark or a continuous illumination (LL) 24 h Light: 0 h Dark since fertilization. The experiment lasted 12 weeks after hatching (wah), the time at which fish were sampled for growth, length, and whole-body cortisol levels. Additionally, 3 heads of fish from each treatment were used to perform a de novo transcriptome analysis using Next-Generation Sequencing. Fish in LL developed the fast growth phenotype with significant differences visible at 4 wah and gained 66% more mass by 12 wah than LD fish. Cortisol levels under LL were below basal levels at all times compared to fish in LD, suggesting circadian dysregulation effects. A strong effect of LL was observed in samples with a generalized down-regulation of genes except for Reactive Oxygen Species responses, genome stability, and growth biological processes. To our knowledge, this work is the first study using a transcriptomic approach to understand environmentally sensitive MRNs that mediate phenotypic plasticity in fish submitted to continuous illumination. This study gives new insights into the plasticity mechanisms of teleost fish under constant illumination.

Evolved transcriptional responses and their trade-offs after long-term adaptation of Bemisia tabaci to a marginally-suitable host

Although generalist insect herbivores can migrate and rapidly adapt to a broad range of host plants, they can face significant difficulties when accidentally migrating to novel and marginally-suitable hosts. What happens, both in performance and gene expression regulation, if these marginally-suitable hosts must be used for multiple generations before migration to a suitable host can take place, largely remains unknown. In this study, we established multigenerational colonies of the whitefly Bemisia tabaci, a generalist phloem-feeding species, adapted to a marginally-suitable host (habanero pepper) or an optimal host (cotton). We used reciprocal host tests to estimate the differences in performance of the populations on both hosts under optimal (30 oC) and mild-stressful (24 oC) temperature conditions, and documented the associated transcriptomic changes. The habanero pepper-adapted population greatly improved its performance on habanero pepper but did not reach its performance level on cotton, the original host. It also showed reduced performance on cotton, relative to the non-adapted population, and an antagonistic effect of the lower-temperature stressor. The transcriptomic data revealed that most of the expression changes, associated with long-term adaptation to habanero pepper, can be categorized as "evolved" with no initial plastic response. Three molecular functions dominated: enhanced formation of cuticle structural constituents, enhanced activity of oxidation-reduction processes involved in neutralization of phytotoxins and reduced production of proteins from the cathepsin B family. Taken together, these findings indicate that generalist insects can adapt to novel host plants by modifying the expression of a relatively small set of specific molecular functions.

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.

How does epigenetics influence the course of evolution?

Epigenetics is the study of changes in gene activity that can be transmitted through cell divisions but cannot be explained by changes in the DNA sequence. Epigenetic mechanisms are central to gene regulation, phenotypic plasticity, development and the preservation of genome integrity. Epigenetic mechanisms are often held to make a minor contribution to evolutionary change because epigenetic states are typically erased and reset at every generation, and are therefore, not heritable. Nonetheless, there is growing appreciation that epigenetic variation makes direct and indirect contributions to evolutionary processes. First, some epigenetic states are transmitted intergenerationally and affect the phenotype of offspring. Moreover, bona fide heritable 'epialleles' exist and are quite common in plants. Such epialleles could, therefore, be subject to natural selection in the same way as conventional DNA sequence-based alleles. Second, epigenetic variation enhances phenotypic plasticity and phenotypic variance and thus can modulate the effect of natural selection on sequence-based genetic variation. Third, given that phenotypic plasticity is central to the adaptability of organisms, epigenetic mechanisms that generate plasticity and acclimation are important to consider in evolutionary theory. Fourth, some genes are under selection to be 'imprinted' identifying the sex of the parent from which they were derived, leading to parent-of-origin-dependent gene expression and effects. These effects can generate hybrid disfunction and contribute to speciation. Finally, epigenetic processes, particularly DNA methylation, contribute directly to DNA sequence evolution, because they act as mutagens on the one hand and modulate genome stability on the other by keeping transposable elements in check. This article is part of the theme issue 'How does epigenetics influence the course of evolution?'

Gene coexpression networks reveal molecular interactions underlying cichlid jaw modularity

Background

The oral and pharyngeal jaw of cichlid fishes are a classic example of evolutionary modularity as their functional decoupling boosted trophic diversification and contributed to the success of cichlid adaptive radiations. Most studies until now have focused on the functional, morphological, or genetic aspects of cichlid jaw modularity. Here we extend this concept to include transcriptional modularity by sequencing whole transcriptomes of the two jaws and comparing their gene coexpression networks.

Results

We show that transcriptional decoupling of gene expression underlies the functional decoupling of cichlid oral and pharyngeal jaw apparatus and the two units are evolving independently in recently diverged cichlid species from Lake Tanganyika. Oral and pharyngeal jaw coexpression networks reflect the common origin of the jaw regulatory program as there is high preservation of gene coexpression modules between the two sets of jaws. However, there is substantial rewiring of genetic architecture within those modules. We define a global jaw coexpression network and highlight jaw-specific and species-specific modules within it. Furthermore, we annotate a comprehensive in silico gene regulatory network linking the Wnt and AHR signalling pathways to jaw morphogenesis and response to environmental cues, respectively. Components of these pathways are significantly differentially expressed between the oral and pharyngeal jaw apparatus.

Conclusion

This study describes the concerted expression of many genes in cichlid oral and pharyngeal jaw apparatus at the onset of the independent life of cichlid fishes. Our findings suggest that – on the basis of an ancestral gill arch network—transcriptional rewiring may have driven the modular evolution of the oral and pharyngeal jaws, highlighting the evolutionary significance of gene network reuse. The gene coexpression and in silico regulatory networks presented here are intended as resource for future studies on the genetics of vertebrate jaw morphogenesis and trophic adaptation.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-021-01787-9.

Evolution of regulatory networks associated with traits under selection in cichlids

BACKGROUND

Seminal studies of vertebrate protein evolution speculated that gene regulatory changes can drive anatomical innovations. However, very little is known about gene regulatory network (GRN) evolution associated with phenotypic effect across ecologically diverse species. Here we use a novel approach for comparative GRN analysis in vertebrate species to study GRN evolution in representative species of the most striking examples of adaptive radiations, the East African cichlids. We previously demonstrated how the explosive phenotypic diversification of East African cichlids can be attributed to diverse molecular mechanisms, including accelerated regulatory sequence evolution and gene expression divergence.

RESULTS

To investigate these mechanisms across species at a genome-wide scale, we develop a novel computational pipeline that predicts regulators for co-extant and ancestral co-expression modules along a phylogeny, and candidate regulatory regions associated with traits under selection in cichlids. As a case study, we apply our approach to a well-studied adaptive trait-the visual system-for which we report striking cases of network rewiring for visual opsin genes, identify discrete regulatory variants, and investigate their association with cichlid visual system evolution. In regulatory regions of visual opsin genes, in vitro assays confirm that transcription factor binding site mutations disrupt regulatory edges across species and segregate according to lake species phylogeny and ecology, suggesting GRN rewiring in radiating cichlids.

CONCLUSIONS

Our approach reveals numerous novel potential candidate regulators and regulatory regions across cichlid genomes, including some novel and some previously reported associations to known adaptive evolutionary traits.

Genetic assimilation and accommodation: Models and mechanisms

Genetic assimilation and genetic accommodation are mechanisms by which novel phenotypes are produced and become established in a population. Novel characters may be fixed and canalized so they are insensitive to environmental variation, or can be plastic and adaptively responsive to environmental variation. In this review we explore the various theories that have been proposed to explain the developmental origin and evolution of novel phenotypes and the mechanisms by which canalization and phenotypic plasticity evolve. These theories and models range from conceptual to mathematical and have taken different views of how genes and environment contribute to the development and evolution of the properties of phenotypes. We will argue that a deeper and more nuanced understanding of genetic accommodation requires a recognition that phenotypes are not static entities but are dynamic system properties with no fixed deterministic relationship between genotype and phenotype. We suggest a mechanistic systems-view of development that allows one to incorporate both genes and environment in a common model, and that enables both quantitative analysis and visualization of the evolution of canalization and phenotypic plasticity.

In a comfort zone and beyond-Ecological plasticity of key marine mediators

Copepods of the genus Calanus are the key components of zooplankton. Understanding their response to a changing climate is crucial to predict the functioning of future warmer high-latitude ecosystems. Although specific Calanus species are morphologically very similar, they have different life strategies and roles in ecosystems. In this study, C. finmarchicus and C. glacialis were thoroughly studied with regard to their plasticity in morphology and ecology both in their preferred original water mass (Atlantic vs. Arctic side of the Polar Front) and in suboptimal conditions (due to, e.g., temperature, turbidity, and competition in Hornsund fjord). Our observations show that "at the same place and time," both species can reach different sizes, take on different pigmentation, be in different states of population development, utilize different reproductive versus lipid accumulation strategies, and thrive on different foods. Size was proven to be a very mutable morphological trait, especially with regard to reduced length of C. glacialis. Both species exhibited pronounced red pigmentation when inhabiting their preferred water mass. In other domains, C. finmarchicus individuals tended to be paler than C. glacialis individuals. Gonad maturation and population development indicated mixed reproductive strategies, although a surprisingly similar population age structure of the two co-occurring species in the fjord was observed. Lipid accumulation was high and not species-specific, and its variability was due to diet differences of the populations. According to the stable isotope composition, both species had a more herbivorous diatom-based diet in their original water masses. While the diet of C. glacialis was rather consistent among the domains studied, C. finmarchicus exhibited much higher variability in its feeding history (based on lipid composition). Our results show that the plasticity of both Calanus species is indeed impressive and may be regulated differently, depending on whether they live in their "comfort zone" or beyond it.

SEGN: Inferring real-time gene networks mediating phenotypic plasticity

The capacity of an organism to alter its phenotype in response to environmental perturbations changes over developmental time and is a process determined by multiple genes that are co-expressed in intricate but organized networks. Characterizing the spatiotemporal change of such gene networks can offer insight into the genomic signatures underlying organismic adaptation, but it represents a major methodological challenge. Here, we integrate the holistic view of systems biology and the interactive notion of evolutionary game theory to reconstruct so-called systems evolutionary game networks (SEGN) that can autonomously detect, track, and visualize environment-induced gene networks along the time axis. The SEGN overcomes the limitations of traditional approaches by inferring context-specific networks, encapsulating bidirectional, signed, and weighted gene-gene interactions into fully informative networks, and monitoring the process of how networks topologically alter across environmental and developmental cues. Based on the design principle of SEGN, we perform a transcriptional plasticity study by culturing Euphrates poplar, a tree that can grow in the saline desert, in saline-free and saline-stress conditions. SEGN characterize previously unknown gene co-regulation that modulates the time trajectories of the trees' response to salt stress. As a marriage of multiple disciplines, SEGN shows its potential to interpret gene interdependence, predict how transcriptional co-regulation responds to various regimes, and provides a hint for exploring the mass, energetic, or signal basis that drives various types of gene interactions.

[Phenotypic plasticity: a brief introduction]

Phenotypic plasticity describes the ability of a given genotype to produce different phenotypes in response to distinct environmental conditions. It has major implications in agronomy, animal husbandry and medicine and is also thought to facilitate evolution. Phenotypic plasticity is widely observed in the wild. It is only relatively recently that the mechanisms involved in phenotypic plasticity have been analysed. Thanks to laboratory experiments we understand better how environmental conditions are involved in phenotypic variations. This article introduces major concepts from the phenotypic plasticity field, presents briefly mechanisms involved in phenotypic plasticity and discusses the links between phenotypic plasticity and evolution.

Mammary tissue transcriptomic analysis for construction of integrated regulatory networks involved in lactogenesis of Ovis aries

The mammary gland experiences vast changes between the onset of lactation and pregnancy. This remodeling involves different functions such as lactation that is controlled by innumerable regulators and various gene networks which are still not completely understood. MicroRNAs (miRNAs) are one of the important non-coding gene regulators which control an extensive range of biological processes. Thus, exploring miRNAs functions is important for solving gene regulation complexity. The main purpose in the present study is to identify the various gene regulative integrated networks involved in lactation progress in mammary gland. We analyzed ovine mammary tissue data sets which included expression profiles of mRNA (genes) and miRNAs related to six ewes in different days of lactation and nutritional treatments. We combined two different types of information: the network that is module inference by mRNAs (RNA-seq data), miRNAs and transcription factors (TFs) expression matrix and prediction of targets via computational methods. To discover the miRNAs regulatory function, 134 modules were predicted by using gene expression data and 14 TFs and 20 miRNAs were allocated to these predicted modules. By applying this integrated computation-based method, 38 miRNA-modules and 35 TF-module interactions were identified from ovine mammary tissue data during lactogenesis. A lot of these modules were involved in lipid and protein metabolism, as well as steroids and vitamin biosynthesis, which would play key roles in mammary tissue and lactation development. These results present new information about the regulatory procedures at the miRNAs and TF levels throughout lactation.

The power and limitations of gene expression pathway analyses toward predicting population response to environmental stressors

Rapid environmental changes impact the global distribution and abundance of species, highlighting the urgency to understand and predict how populations will respond. The analysis of differentially expressed genes has elucidated areas of the genome involved in adaptive divergence to past and present environmental change. Such studies however have been hampered by large numbers of differentially expressed genes and limited knowledge of how these genes work in conjunction with each other. Recent methods (broadly termed "pathway analyses") have emerged that aim to group genes that behave in a coordinated fashion to a factor of interest. These methods aid in functional annotation and uncovering biological pathways, thereby collapsing complex datasets into more manageable units, providing more nuanced understandings of both the organism-level effects of modified gene expression, and the targets of adaptive divergence. Here, we reanalyze a dataset that investigated temperature-induced changes in gene expression in marine-adapted and freshwater-adapted threespine stickleback (Gasterosteus aculeatus), using Weighted Gene Co-expression Network Analysis (WGCNA) with PANTHER Gene Ontology (GO)-Slim overrepresentation and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis. Six modules exhibited a conserved response and six a divergent response between marine and freshwater stickleback when acclimated to 7°C or 22°C. One divergent module showed freshwater-specific response to temperature, and the remaining divergent modules showed differences in height of reaction norms. PPARAa, a transcription factor that regulates fatty acid metabolism and has been implicated in adaptive divergence, was located in a module that had higher expression at 7°C and in freshwater stickleback. This updated methodology revealed patterns that were not found in the original publication. Although such methods hold promise toward predicting population response to environmental stressors, many limitations remain, particularly with regard to module expression representation, database resources, and cross-database integration.

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.

Layers of Cryptic Genetic Variation Underlie a Yeast Complex Trait

To better understand cryptic genetic variation, Lee et al. comprehensively map the genetic basis of a trait that is typically suppressed in a yeast cross. By determining how three different genetic perturbations give rise...

Cryptic genetic variation may be an important contributor to heritable traits, but its extent and regulation are not fully understood. Here, we investigate the cryptic genetic variation underlying a Saccharomyces cerevisiae colony phenotype that is typically suppressed in a cross of the laboratory strain BY4716 (BY) and a derivative of the clinical isolate 322134S (3S). To do this, we comprehensively dissect the trait’s genetic basis in the BYx3S cross in the presence of three different genetic perturbations that enable its expression. This allows us to detect and compare the specific loci that interact with each perturbation to produce the trait. In total, we identify 21 loci, all but one of which interact with just a subset of the perturbations. Beyond impacting which loci contribute to the trait, the genetic perturbations also alter the extent of additivity, epistasis, and genotype–environment interaction among the detected loci. Additionally, we show that the single locus interacting with all three perturbations corresponds to the coding region of the cell surface gene FLO11. While nearly all of the other remaining loci influence FLO11 transcription in cis or trans, the perturbations tend to interact with loci in different pathways and subpathways. Our work shows how layers of cryptic genetic variation can influence complex traits. Here, these layers mainly represent different regulatory inputs into the transcription of a single key gene.

Species-complex diversification and host-plant associations in Bemisia tabaci: A plant-defence, detoxification perspective revealed by RNA-Seq analyses

Insect–plant associations and their role in diversification are mostly studied in specialists. Here, we aimed to identify macroevolution patterns in the relationships between generalists and their host plants that have the potential to promote diversification. We focused on the Bemisia tabaci species complex containing more than 35 cryptic species. Mechanisms for explaining this impressive diversification have focused so far on allopatric forces that assume a common, broad, host range. We conducted a literature survey which indicated that species in the complex differ in their host range, with only few showing a truly broad one. We then selected six species, representing different phylogenetic groups and documented host ranges. We tested whether differences in the species expression profiles of detoxification genes are shaped more by their phylogenetic relationships or by their ability to successfully utilize multiple hosts, including novel ones. Performance assays divided the six species into two groups of three, one showing higher performance on various hosts than the other (the lower performance group). The same grouping pattern appeared when the species were clustered according to their expression profiles. Only species placed in the lower performance group showed a tendency to lower the expression of multiple genes. Taken together, these findings bring evidence for the existence of a common detoxification “machinery,” shared between species that can perform well on multiple hosts. We raise the possibility that this “machinery” might have played a passive role in the diversification of the complex, by allowing successful migration to new/novel environments, leading, in some cases, to fragmentation and speciation.

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.

bric à brac (bab), a central player in the gene regulatory network that mediates thermal plasticity of pigmentation in Drosophila melanogaster

Drosophila body pigmentation has emerged as a major Evo-Devo model. Using two Drosophila melanogaster lines, Dark and Pale, selected from a natural population, we analyse here the interaction between genetic variation and environmental factors to produce this complex trait. Indeed, pigmentation varies with genotype in natural populations and is sensitive to temperature during development. We demonstrate that the bric à brac (bab) genes, that are differentially expressed between the two lines and whose expression levels vary with temperature, participate in the pigmentation difference between the Dark and Pale lines. The two lines differ in a bab regulatory sequence, the dimorphic element (called here bDE). Both bDE alleles are temperature-sensitive, but the activity of the bDE allele from the Dark line is lower than that of the bDE allele from the Pale line. Our results suggest that this difference could partly be due to differential regulation by AbdB. bab has been previously reported to be a repressor of abdominal pigmentation. We show here that one of its targets in this process is the pigmentation gene tan (t), regulated via the tan abdominal enhancer (t_MSE). Furthermore, t expression is strongly modulated by temperature in the two lines. Thus, temperature sensitivity of t expression is at least partly a consequence of bab thermal transcriptional plasticity. We therefore propose that a gene regulatory network integrating both genetic variation and temperature sensitivity modulates female abdominal pigmentation. Interestingly, both bDE and t_MSE were previously shown to have been recurrently involved in abdominal pigmentation evolution in drosophilids. We propose that the environmental sensitivity of these enhancers has turned them into evolutionary hotspots.

Complex traits such as size or disease susceptibility are typically modulated by both genetic variation and environmental conditions. Model organisms such as fruit flies (Drosophila) are particularly appropriate to analyse the interactions between genetic variation and environmental factors during the development of complex phenotypes. Natural populations carry high genetic variation and can be grown in controlled conditions in the laboratory. Here, we use Drosophila melanogaster female abdominal pigmentation, which is both genetically variable and modulated by the environment (temperature) to dissect this kind of interaction. We show that the pigmentation difference between two inbred fly lines is caused by genetic variation in an enhancer of the bab locus, which encodes two transcription factors controlling abdominal pigmentation. Indeed, this enhancer drives differential expression between the two lines. Interestingly, this enhancer is sensitive to temperature in both lines. We show that the effect of bab on pigmentation is mediated by the pigmentation gene tan (t) that is repressed by bab. Thus, the previously reported temperature-sensitive expression of t is a direct consequence of bab transcriptional plasticity.

Towards a gene regulatory network shaping the fins of the Princess cichlid

Variation in fin shape and size contributes to the outstanding morphological diversity of teleost fishes, but the regulation of fin growth has not yet been studied extensively outside the zebrafish model. A previous gene expression study addressing the ornamental elongations of unpaired fins in the African cichlid fish Neolamprologus brichardi identified three genes (cx43, mmp9 and sema3d) with strong and consistent expression differences between short and elongated fin regions. Remarkably, the expression patterns of these genes were not consistent with inferences on their regulatory interactions in zebrafish. Here, we identify a gene expression network (GRN) comprising cx43, mmp9, and possibly also sema3d by a stepwise approach of identifying co-expression modules and predicting their upstream regulators. Among the transcription factors (TFs) predicted as potential upstream regulators of 11 co-expressed genes, six TFs (foxc1, foxp1, foxd3, myc, egr2, irf8) showed expression patterns consistent with their cooperative transcriptional regulation of the gene network. Some of these TFs have already been implicated in teleost fish fin regeneration and formation. We particularly discuss the potential function of foxd3 as driver of the network and its role in the unexpected gene expression correlations observed in N. brichardi.

Similarities in temperature-dependent gene expression plasticity across timescales in threespine stickleback (Gasterosteus aculeatus)

Phenotypic plasticity occurs at a variety of timescales, but little is known about the degree to which plastic responses at different timescales are associated with similar underlying molecular processes, which is critical for assessing the effects of plasticity on evolutionary trajectories. To address this issue, we identified differential gene expression in response to developmental temperature in the muscle transcriptome of adult threespine stickleback (Gasterosteus aculeatus) exposed to 12, 18 and 24°C until hatch and then held at 18°C for 9 months and compared these results to differential gene expression in response to adult thermal acclimation in stickleback developed at 18°C and then acclimated to 5 and 25°C as adults. Adult thermal acclimation affected the expression of 7,940 and 7,015 genes in response to cold and warm acclimation, respectively, and 4,851 of these genes responded in both treatments. In contrast, the expression of only 33 and 29 genes was affected by cold and warm development, respectively. The majority of the genes affected by developmental temperature were also affected by adult acclimation temperature. Many genes that were differentially expressed as a result of adult acclimation were associated with previously identified temperature-dependent effects on DNA methylation patterns, suggesting a role of epigenetic mechanisms in regulating gene expression plasticity during acclimation. Taken together, these results demonstrate similarities between the persistent effects of developmental plasticity on gene expression and the effects of adult thermal acclimation, emphasizing the potential for mechanistic links between plasticity acting at these different life stages.

Phenotypic plasticity, canalization, and the origins of novelty: Evidence and mechanisms from amphibians

A growing number of biologists have begun asking whether environmentally induced phenotypic change--'phenotypic plasticity'--precedes and facilitates the origin and canalization of novel, complex phenotypes. However, such 'plasticity-first evolution' (PFE) remains controversial. Here, we summarize the PFE hypothesis and describe how it can be evaluated in natural systems. We then review the evidence for PFE from amphibians (a group in which phenotypic plasticity is especially widespread) and describe how phenotypic plasticity might have facilitated macroevolutionary change. Finally, we discuss what is known about the proximate mechanisms of PFE in amphibians. We close with suggestions for future research. As we describe, amphibians offer some of the best support for plasticity's role in the origin of evolutionary novelties.

A gene expression study of ornamental fin shape in Neolamprologus brichardi, an African cichlid species

The diversity of fin morphology within and across fish taxa offers great, but still largely unexplored, opportunities to investigate the proximate mechanisms underlying fin shape variation. Relying on available genetic knowledge brought forth mainly by the comprehensive study of the zebrafish caudal fin, we explored candidate molecular mechanisms for the maintenance and formation of the conspicuously elongated filaments adorning the unpaired fins of the East African "princess cichlid" Neolamprologus brichardi. Via qPCR assays, we detected expression differences of candidate genes between elongated and short regions of intact and regenerating fins. The identified genes include skeletogenic and growth factors (igf2b, fgf3, bmp2 and bmp4), components of the WNT pathway (lef1, wnt5b and wnt10) and a regulatory network determining fin ray segment size and junction (cx43, esco2 and sema3d), as well as other genes with different roles (mmp9, msxb and pea3). Interestingly, some of these genes showed fin specific expression differences which are often neglected in studies of model fish that focus on the caudal fin. Moreover, while the observed expression patterns were generally consistent with zebrafish results, we also detected deviating expression correlations and gene functions.

Screening of miRNA profiles and construction of regulation networks in early and late lactation of dairy goat mammary glands

In recent years, studies related to the expression profiles of miRNAs in the dairy goat mammary gland were performed, but regulatory mechanisms in the physiological environment and the dynamic homeostasis of mammary gland development and lactation are not clear. In the present study, sequencing data analysis of early and late lactation uncovered a total of 1,487 unique miRNAs, including 45 novel miRNA candidates and 1,442 known and conserved miRNAs, of which 758 miRNAs were co-expressed and 378 differentially expressed with P < 0.05. Moreover, 76 non-redundant target genes were annotated in 347 GO consortiums, with 3,143 candidate target genes grouped into 33 pathways. Additionally, 18 predicted target genes of 214 miRNAs were directly annotated in mammary gland development and used to construct regulatory networks based on GO annotation and the KEGG pathway. The expression levels of seven known miRNAs and three novel miRNAs were examined using quantitative real-time PCR. The results showed that miRNAs might play important roles in early and late lactation during dairy goat mammary gland development, which will be helpful to obtain a better understanding of the genetic control of mammary gland lactation and development.

Amphibious fishes: evolution and phenotypic plasticity

Amphibious fishes spend part of their life in terrestrial habitats. The ability to tolerate life on land has evolved independently many times, with more than 200 extant species of amphibious fishes spanning 17 orders now reported. Many adaptations for life out of water have been described in the literature, and adaptive phenotypic plasticity may play an equally important role in promoting favourable matches between the terrestrial habitat and behavioural, physiological, biochemical and morphological characteristics. Amphibious fishes living at the interface of two very different environments must respond to issues relating to buoyancy/gravity, hydration/desiccation, low/high O2 availability, low/high CO2 accumulation and high/low NH3 solubility each time they traverse the air-water interface. Here, we review the literature for examples of plastic traits associated with the response to each of these challenges. Because there is evidence that phenotypic plasticity can facilitate the evolution of fixed traits in general, we summarize the types of investigations needed to more fully determine whether plasticity in extant amphibious fishes can provide indications of the strategies used during the evolution of terrestriality in tetrapods.

How plasticity, genetic assimilation and cryptic genetic variation may contribute to adaptive radiations

There is increasing evidence that phenotypic plasticity can promote population divergence by facilitating phenotypic diversification and, eventually, genetic divergence. When a 'plastic' population colonizes a new habitat, it has the possibility to occupy multiple niches by expressing several distinct phenotypes. These initially reflect the population's plastic range but may later become genetically fixed by selection via the process of 'genetic assimilation' (GA). Through this process multiple specialized sister lineages can arise that share a common plastic ancestor - the 'flexible stem'. Here, we review possible molecular mechanisms through which natural selection could fix an initially plastic trait during GA. These mechanisms could also explain how GA may contribute to cryptic genetic variation that can subsequently be coopted into other phenotypes or traits, but also lead to nonadaptive responses. We outline the predicted patterns of genetic and transcriptional divergence accompanying flexible stem radiations. The analysis of such patterns of (retained) adaptive and nonadaptive plastic responses within and across radiating lineages can inform on the state of ongoing GA. We conclude that, depending on the stability of the environment, the molecular architecture underlying plastic traits can facilitate diversification, followed by fixation and consolidation of an adaptive phenotype and degeneration of nonadaptive ones. Additionally, the process of GA may increase the cryptic genetic variation of populations, which on one hand may serve as substrate for evolution, but on another may be responsible for nonadaptive responses that consolidate local allopatry and thus reproductive isolation.

Characterizing phenotypic divergence using three-dimensional geometric morphometrics in four populations of threespine stickleback (Gasterosteus aculeatus; Pisces: Gasterosteidae) in Katmai National Park and Preserve, Alaska

The threespine stickleback (Gasterosteus aculeatus L., 1758) is a vertebrate model for the study of the relationship between phenotype and environment in facilitating rapid evolutionary change. Using four populations from a system of lakes in Katmai National Park and Preserve, Alaska, and microcomputed tomography and three-dimensional geometric morphometrics, we test the hypothesis that stickleback populations inhabiting freshwater environments display cranial phenotypes that are intermediate between the putative ancestral form and the low-plated freshwater populations that demonstrate substantial divergence toward new phenotypic optima. We further test the hypothesis that phenotypic covariance structure is disrupted in the context of such putatively recent adaptive events. We report significant phenotypic differences among all four populations that includes a component of sexual dimorphism. Furthermore, we show evidence of disrupted phenotypic covariance structure among these populations. Taken together, these findings indicate the importance of phenotypic quantification as a key step in elucidating both the ecological processes responsible for rapid adaptive radiations and the role of developmental mechanisms in biasing evolutionary change.

The importance of mechanisms for the evolution of cooperation

Studies aimed at explaining the evolution of phenotypic traits have often solely focused on fitness considerations, ignoring underlying mechanisms. In recent years, there has been an increasing call for integrating mechanistic perspectives in evolutionary considerations, but it is not clear whether and how mechanisms affect the course and outcome of evolution. To study this, we compare four mechanistic implementations of two well-studied models for the evolution of cooperation, the Iterated Prisoner's Dilemma (IPD) game and the Iterated Snowdrift (ISD) game. Behavioural strategies are either implemented by a 1 : 1 genotype-phenotype mapping or by a simple neural network. Moreover, we consider two different scenarios for the effect of mutations. The same set of strategies is feasible in all four implementations, but the probability that a given strategy arises owing to mutation is largely dependent on the behavioural and genetic architecture. Our individual-based simulations show that this has major implications for the evolutionary outcome. In the ISD, different evolutionarily stable strategies are predominant in the four implementations, while in the IPD each implementation creates a characteristic dynamical pattern. As a consequence, the evolved average level of cooperation is also strongly dependent on the underlying mechanism. We argue that our findings are of general relevance for the evolution of social behaviour, pleading for the integration of a mechanistic perspective in models of social evolution.

Genetic assimilation: a review of its potential proximate causes and evolutionary consequences

BACKGROUND

Most, if not all, organisms possess the ability to alter their phenotype in direct response to changes in their environment, a phenomenon known as phenotypic plasticity. Selection can break this environmental sensitivity, however, and cause a formerly environmentally induced trait to evolve to become fixed through a process called genetic assimilation. Essentially, genetic assimilation can be viewed as the evolution of environmental robustness in what was formerly an environmentally sensitive trait. Because genetic assimilation has long been suggested to play a key role in the origins of phenotypic novelty and possibly even new species, identifying and characterizing the proximate mechanisms that underlie genetic assimilation may advance our basic understanding of how novel traits and species evolve.

SCOPE

This review begins by discussing how the evolution of phenotypic plasticity, followed by genetic assimilation, might promote the origins of new traits and possibly fuel speciation and adaptive radiation. The evidence implicating genetic assimilation in evolutionary innovation and diversification is then briefly considered. Next, the potential causes of phenotypic plasticity generally and genetic assimilation specifically are examined at the genetic, molecular and physiological levels and approaches that can improve our understanding of these mechanisms are described. The review concludes by outlining major challenges for future work.

CONCLUSIONS

Identifying and characterizing the proximate mechanisms involved in phenotypic plasticity and genetic assimilation promises to help advance our basic understanding of evolutionary innovation and diversification.

Multi-locus Genotypes Underlying Temperature Sensitivity in a Mutationally Induced Trait

Determining how genetic variation alters the expression of heritable phenotypes across conditions is important for agriculture, evolution, and medicine. Central to this problem is the concept of genotype-by-environment interaction (or 'GxE'), which occurs when segregating genetic variation causes individuals to show different phenotypic responses to the environment. While many studies have sought to identify individual loci that contribute to GxE, obtaining a deeper understanding of this phenomenon may require defining how sets of loci collectively alter the relationship between genotype, environment, and phenotype. Here, we identify combinations of alleles at seven loci that control how a mutationally induced colony phenotype is expressed across a range of temperatures (21, 30, and 37 °C) in a panel of yeast recombinants. We show that five predominant multi-locus genotypes involving the detected loci result in trait expression with varying degrees of temperature sensitivity. By comparing these genotypes and their patterns of trait expression across temperatures, we demonstrate that the involved alleles contribute to temperature sensitivity in different ways. While alleles of the transcription factor MSS11 specify the potential temperatures at which the trait can occur, alleles at the other loci modify temperature sensitivity within the range established by MSS11 in a genetic background- and/or temperature-dependent manner. Our results not only represent one of the first characterizations of GxE at the resolution of multi-locus genotypes, but also provide an example of the different roles that genetic variants can play in altering trait expression across conditions.

Genetic structure is associated with phenotypic divergence in floral traits and reproductive investment in a high-altitude orchid from the Iron Quadrangle, southeastern Brazil

Knowledge of the role of Neotropical montane landscapes in shaping genetic connectivity and local adaptation is essential for understanding the evolutionary processes that have shaped the extraordinary species diversity in these regions. In the present study, we examined the landscape genetics, estimated genetic diversity, and explored genetic relationships with morphological variability and reproductive strategies in seven natural populations of Cattleya liliputana (Orchidaceae). Nuclear microsatellite markers were used for genetic analyses. Spatial Bayesian clustering and population-based analyses revealed significant genetic structuring and high genetic diversity (He = 0.733 ± 0.03). Strong differentiation was found between populations over short spatial scales (FST = 0.138, p < 0.001), reflecting the landscape discontinuity and isolation. Monmonier´s maximum difference algorithm, Bayesian analysis on STRUCTURE and principal component analysis identified one major genetic discontinuity between populations. Divergent genetic groups showed phenotypic divergence in flower traits and reproductive strategies. Increased sexual reproductive effort was associated with rock outcrop type and may be a response to adverse conditions for growth and vegetative reproduction. Here we discuss the effect of restricted gene flow, local adaptation and phenotypic plasticity as drivers of population differentiation in Neotropical montane rock outcrops.