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.

Sex-specific plasticity across generations I: Maternal and paternal effects on sons and daughters.

Intergenerational plasticity or parental effects-when parental environments alter the phenotype of future generations-can influence how organisms cope with environmental change. An intriguing, underexplored possibility is that sex-of both the parent and the offspring-plays an important role in driving the evolution of intergenerational plasticity in both adaptive and non-adaptive ways. Here, we evaluate the potential for sex-specific parental effects in a freshwater population of three-spined sticklebacks Gasterosteus aculeatus by independently and jointly manipulating maternal and paternal experiences and separately evaluating their phenotypic effects in sons versus daughters. We tested the adaptive hypothesis that daughters are more responsive to cues from their mother, whereas sons are more responsive to cues from their father. We exposed mothers, fathers or both parents to visual cues of predation risk and measured offspring antipredator traits and brain gene expression. Predator-exposed fathers produced sons that were more risk-prone, whereas predator-exposed mothers produced more anxious sons and daughters. Furthermore, maternal and paternal effects on offspring survival were non-additive: offspring with a predator-exposed father, but not two predator-exposed parents, had lower survival against live predators. There were also strong sex-specific effects on brain gene expression: exposing mothers versus fathers to predation risk activated different transcriptional profiles in their offspring, and sons and daughters strongly differed in the ways in which their brain gene expression profiles were influenced by parental experience. We found little evidence to support the hypothesis that offspring prioritize their same-sex parent's experience. Parental effects varied with both the sex of the parent and the offspring in complicated and non-additive ways. Failing to account for these sex-specific patterns (e.g. by pooling sons and daughters) would have underestimated the magnitude of parental effects. Altogether, these results draw attention to the potential for sex to influence patterns of intergenerational plasticity and raise new questions about the interface between intergenerational plasticity and sex-specific selective pressures, sexual conflict and sexual selection.

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.

Neurogenomic insights into paternal care and its relation to territorial aggression.

Motherhood is characterized by dramatic changes in brain and behavior, but less is known about fatherhood. Here we report that male sticklebacks-a small fish in which fathers provide care-experience dramatic changes in neurogenomic state as they become fathers. Some genes are unique to different stages of paternal care, some genes are shared across stages, and some genes are added to the previously acquired neurogenomic state. Comparative genomic analysis suggests that some of these neurogenomic dynamics resemble changes associated with pregnancy and reproduction in mammalian mothers. Moreover, gene regulatory analysis identifies transcription factors that are regulated in opposite directions in response to a territorial challenge versus during paternal care. Altogether these results show that some of the molecular mechanisms of parental care might be deeply conserved and might not be sex-specific, and suggest that tradeoffs between opposing social behaviors are managed at the gene regulatory level.

Personal and transgenerational cues are nonadditive at the phenotypic and molecular level.

Organisms can gain information about their environment from their ancestors, their parents or their own personal experience. 'Cue integration' models often start with the simplifying assumption that information from different sources is additive. Here, we test key assumptions and predictions of cue integration theory at both the phenotypic and molecular level in threespined sticklebacks (Gasterosteus aculeatus). We show that regardless of whether cues about predation risk were provided by their father or acquired through personal experience, sticklebacks produced the same set of predator-adapted phenotypes. Moreover, there were nonadditive effects of personal and paternal experience: animals that received cues from both sources resembled animals that received cues from a single source. A similar pattern was detected at the molecular level: there was a core set of genes that were differentially expressed in the brains of offspring regardless of whether risk was experienced by their father, themselves or both. These results provide strong support for cue integration theory because they show that cues provided by parents and personal experience are comparable at both the phenotypic and molecular level, and draw attention to the importance of nonadditive responses to multiple cues.

Temporal dynamics of neurogenomic plasticity in response to social interactions in male threespined sticklebacks.

