Evolution Of Plasticity Research Articles

Phenotypic plasticity in tropical butterflies is linked to climatic seasonality on a macroevolutionary scale.

Phenotypic plasticity can be adaptive in fluctuating environments by providing rapid environment-phenotype matching and this applies particularly in seasonal environments. African Bicyclus butterflies have repeatedly colonized seasonal savannahs from ancestral forests around the late Miocene, and many species now exhibit seasonal polyphenism. On a macroevolutionary scale, it can be expected that savannah species will exhibit higher plasticity because of experiencing stronger environmental seasonality than forest species. We quantified seasonality using environmental niche modeling and surveyed the degree of plasticity in a key wing pattern element (eyespot size) using museum specimens. We showed that species occurring in highly seasonal environments display strong plasticity, while species in less seasonal or aseasonal environments exhibit surprisingly variable degrees of plasticity, including strong to no plasticity. Furthermore, eyespot size plasticity has a moderate phylogenetic signal and the ancestral Bicyclus likely exhibited some degree of plasticity. We propose hypotheses to explain the range of plasticity patterns seen in less seasonal environments and generate testable predictions for the evolution of plasticity in Bicyclus. Our study provides one of the most compelling cases showing links between seasonality and phenotypic plasticity on a macroevolutionary scale and the potential role of plasticity in facilitating the colonization of novel environments.

Fine scale diversity in the lava: genetic and phenotypic diversity in small populations of Arctic charr Salvelinus alpinus.

A major goal in evolutionary biology is to understand the processes underlying phenotypic variation in nature. Commonly, studies have focused on large interconnected populations or populations found along strong environmental gradients. However, studies on small fragmented populations can give strong insight into evolutionary processes in relation to discrete ecological factors. Evolution in small populations is believed to be dominated by stochastic processes, but recent work shows that small populations can also display adaptive phenotypic variation, through for example plasticity and rapid adaptive evolution. Such evolution takes place even though there are strong signs of historical bottlenecks and genetic drift. Here we studied 24 small populations of the freshwater fish Arctic charr (Salvelinus alpinus) found in groundwater filled lava caves. Those populations were found within a few km2-area with no apparent water connections between them. We studied the relative contribution of neutral versus non-neutral evolutionary processes in shaping phenotypic divergence, by contrasting patterns of phenotypic and neutral genetic divergence across populations in relation to environmental measurements. This allowed us to model the proportion of phenotypic variance explained by the environment, taking in to account the observed neutral genetic structure. These populations originated from the nearby Lake Mývatn, and showed small population sizes with low genetic diversity. Phenotypic variation was mostly correlated with neutral genetic diversity with only a small environmental effect. Phenotypic diversity in these cave populations appears to be largely the product of neutral processes, fitting the classical evolutionary expectations. However, the fact that neutral processes did not explain fully the phenotypic patterns suggests that further studies can increase our understanding on how neutral evolutionary processes can interact with other forces of selection at early stages of divergence. The accessibility of these populations has provided the opportunity for long-term monitoring of individual fish, allowing tracking how the environment can influence phenotypic and genetic divergence for shaping and maintaining diversity in small populations. Such studies are important, especially in freshwater, as habitat alteration is commonly breaking populations into smaller units, which may or may not be viable.

Evolutionary divergence of plasticity in brain morphology between ecologically divergent habitats of Trinidadian guppies.

Phenotypic plasticity is critical for organismal performance and can evolve in response to natural selection. Brain morphology is often developmentally plastic, affecting animal performance in a variety of contexts. However, the degree to which the plasticity of brain morphology evolves has rarely been explored. Here, we use Trinidadian guppies (Poecilia reticulata), which are known for their repeated adaptation to high-predation (HP) and low-predation (LP) environments, to examine the evolution and plasticity of brain morphology. We exposed second-generation offspring of individuals from HP and LP sites to 2 different treatments: predation cues and conspecific social environment. Results show that LP guppies had greater plasticity in brain morphology compared to their ancestral HP population, suggesting that plasticity can evolve in response to environmentally divergent habitats. We also show sexual dimorphism in the plasticity of brain morphology, highlighting the importance of considering sex-specific variation in adaptive diversification. Overall, these results may suggest the evolution of brain morphology plasticity as an important mechanism that allows for ecological diversification and adaptation to divergent habitats.

