Caste‐biases in gene expression are specific to developmental stage in the ant Formica exsecta

Understanding how a single genome creates and maintains distinct phenotypes is a central goal in evolutionary biology. Social insects are a striking example of co‐opted genetic backgrounds giving rise to dramatically different phenotypes, such as queen and worker castes. A conserved set of molecular pathways, previously envisioned as a set of ‘toolkit’ genes, has been hypothesized to underlie queen and worker phenotypes in independently evolved social insect lineages. Here, we investigated the toolkit from a developmental point of view, using RNA‐Seq to compare caste‐biased gene expression patterns across three life stages (pupae, emerging adult and old adult) and two female castes (queens and workers) in the ant Formica exsecta. We found that the number of genes with caste‐biased expression increases dramatically from pupal to old adult stages. This result suggests that phenotypic differences between queens and workers at the pupal stage may derive from a relatively low number of caste‐biased genes, compared to higher number of genes required to maintain caste differences at the adult stage. Gene expression patterns were more similar among castes within developmental stages than within castes despite the extensive phenotypic differences between queens and workers. Caste‐biased expression was highly variable among life stages at the level of single genes, but more consistent when gene functions (gene ontology terms) were investigated. Finally, we found that a large part of putative toolkit genes were caste‐biased at least in some life stages in F. exsecta, and the caste‐biases, but not their direction, were more often shared between F. exsecta and other ant species than between F. exsecta and bees. Our results indicate that gene expression should be examined across several developmental stages to fully reveal the genetic basis of polyphenisms.     

Gene duplication and the evolution of phenotypic diversity in insect societies

Gene duplication is an important evolutionary process thought to facilitate the evolution of phenotypic diversity. We investigated if gene duplication was associated with the evolution of phenotypic differences in a highly social insect, the honeybee Apis mellifera. We hypothesized that the genetic redundancy provided by gene duplication could promote the evolution of social and sexual phenotypes associated with advanced societies. We found a positive correlation between sociality and rate of gene duplications across the Apoidea, indicating that gene duplication may be associated with sociality. We also discovered that genes showing biased expression between A. mellifera alternative phenotypes tended to be found more frequently than expected among duplicated genes than singletons. Moreover, duplicated genes had higher levels of caste‐, sex‐, behavior‐, and tissue‐biased expression compared to singletons, as expected if gene duplication facilitated phenotypic differentiation. We also found that duplicated genes were maintained in the A. mellifera genome through the processes of conservation, neofunctionalization, and specialization, but not subfunctionalization. Overall, we conclude that gene duplication may have facilitated the evolution of social and sexual phenotypes, as well as tissue differentiation. Thus this study further supports the idea that gene duplication allows species to evolve an increased range of phenotypic diversity.     

Differential gene expression and phenotypic plasticity in behavioural castes of the primitively eusocial wasp, Polistes canadensis

Understanding how a single genome can produce a variety of different phenotypes is of fundamental importance in evolutionary and developmental biology. One of the most striking examples of phenotypic plasticity is the female caste system found in eusocial insects, where variation in reproductive (queens) and non-reproductive (workers) phenotypes results in a broad spectrum of caste types, ranging from behavioural through to morphological castes. Recent advances in genomic techniques allow novel comparisons on the nature of caste phenotypes to be made at the level of the genes in organisms for which there is little genome information, facilitating new approaches in studying social evolution and behaviour. Using the paper wasp Polistes canadensis as a model system, we investigated for the first time how behavioural castes in primitively eusocial insect societies are associated with differential expression of shared genes. We found that queens and newly emerged females express gene expression patterns that are distinct from each other whilst workers generally expressed intermediate patterns, as predicted by Polistes biology. We compared caste-associated genes in P. canadensis with those expressed in adult queens and workers of more advanced eusocial societies, which represent four independent origins of eusociality. Nine genes were conserved across the four taxa, although their patterns of expression and putative functions varied. Thus, we identify several genes that are putatively of evolutionary importance in the molecular biology that underlies a number of caste systems of independent evolutionary origin.     

DEVELOPMENTAL DECOUPLING OF ALTERNATIVE PHENOTYPES: INSIGHTS FROM THE TRANSCRIPTOMES OF HORN‐POLYPHENIC BEETLES

