Genomic reaction norms: using integrative biology to understand molecular mechanisms of phenotypic plasticity

Genomic determinants, architecture, and constraints in drought-related traits in Corymbia calophylla

BACKGROUND

Drought adaptation is critical to many tree species persisting under climate change, however our knowledge of the genetic basis for trees to adapt to drought is limited. This knowledge gap impedes our fundamental understanding of drought response and application to forest production and conservation. To improve our understanding of the genomic determinants, architecture, and trait constraints, we assembled a reference genome and detected ~ 6.5 M variants in 432 phenotyped individuals for the foundational tree Corymbia calophylla.

RESULTS

We found 273 genomic variants determining traits with moderate heritability (h2SNP = 0.26-0.64). Significant variants were predominantly in gene regulatory elements distributed among several haplotype blocks across all chromosomes. Furthermore, traits were constrained by frequent epistatic and pleiotropic interactions.

CONCLUSIONS

Our results on the genetic basis for drought traits in Corymbia calophylla have several implications for the ability to adapt to climate change: (1) drought related traits are controlled by complex genomic architectures with large haplotypes, epistatic, and pleiotropic interactions; (2) the most significant variants determining drought related traits occurred in regulatory regions; and (3) models incorporating epistatic interactions increase trait predictions. Our findings indicate that despite moderate heritability drought traits are likely constrained by complex genomic architecture potentially limiting trees response to climate change.

Label-Free Quantitative Analysis of Pig Liver Proteome after Hepatitis E Virus Infection

Hepatitis E represents an emerging zoonotic disease caused by the Hepatitis E virus (HEV), for which the main route of transmission is foodborne. In particular, infection in humans has been associated with the consumption of contaminated undercooked meat of pig origin. The aim of this study was to apply comparative proteomics to determine if porcine liver protein profiles could be used to distinguish between pigs seropositive and seronegative for HEV. Preliminarily, an ELISA was used to evaluate the presence of anti-HEV antibodies in the blood serum of 136 animals sent to slaughter. Among the analyzed samples, a seroprevalence of 72.8% was estimated, and it was also possible to identify 10 animals, 5 positive and 5 negative, coming from the same farm. This condition created the basis for the quantitative proteomics comparison between homogeneous animals, in which only the contact with HEV should represent the discriminating factor. The analysis of the proteome in all samples of liver exudate led to the identification of 554 proteins differentially expressed between the two experimental groups, with 293 proteins having greater abundance in positive samples and 261 more represented in negative exudates. The pathway enrichment analysis allowed us to highlight the effect of the interaction between HEV and the host biological system in inducing the potential enrichment of 69 pathways. Among these, carbon metabolism stands out with the involvement of 41 proteins, which were subjected to interactomic analysis. This approach allowed us to focus our attention on three enzymes involved in glycolysis: glucose-6-phosphate isomerase (GPI), glyceraldehyde-3-phosphate dehydrogenase (GAPDH), and fructose-bisphosphate aldolase A (ALDOA). It therefore appears that infection with HEV induced a strengthening of the process, which involves the breakdown of glucose to obtain energy and carbon residues useful for the virus's survival. In conclusion, the label-free LC-MS/MS approach showed effectiveness in highlighting the main differences induced on the porcine liver proteome by the interaction with HEV, providing crucial information in identifying a viral signature on the host metabolism.

Sexual conflict, heterochrony and tissue specificity as evolutionary problems of adaptive plasticity in development

Differential gene expression represents a fundamental cause and manifestation of phenotypic plasticity. Adaptive phenotypic plasticity in gene expression as a trait evolves when alleles that mediate gene regulation serve to increase organismal fitness by improving the alignment of variation in gene expression with variation in circumstances. Among the diverse circumstances that a gene encounters are distinct cell types, developmental stages and sexes, as well as an organism's extrinsic ecological environments. Consequently, adaptive phenotypic plasticity provides a common framework to consider diverse evolutionary problems by considering the shared implications of alleles that produce context-dependent gene expression. From this perspective, adaptive plasticity represents an evolutionary resolution to conflicts of interest that arise from any negatively pleiotropic effects of expression of a gene across ontogeny, among tissues, between the sexes, or across extrinsic environments. This view highlights shared properties within the general relation of fitness, trait expression and context that may nonetheless differ substantively in the grain of selection within and among generations to influence the likelihood of adaptive plasticity as an evolutionary response. Research programmes that historically have focused on these separate issues may use the insights from one another by recognizing their shared dependence on context-dependent gene regulatory evolution.

Responses to predation pressure involve similar sets of genes in two divergent species of Daphnia

Species that are not closely related can express similar inducible traits, but molecular mechanisms underlying the observed responses are often unknown, nor is it known if these mechanisms are shared between such species. Here, we compared transcriptional profiles of two Daphnia species (D. mitsukuri and D. sinensis) from different subgenera, at both juvenile and adult developmental stages. Both species were exposed to the same predation threat (fish kairomones), and both showed similar induced morphological changes (reduced body length). At the early developmental stage, response to predation risk resulted in similar changes in expression levels of 23 orthologues in both species. These orthologues, involved in 107 GO categories, changed in the same direction in both species (over- or underexpressed), in comparison to non-exposed controls. Several of these orthologues were associated with DNA replication, structural constituents of cuticle or innate immune response. In both species, the differentially expressed (DE) genes on average had higher ω (dN /dS ) values than non-DE genes, suggesting that these genes had experienced greater positive selection or lower purifying selection than non-DE genes. Overall, our results suggest that similar suites of genes, responding in similar ways to predation pressure, have been retained in Daphnia for many millions of years.

Relationships between phenotypic plasticity and epigenetic variation in two Caribbean Acropora corals

The plastic ability for a range of phenotypes to be exhibited by the same genotype allows organisms to respond to environmental variation and may modulate fitness in novel environments. Differing capacities for phenotypic plasticity within a population, apparent as genotype by environment interactions (GxE), can therefore have both ecological and evolutionary implications. Epigenetic gene regulation alters gene function in response to environmental cues without changes to the underlying genetic sequence and likely mediates phenotypic variation. DNA methylation is currently the most well described epigenetic mechanism and is related to transcriptional homeostasis in invertebrates. However, evidence quantitatively linking variation in DNA methylation with that of phenotype is lacking in some taxa, including reef-building corals. In this study, spatial and seasonal environmental variation in Bonaire, Caribbean Netherlands was utilized to assess relationships between physiology and DNA methylation profiles within genetic clones across different genotypes of Acropora cervicornis and A. palmata corals. The physiology of both species was highly influenced by environmental variation compared to the effect of genotype. GxE effects on phenotype were only apparent in A. cervicornis. DNA methylation in both species differed between genotypes and seasons and epigenetic variation was significantly related to coral physiological metrics. Furthermore, plastic shifts in physiology across seasons were significantly positively correlated with shifts in DNA methylation profiles in both species. These results highlight the dynamic influence of environmental conditions and genetic constraints on the physiology of two important Caribbean coral species. Additionally, this study provides quantitative support for the role of epigenetic DNA methylation in mediating phenotypic plasticity in invertebrates.

The Epigenetics of Animal Personality