Animals exhibit dramatic immediate behavioral plasticity in response to social interactions, and brief social interactions can shape the future social landscape. However, the molecular mechanisms contributing to behavioral plasticity are unclear. Here, we show that the genome dynamically responds to social interactions with multiple waves of transcription associated with distinct molecular functions in the brain of male threespined sticklebacks, a species famous for its behavioral repertoire and evolution. Some biological functions (e.g., hormone activity) peaked soon after a brief territorial challenge and then declined, while others (e.g., immune response) peaked hours afterwards. We identify transcription factors that are predicted to coordinate waves of transcription associated with different components of behavioral plasticity. Next, using H3K27Ac as a marker of chromatin accessibility, we show that a brief territorial intrusion was sufficient to cause rapid and dramatic changes in the epigenome. Finally, we integrate the time course brain gene expression data with a transcriptional regulatory network, and link gene expression to changes in chromatin accessibility. This study reveals rapid and dramatic epigenomic plasticity in response to a brief, highly consequential social interaction.

Stickleback embryos use ATP-binding cassette transporters as a buffer against exposure to maternally derived cortisol.

Offspring from females that experience stressful conditions during reproduction often exhibit altered phenotypes and many of these effects are thought to arise owing to increased exposure to maternal glucocorticoids. While embryos of placental vertebrates are known to regulate exposure to maternal glucocorticoids via placental steroid metabolism, much less is known about how and whether egg-laying vertebrates can control their steroid environment during embryonic development. We tested the hypothesis that threespine stickleback (Gasterosteus aculeatus) embryos can regulate exposure to maternal steroids via active efflux of maternal steroids from the egg. Embryos rapidly (within 72 h) cleared intact steroids, but blocking ATP-binding cassette (ABC) transporters inhibited cortisol clearance. Remarkably, this efflux of cortisol was sufficient to prevent a transcriptional response of embryos to exogenous cortisol. Taken together, these findings suggest that, much like their placental counterparts, developing fish embryos can actively regulate their exposure to maternal cortisol. These findings highlight the fact that even in egg-laying vertebrates, the realized exposure to maternal steroids is mediated by both maternal and embryonic processes and this has important implications for understanding how maternal stress influences offspring development.

Natural variation in brain gene expression profiles of aggressive and nonaggressive individual sticklebacks.

Within many species, some individuals are consistently more aggressive than others. We examine whether there are differences in brain gene expression between aggressive versus nonaggressive behavioural types of individuals within a natural population of male three-spined sticklebacks (Gasterosteus aculeatus). We compared gene expression profiles of aggressive male sticklebacks to nonaggressive males in four regions of the brain (brainstem, cerebellum, diencephalon and telencephalon). Relatively few genes were differentially expressed between behavioural types in telencephalon, cerebellum and diencephalon, but hundreds of genes were differentially expressed in brainstem, a brain area involved in detecting threats. Six genes that were differentially expressed in response to a territorial intrusion in a previous study were also differentially expressed between behavioural types in this study, implying primarily non-shared but some shared molecular mechanisms. Our findings offer new insights into the molecular causes and correlates of behavioural plasticity and individual variation in behaviour.

Neuromolecular responses to social challenge: common mechanisms across mouse, stickleback fish, and honey bee.

Certain complex phenotypes appear repeatedly across diverse species due to processes of evolutionary conservation and convergence. In some contexts like developmental body patterning, there is increased appreciation that common molecular mechanisms underlie common phenotypes; these molecular mechanisms include highly conserved genes and networks that may be modified by lineage-specific mutations. However, the existence of deeply conserved mechanisms for social behaviors has not yet been demonstrated. We used a comparative genomics approach to determine whether shared neuromolecular mechanisms could underlie behavioral response to territory intrusion across species spanning a broad phylogenetic range: house mouse (Mus musculus), stickleback fish (Gasterosteus aculeatus), and honey bee (Apis mellifera). Territory intrusion modulated similar brain functional processes in each species, including those associated with hormone-mediated signal transduction and neurodevelopment. Changes in chromosome organization and energy metabolism appear to be core, conserved processes involved in the response to territory intrusion. We also found that several homologous transcription factors that are typically associated with neural development were modulated across all three species, suggesting that shared neuronal effects may involve transcriptional cascades of evolutionarily conserved genes. Furthermore, immunohistochemical analyses of a subset of these transcription factors in mouse again implicated modulation of energy metabolism in the behavioral response. These results provide support for conserved genetic "toolkits" that are used in independent evolutions of the response to social challenge in diverse taxa.