Enhanced mobility of dislocation network nodes and its effect on dislocation multiplication and strain hardening

Understanding plastic deformation of crystals in terms of the fundamental physics of dislocations has remained a grand challenge in materials science for decades. To overcome this, the Discrete Dislocation Dynamics (DDD) method has been developed, but its lack of atomistic resolution leaves open the possibility that certain key mechanisms may be overlooked. By comparing large-scale Molecular Dynamics (MD) with DDD simulations performed under identical conditions we uncover significant discrepancies in the predicted strength and microstructure evolution in BCC crystals under high-strain rate conditions. These are traced to unexpected behaviors of dislocation network nodes forming at dislocation intersections, that can move in ways not previously anticipated as revealed by MD. Once these newfound freedoms of nodal motion are incorporated, DDD simulations begin to closely match plastic evolution observed in MD. This additional mechanism of motion whereby non-screw dislocations can change their glide plane profoundly affects fundamental processes of dislocation multiplication, recovery and storage that define strength of metals.

Research on Multi-Physical Field Characteristics of Deep Coal Seam Mining Based on the Rock-Coal-Rock Model

In order to disclose the multi-physical field characteristics of the deep coal seam mining process and their dynamic evolution legislation, based on the “rock-coal-rock” model, during the mining process, the stress field, displacement field, energy field, and plastic zone evolution process are all simulated using FLAC3D6.0. The findings show that stress in the original rock is redistributed as a result of coal seam mining, creating a pressure relief zone in the middle of the goaf and advanced support pressure in the front part of the working face. The roof falls following the termination of coal seam mining. The collapsed blocks fill the middle of the goaf, playing a supporting role. The floor bulges as a new supporting pressure zone forms and builds up high elasticity. The stress reduction zone shifts from a rectangular to an inner circular distribution and an outer square as the working face’s mining distance increases and the range of the fracture field expands accordingly. In addition, a complete model was constructed to verify the correctness of the “rock-coal-rock” model. The stress, displacement, and energy curves of the overlying strata at a distance of 12 m from the bottom of the coal seam in the middle of the goaf obtained by the two methods were basically consistent. Ultimately, the findings of the numerical simulation were compared with the advanced support pressure data that were acquired on-site and they were good. This work can provide a reference for the safe mining of deep coal seams.

Evolution and plasticity of physiological traits in the collembolan Orchesella villosa at fine spatial scales within the city

Abstract Divergent natural selection caused by spatial environmental variation can lead to local adaptation and evolutionary divergence between populations, even those within close proximity to one another. This, however, is only one possibility among other outcomes, such as the evolution of adaptive phenotypic plasticity or local maladaptation. Cities, as fragmented landscapes with high environmental variability across microgeographic scales, present an excellent venue to explore these possible outcomes. Here, we use common garden experiments to explore the potential for microgeographic divergence of physiological traits among populations of a collembolan, Orchesella villosa, living within the footprint of a single city. We assessed if intra-urban variation in the intensity of the heat island effect and soil salinity has led to evolutionary divergence in heat tolerance and salinity tolerance, respectively. While we found little variation in salinity tolerance among populations, there was evidence for both plasticity and evolutionary divergence among populations in response to variation in the urban heat island. Although the adaptive nature of these specific responses is not known, we suggest that cities, and human modified habitats in general, promote evolutionary divergence at small spatial scales by creating barriers to dispersal and imposing divergent selective environments within previously contiguous habitats.

Emergence of phenotypic plasticity through epigenetic mechanisms.

Plasticity is found in all domains of life and is particularly relevant when populations experience variable environmental conditions. Traditionally, evolutionary models of plasticity are non-mechanistic: they typically view reactions norms as the target of selection, without considering the underlying genetics explicitly. Consequently, there have been difficulties in understanding the emergence of plasticity, and in explaining its limits and costs. In this paper, we offer a novel mechanistic approximation for the emergence and evolution of plasticity. We simulate random "epigenetic mutations" in the genotype-phenotype mapping, of the kind enabled by DNA-methylations/demethylations. The frequency of epigenetic mutations at loci affecting the phenotype is sensitive to organism stress (trait-environment mismatch), but is also genetically determined and evolvable. Thus, the "random motion" of epigenetic markers enables developmental learning-like behaviors that can improve adaptation within the limits imposed by the genotypes. However, with random motion being "goal-less," this mechanism is also vulnerable to developmental noise leading to maladaptation. Our individual-based simulations show that epigenetic mutations can hide alleles that are temporarily unfavorable, thus enabling cryptic genetic variation. These alleles can be advantageous at later times, under regimes of environmental change, in spite of the accumulation of genetic loads. Simulations also demonstrate that plasticity is favored by natural selection in constant environments, but more under periodic environmental change. Plasticity also evolves under directional environmental change as long as the pace of change is not too fast and costs are low.

