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

Cross-species functional modules link proteostasis to human normal aging

The evolutionarily conserved nature of the few well-known anti-aging interventions that affect lifespan, such as caloric restriction, suggests that aging-related research in model organisms is directly relevant to human aging. Since human lifespan is a complex trait, a systems-level approach will contribute to a more comprehensive understanding of the underlying aging landscape. Here, we integrate evolutionary and functional information of normal aging across human and model organisms at three levels: gene-level, process-level, and network-level. We identify evolutionarily conserved modules of normal aging across diverse taxa, and notably show proteostasis to be conserved in normal aging. Additionally, we find that mechanisms related to protein quality control network are enriched for genes harboring genetic variants associated with 22 age-related human traits and associated to caloric restriction. These results demonstrate that a systems-level approach, combined with evolutionary conservation, allows the detection of candidate aging genes and pathways relevant to human normal aging.

Biogenic amines and activity levels alter the neural energetic response to aggressive social cues in the honey bee Apis mellifera

Mitochondrial activity is highly dynamic in the healthy brain, and it can reflect both the signaling potential and the signaling history of neural circuits. Recent studies spanning invertebrates to mammals have highlighted a role for neural mitochondrial dynamics in learning and memory processes as well as behavior. In the current study, we investigate the interplay between biogenic amine signaling and neural energetics in the honey bee, Apis mellifera. In this species, aggressive behaviors are regulated by neural energetic state and biogenic amine titers, but it is unclear how these mechanisms are linked to impact behavioral expression. We show that brain mitochondrial number is highest in aggression-relevant brain regions and in individual bees that are most responsive to aggressive cues, emphasizing the importance of energetics in modulating this phenotype. We also show that the neural energetic response to alarm pheromone, an aggression inducing social cue, is activity dependent, modulated by the "fight or flight" insect neurotransmitter octopamine. Two other neuroactive compounds known to cause variation in aggression, dopamine, and serotonin, also modulate neural energetic state in aggression-relevant regions of the brain. However, the effects of these compounds on respiration at baseline and following alarm pheromone exposure are distinct, suggesting unique mechanisms underlying variation in mitochondrial respiration in these circuits. These results motivate new explanations for the ways in which biogenic amines alter sensory perception in the context of aggression. Considering neural energetics improves predictions about the regulation of complex and context-dependent behavioral phenotypes.

Comparative Analysis of Brain and Fat Body Gene Splicing Patterns in the Honey Bee, Apis mellifera

RNA-seq has proven to be a powerful tool to unravel various aspects of the transcriptome, especially the quantification of alternative splicing (AS) that leads to isoform diversity. The honey bee (Apis mellifera) is an important model organism for studying the molecular underpinnings of behavioral plasticity and social behavior, and recent RNA-seq studies of honey bees have revealed AS patterns and their regulation by DNA methylation. However, tissue-specific AS patterns have not been fully explored. In this paper, we characterized AS patterns in two different honey bee tissue types, and also explored their conservation and regulation. We used the RNA-seq data from brain and fat body to improve the existing models of honey bee genes and identified tissue-specific AS patterns. We found that AS genes show high conservation between honey bee and Drosophila melanogaster. We also confirmed and extended previous findings of a correlation between gene body DNA methylation and AS patterns, providing further support for the role of DNA methylation in regulating AS. In addition, our analysis suggests distinct functional roles for tissue-specific alternatively spliced genes. Taken together, our work provides new insights into the conservation and dynamics of AS patterns across different tissue types.

Repeated Evolution of Asexuality Involves Convergent Gene Expression Changes

Asexual reproduction has evolved repeatedly from sexual ancestors across a wide range of taxa. Whereas the costs and benefits associated with asexuality have received considerable attention, the molecular changes underpinning the evolution of asexual reproduction remain relatively unexplored. In particular, it is completely unknown whether the repeated evolution of asexual phenotypes involves similar molecular changes, as previous studies have focused on changes occurring in single lineages. Here, we investigate the extent of convergent gene expression changes across five independent transitions to asexuality in stick insects. We compared gene expression of asexual females to females of close sexual relatives in whole-bodies, reproductive tracts, and legs. We identified a striking amount of convergent gene expression change (up to 8% of genes), greatly exceeding that expected by chance. Convergent changes were also tissue-specific, and most likely driven by selection for functional changes. Genes showing convergent changes in the reproductive tract were associated with meiotic spindle formation and centrosome organization. These genes are particularly interesting as they can influence the production of unreduced eggs, a key barrier to asexual reproduction. Changes in legs and whole-bodies were likely involved in female sexual trait decay, with enrichment in terms such as sperm-storage and pigmentation. By identifying changes occurring across multiple independent transitions to asexuality, our results provide a rare insight into the molecular basis of asexual phenotypes and suggest that the evolutionary path to asexuality is highly constrained, requiring repeated changes to the same key genes.

Evolutionary diversity as a catalyst for biological discovery

The tremendous diversity of animal behaviors has inspired generations of scientists from an array of biological disciplines. To complement investigations of ecological and evolutionary factors contributing to behavioral evolution, modern sequencing, gene editing, computational and neuroscience tools now provide a means to discover the proximate mechanisms upon which natural selection acts to generate behavioral diversity. Social behaviors are motivated behaviors that can differ tremendously between closely related species, suggesting phylogenetic plasticity in their underlying biological mechanisms. In addition, convergent evolution has repeatedly given rise to similar forms of social behavior and mating systems in distantly related species. Social behavioral divergence and convergence provides an entry point for understanding the neurogenetic mechanisms contributing to behavioral diversity. We argue that the greatest strides in discovering mechanisms contributing to social behavioral diversity will be achieved through integration of interdisciplinary comparative approaches with modern tools in diverse species systems. We review recent advances and future potential for discovering mechanisms underlying social behavioral variation; highlighting patterns of social behavioral evolution, oxytocin and vasopressin neuropeptide systems, genetic/transcriptional "toolkits," modern experimental tools, and alternative species systems, with particular emphasis on Microtine rodents and Lake Malawi cichlid fishes.

