Developmental decoupling of alternative phenotypes: insights from the transcriptomes of horn-polyphenic beetles

Origin of division of labor is decoupled from polymorphism in colonial animals

Division of labor, the specialization of sometimes phenotypically divergent cell types or group members, is often associated with ecological success in eukaryotic colonial organisms. Despite its many independent evolutionary origins, how division of labor emerges remains unclear. Conventional hypotheses tend toward an "economic" model, so that biological division of labor may reflect a partitioning of preexisting tasks and morphologies into specialized colony members. Here, we present an alternative model of the origin of division of labor, which can explain the evolution of new functions within a colony. We show that in colonies of the Cretaceous aged (103-96 Ma) fossil bryozoan of the genus Wilbertopora, the first cheilostome bryozoan to evolve polymorphism, preexisting morphologies were not simply partitioned among new members, but instead expanded into novel morphospace as they lost functions, specifically feeding. This expansion occurred primarily during two pulses of heightened morphological disparity, suggesting that the evolution of polymorphism corresponded to relaxed constraints on morphology and perhaps to the exploration of novel functions. Using a simple model of physiological connections, we show that regardless of the functionality of these new colony members, all nonfeeding members could have been supported by neighboring feeding members. This suggests that geometric constraints and physiological connectedness could be prerequisites for evolving both polymorphism and division of labor in modular organisms, and that a classic partitioning model of specialization cannot be broadly applied to biological systems.

Early life nutrient restriction affects hypothalamic-pituitary-interrenal axis gene expression in a diet type-specific manner

Stressful experiences in early life can alter phenotypic expression later in life. For instance, in vertebrates, early life nutrient restriction can modify later life activity of the hypothalamic-pituitary-adrenal/interrenal axis (the HPI in amphibians), including the up- and downstream regulatory components of glucocorticoid signaling. Early life nutrient restriction can also influence later life behavior and metabolism (e.g., fat accumulation). Yet, less is known about whether nutrient stress-induced carryover effects on HPA/HPI axis regulation can vary across environmental contexts, such as the type of diet on which nutrient restriction occurs. Here, we experimentally address this question using the plains spadefoot toad (Spea bombifrons), whose larvae develop in ephemeral habitats that impose intense competition over access to two qualitatively distinct diet types: detritus and live shrimp prey. Consistent with diet type-specific carryover effects of early life nutrient restriction on later life HPI axis regulation, we found that temporary nutrient restriction at the larval stage reduced juvenile (i.e., post-metamorphic) brain gene expression of an upstream glucocorticoid regulator (corticotropin-releasing hormone) and two downstream regulators (glucocorticoid and mineralocorticoid receptors) only on the shrimp diet. These patterns are consistent with known diet type-specific effects of larval nutrient restriction on juvenile corticosterone and behavior. Additionally, larval nutrient restriction increased juvenile body fat levels. Our study indicates that HPA/HPI axis regulatory responses to nutrient restriction can vary remarkably across diet types. Such diet type-specific regulation of the HPA/HPI axis might provide a basis for developmental or evolutionary decoupling of stress-induced carryover effects.

Genetic regulators of a resource polyphenism interact to couple predatory morphology and behaviour

Phenotypic plasticity often requires the coordinated response of multiple traits observed individually as morphological, physiological or behavioural. The integration, and hence functionality, of this response may be influenced by whether and how these component traits share a genetic basis. In the case of polyphenism, or discrete plasticity, at least part of the environmental response is categorical, offering a simple readout for determining whether and to what degree individual components of a plastic response can be decoupled. Here, we use the nematode Pristionchus pacificus, which has a resource polyphenism allowing it to be a facultative predator of other nematodes, to understand the genetic integration of polyphenism. The behavioural and morphological consequences of perturbations to the polyphenism's genetic regulatory network show that both predatory activity and ability are strongly influenced by morphology, different axes of morphological variation are associated with different aspects of predatory behaviour, and rearing environment can decouple predatory morphology from behaviour. Further, we found that interactions between some polyphenism-modifying genes synergistically affect predatory behaviour. Our results show that the component traits of an integrated polyphenic response can be decoupled and, in principle, selected upon individually, and they suggest that multiple routes to functionally comparable phenotypes are possible.

Local adaptation to an altitudinal gradient: the interplay between mean phenotypic trait variation and phenotypic plasticity in Mimulus laciniatus

Organisms can adapt to environmental heterogeneity through two mechanisms: (1) expression of population genetic variation or (2) phenotypic plasticity. In this study we investigated whether patterns of variation in both trait means and phenotypic plasticity along elevational and latitudinal clines in a North American endemic plant, Mimulus laciniatus, were consistent with local adaptation. We grew inbred lines of M. laciniatus from across the species' range in two common gardens varying in day length to measure mean and plastic trait expression in several traits previously shown to be involved in adaptation to M. laciniatus's rocky outcrop microhabitat: flowering time, size-related traits, and leaf shape. We examined correlations between the mean phenotype and phenotypic plasticity, and tested for a relationship between trait variation and population elevation and latitude. We did not find a strong correlation between mean and plastic trait expression at the individual genotype level suggesting that they operate under independent genetic controls. We identified multiple traits that show patterns consistent with local adaptation to elevation: critical photoperiod, flowering time, flower size, mean leaf lobing, and leaf lobing plasticity. These trends occur along multiple geographically independent altitudinal clines indicating that selection is a more likely cause of this pattern than gene flow among nearby populations with similar trait values. We also found that population variation in mean leaf lobing is associated with latitude. Our results indicate that both having more highly lobed leaves and greater leaf shape plasticity may be adaptive at high elevation within M. laciniatus. Our data strongly suggest that traits known to be under divergent selection between M. laciniatus and close relative Mimulus guttatus are also under locally varying selection within M. laciniatus.

A Mediator subunit imparts robustness to a polyphenism decision

