Genetic accommodation via modified endocrine signalling explains phenotypic divergence among spadefoot toad species

The role of phenotypic plasticity and corticosterone in coping with pond drying conditions in yellow-bellied toad (Bombina variegata, Linnaeus 1758) tadpoles

Amphibian larvae inhabiting temporary ponds often exhibit the capacity to accelerate development and undergo metamorphosis in challenging conditions like desiccation. However, not all species exhibit this ability, the yellow-bellied toad (Bombina variegata) is one such example. The underlying mechanisms behind the inability to accelerate development under desiccation remain largely unexplored. The hypothalamic-pituitary-interrenal (HPI) axis and corticosterone (CORT), which act synergistically with thyroid hormone, are thought to facilitate metamorphosis in response to desiccation stress. In this study, we aimed to investigate whether modification in the HPI axis, particularly CORT levels, contributes to the absence of adaptive plasticity in B. variegata under desiccation stress. The study design included four treatments: high water level, high water level with exogenous CORT, low water level, and low water level with metyrapone (a CORT synthesis inhibitor). The main objective was to evaluate the effects of these treatments on whole-body corticosterone levels, life history, morphological traits, and oxidative stress parameters during the prometamorphic and metamorphic climax developmental stages. While low water level had no effect on total corticosterone levels, larval period, body condition index, and metamorphic body shape, it negatively affected metamorph size, mass, and growth rate. Our findings suggest that constant exposure to desiccation stress over generations may have led to modifications in the HPI axis activity in B. variegata, resulting in adaptation to changes in water level, evident through the absence of stress response. Consequently, CORT may not be a relevant stress indicator in desiccation conditions for this species.

Acclimation capacity to global warming of amphibians and freshwater fishes: Drivers, patterns, and data limitations

Amphibians and fishes play a central role in shaping the structure and function of freshwater environments. These organisms have a limited capacity to disperse across different habitats and the thermal buffer offered by freshwater systems is small. Understanding determinants and patterns of their physiological sensitivity across life history is, therefore, imperative to predicting the impacts of climate change in freshwater systems. Based on a systematic literature review including 345 experiments with 998 estimates on 96 amphibian (Anura/Caudata) and 93 freshwater fish species (Teleostei), we conducted a quantitative synthesis to explore phylogenetic, ontogenetic, and biogeographic (thermal adaptation) patterns in upper thermal tolerance (CTmax) and thermal acclimation capacity (acclimation response ratio, ARR) as well as the influence of the methodology used to assess these thermal traits using a conditional inference tree analysis. We found globally consistent patterns in CTmax and ARR, with phylogeny (taxa/order), experimental methodology, climatic origin, and life stage as significant determinants of thermal traits. The analysis demonstrated that CTmax does not primarily depend on the climatic origin but on experimental acclimation temperature and duration, and life stage. Higher acclimation temperatures and longer acclimation times led to higher CTmax values, whereby Anuran larvae revealed a higher CTmax than older life stages. The ARR of freshwater fishes was more than twice that of amphibians. Differences in ARR between life stages were not significant. In addition to phylogenetic differences, we found that ARR also depended on acclimation duration, ramping rate, and adaptation to local temperature variability. However, the amount of data on early life stages is too small, methodologically inconsistent, and phylogenetically unbalanced to identify potential life cycle bottlenecks in thermal traits. We, therefore, propose methods to improve the robustness and comparability of CTmax/ARR data across species and life stages, which is crucial for the conservation of freshwater biodiversity under climate change.

Got milkweed? Genetic assimilation as potential source for the evolution of nonmigratory monarch butterfly wing shape

Monarch butterflies (Danaus plexippus) are well studied for their annual long-distance migration from as far north as Canada to their overwintering grounds in Central Mexico. At the end of the cold season, monarchs start to repopulate North America through short-distance migration over the course of multiple generations. Interestingly, some populations in various tropical and subtropical islands do not migrate and exhibit heritable differences in wing shape and size, most likely an adaptation to island life. Less is known about forewing differences between long- and short-distance migrants in relation to island populations. Given their different migratory behaviors, we hypothesized that these differences would be reflected in wing morphology. To test this, we analyzed forewing shape and size of three different groups: nonmigratory, lesser migratory (migrate short-distances), and migratory (migrate long-distances) individuals. Significant differences in shape appear in all groups using geometric morphometrics. As variation found between migratory and lesser migrants has been shown to be caused by phenotypic plasticity, and lesser migrants develop intermediate forewing shapes between migratory and nonmigratory individuals, we suggest that genetic assimilation might be an important mechanism to explain the heritable variation found between migratory and nonmigratory populations. Additionally, our research confirms previous studies which show that forewing size is significantly smaller in nonmigratory populations when compared to both migratory phenotypes. Finally, we found sexual dimorphism in forewing shape in all three groups, but for size in nonmigratory populations only. This might have been caused by reduced constraints on forewing size in nonmigratory populations.

Larval plastic responses to warming and desiccation delay gonadal maturation in postmetamorphic spadefoot toads

Developmental plasticity allows organisms to adjust life-history traits to varying environmental conditions, which can have concomitant effects across life stages. Many amphibians are suitable model systems to study plasticity because their larvae can adjust growth and differentiation under fluctuating environments. It is unknown, however, whether somatic and gonadal differentiation are equally affected by environmentally induced plasticity or whether their decoupling alters gonadal maturation postmetamorphosis, which may affect fitness. We tested if developmental acceleration in response to warming and desiccation risk results in shifts in gonadal maturation during metamorphosis and postmetamorphic growth in western spadefoot toads (Pelobates cultripes). We found additive effects of increased temperature and desiccation risk on development and growth at metamorphosis, which largely constrained gonadal maturation in metamorphic and postmetamorphic individuals of both sexes. Furthermore, the conditions experienced by larvae incurred sex-specific carryover effects on the gonadal maturation of juveniles 5 months after metamorphosis. In females, high temperature delayed ovarian maturation regardless of the water level. In males, exposure to high temperature and high water levels slightly delayed the testes' maturation. These results highlight the relevance of larval plasticity in the gonadal maturation of species undergoing metamorphosis, which may have implications for population demographics and the evolution of life histories.