Study of nano-broaching properties on nickel-based high-temperature alloys based on molecular dynamics simulation

Broaching technology, as the preferred high-volume machining process for nickel-based high-temperature alloy turbine disk tongue and groove machining, directly determines the reliability and stability of the turbine disk during its service. Although a host of scholars have carried out extensive experimental studies on the broaching process, there are few reports on the exploration of the atomic scale response in the broaching process. In this paper, the plastic deformation, defect evolution and distribution, surface formation, and stress changes during nano broaching of nickel-based single-crystal alloys are analyzed based on the molecular dynamics simulation (MD) method to reveal their broaching properties. The results show that the thermal softening effect induced by the accumulation of cutting temperature is the main reason for the violent plastic deformation within the matrix. The phase transition of the matrix structure, the depth of the subsurface damage layer, and the dislocation density increase with the increase of the cutting distance at the initial stage of the broaching process, and then slow down the upward trend at the stable broaching stage, and fluctuate in a certain range gradually. Meanwhile, the sources of atoms in chips produced by different broaching teeth in the broaching process alternates with the cutting distance. In addition, the diffusion depth of hydrostatic pressure and von Mises stress increases significantly with the increase of broaching distance. The simulation results can provide atomic-scale details for the broaching process of single-crystal Ni Fe Cr alloys and help to optimize the nano-broaching process of nickel-based high-temperature alloys.

Compression-Shear Specimen Stress-State Response and Distribution Characteristics with Wide Stress Triaxiality.

Conventional methods for studying the plastic behavior of materials involve uniaxial tension and uniaxial compression. However, in the metal rolling process, the deformation zone undergoes a complex loading of multidirectional compression and shear. Characterizing the corresponding plastic evolution process online poses challenges, and the existing specimen structures struggle to accurately replicate the deformation-induced loading characteristics. In this study, we aimed to design a compression-shear composite loading specimen that closely mimics the actual processing conditions. The goal was to investigate how the specimen structure influences the stress-strain response in the deformation zone. Using commercial finite element software, a compression-shear composite loading specimen was meticulously designed. Five 304 stainless steel specimens underwent uniaxial compressive loading, with variation angles between the preset notch angle (PNA) of the specimen and compression direction. We employed digital image correlation methods to capture the impact of the PNA on the strain field during compression. Additionally, we aimed to elucidate the plastic response resulting from the stress state of the specimen, particularly in relation to specimen fracture and microstructural evolution.

Plasticity, symbionts and niche construction interact in shaping dung beetle development and evolution.

Developmental plasticity is an important product of evolutionary processes, allowing organisms to maintain high fitness in the face of environmental perturbations. Once evolved, plasticity also has the potential to influence subsequent evolutionary outcomes, for example, by shaping phenotypic variation visible to selection and facilitating the emergence of novel trait variants. Furthermore, organisms may not just respond to environmental conditions through plasticity but may also actively modify the abiotic and (sym)biotic environments to which they themselves respond, causing plasticity to interact in complex ways with niche construction. Here, we explore developmental mechanisms and evolutionary consequences of plasticity in horned dung beetles. First, we discuss how post-invasion evolution of plasticity in an introduced Onthophagus species facilitated rapid range expansion and concurrent local adaptation of life history and morphology to novel climatic conditions. Second, we discuss how, in addition to plastically responding to variation in nutritional conditions, dung beetles engage in behaviors that modify the environment that they themselves respond to during later development. We document that these environment-modifying behaviors mask heritable variation for life history traits within populations, thereby shielding genetic variants from selection. Such cryptic genetic variation may be released and become selectable when these behaviors are compromised. Together, this work documents the complex interactions between plasticity, symbionts and niche construction, and highlights the usefulness of an integrative Eco-Evo-Devo framework to study the varied mechanisms and consequences of plasticity in development and evolution.

