Genome-wide analysis reveals differences in brain gene expression patterns associated with caste and reproductive status in honey bees (Apis mellifera)

A key characteristic of eusocial species is reproductive division of labour. Honey bee colonies typically have a single reproductive queen and thousands of sterile workers. Adult queens differ dramatically from workers in anatomy, physiology, behaviour and lifespan. Young female workers can activate their ovaries and initiate egg laying; these 'reproductive' workers differ from sterile workers in anatomy, physiology, and behaviour. These differences, however, are on a much smaller scale than those observed between the queen and worker castes. Here, we use microarrays to monitor expression patterns of several thousand genes in the brains of same-aged virgin queens, sterile workers, and reproductive workers. We found large differences in expression between queens and both worker groups (~2000 genes), and much smaller differences between sterile and reproductive workers (221 genes). The expression patterns of these 221 genes in reproductive workers are more queen-like, and may represent a core group of genes associated with reproductive physiology. Furthermore, queens and reproductive workers preferentially up-regulate genes associated with the nurse bee behavioural state, which supports the hypothesis of an evolutionary link between worker division of labour and molecular pathways related to reproduction. Finally, several functional groups of genes associated with longevity in other species are significantly up-regulated in queens. Identifying the genes that underlie the differences between queens, sterile workers, and reproductive workers will allow us to begin to characterize the molecular mechanisms underlying the evolution of social behaviour and large-scale remodelling of gene networks associated with polyphenisms.     

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.     

Development and evolution of caste dimorphism in honeybees – a modeling approach

The difference in phenotypes of queens and workers is a hallmark of the highly eusocial insects. The caste dimorphism is often described as a switch‐controlled polyphenism, in which environmental conditions decide an individual's caste. Using theoretical modeling and empirical data from honeybees, we show that there is no discrete larval developmental switch. Instead, a combination of larval developmental plasticity and nurse worker feeding behavior make up a colony‐level social and physiological system that regulates development and produces the caste dimorphism. Discrete queen and worker phenotypes are the result of discrete feeding regimes imposed by nurses, whereas a range of experimental feeding regimes produces a continuous range of phenotypes. Worker ovariole numbers are reduced through feeding‐regime‐mediated reduction in juvenile hormone titers, involving reduced sugar in the larval food. Based on the mechanisms identified in our analysis, we propose a scenario of the evolutionary history of honeybee development and feeding regimes.     

Abdominal microbial communities in ants depend on colony membership rather than caste and are linked to colony productivity

Gut bacteria aid their host in digestion and pathogen defense, and bacterial communities that differ in diversity or composition may vary in their ability to do so. Typically, the gut microbiomes of animals living in social groups converge as members share a nest environment and frequently interact. Social insect colonies, however, consist of individuals that differ in age, physiology, and behavior, traits that could affect gut communities or that expose the host to different bacteria, potentially leading to variation in the gut microbiome within colonies. Here we asked whether bacterial communities in the abdomen of Temnothorax nylanderi ants, composed largely of the gut microbiome, differ between different reproductive and behavioral castes. We compared microbiomes of queens, newly eclosed workers, brood carers, and foragers by high‐throughput 16S rRNA sequencing. Additionally, we sampled individuals from the same colonies twice, in the field and after 2 months of laboratory housing. To disentangle the effects of laboratory environment and season on microbial communities, additional colonies were collected at the same location after 2 months. There were no large differences between ant castes, although queens harbored more diverse microbial communities than workers. Instead, we found effects of colony, environment, and season on the abdominal microbiome. Interestingly, colonies with more diverse communities had produced more brood. Moreover, the queens' microbiome composition was linked to egg production. Although long‐term coevolution between social insects and gut bacteria has been repeatedly evidenced, our study is the first to find associations between abdominal microbiome characteristics and colony productivity in social insects.     

Environmental determinants of population divergence in life‐history traits for an invasive species: climate,seasonality and natural enemies