Promises and limits of an agency perspective in evolutionary developmental biology

An agent-based perspective in the study of complex systems is well established in diverse disciplines, yet is only beginning to be applied to evolutionary developmental biology. In this essay, we begin by defining agency and associated terminology formally. We then explore the assumptions and predictions of an agency perspective, apply these to select processes and key concept areas relevant to practitioners of evolutionary developmental biology, and consider the potential epistemic roles that an agency perspective might play in evo devo. Throughout, we discuss evidence supportive of agential dynamics in biological systems relevant to evo devo and explore where agency thinking may enrich the explanatory reach of research efforts in evolutionary developmental biology.

Early onset of urea synthesis and ammonia detoxification pathways in three terrestrially developing frogs

Frogs evolved terrestrial development multiple times, necessitating mechanisms to avoid ammonia toxicity at early stages. Urea synthesis from ammonia is a key adaptation that reduces water dependence after metamorphosis. We tested for early expression and plasticity of enzymatic mechanisms of ammonia detoxification in three terrestrial-breeding frogs: foam-nest-dwelling larvae of Leptodactylus fragilis (Lf) and arboreal embryos of Hyalinobatrachium fleischmanni (Hf) and Agalychnis callidryas (Ac). Activity of two ornithine-urea cycle (OUC) enzymes, arginase and CPSase, and levels of their products urea and CP in tissues were high in Lf regardless of nest hydration, but reduced in experimental low- vs. high-ammonia environments. High OUC activity in wet and dry nests, comparable to that under experimental high ammonia, suggests terrestrial Lf larvae maintain high capacity for urea excretion regardless of their immediate risk of ammonia toxicity. This may aid survival through unpredictably long waiting periods before rain enables their transition to water. Moderate levels of urea and CP were present in Hf and Ac tissues and enzymatic activities were lower than in Lf. In both species, embryos in drying clutches can hatch and enter the water early, behaviorally avoiding ammonia toxicity. Moreover, glutamine synthetase was active in early stages of all three species, condensing ammonia and glutamate to glutamine as another mechanism of detoxification. Enzyme activity appeared highest in Lf, although substrate and product levels were higher in Ac and Lf. Our results reveal that multiple biochemical mechanisms of ammonia detoxification occur in early life stages of anuran lineages that evolved terrestrial development.

Exploring water-borne corticosterone collection as a non-invasive tool in amphibian conservation physiology: benefits, limitations and future perspectives

Global change exposes wildlife to a variety of environmental stressors and is affecting biodiversity worldwide, with amphibian population declines being at the forefront of the global biodiversity crisis. The use of non-invasive methods to determine the physiological state in response to environmental stressors is therefore an important advance in the field of conservation physiology. The glucocorticoid hormone corticosterone (CORT) is one useful biomarker to assess physiological stress in amphibians, and sampling water-borne (WB) CORT is a novel, non-invasive collection technique. Here, we tested whether WB CORT can serve as a valid proxy of organismal levels of CORT in larvae of the common frog (Rana temporaria). We evaluated the association between tissue and WB CORT levels sampled from the same individuals across ontogenetic stages, ranging from newly hatched larvae to froglets at 10 days after metamorphosis. We also investigated how both tissue and WB CORT change throughout ontogeny. We found that WB CORT is a valid method in pro-metamorphic larvae as values for both methods were highly correlated. In contrast, there was no correlation between tissue and WB CORT in newly hatched, pre-metamorphic larvae, metamorphs or post-metamorphic froglets probably due to ontogenetic changes in respiratory and skin morphology and physiology affecting the transdermal CORT release. Both collection methods consistently revealed a non-linear pattern of ontogenetic change in CORT with a peak at metamorphic climax. Thus, our results indicate that WB CORT sampling is a promising, non-invasive conservation tool for studies on late-stage amphibian larvae. However, we suggest considering that different contexts might affect the reliability of WB CORT and consequently urge future studies to validate this method whenever it is used in new approaches. We conclude proposing some recommendations and perspectives on the use of WB CORT that will aid in broadening its application as a non-invasive tool in amphibian conservation physiology.

Life in plastic, it's not fantastic: Sublethal effects of polyethylene microplastics ingestion throughout amphibian metamorphosis

Microplastics (MP) are an abundant, long-lasting, and widespread type of environmental pollution that is of increasing concern as it might pose a serious threat to ecosystems and species. However, these threats are still largely unknown for amphibians. Here, we used the African clawed frog (Xenopus laevis) as a model species to investigate whether polyethylene MP ingestion affects amphibian growth and development and leads to metabolic changes across two consecutive life stages (larvae and juveniles). Furthermore, we examined whether MP effects were more pronounced at higher rearing temperatures. Larval growth, development, and body condition were recorded, and standard metabolic rate (SMR) and levels of stress hormone (corticosterone, CORT) were measured. We determined variation in size, morphology, and hepatosomatic index in juveniles to identify any potential consequences of MP ingestion across metamorphosis. In both life stages, MP accumulation in the body was assessed. MP ingestion was found to result in sublethal effects on larval growth, development, and metabolism, to lead to allometric carry-over effects on juvenile morphology, and to accumulate in the specimens at both life stages. In larvae, SMR and developmental rate increased in response to MP ingestion; there additionally was a significant interaction of MP ingestion and temperature on development. CORT levels were higher in larvae that ingested MP, except at higher temperature. In juveniles, body was wider, and extremities were longer in animals exposed to MP during the larval stage; a high rearing temperature in combination with MP ingestion counteracted this effect. Our results provide first insights into the effects of MP on amphibians throughout metamorphosis and demonstrate that juvenile amphibians may act as a pathway for MP from freshwater to terrestrial environments. To allow for generalizations across amphibian species, future experiments need to consider the field prevalence and abundance of different MP in amphibians at various life stages.