Animal personality, consistent individual differences in behaviour, is an important concept for understanding how individuals vary in how they cope with environmental challenges. In order to understand the evolutionary significance of animal personality, it is crucial to understand the underlying regulatory mechanisms. Epigenetic marks such as DNA methylation are hypothesised to play a major role in explaining variation in phenotypic changes in response to environmental alterations. Several characteristics of DNA methylation also align well with the concept of animal personality. In this review paper, we summarise the current literature on the role that molecular epigenetic mechanisms may have in explaining personality variation. We elaborate on the potential for epigenetic mechanisms to explain behavioural variation, behavioural development and temporal consistency in behaviour. We then suggest future routes for this emerging field and point to potential pitfalls that may be encountered. We conclude that a more inclusive approach is needed for studying the epigenetics of animal personality and that epigenetic mechanisms cannot be studied without considering the genetic background.

The lack of genetic variation underlying thermal transcriptomic plasticity suggests limited adaptability of the Northern shrimp, Pandalus borealis

Introduction

Genetic variation underlies the populations’ potential to adapt to and persist in a changing environment, while phenotypic plasticity can play a key role in buffering the negative impacts of such change at the individual level.

Methods

We investigated the role of genetic variation in the thermal response of the northern shrimp Pandalus borealis, an ectotherm species distributed in the Arctic and North Atlantic Oceans. More specifically, we estimated the proportion transcriptomic responses explained by genetic variance of female shrimp from three origins after 30 days of exposure to three temperature treatments.

Results

We characterized the P. borealis transcriptome (170,377 transcripts, of which 27.48% were functionally annotated) and then detected a total of 1,607 and 907 differentially expressed transcripts between temperatures and origins, respectively. Shrimp from different origins displayed high but similar level of transcriptomic plasticity in response to elevated temperatures. Differences in transcript expression among origins were not correlated to population genetic differentiation or diversity but to environmental conditions at origin during sampling.

Discussion

The lack of genetic variation explaining thermal plasticity suggests limited adaptability in this species’ response to future environmental changes. These results together with higher mortality observed at the highest temperature indicate that the thermal niche of P. borealis will likely be restricted to higher latitudes in the future. This prediction concurs with current decreases in abundance observed at the southern edge of this species geographical distribution, as it is for other cold-adapted crustaceans.

Integrated omic techniques and their genomic features for invasive weeds

Many emerging invasive weeds display rapid adaptation against different stressful environments compared to their natives. Rapid adaptation and dispersal habits helped invasive populations have strong diversity within the population compared to their natives. Advances in molecular marker techniques may lead to an in-depth understanding of the genetic diversity of invasive weeds. The use of molecular techniques is rapidly growing, and their implications in invasive weed studies are considered powerful tools for genome purposes. Here, we review different approach used multi-omics by invasive weed studies to understand the functional structural and genomic changes in these species under different environmental fluctuations, particularly, to check the accessibility of advance-sequencing techniques used by researchers in genome sequence projects. In this review-based study, we also examine the importance and efficiency of different molecular techniques in identifying and characterizing different genes, associated markers, proteins, metabolites, and key metabolic pathways in invasive and native weeds. Use of these techniques could help weed scientists to further reduce the knowledge gaps in understanding invasive weeds traits. Although these techniques can provide robust insights about the molecular functioning, employing a single omics platform can rarely elucidate the gene-level regulation and the associated real-time expression of weedy traits due to the complex and overlapping nature of biological interactions. We conclude that different multi-omic techniques will provide long-term benefits in launching new genome projects to enhance the understanding of invasive weeds' invasion process.

A Label-Free Quantitative Analysis for the Search of Proteomic Differences between Goat Breeds

The intensification and standardization of livestock farming are causing a decline in the number of animal breeds in many species, such as the goat. The availability of more studies on the potentiality of goat breeds could raise awareness of their importance, conservation and productive possibilities. Label-free quantitative analysis was applied in this study to investigate the proteomic differences between the autochthon Teramana and Saanen goats that could be useful for defining peculiar features of these breeds. A total of 2093 proteins were characterized in the muscle exudate proteome of the Teramana and Saanen breeds. A total of 41 proteins clearly separated the two breeds. Eukaryotic initiation factor proteins and aldehyde-dehydrogenase 7 family-member A1 were up-regulated in the autochthon breed and associated with its resilience, whereas catalase was down-regulated and associated with lower muscular mass. This study is the most detailed report of goat muscle proteome. Several differentially regulated proteins between the two breeds were identified, providing insights into functional pathways that define this organism and its biology.

Disentangling the developmental origins of a novel phenotype: enhancement versus reversal of environmentally induced gene expression

Increasing evidence suggests that many novel traits might have originated via plasticity-led evolution (PLE). Yet, little is known of the developmental processes that underpin PLE, especially in its early stages. One such process is 'phenotypic accommodation', which occurs when, in response to a change in the environment, an organism experiences adjustments across variable parts of its phenotype that improve its fitness. Here, we asked if environmentally induced changes in gene expression are enhanced or reversed during phenotypic accommodation of a novel, complex phenotype in spadefoot toad tadpoles (Spea multiplicata). More genes than expected were affected by both the environment and phenotypic accommodation in the liver and brain. However, although phenotypic accommodation primarily reversed environmentally induced changes in gene expression in liver tissue, it enhanced these changes in brain tissue. Thus, depending on the tissue, phenotypic accommodation may either minimize functional disruption via reversal of gene expression patterns or promote novelty via enhancement of existing expression patterns. Our study thereby provides insights into the developmental origins of a novel phenotype and the incipient stages of PLE.

Transcriptional memory of gene expression across generations participates in transgenerational plasticity of field pennycress in response to cadmium stress

Transgenerational plasticity (TGP) occurs when maternal environments influence the expression of traits in offspring, and in some cases may increase fitness of offspring and have evolutionary significance. However, little is known about the extent of maternal environment influence on gene expression of offspring, and its relationship with trait variations across generations. In this study, we examined TGP in the traits and gene expression of field pennycress (Thlaspi arvense) in response to cadmium (Cd) stress. In the first generation, along with the increase of soil Cd concentration, the total biomass, individual height, and number of seeds significantly decreased, whereas time to flowering, superoxide dismutase (SOD) activity, and content of reduced glutathione significantly increased. Among these traits, only SOD activity showed a significant effect of TGP; the offspring of Cd-treated individuals maintained high SOD activity in the absence of Cd stress. According to the results of RNA sequencing and bioinformatic analysis, 10,028 transcripts were identified as Cd-responsive genes. Among them, only 401 were identified as transcriptional memory genes (TMGs) that maintained the same expression pattern under normal conditions in the second generation as in Cd-treated parents in the first generation. These genes mainly participated in Cd tolerance-related processes such as response to oxidative stress, cell wall biogenesis, and the abscisic acid signaling pathways. The results of weighted correlation network analysis showed that modules correlated with SOD activity recruited more TMGs than modules correlated with other traits. The SOD-coding gene CSD2 was found in one of the modules correlated with SOD activity. Furthermore, several TMGs co-expressed with CSD2 were hub genes that were highly connected to other nodes and critical to the network's topology; therefore, recruitment of TMGs in offspring was potentially related to TGP. These findings indicated that, across generations, transcriptional memory of gene expression played an important role in TGP. Moreover, these results provided new insights into the trait evolution processes mediated by phenotypic plasticity.

