Developmental transcriptomes predict adult social behaviours in the socially flexible sweat bee, Lasioglossum baleicum.

Natural variation can provide important insights into the genetic and environmental factors that shape social behaviour and its evolution. The sweat bee, Lasioglossum baleicum, is a socially flexible bee capable of producing both solitary and eusocial nests. We demonstrate that within a single nesting aggregation, soil temperatures are a strong predictor of the social structure of nests. Sites with warmer temperatures in the spring have a higher frequency of social nests than cooler sites, perhaps because warmer temperatures provide a longer reproductive window for those nests. To identify the molecular correlates of this behavioural variation, we generated a de novo genome assembly for L. baleicum, and we used transcriptomic profiling to compare adults and developing offspring from eusocial and solitary nests. We find that adult, reproductive females have similar expression profiles regardless of social structure in the nest, but that there are strong differences between reproductive females and workers from social nests. We also find substantial differences in the transcriptomic profiles of stage-matched pupae from warmer, social-biased sites compared to cooler, solitary-biased sites. These transcriptional differences are strongly predictive of adult reproductive state, suggesting that the developmental environment may set the stage for adult behaviours in L. baleicum. Together, our results help to characterize the molecular mechanisms shaping variation in social behaviour and highlight a potential role of environmental tuning during development as a factor shaping adult behaviour and physiology in this socially flexible bee.

Developmental transcriptomes predict adult social behaviors in the socially flexible sweat bee, Lasioglossum baleicum.

Natural variation can provide important insights into the genetic and environmental factors that shape social behavior and its evolution. The sweat bee, Lasioglossum baleicum , is a socially flexible bee capable of producing both solitary and eusocial nests. We demonstrate that within a single nesting aggregation, soil temperatures are a strong predictor of the social structure of nests. Sites with warmer temperatures in the spring have a higher frequency of social nests than cooler sites, perhaps because warmer temperatures provide a longer reproductive window for those nests. To identify the molecular correlates of this behavioral variation, we generated a de novo genome assembly for L. baleicum , and we used transcriptomic profiling to compare adults and developing offspring from eusocial and solitary nests. We find that adult, reproductive females have similar expression profiles regardless of social structure in the nest, but that there are strong differences between reproductive females and workers from social nests. We also find substantial differences in the transcriptomic profiles of stage-matched pupae from warmer, social-biased sites compared to cooler, solitary-biased sites. These transcriptional differences are strongly predictive of adult reproductive state, suggesting that the developmental environment may set the stage for adult behaviors in L. baleicum . Together, our results help to characterize the molecular mechanisms shaping variation in social behavior and highlight a potential role of environmental tuning during development as a factor shaping adult behavior and physiology in this socially flexible bee.

Single-cell dissection of aggression in honeybee colonies.

Understanding how genotypic variation results in phenotypic variation is especially difficult for collective behaviour because group phenotypes arise from complex interactions among group members. A genome-wide association study identified hundreds of genes associated with colony-level variation in honeybee aggression, many of which also showed strong signals of positive selection, but the influence of these 'colony aggression genes' on brain function was unknown. Here we use single-cell (sc) transcriptomics and gene regulatory network (GRN) analyses to test the hypothesis that genetic variation for colony aggression influences individual differences in brain gene expression and/or gene regulation. We compared soldiers, which respond to territorial intrusion with stinging attacks, and foragers, which do not. Colony environment showed stronger influences on soldier-forager differences in brain gene regulation compared with brain gene expression. GRN plasticity was strongly associated with colony aggression, with larger differences in GRN dynamics detected between soldiers and foragers from more aggressive relative to less aggressive colonies. The regulatory dynamics of subnetworks composed of genes associated with colony aggression genes were more strongly correlated with each other across different cell types and brain regions relative to other genes, especially in brain regions involved with olfaction and vision and multimodal sensory integration, which are known to mediate bee aggression. These results show how group genetics can shape a collective phenotype by modulating individual brain gene regulatory network architecture.