Does the Glucocorticoid Stress Response Make Toads More Toxic? An Experimental Study on the Regulation of Bufadienolide Toxin Synthesis

Chemical defense is a crucial component of fitness in many organisms, yet the physiological regulation of defensive toxin synthesis is poorly understood, especially in vertebrates. Bufadienolides, the main defensive compounds of toads, are toxic to many predators and other natural enemies, and their synthesis can be upregulated by stressors, including predation risk, high conspecific density, and pollutants. Thus, higher toxin content may be the consequence of a general endocrine stress response in toads. Therefore, we hypothesized that bufadienolide synthesis may be stimulated by elevated levels of corticosterone (CORT), the main glucocorticoid hormone of amphibians, or by upstream regulators that stimulate CORT production. To test these alternatives, we treated common toad tadpoles with exogenous CORT (exoCORT) or metyrapone (MTP, a CORT-synthesis inhibitor that stimulates upstream regulators of CORT by negative feedback) in the presence or absence of predation cues for 2 or 6 days, and subsequently measured their CORT release rates and bufadienolide content. We found that CORT release rates were elevated by exoCORT, and to a lesser extent also by MTP, regardless of treatment length. Bufadienolide content was significantly decreased by treatment with exoCORT for 6 days but was unaffected by exposure to exoCORT for 2 days or to MTP for either 6 or 2 days. The presence or absence of predation cues affected neither CORT release rate nor bufadienolide content. Our results suggest that changes in bufadienolide synthesis in response to environmental challenges are not driven by CORT but may rather be regulated by upstream hormones of the stress response.

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.

Contributions of water-borne corticosterone as one non-invasive biomarker in assessing nitrate pollution stress in tadpoles of Rana temporaria

Among a multitude of stressors to which wildlife is exposed, environmental pollution is a pervasive one that poses a serious threat. The permeable skin of amphibians is likely to increase direct contact of the body with pollutants, making them a group worth studying to access environmental quality. Consequently, finding reliable and complementary biomarkers that will present detectable and predictable changes in response to pollutants is essential to identify pollution sublethal effects on amphibians and to investigate whether these are in part responsible for population declines. The glucocorticoid hormone corticosterone (CORT), involved in many metabolic functions, is often used to measure the physiological stress response to environmental stressors in amphibians. In this study, we evaluated whether water-borne CORT can serve as a non-invasive biomarker for nitrate pollution stress in the European common frog (Rana temporaria) by comparing the effect of nitrate exposure on hormone release rates and on other physiological downstream biomarkers, i.e., ultimate physiological effects of the stressor. Specifically, we investigated the effect of different nitrate concentrations (0, 10, 50, and 100 mg/L) on water-borne CORT release rates, age, size, and body condition. Exposure to nitrate pollution significantly increased age at metamorphosis and water-borne CORT release rates, and led to reduced mass and body condition, but only at higher nitrate concentrations (i.e., 50 and 100 mg/L). Considering this similar sensitivity to other acknowledged biomarkers, water-borne CORT was a reliable biomarker of physiological stress in R. temporaria exposed to nitrate pollution stress in a controlled single-stressor laboratory approach. Thus, water-borne CORT is a promising method to be included in more holistic approaches. We recommend that such approaches keep testing multiple biomarker combinations, as species are exposed to several stressors likely to interact and produce varied outcomes in different biomarkers in their natural habitats.

Chromosome-level assembly, annotation and phylome of Pelobates cultripes, the western spadefoot toad

Genomic resources for amphibians are still hugely under-represented in vertebrate genomic research, despite being a group of major interest for ecology, evolution and conservation. Amphibians constitute a highly threatened group of vertebrates, present a vast diversity in reproductive modes, are extremely diverse in morphology, occupy most ecoregions of the world, and present the widest range in genome sizes of any major group of vertebrates. We combined Illumina, Nanopore and Hi-C sequencing technologies to assemble a chromosome-level genome sequence for an anuran with a moderate genome size (assembly span 3.09 Gb); Pelobates cultripes, the western spadefoot toad. The genome has an N50 length of 330 Mb with 98.6% of the total sequence length assembled into 14 super scaffolds, and 87.7% complete BUSCO genes. We use published transcriptomic data to provide annotations, identifying 32,684 protein-coding genes. We also reconstruct the P. cultripes phylome and identify 2,527 gene expansions. We contribute the first draft of the genome of the western spadefoot toad, P. cultripes. This species represents a relatively basal lineage in the anuran tree with an interesting ecology and a high degree of developmental plasticity, and thus is an important resource for amphibian genomic research.

Reproductive colonization of land by frogs: Embryos and larvae excrete urea to avoid ammonia toxicity