Genome evolution and divergence in cis-regulatory architecture is associated with condition-responsive development in horned dung beetles.

Phenotypic plasticity is thought to be an important driver of diversification and adaptation to environmental variation, yet the genomic mechanisms mediating plastic trait development and evolution remain poorly understood. The Scarabaeinae, or true dung beetles, are a species-rich clade of insects recognized for their highly diversified nutrition-responsive development including that of cephalic horns-evolutionarily novel, secondary sexual weapons that exhibit remarkable intra- and interspecific variation. Here, we investigate the evolutionary basis for horns as well as other key dung beetle traits via comparative genomic and developmental assays. We begin by presenting chromosome-level genome assemblies of three dung beetle species in the tribe Onthophagini (> 2500 extant species) including Onthophagus taurus, O. sagittarius, and Digitonthophagus gazella. Comparing these assemblies to those of seven other species across the order Coleoptera identifies evolutionary changes in coding sequence associated with metabolic regulation of plasticity and metamorphosis. We then contrast chromatin accessibility in developing head horn tissues of high- and low-nutrition O. taurus males and females and identify distinct cis-regulatory architectures underlying nutrition- compared to sex-responsive development, including a large proportion of recently evolved regulatory elements sensitive to horn morph determination. Binding motifs of known and new candidate transcription factors are enriched in these nutrition-responsive open chromatin regions. Our work highlights the importance of chromatin state regulation in mediating the development and evolution of plastic traits, demonstrates gene networks are highly evolvable transducers of environmental and genetic signals, and provides new reference-quality genomes for three species that will bolster future developmental, ecological, and evolutionary studies of this insect group.

Genome Assembly and Annotation of the Dark-Branded Bushbrown Butterfly Mycalesis mineus (Nymphalidae: Satyrinae).

We report a high-quality genome draft assembly of the dark-branded bushbrown, Mycalesis mineus, a member of the Satyrinae subfamily of nymphalid butterflies. This species is emerging as a promising model organism for investigating the evolution and development of phenotypic plasticity. Using 45.99 Gb of long-read data (N50 = 11.11 kb), we assembled a genome size of 497.4 Mb for M. mineus. The assembly is highly contiguous and nearly complete (96.8% of Benchmarking Universal Single-Copy Orthologs lepidopteran genes were complete and single copy). The genome comprises 38.71% of repetitive elements and includes 20,967 predicted protein-coding genes. The assembled genome was super-scaffolded into 28 pseudo-chromosomes using a closely related species, Bicyclus anynana, with a chromosomal-level genome as a template. This valuable genomic tool will advance both ongoing and future research focused on this model organism.

Anisotropic Hardening and Plastic Evolution Characterization on the Pressure-Coupled Drucker Yield Function of ZK61M Magnesium Alloy.

This paper studies the plastic behavior of the ZK61M magnesium alloy through a combination method of experiments and theoretical models. Based on a dog-bone specimen under different loading directions, mechanical tests under uniaxial tension were carried out, and the hardening behavior was characterized by the Swift-Voce hardening law. The von Mises yield function and the pressure-coupled Drucker yield function were used to predict the load-displacement curves of the ZK61M magnesium alloy under various conditions, respectively, where the material parameters were calibrated by using inverse engineering. The experimental results show that the hardening behavior of the ZK61M magnesium alloy has obvious anisotropy, but the effect of the stress state is more important on the strain hardening behavior of the alloy. Compared with the von Mises yield function, the pressure-coupled Drucker yield function is more accurate when characterizing the plastic behavior and strain hardening in different stress states of shear, uniaxial tension, and plane strain tension for the ZK61M alloy.

Convergent and environmentally associated chromatic polymorphism in Bryconops Kner, 1858 (Ostariophysi: Characiformes: Iguanodectidae).

Bryconops Kner, 1858, includes two well defined subgenera based on morphological evidence, with each containing at least one species (B. (Bryconops) caudomaculatus and B. (Creatochanes) melanurus) with a very wide distribution, within which regional populations present color variations. To test if phenotypic variation is related to cladogenetic events, we performed tests for phylogenetic independence and determined the strength of convergence for color characters in relation to water type, as the variation between clear, black and white waters is considered to be one of the major driving forces in the evolution of Amazonian fishes. Color characters for fins above the median line of the body were generally found to be independent from phylogeny and the Wheatsheaf test strongly supports convergence of the dorsal fin color between populations of species in the same type of water, with a similar trend suggested for the color of the dorsal lobe of the caudal fin. This means that simple color characters cannot necessarily be relied upon for taxonomic revisions of the genus as local phenotypic variants may represent environmentally determined plasticity or convergent evolution. Further studies are required to determine the validity of these characters.