Polyphenism is a type of developmental plasticity that translates continuous environmental variability into discontinuous phenotypes. Such discontinuity likely requires a switch between alternative gene-regulatory networks, a principle that has been borne out by mechanisms found to promote morph-specific gene expression. However, whether robustness is required to execute a polyphenism decision has awaited testing at the molecular level. Here, we used a nematode model for polyphenism, Pristionchus pacificus, to identify the molecular regulatory factors that ensure the development of alternative forms. This species has a dimorphism in its adult feeding structures, specifically teeth, which are a morphological novelty that allows predation on other nematodes. Through a forward genetic screen, we determined that a duplicate homolog of the Mediator subunit MDT-15/MED15, P. pacificus MDT-15.1, is necessary for the polyphenism and the robustness of the resulting phenotypes. This transcriptional coregulator, which has a conserved role in metabolic responses to nutritional stress, coordinates these processes with its effects on this diet-induced polyphenism. Moreover, this MED15 homolog genetically interacts with two nuclear receptors, NHR-1 and NHR-40, to achieve dimorphism: Single and double mutants for these three factors result in morphologies that together produce a continuum of forms between the extremes of the polyphenism. In summary, we have identified a molecular regulator that confers discontinuity to a morphological polyphenism, while also identifying a role for MED15 as a plasticity effector.

Sex-specific morphs: the genetics and evolution of intra-sexual variation

Sex-specific morphs exhibit discrete phenotypes, often including many disparate traits, that are observed in only one sex. These morphs have evolved independently in many different animals and are often associated with alternative mating strategies. The remarkable diversity of sex-specific morphs offers unique opportunities to understand the genetic basis of complex phenotypes, as the distinct nature of many morphs makes it easier to both categorize and compare genomes than for continuous traits. Sex-specific morphs also expand the study of sexual dimorphism beyond traditional bimodal comparisons of male and female averages, as they allow for a more expansive range of sexualization. Although ecological and endocrinological studies of sex-specific morphs have been advancing for some time, genomic and transcriptomic studies of morphs are far more recent. These studies reveal not only many different paths to the evolution of sex-specific morphs but also many commonalities, such as the role of sex-determining genes and hormone signalling in morph development, and the mixing of male and female traits within some morphs.

Linking Molecular Mechanisms and Evolutionary Consequences of Resource Polyphenism

Resource polyphenism-the occurrence of environmentally induced, discrete, and intraspecific morphs showing differential niche use-is taxonomically widespread and fundamental to the evolution of ecological function where it has arisen. Despite longstanding appreciation for the ecological and evolutionary significance of resource polyphenism, only recently have its proximate mechanisms begun to be uncovered. Polyphenism switches, especially those influencing and influenced by trophic interactions, offer a route to integrating proximate and ultimate causation in studies of plasticity, and its potential influence on evolution more generally. Here, we use the major events in generalized polyphenic development as a scaffold for linking the molecular mechanisms of polyphenic switching with potential evolutionary outcomes of polyphenism and for discussing challenges and opportunities at each step in this process. Not only does the study of resource polyphenism uncover interesting details of discrete plasticity, it also illuminates and informs general principles at the intersection of development, ecology, and evolution.

Differences in Cell Proliferation and Craniofacial Phenotype of Closely Related Species in the Pupfish Genus Cyprinodon

Understanding the genetic basis for phenotypic differences is fundamental to the study of macroevolutionary patterns of biological diversity. While technological advances in DNA sequencing have made researching genetic variation in wild taxa routine, fully understanding how these variants affect phenotype requires taking the next step to investigate how genetic changes alter cell and tissue interactions that ultimately produce phenotypes. In this article, we investigate a role for cell proliferation as a developmental source of craniofacial diversity in a radiation of 3 species of Cyprinodon from San Salvador Island, Bahamas. Patterns of cell proliferation in the heads of hatching-age fish differ among species of Cyprinodon, and correlate with differences in allometric growth rate among the jaws of 3 distinct species. Regional patterns of cell proliferation in the head are complex, resulting in an unintuitive result in which lower levels of cell proliferation in the posterior head region are associated with the development of relatively larger jaws in one species. We combine these data with previously published morphological and genomic data to show how studying the mechanisms generating phenotype at the cellular and tissue levels of biological organization can help mechanistically link genomic studies with classic morphological studies.

Signals of selection beyond bottlenecks between exotic populations of the bull-headed dung beetle, Onthophagus taurus

Colonization of new environments can lead to population bottlenecks and rapid phenotypic evolution that could be due to neutral and selective processes. Exotic populations of the bull-headed dung beetle (Onthophagus taurus) have differentiated in opposite directions from native beetles in male horn-to-body size allometry and female fecundity. Here we test for genetic and transcriptional differences among two exotic and one native O. taurus populations after three generations in common garden conditions. We sequenced RNA from 24 individuals for each of the three populations including both sexes, and spanning four developmental stages for the two exotic, differentiated populations. Identifying 270,400 high-quality single nucleotide polymorphisms, we revealed a strong signal of genetic differentiation between the three populations, and evidence of recent bottlenecks within and an excess of outlier loci between exotic populations. Differences in gene expression between populations were greatest in prepupae and early adult life stages, stages during which differences in male horn development and female fecundity manifest. Finally, genes differentially expressed between exotic populations also had greater genetic differentiation and performed functions related to chitin biosynthesis and nutrient sensing, possibly underlying allometry and fecundity trait divergences. Our results suggest that beyond bottlenecks, recent introductions have led to genetic and transcriptional differences in genes correlated with observed phenotypic differences.

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.

Molecular evolution of cytochrome C oxidase-I protein of insects living in Saudi Arabia

The study underpins barcode characterization of insect species collected from Saudi Arabia and explored functional constraints during evolution at the DNA and protein levels to expect the possible mechanisms of protein evolution in insects. Codon structure designated AT-biased insect barcode of the cytochrome C oxidase I (COI). In addition, the predicted 3D structure of COI protein indicated tyrosine in close proximity with the heme ligand, depicted substitution to phenylalanine in two Hymenopteran species. This change resulted in the loss of chemical bonding with the heme ligand. The estimated nucleotide substitution matrices in insect COI barcode generally showed a higher probability of transversion compared with the transition. Computations of codon-by-codon nonsynonymous substitutions in Hymenopteran and Hemipteran species indicated that almost half of the codons are under positive evolution. Nevertheless, codons of COI barcode of Coleoptera, Lepidoptera and Diptera are mostly under purifying selection. The results reinforce that codons in helices 2, 5 and 6 and those in loops 2–3 and 5–6 are mostly conserved and approach strong purifying selection. The overall results argue the possible evolutionary position of Hymenopteran species among those of other insects.