Convergent and complementary selection shaped gains and losses of eusociality in sweat bees.

Sweat bees have repeatedly gained and lost eusociality, a transition from individual to group reproduction. Here we generate chromosome-length genome assemblies for 17 species and identify genomic signatures of evolutionary trade-offs associated with transitions between social and solitary living. Both young genes and regulatory regions show enrichment for these molecular patterns. We also identify loci that show evidence of complementary signals of positive and relaxed selection linked specifically to the convergent gains and losses of eusociality in sweat bees. This includes two pleiotropic proteins that bind and transport juvenile hormone (JH)-a key regulator of insect development and reproduction. We find that one of these proteins is primarily expressed in subperineurial glial cells that form the insect blood-brain barrier and that brain levels of JH vary by sociality. Our findings are consistent with a role of JH in modulating social behaviour and suggest that eusocial evolution was facilitated by alteration of the proteins that bind and transport JH, revealing how an ancestral developmental hormone may have been co-opted during one of life's major transitions. More broadly, our results highlight how evolutionary trade-offs have structured the molecular basis of eusociality in these bees and demonstrate how both directional selection and release from constraint can shape trait evolution.

NAPS: Integrating pose estimation and tag-based tracking.

1. Significant advances in computational ethology have allowed the quantification of behaviour in unprecedented detail. Tracking animals in social groups, however, remains challenging as most existing methods can either capture pose or robustly retain individual identity over time but not both. 2. To capture finely resolved behaviours while maintaining individual identity, we built NAPS (NAPS is ArUco Plus SLEAP), a hybrid tracking framework that combines state-of-the-art, deep learning-based methods for pose estimation (SLEAP) with unique markers for identity persistence (ArUco). We show that this framework allows the exploration of the social dynamics of the common eastern bumblebee (Bombus impatiens). 3. We provide a stand-alone Python package for implementing this framework along with detailed documentation to allow for easy utilization and expansion. We show that NAPS can scale to long timescale experiments at a high frame rate and that it enables the investigation of detailed behavioural variation within individuals in a group. 4. Expanding the toolkit for capturing the constituent behaviours of social groups is essential for understanding the structure and dynamics of social networks. NAPS provides a key tool for capturing these behaviours and can provide critical data for understanding how individual variation influences collective dynamics.

SimiC enables the inference of complex gene regulatory dynamics across cell phenotypes.

Single-cell RNA-Sequencing has the potential to provide deep biological insights by revealing complex regulatory interactions across diverse cell phenotypes at single-cell resolution. However, current single-cell gene regulatory network inference methods produce a single regulatory network per input dataset, limiting their capability to uncover complex regulatory relationships across related cell phenotypes. We present SimiC, a single-cell gene regulatory inference framework that overcomes this limitation by jointly inferring distinct, but related, gene regulatory dynamics per phenotype. We show that SimiC uncovers key regulatory dynamics missed by previously proposed methods across a range of systems, both model and non-model alike. In particular, SimiC was able to uncover CAR T cell dynamics after tumor recognition and key regulatory patterns on a regenerating liver, and was able to implicate glial cells in the generation of distinct behavioral states in honeybees. SimiC hence establishes a new approach to quantitating regulatory architectures between distinct cellular phenotypes, with far-reaching implications for systems biology.

Context-dependent influence of threat on honey bee social network dynamics and brain gene expression.

Adverse social experience affects social structure by modifying the behavior of individuals, but the relationship between an individual's behavioral state and its response to adversity is poorly understood. We leveraged naturally occurring division of labor in honey bees and studied the biological embedding of environmental threat using laboratory assays and automated behavioral tracking of whole colonies. Guard bees showed low intrinsic levels of sociability compared with foragers and nurse bees, but large increases in sociability following exposure to a threat. Threat experience also modified the expression of caregiving-related genes in a brain region called the mushroom bodies. These results demonstrate that the biological embedding of environmental experience depends on an individual's societal role and, in turn, affects its future sociability.

