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.

Automated monitoring of honey bees with barcodes and artificial intelligence reveals two distinct social networks from a single affiliative behavior.

Barcode-based tracking of individuals is revolutionizing animal behavior studies, but further progress hinges on whether in addition to determining an individual's location, specific behaviors can be identified and monitored. We achieve this goal using information from the barcodes to identify tightly bounded image regions that potentially show the behavior of interest. These image regions are then analyzed with convolutional neural networks to verify that the behavior occurred. When applied to a challenging test case, detecting social liquid transfer (trophallaxis) in the honey bee hive, this approach yielded a 67% higher sensitivity and an 11% lower error rate than the best detector for honey bee trophallaxis so far. We were furthermore able to automatically detect whether a bee donates or receives liquid, which previously required manual observations. By applying our trophallaxis detector to recordings from three honey bee colonies and performing simulations, we discovered that liquid exchanges among bees generate two distinct social networks with different transmission capabilities. Finally, we demonstrate that our approach generalizes to detecting other specific behaviors. We envision that its broad application will enable automatic, high-resolution behavioral studies that address a broad range of previously intractable questions in evolutionary biology, ethology, neuroscience, and molecular biology.

Individual differences in honey bee behavior enabled by plasticity in brain gene regulatory networks.

Understanding the regulatory architecture of phenotypic variation is a fundamental goal in biology, but connections between gene regulatory network (GRN) activity and individual differences in behavior are poorly understood. We characterized the molecular basis of behavioral plasticity in queenless honey bee (Apis mellifera) colonies, where individuals engage in both reproductive and non-reproductive behaviors. Using high-throughput behavioral tracking, we discovered these colonies contain a continuum of phenotypes, with some individuals specialized for either egg-laying or foraging and 'generalists' that perform both. Brain gene expression and chromatin accessibility profiles were correlated with behavioral variation, with generalists intermediate in behavior and molecular profiles. Models of brain GRNs constructed for individuals revealed that transcription factor (TF) activity was highly predictive of behavior, and behavior-associated regulatory regions had more TF motifs. These results provide new insights into the important role played by brain GRN plasticity in the regulation of behavior, with implications for social evolution.

Brain microRNAs among social and solitary bees.

Evolutionary transitions to a social lifestyle in insects are associated with lineage-specific changes in gene expression, but the key nodes that drive these regulatory changes are unknown. We examined the relationship between social organization and lineage-specific microRNAs (miRNAs). Genome scans across 12 bee species showed that miRNA copy-number is mostly conserved and not associated with sociality. However, deep sequencing of small RNAs in six bee species revealed a substantial proportion (20-35%) of detected miRNAs had lineage-specific expression in the brain, 24-72% of which did not have homologues in other species. Lineage-specific miRNAs disproportionately target lineage-specific genes, and have lower expression levels than shared miRNAs. The predicted targets of lineage-specific miRNAs are not enriched for genes with caste-biased expression or genes under positive selection in social species. Together, these results suggest that novel miRNAs may coevolve with novel genes, and thus contribute to lineage-specific patterns of evolution in bees, but do not appear to have significant influence on social evolution. Our analyses also support the hypothesis that many new miRNAs are purged by selection due to deleterious effects on mRNA targets, and suggest genome structure is not as influential in regulating bee miRNA evolution as has been shown for mammalian miRNAs.

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.

Developmental plasticity shapes social traits and selection in a facultatively eusocial bee.

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.

Honey bee virus causes context-dependent changes in host social behavior.

Anthropogenic changes create evolutionarily novel environments that present opportunities for emerging diseases, potentially changing the balance between host and pathogen. Honey bees provide essential pollination services, but intensification and globalization of honey bee management has coincided with increased pathogen pressure, primarily due to a parasitic mite/virus complex. Here, we investigated how honey bee individual and group phenotypes are altered by a virus of concern, Israeli acute paralysis virus (IAPV). Using automated and manual behavioral monitoring of IAPV-inoculated individuals, we find evidence for pathogen manipulation of worker behavior by IAPV, and reveal that this effect depends on social context; that is, within versus between colony interactions. Experimental inoculation reduced social contacts between honey bee colony members, suggesting an adaptive host social immune response to diminish transmission. Parallel analyses with double-stranded RNA (dsRNA)-immunostimulated bees revealed these behaviors are part of a generalized social immune defensive response. Conversely, inoculated bees presented to groups of bees from other colonies experienced reduced aggression compared with dsRNA-immunostimulated bees, facilitating entry into susceptible colonies. This reduction was associated with a shift in cuticular hydrocarbons, the chemical signatures used by bees to discriminate colony members from intruders. These responses were specific to IAPV infection, suggestive of pathogen manipulation of the host. Emerging bee pathogens may thus shape host phenotypes to increase transmission, a strategy especially well-suited to the unnaturally high colony densities of modern apiculture. These findings demonstrate how anthropogenic changes could affect arms races between human-managed hosts and their pathogens to potentially affect global food security.