Developmental mechanisms play an important role in determining the costs, limits, and evolutionary consequences of phenotypic plasticity. One issue central to these claims is the hypothesis of developmental decoupling, where alternate morphs result from evolutionarily independent developmental pathways. We address this assumption through a microarray study that tests whether differences in gene expression between alternate morphs are as divergent as those between sexes, a classic example of developmental decoupling. We then examine whether genes with morph‐biased expression are less conserved than genes with shared expression between morphs, as predicted if developmental decoupling relaxes pleiotropic constraints on divergence. We focus on the developing horns and brains of two species of horned beetles with impressive sexual‐ and morph‐dimorphism in the expression of horns and fighting behavior. We find that patterns of gene expression were as divergent between morphs as they were between sexes. However, overall patterns of gene expression were also highly correlated across morphs and sexes. Morph‐biased genes were more evolutionarily divergent, suggesting a role of relaxed pleiotropic constraints or relaxed selection. Together these results suggest that alternate morphs are to some extent developmentally decoupled, and that this decoupling has significant evolutionary consequences. However, alternative morphs may not be as developmentally decoupled as sometimes assumed and such hypotheses of development should be revisited and refined.     

Robustness to noise in gene expression evolves despite epistatic constraints in a model of gene networks

Stochastic noise in gene expression causes variation in the development of phenotypes, making such noise a potential target of stabilizing selection. Here, we develop a new simulation model of gene networks to study the adaptive landscape underlying the evolution of robustness to noise. We find that epistatic interactions between the determinants of the expression of a gene and its downstream effect impose significant constraints on evolution, but these interactions do allow the gradual evolution of increased robustness. Despite strong sign epistasis, adaptation rarely proceeds via deleterious intermediate steps, but instead occurs primarily through small beneficial mutations. A simple mathematical model captures the relevant features of the single‐gene fitness landscape and explains counterintuitive patterns, such as a correlation between the mean and standard deviation of phenotypes. In more complex networks, mutations in regulatory regions provide evolutionary pathways to increased robustness. These results chart the constraints and possibilities of adaptation to reduce expression noise and demonstrate the potential of a novel modeling framework for gene networks.     

The contribution of the mitochondrial genome to sex‐specific fitness variance

Maternal inheritance of mitochondrial DNA (mtDNA) facilitates the evolutionary accumulation of mutations with sex‐biased fitness effects. Whereas maternal inheritance closely aligns mtDNA evolution with natural selection in females, it makes it indifferent to evolutionary changes that exclusively benefit males. The constrained response of mtDNA to selection in males can lead to asymmetries in the relative contributions of mitochondrial genes to female versus male fitness variation. Here, we examine the impact of genetic drift and the distribution of fitness effects (DFE) among mutations—including the correlation of mutant fitness effects between the sexes—on mitochondrial genetic variation for fitness. We show how drift, genetic correlations, and skewness of the DFE determine the relative contributions of mitochondrial genes to male versus female fitness variance. When mutant fitness effects are weakly correlated between the sexes, and the effective population size is large, mitochondrial genes should contribute much more to male than to female fitness variance. In contrast, high fitness correlations and small population sizes tend to equalize the contributions of mitochondrial genes to female versus male variance. We discuss implications of these results for the evolution of mitochondrial genome diversity and the genetic architecture of female and male fitness.     

Soldier‐biased gene expression in a subterranean termite implies functional specialization of the defensive caste

In a termite colony, reproduction is typically monopolized by a small number of sexuals that are supported by reproductively altruistic soldiers and workers. We expect caste differentiation to be associated with clear‐cut differences in gene expression, and for these differences to reflect caste function and development. Here, we use RNA‐Sequencing to compare the gene expression profiles of sexual nymphs and two non‐reproductive helper castes (i.e., workers and soldiers) of the Eastern subterranean termite Reticulitermes flavipes. We found that of n = 93 genes that are strictly expressed as a function of caste, a majority (78%) show a soldier‐specific pattern. This conspicuous soldier‐bias in genome‐wide expression suggests that this defensively specialized caste is functionally well‐differentiated from both the reproductive and the other non‐reproductive caste of this species, despite a shared developmental program with workers. Gene ontology analysis supports the notion of functional specialization by soldiers, as soldier‐biased gene sets are enriched for novel biological processes. Whether this pattern reflects ancient or more recent bouts of selection for caste novelty at the gene‐regulatory level is not known, but because soldiers are sterile and thus have no direct fitness, any selection for novelty must have been mediated indirectly, through reproducing relatives.     