Parallel mechanisms of visual memory formation across distinct regions of the honey bee brain.

Visual learning is vital to the behavioral ecology of the Western honey bee (Apis mellifera). Honey bee workers forage for floral resources, a behavior that requires the learning and long-term memory of visual landmarks, but how these memories are mapped to the brain remains poorly understood. To address this gap in our understanding, we collected bees that successfully learned visual associations in a conditioned aversion paradigm and compared gene expression correlates of memory formation in the mushroom bodies, a higher-order sensory integration center classically thought to contribute to learning, as well as the optic lobes, the primary visual neuropil responsible for sensory transduction of visual information. We quantified expression of CREB and CaMKII, two classical genetic markers of learning, and fen-1, a gene specifically associated with punishment learning in vertebrates. As expected, we found substantial involvement of the mushroom bodies for all three markers but additionally report the involvement of the optic lobes across a similar time course. Our findings imply the molecular involvement of a sensory neuropil during visual associative learning parallel to a higher-order brain region, furthering our understanding of how a tiny brain processes environmental signals.

Neural and Molecular Mechanisms of Biological Embedding of Social Interactions.

Animals operate in complex environments, and salient social information is encoded in the nervous system and then processed to initiate adaptive behavior. This encoding involves biological embedding, the process by which social experience affects the brain to influence future behavior. Biological embedding is an important conceptual framework for understanding social decision-making in the brain, as it encompasses multiple levels of organization that regulate how information is encoded and used to modify behavior. The framework we emphasize here is that social stimuli provoke short-term changes in neural activity that lead to changes in gene expression on longer timescales. This process, simplified-neurons are for today and genes are for tomorrow-enables the assessment of the valence of a social interaction, an appropriate and rapid response, and subsequent modification of neural circuitry to change future behavioral inclinations in anticipation of environmental changes. We review recent research on the neural and molecular basis of biological embedding in the context of social interactions, with a special focus on the honeybee.

Neurodevelopmental and transcriptomic effects of CRISPR/Cas9-induced somatic orco mutation in honey bees.

In insects, odorant receptors facilitate olfactory communication and require the functionality of the highly conserved co-receptor gene orco. Genome editing studies in a few species of ants and moths have revealed that orco can also have a neurodevelopmental function, in addition to its canonical role in adult olfaction, discovered first in Drosophila melanogaster. To extend this analysis, we determined whether orco mutations also affect the development of the adult brain of the honey bee Apis mellifera, an important model system for social behavior and chemical communication. We used CRISPR/Cas9 to knock out orco and examined anatomical and molecular consequences. To increase efficiency, we coupled embryo microinjection with a laboratory egg collection and in vitro rearing system. This new workflow advances genomic engineering technologies in honey bees by overcoming restrictions associated with field studies. We used Sanger sequencing to quickly select individuals with complete orco knockout for neuroanatomical analyses and later validated and described the mutations with amplicon sequencing. Mutant bees had significantly fewer glomeruli, smaller total volume of all the glomeruli, and higher mean individual glomerulus volume in the antennal lobe compared to wild-type controls. RNA-Sequencing revealed that orco knockout also caused differential expression of hundreds of genes in the antenna, including genes related to neural development and genes encoding odorant receptors. The expression of other types of chemoreceptor genes was generally unaffected, reflecting specificity of CRISPR activity in this study. These results suggest that neurodevelopmental effects of orco are related to specific insect life histories.

Behavior-related gene regulatory networks: A new level of organization in the brain.