Vertebrate Phenological Plasticity: from Molecular Mechanisms to Ecological and Evolutionary Implications

Seasonal variation in the availability of essential resources is one of the most important drivers of natural selection on the phasing and duration of annually recurring life-cycle events. Shifts in seasonal timing are among the most commonly reported responses to climate change and the capacity of organisms to adjust their timing, either through phenotypic plasticity or evolution, is a critical component of resilience. Despite growing interest in documenting and forecasting the impacts of climate change on phenology, our ability to predict how individuals, populations, and species might alter their seasonal timing in response to their changing environments is constrained by limited knowledge regarding the cues animals use to adjust timing, the endogenous genetic and molecular mechanisms that transduce cues into neural and endocrine signals, and the inherent capacity of animals to alter their timing and phasing within annual cycles. Further, the fitness consequences of phenological responses are often due to biotic interactions within and across trophic levels, rather than being simple outcomes of responses to changes in the abiotic environment. Here, we review the current state of knowledge regarding the mechanisms that control seasonal timing in vertebrates, as well as the ecological and evolutionary consequences of individual, population, and species-level variation in phenological responsiveness. Understanding the causes and consequences of climate-driven phenological shifts requires combining ecological, evolutionary, and mechanistic approaches at individual, populational, and community scales. Thus, to make progress in forecasting phenological responses and demographic consequences, we need to further develop interdisciplinary networks focused on climate change science.

Phenotypic and transcriptional response of Daphnia pulicaria to the combined effects of temperature and predation

Daphnia, an ecologically important zooplankton species in lakes, shows both genetic adaptation and phenotypic plasticity in response to temperature and fish predation, but little is known about the molecular basis of these responses and their potential interactions. We performed a factorial experiment exposing laboratory-propagated Daphnia pulicaria clones from two lakes in the Sierra Nevada mountains of California to normal or high temperature (15°C or 25°C) in the presence or absence of fish kairomones, then measured changes in life history and gene expression. Exposure to kairomones increased upper thermal tolerance limits for physiological activity in both clones. Cloned individuals matured at a younger age in response to higher temperature and kairomones, while size at maturity, fecundity and population intrinsic growth were only affected by temperature. At the molecular level, both clones expressed more genes differently in response to temperature than predation, but specific genes involved in metabolic, cellular, and genetic processes responded differently between the two clones. Although gene expression differed more between clones from different lakes than experimental treatments, similar phenotypic responses to predation risk and warming arose from these clone-specific patterns. Our results suggest that phenotypic plasticity responses to temperature and kairomones interact synergistically, with exposure to fish predators increasing the tolerance of Daphnia pulicaria to stressful temperatures, and that similar phenotypic responses to temperature and predator cues can be produced by divergent patterns of gene regulation.

Extent and complexity of RNA processing in honey bee queen and worker caste development

The distinct honeybee (Apis mellifera) worker and queen castes have become a model for the study of genomic mechanisms of phenotypic plasticity. Here we performed a nanopore-based direct RNA sequencing with exceptionally long reads to compare the mRNA transcripts between queen and workers at three points during their larval development. We found thousands of significantly differentially expressed transcript isoforms (DEIs) between queen and worker larvae. These DEIs were formatted by a flexible splicing system. We showed that poly(A) tails participated in this caste differentiation by negatively regulating the expression of DEIs. Hundreds of isoforms uniquely expressed in either queens or workers during their larval development, and isoforms were expressed at different points in queen and worker larval development demonstrating a dynamic relationship between isoform expression and developmental mechanisms. These findings show the full complexity of RNA processing and transcript expression in honey bee phenotypic plasticity.

•

Honeybee caste differentiation has a complexity of RNA processing

•

Isoforms differentially express between queens and workers during larval development

•

Isoforms are formatted by a flexible alternative splicing system

•

Poly(A) tails are negatively correlated with isoform expression

Genomic reaction norms inform predictions of plastic and adaptive responses to climate change

Genomic reaction norms represent the range of gene expression phenotypes (usually mRNA transcript levels) expressed by a genotype along an environmental gradient. Reaction norms derived from common‐garden experiments are powerful approaches for disentangling plastic and adaptive responses to environmental change in natural populations. By treating gene expression as a phenotype in itself, genomic reaction norms represent invaluable tools for exploring causal mechanisms underlying organismal responses to climate change across multiple levels of biodiversity.

Our goal is to provide the context, framework and motivation for applying genomic reaction norms to study the responses of natural populations to climate change.

Here, we describe the utility of integrating genomics with common‐garden‐gradient experiments under a reaction norm analytical framework to answer fundamental questions about phenotypic plasticity, local adaptation, their interaction (i.e. genetic variation in plasticity) and future adaptive potential.

An experimental and analytical framework for constructing and analysing genomic reaction norms is presented within the context of polygenic climate change responses of structured populations with gene flow. Intended for a broad eco‐evo readership, we first briefly review adaptation with gene flow and the importance of understanding the genomic basis and spatial scale of adaptation for conservation and management of structured populations under anthropogenic change. Then, within a high‐dimensional reaction norm framework, we illustrate how to distinguish plastic, differentially expressed (difference in reaction norm intercepts) and differentially plastic (difference in reaction norm slopes) genes, highlighting the areas of opportunity for applying these concepts.

We conclude by discussing how genomic reaction norms can be incorporated into a holistic framework to understand the eco‐evolutionary dynamics of climate change responses from molecules to ecosystems. We aim to inspire researchers to integrate gene expression measurements into common‐garden experimental designs to investigate the genomics of climate change responses as sequencing costs become increasingly accessible.

Strong parallel differential gene expression induced by hatchery rearing weakly associated with methylation signals in adult Coho Salmon (O. kisutch)

Human activities and resource exploitation led to a massive decline of wild salmonid populations, consequently, numerous conservation programs have been developed to supplement wild populations. However, many studies documented reduced fitness of hatchery-born relative to wild fish. Here, by using both RNA sequencing and Whole Genome Bisulfite Sequencing of hatchery and wild-born adult Coho salmon (Oncorhynchus kisutch) originating from two previously studied river systems, we show that early-life hatchery-rearing environment-induced significant and parallel gene expression differentiation is maintained until Coho come back to their natal river for reproduction. A total of 3,643 genes differentially expressed and 859 coexpressed genes were downregulated in parallel in hatchery-born fish from both rivers relative to their wild congeners. Among those genes, 26 displayed a significant relationship between gene expression and the median gene body methylation and 669 single CpGs displayed a significant correlation between methylation level and the associated gene expression. The link between methylation and gene expression was weak suggesting that DNA methylation is not the only player in mediating hatchery-related expression differences. Yet, significant gene expression differentiation was observed despite 18 months spent in a common environment (i.e., the sea). Finally, the differentiation is observed in parallel in two different river systems, highlighting the fact that early-life environment may account for at least some of the reduced fitness of the hatchery salmon in the wild. These results illustrate the relevance and importance of considering both epigenome and transcriptome to evaluate the costs and benefits of large-scale supplementation programs.