Invasive species cope with novel environments through both phenotypic plasticity and evolutionary change. However, the environmental factors that cause evolutionary divergence in invasive species are poorly understood. We developed predictions for how different life‐history traits, and plasticity in those traits, may respond to environmental gradients in seasonal temperatures, season length and natural enemies. We then tested these predictions in four geographic populations of the invasive cabbage white butterfly (Pieris rapae) from North America. We examined the influence of two rearing temperatures (20 and 26.7 °C) on pupal mass, pupal development time, immune function and fecundity. As predicted, development time was shorter and immune function was greater in populations adapted to longer season length. Also, phenotypic plasticity in development time was greater in regions with shorter growing seasons. Populations differed significantly in mean and plasticity of body mass and fecundity, but these differences were not associated with seasonal temperatures or season length. Our study shows that some life‐history traits, such as development time and immune function, can evolve rapidly in response to latitudinal variation in season length and natural enemies, whereas others traits did not. Our results also indicate that phenotypic plasticity in development time can also diverge rapidly in response to environmental conditions for some traits.     

Sexually dimorphic gene expression and transcriptome evolution provide mixed evidence for a fast‐Z effect in Heliconius

Sex chromosomes have different evolutionary properties compared to autosomes due to their hemizygous nature. In particular, recessive mutations are more readily exposed to selection, which can lead to faster rates of molecular evolution. Here, we report patterns of gene expression and molecular evolution for a group of butterflies. First, we improve the completeness of the Heliconius melpomene reference annotation, a neotropical butterfly with a ZW sex determination system. Then, we analyse RNA from male and female whole abdomens and sequence female ovary and gut tissue to identify sex‐ and tissue‐specific gene expression profiles in H. melpomene. Using these expression profiles, we compare (a) sequence divergence and polymorphism; (b) the strength of positive and negative selection; and (c) rates of adaptive evolution, for Z and autosomal genes between two species of Heliconius butterflies, H. melpomene and H. erato. We show that the rate of adaptive substitutions is higher for Z than autosomal genes, but contrary to expectation, it is also higher for male‐biased than female‐biased genes. Additionally, we find no significant increase in the rate of adaptive evolution or purifying selection on genes expressed in ovary tissue, a heterogametic‐specific tissue. Our results contribute to a growing body of literature from other ZW systems that also provide mixed evidence for a fast‐Z effect where hemizygosity influences the rate of adaptive substitutions.     

Cross‐generational effects of temperature on flight performance,and associated life‐history traits in an insect

Nongenetic parental effects may affect offspring phenotype, and in species with multiple generations per year, these effects may cause life‐history traits to vary over the season. We investigated the effects of parental, offspring developmental and offspring adult temperatures on a suite of life‐history traits in the globally invasive agricultural pest Grapholita molesta. A low parental temperature resulted in female offspring that developed faster at low developmental temperature compared with females whose parents were reared at high temperature. Furthermore, females whose parents were reared at low temperature were heavier and more fecund and had better flight abilities than females whose parents were reared at high temperature. In addition to these cross‐generational effects, females developed at low temperature had similar flight abilities at low and high ambient temperatures, whereas females developed at high temperature had poorer flight abilities at low than at high ambient temperature. Our findings demonstrate a pronounced benefit of low parental temperature on offspring performance, as well as between‐ and within‐generation effects of acclimation to low temperature. In cooler environments, the offspring generation is expected to develop more rapidly than the parental generation and to comprise more fecund and more dispersive females. By producing phenotypes that are adaptive to the conditions inducing them as well as heritable, cross‐generational plasticity can influence the evolutionary trajectory of populations. The potential for short‐term acclimation to low temperature may allow expanding insect populations to better cope with novel environments and may help to explain the spread and establishment of invasive species.     

Temperature‐mediated trade‐offs and changes in life‐history integration in two slipper limpets (Gastropoda: Calyptraeidae) with planktotrophic development