Cell fate regulation governed by p53: Friends or reversible foes in cancer therapy.

Cancer is a leading cause of death worldwide. Targeted therapies aimed at key oncogenic driver mutations in combination with chemotherapy and radiotherapy as well as immunotherapy have benefited cancer patients considerably. Tumor protein p53 (TP53), a crucial tumor suppressor gene encoding p53, regulates numerous downstream genes and cellular phenotypes in response to various stressors. The affected genes are involved in diverse processes, including cell cycle arrest, DNA repair, cellular senescence, metabolic homeostasis, apoptosis, and autophagy. However, accumulating recent studies have continued to reveal novel and unexpected functions of p53 in governing the fate of tumors, for example, functions in ferroptosis, immunity, the tumor microenvironment and microbiome metabolism. Among the possibilities, the evolutionary plasticity of p53 is the most controversial, partially due to the dizzying array of biological functions that have been attributed to different regulatory mechanisms of p53 signaling. Nearly 40 years after its discovery, this key tumor suppressor remains somewhat enigmatic. The intricate and diverse functions of p53 in regulating cell fate during cancer treatment are only the tip of the iceberg with respect to its equally complicated structural biology, which has been painstakingly revealed. Additionally, TP53 mutation is one of the most significant genetic alterations in cancer, contributing to rapid cancer cell growth and tumor progression. Here, we summarized recent advances that implicate altered p53 in modulating the response to various cancer therapies, including chemotherapy, radiotherapy, and immunotherapy. Furthermore, we also discussed potential strategies for targeting p53 as a therapeutic option for cancer.

Influence of principal stress orientation on stress distribution and plastic zone evolution of rock surrounding tunnels

During the excavation of tunnels, the principal stress orientation changes, with a significant impact on the stress distribution, mechanical properties, and plastic zone evolution of rocks surrounding tunnels, causing severe deformation control, and monitoring problems in the stability of tunnels. Currently, biaxial compression tests were conducted to explore the influence of principal stress orientation on mechanical properties of rocks surrounding tunnels. The analytical solution of stress and the model of plastic zone of rocks considering the principal stress orientation and the distance from the excavation boundary were established to reveal the failure mechanism of surrounding rock under different principal stress orientations. With an increase in the angle between the principal stress orientation and the long axis of tunnel, the maximum tangential stress around the tunnel gradually changed to the minimum horizontal principal stress direction, and its value gradually increased, leading to quicker failure of surrounding rocks, and reducing the strength enhancement effect of the same section size of the tunnel. However, the increase in the angle reduced the damage range and the range of the plastic zone around the tunnel and caused the plastic zone to gradually approach the bottom and roof. The maximum depth of the plastic zone remained parallel to or nearly parallel to the minimum horizontal principal stress direction. When the principal stress orientation was kept constant, the maximum depth of the plastic zone shifted to the minimum horizontal principal stress direction with an increase in the vertical principal stress.

Modeling stress-induced responses: plasticity in continuous state space and gradual clonal evolution.

Mathematical models of cancer and bacterial evolution have generally stemmed from a gene-centric framework, assuming clonal evolution via acquisition of resistance-conferring mutations and selection of their corresponding subpopulations. More recently, the role of phenotypic plasticity has been recognized and models accounting for phenotypic switching between discrete cell states (e.g., epithelial and mesenchymal) have been developed. However, seldom do models incorporate both plasticity and mutationally driven resistance, particularly when the state space is continuous and resistance evolves in a continuous fashion. In this paper, we develop a framework to model plastic and mutational mechanisms of acquiring resistance in a continuous gradual fashion. We use this framework to examine ways in which cancer and bacterial populations can respond to stress and consider implications for therapeutic strategies. Although we primarily discuss our framework in the context of cancer and bacteria, it applies broadly to any system capable of evolving via plasticity and genetic evolution.

Instability mechanism and control technology of soft rock roadways affected by mining in karst mountain area