Brain transcriptome analysis reveals gene expression differences associated with dispersal behaviour between range-front and range-core populations of invasive cane toads in Australia

Understanding the mechanisms allowing invasive species to adapt to novel environments is a challenge in invasion biology. Many invaders demonstrate rapid evolution of behavioural traits involved in range expansion such as locomotor activity, exploration and risk‐taking. However, the molecular mechanisms that underpin these changes are poorly understood. In 86 years, invasive cane toads (Rhinella marina) in Australia have drastically expanded their geographic range westward from coastal Queensland to Western Australia. During their range expansion, toads have undergone extensive phenotypic changes, particularly in behaviours that enhance the toads’ dispersal ability. Common‐garden experiments have shown that some changes in behavioural traits related to dispersal are heritable. At the molecular level, it is currently unknown whether these changes in dispersal‐related behaviour are underlain by small or large differences in gene expression, nor is known the biological function of genes showing differential expression. Here, we used RNA‐seq to gain a better understanding of the molecular mechanisms underlying dispersal‐related behavioural changes. We compared the brain transcriptomes of toads from the Hawai'ian source population, as well as three distinct populations from across the Australian invasive range. We found markedly different gene expression profiles between the source population and Australian toads. By contrast, toads from across the Australian invasive range had very similar transcriptomic profiles. Yet, key genes with functions putatively related to dispersal behaviour showed differential expression between populations located at each end of the invasive range. These genes could play an important role in the behavioural changes characteristic of range expansion in Australian cane toads.

The role of epigenetic regulation in adaptive phenotypic plasticity of plants

In recent decades, knowledge about the role of epigenetic regulation of gene expression in plant responses to external stimuli and in adaptation of plants to adverse environmental fluctuations have extended significantly. DNA methylation is considered as the main molecular mechanism that provides genomic information and contributes to the understanding of the molecular basis of phenotypic variations based on epigenetic modifications. Unfortunately, the vast majority of research in this area has been performed on the model species Arabidopsis thaliana. The development of the methylation-sensitive amplified polymorphism (MSAP) method has made it possible to implement the large-scale detection of DNA methylation alterations in wild non-model and agricultural plants with large and highly repetitive genomes in natural and manipulated habitats. The article presents current information on DNA methylation in species of natural communities and crops and its importance in plant development and adaptive phenotypic plasticity, along with brief reviews of current ideas about adaptive phenotypic plasticity and epigenetic regulation of gene expression. The great potential of further studies of the epigenetic role in phenotypic plasticity of a wide range of non-model species in natural populations and agrocenoses for understanding the molecular mechanisms of plant existence in the changing environment in onto- and phylogeny, directly related to the key tasks of forecasting the effects of global warming and crop selection, is emphasized. Specific taxa of the Ukrainian flora, which, in authors’ opinion, are promising and interesting for this type of research, are recommended.

A neglected conceptual problem regarding phenotypic plasticity's role in adaptive evolution: The importance of genetic covariance and social drive

There is tantalizing evidence that phenotypic plasticity can buffer novel, adaptive genetic variants long enough to permit their evolutionary spread, and this process is often invoked in explanations for rapid adaptive evolution. However, the strength and generality of evidence for it is controversial. We identify a conceptual problem affecting this debate: recombination, segregation, and independent assortment are expected to quickly sever associations between genes controlling novel adaptations and genes contributing to trait plasticity that facilitates the novel adaptations by reducing their indirect fitness costs. To make clearer predictions about this role of plasticity in facilitating genetic adaptation, we describe a testable genetic mechanism that resolves the problem: genetic covariance between new adaptive variants and trait plasticity that facilitates their persistence within populations. We identify genetic architectures that might lead to such a covariance, including genetic coupling via physical linkage and pleiotropy, and illustrate the consequences for adaptation rates using numerical simulations. Such genetic covariances may also arise from the social environment, and we suggest the indirect genetic effects that result could further accentuate the process of adaptation. We call the latter mechanism of adaptation social drive, and identify methods to test it. We suggest that genetic coupling of plasticity and adaptations could promote unusually rapid ‘runaway’ evolution of novel adaptations. The resultant dynamics could facilitate evolutionary rescue, adaptive radiations, the origin of novelties, and other commonly studied processes.

Epigenomic modifications induced by hatchery rearing persist in germ line cells of adult salmon after their oceanic migration

Human activities induce direct or indirect selection pressure on natural population and may ultimately affect population's integrity. While numerous conservation programs aimed to minimize human-induced genomic variation, human-induced environmental variation may generate epigenomic variation potentially affecting fitness through phenotypic modifications. Major questions remain pertaining to how much epigenomic variation arises from environmental heterogeneity, whether this variation can persist throughout life, and whether it can be transmitted across generations. We performed whole genome bisulfite sequencing (WGBS) on the sperm of genetically indistinguishable hatchery and wild-born migrating adults of Coho salmon (Oncorhynchus kisutch) from two geographically distant rivers at different epigenome scales. Our results showed that coupling WGBS with fine-scale analyses (local and chromosomal) allows the detection of parallel early-life hatchery-induced epimarks that differentiate wild from hatchery-reared salmon. Four chromosomes and 183 differentially methylated regions (DMRs) displayed a significant signal of methylation differentiation between hatchery and wild-born Coho salmon. Moreover, those early-life epimarks persisted in germ line cells despite about 1.5 year spent in the ocean following release from hatchery, opening the possibility for transgenerational inheritance. Our results strengthen the hypothesis that epigenomic modifications environmentally induced during early-life development persist in germ cells of adults until reproduction, which could potentially impact their fitness.

Ganoderma (Ganodermataceae, Basidiomycota) Species from the Greater Mekong Subregion

The cosmopolitan fungal genus Ganoderma is an important pathogen on arboreal plant hosts, particularly in tropical and temperate regions. It has long been used as a traditional medicine because of its medicinal properties and chemical constituents. In this study, Ganoderma collections were made in the Greater Mekong Subregion (GMS), encompassing tropical parts of Laos, Myanmar, Thailand, Vietnam, and temperate areas in Yunnan Province, China. The specimens used in this study are described based on micro-macro-characteristics and phylogenetic analysis of combined ITS, LSU, TEF1α, and RPB2 sequence data. In this comprehensive study, we report 22 Ganoderma species from the GMS, namely, G. adspersum, G. applanatum, G. australe, G. calidophilum, G. ellipsoideum, G. flexipes, G. gibbosum, G. heohnelianum, G. hochiminhense, G. leucocontextum, G. lucidum, G. multiplicatum, G. multipileum, G. myanmarense, G. orbiforme, G. philippii, G. resinaceum, G. sichuanense, G. sinense, G. subresinosum, G. williamsianum, and G. tsugae. Some of these species were reported in more than one country within the GMS. Of these 22 species, 12 were collected from Yunnan Province, China; three were collected from Laos; three species, two new records, and one new species were collected from Myanmar; 15 species and four new records were collected from Thailand, and one new species was collected from Vietnam. Comprehensive descriptions, color photographs of macro- and micro-characteristics, the distribution of Ganoderma within the GMS, as well as a phylogenetic tree showing the placement of all reported Ganoderma from the GMS are provided.

Epigenetic Regulation and Environmental Sex Determination in Cichlid Fishes