Morph-specific artificial selection reveals a constraint on the evolution of polyphenisms

Theory predicts that the evolution of polyphenic variation is facilitated where morphs are genetically uncoupled and free to evolve towards their phenotypic optima. However, the assumption that developmentally plastic morphs can evolve independently has not been tested directly. Using morph-specific artificial selection, we investigated correlated evolution between the sexes and male morphs of the bulb mite Rhizoglyphus echinopus Large 'fighter' males have a thick and sharply terminating pair of legs used to kill rival males, while small 'scrambler' males have unmodified legs, and search for unguarded females, avoiding fights. We selected on the relative leg width of only the fighter male morph, tracked the evolutionary responses in fighters and the correlated evolutionary responses in scramblers and females that were untouched by direct selection. Fighters diverged in relative leg thickness after six generations; assaying scramblers and females at the ninth generation we observed correlated responses in relative leg width in both. Our results represent strong evidence for the evolution of intraspecific phenotypic diversity despite correlated evolution between morphs and sexes, challenging the idea that male morphs are genetically uncoupled and free to independently respond to selection. We therefore question the perceived necessity for genetic independence in traits with extreme phenotypic plasticity.

On the Reciprocally Causal and Constructive Nature of Developmental Plasticity and Robustness

Exposure to environmental variation is a characteristic feature of normal development, one that organisms can respond to during their lifetimes by actively adjusting or maintaining their phenotype in order to maximize fitness. Plasticity and robustness have historically been studied by evolutionary biologists through quantitative genetic and reaction norm approaches, while more recent efforts emerging from evolutionary developmental biology have begun to characterize the molecular and developmental genetic underpinnings of both plastic and robust trait formation. In this review, we explore how our growing mechanistic understanding of plasticity and robustness is beginning to force a revision of our perception of both phenomena, away from our conventional view of plasticity and robustness as opposites along a continuum and toward a framework that emphasizes their reciprocal, constructive, and integrative nature. We do so in three sections. Following an introduction, the first section looks inward and reviews the genetic, epigenetic, and developmental mechanisms that enable organisms to sense and respond to environmental conditions, maintaining and adjusting trait formation in the process. In the second section, we change perspective and look outward, exploring the ways in which organisms reciprocally shape their environments in ways that influence trait formation, and do so through the lens of behavioral plasticity, niche construction, and host-microbiota interactions. In the final section, we revisit established plasticity and robustness concepts in light of these findings, and highlight research opportunities to further advance our understanding of the causes, mechanisms, and consequences of these ubiquitous, and interrelated, phenomena.

Sexual dimorphism and heightened conditional expression in a sexually selected weapon in the Asian rhinoceros beetle

Among the most dramatic examples of sexual selection are the weapons used in battles between rival males over access to females. As with ornaments of female choice, the most "exaggerated" sexually selected weapons vary from male to male more widely than other body parts (hypervariability), and their growth tends to be more sensitive to nutritional state or physiological condition compared with growth of other body parts ("heightened" conditional expression). Here, we use RNAseq analysis to build on recent work exploring these mechanisms in the exaggerated weapons of beetles, by examining patterns of differential gene expression in exaggerated (head and thorax horns) and non-exaggerated (wings, genitalia) traits in the Asian rhinoceros beetle, Trypoxylus dichotomus. Our results suggest that sexually dimorphic expression of weaponry involves large-scale changes in gene expression, relative to other traits, while nutrition-driven changes in gene expression in these same weapons are less pronounced. However, although fewer genes overall were differentially expressed in high- vs. low-nutrition individuals, the number of differentially expressed genes varied predictably according to a trait's degree of condition dependence (head horn > thorax horn > wings > genitalia). Finally, we observed a high degree of similarity in direction of effects (vectors) for subsets of differentially expressed genes across both sexually dimorphic and nutritionally responsive growth. Our results are consistent with a common set of mechanisms governing sexual size dimorphism and condition dependence.

Mechanisms of Plastic Rescue in Novel Environments

Adaptive phenotypic plasticity provides a mechanism of developmental rescue in novel and rapidly changing environments. Understanding the underlying mechanism of plasticity is important for predicting both the likelihood that a developmental response is adaptive and associated life-history trade-offs that could influence patterns of subsequent evolutionary rescue. Although evolved developmental switches may move organisms toward a new adaptive peak in a novel environment, such mechanisms often result in maladaptive responses. The induction of generalized physiological mechanisms in new environments is relatively more likely to result in adaptive responses to factors such as novel toxins, heat stress, or pathogens. Developmental selection forms of plasticity, which rely on within-individual selective processes, such as shaping of tissue architecture, trial-and-error learning, or acquired immunity, are particularly likely to result in adaptive plasticity in a novel environment. However, both the induction of plastic responses and the ability to be plastic through developmental selection come with significant costs, resulting in delays in reproduction, increased individual investment, and reduced fecundity. Thus, we might expect complex interactions between plastic responses that allow survival in novel environments and subsequent evolutionary responses at the population level.

Masculinization of gene expression is associated with male quality in Drosophila melanogaster

The signature of sexual selection has been revealed through the study of differences in patterns of genome-wide gene expression, both between the sexes and between alternative reproductive morphs within a single sex. What remains unclear, however, is whether differences in gene expression patterns between individuals of a given sex consistently map to variation in individual quality. Such a pattern, particularly if found in males, would provide unambiguous evidence that the phenotypic response to sexual selection is shaped through sex-specific alterations to the transcriptome. To redress this knowledge gap, we explored whether patterns of sex-biased gene expression are associated with variation in male reproductive quality in Drosophila melanogaster. We measured two male reproductive phenotypes, and their association with sex-biased gene expression, across a selection of inbred lines from the Drosophila Genetic Reference Panel. Genotypes with higher expression of male-biased genes produced males exhibiting shorter latencies to copulation, and higher capacity to inseminate females. Conversely, female-biased genes tended to show negative associations with these male reproductive traits across genotypes. We uncovered similar patterns, by reanalyzing a published dataset from a second D. melanogaster population. Our results reveal the footprint of sexual selection in masculinising the male transcriptome.

Recent advances in understanding the mechanisms of sexually dimorphic plasticity: insights from beetle weapons and future directions