Relaxed genetic constraint is ancestral to the evolution of phenotypic plasticity

Phenotypic plasticity--the capacity of a single genotype to produce different phenotypes in response to varying environmental conditions--is widespread. Yet, whether, and how, plasticity impacts evolutionary diversification is unclear. According to a widely discussed hypothesis, plasticity promotes rapid evolution because genes expressed differentially across different environments (i.e., genes with "biased" expression) experience relaxed genetic constraint and thereby accumulate variation faster than do genes with unbiased expression. Indeed, empirical studies confirm that biased genes evolve faster than unbiased genes in the same genome. An alternative hypothesis holds, however, that the relaxed constraint and faster evolutionary rates of biased genes may be a precondition for, rather than a consequence of, plasticity's evolution. Here, we evaluated these alternative hypotheses by characterizing evolutionary rates of biased and unbiased genes in two species of frogs that exhibit a striking form of phenotypic plasticity. We also characterized orthologs of these genes in four species of frogs that had diverged from the two plastic species before the plasticity evolved. We found that the faster evolutionary rates of biased genes predated the evolution of the plasticity. Furthermore, biased genes showed greater expression variance than did unbiased genes, suggesting that they may be more dispensable. Phenotypic plasticity may therefore evolve when dispensable genes are co-opted for novel function in environmentally induced phenotypes. Thus, relaxed genetic constraint may be a cause--not a consequence--of the evolution of phenotypic plasticity, and thereby contribute to the evolution of novel traits.     

Synthetic circuits reveal how mechanisms of gene regulatory networks constrain evolution

Phenotypic variation is the raw material of adaptive Darwinian evolution. The phenotypic variation found in organismal development is biased towards certain phenotypes, but the molecular mechanisms behind such biases are still poorly understood. Gene regulatory networks have been proposed as one cause of constrained phenotypic variation. However, most pertinent evidence is theoretical rather than experimental. Here, we study evolutionary biases in two synthetic gene regulatory circuits expressed in Escherichia coli that produce a gene expression stripe—a pivotal pattern in embryonic development. The two parental circuits produce the same phenotype, but create it through different regulatory mechanisms. We show that mutations cause distinct novel phenotypes in the two networks and use a combination of experimental measurements, mathematical modelling and DNA sequencing to understand why mutations bring forth only some but not other novel gene expression phenotypes. Our results reveal that the regulatory mechanisms of networks restrict the possible phenotypic variation upon mutation. Consequently, seemingly equivalent networks can indeed be distinct in how they constrain the outcome of further evolution.     

The microevolutionary response to male‐limited X‐chromosome evolution in Drosophila melanogaster reflects macroevolutionary patterns

Due to its hemizygous inheritance and role in sex determination, the X‐chromosome is expected to play an important role in the evolution of sexual dimorphism and to be enriched for sexually antagonistic genetic variation. By forcing the X‐chromosome to only be expressed in males over >40 generations, we changed the selection pressures on the X to become similar to those experienced by the Y. This releases the X from any constraints arising from selection in females and should lead to specialization for male fitness, which could occur either via direct effects of X‐linked loci or trans‐regulation of autosomal loci by the X. We found evidence of masculinization via up‐regulation of male‐benefit sexually antagonistic genes and down‐regulation of X‐linked female‐benefit genes. Potential artefacts of the experimental evolution protocol are discussed and cannot be wholly discounted, leading to several caveats. Interestingly, we could detect evidence of microevolutionary changes consistent with previously documented macroevolutionary patterns, such as changes in expression consistent with previously established patterns of sexual dimorphism, an increase in the expression of metabolic genes related to mito‐nuclear conflict and evidence that dosage compensation effects can be rapidly altered. These results confirm the importance of the X in the evolution of sexual dimorphism and as a source for sexually antagonistic genetic variation and demonstrate that experimental evolution can be a fruitful method for testing theories of sex chromosome evolution.     

Delusions of grandeur: Seed count is not a good fitness proxy under individual variation in phenology