Changes of gene expression but not cytosine methylation are associated with male parental care reflecting behavioural state, social context and individual flexibility

Behaviour is often a front line response to changing environments. Recent studies show behavioural changes are associated with changes of gene expression; however, these studies have primarily focused on discrete behavioural states. We build on these studies by addressing additional contexts that produce qualitatively similar behavioural changes. We measured levels of gene expression and cytosine methylation, which is hypothesized to regulate the transcriptional architecture of behavioural transitions, within the brain during male parental care of the burying beetle Nicrophorus vespilloides in a factorial design. Male parenting is a suitably plastic behaviour because although male N. vespilloides typically do not provide direct care (i.e. feed offspring) when females are present, levels of feeding by a male equivalent to the female can be induced by removing the female. We examined three different factors: behavioural state (caring versus non-caring), social context (with or without a female mate) and individual flexibility (if a male switched to direct care after his mate was removed). The greatest number of differentially expressed genes were associated with behavioural state, followed by social context and individual flexibility. Cytosine methylation was not associated with changes of gene expression in any of the factors. Our results suggest a hierarchical association between gene expression and the different factors, but that this process is not controlled by cytosine methylation. Our results further suggest that the extent a behaviour is transient plays an underappreciated role in determining its underpinning molecular mechanisms.

Conserved transcriptomic profiles underpin monogamy across vertebrates

Social monogamy, typically characterized by the formation of a pair bond, increased territorial defense, and often biparental care, has independently evolved multiple times in animals. Despite the independent evolutionary origins of monogamous mating systems, several homologous brain regions and neuropeptides and their receptors have been shown to play a conserved role in regulating social affiliation and parental care, but little is known about the neuromolecular mechanisms underlying monogamy on a genomic scale. Here, we compare neural transcriptomes of reproductive males in monogamous and nonmonogamous species pairs of Peromyscus mice, Microtus voles, parid songbirds, dendrobatid frogs, and Xenotilapia species of cichlid fishes. We find that, while evolutionary divergence time between species or clades did not explain gene expression similarity, characteristics of the mating system correlated with neural gene expression patterns, and neural gene expression varied concordantly across vertebrates when species transition to monogamy. Our study provides evidence of a universal transcriptomic mechanism underlying the evolution of monogamy in vertebrates.

Dominance rank-associated gene expression is widespread, sex-specific, and a precursor to high social status in wild male baboons

In humans and other hierarchical species, social status is tightly linked to variation in health and fitness-related traits. Experimental manipulations of social status in female rhesus macaques suggest that this relationship is partially explained by status effects on immune gene regulation. However, social hierarchies are established and maintained in different ways across species: While some are based on kin-directed nepotism, others emerge from direct physical competition. We investigated how this variation influences the relationship between social status and immune gene regulation in wild baboons, where hierarchies in males are based on fighting ability but female hierarchies are nepotistic. We measured rank-related variation in gene expression levels in adult baboons of both sexes at baseline and in response to ex vivo stimulation with the bacterial endotoxin lipopolysaccharide (LPS). We identified >2,000 rank-associated genes in males, an order of magnitude more than in females. In males, high status predicted increased expression of genes involved in innate immunity and preferential activation of the NF-κB-mediated proinflammatory pathway, a pattern previously associated with low status in female rhesus macaques. Using Mendelian randomization, we reconcile these observations by demonstrating that high status-associated gene expression patterns are precursors, not consequences, of high social status in males, in support of the idea that physiological condition determines who attains high rank. Together, our work provides a test of the relationship between social status and immune gene regulation in wild primates. It also emphasizes the importance of social context in shaping the relationship between social status and immune function.

Social status alters chromatin accessibility and the gene regulatory response to glucocorticoid stimulation in rhesus macaques

Low social status is an important predictor of disease susceptibility and mortality risk in humans and other social mammals. These effects are thought to stem in part from dysregulation of the glucocorticoid (GC)-mediated stress response. However, the molecular mechanisms that connect low social status and GC dysregulation to downstream health outcomes remain elusive. Here, we used an in vitro GC challenge to investigate the consequences of experimentally manipulated social status (i.e., dominance rank) for immune cell gene regulation in female rhesus macaques, using paired control and GC-treated peripheral blood mononuclear cell samples. We show that social status not only influences immune cell gene expression but also chromatin accessibility at hundreds of regions in the genome. Social status effects on gene expression were less pronounced following GC treatment than under control conditions. In contrast, social status effects on chromatin accessibility were stable across conditions, resulting in an attenuated relationship between social status, chromatin accessibility, and gene expression after GC exposure. Regions that were more accessible in high-status animals and regions that become more accessible following GC treatment were enriched for a highly concordant set of transcription factor binding motifs, including motifs for the GC receptor cofactor AP-1. Together, our findings support the hypothesis that social status alters the dynamics of GC-mediated gene regulation and identify chromatin accessibility as a mechanism involved in social stress-driven GC resistance. More broadly, they emphasize the context-dependent nature of social status effects on gene regulation and implicate epigenetic remodeling of chromatin accessibility as a contributing factor.

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.

Conserved Genes Underlie Phenotypic Plasticity in an Incipiently Social Bee