Origin and evolution of protein fold designs inferred from phylogenomic analysis of CATH domain structures in proteomes.

The spatial arrangements of secondary structures in proteins, irrespective of their connectivity, depict the overall shape and organization of protein domains. These features have been used in the CATH and SCOP classifications to hierarchically partition fold space and define the architectural make up of proteins. Here we use phylogenomic methods and a census of CATH structures in hundreds of genomes to study the origin and diversification of protein architectures (A) and their associated topologies (T) and superfamilies (H). Phylogenies that describe the evolution of domain structures and proteomes were reconstructed from the structural census and used to generate timelines of domain discovery. Phylogenies of CATH domains at T and H levels of structural abstraction and associated chronologies revealed patterns of reductive evolution, the early rise of Archaea, three epochs in the evolution of the protein world, and patterns of structural sharing between superkingdoms. Phylogenies of proteomes confirmed the early appearance of Archaea. While these findings are in agreement with previous phylogenomic studies based on the SCOP classification, phylogenies unveiled sharing patterns between Archaea and Eukarya that are recent and can explain the canonical bacterial rooting typically recovered from sequence analysis. Phylogenies of CATH domains at A level uncovered general patterns of architectural origin and diversification. The tree of A structures showed that ancient structural designs such as the 3-layer (αßα) sandwich (3.40) or the orthogonal bundle (1.10) are comparatively simpler in their makeup and are involved in basic cellular functions. In contrast, modern structural designs such as prisms, propellers, 2-solenoid, super-roll, clam, trefoil and box are not widely distributed and were probably adopted to perform specialized functions. Our timelines therefore uncover a universal tendency towards protein structural complexity that is remarkable.

POLYAR, a new computer program for prediction of poly(A) sites in human sequences.

BACKGROUND: mRNA polyadenylation is an essential step of pre-mRNA processing in eukaryotes. Accurate prediction of the pre-mRNA 3'-end cleavage/polyadenylation sites is important for defining the gene boundaries and understanding gene expression mechanisms. RESULTS: 28761 human mapped poly(A) sites have been classified into three classes containing different known forms of polyadenylation signal (PAS) or none of them (PAS-strong, PAS-weak and PAS-less, respectively) and a new computer program POLYAR for the prediction of poly(A) sites of each class was developed. In comparison with polya_svm (till date the most accurate computer program for prediction of poly(A) sites) while searching for PAS-strong poly(A) sites in human sequences, POLYAR had a significantly higher prediction sensitivity (80.8% versus 65.7%) and specificity (66.4% versus 51.7%) However, when a similar sort of search was conducted for PAS-weak and PAS-less poly(A) sites, both programs had a very low prediction accuracy, which indicates that our knowledge about factors involved in the determination of the poly(A) sites is not sufficient to identify such polyadenylation regions. CONCLUSIONS: We present a new classification of polyadenylation sites into three classes and a novel computer program POLYAR for prediction of poly(A) sites/regions of each of the class. In tests, POLYAR shows high accuracy of prediction of the PAS-strong poly(A) sites, though this program's efficiency in searching for PAS-weak and PAS-less poly(A) sites is not very high but is comparable to other available programs. These findings suggest that additional characteristics of such poly(A) sites remain to be elucidated. POLYAR program with a stand-alone version for downloading is available at http://cub.comsats.edu.pk/polyapredict.htm.