The concept of fitness is central to evolutionary biology, yet it is difficult to define and to measure. In plant biology, fitness is often measured as seed count. However, under an array of circumstances, seed count may be a biased proxy of fitness, for example when individuals vary in allocation to sexual versus asexual reproduction. A more subtle example, but also likely to be important in natural populations, is when interindividual variation in conditions during development results in variation in offspring quality among seed parents. In monocarpic (semelparous) plants, this is expected to result from variation in effective season length experienced among individuals that reach reproductive maturity at different times. Here, we manipulate growing season length to ask whether seed count is an accurate representation of parental fitness in the monocarpic herb Lobelia inflata. Simple seed count suggests a paradoxical fitness advantage under constrained‐season length. However, we find that the apparent fitness advantage of a constrained‐season length is overridden by low relative per‐seed fitness. Furthermore, the fitness deficit in the constrained environment is associated primarily with an accelerating decrease in viability and seedling survival in seeds derived from fruits produced progressively later in the season. In this study, the overall fitness value of a seed under a constrained season is 0.774 of that observed under a long season.     

Conserved roles of Osiris genes in insect development,polymorphism and protection

Much of the variation among insects is derived from the different ways that chitin has been moulded to form rigid structures, both internal and external. In this study, we identify a highly conserved expression pattern in an insect‐only gene family, the Osiris genes, that is essential for development, but also plays a significant role in phenotypic plasticity and in immunity/toxicity responses. The majority of Osiris genes exist in a highly syntenic cluster, and the cluster itself appears to have arisen very early in the evolution of insects. We used developmental gene expression in the fruit fly, Drosophila melanogaster, the bumble bee, Bombus terrestris, the harvester ant, Pogonomyrmex barbatus, and the wood ant, Formica exsecta, to compare patterns of Osiris gene expression both during development and between alternate caste phenotypes in the polymorphic social insects. Developmental gene expression of Osiris genes is highly conserved across species and correlated with gene location and evolutionary history. The social insect castes are highly divergent in pupal Osiris gene expression. Sets of co‐expressed genes that include Osiris genes are enriched in gene ontology terms related to chitin/cuticle and peptidase activity. Osiris genes are essential for cuticle formation in both embryos and pupae, and genes co‐expressed with Osiris genes affect wing development. Additionally, Osiris genes and those co‐expressed seem to play a conserved role in insect toxicology defences and digestion. Given their role in development, plasticity, and protection, we propose that the Osiris genes play a central role in insect adaptive evolution.     

Genetic accommodation in the wild: evolution of gene expression plasticity during character displacement

Ecological character displacement is considered crucial in promoting diversification, yet relatively little is known of its underlying mechanisms. We examined whether evolutionary shifts in gene expression plasticity (‘genetic accommodation’) mediate character displacement in spadefoot toads. Where Spea bombifrons and S. multiplicata occur separately in allopatry (the ancestral condition), each produces alternative, diet‐induced, larval ecomorphs: omnivores, which eat detritus, and carnivores, which specialize on shrimp. By contrast, where these two species occur together in sympatry (the derived condition), selection to minimize competition for detritus has caused S. bombifrons to become nearly fixed for producing only carnivores, suggesting that character displacement might have arisen through an extreme form of genetic accommodation (‘genetic assimilation’) in which plasticity is lost. Here, we asked whether we could infer a signature of this process in regulatory changes of specific genes. In particular, we investigated whether genes that are normally expressed more highly in one morph (‘biased’ genes) have evolved reduced plasticity in expression levels among S. bombifrons from sympatry compared to S. bombifrons from allopatry. We reared individuals from sympatry vs. allopatry on detritus or shrimp and measured the reaction norms of nine biased genes. Although different genes displayed different patterns of gene regulatory evolution, the combined gene expression profiles revealed that sympatric individuals had indeed lost the diet‐induced gene expression plasticity present in allopatric individuals. Our data therefore provide one of the few examples from natural populations in which genetic accommodation/assimilation can be traced to regulatory changes of specific genes. Such genetic accommodation might mediate character displacement in many systems.     

Sex‐biased transcriptome divergence along a latitudinal gradient