Despite a strong history of theoretical work on the mechanisms of social evolution, relatively little is known of the molecular genetic changes that accompany transitions from solitary to eusocial forms. Here, we provide the first genome of an incipiently social bee that shows both solitary and social colony organization in sympatry, the Australian carpenter bee Ceratina australensis. Through comparative analysis, we provide support for the role of conserved genes and cis-regulation of gene expression in the phenotypic plasticity observed in nest-sharing, a rudimentary form of sociality. Additionally, we find that these conserved genes are associated with caste differences in advanced eusocial species, suggesting these types of mechanisms could pave the molecular pathway from solitary to eusocial living. Genes associated with social nesting in this species show signatures of being deeply conserved, in contrast to previous studies in other bees showing novel and faster-evolving genes are associated with derived sociality. Our data provide support for the idea that the earliest social transitions are driven by changes in gene regulation of deeply conserved genes.

Honey bee neurogenomic responses to affiliative and agonistic social interactions

Social interactions can be divided into two categories, affiliative and agonistic. How neurogenomic responses reflect these opposing valences is a central question in the biological embedding of experience. To address this question, we exposed honey bees to a queen larva, which evokes nursing, an affiliative alloparenting interaction, and measured the transcriptomic response of the mushroom body brain region at different times after exposure. Hundreds of genes were differentially expressed at distinct time points, revealing a dynamic temporal patterning of the response. Comparing these results to our previously published research on agonistic aggressive interactions, we found both shared and unique transcriptomic responses to each interaction. The commonly responding gene set was enriched for nuclear receptor signaling, the set specific to nursing was enriched for olfaction and neuron differentiation, and the set enriched for aggression was enriched for cytoskeleton, metabolism, and chromosome organization. Whole brain histone profiling after the affiliative interaction revealed few changes in chromatin accessibility, suggesting that the transcriptomic changes derive from already accessible areas of the genome. Although only one stimulus of each type was studied, we suggest that elements of the observed transcriptomic responses reflect molecular encoding of stimulus valence, thus priming individuals for future encounters. This hypothesis is supported by behavioral analyses showing that bees responding to either the affiliative or agonistic stimulus exhibited a higher probability of repeating the same behavior but a lower probability of performing the opposite behavior. These findings add to our understanding of the biological embedding at the molecular level.

The effect of maternal care on gene expression and DNA methylation in a subsocial bee

Developmental plasticity describes the influence of environmental factors on phenotypic variation. An important mediator of developmental plasticity in many animals is parental care. Here, we examine the consequences of maternal care on offspring after the initial mass provisioning of brood in the small carpenter bee, Ceratina calcarata. Removal of the mother during larval development leads to increased aggression and avoidance in adulthood. This corresponds with changes in expression of over one thousand genes, alternative splicing of hundreds of genes, and significant changes to DNA methylation. We identify genes related to metabolic and neuronal functions that may influence developmental plasticity and aggression. We observe no genome-wide association between differential DNA methylation and differential gene expression or splicing, though indirect relationships may exist between these factors. Our results provide insight into the gene regulatory context of DNA methylation in insects and the molecular avenues through which variation in maternal care influences developmental plasticity.

Cross-species systems analysis of evolutionary toolkits of neurogenomic response to social challenge

Social challenges like territorial intrusions evoke behavioral responses in widely diverging species. Recent work has showed that evolutionary "toolkits"-genes and modules with lineage-specific variations but deep conservation of function-participate in the behavioral response to social challenge. Here, we develop a multispecies computational-experimental approach to characterize such a toolkit at a systems level. Brain transcriptomic responses to social challenge was probed via RNA-seq profiling in three diverged species-honey bees, mice and three-spined stickleback fish-following a common methodology, allowing fair comparisons across species. Data were collected from multiple brain regions and multiple time points after social challenge exposure, achieving anatomical and temporal resolution substantially greater than previous work. We developed statistically rigorous analyses equipped to find homologous functional groups among these species at the levels of individual genes, functional and coexpressed gene modules, and transcription factor subnetworks. We identified six orthogroups involved in response to social challenge, including groups represented by mouse genes Npas4 and Nr4a1, as well as common modulation of systems such as transcriptional regulators, ion channels, G-protein-coupled receptors and synaptic proteins. We also identified conserved coexpression modules enriched for mitochondrial fatty acid metabolism and heat shock that constitute the shared neurogenomic response. Our analysis suggests a toolkit wherein nuclear receptors, interacting with chaperones, induce transcriptional changes in mitochondrial activity, neural cytoarchitecture and synaptic transmission after social challenge. It shows systems-level mechanisms that have been repeatedly co-opted during evolution of analogous behaviors, thus advancing the genetic toolkit concept beyond individual genes.

The evolution of a series of behavioral traits is associated with autism-risk genes in cavefish

Background

An essential question in evolutionary biology is whether shifts in a set of polygenic behaviors share a genetic basis across species. Such a behavioral shift is seen in the cave-dwelling Mexican tetra, Astyanax mexicanus. Relative to surface-dwelling conspecifics, cavefish do not school (asocial), are hyperactive and sleepless, adhere to a particular vibration stimulus (imbalanced attention), behave repetitively, and show elevated stress hormone levels. Interestingly, these traits largely overlap with the core symptoms of human autism spectrum disorder (ASD), raising the possibility that these behavioral traits are underpinned by a similar set of genes (i.e. a repeatedly used suite of genes).

Result

Here, we explored whether modification of ASD-risk genes underlies cavefish evolution. Transcriptomic analyses revealed that > 58.5% of 3152 cavefish orthologs to ASD-risk genes are significantly up- or down-regulated in the same direction as genes in postmortem brains from ASD patients. Enrichment tests suggest that ASD-risk gene orthologs in A. mexicanus have experienced more positive selection than other genes across the genome. Notably, these positively selected cavefish ASD-risk genes are enriched for pathways involved in gut function, inflammatory diseases, and lipid/energy metabolism, similar to symptoms that frequently coexist in ASD patients. Lastly, ASD drugs mitigated cavefish’s ASD-like behaviors, implying shared aspects of neural processing.