Many traits that are sexually dimorphic, appearing either differently or uniquely in one sex, are also sensitive to an organism's condition. This phenomenon seems to have evolved to limit genetic conflict between traits that are under different selective pressures in each sex. Recent work has shed light on the molecular and developmental mechanisms that govern this condition sensitive growth, and this work has now expanded to encompass both sexual dimorphism as well as conditionally plastic growth, as it seems the two phenomena are linked on a molecular level. In all cases studied the gene doublesex, a conserved regulator of sex differentiation, controls both sexual dimorphism as well as the condition-dependent plastic responses common to these traits. However, the advent of next-generation -omics technologies has allowed researchers to decipher the common and diverged mechanisms of sexually dimorphic plasticity and expand investigations beyond the foundation laid by studies utilizing beetle weapons.

The distinct roles of insulin signaling in polyphenic development

Many insects have the ability to develop alternative morphologies in response to specific environmental signals such as photoperiod, temperature, nutrition and crowding. These signals are integrated by the brain and result in alternative patterns of secretion of developmental hormones like ecdysone, juvenile hormone and insulin-like growth factors, which, in turn, direct alternative developmental trajectories. Insulin signaling appears to be particularly important when the polyphenism involves differences in the sizes of the body, appendages and other structures, such as wings, mandibles and horns. Here we review recent advances in understanding the role of insulin signaling, and its interaction with other hormones, in the development of polyphenisms.

Temporal variation in brain transcriptome is associated with the expression of female mimicry as a sequential male alternative reproductive tactic in fish

Distinct patterns of gene expression often underlie intra- and intersexual differences, and the study of this set of coregulated genes is essential to understand the emergence of complex behavioural phenotypes. Here, we describe the development of a de novo transcriptome and brain gene expression profiles of wild-caught peacock blenny, Salaria pavo, an intertidal fish with sex-role reversal in courtship behaviour (i.e., females are the courting sex) and sequential alternative reproductive tactics in males (i.e., larger and older nest-holder males and smaller and younger sneaker males occur). Sneakers mimic both female's courtship behaviour and nuptial coloration to get access to nests and sneak fertilizations, and later in life transition into nest-holder males. Thus, this species offers the unique opportunity to study how the regulation of gene expression can contribute to intersex phenotypes and to the sequential expression of male and female behavioural phenotypes by the same individual. We found that at the whole brain level, expression of the sneaker tactic was paralleled by broader and divergent gene expression when compared to either females or nest-holder males, which were more similar between themselves. When looking at sex-biased transcripts, sneaker males are intersex rather than being either nest-holder or female-like, and their transcriptome is simultaneously demasculinized for nest-holder-biased transcripts and feminized for female-biased transcripts. These results indicate that evolutionary changes in reproductive plasticity can be achieved through regulation of gene expression, and in particular by varying the magnitude of expression of sex-biased genes, throughout the lifetime of the same individual.

Accelerated Evolution of Developmentally Biased Genes in the Tetraphenic Ant Cardiocondyla obscurior

Plastic gene expression underlies phenotypic plasticity and plastically expressed genes evolve under different selection regimes compared with ubiquitously expressed genes. Social insects are well-suited models to elucidate the evolutionary dynamics of plastic genes for their genetically and environmentally induced discrete polymorphisms. Here, we study the evolution of plastically expressed genes in the ant Cardiocondyla obscurior—a species that produces two discrete male morphs in addition to the typical female polymorphism of workers and queens. Based on individual-level gene expression data from 28 early third instar larvae, we test whether the same evolutionary dynamics that pertain to plastically expressed genes in adults also pertain to genes with plastic expression during development. In order to quantify plasticity of gene expression over multiple contrasts, we develop a novel geometric measure. For genes expressed during development, we show that plasticity of expression is positively correlated with evolutionary rates. We furthermore find a strong correlation between expression plasticity and expression variation within morphs, suggesting a close link between active and passive plasticity of gene expression. Our results support the notion of relaxed selection and neutral processes as important drivers in the evolution of adaptive plasticity.

Signals of selection in conditionally expressed genes in the diversification of three horned beetle species

Species radiations may be facilitated by phenotypic differences already present within populations, such as those arising through sex-specific development or developmental processes biased towards particular reproductive or trophic morphs. We sought to test this hypothesis by utilizing a comparative transcriptomic approach to contrast among- and within-species differentiation using three horned beetle species in the genus Onthophagus. These three species exhibit differences along three phenotypic axes reflective of much of the interspecific diversity present within the genus: horn location, polarity of sexual dimorphism and degree of nutritional sensitivity. Our approach combined de novo transcript assembly, assessment of amino acid substitutions (dN/dS) across orthologous gene pairs and integration of gene function and conditional gene expression data. We identified 17 genes across the three species pairs related to axis patterning, development and metabolism with dN/dS > 1 and detected elevated dN/dS in genes related to metabolism and biosynthesis in the most closely related species pair, which is characterized by a loss of nutritional polyphenism and a reversal of sexual dimorphism. Further, we found that genes that are conditionally expressed (i.e. as a function of sex, nutrition or body region) within one of our focal species also showed significantly stronger signals of positive or relaxed purifying selection between species divergent along the same morphological axis (i.e. polarity of sexual dimorphism, degree of nutritional sensitivity or location of horns). Our findings thus reveal a positive relationship between intraspecific differentiation due to condition-specific development and genetic divergences among species.

Mixed evidence for the erosion of intertactical genetic correlations through intralocus tactical conflict

Alternative reproductive tactics, whereby members of the same sex use different tactics to secure matings, are often associated with conditional intrasexual dimorphisms. Given the different selective pressures on males adopting each mating tactic, intrasexual dimorphism is more likely to arise if phenotypes are genetically uncoupled and free to evolve towards their phenotypic optima. However, in this context, genetic correlations between male morphs could result in intralocus tactical conflict (ITC). We investigated the genetic architecture of male dimorphism in bulb mites (Rhizoglyphus echinopus) and earwigs (Forficula auricularia). We used half-sibling breeding designs to assess the heritability and intra/intersexual genetic correlations of dimorphic and monomorphic traits in each species. We found two contrasting patterns; F. auricularia exhibited low intrasexual genetic correlations for the dimorphic trait, suggesting that the ITC is moving towards a resolution. Meanwhile, R. echinopus exhibited high and significant intrasexual genetic correlations for most traits, suggesting that morphs in the bulb mite may be limited in evolving to their optima. This also shows that intrasexual dimorphisms can evolve despite strong genetic constraints, contrary to current predictions. We discuss the implications of this genetic constraint and emphasize the potential importance of ITC for our understanding of intrasexual dimorphisms.