Vertebrate colonization of land has occurred multiple times, including over 50 origins of terrestrial eggs in frogs. Some environmental factors and phenotypic responses that facilitated these transitions are known, but responses to water constraints and risk of ammonia toxicity during early development are poorly understood. We tested if ammonia accumulation and dehydration risk induce a shift from ammonia to urea excretion during early stages of four anurans, from three origins of terrestrial development. We quantified ammonia and urea concentrations during early development on land, under well‐hydrated and dry conditions. Where we found urea excretion, we tested for a plastic increase under dry conditions and with ammonia accumulation in developmental environments. We assessed the potential adaptive role of urea excretion by comparing ammonia tolerance measured in 96h‐LC₅₀ tests with ammonia levels in developmental environments. Ammonia accumulated in foam nests and perivitelline fluid, increasing over development and reaching higher concentrations under dry conditions. All four species showed high ammonia tolerance, compared to fishes and aquatic‐breeding frogs. Both nest‐dwelling larvae of Leptodactylus fragilis and late embryos of Hyalinobatrachium fleischmanni excreted urea, showing a plastic increase under dry conditions. These two species can develop the longest on land and urea excretion appears adaptive, preventing their exposure to potentially lethal levels of ammonia. Neither late embryos of Agalychnis callidryas nor nest‐dwelling larvae of Engystomops pustulosus experienced toxic ammonia levels under dry conditions, and neither excreted urea. Our results suggest that an early onset of urea excretion, its increase under dry conditions, and elevated ammonia tolerance can all help prevent ammonia toxicity during terrestrial development. High ammonia represents a general risk for development which may be exacerbated as climate change increases dehydration risk for terrestrial‐breeding frogs. It may also be a cue that elicits adaptive physiological responses during early development.

Context dependent variation in corticosterone and phenotypic divergence of Rana arvalis populations along an acidification gradient

Background

Physiological processes, as immediate responses to the environment, are important mechanisms of phenotypic plasticity and can influence evolution at ecological time scales. In stressful environments, physiological stress responses of individuals are initiated and integrated via the release of hormones, such as corticosterone (CORT). In vertebrates, CORT influences energy metabolism and resource allocation to multiple fitness traits (e.g. growth and morphology) and can be an important mediator of rapid adaptation to environmental stress, such as acidification. The moor frog, Rana arvalis, shows adaptive divergence in larval life-histories and predator defense traits along an acidification gradient in Sweden. Here we take a first step to understanding the role of CORT in this adaptive divergence. We conducted a fully factorial laboratory experiment and reared tadpoles from three populations (one acidic, one neutral and one intermediate pH origin) in two pH treatments (Acid versus Neutral pH) from hatching to metamorphosis. We tested how the populations differ in tadpole CORT profiles and how CORT is associated with tadpole life-history and morphological traits.

Results

We found clear differences among the populations in CORT profiles across different developmental stages, but only weak effects of pH treatment on CORT. Tadpoles from the acid origin population had, on average, lower CORT levels than tadpoles from the neutral origin population, and the intermediate pH origin population had intermediate CORT levels. Overall, tadpoles with higher CORT levels developed faster and had shorter and shallower tails, as well as shallower tail muscles.

Conclusions

Our common garden results indicate among population divergence in CORT levels, likely reflecting acidification mediated divergent selection on tadpole physiology, concomitant to selection on larval life-histories and morphology. However, CORT levels were highly environmental context dependent. Jointly these results indicate a potential role for CORT as a mediator of multi-trait divergence along environmental stress gradients in natural populations. At the same time, the population level differences and high context dependency in CORT levels suggest that snapshot assessment of CORT in nature may not be reliable bioindicators of stress.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-022-01967-1.

Constraints and Opportunities for the Evolution of Metamorphic Organisms in a Changing Climate

We argue that developmental hormones facilitate the evolution of novel phenotypic innovations and timing of life history events by genetic accommodation. Within an individual’s life cycle, metamorphic hormones respond readily to environmental conditions and alter adult phenotypes. Across generations, the many effects of hormones can bias and at times constrain the evolution of traits during metamorphosis; yet, hormonal systems can overcome constraints through shifts in timing of, and acquisition of tissue specific responses to, endocrine regulation. Because of these actions of hormones, metamorphic hormones can shape the evolution of metamorphic organisms. We present a model called a developmental goblet, which provides a visual representation of how metamorphic organisms might evolve. In addition, because developmental hormones often respond to environmental changes, we discuss how endocrine regulation of postembryonic development may impact how organisms evolve in response to climate change. Thus, we propose that developmental hormones may provide a mechanistic link between climate change and organismal adaptation.

Interdependence of Thyroid and Corticosteroid Signaling in Vertebrate Developmental Transitions

Post-embryonic acute developmental processes mainly allow the transition from one life stage in a specific ecological niche to the next life stage in a different ecological niche. Metamorphosis, an emblematic type of these post-embryonic developmental processes, has occurred repeatedly and independently in various phylogenetic groups throughout metazoan evolution, such as in cnidarian, insects, molluscs, tunicates, or vertebrates. This review will focus on metamorphoses and developmental transitions in vertebrates, including typical larval metamorphosis in anuran amphibians, larval and secondary metamorphoses in teleost fishes, egg hatching in sauropsids and birth in mammals. Two neuroendocrine axes, the hypothalamic-pituitary-thyroid and the hypothalamic-pituitary-adrenal/interrenal axes, are central players in the regulation of these life transitions. The review will address the molecular and functional evolution of these axes and their interactions. Mechanisms of integration of internal and environmental cues, and activation of these neuroendocrine axes represent key questions in an “eco-evo-devo” perspective of metamorphosis. The roles played by developmental transitions in the innovation, adaptation, and plasticity of life cycles throughout vertebrates will be discussed. In the current context of global climate change and habitat destruction, the review will also address the impact of environmental factors, such as global warming and endocrine disruptors on hypothalamic-pituitary-thyroid and hypothalamic-pituitary-adrenal/interrenal axes, and regulation of developmental transitions.

