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

Pathophysiology in cortico-amygdala circuits and excessive aversion processing: the role of oligodendrocytes and myelination

Stress-related psychiatric illnesses, such as major depressive disorder, anxiety and post-traumatic stress disorder, present with alterations in emotional processing, including excessive processing of negative/aversive stimuli and events. The bidirectional human/primate brain circuit comprising anterior cingulate cortex and amygdala is of fundamental importance in processing emotional stimuli, and in rodents the medial prefrontal cortex-amygdala circuit is to some extent analogous in structure and function. Here, we assess the comparative evidence for: (i) Anterior cingulate/medial prefrontal cortex<->amygdala bidirectional neural circuits as major contributors to aversive stimulus processing; (ii) Structural and functional changes in anterior cingulate cortex<->amygdala circuit associated with excessive aversion processing in stress-related neuropsychiatric disorders, and in medial prefrontal cortex<->amygdala circuit in rodent models of chronic stress-induced increased aversion reactivity; and (iii) Altered status of oligodendrocytes and their oligodendrocyte lineage cells and myelination in anterior cingulate/medial prefrontal cortex<->amygdala circuits in stress-related neuropsychiatric disorders and stress models. The comparative evidence from humans and rodents is that their respective anterior cingulate/medial prefrontal cortex<->amygdala circuits are integral to adaptive aversion processing. However, at the sub-regional level, the anterior cingulate/medial prefrontal cortex structure-function analogy is incomplete, and differences as well as similarities need to be taken into account. Structure-function imaging studies demonstrate that these neural circuits are altered in both human stress-related neuropsychiatric disorders and rodent models of stress-induced increased aversion processing. In both cases, the changes include altered white matter integrity, albeit the current evidence indicates that this is decreased in humans and increased in rodent models. At the cellular-molecular level, in both humans and rodents, the current evidence is that stress disorders do present with changes in oligodendrocyte lineage, oligodendrocytes and/or myelin in these neural circuits, but these changes are often discordant between and even within species. Nonetheless, by integrating the current comparative evidence, this review provides a timely insight into this field and should function to inform future studies-human, monkey and rodent-to ascertain whether or not the oligodendrocyte lineage and myelination are causally involved in the pathophysiology of stress-related neuropsychiatric disorders.

SPREd: a simulation-supervised neural network tool for gene regulatory network reconstruction

Summary