Studying environmental sex determination (ESD) in cichlids provides a phylogenetic and comparative approach to understand the evolution of the underlying mechanisms, their impact on the evolution of the overlying systems, and the neuroethology of life history strategies. Natural selection normally favors parents who invest equally in the development of male and female offspring, but evolution may favor deviations from this 50:50 ratio when environmental conditions produce an advantage for doing so. Many species of cichlids demonstrate ESD in response to water chemistry (temperature, pH, and oxygen concentration). The relative strengths of and the exact interactions between these factors vary between congeners, demonstrating genetic variation in sensitivity. The presence of sizable proportions of the less common sex towards the environmental extremes in most species strongly suggests the presence of some genetic sex-determining loci acting in parallel with the ESD factors. Sex determination and differentiation in these species does not seem to result in the organization of a final and irreversible sexual fate, so much as a life-long ongoing battle between competing male- and female-determining genetic and hormonal networks governed by epigenetic factors. We discuss what is and is not known about the epigenetic mechanism behind the differentiation of both gonads and sex differences in the brain. Beyond the well-studied tilapia species, the 2 best-studied dwarf cichlid systems showing ESD are the South American genus Apistogramma and the West African genus Pelvicachromis. Both species demonstrate male morphs with alternative reproductive tactics. We discuss the further neuroethology opportunities such systems provide to the study of epigenetics of alternative life history strategies and other behavioral variation.

Polyphenism of a Novel Trait Integrated Rapidly Evolving Genes into Ancestrally Plastic Networks

Developmental polyphenism, the ability to switch between phenotypes in response to environmental variation, involves the alternating activation of environmentally sensitive genes. Consequently, to understand how a polyphenic response evolves requires a comparative analysis of the components that make up environmentally sensitive networks. Here, we inferred coexpression networks for a morphological polyphenism, the feeding-structure dimorphism of the nematode Pristionchus pacificus. In this species, individuals produce alternative forms of a novel trait—moveable teeth, which in one morph enable predatory feeding—in response to environmental cues. To identify the origins of polyphenism network components, we independently inferred coexpression modules for more conserved transcriptional responses, including in an ancestrally nonpolyphenic nematode species. Further, through genome-wide analyses of these components across the nematode family (Diplogastridae) in which the polyphenism arose, we reconstructed how network components have changed. To achieve this, we assembled and resolved the phylogenetic context for five genomes of species representing the breadth of Diplogastridae and a hypothesized outgroup. We found that gene networks instructing alternative forms arose from ancestral plastic responses to environment, specifically starvation-induced metabolism and the formation of a conserved diapause (dauer) stage. Moreover, loci from rapidly evolving gene families were integrated into these networks with higher connectivity than throughout the rest of the P. pacificus transcriptome. In summary, we show that the modular regulatory outputs of a polyphenic response evolved through the integration of conserved plastic responses into networks with genes of high evolutionary turnover.

Thyroid hormones regulate the formation and environmental plasticity of white bars in clownfishes

Determining how plasticity of developmental traits responds to environmental conditions is a challenge that must combine evolutionary sciences, ecology, and developmental biology. During metamorphosis, fish alter their morphology and color pattern according to environmental cues. We observed that juvenile clownfish (Amphiprion percula) modulate the developmental timing of their adult white bar formation during metamorphosis depending on the sea anemone species in which they are recruited. We observed an earlier formation of white bars when clownfish developed with Stichodactyla gigantea (Sg) than with Heteractis magnifica (Hm). As these bars, composed of iridophores, form during metamorphosis, we hypothesized that timing of their development may be thyroid hormone (TH) dependent. We treated clownfish larvae with TH and found that white bars developed earlier than in control fish. We further observed higher TH levels, associated with rapid white bar formation, in juveniles recruited in Sg than in Hm, explaining the faster white bar formation. Transcriptomic analysis of Sg recruits revealed higher expression of duox, a dual oxidase implicated in TH production as compared to Hm recruits. Finally, we showed that duox is an essential regulator of iridophore pattern timing in zebrafish. Taken together, our results suggest that TH controls the timing of adult color pattern formation and that shifts in duox expression and TH levels are associated with ecological differences resulting in divergent ontogenetic trajectories in color pattern development.

Cross-species transcriptomics uncovers genes underlying genetic accommodation of developmental plasticity in spadefoot toads

That hardcoded genomes can manifest as plastic phenotypes responding to environmental perturbations is a fascinating feature of living organisms. How such developmental plasticity is regulated at the molecular level is beginning to be uncovered aided by the development of -omic techniques. Here, we compare the transcriptome-wide responses of two species of spadefoot toads with differing capacity for developmental acceleration of their larvae in the face of a shared environmental risk: pond drying. By comparing gene expression profiles over time and performing cross-species network analyses, we identified orthologues and functional gene pathways whose environmental sensitivity in expression have diverged between species. Genes related to lipid, cholesterol and steroid biosynthesis and metabolism make up most of a module of genes environmentally responsive in one species, but canalized in the other. The evolutionary changes in the regulation of the genes identified through these analyses may have been key in the genetic accommodation of developmental plasticity in this system.

Fine-scale empirical data on niche divergence and homeolog expression patterns in an allopolyploid and its diploid progenitor species

Polyploidization is pervasive in plants, but little is known about the niche divergence of wild allopolyploids (species that harbor polyploid genomes originating from different diploid species) relative to their diploid progenitor species and the gene expression patterns that may underlie such ecological divergence. We conducted a fine-scale empirical study on habitat and gene expression of an allopolyploid and its diploid progenitors. We quantified soil properties and light availability of habitats of an allotetraploid Cardamine flexuosa and its diploid progenitors Cardamine amara and Cardamine hirsuta in two seasons. We analyzed expression patterns of genes and homeologs (homeologous gene copies in allopolyploids) using RNA sequencing. We detected niche divergence between the allopolyploid and its diploid progenitors along water availability gradient at a fine scale: the diploids in opposite extremes and the allopolyploid in a broader range between diploids, with limited overlap with diploids at both ends. Most of the genes whose homeolog expression ratio changed among habitats in C. flexuosa varied spatially and temporally. These findings provide empirical evidence for niche divergence between an allopolyploid and its diploid progenitor species at a fine scale and suggest that divergent expression patterns of homeologs in an allopolyploid may underlie its persistence in diverse habitats.

Genotype by environment interaction for fat and protein yields via reaction norms in Holstein cattle of southern Brazil

Our objective was to evaluate the genetic merit of Holstein cattle population in southern Brazil in response to variations in the regional temperature by analyzing the genotype by environment interaction using reaction norms. Fat yield (FY) and protein yield (PY) data of 67 360 primiparous cows were obtained from the database of the Paraná Holstein Breeders Association, Brazil (APCBRH). The regional average annual temperature was used as the environmental variable. A random regression model was adopted applying mixed models with Restricted Maximum Likelihood (REML) algorithm using WOMBAT software. The genetic merit of the 15 most representative bulls, depending on the temperature gradient, was evaluated. Heritability ranged from 0.21 to 0.27 for FY and from 0.14 to 0.20 for PY. The genetic correlation observed among the environmental gradients proved to be higher than 0.80 for both traits. Slight reranking of bulls for both traits was detected, demonstrating that non-relevant genotype by environment interaction for FY and PY were observed. Consequently, no inclusion of the temperature effect in the model of genetic evaluation in southern Brazilian Holstein breed is required.