Intraspecific variation in egg size and hatching size, and the genetic and environmental trade‐offs that contribute to variation, are the basis of the evolution of life histories. The present study examined both univariate and multivariate temperature‐mediated plasticity of life‐history traits, as well as temperature‐mediated trade‐offs in egg size and clutch size, in two planktotrophic species of marine slipper limpets, Crepidula. Previous work with two species of Crepidula with large eggs and lecithotrophic development has shown a significant effect of temperature on egg size and hatching size. To further examine the effect of temperature on egg size in Crepidula, the effects of temperature on egg size and hatching size, as well as the possible trade‐offs with other the life‐history features, were examined for two planktotrophic species: Crepidula incurva and Crepidula cf. marginalis. Field‐collected juveniles were raised at 23 or 28 °C and egg size, hatching size, capsules/brood, eggs/capsule, time to hatch, interbrood interval, and final body weight were recorded. Consistent with results for the lecithotrophic Crepidula, egg size and hatching size decreased with temperature in the planktotrophic species. The affects of maternal identity and individual brood account for more than half of the intraspecific variation in egg size and hatching size. Temperature also showed a significant effect on reproductive rate, with time to hatch and interbrood interval both decreasing with increasing temperature. However, temperature had contrasting effects on the number of offspring. Crepidula cf. marginalis has significantly more eggs/capsule and therefore more eggs per brood at 28 °C compared to 23 °C, although capsules/brood did not vary with temperature. Crepidula incurva, on the other hand, produced significantly more capsules/brood and more eggs per brood at the lower temperature, whereas the number of eggs/capsule did not vary with temperature. The phenotypic variance–covariance matrix of life‐history variables showed a greater response to temperature in C. incurva than in C. cf. marginalis, and temperature induced trade‐offs between offspring size and number differ between the species. These differences suggest that temperature changes as a result of seasonal upwelling along the coast of Panama will effect the reproduction and evolution of life histories of these two co‐occurring species differently. © 2012 The Linnean Society of London, Biological Journal of the Linnean Society, 2012, ?? , ??–??.     

Evolution and plasticity: Divergence of song discrimination is faster in birds with innate song than in song learners in Neotropical passerine birds

Plasticity is often thought to accelerate trait evolution and speciation. For example, plasticity in birdsong may partially explain why clades of song learners are more diverse than related clades with innate song. This “song learning” hypothesis predicts that (1) differences in song traits evolve faster in song learners, and (2) behavioral discrimination against allopatric song (a proxy for premating reproductive isolation) evolves faster in song learners. We tested these predictions by analyzing acoustic traits and conducting playback experiments in allopatric Central American sister pairs of song learning oscines (N = 42) and nonlearning suboscines (N = 27). We found that nonlearners evolved mean acoustic differences slightly faster than did leaners, and that the mean evolutionary rate of song discrimination was 4.3 times faster in nonlearners than in learners. These unexpected results may be a consequence of significantly greater variability in song traits in song learners (by 54–79%) that requires song‐learning oscines to evolve greater absolute differences in song before achieving the same level of behavioral song discrimination as nonlearning suboscines. This points to “a downside of learning” for the evolution of species discrimination, and represents an important example of plasticity reducing the rate of evolution and diversification by increasing variability.     

The population genetic test Tajima's D identifies genes encoding pathogen‐associated molecular patterns and other virulence‐related genes in Ralstonia solanacearum

The detection of pathogen‐associated molecular patterns (PAMPs) by plant pattern recognition receptors (PRRs) is an essential part of plant immunity. Until recently, elf18, an epitope of elongation factor‐Tu (EF‐Tu), was the sole confirmed PAMP of Ralstonia solanacearum, the causal agent of bacterial wilt disease, limiting our understanding of R. solanacearum–plant interactions. Therefore, we set out to identify additional R. solanacearum PAMPs based on the hypothesis that genes encoding PAMPs are under selection to avoid recognition by plant PRRs. We calculated Tajima's D, a population genetic test statistic which identifies genes that do not evolve neutrally, for 3003 genes conserved in 37 R. solanacearum genomes. The screen flagged 49 non‐neutrally evolving genes, including not only EF‐Tu but also the gene for Cold Shock Protein C, which encodes the PAMP csp22. Importantly, an R. solanacearum allele of this PAMP was recently identified in a parallel independent study. Genes coding for efflux pumps, some with known roles in virulence, were also flagged by Tajima's D. We conclude that Tajima's D is a straightforward test to identify genes encoding PAMPs and other virulence‐related genes in plant pathogen genomes.