Evolutionary constraints shape caste-specific gene expression across 15 ant species

Development of polymorphic phenotypes from similar genomes requires gene expression differences. However, little is known about how morph-specific gene expression patterns vary on a broad phylogenetic scale. We hypothesize that evolution of morph-specific gene expression, and consequently morph-specific phenotypic evolution, may be constrained by gene essentiality and the amount of pleiotropic constraints. Here, we use comparative transcriptomics of queen and worker morphs, that is, castes, from 15 ant species to understand the constraints of morph-biased gene expression. In particular, we investigate how measures of evolutionary constraints at the sequence level (expression level, connectivity, and number of gene ontology [GO] terms) correlate with morph-biased expression. Our results show that genes indeed vary in their potential to become morph-biased. The existence of genes that are constrained in becoming caste-biased potentially limits the evolutionary decoupling of the caste phenotypes, that is, it might result in "caste load" occasioning from antagonistic fitness variation, similarly to sexually antagonistic fitness variation between males and females. On the other hand, we suggest that genes under low constraints are released from antagonistic variation and thus more likely to be co-opted for morph specific use. Overall, our results suggest that the factors that affect sequence evolutionary rates and evolution of plastic expression may largely overlap.

Asymmetric interactions between doublesex and tissue- and sex-specific target genes mediate sexual dimorphism in beetles

Sexual dimorphisms fuel significant intraspecific variation and evolutionary diversification. Yet the developmental-genetic mechanisms underlying sex-specific development remain poorly understood. Here, we focus on the conserved sex-determination gene doublesex (dsx) and the mechanisms by which it mediates sex-specific development in a horned beetle species by combining systemic dsx knockdown, high-throughput sequencing of diverse tissues and a genome-wide analysis of Dsx-binding sites. We find that Dsx regulates sex-biased expression predominantly in males, that Dsx's target repertoires are highly sex- and tissue-specific and that Dsx can exercise its regulatory role via two distinct mechanisms: as a sex-specific modulator by regulating strictly sex-specific targets, or as a switch by regulating the same genes in males and females in opposite directions. More generally, our results suggest Dsx can rapidly acquire new target gene repertoires to accommodate evolutionarily novel traits, evidenced by the large and unique repertoire identified in head horns, a recent morphological innovation.

Caste-biased gene expression in a facultatively eusocial bee suggests a role for genetic accommodation in the evolution of eusociality

Developmental plasticity may accelerate the evolution of phenotypic novelty through genetic accommodation, but studies of genetic accommodation often lack knowledge of the ancestral state to place selected traits in an evolutionary context. A promising approach for assessing genetic accommodation involves using a comparative framework to ask whether ancestral plasticity is related to the evolution of a particular trait. Bees are an excellent group for such comparisons because caste-based societies (eusociality) have evolved multiple times independently and extant species exhibit different modes of eusociality. We measured brain and abdominal gene expression in a facultatively eusocial bee, Megalopta genalis, and assessed whether plasticity in this species is functionally linked to eusocial traits in other bee lineages. Caste-biased abdominal genes in M. genalis overlapped significantly with caste-biased genes in obligately eusocial bees. Moreover, caste-biased genes in M. genalis overlapped significantly with genes shown to be rapidly evolving in multiple studies of 10 bee species, particularly for genes in the glycolysis pathway and other genes involved in metabolism. These results provide support for the idea that eusociality can evolve via genetic accommodation, with plasticity in facultatively eusocial species like M. genalis providing a substrate for selection during the evolution of caste in obligately eusocial lineages.

Transcriptomics of Intralocus Sexual Conflict: Gene Expression Patterns in Females Change in Response to Selection on a Male Secondary Sexual Trait in the Bulb Mite

Intralocus sexual conflict (IASC) prevents males and females from reaching their disparate phenotypic optima and is widespread, but little is known about its genetic underpinnings. In Rhizoglyphus robini, a mite species with alternative male morphs, elevated sexual dimorphism of the armored fighter males (compared to more feminized scramblers males) was previously reported to be associated with increased IASC. Because IASC persists if gene expression patterns are correlated between sexes, we compared gene expression patterns of males and females from the replicate lines selected for increased proportion of fighter or scrambler males (F- and S-lines, respectively). Specifically, we tested the prediction that selection for fighter morph caused correlated changes in gene expression patterns in females. We identified 532 differentially expressed genes (FDR < 0.05) between the F-line and S-line males. Consistent with the prediction, expression levels of these genes also differed between females from respective lines. Thus, significant proportion of genes differentially expressed between sexually selected male phenotypes showed correlated expression levels in females, likely contributing to elevated IASC in F-lines reported in a previous study.

Gene networks and developmental context: the importance of understanding complex gene expression patterns in evolution

Animal development is the product of distinct components and interactions-genes, regulatory networks, and cells-and it exhibits emergent properties that cannot be inferred from the components in isolation. Often the focus is on the genotype-to-phenotype map, overlooking the process of development that turns one into the other. We propose a move toward micro-evolutionary analysis of development, incorporating new tools that enable cell type resolution and single-cell microscopy. Using the sex determination pathway in Drosophila to illustrate potential avenues of research, we highlight some of the questions that these emerging technologies can address. For example, they provide an unprecedented opportunity to study heterogeneity within cell populations, and the potential to add the dimension of time to gene regulatory network analysis. Challenges still remain in developing methods to analyze this data and to increase the throughput. However this line of research has the potential to bridge the gaps between previously more disparate fields, such as population genetics and development, opening up new avenues of research.

Hedgehog signaling enables nutrition-responsive inhibition of an alternative morph in a polyphenic beetle