Plasticity of Leaf Traits of Juglans regia L. f. luodianense Liu et Xu Seedlings Under Different Light Conditions in Karst Habitats

This study examined the effects of light intensity on the plasticity of the leaves of Juglans regia f. luodianense seedlings in karst habitat and how they respond to changes in light intensity. The light intensity of 1-year-old seedlings of J. regia f. luodianense in different niches in a karst area was set as 100% (bare land), 75% (forest margin), 50% (forest gap), and 25% (under forest) of natural light. The material harvested after four months was compared to analyze the differences in various morphological characteristics, biomass allocation, and physiological characteristics of the leaves of seedlings of J. regia f. luodianense, and a comprehensive evaluation of the plasticity indexes was conducted. The results showed that under moderate (50%) full light intensity, the leaf area, specific leaf area, leaf biomass, and chlorophyll content increased, and improved photosynthesis and promoted the accumulation of free proline content and peroxidase (POD) activity. The accumulation of malondialdehyde was also the lowest in this treatment, indicating that the plants had the strongest adaptability under this light intensity. Moreover, under high (75%) full light intensity, the above functional characteristics of plants showed good performance. Under low (25%) full light intensity, plants also had higher specific leaf area, leaf biomass, and photosynthetic parameters. However, under full light, the cell membrane permeability decreased, the chlorophyll accumulation was the lowest, and the photosynthetic index was seriously inhibited. Our results showed that the plasticity of morphological characters was greater than that of biomass allocation and physiological characters; POD activity and stomatal conductance were the highest, followed by leaf area and chlorophyll b, whereas the plasticity of palisade tissue/sponge tissue thickness and lower-epidermis thickness were the lowest. In summary, there are evident differences in the sensitivity and regulation mechanisms of morphological characteristics, biomass allocation, and physiological indices of the seedling leaves of J. regia f. luodianense in response to light intensity. During the stage of seedling establishment, only the plants in the bare ground under full light can be induced to show obvious inhibition of phenotypic traits. In contrast, the plants in the forest margins and gaps and under the forest habitats under light intensity can regulate their own characteristics to maintain their growth and development. The wide light range and strong plasticity of the species might be two of the important reasons for its existence in a highly heterogeneous karst habitat.

Gene expression remodelling and immune response during adaptive divergence in an African cichlid fish

Variation in gene expression contributes to ecological speciation by facilitating population persistence in novel environments. Likewise, immune responses can be of relevance in speciation driven by adaptation to different environments. Previous studies examining gene expression differences between recently diverged ecotypes have often relied on only one pair of populations, targeted the expression of only a subset of genes or used wild-caught individuals. Here, we investigated the contribution of habitat-specific parasites and symbionts and the underlying immunological abilities of ecotype hosts to adaptive divergence in lake-river population pairs of the cichlid fish Astatotilapia burtoni. To shed light on the role of phenotypic plasticity in adaptive divergence, we compared parasite and microbiota communities, immune response, and gene expression patterns of fish from natural habitats and a lake-like pond set-up. In all investigated population pairs, lake fish were more heavily parasitized than river fish, in terms of both parasite taxon composition and infection abundance. The innate immune response in the wild was higher in lake than in river populations and was elevated in a river population exposed to lake parasites in the pond set-up. Environmental differences between lake and river habitat and their distinct parasite communities have shaped differential gene expression, involving genes functioning in osmoregulation and immune response. Most changes in gene expression between lake and river samples in the wild and in the pond set-up were based on a plastic response. Finally, gene expression and bacterial communities of wild-caught individuals and individuals acclimatized to lake-like pond conditions showed shifts underlying adaptive phenotypic plasticity.

Host behaviour alteration by its parasite: from brain gene expression to functional test

Many parasites with complex life cycles modify their intermediate hosts' behaviour, presumably to increase transmission to their final host. The threespine stickleback (Gasterosteus aculeatus) is an intermediate host in the cestode Schistocephalus solidus life cycle, which ends in an avian host, and shows increased risky behaviours when infected. We studied brain gene expression profiles of sticklebacks infected with S. solidus to determine the proximal causes of these behavioural alterations. We show that infected fish have altered expression levels in genes involved in the inositol pathway. We thus tested the functional implication of this pathway and successfully rescued normal behaviours in infected sticklebacks using lithium exposure. We also show that exposed but uninfected fish have a distinct gene expression profile from both infected fish and control individuals, allowing us to separate gene activity related to parasite exposure from consequences of a successful infection. Finally, we find that selective serotonin reuptake inhibitor-treated sticklebacks and infected fish do not have similarly altered gene expression, despite their comparable behaviours, suggesting that the serotonin pathway is probably not the main driver of phenotypic changes in infected sticklebacks. Taken together, our results allow us to predict that if S. solidus directly manipulates its host, it could target the inositol pathway.

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.

Genotype by environment interaction for somatic cell score in Holstein cattle of southern Brazil via reaction norms

Objective

The objective of this study was to evaluate the genetic behavior of a population of Holstein cattle in response to the variation of environmental temperature by analyzing the effects of genotype by environment interaction (GEI) through reaction norms for the somatic cell score (SCS).

Methods

Data was collected for 67,206 primiparous cows from the database of the Paraná Holstein Breeders Association in Brazil, with the aim of evaluating the temperature effect, considered as an environmental variable, distinguished under six gradients, with the variation range found being 17°C to 19.5°C, over the region. A reaction norm model was adopted utilizing the fourth order under the Legendre polynomials, using the mixed models of analysis by the restricted maximum likelihood method by the WOMBAT software. Additionally, the genetic behavior of the 15 most representative bulls was assessed, in response to the changes in the temperature gradient.

Results

A mean score of 2.66 and a heritability variation from 0.17 to 0.23 was found in the regional temperature increase. The correlation between the environmental gradients proved to be higher than 0.80. Distinctive genetic behaviors were observed according to the increase in regional temperature, with an observed increase of up to 0.258 in the breeding values of some animals, as well as a reduction in the breeding of up to 0.793, with occasional reclassifications being observed as the temperature increased.

Conclusion

Non-relevant GEI for SCS were observed in Holstein cattle herds of southern Brazil. Thus, the inclusion of the temperature effect in the model of genetic evaluation of SCS for the southern Brazilian Holstein breed is not required.

Gene Expression and Photophysiological Changes in Pocillopora acuta Coral Holobiont Following Heat Stress and Recovery