Transcriptome and Methylome Analysis Reveal Complex Cross-Talks between Thyroid Hormone and Glucocorticoid Signaling at Xenopus Metamorphosis

Background: Most work in endocrinology focus on the action of a single hormone, and very little on the cross-talks between two hormones. Here we characterize the nature of interactions between thyroid hormone and glucocorticoid signaling during Xenopus tropicalis metamorphosis. Methods: We used functional genomics to derive genome wide profiles of methylated DNA and measured changes of gene expression after hormonal treatments of a highly responsive tissue, tailfin. Clustering classified the data into four types of biological responses, and biological networks were modeled by system biology. Results: We found that gene expression is mostly regulated by either T₃ or CORT, or their additive effect when they both regulate the same genes. A small but non-negligible fraction of genes (12%) displayed non-trivial regulations indicative of complex interactions between the signaling pathways. Strikingly, DNA methylation changes display the opposite and are dominated by cross-talks. Conclusion: Cross-talks between thyroid hormones and glucocorticoids are more complex than initially envisioned and are not limited to the simple addition of their individual effects, a statement that can be summarized with the pseudo-equation: TH ∙ GC > TH + GC. DNA methylation changes are highly dynamic and buffered from genome expression.

Recombination facilitates genetic assimilation of new traits in gene regulatory networks

A new phenotypic variant may appear first in organisms through plasticity, that is, as a response to an environmental signal or other nongenetic perturbation. If such trait is beneficial, selection may increase the frequency of alleles that enable and facilitate its development. Thus, genes may take control of such traits, decreasing dependence on nongenetic disturbances, in a process called genetic assimilation. Despite an increasing amount of empirical studies supporting genetic assimilation, its significance is still controversial. Whether genetic assimilation is widespread depends, to a great extent, on how easily mutation and recombination reduce the trait's dependence on nongenetic perturbations. Previous research suggests that this is the case for mutations. Here we use simulations of gene regulatory network dynamics to address this issue with respect to recombination. We find that recombinant offspring of parents that produce a new phenotype through plasticity are more likely to produce the same phenotype without requiring any perturbation. They are also prone to preserve the ability to produce that phenotype after genetic and nongenetic perturbations. Our work also suggests that ancestral plasticity can play an important role for setting the course that evolution takes. In sum, our results indicate that the manner in which phenotypic variation maps unto genetic variation facilitates evolution through genetic assimilation in gene regulatory networks. Thus, we contend that the importance of this evolutionary mechanism should not be easily neglected.

Adaptive and nonadaptive plasticity in changing environments: Implications for sexual species with different life history strategies

Populations adapt to novel environmental conditions by genetic changes or phenotypic plasticity. Plastic responses are generally faster and can buffer fitness losses under variable conditions. Plasticity is typically modeled as random noise and linear reaction norms that assume simple one-to-one genotype-phenotype maps and no limits to the phenotypic response. Most studies on plasticity have focused on its effect on population viability. However, it is not clear, whether the advantage of plasticity depends solely on environmental fluctuations or also on the genetic and demographic properties (life histories) of populations. Here we present an individual-based model and study the relative importance of adaptive and nonadaptive plasticity for populations of sexual species with different life histories experiencing directional stochastic climate change. Environmental fluctuations were simulated using differentially autocorrelated climatic stochasticity or noise color, and scenarios of directional climate change. Nonadaptive plasticity was simulated as a random environmental effect on trait development, while adaptive plasticity as a linear, saturating, or sinusoidal reaction norm. The last two imposed limits to the plastic response and emphasized flexible interactions of the genotype with the environment. Interestingly, this assumption led to (a) smaller phenotypic than genotypic variance in the population (many-to-one genotype-phenotype map) and the coexistence of polymorphisms, and (b) the maintenance of higher genetic variation-compared to linear reaction norms and genetic determinism-even when the population was exposed to a constant environment for several generations. Limits to plasticity led to genetic accommodation, when costs were negligible, and to the appearance of cryptic variation when limits were exceeded. We found that adaptive plasticity promoted population persistence under red environmental noise and was particularly important for life histories with low fecundity. Populations producing more offspring could cope with environmental fluctuations solely by genetic changes or random plasticity, unless environmental change was too fast.

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.

Stress hormones mediate developmental plasticity in vertebrates with complex life cycles

The environment experienced by developing organisms can shape the timing and character of developmental processes, generating different phenotypes from the same genotype, each with different probabilities of survival and performance as adults. Chordates have two basic modes of development, indirect and direct. Species with indirect development, which includes most fishes and amphibians, have a complex life cycle with a free-swimming larva that is typically a growth stage, followed by a metamorphosis into the adult form. Species with direct development, which is an evolutionarily derived developmental mode, develop directly from embryo to the juvenile without an intervening larval stage. Among the best studied species with complex life cycles are the amphibians, especially the anurans (frogs and toads). Amphibian tadpoles are exposed to diverse biotic and abiotic factors in their developmental habitat. They have extensive capacity for developmental plasticity, which can lead to the expression of different, adaptive morphologies as tadpoles (polyphenism), variation in the timing of and size at metamorphosis, and carry-over effects on the phenotype of the juvenile/adult. The neuroendocrine stress axis plays a pivotal role in mediating environmental effects on amphibian development. Before initiating metamorphosis, if tadpoles are exposed to predators they upregulate production of the stress hormone corticosterone (CORT), which acts directly on the tail to cause it to grow, thereby increasing escape performance. When tadpoles reach a minimum body size to initiate metamorphosis they can vary the timing of transformation in relation to growth opportunity or mortality risk in the larval habitat. They do this by modulating the production of thyroid hormone (TH), the primary inducer of metamorphosis, and CORT, which synergizes with TH to promote tissue transformation. Hypophysiotropic neurons that release the stress neurohormone corticotropin-releasing factor (CRF) are activated in response to environmental stress (e.g., pond drying, food restriction, etc.), and CRF accelerates metamorphosis by directly inducing secretion of pituitary thyrotropin and corticotropin, thereby increasing secretion of TH and CORT. Although activation of the neuroendocrine stress axis promotes immediate survival in a deteriorating larval habitat, costs may be incurred such as reduced tadpole growth and size at metamorphosis. Small size at transformation can impair performance of the adult, reducing probability of survival in the terrestrial habitat, or fecundity. Furthermore, elevations in CORT in the tadpole caused by environmental stressors cause long term, stable changes in neuroendocrine function, behavior and physiology of the adult, which can affect fitness. Comparative studies show that the roles of stress hormones in developmental plasticity are conserved across vertebrate taxa including humans.

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.