Conclusion

Overall, our study indicates that ASD-risk genes and associated pathways (especially digestive, immune and metabolic pathways) may be repeatedly used for shifts in polygenic behaviors across evolutionary time.

Electronic supplementary material

The online version of this article (10.1186/s12862-018-1199-9) contains supplementary material, which is available to authorized users.

Defense against territorial intrusion is associated with DNA methylation changes in the honey bee brain

BACKGROUND

Aggression is influenced by individual variation in temperament as well as behavioral plasticity in response to adversity. DNA methylation is stably maintained over time, but also reversible in response to specific environmental conditions, and may thus be a neuromolecular regulator of both of these processes. A previous study reported DNA methylation differences between aggressive Africanized and gentle European honey bees. We investigated whether threat-induced aggression altered DNA methylation profiles in the honey bee brain in response to a behavioral stimulus (aggression-provoking intruder bee or inert control). We sampled five minutes and two hours after stimulus exposure to examine the effect of time on epigenetic profiles of aggression.

RESULTS

There were DNA methylation differences between aggressive and control bees for individual cytosine-guanine dinucleotides (CpGs) across the genome. Eighteen individual CpG sites showed significant difference between aggressive and control bees 120 min post stimulus. For clusters of CpGs, we report four genomic regions differentially methylated between aggressive and control bees at the 5-min time point, and 50 regions differentially methylated at the120-minute time point following intruder exposure. Differential methylation occurred at genes involved in neural plasticity, chromatin remodeling and hormone signaling. Additionally, there was a significant overlap of differential methylation with previously published epigenetic differences that distinguish aggressive Africanized and gentle European honey bees, suggesting an evolutionarily conserved use of brain DNA methylation in the regulation of aggression. Lastly, we identified individually statistically suggestive CpGs that as a group were significantly associated with differentially expressed genes underlying aggressive behavior and also co-localize with binding sites of transcription factors involved in neuroplasticity or neurodevelopment.

CONCLUSIONS

There were DNA methylation differences in the brain associated with response to an intruder. These differences increased in number a few hours after the initial exposure and overlap with previously reported aggression-associated genes and neurobiologically relevant transcription factor binding sites. Many DNA methylation differences that occurred in association with the expression of aggression in real time also exist between Africanized bees and European bees, suggesting an evolutionarily conserved role for epigenetic regulation in aggressive behavior.

Dissecting the Transcriptional Patterns of Social Dominance across Teleosts

In many species, under varying ecological conditions, social interactions among individuals result in the formation of dominance hierarchies. Despite general similarities, there are robust differences among dominance hierarchies across species, populations, environments, life stages, sexes, and individuals. Understanding the proximate mechanisms underlying the variation is an important step toward understanding the evolution of social behavior. However, physiological changes associated with dominance, such as gonadal maturation and somatic growth, often complicate efforts to identify the specific underlying mechanisms. Traditional gene expression analyses are useful for generating candidate gene lists, but are biased by choice of significance cut-offs and difficult to use for between-study comparisons. In contrast, complementary analysis tools allow one to both test a priori hypotheses and generate new hypotheses. Here we employ a meta-analysis of high-throughput expression profiling experiments to investigate the gene expression patterns that underlie mechanisms and evolution of behavioral social phenotypes. Specifically, we use a collection of datasets on social dominance in fish across social contexts, sex, and species. Using experimental manipulation to produce female dominance hierarchies in the cichlid Astatotilapia burtoni, heralded as a genomic model of social dominance, we generate gene lists, and assess molecular gene modules. In the dominant female gene expression profile, we demonstrate a strong pattern of up-regulation of genes previously identified as having male-biased expression and furthermore, compare expression biases between male and female dominance phenotypes. Using a threshold-free approach to identify correlation throughout ranked gene lists, we query previously published datasets associated with maternal behavior, alternative reproductive tactics, cooperative breeding, and sex-role reversal to describe correlations among these various neural gene expression profiles associated with different instances of social dominance. These complementary approaches capitalize on the high-throughput gene expression profiling from similar behavioral phenotypes in order to address the mechanisms associated with social dominance behavioral phenotypes.

Candidate genes for cooperation and aggression in the social wasp Polistes dominula

Cooperation and aggression are ubiquitous in social groups, and the genetic mechanisms underlying these behaviours are of great interest for understanding how social group formation is regulated and how it evolves. In this study, we used a candidate gene approach to investigate the patterns of expression of key genes for cooperation and aggression in the brain of a primitively eusocial wasp, Polistes dominula, during colony founding, when multiple foundresses can join the same nest and establish subtle hierarchies of dominance. We used a comparative approach to select candidate genes for cooperation and aggression looking at two previously published studies on global gene expression in wasps and ants. We tested the expression of these genes in P. dominula wasps that were either displaying aggressive behaviour (dominant and single foundresses) or cooperation (subordinate foundresses and workers) towards nestmates. One gene in particular, the egg yolk protein vitellogenin, known for its reproductive role in insects, displayed patterns of expression that strongly matched wasp social rank. We characterize the genomic context of vitellogenin by building a head co-expression gene network for P. dominula, and we discuss a potential role for vitellogenin as a mediator of social interactions in wasps.

Genomic tools for behavioural ecologists to understand repeatable individual differences in behaviour

Behaviour is a key interface between an animal's genome and its environment. Repeatable individual differences in behaviour have been extensively documented in animals, but the molecular underpinnings of behavioural variation among individuals within natural populations remain largely unknown. Here, we offer a critical review of when molecular techniques may yield new insights, and we provide specific guidance on how and whether the latest tools available are appropriate given different resources, system and organismal constraints, and experimental designs. Integrating molecular genetic techniques with other strategies to study the proximal causes of behaviour provides opportunities to expand rapidly into new avenues of exploration. Such endeavours will enable us to better understand how repeatable individual differences in behaviour have evolved, how they are expressed and how they can be maintained within natural populations of animals.