The recruitment of modular developmental genetic components into new developmental contexts has been proposed as a central mechanism enabling the origin of novel traits and trait functions without necessitating the origin of novel pathways. Here, we investigate the function of the hedgehog (Hh) signaling pathway, a highly conserved pathway best understood for its role in patterning anterior/posterior (A/P) polarity of diverse traits, in the developmental evolution of beetle horns, an evolutionary novelty, and horn polyphenisms, a highly derived form of environment-responsive trait induction. We show that interactions among pathway members are conserved during development of Onthophagus horned beetles and have retained the ability to regulate A/P polarity in traditional appendages, such as legs. At the same time, the Hh signaling pathway has acquired a novel and highly unusual role in the nutrition-dependent regulation of horn polyphenisms by actively suppressing horn formation in low-nutrition males. Down-regulation of Hh signaling lifts this inhibition and returns a highly derived sigmoid horn body size allometry to its presumed ancestral, linear state. Our results suggest that recruitment of the Hh signaling pathway may have been a key step in the evolution of trait thresholds, such as those involved in horn polyphenisms and the corresponding origin of alternative phenotypes and complex allometries.

De novo transcriptome assembly of the grapevine phylloxera allows identification of genes differentially expressed between leaf- and root-feeding forms

Background

Grapevine phylloxera, an insect related to true aphids, is a major historic pest of viticulture only controlled through the selection of resistant rootstocks or through quarantine regulations where grapevine is cultivated own-rooted. Transcriptomic data could help understand the bases of its original life-traits, including a striking case of polyphenism, with forms feeding on roots and forms feeding in leaf-galls. Comparisons with true aphids (for which complete genomes have been sequenced) should also allow to link differences in life-traits of the two groups with changes in gene repertoires or shifts in patterns of expression.

Results

We sequenced transcriptomes of the grapevine phylloxera (Illumina technology), choosing three life-stages (adults on roots or on leaf galls, and eggs) to cover a large catalogue of transcripts, and performed a de novo assembly. This resulted in 105,697 contigs, which were annotated: most contigs had a best blastx hit to the pea aphid (phylogenetically closest complete genome), while very few bacterial hits were recorded (except for Probionibacterium acnes). Coding sequences were predicted from this data set (17,372 sequences), revealing an extremely high AT-bias (at the third codon position). Differential expression (DE) analysis among root-feeding and gall-feeding showed that i) the root-feeding form displayed a much larger number of differentially expressed transcripts ii) root-feeding biased genes were enriched in some categories, for example cuticular proteins and genes associated with cell-cell signaling iii) leaf-galling-biased genes were enriched in genes associated with the nucleus and DNA-replication, suggesting a metabolism more oriented towards fast and active multiplication. We also identified a gene family with a very high expression level (copies totaling nearly 10 % of the reads) in the grapevine phylloxera (both in root and leaf galling forms), but usually expressed at very low levels in true aphids (except in sexual oviparous females). These transcripts thus appear to be associated with oviparity.

Conclusions

Our study illustrated major intraspecific changes in transcriptome profiles, related with different life-styles (and the feeding on roots versus in leaf-galls). At a different scale, we could also illustrate one major shift in expression levels associated with changes in life-traits that occurred along evolution and that respectively characterize (strictly oviparous) grapevine phylloxera and (mostly viviparous) true aphids.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2530-8) contains supplementary material, which is available to authorized users.

The transcriptomic basis of tissue- and nutrition-dependent sexual dimorphism in the beetle Onthophagus taurus

Sexual dimorphism accounts for a large fraction of intraspecific diversity. However, not all traits are equally sexually dimorphic; instead, individuals are mosaics of tissues that vary in their ability to exhibit dimorphism. Furthermore, the degree of a trait's sexual dimorphism is frequently environment‐dependent, with elaborate sexual dimorphism commonly being restricted to high nutritional conditions. Understanding the developmental basis and evolution of condition‐dependent sexual dimorphism can be critically informed by determining – across tissues and nutritional conditions – what sex‐biased genes are deployed and how they interact and translate into functional processes. Indeed, key theories concerning the evolution of condition‐dependent sexually dimorphic traits rest on assumptions regarding their developmental genetic underpinnings, yet, have largely gone unexamined by empirical studies. Here, we provide such evidence by investigating the transcriptomic basis of tissue‐ and nutrition‐dependent sexual dimorphism in the bull‐headed dung beetle Onthophagus taurus. Our findings suggest (1) that generating morphological sexual dimorphism requires sex‐biased gene expression in and developmental remodeling of both sexes, regardless of which sex exhibits externally visible trait exaggeration, (2) that although sexually dimorphic phenotypes are comprised of traits underlain by independent repertoires of sex‐biased gene expression, they act similarly at a functional level, and (3) that sexual dimorphism and condition‐dependence share common genetic underpinnings specifically in sexually‐selected traits.

Constraints on the evolution of phenotypic plasticity: limits and costs of phenotype and plasticity

Phenotypic plasticity is ubiquitous and generally regarded as a key mechanism for enabling organisms to survive in the face of environmental change. Because no organism is infinitely or ideally plastic, theory suggests that there must be limits (for example, the lack of ability to produce an optimal trait) to the evolution of phenotypic plasticity, or that plasticity may have inherent significant costs. Yet numerous experimental studies have not detected widespread costs. Explicitly differentiating plasticity costs from phenotype costs, we re-evaluate fundamental questions of the limits to the evolution of plasticity and of generalists vs specialists. We advocate for the view that relaxed selection and variable selection intensities are likely more important constraints to the evolution of plasticity than the costs of plasticity. Some forms of plasticity, such as learning, may be inherently costly. In addition, we examine opportunities to offset costs of phenotypes through ontogeny, amelioration of phenotypic costs across environments, and the condition-dependent hypothesis. We propose avenues of further inquiry in the limits of plasticity using new and classic methods of ecological parameterization, phylogenetics and omics in the context of answering questions on the constraints of plasticity. Given plasticity's key role in coping with environmental change, approaches spanning the spectrum from applied to basic will greatly enrich our understanding of the evolution of plasticity and resolve our understanding of limits.

Eco-Evo-Devo: developmental symbiosis and developmental plasticity as evolutionary agents