Rate variation in the evolution of non-coding DNA associated with social evolution in bees.

The evolutionary origins of eusociality represent increases in complexity from individual to caste-based, group reproduction. These behavioural transitions have been hypothesized to go hand in hand with an increased ability to regulate when and where genes are expressed. Bees have convergently evolved eusociality up to five times, providing a framework to test this hypothesis. To examine potential links between putative gene regulatory elements and social evolution, we compare alignable, non-coding sequences in 11 diverse bee species, encompassing three independent origins of reproductive division of labour and two elaborations of eusocial complexity. We find that rates of evolution in a number of non-coding sequences correlate with key social transitions in bees. Interestingly, while we find little evidence for convergent rate changes associated with independent origins of social behaviour, a number of molecular pathways exhibit convergent rate changes in conjunction with subsequent elaborations of social organization. We also present evidence that many novel non-coding regions may have been recruited alongside the origin of sociality in corbiculate bees; these loci could represent gene regulatory elements associated with division of labour within this group. Thus, our findings are consistent with the hypothesis that gene regulatory innovations are associated with the evolution of eusociality and illustrate how a thorough examination of both coding and non-coding sequence can provide a more complete understanding of the molecular mechanisms underlying behavioural evolution. This article is part of the theme issue 'Convergent evolution in the genomics era: new insights and directions'.

Draft Genome Assembly and Population Genetics of an Agricultural Pollinator, the Solitary Alkali Bee (Halictidae: Nomia melanderi).

Alkali bees (Nomia melanderi) are solitary relatives of the halictine bees, which have become an important model for the evolution of social behavior, but for which few solitary comparisons exist. These ground-nesting bees defend their developing offspring against pathogens and predators, and thus exhibit some of the key traits that preceded insect sociality. Alkali bees are also efficient native pollinators of alfalfa seed, which is a crop of major economic value in the United States. We sequenced, assembled, and annotated a high-quality draft genome of 299.6 Mbp for this species. Repetitive content makes up more than one-third of this genome, and previously uncharacterized transposable elements are the most abundant type of repetitive DNA. We predicted 10,847 protein coding genes, and identify 479 of these undergoing positive directional selection with the use of population genetic analysis based on low-coverage whole genome sequencing of 19 individuals. We found evidence of recent population bottlenecks, but no significant evidence of population structure. We also identify 45 genes enriched for protein translation and folding, transcriptional regulation, and triglyceride metabolism evolving slower in alkali bees compared to other halictid bees. These resources will be useful for future studies of bee comparative genomics and pollinator health research.

Genetic accommodation and the role of ancestral plasticity in the evolution of insect eusociality.

For over a century, biologists have proposed a role for phenotypic plasticity in evolution, providing an avenue for adaptation in addition to 'mutation-first' models of evolutionary change. According to the various versions of this idea, the ability of organisms to respond adaptively to their environment through phenotypic plasticity may lead to novel phenotypes that can be screened by natural selection. If these initially environmentally induced phenotypes increase fitness, then genetic accommodation can lead to allele frequency change, influencing the expression of those phenotypes. Despite the long history of 'plasticity-first' models, the importance of genetic accommodation in shaping evolutionary change has remained controversial - it is neither fully embraced nor completely discarded by most evolutionary biologists. We suggest that the lack of acceptance of genetic accommodation in some cases is related to a lack of information on its molecular mechanisms. However, recent reports of epigenetic transgenerational inheritance now provide a plausible mechanism through which genetic accommodation may act, and we review this research here. We also discuss current evidence supporting a role for genetic accommodation in the evolution of eusociality in social insects, which have long been models for studying the influence of the environment on phenotypic variation, and may be particularly good models for testing hypotheses related to genetic accommodation. Finally, we introduce 'eusocial engineering', a method by which novel social phenotypes are first induced by environmental modification and then studied mechanistically to understand how environmentally induced plasticity may lead to heritable changes in social behavior. We believe the time is right to incorporate genetic accommodation into models of the evolution of complex traits, armed with new molecular tools and a better understanding of non-genetic heritable elements.

Caste-biased gene expression in a facultatively eusocial bee suggests a role for genetic accommodation in the evolution of eusociality.