Neuronal networks are the standard heuristic model today for describing brain activity associated with animal behavior. Recent studies have revealed an extensive role for a completely distinct layer of networked activities in the brain-the gene regulatory network (GRN)-that orchestrates expression levels of hundreds to thousands of genes in a behavior-related manner. We examine emerging insights into the relationships between these two types of networks and discuss their interplay in spatial as well as temporal dimensions, across multiple scales of organization. We discuss properties expected of behavior-related GRNs by drawing inspiration from the rich literature on GRNs related to animal development, comparing and contrasting these two broad classes of GRNs as they relate to their respective phenotypic manifestations. Developmental GRNs also represent a third layer of network biology, playing out over a third timescale, which is believed to play a crucial mediatory role between neuronal networks and behavioral GRNs. We end with a special emphasis on social behavior, discuss whether unique GRN organization and cis-regulatory architecture underlies this special class of behavior, and review literature that suggests an affirmative answer.

Meta-analysis of honey bee neurogenomic response links Deformed wing virus type A to precocious behavioral maturation.

Crop pollination by the western honey bee Apis mellifera is vital to agriculture but threatened by alarmingly high levels of colony mortality, especially in Europe and North America. Colony loss is due, in part, to the high viral loads of Deformed wing virus (DWV), transmitted by the ectoparasitic mite Varroa destructor, especially throughout the overwintering period of a honey bee colony. Covert DWV infection is commonplace and has been causally linked to precocious foraging, which itself has been linked to colony loss. Taking advantage of four brain transcriptome studies that unexpectedly revealed evidence of covert DWV-A infection, we set out to explore whether this effect is due to DWV-A mimicking naturally occurring changes in brain gene expression that are associated with behavioral maturation. Consistent with this hypothesis, we found that brain gene expression profiles of DWV-A infected bees resembled those of foragers, even in individuals that were much younger than typical foragers. In addition, brain transcriptional regulatory network analysis revealed a positive association between DWV-A infection and transcription factors previously associated with honey bee foraging behavior. Surprisingly, single-cell RNA-Sequencing implicated glia, not neurons, in this effect; there are relatively few glial cells in the insect brain and they are rarely associated with behavioral plasticity. Covert DWV-A infection also has been linked to impaired learning, which together with precocious foraging can lead to increased occurrence of infected bees from one colony mistakenly entering another colony, especially under crowded modern apiary conditions. These findings provide new insights into the mechanisms by which DWV-A affects honey bee health and colony survival.

Valence of social information is encoded in different subpopulations of mushroom body Kenyon cells in the honeybee brain.

Over 600 Myr of evolutionary divergence between vertebrates and invertebrates is associated with considerable neuroanatomical variation both across and within these lineages. By contrast, valence encoding is an important behavioural trait that is evolutionarily conserved across vertebrates and invertebrates, and enables individuals to distinguish between positive (potentially beneficial) and negative (potentially harmful) situations. We tested the hypothesis that social interactions of positive and negative valence are modularly encoded in the honeybee brain (i.e. encoded in different cellular subpopulations) as in vertebrate brains. In vertebrates, neural activation patterns are distributed across distinct parts of the brain, suggesting that discrete circuits encode positive or negative stimuli. Evidence for this hypothesis would suggest a deep homology of neural organization between insects and vertebrates for valence encoding, despite vastly different brain sizes. Alternatively, overlapping localization of valenced social information in the brain would imply a 're-use' of circuitry in response to positive and negative social contexts, potentially to overcome the energetic constraints of a tiny brain. We used immediate early gene expression to map positively and negatively valenced social interactions in the brain of the western honeybee Apis mellifera. We found that the valence of a social signal is represented by distinct anatomical subregions of the mushroom bodies, an invertebrate sensory neuropil associated with social behaviour, multimodal sensory integration, learning and memory. Our results suggest that the modularization of valenced social information in the brain is a fundamental property of neuroanatomical organization.

Division of labor in honey bees is associated with transcriptional regulatory plasticity in the brain.