The integration of research from developmental biology and ecology into evolutionary theory has given rise to a relatively new field, ecological evolutionary developmental biology (Eco-Evo-Devo). This field integrates and organizes concepts such as developmental symbiosis, developmental plasticity, genetic accommodation, extragenic inheritance and niche construction. This Review highlights the roles that developmental symbiosis and developmental plasticity have in evolution. Developmental symbiosis can generate particular organs, can produce selectable genetic variation for the entire animal, can provide mechanisms for reproductive isolation, and may have facilitated evolutionary transitions. Developmental plasticity is crucial for generating novel phenotypes, facilitating evolutionary transitions and altered ecosystem dynamics, and promoting adaptive variation through genetic accommodation and niche construction. In emphasizing such non-genomic mechanisms of selectable and heritable variation, Eco-Evo-Devo presents a new layer of evolutionary synthesis.

The nutritionally responsive transcriptome of the polyphenic beetle Onthophagus taurus and the importance of sexual dimorphism and body region

Developmental responses to nutritional variation represent one of the ecologically most important classes of adaptive plasticity. However, knowledge of genome-wide patterns of nutrition-responsive gene expression is limited. Here, we studied genome-wide transcriptional responses to nutritional variation and their dependency on trait and sex in the beetle Onthophagus taurus. We find that averaged across the transcriptome, nutrition contributes less to overall variation in gene expression than do sex or body region, but that for a modest subset of genes nutrition is by far the most important determinant of expression variation. Furthermore, our results reject the hypothesis that a common machinery may underlie nutrition-sensitive development across body regions. Instead, we find that magnitude (measured by number of differentially expressed contigs), composition (measured by functional enrichment) and evolutionary consequences (measured by patterns of sequence variation) are heavily dependent on exactly which body region is considered and the degree of sexual dimorphism observed on a morphological level. More generally, our findings illustrate that studies into the developmental mechanisms and evolutionary consequences of nutrition-biased gene expression must take into account the dynamics and complexities imposed by other sources of variation in gene expression such as sexual dimorphism and trait type.

The locus of sexual selection: moving sexual selection studies into the post-genomics era

Sexual selection drives fundamental evolutionary processes such as trait elaboration and speciation. Despite this importance, there are surprisingly few examples of genes unequivocally responsible for variation in sexually selected phenotypes. This lack of information inhibits our ability to predict phenotypic change due to universal behaviours, such as fighting over mates and mate choice. Here, we discuss reasons for this apparent gap and provide recommendations for how it can be overcome by adopting contemporary genomic methods, exploiting underutilized taxa that may be ideal for detecting the effects of sexual selection and adopting appropriate experimental paradigms. Identifying genes that determine variation in sexually selected traits has the potential to improve theoretical models and reveal whether the genetic changes underlying phenotypic novelty utilize common or unique molecular mechanisms. Such a genomic approach to sexual selection will help answer questions in the evolution of sexually selected phenotypes that were first asked by Darwin and can furthermore serve as a model for the application of genomics in all areas of evolutionary biology.

Intralocus tactical conflict: genetic correlations between fighters and sneakers of the dung beetle Onthophagus taurus

Males and females differ in their phenotypic optima for many traits, and as the majority of genes are expressed in both sexes, some alleles can be beneficial to one sex but harmful to the other (intralocus sexual conflict; ISC). ISC theory has recently been extended to intrasexual dimorphisms, where certain alleles may have opposite effects on the fitness of males of different morphs that employ alternative reproductive tactics (intralocus tactical conflict; ITC). Here, we use a half-sib breeding design to investigate the genetic basis for ISC and ITC in the dung beetle Onthophagus taurus. We found positive heritabilities and intersexual genetic correlations for almost all traits investigated. Next, we calculated the intrasexual genetic correlation between males of different morphs for horn length, a sexually selected trait, and compared it to intrasexual correlations for naturally selected traits in both sexes. Intrasexual genetic correlations did not differ significantly between the sexes or between naturally and sexually selected traits, failing to support the hypothesis that horns present a reduction of intrasexual genetic correlations due to ITC. We discuss the implications for the idea of developmental reprogramming between male morphs and emphasize the importance of genetic correlations as constraints for the evolution of dimorphisms.

Neutral and adaptive explanations for an association between caste-biased gene expression and rate of sequence evolution

The castes of social insects provide outstanding opportunities to address the causes and consequences of evolution of discrete phenotypes, i.e., polymorphisms. Here we focus on recently described patterns of a positive association between the degree of caste-specific gene expression and the rate of sequence evolution. We outline how neutral and adaptive evolution can cause genes that are morph-biased in their expression profiles to exhibit historical signatures of faster or slower sequence evolution compared to unbiased genes. We conclude that evaluation of different hypotheses will benefit from (i) reconstruction of the phylogenetic origin of biased expression and changes in rates of sequence evolution, and (ii) replicated data on gene expression variation within versus between morphs. Although the data are limited at present, we suggest that the observed phylogenetic and intra-population variation in gene expression lends support to the hypothesis that the association between caste-biased expression and rate of sequence evolution largely is a result of neutral processes.

Development and evolution of insect polyphenisms: novel insights through the study of sex determination mechanisms

Polyphenism is defined as a single individual's ability to develop into two or more alternative phenotypes. Polyphenic development is taxonomically widespread among insects, cued by diverse environmental factors, and enables single genotypes to accommodate breathtaking phenotypic diversity. Most research on the developmental regulation and evolution of insect polyphenisms has focused on endocrine control mechanisms, in particular the role of the sesquiterpenoid juvenile hormone. Here we review recent finding that suggest additional and previously overlooked mechanisms that underlie the developmental regulation and rapid evolution of polyphenic development. Specifically, we focus on the role of somatic sex determination mechanisms, which mediate body-region and tissue-specific differentiation as a function of sex across insects. Recent work on Onthophagus horned beetles suggest that the same mechanisms have been co-opted to regulate the development of nutritionally cued, alternative male morphs, and that rapid changes in these mechanisms underlie rapid evolutionary changes in patterns of phenotype expression across Onthophagus species, including the loss of old and gain of novel locations for horn development, the evolution of reversed sexual dimorphisms, and the secondary loss of male polyphenism. We discuss how these findings expand the way we think about the origins and diversification of polyphenisms, and close by briefly highlighting potentially fruitful avenues for future research.

Accelerated evolution of morph-biased genes in pea aphids