The ability of corals to withstand changes in their surroundings is a critical survival mechanism for coping with environmental stress. While many studies have examined responses of the coral holobiont to stressful conditions, its capacity to reverse responses and recover when the stressor is removed is not well-understood. In this study, we investigated among-colony responses of Pocillopora acuta from two sites with differing distance to the mainland (Kusu (closer to the mainland) and Raffles Lighthouse (further from the mainland)) to heat stress through differential expression analysis of target genes and quantification of photophysiological metrics. We then examined how these attributes were regulated after the stressor was removed to assess the recovery potential of P. acuta. The fragments that were subjected to heat stress (2 °C above ambient levels) generally exhibited significant reduction in their endosymbiont densities, but the extent of recovery following stress removal varied depending on natal site and colony. There were minimal changes in chl a concentration and maximum quantum yield (Fv/Fm, the proportion of variable fluorescence (Fv) to maximum fluorescence (Fm)) in heat-stressed corals, suggesting that the algal endosymbionts’ Photosystem II was not severely compromised. Significant changes in gene expression levels of selected genes of interest (GOI) were observed following heat exposure and stress removal among sites and colonies, including Actin, calcium/calmodulin-dependent protein kinase type IV (Camk4), kinesin-like protein (KIF9), and small heat shock protein 16.1 (Hsp16.1). The most responsive GOIs were Actin, a major component of the cytoskeleton, and the adaptive immune-related Camk4 which both showed significant reduction following heat exposure and subsequent upregulation during the recovery phase. Our findings clearly demonstrate specific responses of P. acuta in both photophysiological attributes and gene expression levels, suggesting differential capacity of P. acuta corals to tolerate heat stress depending on the colony, so that certain colonies may be more resilient than others.

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.

Brain transcriptomics of agonistic behaviour in the weakly electric fish Gymnotus omarorum, a wild teleost model of non-breeding aggression

Differences in social status are often mediated by agonistic encounters between competitors. Robust literature has examined social status-dependent brain gene expression profiles across vertebrates, yet social status and reproductive state are often confounded. It has therefore been challenging to identify the neuromolecular mechanisms underlying social status independent of reproductive state. Weakly electric fish, Gymnotus omarorum, display territorial aggression and social dominance independent of reproductive state. We use wild-derived G. omarorum males to conduct a transcriptomic analysis of non-breeding social dominance relationships. After allowing paired rivals to establish a dominance hierarchy, we profiled the transcriptomes of brain sections containing the preoptic area (region involved in regulating aggressive behaviour) in dominant and subordinate individuals. We identified 16 differentially expressed genes (FDR < 0.05) and numerous genes that co-varied with behavioural traits. We also compared our results with previous reports of differential gene expression in other teleost species. Overall, our study establishes G. omarorum as a powerful model system for understanding the neuromolecular bases of social status independent of reproductive state.

Forebrain Transcriptional Response to Transient Changes in Circulating Androgens in a Cichlid Fish

It has been hypothesized that androgens respond to the social interactions as a way to adjust the behavior of individuals to the challenges of the social environment in an adaptive manner. Therefore, it is expected that transient changes in circulating androgen levels within physiological scope should impact the state of the brain network that regulates social behavior, which should translate into adaptive behavioral changes. Here, we examined the effect that a transient peak in androgen circulating levels, which mimics socially driven changes in androgen levels, has on the forebrain state, which harbors most nuclei of the social decision-making network. For this purpose, we successfully induced transient changes in circulating androgen levels in an African cichlid fish (Mozambique tilapia, Oreochromis mossambicus) commonly used as a model in behavioral neuroendocrinology by injecting 11-ketotestosterone or testosterone, and compared the forebrain transcriptome of these individuals to control fish injected with vehicle. Forebrain samples were collected 30 min and 60 min after injection and analyzed using RNAseq. Our results showed that a transient peak in 11-ketotestosterone drives more accentuated changes in forebrain transcriptome than testosterone, and that transcriptomic impact was greater at the 30 min than at the 60 min post-androgen administration. Several genes involved in the regulation of translation, steroid metabolism, ion channel membrane receptors, and genes involved in epigenetic mechanisms were differentially expressed after 11-ketotestosterone or testosterone injection. In summary, this study identified specific candidate genes that may regulate socially driven changes in behavioral flexibility mediated by androgens.

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.

Phenotypic plasticity as a mechanism of cave colonization and adaptation

A widely accepted model for the evolution of cave animals posits colonization by surface ancestors followed by the acquisition of adaptations over many generations. However, the speed of cave adaptation in some species suggests mechanisms operating over shorter timescales. To address these mechanisms, we used Astyanax mexicanus, a teleost with ancestral surface morphs (surface fish, SF) and derived cave morphs (cavefish, CF). We exposed SF to completely dark conditions and identified numerous altered traits at both the gene expression and phenotypic levels. Remarkably, most of these alterations mimicked CF phenotypes. Our results indicate that many cave-related traits can appear within a single generation by phenotypic plasticity. In the next generation, plasticity can be further refined. The initial plastic responses are random in adaptive outcome but may determine the subsequent course of evolution. Our study suggests that phenotypic plasticity contributes to the rapid evolution of cave-related traits in A. mexicanus.

The Mexican tetra is a fish that has two forms: a surface-dwelling form, which has eyes and silvery grey appearance, and a cave-dwelling form, which is blind and has lost its pigmentation. Recent studies have shown that the cave-dwelling form evolved rapidly within the last 200,000 years from an ancestor that lived at the surface. The recent evolution of the cave-dwelling form of the tetra poses an interesting evolutionary question: how did the surface-dwelling ancestor of the tetra quickly adapt to the new and challenging environment found in the caves?

‘Phenotypic plasticity’ is a phenomenon through which a single set of genes can produce different observable traits depending on the environment. An example of phenotypic plasticity occurs in response to diet: in animals, poor diets can lead to an increase in the size of the digestive organs and to the animals eating more. To see if surface-dwelling tetras can quickly adapt to cave environments through phenotypic plasticity, Bilandžija et al. have exposed these fish to complete darkness (the major feature of the cave environment) for two years. After spending up to two years in the dark, these fish were compared to normal surface-dwelling and cave-dwelling tetras.

Results revealed that surface-dwelling tetras raised in the dark exhibited traits associated with cave-dwelling tetras. These traits included changes in the activity of many genes involved in diverse processes, resistance to starvation, metabolism, and levels of hormones and molecules involved in neural signaling, which could lead to changes in behavior. However, the fish also exhibited traits, including an increase in the cells responsible for pigmentation, that would have no obvious benefit in the darkness. Even though the changes observed require no genetic mutations, they can help or hinder the fish’s survival once they occur, possibly determining subsequent evolution. Thus, a trait beneficial for surviving in the dark that appears simply through phenotypic plasticity may eventually be selected for and genetic mutations that encode it more reliably may appear too.

These results shed light on how species may quickly adapt to new environments without accumulating genetic mutations, which can take hundreds of thousands of years. They also may help to explain how colonizer species succeed in challenging environments. The principles described by Bilandžija et al. can be applied to different organisms adapting to new environments, and may help understand the role of phenotypic plasticity in evolution.

Evolutionary genomics can improve prediction of species' responses to climate change