Advancing behavioural genomics by considering timescale

Animal behavioural traits often covary with gene expression, pointing towards a genomic constraint on organismal responses to environmental cues. This pattern highlights a gap in our understanding of the time course of environmentally responsive gene expression, and moreover, how these dynamics are regulated. Advances in behavioural genomics explore how gene expression dynamics are correlated with behavioural traits that range from stable to highly labile. We consider the idea that certain genomic regulatory mechanisms may predict the timescale of an environmental effect on behaviour. This temporally minded approach could inform both organismal and evolutionary questions ranging from the remediation of early life social trauma to understanding the evolution of trait plasticity.

Brain mitochondrial bioenergetics change with rapid and prolonged shifts in aggression in the honey bee, Apis mellifera

Neuronal function demands high-level energy production, and as such, a decline in mitochondrial respiration characterizes brain injury and disease. A growing number of studies, however, link brain mitochondrial function to behavioral modulation in non-diseased contexts. In the honey bee, we show for the first time that an acute social interaction, which invokes an aggressive response, may also cause a rapid decline in brain mitochondrial bioenergetics. The degree and speed of this decline has only been previously observed in the context of brain injury. Furthermore, in the honey bee, age-related increases in aggressive tendency are associated with increased baseline brain mitochondrial respiration, as well as increased plasticity in response to metabolic fuel type in vitro. Similarly, diet restriction and ketone body feeding, which commonly enhance mammalian brain mitochondrial function in vivo, cause increased aggression. Thus, even in normal behavioral contexts, brain mitochondria show a surprising degree of variation in function over both rapid and prolonged timescales, with age predicting both baseline function and plasticity in function. These results suggest that mitochondrial function is integral to modulating aggression-related neuronal signaling. We hypothesize that variation in function reflects mitochondrial calcium buffering activity, and that shifts in mitochondrial function signal to the neuronal soma to regulate gene expression and neural energetic state. Modulating brain energetic state is emerging as a critical component of the regulation of behavior in non-diseased contexts.

Vasopressin and the Neurogenetics of Parental Care

Making robust connections between genetic variation, neurophysiology, and social behavior remains a challenge. A study by Bendesky et al. (2017) tackles this challenge by dissecting the genetic architecture of parental care in deer mice to discover an important contribution of vasopressin signaling to the evolution of nest building.

Deep evolutionary conservation of autism-related genes

E. O. Wilson proposed in Sociobiology that similarities between human and animal societies reflect common mechanistic and evolutionary roots. When introduced in 1975, this controversial hypothesis was beyond science's ability to test. We used genomic analyses to determine whether superficial behavioral similarities in humans and the highly social honey bee reflect common molecular mechanisms. Here, we report that gene expression signatures for individual bees unresponsive to various salient social stimuli are significantly enriched for autism spectrum disorder-related genes. These signatures occur in the mushroom bodies, a high-level integration center of the insect brain. Furthermore, our finding of enrichment was unique to autism spectrum disorders; brain gene expression signatures from other honey bee behaviors do not show this enrichment, nor do datasets from other human behavioral and health conditions. These results demonstrate deep conservation for genes associated with a human social pathology and individual differences in insect social behavior, thus providing an example of how comparative genomics can be used to test sociobiological theory.

Divergence in Morris Water Maze-Based Cognitive Performance under Chronic Stress Is Associated with the Hippocampal Whole Transcriptomic Modification in Mice

Individual susceptibility determines the magnitude of stress effects on cognitive function. The hippocampus, a brain region of memory consolidation, is vulnerable to stressful environments, and the impact of stress on hippocampus may determine individual variability in cognitive performance. Therefore, the purpose of this study was to define the relationship between the divergence in spatial memory performance under chronically unpredictable stress and an associated transcriptomic alternation in hippocampus, the brain region of spatial memory consolidation. Multiple strains of BXD (B6 × D2) recombinant inbred mice went through a 4-week chronic variable stress (CVS) paradigm, and the Morris water maze (MWM) test was conducted during the last week of CVS to assess hippocampal-dependent spatial memory performance and grouped animals into low and high performing groups based on the cognitive performance. Using hippocampal whole transcriptome RNA-sequencing data, differential expression, PANTHER analysis, WGCNA, Ingenuity's upstream regulator analysis in the Ingenuity Pathway Analysis® and phenotype association analysis were conducted. Our data identified multiple genes and pathways that were significantly associated with chronic stress-associated cognitive modification and the divergence in hippocampal dependent memory performance under chronic stress. Biological pathways associated with memory performance following chronic stress included metabolism, neurotransmitter and receptor regulation, immune response and cellular process. The Ingenuity's upstream regulator analysis identified 247 upstream transcriptional regulators from 16 different molecule types. Transcripts predictive of cognitive performance under high stress included genes that are associated with a high occurrence of Alzheimer's and cognitive impairments (e.g., Ncl, Eno1, Scn9a, Slc19a3, Ncstn, Fos, Eif4h, Copa, etc.). Our results show that the variable effects of chronic stress on the hippocampal transcriptome are related to the ability to complete the MWM task and that the modulations of specific pathways are indicative of hippocampal dependent memory performance. Thus, the divergence in spatial memory performance following chronic stress is related to the unique pattern of gene expression within the hippocampus.

Stress Induced Polarization of Immune-Neuroendocrine Phenotypes in Gallus gallus

