Single-cell dissection of a collective behaviour in honeybees

Understanding how genotypic variation results in phenotypic variation, a major challenge in biology, is especially difficult for collective behaviour because collective group phenotypes arise from complex interactions between group members1. Honeybees aggressively defend their colony from attacks with highly integrated collective behaviour in which different groups of bees play specific roles, giving rise to distinct colony-level differences in aggression. A previous genome-wide association study of a population of Africanized honeybees (Apis mellifera scutellata) from Puerto Rico that recently evolved decreased aggression identified hundreds of genes with single nucleotide polymorphisms (SNPs) that associated with colony-level variation in aggression2. Many of these SNPs also showed strong signals of selection for decreased aggression2,3, but their influence on brain function was unknown. Using brain single-cell (sc) transcriptomics and sc gene regulatory network analysis, we show here that variants of these genes give rise to genetic differences in transcription factor-target gene relationships. These differences involved the activity of several TFs, some that have been previously associated with aggression, like single stranded-binding protein c31A, and some that have been associated with tissue morphogenesis but not behaviour, like apontic. The activity of these and other TFs was located in specific brain cell populations related to olfaction and vision, the two sensory modalities that bees use in colony defence. They also implicate metabolism of serotonin, a neurochemical already known to influence honeybee aggression, but not from a genetic perspective. Surprisingly, genetic differences were more pronounced in the brains of forager bees than in similarly aged but more aggressive soldier bees, pointing to an evolutionary change in division of labour for colony defence. Our results demonstrate how group genetics can shape a collective phenotype by modulating individual brain gene regulatory network architecture.