Corticosterone Is Essential for Survival Through Frog Metamorphosis

Thyroid hormone (TH) is required for frog metamorphosis, and corticosterone (CORT) increases TH signaling to accelerate metamorphic progression. However, a requirement for CORT in metamorphosis has been difficult to assess prior to the recent development of gene-editing technologies. We addressed this long-standing question using transcription activator-like effector nuclease (TALEN) gene disruption to knock out proopiomelanocortin (pomc) and disrupt CORT production in Xenopus tropicalis. As expected, mutant tadpoles had a reduced peak of plasma CORT at metamorphosis with correspondingly reduced expression of the CORT-response gene Usher syndrome type-1G (ush1g). Mutants had reduced rates of growth and development and exhibited lower expression levels of 2 TH response genes, Krüppel-like factor 9 (klf9) and TH receptor β (thrb). In response to exogenous TH, mutants had reduced TH response gene induction and slower morphological change. Importantly, death invariably occurred during tail resorption, unless rescued by exogenous CORT and, remarkably, by exogenous TH. The ability of exogenous TH by itself to overcome death in pomc mutants indicates that the CORT-dependent increase in TH signaling may ensure functional organ transformation required for survival through metamorphosis and/or may shorten the nonfeeding metamorphic transition to avoid lethal inanition.

Shifts in the developmental rate of spadefoot toad larvae cause decreased complexity of post-metamorphic pigmentation patterns

Amphibian larvae are plastic organisms that can adjust their growth and developmental rates to local environmental conditions. The consequences of such developmental alterations have been studied in detail, both at the phenotypic and physiological levels. While largely unknown, it is of great importance to assess how developmental alterations affect the pigmentation pattern of the resulting metamorphs, because pigmentation is relevant for communication, mate choice, and camouflage and hence influences the overall fitness of the toads. Here we quantify the variation in several aspects of the pigmentation pattern of juvenile spadefoot toads experimentally induced to accelerate their larval development in response to decreased water level. It is known that induced developmental acceleration comes at the cost of reduced size at metamorphosis, higher metabolic rate, and increased oxidative stress. In this study, we show that spadefoot toads undergoing developmental acceleration metamorphosed with a less complex, more homogeneous, darker dorsal pattern consisting of continuous blotches, compared to the more contrasted pattern with segregated blotches and higher fractal dimension in normally developing individuals, and at a smaller size. We also observed a marked effect of population of origin in the complexity of the pigmentation pattern. Complexity of the post-metamorphic dorsal pigmentation could therefore be linked to pre-metamorphic larval growth and development.

Identification of candidate loci for adaptive phenotypic plasticity in natural populations of spadefoot toads

Phenotypic plasticity allows organisms to alter their phenotype in direct response to changes in the environment. Despite growing recognition of plasticity's role in ecology and evolution, few studies have probed plasticity's molecular bases-especially using natural populations. We investigated the genetic basis of phenotypic plasticity in natural populations of spadefoot toads (Spea multiplicata). Spea tadpoles normally develop into an "omnivore" morph that is favored in long-lasting, low-density ponds. However, if tadpoles consume freshwater shrimp or other tadpoles, they can alternatively develop (via plasticity) into a "carnivore" morph that is favored in ephemeral, high-density ponds. By combining natural variation in pond ecology and morph production with population genetic approaches, we identified candidate loci associated with each morph (carnivores vs. omnivores) and loci associated with adaptive phenotypic plasticity (adaptive vs. maladaptive morph choice). Our candidate morph loci mapped to two genes, whereas our candidate plasticity loci mapped to 14 genes. In both cases, the identified genes tended to have functions related to their putative role in spadefoot tadpole biology. Our results thereby form the basis for future studies into the molecular mechanisms that mediate plasticity in spadefoots. More generally, these results illustrate how diverse loci might mediate adaptive plasticity.

Nutrition-responsive gene expression and the developmental evolution of insect polyphenism

Nutrition-responsive development is a ubiquitous and highly diversified example of phenotypic plasticity, yet its underlying molecular and developmental mechanisms and modes of evolutionary diversification remain poorly understood. We measured genome-wide transcription in three closely related species of horned beetles exhibiting strikingly diverse degrees of nutrition responsiveness in the development of male weaponry. We show that (1) counts of differentially expressed genes between low- and high-nutritional backgrounds mirror species-specific degrees of morphological nutrition responsiveness; (2) evolutionary exaggeration of morphological responsiveness is underlain by both amplification of ancestral nutrition-responsive gene expression and recruitment of formerly low nutritionally responsive genes; and (3) secondary loss of morphological responsiveness to nutrition coincides with a dramatic reduction in gene expression plasticity. Our results further implicate genetic accommodation of ancestrally high variability of gene expression plasticity in both exaggeration and loss of nutritional plasticity, yet reject a major role of taxon-restricted genes in the developmental regulation and evolution of nutritional plasticity.