Immune-neuroendocrine phenotypes (INPs) stand for population subgroups differing in immune-neuroendocrine interactions. While mammalian INPs have been characterized thoroughly in rats and humans, avian INPs were only recently described in Coturnix coturnix (quail). To assess the scope of this biological phenomenon, herein we characterized INPs in Gallus gallus (a domestic hen strain submitted to a very long history of strong selective breeding pressure) and evaluated whether a social chronic stress challenge modulates the individuals' interplay affecting the INP subsets and distribution. Evaluating plasmatic basal corticosterone, interferon-γ and interleukin-4 concentrations, innate/acquired leukocyte ratio, PHA-P skin-swelling and induced antibody responses, two opposite INP profiles were found: LEWIS-like (15% of the population) and FISCHER-like (16%) hens. After chronic stress, an increment of about 12% in each polarized INP frequency was found at expenses of a reduction in the number of birds with intermediate responses. Results show that polarized INPs are also a phenomenon occurring in hens. The observed inter-individual variation suggest that, even after a considerable selection process, the population is still well prepared to deal with a variety of immune-neuroendocrine challenges. Stress promoted disruptive effects, leading to a more balanced INPs distribution, which represents a new substrate for challenging situations.

XGSA: A statistical method for cross-species gene set analysis

MOTIVATION

Gene set analysis is a powerful tool for determining whether an experimentally derived set of genes is statistically significantly enriched for genes in other pre-defined gene sets, such as known pathways, gene ontology terms, or other experimentally derived gene sets. Current gene set analysis methods do not facilitate comparing gene sets across different organisms as they do not explicitly deal with homology mapping between species. There lacks a systematic investigation about the effect of complex gene homology on cross-species gene set analysis.

RESULTS

In this study, we show that not accounting for the complex homology structure when comparing gene sets in two species can lead to false positive discoveries, especially when comparing gene sets that have complex gene homology relationships. To overcome this bias, we propose a straightforward statistical approach, called XGSA, that explicitly takes the cross-species homology mapping into consideration when doing gene set analysis. Simulation experiments confirm that XGSA can avoid false positive discoveries, while maintaining good statistical power compared to other ad hoc approaches for cross-species gene set analysis. We further demonstrate the effectiveness of XGSA with two real-life case studies that aim to discover conserved or species-specific molecular pathways involved in social challenge and vertebrate appendage regeneration.

AVAILABILITY AND IMPLEMENTATION

The R source code for XGSA is available under a GNU General Public License at http://github.com/VCCRI/XGSA CONTACT: jho@victorchang.edu.au.

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.

Social interactions provoke changes in the brain and behavior but their underlying molecular mechanisms remain obscure. Male sticklebacks are small fish whose fitness depends on their ability to defend a territory. Here, by measuring the time course of gene expression in response to a territorial challenge in two brain regions, we show that a single brief territorial intrusion provoked waves of gene expression that persisted for hours afterwards, with waves of transcription associated with distinct biological processes. Moreover, a single territorial challenge caused dramatic changes to the epigenome. Changes in chromatin accessibility corresponded to changes in gene expression, and to the activity of transcription factors operating within gene regulatory networks. This study reveals rapid and dramatic epigenomic plasticity in response to a brief, highly consequential social interaction. These results suggest that meaningful social interactions (even brief ones) can provoke waves of transcription and changes to the epigenome which lead to changes in neural functioning, and those changes are a mechanism by which animals update their assessment of their social world.

Genetic editing of the androgen receptor contributes to impaired male courtship behavior in zebrafish

Elucidating the genes that contribute to behavioral variation has become an important endeavor in behavioral studies. While advances in genomics have narrowed down the list of candidate genes, functional validation of them has lagged behind, partly because of challenges associated with rapid gene manipulations. Consequently, few studies have demonstrated causal genetic changes linked to behaviors. The 'gene editing revolution' has offered unprecedented opportunities to investigate candidate genes responsible for critical behaviors. Here, we edited the androgen receptor gene (AR), which is associated with male reproductive behavior in zebrafish, using TAL effector nucleases (TALENs), and tested whether modifications at the AR impacted courtship during mating trials. We reveal that males lacking AR courted females significantly less, showing reduced levels of stereotypic behaviors. Consistent with previous studies, disrupting androgen mechanisms can lead to behavioral changes with potential fitness consequences. Our study highlights the possibility of genetically altering a reproductive behavior, further solidifying the link between genotype and behavior.

Transcriptional regulatory dynamics drive coordinated metabolic and neural response to social challenge in mice

Agonistic encounters are powerful effectors of future behavior, and the ability to learn from this type of social challenge is an essential adaptive trait. We recently identified a conserved transcriptional program defining the response to social challenge across animal species, highly enriched in transcription factor (TF), energy metabolism, and developmental signaling genes. To understand the trajectory of this program and to uncover the most important regulatory influences controlling this response, we integrated gene expression data with the chromatin landscape in the hypothalamus, frontal cortex, and amygdala of socially challenged mice over time. The expression data revealed a complex spatiotemporal patterning of events starting with neural signaling molecules in the frontal cortex and ending in the modulation of developmental factors in the amygdala and hypothalamus, underpinned by a systems-wide shift in expression of energy metabolism-related genes. The transcriptional signals were correlated with significant shifts in chromatin accessibility and a network of challenge-associated TFs. Among these, the conserved metabolic and developmental regulator ESRRA was highlighted for an especially early and important regulatory role. Cell-type deconvolution analysis attributed the differential metabolic and developmental signals in this social context primarily to oligodendrocytes and neurons, respectively, and we show that ESRRA is expressed in both cell types. Localizing ESRRA binding sites in cortical chromatin, we show that this nuclear receptor binds both differentially expressed energy-related and neurodevelopmental TF genes. These data link metabolic and neurodevelopmental signaling to social challenge, and identify key regulatory drivers of this process with unprecedented tissue and temporal resolution.

Sequential social experiences interact to modulate aggression but not brain gene expression in the honey bee (Apis mellifera)

BACKGROUND