Phenotypic plasticity, the production of alternative phenotypes (or morphs) from the same genotype due to environmental factors, results in some genes being expressed in a morph-biased manner. Theoretically, these morph-biased genes experience relaxed selection, the consequence of which is the buildup of slightly deleterious mutations at these genes. Over time, this is expected to result in increased protein divergence at these genes between species and a signature of relaxed purifying selection within species. Here we test these theoretical expectations using morph-biased genes in the pea aphid, a species that produces multiple morphs via polyphenism. We find that morph-biased genes exhibit faster rates of evolution (in terms of dN/dS) relative to unbiased genes and that divergence generally increases with increasing morph bias. Further, genes with expression biased toward rarer morphs (sexual females and males) show faster rates of evolution than genes expressed in the more common morph (asexual females), demonstrating that the amount of time a gene spends being expressed in a morph is associated with its rate of evolution. And finally, we show that genes expressed in the rarer morphs experience decreased purifying selection relative to unbiased genes, suggesting that it is a relaxation of purifying selection that contributes to their faster rates of evolution. Our results provide an important empirical look at the impact of phenotypic plasticity on gene evolution.

Phenotypic and genomic plasticity of alternative male reproductive tactics in sailfin mollies

A major goal of modern evolutionary biology is to understand the causes and consequences of phenotypic plasticity, the ability of a single genotype to produce multiple phenotypes in response to variable environments. While ecological and quantitative genetic studies have evaluated models of the evolution of adaptive plasticity, some long-standing questions about plasticity require more mechanistic approaches. Here, we address two of those questions: does plasticity facilitate adaptive evolution? And do physiological costs place limits on plasticity? We examine these questions by comparing genetically and plastically regulated behavioural variation in sailfin mollies (Poecilia latipinna), which exhibit striking variation in plasticity for male mating behaviour. In this species, some genotypes respond plastically to a change in the social environment by switching between primarily courting and primarily sneaking behaviour. In contrast, other genotypes have fixed mating strategies (either courting or sneaking) and do not display plasticity. We found that genetic and plastic variation in behaviour were accompanied by partially, but not completely overlapping changes in brain gene expression, in partial support of models that predict that plasticity can facilitate adaptive evolution. We also found that behavioural plasticity was accompanied by broader and more robust changes in brain gene expression, suggesting a substantial physiological cost to plasticity. We also observed that sneaking behaviour, but not courting, was associated with upregulation of genes involved in learning and memory, suggesting that sneaking is more cognitively demanding than courtship.

Insights into the development and evolution of exaggerated traits using de novo transcriptomes of two species of horned scarab beetles

Scarab beetles exhibit an astonishing variety of rigid exo-skeletal outgrowths, known as “horns”. These traits are often sexually dimorphic and vary dramatically across species in size, shape, location, and allometry with body size. In many species, the horn exhibits disproportionate growth resulting in an exaggerated allometric relationship with body size, as compared to other traits, such as wings, that grow proportionately with body size. Depending on the species, the smallest males either do not produce a horn at all, or they produce a disproportionately small horn for their body size. While the diversity of horn shapes and their behavioural ecology have been reasonably well studied, we know far less about the proximate mechanisms that regulate horn growth. Thus, using 454 pyrosequencing, we generated transcriptome profiles, during horn growth and development, in two different scarab beetle species: the Asian rhinoceros beetle, Trypoxylus dichotomus, and the dung beetle, Onthophagus nigriventris. We obtained over half a million reads for each species that were assembled into over 6,000 and 16,000 contigs respectively. We combined these data with previously published studies to look for signatures of molecular evolution. We found a small subset of genes with horn-biased expression showing evidence for recent positive selection, as is expected with sexual selection on horn size. We also found evidence of relaxed selection present in genes that demonstrated biased expression between horned and horn-less morphs, consistent with the theory of developmental decoupling of phenotypically plastic traits.

Gene expression patterns associated with caste and reproductive status in ants: worker-specific genes are more derived than queen-specific ones

Variation in gene expression leads to phenotypic diversity and plays a central role in caste differentiation of eusocial insect species. In social Hymenoptera, females with the same genetic background can develop into queens or workers, which are characterized by divergent morphologies, behaviours and lifespan. Moreover, many social insects exhibit behaviourally distinct worker castes, such as brood-tenders and foragers. Researchers have just started to explore which genes are differentially expressed to achieve this remarkable phenotypic plasticity. Although the queen is normally the only reproductive individual in the nest, following her removal, young brood-tending workers often develop ovaries and start to reproduce. Here, we make use of this ability in the ant Temnothorax longispinosus and compare gene expression patterns in the queens and three worker castes along a reproductive gradient. We found the largest expression differences between the queen and the worker castes (~2500 genes) and the smallest differences between infertile brood-tenders and foragers (~300 genes). The expression profile of fertile workers is more worker-like, but to a certain extent intermediate between the queen and the infertile worker castes. In contrast to the queen, a high number of differentially expressed genes in the worker castes are of unknown function, pointing to the derived status of hymenopteran workers within insects.

Transcriptome analyses of primitively eusocial wasps reveal novel insights into the evolution of sociality and the origin of alternative phenotypes

BACKGROUND

Understanding how alternative phenotypes arise from the same genome is a major challenge in modern biology. Eusociality in insects requires the evolution of two alternative phenotypes - workers, who sacrifice personal reproduction, and queens, who realize that reproduction. Extensive work on honeybees and ants has revealed the molecular basis of derived queen and worker phenotypes in highly eusocial lineages, but we lack equivalent deep-level analyses of wasps and of primitively eusocial species, the latter of which can reveal how phenotypic decoupling first occurs in the early stages of eusocial evolution.

RESULTS

We sequenced 20 Gbp of transcriptomes derived from brains of different behavioral castes of the primitively eusocial tropical paper wasp Polistes canadensis. Surprisingly, 75% of the 2,442 genes differentially expressed between phenotypes were novel, having no significant homology with described sequences. Moreover, 90% of these novel genes were significantly upregulated in workers relative to queens. Differential expression of novel genes in the early stages of sociality may be important in facilitating the evolution of worker behavioral complexity in eusocial evolution. We also found surprisingly low correlation in the identity and direction of expression of differentially expressed genes across similar phenotypes in different social lineages, supporting the idea that social evolution in different lineages requires substantial de novo rewiring of molecular pathways.

CONCLUSIONS

These genomic resources for aculeate wasps and first transcriptome-wide insights into the origin of castes bring us closer to a more general understanding of eusocial evolution and how phenotypic diversity arises from the same genome.