Reconstruction of gene regulatory networks (GRNs) from expression data is a significant open problem. Common approaches train a machine learning (ML) model to predict a gene's expression using transcription factors' (TFs') expression as features and designate important features/TFs as regulators of the gene. Here, we present an entirely different paradigm, where GRN edges are directly predicted by the ML model. The new approach, named "SPREd," is a simulation-supervised neural network for GRN inference. Its inputs comprise expression relationships (e.g. correlation, mutual information) between the target gene and each TF and between pairs of TFs. The output includes binary labels indicating whether each TF regulates the target gene. We train the neural network model using synthetic expression data generated by a biophysics-inspired simulation model that incorporates linear as well as non-linear TF-gene relationships and diverse GRN configurations. We show SPREd to outperform state-of-the-art GRN reconstruction tools GENIE3, ENNET, PORTIA, and TIGRESS on synthetic datasets with high co-expression among TFs, similar to that seen in real data. A key advantage of the new approach is its robustness to relatively small numbers of conditions (columns) in the expression matrix, which is a common problem faced by existing methods. Finally, we evaluate SPREd on real data sets in yeast that represent gold-standard benchmarks of GRN reconstruction and show it to perform significantly better than or comparably to existing methods. In addition to its high accuracy and speed, SPREd marks a first step toward incorporating biophysics principles of gene regulation into ML-based approaches to GRN reconstruction.

Availability and implementation

Data and code are available from https://github.com/iiiime/SPREd.

Isolated loss of the AUTS2 long isoform, brain-wide or targeted to Calbindin-lineage cells, generates a specific suite of brain, behavioral, and molecular pathologies

Rearrangements within the AUTS2 region are associated with a rare syndromic disorder with intellectual disability, developmental delay, and behavioral abnormalities as core features. In addition, smaller regional variants are linked to wide range of neuropsychiatric disorders, underscoring the gene's essential role in brain development. Like many essential neurodevelopmental genes, AUTS2 is large and complex, generating distinct long (AUTS2-l) and short (AUTS2-s) protein isoforms from alternative promoters. Although evidence suggests unique isoform functions, the contributions of each isoform to specific AUTS2-linked phenotypes have not been clearly resolved. Furthermore, Auts2 is widely expressed across the developing brain, but cell populations most central to disease presentation have not been determined. In this study, we focused on the specific roles of AUTS2-l in brain development, behavior, and postnatal brain gene expression, showing that brain-wide AUTS2-l ablation leads to specific subsets of the recessive pathologies associated with mutations in 3' exons (exons 8-19) that disrupt both major isoforms. We identify downstream genes that could explain expressed phenotypes including hundreds of putative direct AUTS2-l target genes. Furthermore, in contrast to 3' Auts2 mutations which lead to dominant hypoactivity, AUTS2-l loss-of-function is associated with dominant hyperactivity and repetitive behaviors, phenotypes exhibited by many human patients. Finally, we show that AUTS2-l ablation in Calbindin 1-expressing cell lineages is sufficient to yield learning/memory deficits and hyperactivity with abnormal dentate gyrus granule cell maturation, but not other phenotypic effects. These data provide new clues to in vivo AUTS2-l functions and novel information relevant to genotype-phenotype correlations in the human AUTS2 region.

SEAMoD: A fully interpretable neural network for cis-regulatory analysis of differentially expressed genes

A common way to investigate gene regulatory mechanisms is to identify differentially expressed genes using transcriptomics, find their candidate enhancers using epigenomics, and search for over-represented transcription factor (TF) motifs in these enhancers using bioinformatics tools. A related follow-up task is to model gene expression as a function of enhancer sequences and rank TF motifs by their contribution to such models, thus prioritizing among regulators. We present a new computational tool called SEAMoD that performs the above tasks of motif finding and sequence-to-expression modeling simultaneously. It trains a convolutional neural network model to relate enhancer sequences to differential expression in one or more biological conditions. The model uses TF motifs to interpret the sequences, learning these motifs and their relative importance to each biological condition from data. It also utilizes epigenomic information in the form of activity scores of putative enhancers and automatically searches for the most promising enhancer for each gene. Compared to existing neural network models of non-coding sequences, SEAMoD uses far fewer parameters, requires far less training data, and emphasizes biological interpretability. We used SEAMoD to understand regulatory mechanisms underlying the differentiation of neural stem cell (NSC) derived from mouse forebrain. We profiled gene expression and histone modifications in NSC and three differentiated cell types and used SEAMoD to model differential expression of nearly 12,000 genes with an accuracy of 81%, in the process identifying the Olig2, E2f family TFs, Foxo3, and Tcf4 as key transcriptional regulators of the differentiation process.

SPREd: A simulation-supervised neural network tool for gene regulatory network reconstruction

Reconstruction of gene regulatory networks (GRNs) from expression data is a significant open problem. Common approaches train a machine learning (ML) model to predict a gene's expression using transcription factors' (TFs') expression as features and designate important features/TFs as regulators of the gene. Here, we present an entirely different paradigm, where GRN edges are directly predicted by the ML model. The new approach, named "SPREd" is a simulation-supervised neural network for GRN inference. Its inputs comprise expression relationships (e.g., correlation, mutual information) between the target gene and each TF and between pairs of TFs. The output includes binary labels indicating whether each TF regulates the target gene. We train the neural network model using synthetic expression data generated by a biophysics-inspired simulation model that incorporates linear as well as non-linear TF-gene relationships and diverse GRN configurations. We show SPREd to outperform state-of-the-art GRN reconstruction tools GENIE3, ENNET, PORTIA and TIGRESS on synthetic datasets with high co-expression among TFs, similar to that seen in real data. A key advantage of the new approach is its robustness to relatively small numbers of conditions (columns) in the expression matrix, which is a common problem faced by existing methods. Finally, we evaluate SPREd on real data sets in yeast that represent gold standard benchmarks of GRN reconstruction and show it to perform significantly better than or comparably to existing methods. In addition to its high accuracy and speed, SPREd marks a first step towards incorporating biophysics principles of gene regulation into ML-based approaches to GRN reconstruction.

The altered hippocampal functional connectivity and serum brain-derived neurotrophic factor level predict cognitive decline in patients with knee osteoarthritis

Patients with knee osteoarthritis (KOA) often suffer from cognitive decline and increased dementia risk, but the neurobiological mechanisms remain unclear. In this study, we evaluated cognitive performance and collected brain magnetic resonance imaging (MRI) data and blood samples from cognitively normal KOA patients at baseline sessions and reevaluated their cognition after 5 years. We also collected MRI data from matched healthy controls. Results showed that KOA patients exhibited dysregulated functional connectivities between the hippocampus and thalamus/superior frontal gyrus compared with healthy controls. The altered hippocampal functional connectivities were associated with serum brain-derived neurotrophic factor (BDNF) levels and spatial expression of genes enriched in synaptic plasticity. The hippocampus-thalamus functional connectivity was significantly correlated with patients' memory scores. Moreover, the baseline hippocampus-thalamus functional connectivity and BDNF levels significantly predicted the development of cognitive decline in KOA patients in the follow-up session. Our findings provide insight into the neurobiological underpinnings of KOA and cognitive decline.

Isolated loss of the AUTS2 long isoform, brain-wide or targeted to Calbindin -lineage cells, generates a specific suite of brain, behavioral and molecular pathologies

Rearrangements within the AUTS2 region are associated with a rare syndromic disorder with intellectual disability, developmental delay and behavioral abnormalities as core features. In addition, smaller regional variants are linked to wide range of neuropsychiatric disorders, underscoring the gene's essential role in brain development. Like many essential neurodevelopmental genes, AUTS2 is large and complex, generating distinct long (AUTS2-l) and short (AUTS2-s) protein isoforms from alternative promoters. Although evidence suggests unique isoform functions, the contributions of each isoform to specific AUTS2- linked phenotypes have not been clearly resolved. Furthermore, Auts2 is widely expressed across the developing brain, but cell populations most central to disease presentation have not been determined. In this study, we focused on the specific roles of AUTS2-l in brain development, behavior, and postnatal brain gene expression, showing that brain-wide AUTS2-l ablation leads to specific subsets of the recessive pathologies associated with C-terminal mutations that disrupt both isoforms. We identify downstream genes that could explain expressed phenotypes including hundreds of putative direct AUTS2- l target genes. Furthermore, in contrast to C-terminal Auts2 mutations which lead to dominant hypoactivity, AUTS2-l loss-of-function is associated with dominant hyperactivity, a phenotype exhibited by many human patients. Finally, we show that AUTS2-l ablation in Calbindin 1 -expressing cell lineages is sufficient to yield learning/memory deficits and hyperactivity with abnormal dentate gyrus granule cell maturation, but not other phenotypic effects. These data provide new clues to in vivo AUTS2-l functions and novel information relevant to genotype-phenotype correlations in the human AUTS2 region.

Multifaceted Transcriptional Network of Estrogen-Related Receptor Alpha in Health and Disease

Estrogen-related receptors (ERRα, β and γ in mammals) are orphan members of the nuclear receptor superfamily acting as transcription factors. ERRs are expressed in several cell types and they display various functions in normal and pathological contexts. Amongst others, they are notably involved in bone homeostasis, energy metabolism and cancer progression. In contrast to other nuclear receptors, the activities of the ERRs are apparently not controlled by a natural ligand but they rely on other means such as the availability of transcriptional co-regulators. Here we focus on ERRα and review the variety of co-regulators that have been identified by various means for this receptor and their reported target genes. ERRα cooperates with distinct co-regulators to control the expression of distinct sets of target genes. This exemplifies the combinatorial specificity of transcriptional regulation that induces discrete cellular phenotypes depending on the selected coregulator. We finally propose an integrated view of the ERRα transcriptional network.

Galnt17 loss-of-function leads to developmental delay and abnormal coordination, activity, and social interactions with cerebellar vermis pathology

GALNT17 encodes a N-acetylgalactosaminyltransferase (GalNAc-T) protein specifically involved in mucin-type O-linked glycosylation of target proteins, a process important for cell adhesion, cell signaling, neurotransmitter activity, neurite outgrowth, and neurite sensing. GALNT17, also known as WBSCR17, is located at the edge of the Williams-Beuren Syndrome (WBS) critical region and adjacent to the AUTS2 locus, genomic regions associated with neurodevelopmental phenotypes that are thought to be co-regulated. Although previous data have implicated Galnt17 in neurodevelopment, the in vivo functions of this gene have not been investigated. In this study, we have analyzed behavioral, brain pathology, and molecular phenotypes exhibited by Galnt17 knockout (Galnt17-/-) mice. We show that Galnt17-/- mutants exhibit developmental neuropathology within the cerebellar vermis, along with abnormal activity, coordination, and social interaction deficits. Transcriptomic and protein analysis revealed reductions in both mucin type O-glycosylation and heparan sulfate synthesis in the developing mutant cerebellum along with disruption of pathways central to neuron differentiation, axon pathfinding, and synaptic signaling, consistent with the mutant neuropathology. These brain and behavioral phenotypes and molecular data confirm a specific role for Galnt17 in brain development and suggest new clues to factors that could contribute to phenotypes in certain WBS and AUTS2 syndrome patients.

A compound decision approach to covariance matrix estimation

Covariance matrix estimation is a fundamental statistical task in many applications, but the sample covariance matrix is sub-optimal when the sample size is comparable to or less than the number of features. Such high-dimensional settings are common in modern genomics, where covariance matrix estimation is frequently employed as a method for inferring gene networks. To achieve estimation accuracy in these settings, existing methods typically either assume that the population covariance matrix has some particular structure, for example sparsity, or apply shrinkage to better estimate the population eigenvalues. In this paper, we study a new approach to estimating high-dimensional covariance matrices. We first frame covariance matrix estimation as a compound decision problem. This motivates defining a class of decision rules and using a nonparametric empirical Bayes g-modeling approach to estimate the optimal rule in the class. Simulation results and gene network inference in an RNA-seq experiment in mouse show that our approach is comparable to or can outperform a number of state-of-the-art proposals. This article is protected by copyright. All rights reserved.

An epigenomic shift in amygdala marks the transition to maternal behaviors in alloparenting virgin female mice

Adults of many species will care for young offspring that are not their own, a phenomenon called alloparenting. However, in many cases, nonparental adults must be sensitized by repeated or extended exposures to newborns before they will robustly display parental-like behaviors. To capture neurogenomic events underlying the transition to active parental caring behaviors, we analyzed brain gene expression and chromatin profiles of virgin female mice co-housed with pregnant dams during pregnancy and after birth. After an initial display of antagonistic behaviors and a surge of defense-related gene expression, we observed a dramatic shift in the chromatin landscape specifically in amygdala of the pup-exposed virgin females compared to females co-housed with mother before birth, accompanied by a dampening of anxiety-related gene expression. This epigenetic shift coincided with hypothalamic expression of the oxytocin gene and the emergence of behaviors and gene expression patterns classically associated with maternal care. The results outline a neurogenomic program associated with dramatic behavioral changes and suggest molecular networks relevant to human postpartum mental health.

Brain transcriptome analysis reveals gene expression differences associated with dispersal behaviour between range-front and range-core populations of invasive cane toads in Australia

Understanding the mechanisms allowing invasive species to adapt to novel environments is a challenge in invasion biology. Many invaders demonstrate rapid evolution of behavioural traits involved in range expansion such as locomotor activity, exploration and risk‐taking. However, the molecular mechanisms that underpin these changes are poorly understood. In 86 years, invasive cane toads (Rhinella marina) in Australia have drastically expanded their geographic range westward from coastal Queensland to Western Australia. During their range expansion, toads have undergone extensive phenotypic changes, particularly in behaviours that enhance the toads’ dispersal ability. Common‐garden experiments have shown that some changes in behavioural traits related to dispersal are heritable. At the molecular level, it is currently unknown whether these changes in dispersal‐related behaviour are underlain by small or large differences in gene expression, nor is known the biological function of genes showing differential expression. Here, we used RNA‐seq to gain a better understanding of the molecular mechanisms underlying dispersal‐related behavioural changes. We compared the brain transcriptomes of toads from the Hawai'ian source population, as well as three distinct populations from across the Australian invasive range. We found markedly different gene expression profiles between the source population and Australian toads. By contrast, toads from across the Australian invasive range had very similar transcriptomic profiles. Yet, key genes with functions putatively related to dispersal behaviour showed differential expression between populations located at each end of the invasive range. These genes could play an important role in the behavioural changes characteristic of range expansion in Australian cane toads.

Estrogen-related Receptor Alpha (ERRα) is Required for PGC-1α-dependent Gene Expression in the Mouse Brain

Deficiency in peroxisome proliferator-activated receptor gamma coactivator 1-alpha. (PGC-1α) expression or function is implicated in numerous neurological and psychiatric disorders. PGC-1α is required for the expression of genes involved in synchronous neurotransmitter release, axonal integrity, and metabolism, especially in parvalbumin-positive interneurons. As a transcriptional coactivator, PGC-1α requires transcription factors to specify cell-type-specific gene programs; while much is known about these factors in peripheral tissues, it is unclear if PGC-1α utilizes these same factors in neurons. Here, we identified putative transcription factors controlling PGC-1α-dependent gene expression in the brain using bioinformatics and then validated the role of the top candidate in a knockout mouse model. We transcriptionally profiled cells overexpressing PGC-1α and searched for over-represented binding motifs in the promoters of upregulated genes. Binding sites of the estrogen-related receptor (ERR) family of transcription factors were enriched, and blockade of ERRα attenuated PGC-1α-mediated induction of mitochondrial and synaptic genes in cell culture. Localization in the mouse brain revealed enrichment of ERRα expression in parvalbumin-expressing neurons with tight correlation of expression with PGC-1α across brain regions. In ERRα null mice, PGC-1α-dependent genes were reduced in multiple regions, including neocortex, hippocampus, and cerebellum, though not to the extent observed in PGC-1α null mice. Behavioral assessment revealed ambulatory hyperactivity in response to amphetamine and impairments in sensorimotor gating without the overt motor impairment characteristic of PGC-1α null mice. These data suggest that ERRα is required for normal levels of expression of PGC-1α-dependent genes in neurons but that additional factors may be involved in their regulation.

Dynamic neurogenomic responses to social interactions and dominance outcomes in female paper wasps

Social interactions have large effects on individual physiology and fitness. In the immediate sense, social stimuli are often highly salient and engaging. Over longer time scales, competitive interactions often lead to distinct social ranks and differences in physiology and behavior. Understanding how initial responses lead to longer-term effects of social interactions requires examining the changes in responses over time. Here we examined the effects of social interactions on transcriptomic signatures at two times, at the end of a 45-minute interaction and 4 hours later, in female Polistes fuscatus paper wasp foundresses. Female P. fuscatus have variable facial patterns that are used for visual individual recognition, so we separately examined the transcriptional dynamics in the optic lobe and the non-visual brain. Results demonstrate much stronger transcriptional responses to social interactions in the non-visual brain compared to the optic lobe. Differentially regulated genes in response to social interactions are enriched for memory-related transcripts. Comparisons between winners and losers of the encounters revealed similar overall transcriptional profiles at the end of an interaction, which significantly diverged over the course of 4 hours, with losers showing changes in expression levels of genes associated with aggression and reproduction in paper wasps. On nests, subordinate foundresses are less aggressive, do more foraging and lay fewer eggs compared to dominant foundresses and we find losers shift expression of many genes in the non-visual brain, including vitellogenin, related to aggression, worker behavior, and reproduction within hours of losing an encounter. These results highlight the early neurogenomic changes that likely contribute to behavioral and physiological effects of social status changes in a social insect.

Neural and Molecular Mechanisms of Biological Embedding of Social Interactions

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

A unique neurogenomic state emerges after aggressive confrontations in males of the fish Betta splendens

Territorial defense involves frequent aggressive confrontations with competitors, but little is known about how brain-transcriptomic profiles change between individuals competing for territory establishment. Our previous study elucidated that when two fish Betta splendens males interact, transcriptomes across their brains synchronize in a way that reflects a mutual assessment process between them at the gene expression level. Here we aim to evaluate how the brain-transcriptomic profiles of opponents change immediately after shifting their social status (i.e., the winner/loser has emerged) and 30 min after this shift. We showed that changes in the expression of certain genes are unique to different fighting stages and the expression patterns of certain genes are transiently or persistently changed across all fighting stages. These brain transcriptomic responses are in accordance with behavioral changes across the fight. Strikingly, the specificity of the brain-transcriptomic synchronization of a pair during fighting was gradually lost after fighting ceased, leading to the emergence of a basal neurogenomic state in which the changes in gene expression were reduced to minimum and consistent across all individuals. This state shares common characteristics with the hibernation state that animals adopt to minimize their metabolic rates to save energy. Interestingly, expression changes for genes related to metabolism, autism spectrum disorder, and long-term memory still differentiated losers from winners. Together, the fighting system using male B. splendens provides a promising platform for investigating neurogenomic states of aggression in vertebrates.

Transcriptional Regulation of ROS Homeostasis by the ERR Subfamily of Nuclear Receptors

Reactive oxygen species (ROS) such as superoxide anion (O₂^•-) and hydrogen peroxide (H₂O₂) are generated endogenously by processes such as mitochondrial oxidative phosphorylation, or they may arise from exogenous sources like bacterial invasion. ROS can be beneficial (oxidative eustress) as signaling molecules but also harmful (oxidative distress) to cells when ROS levels become unregulated in response to physiological, pathological or pharmacological insults. Indeed, abnormal ROS levels have been shown to contribute to the etiology of a wide variety of diseases. Transcriptional control of metabolic genes is a crucial mechanism to coordinate ROS homeostasis. Therefore, a better understanding of how ROS metabolism is regulated by specific transcription factors can contribute to uncovering new therapeutic strategies. A large body of work has positioned the estrogen-related receptors (ERRs), transcription factors belonging to the nuclear receptor superfamily, as not only master regulators of cellular energy metabolism but, most recently, of ROS metabolism. Herein, we will review the role played by the ERRs as transcriptional regulators of ROS generation and antioxidant mechanisms and also as ROS sensors. We will assess how the control of ROS homeostasis by the ERRs can be linked to physiology and disease and the possible contribution of manipulating ERR activity in redox medicine.

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 α-Tocopherol Concentration and Stereoisomer Profile Alter Hippocampal Gene Expression in Weanling Mice

BACKGROUND

Alpha-tocopherol (αT), the bioactive constituent of vitamin E, is essential for fertility and neurological development. Synthetic αT (8 stereoisomers; all rac-αT) is added to infant formula at higher concentrations than natural αT (RRR-αT only) to adjust for bio-potency differences, but its effects on brain development are poorly understood.

OBJECTIVES

The objective was to determine the impact of bio-potency-adjusted dietary all rac-αT versus RRR-αT, fed to dams, on the hippocampal gene expression in weanling mice.

METHODS

Male/female pairs of C57BL/6J mice were fed AIN 93-G containing RRR-αT (NAT) or all rac-αT (SYN) at 37.5 or 75 IU/kg (n = 10/group) throughout gestation and lactation. Male pups were euthanized at 21 days. Half the brain was evaluated for the αT concentration and stereoisomer distribution. The hippocampus was dissected from the other half, and RNA was extracted and sequenced. Milk αT was analyzed in separate dams.

RESULTS

A total of 797 differentially expressed genes (DEGs) were identified in the hippocampi across the 4 dietary groups, at a false discovery rate of 10%. Comparing the NAT-37.5 group to the NAT-75 group or the SYN-37.5 group to the SYN-75 group, small differences in brain αT concentrations (10%; P < 0.05) led to subtle changes (<10%) in gene expression of 600 (NAT) or 487 genes (SYN), which were statistically significant. Marked differences in brain αT stereoisomer profiles (P < 0.0001) had a small effect on fewer genes (NAT-37.5 vs. SYN-37.5, 179; NAT-75 vs. SYN-75, 182). Most of the DEGs were involved in transcription regulation and synapse formation. A network analysis constructed around known vitamin E interacting proteins (VIPs) revealed a group of 32 DEGs between NAT-37.5 vs. SYN-37.5, explained by expression of the gene for the VIP, protein kinase C zeta (Pkcz).

CONCLUSIONS

In weanling mouse hippocampi, a network of genes involved in transcription regulation and synapse formation was differentially affected by dam diet αT concentration and source: all rac-αT or RRR-αT.

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.

Biological embedding of experience: A primer on epigenetics

Biological embedding occurs when life experience alters biological processes to affect later life health and well-being. Although extensive correlative data exist supporting the notion that epigenetic mechanisms such as DNA methylation underlie biological embedding, causal data are lacking. We describe specific epigenetic mechanisms and their potential roles in the biological embedding of experience. We also consider the nuanced relationships between the genome, the epigenome, and gene expression. Our ability to connect biological embedding to the epigenetic landscape in its complexity is challenging and complicated by the influence of multiple factors. These include cell type, age, the timing of experience, sex, and DNA sequence. Recent advances in molecular profiling and epigenome editing, combined with the use of comparative animal and human longitudinal studies, should enable this field to transition from correlative to causal analyses.

Geometric regulation of histone state directs melanoma reprogramming

Malignant melanoma displays a high degree of cellular plasticity during disease progression. Signals in the tumor microenvironment are believed to influence melanoma plasticity through changes in the epigenetic state to guide dynamic differentiation and de-differentiation. Here we uncover a relationship between geometric features at perimeter regions of melanoma aggregates, and reprogramming to a stem cell-like state through histone marks H3K4Me2 and H3K9Ac. Using an in vitro tumor microengineering approach, we find spatial enrichment of these histone modifications with concurrent expression of stemness markers. The epigenetic modifier PRDM14 overlaps with H3K9Ac and shows elevated expression in cells along regions of perimeter curvature. siRNA knockdown of PRDM14 abolishes the MIC phenotype suggesting a role in regulating melanoma heterogeneity. Our results suggest mechanotransduction at the periphery of melanoma aggregates may orchestrate the activity of epigenetic modifiers to regulate histone state, cellular plasticity, and tumorigenicity.

Junmin Lee et al. study the role of geometric features at the perimeter regions of melanoma aggregates in programming stem cell-like state through histone marks. They use a tumor microengineering approach in vitro and report a spatial enrichment of histone modifications with stemness markers. Their work uncovers a mechanotransduction signaling that regulates epigenetic modifiers to regulate tumorigenicity.

Behavioral and brain- transcriptomic synchronization between the two opponents of a fighting pair of the fish Betta splendens

Conspecific male animals fight for resources such as food and mating opportunities but typically stop fighting after assessing their relative fighting abilities to avoid serious injuries. Physiologically, how the fighting behavior is controlled remains unknown. Using the fighting fish Betta splendens, we studied behavioral and brain-transcriptomic changes during the fight between the two opponents. At the behavioral level, surface-breathing, and biting/striking occurred only during intervals between mouth-locking. Eventually, the behaviors of the two opponents became synchronized, with each pair showing a unique behavioral pattern. At the physiological level, we examined the expression patterns of 23,306 brain transcripts using RNA-sequencing data from brains of fighting pairs after a 20-min (D20) and a 60-min (D60) fight. The two opponents in each D60 fighting pair showed a strong gene expression correlation, whereas those in D20 fighting pairs showed a weak correlation. Moreover, each fighting pair in the D60 group showed pair-specific gene expression patterns in a grade of membership analysis (GoM) and were grouped as a pair in the heatmap clustering. The observed pair-specific individualization in brain-transcriptomic synchronization (PIBS) suggested that this synchronization provides a physiological basis for the behavioral synchronization. An analysis using the synchronized genes in fighting pairs of the D60 group found genes enriched for ion transport, synaptic function, and learning and memory. Brain-transcriptomic synchronization could be a general phenomenon and may provide a new cornerstone with which to investigate coordinating and sustaining social interactions between two interacting partners of vertebrates.

Agonistic encounters induce changes in the brain and behavior, but their underlying molecular mechanisms remain poorly understood. The fighting fish Betta splendens are small freshwater fish that are well known for their aggressiveness and are widely used to study aggression. Here, by measuring aggressive behavior displays (bite/strike/surface-breathing) between two opponents during fighting, we demonstrate that the two opponents in each fighting pair showed similar fighting configurations by influencing each other. In addition, we compared brain gene expression between opponents and showed synchronization of gene expression within a fighting pair, leading to pair-specific synchronization in genes associated with ion transport, synapse function, and learning and memory. This study presents the possibility that similar behaviors in pairs of animals under similar conditions may trigger synchronizing waves of transcription between the individuals, providing a hint to support the idea that fighting behaviors contain cooperative aspects at the molecular level.

The Role of Chromatin Accessibility in cis-Regulatory Evolution

Transcription factor (TF) binding is determined by sequence as well as chromatin accessibility. Although the role of accessibility in shaping TF-binding landscapes is well recorded, its role in evolutionary divergence of TF binding, which in turn can alter cis-regulatory activities, is not well understood. In this work, we studied the evolution of genome-wide binding landscapes of five major TFs in the core network of mesoderm specification, between Drosophila melanogaster and Drosophila virilis, and examined its relationship to accessibility and sequence-level changes. We generated chromatin accessibility data from three important stages of embryogenesis in both Drosophila melanogaster and Drosophila virilis and recorded conservation and divergence patterns. We then used multivariable models to correlate accessibility and sequence changes to TF-binding divergence. We found that accessibility changes can in some cases, for example, for the master regulator Twist and for earlier developmental stages, more accurately predict binding change than is possible using TF-binding motif changes between orthologous enhancers. Accessibility changes also explain a significant portion of the codivergence of TF pairs. We noted that accessibility and motif changes offer complementary views of the evolution of TF binding and developed a combined model that captures the evolutionary data much more accurately than either view alone. Finally, we trained machine learning models to predict enhancer activity from TF binding and used these functional models to argue that motif and accessibility-based predictors of TF-binding change can substitute for experimentally measured binding change, for the purpose of predicting evolutionary changes in enhancer activity.

A Mouse Mutation That Dysregulates Neighboring Galnt17 and Auts2 Genes Is Associated with Phenotypes Related to the Human AUTS2 Syndrome

AUTS2 was originally discovered as the gene disrupted by a translocation in human twins with Autism spectrum disorder, intellectual disability, and epilepsy. Since that initial finding, AUTS2-linked mutations and variants have been associated with a very broad array of neuropsychiatric disorders, sugg esting that AUTS2 is required for fundamental steps of neurodevelopment. However, genotype-phenotype correlations in this region are complicated, because most mutations could also involve neighboring genes. Of particular interest is the nearest downstream neighbor of AUTS2, GALNT17, which encodes a brain-expressed N-acetylgalactosaminyltransferase of unknown brain function. Here we describe a mouse (Mus musculus) mutation, T(5G2;8A1)GSO (abbreviated 16Gso), a reciprocal translocation that breaks between Auts2 and Galnt17 and dysregulates both genes. Despite this complex regulatory effect, 16Gso homozygotes model certain human AUTS2-linked phenotypes very well. In addition to abnormalities in growth, craniofacial structure, learning and memory, and behavior, 16Gso homozygotes display distinct pathologies of the cerebellum and hippocampus that are similar to those associated with autism and other types of AUTS2-linked neurological disease. Analyzing mutant cerebellar and hippocampal transcriptomes to explain this pathology, we identified disturbances in pathways related to neuron and synapse maturation, neurotransmitter signaling, and cellular stress, suggesting possible cellular mechanisms. These pathways, coupled with the translocation's selective effects on Auts2 isoforms and coordinated dysregulation of Galnt17, suggest novel hypotheses regarding the etiology of the human "AUTS2 syndrome" and the wide array of neurodevelopmental disorders linked to variance in this genomic region.

Subregion-specific rules govern the distribution of neuronal immediate-early gene induction

The induction of immediate-early gene (IEG) expression in brain nuclei in response to an experience is necessary for the formation of long-term memories. Additionally, the rapid dynamics of IEG induction and decay motivates the common use of IEG expression as markers for identification of neuronal assemblies ("ensembles") encoding recent experience. However, major gaps remain in understanding the rules governing the distribution of IEGs within neuronal assemblies. Thus, the extent of correlation between coexpressed IEGs, the cell specificity of IEG expression, and the spatial distribution of IEG expression have not been comprehensively studied. To address these gaps, we utilized quantitative multiplexed single-molecule fluorescence in situ hybridization (smFISH) and measured the expression of IEGs (Arc, Egr2, and Nr4a1) within spiny projection neurons (SPNs) in the dorsal striatum of mice following acute exposure to cocaine. Exploring the relevance of our observations to other brain structures and stimuli, we also analyzed data from a study of single-cell RNA sequencing of mouse cortical neurons. We found that while IEG expression is graded, the expression of multiple IEGs is tightly correlated at the level of individual neurons. Interestingly, we observed that region-specific rules govern the induction of IEGs in SPN subtypes within striatal subdomains. We further observed that IEG-expressing assemblies form spatially defined clusters within which the extent of IEG expression correlates with cluster size. Together, our results suggest the existence of IEG-expressing neuronal "superensembles," which are associated in spatial clusters and characterized by coherent and robust expression of multiple IEGs.

Social history and exposure to pathogen signals modulate social status effects on gene regulation in rhesus macaques

Social experience is an important predictor of disease susceptibility and survival in humans and other social mammals. Chronic social stress is thought to generate a proinflammatory state characterized by elevated antibacterial defenses and reduced investment in antiviral defense. Here we manipulated long-term social status in female rhesus macaques to show that social subordination alters the gene expression response to ex vivo bacterial and viral challenge. As predicted by current models, bacterial lipopolysaccharide polarizes the immune response such that low status corresponds to higher expression of genes in NF-κB-dependent proinflammatory pathways and lower expression of genes involved in the antiviral response and type I IFN signaling. Counter to predictions, however, low status drives more exaggerated expression of both NF-κB- and IFN-associated genes after cells are exposed to the viral mimic Gardiquimod. Status-driven gene expression patterns are linked not only to social status at the time of sampling, but also to social history (i.e., past social status), especially in unstimulated cells. However, for a subset of genes, we observed interaction effects in which females who fell in rank were more strongly affected by current social status than those who climbed the social hierarchy. Taken together, our results indicate that the effects of social status on immune cell gene expression depend on pathogen exposure, pathogen type, and social history-in support of social experience-mediated biological embedding in adulthood, even in the conventionally memory-less innate immune system.

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.

Compared to motherhood, the molecular changes associated with fatherhood are less understood. Here, the authors investigate gene expression changes associated with paternal care in male stickleback fish, and compare them with patterns in territorial aggression.

Hydrogel-Tissue Chemistry: Principles and Applications

Over the past five years, a rapidly developing experimental approach has enabled high-resolution and high-content information retrieval from intact multicellular animal (metazoan) systems. New chemical and physical forms are created in the hydrogel-tissue chemistry process, and the retention and retrieval of crucial phenotypic information regarding constituent cells and molecules (and their joint interrelationships) are thereby enabled. For example, rich data sets defining both single-cell-resolution gene expression and single-cell-resolution activity during behavior can now be collected while still preserving information on three-dimensional positioning and/or brain-wide wiring of those very same neurons-even within vertebrate brains. This new approach and its variants, as applied to neuroscience, are beginning to illuminate the fundamental cellular and chemical representations of sensation, cognition, and action. More generally, reimagining metazoans as metareactants-or positionally defined three-dimensional graphs of constituent chemicals made available for ongoing functionalization, transformation, and readout-is stimulating innovation across biology and medicine.

Social Environment and Epigenetics

Our social environment, from the microscopic to the macro-social, affects us for the entirety of our lives. One integral line of research to examine how interpersonal and societal environments can get "under the skin" is through the lens of epigenetics. Epigenetic mechanisms are adaptations made to our genome in response to our environment which include tags placed on and removed from the DNA itself to how our DNA is packaged, affecting how our genes are read, transcribed, and interact. These tags are affected by social environments and can persist over time; this may aid us in responding to experiences and exposures, both the enriched and the disadvantageous. From memory formation to immune function, the experience-dependent plasticity of epigenetic modifications to micro- and macro-social environments may contribute to the process of learning from comfort, pain, and stress to better survive in whatever circumstances life has in store.

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.

L-Lactate Regulates the Expression of Synaptic Plasticity and Neuroprotection Genes in Cortical Neurons: A Transcriptome Analysis

Lactate, a product of aerobic glycolysis in astrocytes, is required for memory formation and consolidation, and has recently emerged as a signaling molecule for neurons and various cell types in peripheral tissues. In particular lactate stimulates mRNA expression of a few plasticity-related genes. Here, we describe a RNA-seq study that unravels genome-wide transcriptomic responses to this energy metabolite in cortical neurons. Our results show that mRNA expression of 20 immediate-early genes involved in the MAPK signaling pathway and in synaptic plasticity were increased by more than twofold following 1 h of lactate stimulation. This effect was dependent on NMDA receptor (NMDAR) activity since it was prevented by pre-treatment with MK-801. Comparison with published datasets showed that a significant proportion of genes modulated by lactate were similarly regulated by a stimulation protocol activating specifically synaptic NMDARs known to result in upregulation of pro-survival and downregulation of pro-death genes. Remarkably, transcriptional responses to lactate were reproduced by NADH (for 74 of the 113 genes, FDR < 0.05), suggesting a redox-dependent mechanism of action. Longer-term gene expression changes observed after 6 h of lactate treatment affected genes involved in regulating neuronal excitability and genes coding for proteins localized at synapses. Gene set enrichment analyses performed with ranked lists of expressed genes revealed effects on molecular functions involved in epigenetic modulation, and on processes relevant to sleep physiology and behavioral phenotypes such as anxiety and hyperactivity. Overall, these results strengthen the notion that lactate effectively regulates activity-dependent and synaptic genes, and highlight new signaling effects of lactate in plasticity and neuroprotection.

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.

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.

Emerging Roles of Estrogen-Related Receptors in the Brain: Potential Interactions with Estrogen Signaling

In addition to their well-known role in the female reproductive system, estrogens can act in the brain to regulate a wide range of behaviors and physiological functions in both sexes. Over the past few decades, genetically modified animal models have greatly increased our knowledge about the roles of estrogen receptor (ER) signaling in the brain in behavioral and physiological regulations. However, less attention has been paid to the estrogen-related receptors (ERRs), the members of orphan nuclear receptors whose sequences are homologous to ERs but lack estrogen-binding ability. While endogenous ligands of ERRs remain to be determined, they seemingly share transcriptional targets with ERs and their expression can be directly regulated by ERs through the estrogen-response element embedded within the regulatory region of the genes encoding ERRs. Despite the broad expression of ERRs in the brain, we have just begun to understand the fundamental roles they play at molecular, cellular, and circuit levels. Here, we review recent research advancement in understanding the roles of ERs and ERRs in the brain, with particular emphasis on ERRs, and discuss possible cross-talk between ERs and ERRs in behavioral and physiological regulations.

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.

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.