In highly structured societies, individuals behave flexibly and cooperatively in order to achieve a particular group-level outcome. However, even in social species, environmental inputs can have long lasting effects on individual behavior, and variable experiences can even result in consistent individual differences and constrained behavioral flexibility. Despite the fact that such constraints on behavior could have implications for behavioral optimization at the social group level, few studies have explored how social experiences accumulate over time, and the mechanistic basis of these effects. In the current study, I evaluate how sequential social experiences affect individual and group level aggressive phenotypes, and individual brain gene expression, in the highly social honey bee (Apis mellifera). To do this, I combine a whole colony chronic predator disturbance treatment with a lab-based manipulation of social group composition.

RESULTS

Compared to the undisturbed control, chronically disturbed individuals show lower aggression levels overall, but also enhanced behavioral flexibility in the second, lab-based social context. Disturbed bees display aggression levels that decline with increasing numbers of more aggressive, undisturbed group members. However, group level aggressive phenotypes are similar regardless of the behavioral tendencies of the individuals that make up the group, suggesting a combination of underlying behavioral tendency and negative social feedback influences the aggressive behaviors displayed, particularly in the case of disturbed individuals. An analysis of brain gene expression showed that aggression related biomarker genes reflect an individual's disturbance history, but not subsequent social group experience or behavioral outcomes.

CONCLUSIONS

In highly social animals with collective behavioral phenotypes, social context may mask underlying variation in individual behavioral tendencies. Moreover, gene expression patterns may reflect behavioral tendency, while behavioral outcomes are further regulated by social cues perceived in real-time.

Neuroendocrine-immune circuits, phenotypes, and interactions

Multidirectional interactions among the immune, endocrine, and nervous systems have been demonstrated in humans and non-human animal models for many decades by the biomedical community, but ecological and evolutionary perspectives are lacking. Neuroendocrine-immune interactions can be conceptualized using a series of feedback loops, which culminate into distinct neuroendocrine-immune phenotypes. Behavior can exert profound influences on these phenotypes, which can in turn reciprocally modulate behavior. For example, the behavioral aspects of reproduction, including courtship, aggression, mate selection and parental behaviors can impinge upon neuroendocrine-immune interactions. One classic example is the immunocompetence handicap hypothesis (ICHH), which proposes that steroid hormones act as mediators of traits important for female choice while suppressing the immune system. Reciprocally, neuroendocrine-immune pathways can promote the development of altered behavioral states, such as sickness behavior. Understanding the energetic signals that mediate neuroendocrine-immune crosstalk is an active area of research. Although the field of psychoneuroimmunology (PNI) has begun to explore this crosstalk from a biomedical standpoint, the neuroendocrine-immune-behavior nexus has been relatively underappreciated in comparative species. The field of ecoimmunology, while traditionally emphasizing the study of non-model systems from an ecological evolutionary perspective, often under natural conditions, has focused less on the physiological mechanisms underlying behavioral responses. This review summarizes neuroendocrine-immune interactions using a comparative framework to understand the ecological and evolutionary forces that shape these complex physiological interactions.

Putting the "Biology" Back into "Neurobiology": The Strength of Diversity in Animal Model Systems for Neuroscience Research

Current trends in neuroscience research have moved toward a reliance on rodent animal models to study most aspects of brain function. Such laboratory-reared animals are highly inbred, have been disengaged from their natural environments for generations and appear to be of limited predictive value for successful clinical outcomes. In this Perspective article, we argue that research on a rich diversity of animal model systems is fundamental to new discoveries in evolutionarily conserved core physiological and molecular mechanisms that are the foundation of human brain function. Analysis of neural circuits across phyla will reveal general computational solutions that form the basis for adaptive behavioral responses. Further, we stress that development of ethoexperimental approaches to improve our understanding of behavioral nuance will help to realign our research strategies with therapeutic goals and improve the translational validity of specific animal models. Finally, we suggest that neuroscience has a role in environmental conservation of habitat and fauna that will preserve and protect the ecological settings that drive species-specific behavioral adaptations. A rich biodiversity will enhance our understanding of human brain function and lead in unpredicted directions for development of therapeutic treatments for neurological disorders.

Conservation in Mammals of Genes Associated with Aggression-Related Behavioral Phenotypes in Honey Bees

The emerging field of sociogenomics explores the relations between social behavior and genome structure and function. An important question is the extent to which associations between social behavior and gene expression are conserved among the Metazoa. Prior experimental work in an invertebrate model of social behavior, the honey bee, revealed distinct brain gene expression patterns in African and European honey bees, and within European honey bees with different behavioral phenotypes. The present work is a computational study of these previous findings in which we analyze, by orthology determination, the extent to which genes that are socially regulated in honey bees are conserved across the Metazoa. We found that the differentially expressed gene sets associated with alarm pheromone response, the difference between old and young bees, and the colony influence on soldier bees, are enriched in widely conserved genes, indicating that these differences have genomic bases shared with many other metazoans. By contrast, the sets of differentially expressed genes associated with the differences between African and European forager and guard bees are depleted in widely conserved genes, indicating that the genomic basis for this social behavior is relatively specific to honey bees. For the alarm pheromone response gene set, we found a particularly high degree of conservation with mammals, even though the alarm pheromone itself is bee-specific. Gene Ontology identification of human orthologs to the strongly conserved honey bee genes associated with the alarm pheromone response shows overrepresentation of protein metabolism, regulation of protein complex formation, and protein folding, perhaps associated with remodeling of critical neural circuits in response to alarm pheromone. We hypothesize that such remodeling may be an adaptation of social animals to process and respond appropriately to the complex patterns of conspecific communication essential for social organization.

Routine handling methods affect behaviour of three-spined sticklebacks in a novel test of anxiety