Hormones, developmental plasticity, and adaptive evolution: Endocrine flexibility as a catalyst for 'plasticity-first' phenotypic divergence

Explaining how populations adapt to environments is among the foremost objectives of evolutionary theory. Over generations, natural selection impels the phenotypic distribution of a population based on individual variation in phenotype and fitness. However, environmental conditions can also shape how individuals develop within their lifetime to influence which phenotypes are expressed in a population. It has been proposed that such environmentally-initiated phenotypic variation - also termed developmental plasticity - may enable adaptive evolution under some scenarios. As dynamic regulators of development and phenotypic expression, hormones are important physiological mediators of developmental plasticity. Patterns of hormone secretion, hormone transport, and the sensitivity of tissues to hormones can each be altered by environmental conditions, and understanding how endocrine regulation shapes phenotypic development in an ecologically-relevant context has much to contribute toward clarifying the role of plasticity in evolutionary adaptation. This article explores how the environmental sensitivity of endocrine regulation may facilitate 'plasticity-first' evolution by generating phenotypic variants that precede adaptation to altered or novel environments. Predictions arising from 'plasticity-first' evolution are examined in the context of thyroid hormone mediation of morphological plasticity in Cyprinodon pupfishes from the Death Valley region of California and Nevada, USA. This clade of extremophile fishes diversified morphologically over the last ~20,000 years, and observations that some populations experienced contemporary phenotypic differentiation under recent habitat change provide evidence that hormone-mediate plasticity preceded genetic assimilation of morphology in one of the region's species: the Devils Hole pupfish, Cyprinodon diabolis. This example illustrates how conceptualizing hormones not only as regulators of homeostasis, but also as developmental intermediaries between environment conditions and phenotypic variation at the individual-, population-, and species-levels can enrich our understanding of endocrine regulation both as a facilitator of phenotypic change under shifting environments, and as important proximate mechanisms that may initiate 'plasticity-first' evolutionary adaptation.

Cellular Innovation of the Cyanobacterial Heterocyst by the Adaptive Loss of Plasticity

Cellular innovation is central to biological diversification, yet its underlying mechanisms remain poorly understood [1]. One potential source of new cellular traits is environmentally induced phenotypic variation, or phenotypic plasticity. The plasticity-first hypothesis [2-4] proposes that natural selection can improve upon an ancestrally plastic phenotype to produce a locally adaptive trait, but the role of plasticity for adaptive evolution is still unclear [5-10]. Here, we show that a structurally novel form of the heterocyst, the specialized nitrogen-fixing cell of the multicellular cyanobacterium Fischerella thermalis, has evolved multiple times from ancestrally plastic developmental variation during adaptation to high temperature. Heterocyst glycolipids (HGs) provide an extracellular gas diffusion barrier that protects oxygen-sensitive nitrogenase [11, 12], and cyanobacteria typically exhibit temperature-induced plasticity in HG composition that modulates heterocyst permeability [13, 14]. By contrast, high-temperature specialists of F. thermalis constitutively overproduce glycolipid isomers associated with high temperature to levels unattained by plastic strains. This results in a less-permeable heterocyst, which is advantageous at high temperature but deleterious at low temperature for both nitrogen fixation activity and fitness. Our study illustrates how the origin of a novel cellular phenotype by the genetic assimilation and adaptive refinement of a plastic trait can be a source of biological diversity and contribute to ecological specialization.

De Novo Assembly and Annotation of the Larval Transcriptome of Two Spadefoot Toads Widely Divergent in Developmental Rate

Amphibians are highly vulnerable and diverse vertebrates for which we still have modest genomic resources. Amphibian larvae are key components of continental wetlands, where they have strong influences on energy fluxes, nutrient cycling, and community structure. Amphibian larvae are highly sensitive to environmental conditions and can often alter their physiology, behavior and even morphology in response to the local conditions experienced, although we still know relatively little about the transcriptomic changes that enable such plasticity. Here we contribute the larval transcriptomes of two spadefoot toad species with divergent developmental rates and degree of developmental plasticity in response to pond drying.

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.

Insufficiency of Thyroid Hormone in Frog Metamorphosis and the Role of Glucocorticoids

Thyroid hormone (TH) is the most important hormone in frog metamorphosis, a developmental process which will not occur in the absence of TH but can be induced precociously by exogenous TH. However, such treatments including in-vitro TH treatments often do not replicate the events of natural metamorphosis in many organs, including lung, brain, blood, intestine, pancreas, tail, and skin. A potential explanation for the discrepancy between natural and TH-induced metamorphosis is the involvement of glucocorticoids (GCs). GCs are not able to advance development by themselves but can modulate the rate of developmental progress induced by TH via increased tissue sensitivity to TH. Global gene expression analyses and endocrine experiments suggest that GCs may also have direct actions required for completion of metamorphosis independent of their effects on TH signaling. Here, we provide a new review and analysis of the requirement and necessity of TH signaling in light of recent insights from gene knockout frogs. We also examine the independent and interactive roles GCs play in regulating morphological and molecular metamorphic events dependent upon TH.

How stabilizing selection and nongenetic inheritance combine to shape the evolution of phenotypic plasticity