Global climate change (GCC) increasingly threatens biodiversity through the loss of species, and the transformation of entire ecosystems. Many species are challenged by the pace of GCC because they might not be able to respond fast enough to changing biotic and abiotic conditions. Species can respond either by shifting their range, or by persisting in their local habitat. If populations persist, they can tolerate climatic changes through phenotypic plasticity, or genetically adapt to changing conditions depending on their genetic variability and census population size to allow for de novo mutations. Otherwise, populations will experience demographic collapses and species may go extinct. Current approaches to predicting species responses to GCC begin to combine ecological and evolutionary information for species distribution modelling. Including an evolutionary dimension will substantially improve species distribution projections which have not accounted for key processes such as dispersal, adaptive genetic change, demography, or species interactions. However, eco-evolutionary models require new data and methods for the estimation of a species' adaptive potential, which have so far only been available for a small number of model species. To represent global biodiversity, we need to devise large-scale data collection strategies to define the ecology and evolutionary potential of a broad range of species, especially of keystone species of ecosystems. We also need standardized and replicable modelling approaches that integrate these new data to account for eco-evolutionary processes when predicting the impact of GCC on species' survival. Here, we discuss different genomic approaches that can be used to investigate and predict species responses to GCC. This can serve as guidance for researchers looking for the appropriate experimental setup for their particular system. We furthermore highlight future directions for moving forward in the field and allocating available resources more effectively, to implement mitigation measures before species go extinct and ecosystems lose important functions.

Divergent Allometric Trajectories in Gene Expression and Coexpression Produce Species Differences in Sympatrically Speciating Midas Cichlid Fish

The mechanisms of speciation without geographic isolation (i.e., sympatric speciation) remain debated. This is due in part to the fact that the genomic landscape that could promote or hinder species divergence in the presence of gene flow is still largely unknown. However, intensive research is now centered on understanding the genetic architecture of adaptive traits associated with this process as well as how gene expression might affect these traits. Here, using RNA-Seq data, we investigated gene expression of sympatrically speciating benthic and limnetic Neotropical cichlid fishes at two developmental stages. First, we identified groups of coexpressed genes (modules) at each stage. Although there are a few large and well-preserved modules, most of the other modules are not preserved across life stages. Second, we show that later in development more and larger coexpression modules are associated with divergence between benthic and limnetic fish compared with the earlier life stage. This divergence between benthic and limnetic fish in coexpression mirrors divergence in overall expression between benthic and limnetic fish, which is more pronounced later in life. Our results reveal that already at 1-day posthatch benthic and limnetic fish diverge in (co)expression, and that this divergence becomes more substantial when fish are free-swimming but still unlikely to have divergent swimming and feeding habits. More importantly, our study describes how the coexpression of several genes through development, as opposed to individual genes, is associated with benthic–limnetic species differences, and how two morphogenetic trajectories diverge as fish grow older.

Tissue-specific expression profiles and positive selection analysis in the tree swallow (Tachycineta bicolor) using a de novo transcriptome assembly

Tree swallows (Tachycineta bicolor) are one of the most commonly studied wild birds in North America. They have advanced numerous research areas, including life history, physiology, and organismal responses to global change; however, transcriptomic resources are scarce. To further advance the utility of this system for biologists across disciplines, we generated a transcriptome for the tree swallow using six tissues (brain, blood, ovary, spleen, liver, and muscle) collected from breeding females. We de novo assembled 207,739 transcripts, which we aligned to 14,717 high confidence protein-coding genes. We then characterized each tissue with regard to its unique genes and processes and applied this transcriptome to two fundamental questions in evolutionary biology and endocrinology. First, we analyzed 3,015 single-copy orthologs and identified 46 genes under positive selection in the tree swallow lineage, including those with putative links to adaptations in this species. Second, we analyzed tissue-specific expression patterns of genes involved in sex steroidogenesis and processing. Enzymes capable of synthesizing these behaviorally relevant hormones were largely limited to the ovary, whereas steroid binding genes were found in nearly all other tissues, highlighting the potential for local regulation of sex steroid-mediated traits. These analyses provide new insights into potential sources of phenotypic variation in a free-living female bird and advance our understanding of fundamental questions in evolutionary and organismal biology.

Altering social cue perception impacts honey bee aggression with minimal impacts on aggression-related brain gene expression

Gene expression changes resulting from social interactions may give rise to long term behavioral change, or simply reflect the activity of neural circuitry associated with behavioral expression. In honey bees, social cues broadly modulate aggressive behavior and brain gene expression. Previous studies suggest that expression changes are limited to contexts in which social cues give rise to stable, relatively long-term changes in behavior. Here we use a traditional beekeeping approach that inhibits aggression, smoke exposure, to deprive individuals of aggression-inducing olfactory cues and evaluate whether behavioral changes occur in absence of expression variation in a set of four biomarker genes (drat, cyp6g1/2, GB53860, inos) associated with aggression in previous studies. We also evaluate two markers of a brain hypoxic response (hif1α, hsf) to determine whether smoke induces molecular changes at all. We find that bees with blocked sensory perception as a result of smoke exposure show a strong, temporary inhibition of aggression relative to bees allowed to perceive normal social cues. However, blocking sensory perception had minimal impacts on aggression-relevant gene expression, althought it did induce a hypoxic molecular response in the brain. Results suggest that certain genes differentiate social cue-induced changes in aggression from long-term modulation of this phenotype.

Physiological status of silver carp (Hypophthalmichthys molitrix) in the Illinois River: An assessment of fish at the leading edge of the invasion front

Silver carp (Hypophthalmichthys molitrix) are invasive to North America, and their range has expanded within the Mississippi River Basin, seemingly unchecked, since their introduction in the late 1970s, with the exception of the upper reaches of the Illinois River. With the imminent threat of their movement into the Great Lakes, the goal of the present study was to assess whether differences in the physiological status between silver carp at the leading edge of their invasion front and core population sites could explain their lack of expansion upstream toward Lake Michigan over the past decade. A transcriptomic approach using RNA sequencing and analysis of plasma variables were used to quantify differences among fish at the leading edge and two downstream core population sites. Leading-edge fish exhibited upregulation of genes associated with xenobiotic defense (e.g., ATP-binding cassette C1 [abcc1], abcc2, abcc6), decreased cell integrity (i.e., macroautophagy and apoptosis; autophagy-related protein 9A [atg9a], caspase 3b [casp3b]), and cholesterol metabolism (e.g., abca1, apolipoprotein A1 [apoa1], sterol O-acyltransferase [soat1]) and downregulation of genes associated with DNA repair (e.g., tumor suppressor p53-binding protein 1 [tp53bp1]) compared to core population sites. Transcriptomic profiles of leading-edge fish were consistent with fish inhabiting a polluted environment and suggest that poorer water quality conditions upstream of the leading edge may represent a non-permanent barrier to silver carp range expansion. The present study provides potential molecular targets for monitoring the physiological status of silver carp over time and in response to future improvements in water quality upstream of their leading edge.

Genomics of Developmental Plasticity in Animals

Developmental plasticity refers to the property by which the same genotype produces distinct phenotypes depending on the environmental conditions under which development takes place. By allowing organisms to produce phenotypes adjusted to the conditions that adults will experience, developmental plasticity can provide the means to cope with environmental heterogeneity. Developmental plasticity can be adaptive and its evolution can be shaped by natural selection. It has also been suggested that developmental plasticity can facilitate adaptation and promote diversification. Here, we summarize current knowledge on the evolution of plasticity and on the impact of plasticity on adaptive evolution, and we identify recent advances and important open questions about the genomics of developmental plasticity in animals. We give special attention to studies using transcriptomics to identify genes whose expression changes across developmental environments and studies using genetic mapping to identify loci that contribute to variation in plasticity and can fuel its evolution.