For social animals, the genotypes of group members affect the social environment, and thus individual behavior, often indirectly. We used genome-wide association studies (GWAS) to determine the influence of individual vs. group genotypes on aggression in honey bees. Aggression in honey bees arises from the coordinated actions of colony members, primarily nonreproductive “soldier” bees, and thus, experiences evolutionary selection at the colony level. Here, we show that individual behavior is influenced by colony environment, which in turn, is shaped by allele frequency within colonies. Using a population with a range of aggression, we sequenced individual whole genomes and looked for genotype–behavior associations within colonies in a common environment. There were no significant correlations between individual aggression and specific alleles. By contrast, we found strong correlations between colony aggression and the frequencies of specific alleles within colonies, despite a small number of colonies. Associations at the colony level were highly significant and were very similar among both soldiers and foragers, but they covaried with one another. One strongly significant association peak, containing an ortholog of the Drosophila sensory gene dpr4 on linkage group (chromosome) 7, showed strong signals of both selection and admixture during the evolution of gentleness in a honey bee population. We thus found links between colony genetics and group behavior and also, molecular evidence for group-level selection, acting at the colony level. We conclude that group genetics dominates individual genetics in determining the fatal decision of honey bees to sting. Significance Honey bee colony defense is an emergent trait composed of individual aggressive responses. Here, we investigated the relationship between individual genotype, colony allele frequency, and aggression in individual bees. Our findings show that the colony-level defense response strongly correlates with colony-level allele frequency in a way that can be used to identify causative genomic regions. Importantly, we were able to validate a key associated region as also being under selection. As very similar allele frequency correlations are observed in both soldier and forager bees, we conclude that group genetics is more important than individual genetics in this case, giving further insight into the relationships between “nature,” “nurture,” and behavioral evolution.

Developmental plasticity generates phenotypic variation, but how it contributes to evolutionary change is unclear. Phenotypes of individuals in caste-based (eusocial) societies are particularly sensitive to developmental processes, and the evolutionary origins of eusociality may be rooted in developmental plasticity of ancestral forms. We used an integrative genomics approach to evaluate the relationships among developmental plasticity, molecular evolution, and social behavior in a bee species (Megalopta genalis) that expresses flexible sociality, and thus provides a window into the factors that may have been important at the evolutionary origins of eusociality. We find that differences in social behavior are derived from genes that also regulate sex differentiation and metamorphosis. Positive selection on social traits is influenced by the function of these genes in development. We further identify evidence that social polyphenisms may become encoded in the genome via genetic changes in regulatory regions, specifically in transcription factor binding sites. Taken together, our results provide evidence that developmental plasticity provides the substrate for evolutionary novelty and shapes the selective landscape for molecular evolution in a major evolutionary innovation: Eusociality.

It is often assumed that the mutation rate is an evolutionarily optimized property of a taxon. The relevant mutation rate is for mutations that affect fitness, U, but the strength of selection on the mutation rate depends on the average effect of a mutation. Determination of U is complicated by the possibility that mutational effects depend on the particular environmental context in which the organism exists. It has been suggested that the effects of deleterious mutations are typically magnified in stressful environments, but most studies confound genotype with environment, so it is unclear to what extent environmental specificity of mutations is specific to a particular starting genotype. We report a study designed to separate effects of species, genotype, and environment on the degradation of fitness resulting from new mutations. Mutations accumulated for >200 generations at 20° in two strains of two species of nematodes that differ in thermal sensitivity. Caenorhabditis briggsae and C. elegans have similar demography at 20°, but C. elegans suffers markedly reduced fitness at 25°. We find little evidence that mutational properties differ depending on environmental conditions and mutational correlations between environments are close to those expected if effects were identical in both environments.

Circadian rhythms in social insects are highly plastic and are modulated by multiple factors. In addition, complex behaviors such as sun-compass orientation and time learning are clearly regulated by the circadian system in these organisms. Despite these unique features of social insect clocks, the mechanisms as well as the functional and evolutionary relevance of these traits remain largely unknown. Here we show a modification of the Drosophila activity monitoring (DAM) system that allowed us to measure locomotor rhythms of the honey bee, Apis mellifera (three variants; gAHB, carnica and caucasica), and two paper wasps (Polistes crinitus and Mischocyttarus phthisicus). A side-by-side comparison of the endogenous period under constant darkness (free-running period) led us to the realization that these social insects exhibit significant deviations from the Earth9s 24 h rotational period as well as a large degree of inter-individual variation compared with Drosophila. Experiments at different temperatures, using honey bees as a model, revealed that testing the endogenous rhythm at 35°C, which is the hive9s core temperature, results in average periods closer to 24 h compared with 25°C (23.8 h at 35°C versus 22.7 h at 25°C). This finding suggests that the degree of tuning of circadian temperature compensation varies among different organisms. We expect that the commercial availability, cost-effectiveness and integrated nature of this monitoring system will facilitate the growth of the circadian field in these social insects and catalyze our understanding of the mechanisms as well as the functional and evolutionary relevance of circadian rhythms.

Highly aggressive Africanized honeybees (AHB) invaded Puerto Rico (PR) in 1994, displacing gentle European honeybees (EHB) in many locations. Gentle AHB (gAHB), unknown anywhere else in the world, subsequently evolved on the island within a few generations. Here we sequence whole genomes from gAHB and EHB populations, as well as a North American AHB population, a likely source of the founder AHB on PR. We show that gAHB retains high levels of genetic diversity after evolution of gentle behaviour, despite selection on standing variation. We observe multiple genomic loci with significant signatures of selection. Rapid evolution during colonization of novel habitats can generate major changes to characteristics such as morphological or colouration traits, usually controlled by one or more major genetic loci. Here we describe a soft selective sweep, acting at multiple loci across the genome, that occurred during, and may have mediated, the rapid evolution of a behavioural trait.

SUMMARY: Honey bees provide a model system to elucidate the relationship between sociality and complex behaviors within the same species, as females (workers) are highly social and males (drones) are more solitary. We report on aversive learning studies in drone and worker honey bees (Apis mellifera anatolica) in escape, punishment and discriminative punishment situations. In all three experiments, a newly developed electric shock avoidance assay was used. The comparisons of expected and observed responses were performed with conventional statistical methods and a systematic randomization modeling approach called object oriented modeling. The escape experiment consisted of two measurements recorded in a master–yoked paradigm: frequency of response and latency to respond following administration of shock. Master individuals could terminate an unavoidable shock triggered by a decrementing 30 s timer by crossing the shuttlebox centerline following shock activation. Across all groups, there was large individual response variation. When assessing group response frequency and latency, master subjects performed better than yoked subjects for both workers and drones. In the punishment experiment, individuals were shocked upon entering the shock portion of a bilaterally wired shuttlebox. The shock portion was spatially static and unsignalled. Only workers effectively avoided the shock. The discriminative punishment experiment repeated the punishment experiment but included a counterbalanced blue and yellow background signal and the side of shock was manipulated. Drones correctly responded less than workers when shock was paired with blue. However, when shock was paired with yellow there was no observable difference between drones and workers.