Developmental plasticity may accelerate the evolution of phenotypic novelty through genetic accommodation, but studies of genetic accommodation often lack knowledge of the ancestral state to place selected traits in an evolutionary context. A promising approach for assessing genetic accommodation involves using a comparative framework to ask whether ancestral plasticity is related to the evolution of a particular trait. Bees are an excellent group for such comparisons because caste-based societies (eusociality) have evolved multiple times independently and extant species exhibit different modes of eusociality. We measured brain and abdominal gene expression in a facultatively eusocial bee, Megalopta genalis, and assessed whether plasticity in this species is functionally linked to eusocial traits in other bee lineages. Caste-biased abdominal genes in M. genalis overlapped significantly with caste-biased genes in obligately eusocial bees. Moreover, caste-biased genes in M. genalis overlapped significantly with genes shown to be rapidly evolving in multiple studies of 10 bee species, particularly for genes in the glycolysis pathway and other genes involved in metabolism. These results provide support for the idea that eusociality can evolve via genetic accommodation, with plasticity in facultatively eusocial species like M. genalis providing a substrate for selection during the evolution of caste in obligately eusocial lineages.

Transcriptional, translational, and physiological signatures of undernourished honey bees (Apis mellifera) suggest a role for hormonal factors in hypopharyngeal gland degradation.

Honey bee colonies function as a superorganism, where facultatively sterile female workers perform various tasks that support the hive. Nurse workers undergo numerous anatomical and physiological changes in preparation for brood rearing, including the growth of hypopharyngeal glands (HGs). These glands produce the major protein fraction of a protein- and lipid-rich jelly used to sustain developing larvae. Pollen intake is positively correlated with HG growth, but growth in the first three days is similar regardless of diet, suggesting that initial growth is a pre-determined process while later HG development depends on nutrient availability during a critical window in early adulthood (>3 d). It is unclear whether the resultant size differences in nurse HG are simply due to growth arrest or active degradation of the tissue. To determine what processes cause such differences in HG size, we catalogued the differential expression of both gene transcripts and proteins in the HGs of 8 d old bees that were fed diets containing pollen or no pollen. 3438 genes and 367 proteins were differentially regulated due to nutrition. Of the genes and proteins differentially expressed, undernourished bees exhibited more gene and protein up-regulation compared to well-nourished bees, with the affected processes including salivary gland apoptosis, oogenesis, and hormone signaling. Protein secretion was virtually the only process up-regulated in well-nourished bees. Further assays demonstrated that inhibition of ultraspiracle, one component of the ecdysteroid receptor, in the fat body caused larger HGs. Undernourished bees also had higher acid phosphatase activity, a physiological marker of cell death, compared to well-nourished bees. These results support a connection between poor nutrition, hormonal signaling, and HG degradation.

Developmental Transcriptome for a Facultatively Eusocial Bee, Megalopta genalis.

Transcriptomes provide excellent foundational resources for mechanistic and evolutionary analyses of complex traits. We present a developmental transcriptome for the facultatively eusocial bee Megalopta genalis, which represents a potential transition point in the evolution of eusociality. A de novo transcriptome assembly of Megalopta genalis was generated using paired-end Illumina sequencing and the Trinity assembler. Males and females of all life stages were aligned to this transcriptome for analysis of gene expression profiles throughout development. Gene Ontology analysis indicates that stage-specific genes are involved in ion transport, cell-cell signaling, and metabolism. A number of distinct biological processes are upregulated in each life stage, and transitions between life stages involve shifts in dominant functional processes, including shifts from transcriptional regulation in embryos to metabolism in larvae, and increased lipid metabolism in adults. We expect that this transcriptome will provide a useful resource for future analyses to better understand the molecular basis of the evolution of eusociality and, more generally, phenotypic plasticity.

Transcriptional markers of sub-optimal nutrition in developing Apis mellifera nurse workers.