Fish are increasingly popular subjects in behavioural and neurobiological research. It is therefore important that they are housed and handled appropriately to ensure good welfare and reliable scientific findings, and that species-appropriate behavioural tests (e.g. of cognitive/affective states) are developed. Routine handling of captive animals may cause physiological stress responses that lead to anxiety-like states (e.g. increased perception of danger). In fish, these may be particularly pronounced when handling during tank-to-tank transfer involves removal from water into air. Here we develop and use a new combined scototaxis (preference for dark over light areas) and novel-tank-diving test, alongside conventional open-field and novel-object tests, to measure the effects of transferring three-spined sticklebacks (Gasterosteus aculeatus) between tanks using a box or net (in and out of water respectively). Preference tests for dark over light areas confirmed the presence of scototaxis in this species. Open-field and novel-object tests failed to detect any significant differences between net and box-handled fish. However, the combined diving and scototaxis detected consistent differences between the treatments. Net-handled fish spent less time on the dark side of the tank, less time in the bottom third, and kept a greater distance from the 'safe' bottom dark area than box-handled fish. Possible explanations for this reduction in anxiety-like behaviour in net-handled fish are discussed. The combined diving and scototaxis test may be a sensitive and taxon-appropriate method for measuring anxiety-like states in fish.

Characterizing gene sets using discriminative random walks with restart on heterogeneous biological networks

Motivation: Analysis of co-expressed gene sets typically involves testing for enrichment of different annotations or ‘properties’ such as biological processes, pathways, transcription factor binding sites, etc., one property at a time. This common approach ignores any known relationships among the properties or the genes themselves. It is believed that known biological relationships among genes and their many properties may be exploited to more accurately reveal commonalities of a gene set. Previous work has sought to achieve this by building biological networks that combine multiple types of gene–gene or gene–property relationships, and performing network analysis to identify other genes and properties most relevant to a given gene set. Most existing network-based approaches for recognizing genes or annotations relevant to a given gene set collapse information about different properties to simplify (homogenize) the networks.

Results: We present a network-based method for ranking genes or properties related to a given gene set. Such related genes or properties are identified from among the nodes of a large, heterogeneous network of biological information. Our method involves a random walk with restarts, performed on an initial network with multiple node and edge types that preserve more of the original, specific property information than current methods that operate on homogeneous networks. In this first stage of our algorithm, we find the properties that are the most relevant to the given gene set and extract a subnetwork of the original network, comprising only these relevant properties. We then re-rank genes by their similarity to the given gene set, based on a second random walk with restarts, performed on the above subnetwork. We demonstrate the effectiveness of this algorithm for ranking genes related to Drosophila embryonic development and aggressive responses in the brains of social animals.

Availability and Implementation: DRaWR was implemented as an R package available at veda.cs.illinois.edu/DRaWR.

Contact: blatti@illinois.edu

Supplementary information:Supplementary data are available at Bioinformatics online.

Brain regions and molecular pathways responding to food reward type and value in honey bees

The ability of honey bees to evaluate differences in food type and value is crucial for colony success, but these assessments are made by individuals who bring food to the hive, eating little, if any, of it themselves. We tested the hypothesis that responses to food type (pollen or nectar) and value involve different subsets of brain regions, and genes responsive to food. mRNA in situ hybridization of c-jun revealed that brain regions responsive to differences in food type were mostly different from regions responsive to differences in food value, except those dorsal and lateral to the mushroom body calyces, which responded to all three. Transcriptomic profiles of the mushroom bodies generated by RNA sequencing gave the following results: (1) responses to differences in food type or value included a subset of molecular pathways involved in the response to food reward; (2) genes responsive to food reward, food type and food value were enriched for (the Gene Ontology categories) mitochondrial and endoplasmic reticulum activity; (3) genes responsive to only food and food type were enriched for regulation of transcription and translation; and (4) genes responsive to only food and food value were enriched for regulation of neuronal signaling. These results reveal how activities necessary for colony survival are channeled through the reward system of individual honey bees.

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.

A fluorescence in situ hybridization (FISH) protocol for stickleback tissue

BACKGROUND

Threespine stickleback are an important model for behaviour and evolutionary studies. A growing number of quantitative trait loci (QTL) and gene expression studies are identifying genes related to ecologically important traits in sticklebacks. In order to visualize the expression of candidate genes, we developed a fluorescence in situ hybridization (FISH) protocol.

METHODS

We present a protocol for FISH on fresh or flash-frozen dissected tissue, using either cryo- or paraffin embedding. The protocol covers probe design guidelines and synthesis, sample embedding, sectioning, and the hybridization process. The protocol is optimized for brain tissue. Key steps for modifying the protocol for other tissues are noted.

RESULTS

The FISH protocol resulted in specific labelling under all combinations of dissection and embedding conditions. Paraffin embedding preserved morphology better than cryo-embedding. We provide representative results showing the expression of glial fibrillary acidic protein (GFAP), oxytocin receptor (OXTR), and tyrosine hydroxylase (TH) in the brain.

Integrating molecular mechanisms into quantitative genetics to understand consistent individual differences in behavior

It is now well appreciated that individual animals behave differently from one another and that individual differences in behaviors-personality differences-are maintained through time and across situations. Quantitative genetics has emerged as a conceptual basis for understanding the key ingredients of personality: (co)variation and plasticity. However, the results from quantitative genetic analyses are often divorced from underlying molecular or other proximate mechanisms. This disconnect has the potential to impede an integrated understanding of behavior and is a disconnect present throughout evolutionary ecology. Here we discuss some of the main conceptual connections between personality and quantitative genetics, the relationship of both with genomic tools, and areas that require integration. With its consideration of both trait variation and plasticity, the study of animal personality offers new opportunities to incorporate molecular mechanisms into both the trait partitioning and reaction norm frameworks provided by quantitative genetics.