The deformation and damage of soft rock roadways in karst mountainous area is serious, which greatly restricts the safe and efficient mining of coal resources. To explore the instability mechanism of soft rock roadway under the influence of mining, the 226 rail cross-cut of Shanjiaoshu Coal Mine in China was taken as an engineering background. The methods of field monitoring and drilling imaging were adopted to analyze the roadway failure characteristics. In addition, the distributions of the plastic zone, surface displacement, and fracture fractal evolution characteristics of the surrounding rock for different stress wave peak and stress wave frequency are compared and analyzed using UDEC software. Then, the combined support scheme of “secondary strong bolting and grouting” is proposed to control the severe deformation of the roadway. The primary support mainly adopts short hollow grouting anchor cable to strengthen the shallow broken surrounding rock within 0–3 m. The secondary reinforced support adopts long hollow grouting anchor cable or constant resistance large deformation anchor cable to fully mobilize the bearing capacity of deep surrounding rock and to ensure the long-term stability of the roadway. On-site monitoring showed that the scheme could effectively control the severe deformation of the roadway, which also provided a reference for other similar roadway support.

Plastic and damage energy dissipation characteristics and damage evolution of Beishan granite under triaxial cyclic loading

The damage and failure of rock under triaxial stress is an irreversible process of energy dissipation, which mainly includes plastic energy dissipation and damage energy dissipation. However, it is difficult and rare to separate plastic energy dissipation and damage energy dissipation and to conduct relevant research on this basis. In this work, a loading control method with a combination of load (0.5 kN/s) and circumferential deformation (0.01 mm/min) was used, and triaxial cyclic loading and unloading tests under five different confining pressures were performed on granite from the Beishan site area of China's high-level radioactive waste repository. The results indicate that the angle between the main fracture surface and the longitudinal axis increases with the increase of confining pressure. According to the characteristics of the stress–strain curve and energy dissipation, the test process can be divided into five stages: the initial compaction and elastic stage, the stable crack growth stage, the unstable crack growth stage, the postpeak unstable fracture stage and the residual stress stage. Importantly, a separate determination method of plastic energy dissipation and damage energy dissipation was proposed, and the evolution characteristics of plastic energy dissipation and damage energy dissipation were revealed. The ratio of plastic energy dissipation to input energy decreases slightly in the initial compaction and elastic stage, stabilizes at a low level in the stable and unstable crack growth stages, increases significantly in the postpeak unstable fracture stage and remains at a high level in the residual stress stage. Meanwhile, the ratio of damage energy dissipation to input energy increases slightly before the unstable crack growth stage, increases abruptly at the end of the unstable crack growth stage and the initial part of the postpeak unstable fracture stage, then declines rapidly, and finally remains at a low level. Furthermore, the qualitative and quantitative research results of damage evolution show that the plastic shear strain and damage variable established by damage energy dissipation can qualitatively and quantitatively describe the damage evolution of Beishan granite respectively. Finally, the verification results of Beishan granite and two other rock types show that the separation method for calculating plastic energy dissipation and damage energy dissipation proposed in this work has good adaptability, and the qualitative damage development analysis based on plastic shear strain can be verified and supplemented with the quantitative damage development analysis based on energy theory.

Molecular basis of phenotypic plasticity in a marine ciliate.

Phenotypic plasticity, which involves phenotypic transformation in the absence of genetic change, may serve as a strategy for organisms to survive in complex and highly fluctuating environments. However, its reaction norm, molecular basis, and evolution remain unclear in most organisms, especially microbial eukaryotes. In this study, we explored these questions by investigating the reaction norm, regulation, and evolution of phenotypic plasticity in the cosmopolitan marine free-living ciliates Glauconema spp., which undergo significant phenotypic changes in response to food shortages. This study led to the de novo assembly of macronuclear genomes using long-read sequencing, identified hundreds of differentially expressed genes associated with phenotypic plasticity in different life stages, validated the function of two of these genes, and revealed that the reaction norm of body shape in response to food density follows a power-law distribution. Purifying selection may be the dominant evolutionary force acting on the genes associated with phenotypic plasticity, and the overall data support the hypothesis that phenotypic plasticity is a trait maintained by natural selection. This study provides novel insight into the developmental genetics of phenotypic plasticity in non-model unicellular eukaryotes and sheds light on the complexity and long evolutionary history of this important survival strategy.