Sex‐dependent gene expression is likely an important genomic mechanism that allows sex‐specific adaptation to environmental changes. Among Drosophila species, sex‐biased genes display remarkably consistent evolutionary patterns; male‐biased genes evolve faster than unbiased genes in both coding sequence and expression level, suggesting sex differences in selection through time. However, comparatively little is known of the evolutionary process shaping sex‐biased expression within species. Latitudinal clines offer an opportunity to examine how changes in key ecological parameters also influence sex‐specific selection and the evolution of sex‐biased gene expression. We assayed male and female gene expression in Drosophila serrata along a latitudinal gradient in eastern Australia spanning most of its endemic distribution. Analysis of 11 631 genes across eight populations revealed strong sex differences in the frequency, mode and strength of divergence. Divergence was far stronger in males than females and while latitudinal clines were evident in both sexes, male divergence was often population specific, suggesting responses to localized selection pressures that do not covary predictably with latitude. While divergence was enriched for male‐biased genes, there was no overrepresentation of X‐linked genes in males. By contrast, X‐linked divergence was elevated in females, especially for female‐biased genes. Many genes that diverged in D. serrata have homologs also showing latitudinal divergence in Drosophila simulans and Drosophila melanogaster on other continents, likely indicating parallel adaptation in these distantly related species. Our results suggest that sex differences in selection play an important role in shaping the evolution of gene expression over macro‐ and micro‐ecological spatial scales.     

Time constraints related to sexual maturation and prolonged copulation in the female‐dimorphic damselfly Ischnura senegalensis

Time constraints are critical for reproductive success. To understand the spatiotemporal dynamics of morph frequency in the female‐dimorphic damselfly Ischnura senegalensis, we compared two different morphs for two important time constraints on female reproductive output, i.e. post‐emergence sexual maturation and prolonged copulation. The females of both morphs achieved sexual maturation 4–5 days after emergence, suggesting that the rate of sexual mutation does not result in morph‐specific fitness. The copulation durations declined with the time of onset of copulation in both morphs. Consequently, all copulations terminated at approximately 12:00 hours. Because females show foraging and oviposition activity only after copulation, the copulation duration does not result in morph‐specific time constraints. These two important time constraints do not account for morph‐specific reproductive success and do not affect the evolutionary equilibrium of morph frequency in I. senegalensis.     

Building a new research framework for social evolution: intralocus caste antagonism

The breeding and non‐breeding ‘castes’ of eusocial insects provide a striking example of role‐specific selection, where each caste maximises fitness through different morphological, behavioural and physiological trait values. Typically, queens are long‐lived egg‐layers, while workers are short‐lived, largely sterile foragers. Remarkably, the two castes are nevertheless produced by the same genome. The existence of inter‐caste genetic correlations is a neglected consequence of this shared genome, potentially hindering the evolution of caste dimorphism: alleles that increase the productivity of queens may decrease the productivity of workers and vice versa, such that each caste is prevented from reaching optimal trait values. A likely consequence of this ‘intralocus caste antagonism’ should be the maintenance of genetic variation for fitness and maladaptation within castes (termed ‘caste load’), analogous to the result of intralocus sexual antagonism. The aim of this review is to create a research framework for understanding caste antagonism, drawing in part upon conceptual similarities with sexual antagonism. By reviewing both the social insect and sexual antagonism literature, we highlight the current empirical evidence for caste antagonism, discuss social systems of interest, how antagonism might be resolved, and challenges for future research. We also introduce the idea that sexual and caste antagonism could interact, creating a three‐way antagonism over gene expression. This includes unpacking the implications of haplodiploidy for the outcome of this complex interaction.     

Morph‐dependent fitness and directional change of morph frequencies over time in a Dutch population of Common buzzards Buteo buteo

How genetic polymorphisms are maintained in a population is a key question in evolutionary ecology. Previous work on a plumage colour polymorphism in the common buzzard Buteo buteo suggested heterozygote advantage as the mechanism maintaining the co‐existence of three morphs (light, intermediate and dark). We took advantage of 20 years of life‐history data collected in a Dutch population to replicate earlier studies on the relationship between colour morph and fitness in this species. We examined differences between morphs in adult apparent survival, breeding success, annual number of fledglings produced and cumulative reproductive success. We found that cumulative reproductive success differed among morphs, with the intermediate morph having highest fitness. We also found assortative mating for colour morph, whereby assortative pairs were more likely to produce offspring and had longer‐lasting pair bonds than disassortative pairs. Over the 20‐year study period, the proportion of individuals with an intermediate morph increased. This apparent evolutionary change did not just arise from selection on individual phenotypes, but also from fitness benefits of assortative mating. The increased frequency of intermediates might also be due to immigration or drift. We hypothesize that genetic variation is maintained through spatial variation in selection pressures. Further studies should investigate morph‐dependent dispersal behaviour and habitat choice.