Relatively little is known about whether and how nongenetic inheritance interacts with selection to impact the evolution of phenotypic plasticity. Here, we empirically evaluated how stabilizing selection and a common form of nongenetic inheritance-maternal environmental effects-jointly influence the evolution of phenotypic plasticity in natural populations of spadefoot toads. We compared populations that previous fieldwork has shown to have evolved conspicuous plasticity in resource-use phenotypes ("resource polyphenism") with those that, owing to stabilizing selection favouring a narrower range of such phenotypes, appear to have lost this plasticity. We show that: (a) this apparent loss of plasticity in nature reflects a condition-dependent maternal effect and not a genetic loss of plasticity, that is "genetic assimilation," and (b) this plasticity is not costly. By shielding noncostly plasticity from selection, nongenetic inheritance generally, and maternal effects specifically, can preclude genetic assimilation from occurring and consequently impede adaptive (genetic) evolution.

Comparison the sensitivity of amphibian metamorphosis assays with NF 48 stage and NF 51 stage Xenopus laevis tadpoles

The amphibian metamorphosis assay (AMA) was proposed by the Organization for Economic Cooperation and Development (OECD) to screen thyroid disruptors of vertebrate species. The general experimental design of the AMA exposes Nieuwkoop and Faber (NF) stage 51 Xenopus laevis tadpoles to test chemical concentrations for 21 d. However, recent studies demonstrated that thyroid gland began to function after NF stage 45 in X. laevis. Thus, in this study, we initiated exposure with NF stage 48 tadpoles when the thyroid gland is still in a preliminary development period, to compare the sensitivity of the AMA with NF 48 stage and NF 51 stage tadpoles. Further, the application and sensitivity of the optimized AMA were evaluated and validated by two known thyroid toxicants methimazole (MMI) and sodium perchlorate (SP). The observational endpoints are developmental stage, hind limb length (HLL), snout-vent length (SVL), wet weight, and daily observations of mortality. The results were as follows. Although the sensitivity to endpoint of growth, such as wet weight and SVL was similar between the two assays, our optimized AMA detected delaying effects of 1 mg/L MMI and 32 μg/L SP on metamorphosis development both on day 7 and at test termination, which were lower than those in AMA. Additionally, it is easier to get a large number of animals at NF stage 48 than NF stage 51 in a short time. Thus, it is suggested that the NF stage 48 tadpoles might be applied to the AMA for efficiently screening the thyroid-active substances.

Plasticity-led evolution: evaluating the key prediction of frequency-dependent adaptation

Plasticity-led evolution occurs when a change in the environment triggers a change in phenotype via phenotypic plasticity, and this pre-existing plasticity is subsequently refined by selection into an adaptive phenotype. A critical, but largely untested prediction of plasticity-led evolution (and evolution by natural selection generally) is that the rate and magnitude of evolutionary change should be positively associated with a phenotype's frequency of expression in a population. Essentially, the more often a phenotype is expressed and exposed to selection, the greater its opportunity for adaptive refinement. We tested this prediction by competing against each other spadefoot toad tadpoles from different natural populations that vary in how frequently they express a novel, environmentally induced carnivore ecomorph. As expected, laboratory-reared tadpoles whose parents were derived from populations that express the carnivore ecomorph more frequently were superior competitors for the resource for which this ecomorph is specialized-fairy shrimp. These tadpoles were better at using this resource both because they were more efficient at capturing and consuming shrimp and because they produced more exaggerated carnivore traits. Moreover, they exhibited these more carnivore-like features even without experiencing the inducing cue, suggesting that this ecomorph has undergone an extreme form of plasticity-led evolution-genetic assimilation. Thus, our findings provide evidence that the frequency of trait expression drives the magnitude of adaptive refinement, thereby validating a key prediction of plasticity-led evolution specifically and adaptive evolution generally.

Developmental Bias and Evolution: A Regulatory Network Perspective

Phenotypic variation is generated by the processes of development, with some variants arising more readily than others-a phenomenon known as "developmental bias." Developmental bias and natural selection have often been portrayed as alternative explanations, but this is a false dichotomy: developmental bias can evolve through natural selection, and bias and selection jointly influence phenotypic evolution. Here, we briefly review the evidence for developmental bias and illustrate how it is studied empirically. We describe recent theory on regulatory networks that explains why the influence of genetic and environmental perturbation on phenotypes is typically not uniform, and may even be biased toward adaptive phenotypic variation. We show how bias produced by developmental processes constitutes an evolving property able to impose direction on adaptive evolution and influence patterns of taxonomic and phenotypic diversity. Taking these considerations together, we argue that it is not sufficient to accommodate developmental bias into evolutionary theory merely as a constraint on evolutionary adaptation. The influence of natural selection in shaping developmental bias, and conversely, the influence of developmental bias in shaping subsequent opportunities for adaptation, requires mechanistic models of development to be expanded and incorporated into evolutionary theory. A regulatory network perspective on phenotypic evolution thus helps to integrate the generation of phenotypic variation with natural selection, leaving evolutionary biology better placed to explain how organisms adapt and diversify.

Morphological novelty emerges from pre-existing phenotypic plasticity

Plasticity-first evolution (PFE) posits that novel features arise when selection refines pre-existing phenotypic plasticity into an adaptive phenotype. However, PFE is controversial because few tests have been conducted in natural populations. Here we present evidence that PFE fostered the origin of an evolutionary novelty that allowed certain amphibians to invade a new niche-a distinctive carnivore morph. We compared morphology, gene expression and growth of three species of spadefoot toad tadpoles when reared on alternative diets: Scaphiopus holbrookii, which (like most frogs) never produce carnivores; Spea multiplicata, which sometimes produce carnivores, but only through diet-induced plasticity; and Spea bombifrons, which often produce carnivores regardless of diet. Consistent with PFE, we found diet-induced plasticity-in morphology and gene expression-in Sc. holbrookii, adaptive refinement of this plasticity in Sp. multiplicata, and further refinement of the carnivore phenotype in Sp. bombifrons. Generally, phenotypic plasticity might play a significant, if underappreciated, role in evolutionary innovation.

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.