BACKGROUND: Honey bees (Apis mellifera) contribute substantially to the worldwide economy and ecosystem health as pollinators. Pollen is essential to the bee's diet, providing protein, lipids, and micronutrients. The dramatic shifts in physiology, anatomy, and behavior that accompany normal worker development are highly plastic and recent work demonstrates that development, particularly the transition from nurse to foraging roles, is greatly impacted by diet. However, the role that diet plays in the developmental transition of newly eclosed bees to nurse workers is poorly understood. To further understand honey bee nutrition and the role of diet in nurse development, we used a high-throughput screen of the transcriptome of 3 day and 8 day old worker bees fed either honey and stored pollen (rich diet) or honey alone (poor diet) within the hive. We employed a three factor (age, diet, age x diet) analysis of the transcriptome to determine whether diet affected nurse worker physiology and whether poor diet altered the developmental processes normally associated with aging. RESULTS: Substantial changes in gene expression occurred due to starvation. Diet-induced changes in gene transcription occurring in younger bees were largely a subset of those occurring in older bees, but certain signatures of starvation were only evident 8 day old workers. Of the 18,542 annotated transcripts in the A. mellifera genome, 150 transcripts exhibited differential expression due to poor diet at 3d of age compared with 17,226 transcripts that differed due to poor diet at 8d of age, and poor diet caused more frequent down-regulation of gene expression in younger bees compared to older bees. In addition, the age-related physiological changes that accompanied early adult development differed due to the diet these young adult bees were fed. More frequent down-regulation of gene expression was observed in developing bees fed a poor diet compared to those fed an adequate diet. Functional analyses also suggest that the physiological and developmental processes occurring in well-fed bees are vastly different than those occurring in pollen deprived bees. Our data support the hypothesis that poor diet causes normal age-related development to go awry. CONCLUSION: Poor nutrition has major consequences for the expression of genes underlying the physiology and age-related development of nurse worker bees. More work is certainly needed to fully understand the consequences of starvation and the complex biology of nutrition and development in this system, but the genes identified in the present study provide a starting point for understanding the consequences of poor diet and for mitigating the economic costs of colony starvation.

Microbial ecology of the hive and pollination landscape: bacterial associates from floral nectar, the alimentary tract and stored food of honey bees (Apis mellifera).

Nearly all eukaryotes are host to beneficial or benign bacteria in their gut lumen, either vertically inherited, or acquired from the environment. While bacteria core to the honey bee gut are becoming evident, the influence of the hive and pollination environment on honey bee microbial health is largely unexplored. Here we compare bacteria from floral nectar in the immediate pollination environment, different segments of the honey bee (Apis mellifera) alimentary tract, and food stored in the hive (honey and packed pollen or "beebread"). We used cultivation and sequencing to explore bacterial communities in all sample types, coupled with culture-independent analysis of beebread. We compare our results from the alimentary tract with both culture-dependent and culture-independent analyses from previous studies. Culturing the foregut (crop), midgut and hindgut with standard media produced many identical or highly similar 16S rDNA sequences found with 16S rDNA clone libraries and next generation sequencing of 16S rDNA amplicons. Despite extensive culturing with identical media, our results do not support the core crop bacterial community hypothesized by recent studies. We cultured a wide variety of bacterial strains from 6 of 7 phylogenetic groups considered core to the honey bee hindgut. Our results reveal that many bacteria prevalent in beebread and the crop are also found in floral nectar, suggesting frequent horizontal transmission. From beebread we uncovered a variety of bacterial phylotypes, including many possible pathogens and food spoilage organisms, and potentially beneficial bacteria including Lactobacillus kunkeei, Acetobacteraceae and many different groups of Actinobacteria. Contributions of these bacteria to colony health may include general hygiene, fungal and pathogen inhibition and beebread preservation. Our results are important for understanding the contribution to pollinator health of both environmentally vectored and core microbiota, and the identification of factors that may affect bacterial detection and transmission, colony food storage and disease susceptibility.

Plasticity of the worker bumblebee brain in relation to age and rearing environment.

The environment experienced during development can dramatically affect the brain, with possible implications for sensory processing, learning, and memory. Although the effects of single sensory modalities on brain development have been repeatedly explored, the additive or interactive effects of multiple modalities have been less thoroughly investigated. We asked how experience with multisensory stimuli affected brain development in the bumblebee Bombus impatiens. First, to establish the timeline of brain development during early adulthood, we estimated regional brain volumes across a range of ages. We discovered significant age-related volume changes in nearly every region of the brain. Next, to determine whether these changes were dependent upon certain environmental stimuli, we manipulated the visual and olfactory stimuli available to newly emerged bumblebee workers in a factorial manner. Newly emerged bumblebees were maintained in the presence or absence of supplemental visual and/or olfactory stimuli for 7 days, after which the volumes of several brain regions were estimated. We found that the volumes of the mushroom body lobes and calyces were larger in the absence of visual stimuli. Additionally, visual deprivation was associated with the expression of larger antennal lobes, the primary olfactory processing regions of the brain. In contrast, exposure to plant-derived olfactory stimuli did not have a significant effect on brain region volumes. This study is the first to explore the separate and interactive effects of visual and olfactory stimuli on bee brain development. Assessing the timing and sensitivity of brain development is a first step toward understanding how different rearing environments differentially affect regional brain volumes in this species. Our findings suggest that environmental factors experienced during the first week of adulthood can modify bumblebee brain development in many subtle ways.