Studies in evolutionary and developmental biology show that relationships between transcription factors (TFs) and their target genes can be altered to result in novel regulatory relationships that generate phenotypic plasticity. We hypothesized that context-dependent shifts in the nervous system associated with behavior may also be linked to changes in TF-target relationships over physiological time scales. We tested this hypothesis using honey bee (Apis mellifera) division of labor as a model system by performing bioinformatic analyses of previously published brain transcriptomic profiles together with new RNAi and behavioral experiments. The bioinformatic analyses identified five TFs that exhibited strong signatures of regulatory plasticity as a function of division of labor. RNAi targeting of one of these TFs (broad complex) and a related TF that did not exhibit plasticity (fushi tarazu transcription factor 1) was administered in conjunction with automated analyses of foraging behavior in the field, laboratory assays of aggression and brood care behavior, and endocrine treatments. The results showed that changes in the regulatory relationships of these TFs were associated with behavioral state, social context and endocrine state. These findings provide the first empirical evidence that TF-target relationships in the brain are altered in conjunction with behavior and social context. They also suggest that one mechanism for this plasticity involves pleiotropic TFs high up in regulatory hierarchies producing behavior-specific transcriptional responses by activating different downstream TFs to induce discrete context-dependent transcriptional cascades. These findings provide new insights into the dynamic nature of the transcriptional regulatory architecture underlying behavior in the brain.

Adult stem cells in the knifefish cerebellum.

Adult neurogenesis has been described in dozens of brain regions in teleost fish, with the largest number of new neurons being generated in the cerebellum. Here, we characterized the cerebellar neural stem/progenitor cells (NSPCs) in the brown ghost knifefish (Apteronotus leptorhynchus), an established model system of adult neurogenesis. The majority of the new cerebellar cells arise from neurogenic niches located medially, at the interface of the dorsal/ventral molecular layers and the granular layer. NSPCs within these niches give rise to transit-amplifying progenitors which populate the molecular layer, where they continue to proliferate during their migration toward target areas in the granular layer. At any given time, the majority of proliferating cells are located in the molecular layer. Immunohistochemical staining revealed that the stem cell markers Sox2, Meis1/2/3, Islet1, and, to a lesser extent, Pax6, are widely expressed in all regions of the adult cerebellum. A large subpopulation of these NSPCs coexpress S100, GFAP, and/or vimentin, indicating astrocytic identity. This is further supported by the specific effect of the gliotoxin l-methionine sulfoximine, which leads to a targeted decrease in the number of GFAP+ cells that coexpress Sox2 or the proliferation marker PCNA. Pulse-chase analysis of the label size associated with new cells after administration of 5-bromo-2'-deoxyuridine demonstrated that, on average, two additional cell divisions occur after completion of the initial mitotic cycle. Overall numbers of NSPCs in the cerebellum niches increase consistently over time, presumably in parallel with the continuous growth of the brain.

Age-related changes in stem cell dynamics, neurogenesis, apoptosis, and gliosis in the adult brain: a novel teleost fish model of negligible senescence.

Adult neurogenesis, the generation of new neurons in the adult central nervous system, is a reported feature of all examined vertebrate species. However, a dramatic decline in the rates of cell proliferation and neuronal differentiation occurs in mammals, typically starting near the onset of sexual maturation. In the present study, we examined possible age-related changes associated with adult neurogenesis in the brain of brown ghost knifefish (Apteronotus leptorhynchus), a teleost fish distinguished by its enormous neurogenic potential. Contrary to the well-established alterations in the mammalian brain during aging, in the brain of this teleostean species we could not find evidence for any significant age-related decline in the absolute levels of stem/progenitor cell proliferation, neuronal and glial differentiation, or long-term survival of newly generated cells. Moreover, there was no indication that the amount of glial fibrillary acidic protein or the number of apoptotic cells in the brain was altered significantly over the course of adult life. We hypothesize that this first demonstration of negligible cellular senescence in the vertebrate brain is related to the continued growth of this species and to the lack of reproductive senescence during adulthood. The establishment of the adult brain of this species as a novel model of negligible senescence provides new opportunities for the advancement of our understanding of the biology of aging and the fundamental mechanisms that underlie senescence in the brain.