#PAGE_PARAMS# #ADS_HEAD_SCRIPTS# #MICRODATA#

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


Authors: Trieu-Duc Vu aff001;  Yuki Iwasaki aff002;  Shuji Shigenobu aff007;  Akiko Maruko aff001;  Kenshiro Oshima aff001;  Erica Iioka aff002;  Chao-Li Huang aff008;  Takashi Abe aff009;  Satoshi Tamaki aff010;  Yi-Wen Lin aff005;  Chih-Kuan Chen aff004;  Mei-Yeh Lu aff004;  Masaru Hojo aff005;  Hao-Ven Wang aff005;  Shun-Fen Tzeng aff005;  Hao-Jen Huang aff005;  Akio Kanai aff010;  Takashi Gojobori aff011;  Tzen-Yuh Chiang aff005;  H. Sunny Sun aff012;  Wen-Hsiung Li aff004;  Norihiro Okada aff001
Authors place of work: School of Pharmacy, Kitasato University, Tokyo, Japan aff001;  Foundation for Advancement of International Science, Tsukuba, Japan aff002;  Life Sciences and Biotechnology Dept, Tokyo Institute of Technology, Tokyo, Japan aff003;  Biodiversity Research Center, Academia Sinica, Taipei, Taiwan aff004;  Department of Life Sciences, National Cheng Kung University, Tainan, Taiwan aff005;  Nagahama Institute of Bio-Science and Technology, Nagahama, Japan aff006;  National Institute for Basic Biology, Okazaki, Japan aff007;  Institute of Tropical Plant Sciences, National Cheng Kung University, Tainan, Taiwan aff008;  Graduate School of Science and Technology, Niigata University, Niigata, Japan aff009;  Institute for Advanced Biosciences, Keio University, Yamagata, Japan aff010;  Computational Bioscience Research Center, King Abdullah University of Science and Technology, Thuwal, Kingdom of Saudi Arabia aff011;  Institute of Molecular Medicine, National Cheng Kung University, Tainan, Taiwan aff012;  Department of Ecology and Evolution, University of Chicago, IL, United States of America aff013
Published in the journal: Behavioral and brain- transcriptomic synchronization between the two opponents of a fighting pair of the fish Betta splendens. PLoS Genet 16(6): e32767. doi:10.1371/journal.pgen.1008831
Category: Research Article
doi: https://doi.org/10.1371/journal.pgen.1008831

Summary

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.

Keywords:

Gene expression – Learning and memory – Animal behavior – Transcriptome analysis – Zebrafish – Fish physiology – RNA sequencing – MAPK signaling cascades

Introduction

Animal contests result in an unequal division of resources and thus are a principal driver of natural selection [1]. Conspecifics animals fight for resources such as food and mating opportunities, but they usually stop fighting after assessing their relative fighting abilities to avoid serious injuries [24]. The contest outcomes have been commonly explained by the mutual assessment model [5], which, for example, was demonstrated in the case of zebrafish [6] and in the fish Betta splendens [7]. Yet, how such mutual assessments are carried out at the molecular level during fighting remains unknown.

Aggression is ubiquitous in animals because it may increase an individual’s chance to survive and to transmit its genes to the next generation [8]. Various types of aggression in animals have been described [9]. The fighting fish B. splendens is well known for its long fighting duration. It belongs to the suborder Anabantoidei of the perciform ray-finned freshwater fish, which are distinguished by their possession of a lung-like labyrinth organ that helps them to carry out surface-breathing [1012]. This fish has been commonly used to study the biological mechanism of aggression for various reasons. B. splendens males are extremely aggressive and have stereotypical social displays [13, 14] (S1A Fig). In their natural habitats, males fiercely protect their territories, where they build a bubble nest to hold fertilized eggs [15]. In the laboratory, a male aggressively attacks any intruder or its own mirror image to maintain its territory [14, 16]. Previous studies on B. splendens and other fish focused on aggressive displays between males but usually described only one particular behavior (e.g., displaying, biting, or striking) [17, 18] or compared global differences in aggressiveness between the loser and the winner [1923].

More recently, physiological and neurological studies have been conducted to examine aggressive behaviors using B. splendens in combination with pharmacological treatments, and some important signaling molecules linked to aggression have been discovered, such as serotonin (5-HT) [24], dopamine [25], and GABA [26]. Additionally, a suite of genes linked to aggressiveness has been identified in different fish such as cichlid fish [27], zebrafish [28], and rainbow trout [29]. However, these analyses focused on differentially expressed genes (DEGs) in the brains of the eventual winner and loser, revealing only the differences in the brain-transcriptomic state between them without comparing the brain-transcriptomic responses of the two opponents in a fighting pair during a fight. As a result, no comprehensive molecular analysis of the brain-transcriptomic state of male-male competition during a fight in B. splendens has been attempted at the genomic level. Two major obstacles to such studies were the large number of genes known to influence aggression in animals [30, 31] and the lack of a reference genome for this fish, which was released only recently [32].

Behavioral and physiological synchronization have long been studied by scientists given its importance in reproductive/mating behavior, cooperative behavior, aggressive behavior, etc. [33, 34]. Crucial signaling molecules such as 5-HT and arginine-vasotocin (AVT) modulate the levels of such cooperation at the physiological level in fish [3537]. Furthermore, recent studies have provided physiological evidence for such synchronization at the neuronal level between two interacting bats or mice [38, 39]. Given that clear behavioral evidence over the course of several decades as well as the recent physiological findings described above, we hypothesized that synchronization at the gene expression level between two interacting partners, especially in brain genes associated with synaptic function, ion transport, and aggression, might also occur. We specifically ask the following questions 1) what is the relationship between brain-transcriptomic activity in two individuals that are engaged in social interactions? 2) how does this relationship change across different timescales to facilitate social interactions?

We evaluated this possibility by examining the behavioral and brain-transcriptomic changes between two interacting opponents of B. splendens during fighting. Males of this fish undergo a stress response during a 15-min fight [14, 16]. Over a similar time period, differences in the expression of early genes as well as other socially responsive genes have been detected in various fish [4047]. Accordingly, we examined a 20-min and a 60-min fighting period to correlate the brain transcriptome responses with behavioral changes [4447]. Whereas the 20-min fight allowed us to capture immediate early genes (IEGs)-known as genomic markers for brain activity and markedly new experiences for fish, the 60-min fight permitted us to evaluate brain-transcriptomic responses for the full range of volitional behavior in response to an opponent over the course of the encounter. We found the presence of a unique temporal fighting structure by comparing the frequencies of each behavior (i.e., biting/striking, mouth-locking, and surface-breathing) of the two opponents in the intervals between mouth-lockings over the time course of a fight. Also, we discovered the synchronization of the brain transcriptomes between the opponents in a fighting pair by testing for correlation values of gene expression between paired fish and unpaired fish. Furthermore, we identified and examined the synchronized gene sets by analyzing DEGs, gene ontology and gene pathway enrichment to understand what kinds of brain activity might underlie the transcriptomic synchronization during competitive interaction.

Results

Behavioral synchronization achieved through fighting interaction

To examine the behavioral changes of the two opponents of a fighting pair, we selected male fish and allowed them to fight in pairs for a duration of 20 min (D20) or 60 min (D60) or until the conflicts were resolved (typical fighting pairs) (Fig 1A and 1B). The data obtained from 17 typical fighting pairs revealed that the fighting process in B. splendens proceeded in a stereotypic manner: surface-breathing (started at 0.94 ± 0.24 min), biting/striking (3.49 ± 0.60 min) and mouth-locking (10.60 ± 1.12 min) (S1 Table, S1 Video). During a 60-min fight, the mouth-locking events occurred an average of 12.12 ± 2.40 times, and each event lasted for 1.29 ± 0.24 min; therefore, mouth-locking events accounted for ~26.88% of the 60 min. Although the fighting duration varied dramatically among the pairs, surface-breathing, biting/striking, and mouth-locking occurred in this order among all fighting pairs, showing a temporal fighting structure.

Fig. 1. Schematic representation and behavioral synchronization.
Schematic representation and behavioral synchronization.
(A) Setup of the fighting experiment. (B) Schematic illustration of dynamic fighting behaviors between two male fish. The fighting duration and the first occurrence of each specific behavior are recorded. Abbreviations: Br, surface-breathing event; B–S, biting and striking attack; and ML, mouth-locking. (C) Behavioral differences among all fighting pairs with respect to biting/striking attack frequency. (D) Behavioral differences among all fighting pairs with respect to surface-breathing frequency. Behavioral analysis of two opponents in a pair; the data were obtained from three pairs, M5 vs. M8, E13 vs. E14, and E1 vs. E2 with respect to biting/striking attacks and surface-breathing events. Blue lines indicate the frequencies of behaviors of the ultimate winner, and dark-yellow lines indicate those of the ultimate loser in a fighting pair. Black bars indicate mouth-locking events and durations, during which biting/striking and surface-breathing behaviors are restricted, n = 6.

To understand in detail the fighting configuration, a 60-min fighting period for eight of the 17 typical fighting pairs was analyzed in 2-min windows with a 1-min overlap in terms of the frequency and duration of surface-breathing and biting/striking (see Methods). Among pairs, the patterns of surface-breathing and biting/striking were very similar between the opponents of each fighting pair (Fig 1C and 1D, S2A–S2J Fig), because these behaviors occur only between two consecutive mouth-lockings. Accordingly, the frequency, duration, and distribution of mouth-lockings during the fighting period determined the fight structure, and each fighting pair showed a pair-specific fighting pattern. More specifically, the frequencies of biting/striking and surface-breathing behaviors of the two opponents showed a high correlation for each of the eight pairs (S2 Table and S3 Table).

Brain-transcriptomic response triggered by the fighting experience

To examine the transcriptome changes during a fight, we conducted RNA-seq analysis of 25 whole-brain RNA samples, including those from five non-fighting individuals (i.e., before fighting group, B), from five pairs after a 20-min fight (the D20 group), and from five pairs after a 60-min fight (the D60 group) (Fig 2A; S4 Table). A principle component analysis (PCA) and a linear discriminant function analysis (LDA) from the 25 libraries were done on the same dataset—top 50% most variable gene transcripts (11,653 genes) to differentiate the fighting groups (see Methods). While the PCA clearly illustrated that paired fish in the D60 group were closely clustered together (S4E Fig), the LDA nicely separated the fighting groups with two significant functions: Function 1 (LD1), which separated the D60 group from the D20 and non-fighting group, explained 69.1% of the variance, and Function 2 (LD2), which separated the non-fighting group from the D60 and D20 groups, explained 30.9% of the variance (Fig 2B). A high consistency among the brain transcriptome profiles induced by each social experience (B, D20 and D60) indicates that a particular brain-transcriptomic state (co-regulated gene sets) underlies each behavioral state.

Fig. 2. Brain-transcriptomic response.
Brain-transcriptomic response.
(A) Overview of experimental design. (B) Clustering of all 25 brain samples using the linear discriminant function analysis (LDA); blue, non-fighting group (B); red, D20 group; and green, D60 group. (C) Bar graph showing the number of DEGs generated from the comparisons of B vs. D60 and B vs. D20. Significantly upregulated and downregulated DEGs are represented in red and green, respectively. (D) Volcano plot of the DEGs obtained from the B vs. D20 and B vs. D60 comparisons. Vertical lines indicate the threshold for a relative expression fold change (FC) of 2 or −2 as compared with controls (B). The horizontal line represents the threshold of a 0.05 FDR value. The red points were significantly upregulated or downregulated in B vs. D20 and B vs. D60. (E) A plot in which the D20 DEGs are projected onto the D60 DEGs and vice versa; ‘ns’ means not significant, ‘both’ means significant in both D20 and D60 groups. Blue, significant only in the D20; green, significant only in the D60.

To evaluate the brain transcriptomic differences between the D20 and D60 groups, we tested for DEGs. Comparing each social treatment (i.e., D20 or D60) with the non-fighting group (B) revealed a total of 1,082 DEGs (FDR < 0.05, |log FC| > 2), with a greater number of up-regulated genes than down-regulated genes in both fighting groups (Fig 2C). Noticeably, FDR means false discovery rate and FC refers to fold change. Whereas 859 and 518 genes showed significant increases, only 40 and 46 genes showed significant decreases in the D60 and D20 groups, respectively (Fig 2D). A plot in which the D20 DEGs are projected onto the D60 DEGs and vice versa revealed that most of the genes with down-regulated expression in the D20 group were also down-regulated in the D60 group and that those with up-regulated expression in the D20 group remained up-regulated in the D60 group. There were only three genes (unannotated) that were down-regulated in the D20 but were up-regulated in the D60, whereas just one gene (unannotated) was up-regulated in the D20 but was down-regulated in the D60 (Fig 2E).

Among 1,082 DEGs, 564 genes were specifically expressed in the D60 group (D60-specific genes), 223 genes were specifically expressed in the D20 group (D20-specific genes) and 295 genes were expressed in both the D20 and the D60 groups (common genes) (Fig 3A, S5 Table). We found that the probability of 295 overlapping DEGs between the D20 and D60 groups was greater than that expected by chance (p < 7.2e-298, hypergeometric test). In addition to the attributes of their expression across fighting stages, several of these DEGs genes have previously known to associate with neuronal activity (e.g., c-fos, fosb, egr4), neural plasticity processes related to learning and memory (e.g., bdnf, npas4, nr4a1), neurotransmitter (e.g., grin2a, grik2) and epigenetic function (rnf44) (Fig 3B and see Discussion) [6, 48]. The heatmap generated from all 1,082 DEGs illustrated three main patterns of gene expression that were distinguishable among the fighting groups (B, D20, and D60) (Fig 3C). We also note that a large set of unannotated DEGs was uncovered in this analysis (S5 Table).

Fig. 3. Characterization of DEGs.
Characterization of DEGs.
(A) Venn diagram displaying the number of common and specific DEGs identified in the B (non-fighting group) vs. D60 and B vs. D20 comparisons. (B) Representative genes among the DEGs that are related to neuronal functions; the numbers in the parentheses indicate which clusters in the heatmap these genes belong to; hyphen ‘-’ means not found. (C) Heatmap using all 1,082 DEGs illustrates four main patterns of gene expression within each of the three groups (B, D20, and D60). Similarities between the fighting groups as shown by hierarchical clustering can be seen above the heatmap. The bootstrap value at the node was obtained by the hclust function in R. Intensity of color indicates expression levels: red, high expression; blue, low expression. The numbers (1−3) on the right side refer to the three main patterns of gene expression.

Genes associated with each fighting group (D20 and D60)

A weighted gene co-expression network analysis (WGCNA) was used to find gene clusters (modules) among all 23,306 genes across all 25 samples (see Methods). Our goal was to identify modules associated with the fighting groups (D20 and D60). This analysis identified 37 co-expression gene modules, which are represented by different colors as named by the program (Fig 4A, S4D Fig). To relate these gene modules to the fighting groups (D20 or D60), associations between the eigengene of each gene module and the fighting groups were computed. The results showed that the two fighting groups were associated with different modules. Particularly, the D20 samples were characterized by an up-regulation of the MEmidnightblue gene module, which included 186 genes that are associated with transcription regulation (Fig 4B, S3A Fig and S3F Fig). The D60 group was characterized by three modules: an up-regulation of the MEgrey60 module (168 genes), which is enriched in genes related to stress response, transcription regulation, and MAPK signaling (Fig 4B, S3B Fig and S3G Fig); an up-regulation of the MEturquoise module (6,390 genes), which is enriched in genes related to phosphorylation, ion transport, protein transport, and MAPK signaling (Fig 4B, S3D Fig and S3I Fig); and an down-regulation of the MEblue module (4,256 genes), which is enriched in genes related to translation, metabolic process, DNA-repair, and MAPK signaling (Fig 4B, S3C Fig and S3H Fig).

Fig. 4. Co-expression gene modules for different fighting groups (B, D20, and D60).
Co-expression gene modules for different fighting groups (B, D20, and D60).
(A) Associations between patterns of expression in the 37 identified modules across all three groups (B, D20, and D60). The colors of the boxes are scaled with the value of the correlation coefficients, ranging from red (r = −1) to green (r = 1). (B) Eigengene values of samples separated by group (B, D20, and D60) for gene modules significantly associated with four groups (blue, midnight blue, turquoise, and grey60). Representative genes enriched in each module are shown.

We found that the D20-specific genes, D60-specific genes, and common genes in the DEG analysis appeared in different modules in this WGCNA analysis. The D20-specific genes were present in the MEbrown (104 genes), MEred (86 genes), MEgreen (5 genes), MEdarkgreen (4 genes), and other gene modules (S6 Table). The D60-specifc genes were present in the MEturquoise (444 genes), MEblue (29 genes), MEgrey60 (18 genes), MElightgreen (17 genes), and other gene modules (S6 Table). The common genes were present in the MEbrown (170 genes), MEturquoise (79 genes), MEred (20 genes), and other modules (S6 Table). The DEG and WGCNA analyses used different approaches but showed a high agreement in the characterization of the transcriptome changes associated with each fighting group across the fight.

Brain-transcriptomic synchronization occurs in a pair-specific manner

We discovered the presence of brain-transcriptomic synchronization by first testing for the correlation coefficients (r values) of expressed gene transcripts between the paired fish and unpaired fish over the 23,306 gene transcripts for the D20 (D20-paired fish vs. D20-unpaired fish) and D60 (D60-paired fish vs. D60-unpaired fish) fighting groups (S4A Fig and S4B Fig). The results revealed that the r values for paired fish were significantly higher than those for unpaired fish in both groups (permutation test, p < 0.05), especially for the D60 group (permutation test, p = 0.001) (Fig 5A). Furthermore, the clustering analysis in the heatmap using the top 10% most variable genes (2,330 genes) showed that all five pairs of the D60 group were clustered together, in sharp contrast to the single pair that clustered together in the D20 group (Fig 5B). It is interesting to note that similar results were obtained when individual profiles of all the 25 samples using all gene set (23,306 genes) were clustered (S5C Fig).

Fig. 5. The structure of brain-transcriptomic synchronization.
The structure of brain-transcriptomic synchronization.
(A) Boxplot of correlation coefficient values (r) for comparisons of the TMMs for 23,306 gene transcripts between paired fish and unpaired fish from the same treatments (D20-paired vs. D20-unpaired and D60-paired vs. D60-unpaired, permutation test). *p < 0.05; ***p < 0.001. (B) Heatmap using the 25 cDNA libraries with each of the 2,330 gene contigs (top 10% most variable genes). Intensity of color indicates expression level (dark orange, high expression; sky blue, low expression). Similarities between individuals within fighting pairs and between fighting pairs as shown by hierarchical clustering can be seen above the heatmap. Bootstrap values at the nodes were obtained by hclust function in R. Bootstrap values in red indicate the two opponents in a pair were clustered together. (C) GoM analysis using TMM values of all 23,306 genes generated from the 25 libraries with a set of seven clusters represented by different colors.

Next, we subsequently examined the fine structure of brain-transcriptomic synchronization in each fighting pair by applying the grade of membership (GoM) model to the 25 brain samples, each of which included the 23,306 gene transcripts. This GoM model allowed us to highlight similarities among the samples by determining their similar membership (cluster) proportions (see Methods). Interestingly, the composition of clusters for individuals in each pair from the D60 group was very similar, and moreover the pattern from any given pair of the D60 group was unique relative to that of the other pairs. This observation suggested that brain-transcriptomic synchronization is achieved according to the specific interactions between the two opponents of a pair, such that a particular set of genes was upregulated or downregulated and synchronized in a pair-specific manner (Fig 5C). The top five driving genes (i.e., the genes that were most distinctively differentially expressed) in each cluster were extracted from the total gene transcripts and were found to be unique for each cluster (S7 Table). Altogether, these findings suggest that the brain-transcriptomic state of the two opponents in a pair were mutually influenced during fighting and became synchronized in a pair-specific manner during the fighting process. This synchronization was weak in the D20 group but became strong in the D60 group.

Characterization of synchronized genes

To identify synchronized genes, we conducted a two-step analysis (see Methods). First, we obtained 868 and 2,409 DEGs from the comparison of the D20 and D60 groups relative to the non-fighting group, respectively (FDR ≤ 0.05 and log FC > 0) (Fig 6A, S8 Table). Here we focused on the up-regulated gene sets because the number of these genes was much higher than that of the down-regulated genes among the DEGs in both groups (Fig 2C), and we used a less stringent criteria than was used in the analysis of Fig 2C (log FC > 0 instead of |log FC| >2) so that more DEGs will be included for the downstream analysis. Next, from these DEGs, we extracted 172 and 1,522 genes that were synchronized in at least one pair from the D20 and the D60 group, respectively, according to the criteria defined in the Methods (Fig 6A, S9 Table). Comparisons between the two synchronized gene sets showed that the number of synchronized genes within one pair, two pairs, three pairs, four pairs, and five pairs was much higher in the D60 group than the D20 group (Fig 6B). The composition of these synchronized gene sets confirmed that brain-transcriptomic synchronization is a major process because more than half (63.2%) of the total DEGs were involved in this process. Moreover, almost half of the synchronized genes were synchronized in five fighting pairs (Fig 6C). It is interesting to note, however, that the synchrony levels were different among each particular fighting pair as shown for neuronal activity−associated genes such as bdnf, calm2a, and several other genes (Fig 6D, S4C Fig). The overlapping synchronized genes among the five pairs in the Venn diagrams show the number of genes involved in brain-transcriptomic synchronization across all five fighting pairs as well as those that were shared between each particular fighting pair in both the D20 (S5A Fig) and D60 (S5B Fig) fighting groups.

Fig. 6. Synchronization associated with learning and memory genes.
Synchronization associated with learning and memory genes.
(A) Workflow illustrating the GO and KEGG pathway analyses. (B) Bar graph showing the number of synchronized genes in one, two, three, four, and five pairs in the D60 and D20 groups. (C) Pie charts showing proportions of synchronized genes. (Left), the percentage of synchronized vs. non-synchronized genes; (right), synchronized genes in one, two, three, four, and five pairs; the data here represent pooled data for both groups (D) Examples of genes showing pair-specific individualization of brain-transcriptomic synchronization. (Left), representation for identification a baseline (Bi) of a particular gene i in pair (s) j. The gene i demonstrated here is camk1da from the D60 group. Each dot represents a pairwise (the x-axis shows the total of 45 pairwise from the D60 group); blue represents the D60 paired fish; orange represents the D60 unpaired fish; the green line refers to the 25th percentile rank. (Center), bdnf represents IEGs and (right), calm2a represents memory related genes. (E), Significantly enriched biological processes among the synchronized genes in the D60. (F) Simplified calcium signaling pathway map generated from the D60 synchronized genes (see S6A Fig and Discussions for detailed information). The up-regulated synchronized genes involved in memory consolidation are shown in red. GPCR, G protein−coupled receptor; CaV1, L-type calcium channel.

To gain insight into the biological functions of synchronized genes involved in brain-transcriptomic synchronization, we first examined the overrepresentation of the Gene Ontology (GO) terms related to biological process (BP) within the 1,282 synchronized genes (all of which were well annotated) in the D60 using DAVID (p < 0.05) (Fig 6A). Accordingly, 49 overrepresented BP terms were enriched (S10 Table); however, as the GO terms obtained from DAVID seemed to be redundant, we next used the REVIGO tool to filter these redundant BP terms. Subsequently, 31 overrepresented BP terms remained (S10 Table). Interestingly, this synchronized gene set was mainly enriched in terms related to neurotransmitter secretion, ion transport, and synapse function (Fig 6E, S10 Table). We also looked for overrepresentations among molecular function (MF) and cellular component (CC) terms. Enriched MF terms were associated with transcription factor activity, receptor activity, and ion channel activity (S11 Table), whereas enriched CC terms were associated with membrane, synapse, and AMPA glutamate receptor complex (S12 Table). Furthermore, overrepresentation analysis of KEGG pathways (p < 0.05) within 1,282 synchronized genes revealed 19 pathways that were significantly enriched and are known to have causal effects on social behavior such as calcium signaling [49], MAPK signaling [50], GnRH [51], and others (S13 Table).

Brain-transcriptomic synchronization associated with long-term memory genes

In mice and bats, social interactions between two interacting partners induce inter-brain correlation of neural activity in the prefrontal cortex of their brains [38, 39], this activity has been implicated in supporting working memory in human [52]. Learning and memory are crucial processes in the winner or loser effect phenomenon, in which fish increase or decrease their likelihood of winning or losing in a future fight after experiencing previous wins or losses, respectively [53]. This motivated us to examine whether genes related to long-term memory were involved in the brain-transcriptomic synchronization in the D60 group.

We first extracted 378 genes from the 1,522 synchronized genes of the D60 group that were linked to learning or memory (see Methods). Next, these 378 genes were used for a KEGG pathway enrichment analysis using DAVID (Fig 6A). As a result, 19 pathways were significantly enriched such as insulin signaling pathway, calcium signaling pathway, adrenergic signaling pathway, etc. (p < 0.05, S13 Table). Among the top significantly enriched pathways shown in S13 Table, the calcium signaling pathway, which was ranked second when pathways were sorted by p-values, is well known for its involvement in learning and long-term memory formation [49]. In support of this finding, several genes in this pathway were up-regulated in the route leading to learning and memory (Fig 6F, S6A Fig).

Next, we used the calcium pathway as a proxy for learning and memory examination for other gene sets including 104 synchronized genes in the D20, 567 non-synchronized genes in the D20, and 791 non-synchronized genes in the D60 to see if they were involved in pathways related to learning and memory (Fig 6A). By comparing the gene expression of those involved in the calcium pathways generated from the KEGG analysis, we found that the D20 synchronized gene set was not enriched for the calcium pathway, and none of the genes involved in calcium pathways generated from the non-synchronized gene sets were included in the pathway of learning and memory (S6B Fig and S6C Fig). Collectively, these results suggest that learning and memory may play a key role in brain transcriptomic synchronization after the fish experience a long period of fighting (see Discussion).

Discussion

Behavioral synchronization in vertebrates

In the case of competition among conspecific male animals for resources and mates, they typically stop fighting after assessing their relative fighting abilities to avoid serious injuries. How animals assess another’s ability remains poorly understood. Here, we demonstrated behavioral synchronization in each B. splendens pair during fighting. Although the behavior sequences (displaying, biting, and striking) in B. splendens are in accordance with other fish such as zebrafish (Danio rerio) [53] and cichlid fish (Nannacara anomala) [54], characteristic pair-specific behavioral synchronization as defined by the pattern of mouth-locking behavior is observed only in B. splendens, underscoring the unique behavior of this fish. During fighting, it seems that each B. splendens male in a fighting pair attempts to increase its chance of winning by imposing its aggressive behaviors upon its opponent, leading to tightly synchronized behaviors. A few previous studies have found synchronization in other animal behaviors such as foraging (seeking food) and cooperative hunting [5558]. For example, during foraging by a pair of sticklebacks or Atlantic salmon fish, if there is an advantage to foraging together, the behavior of both individuals becomes synchronized; however, when one of the two fish receive additional food, they tend to behave independently.

Our behavior analysis showed the presence of a specific behavioral synchronization pattern for each fighting pair. It is likely that each individual of B. splendens has its own behavioral characteristics, referred to as “personality” (e.g. shy, bold, proactive, or reactive) [59]. Accordingly, the characteristics of the synchronization pattern for each pair might be influenced by the personality of the individual fish involved (e.g., how shy or bold they are), as each B. splendens individual needs to modify its behavior to match the ability of its competitor. This explanation is consistent with previous studies in other fish such as guppies and sticklebacks and in pair-bonded birds such as the great tits (Parus major), in which the frequency and distribution of their interactions depend on their body length, color, and proactive or passive exploration personality [60, 61].

A model for behavioral and brain-transcriptomic synchronization

To explain how behavioral synchronization and brain-transcriptomic synchronization occur, we propose the model presented in Fig 7, in that the mutual assessment could be involved. The repeated attempts of each opponent to match the actions of the other opponent lead to behavioral and then transcriptomic synchronization. At the behavioral level, the two fish initiate and escalate the fight (Fig 7A). Their interactions provide signals for the mutual assessment process to assess the opponent’s fighting ability. Moreover, to match the attacks of its opponent, each fish adjusts its fighting behavior to achieve synchronization.

Fig. 7. A model for behavioral and brain-transcriptomic synchronization.
A model for behavioral and brain-transcriptomic synchronization.
(A) Two fish, A and B, initiate aggressive behaviors (opercular display, fin spread, bite/strike). The repeated attempts of each fighting opponent to match the actions of the other opponent lead to behavioral synchronization. (B) Continued physical interaction provides repeated signals for the brains of each fish to adjust the expression of genes in the neurocircuits as mediated by mutual assessment; these genes include immediate-early genes (IEGs), synchronized genes, and learning/memory-associated genes, gradually leading to brain-transcriptomic synchronization. From 0 min to 20 min, different colors are used to refer to transcriptomic plasticity and to partially reflect the gene expression patterns of the DEGs in groups 1 & 2 in Fig 3. From 20 min to 60 min, a single color was used to represent transcriptomic synchronization and to partially reflect the gene expression patterns of the DEGs in groups 3 in Fig 3. (C) Learned association through mutual assessment leads to changes in the physiology of the neurocircuits as mediated by second messengers and intracellular signaling pathways to convey signals for brain-transcriptomic synchronization.

At the physiological level (Fig 7B), the signals of mutual assessment in each fighting individual trigger a transcriptomic response in neurons leading to changes in the physiology of neurocircuits. The first transcriptional wave after neuronal activation induces rapid transcription of immediate-early genes (IEGs) such as c-fos [43]; brain-derived neurotrophic factor, bdnf [62]; and nuclear receptor transcription factor subfamily nr4a(1,3) [63]. Then, a downstream genomic response cascade takes place with the involvement of synchronized genes (S9 Table) and learning- and/or memory-associated genes such as calm2a, camk2b1, ppp3r1b, and others [64], eventually leading to the pair-specific individualization of brain-transcriptomic synchronization (PIBS). These processes are mediated by second messengers, intracellular pathways such as the MAPK pathway [50], ligand-receptor pathways, and/or calcium signaling pathways [64]. Collectively, PIBS presents a mechanism for coordinating and sustaining interactions between opponents.

Although our data are correlative and many details of the proposed model at the molecular level remain unclear, the model may provide a useful basis for future study. It will be interesting to see if both behavioral and brain-transcriptomic synchronization also occur during other interactive behaviors.

Insights into brain-transcriptomic synchronization

We demonstrated the presence of brain-transcriptomic synchronization, which was observable to a lesser extent after a 20-min fight and was strengthened after a 60-min fight in accordance with behavioral synchronization between fighting pairs of B. splendens. We characterized it as follows.

First, the DEG analysis illustrated three main patterns of gene expression that reflected changes in the transcriptomic landscapes between the social groups (B, D20, and D60; Fig 3C), with the involvement of IEGs in both fighting groups. As mentioned above, some IEGs exhibited strong upregulation in the D60 group such as cfos, bdnf, nr4a3, and neurotransmitter transporter glycine member 9 (slc6a9) [65] (Fig 3B), which are involved in long-term memory formation. Thus, a fighting B. splendens seems to activate molecular pathways involved in memory formation during fighting. This observation is further supported by the detected enrichment of the MAPK signaling pathway in the D60-specific gene sets, which is known to be associated with cognitive processes [50]. Several dual specificity phosphatase (dusp) subunits (dusp2, dusp4, dusp3a, and dusp7), modulators of the MAPK signaling pathway [66], were upregulated and synchronized. In particular, the involvement of the calcium signaling pathway appears to be crucial. Different routes of Ca2+ entry lead to the activation of different downstream pathways (Fig 6F; S6C Fig and S6D Fig). The activation of genes encoding the NMDAR subunits (grinaa, grin2aa, and grin2b), AMPAR subunit (grm8a), and CaV subunit (cav3, T-type calcium channels and cav1, L-type calcium channels) leads to increases in nuclear Ca2+ [67], which in turn binds the calmodulin (CALM) subunit (calm2a) to activate the α-calcium calmodulin-dependent protein kinase II (CAMK2) subunits (camk1da, camk2, and camk2b1) and calcineurin (CaN) subunit (ppp3r1b), eventually triggering long-term memory formation [68]. Noticeably, each fighting pair activates a unique route within the Ca2+ signaling pathway (S14 Table).

Interestingly, the number of up-regulated genes was much higher than that of down-regulated genes among DEGs in both fighting groups (Fig 2C). Our gene enrichment analysis showed that the DEGs in the D20 and D60 were associated with ion transport, synaptic function, and long-term memory. The gene with the highest expression was thoc2, a subunit of three prime repair exonuclease (TREX), which is required for efficient export of polyadenylated RNA and spliced mRNA [69]. This asymmetry in gene expression variation between the B and D20/D60 groups was quite striking and is consistent with previous studies in sticklebacks [70], mice [71], and honeybees [72] to some extent. These studies demonstrated that up-regulated DEGs were involved in insulin signaling, oxidative phosphorylation, mitochondrial function, and regulation of metabolic processes, showing deeply conserved metabolic signal after challenge and representing a shift in the balance of metabolism [71]. Together, highly expressed genes from the D20 and D60 groups suggest the possibility of a shift in metabolism that prefigures later neuronal plasticity in the brain. It is also possible that the asymmetry in DEGs between baseline and fighting groups could be a molecular manifestation of increased inhibitory control of gene expression on aggression so as to avoid an over-escalation [73].

Second, the heatmap and the GoM analyses of the 25 brain samples collected at three different fighting time points together with DEG analysis revealed the presence of brain-transcriptomic synchronization, with more than a half (63.2%) of all DEGs (B vs. D60) involved in this process (Figs 5 and 6). The GoM analysis demonstrated that the composition of multiple clusters from individuals in a pair, especially D60 pairs, was very similar, illustrating nicely the presence of brain-transcriptomic synchronization in a pair-specific manner (PIBS). The GoM model approach has been widely used in population genetics to cluster admixed individuals who have ancestry from multiple populations. Recently, it has been used for assessing global structure in RNA sequencing data in mice and humans [71, 74]. Future studies using this approach along with existing and widely used visualization methods such as LDA and hierarchical clustering (heatmap) will provide a richer summary of the information contained within RNA-seq data [75]. In addition, it will be interesting to find a way to quantify the strength of association between the GoM results and behavioral synchronization of any given pair. Currently, however, as our behavioral analysis here did not include all sets of behaviors requiring synchrony, such as tail-beating, body-orientation, and body-coloration, such quantification will require future studies.

Third, the WGCNA showed that each of the two groups, D20 and D60, was associated with distinct gene modules: MEgrey60, MEblue, and MEturquoise for the D60 group and MEmidnightblue for the D20 group (Fig 4). Although the DEG analysis and WGCNA represent different approaches, they gave consistent results with respect to the characterization of transcriptome changes associated with each fighting group. Of the 564 D60-specific genes identified by the DEG analysis, 491 (87.06%) were present in the three gene modules associated with D60 based on the WGCNA (S6 Table). Our results are consistent with a study that examined DEGs between the winner, loser, and mirror phenotypes (i.e., when a fish fights against its own image in a mirror) in zebrafish using both DEG and WGCNA approaches [6]. Although this study identified fewer gene modules, it showed extensive agreement for the DEG numbers between the two methods [6]. Other recent studies using both DEGs and WGCNA methods have demonstrated fruitful results [71, 72]. Thus, this integrative approach is useful for examining the relationship between gene expression and targeted phenotypes, and one can corroborate the other.

Transcriptome synchronization using whole brains vs. brain regions

We used the whole brain in this study given that social influences and behavioral changes induce modifications in brain gene expression [76], and that external information can trigger brief variation from the baseline of brain gene expression [77]. Although it is uncertain how these changes are conveyed to other organs, clearly the brain plays a key role in where social experience is integrated with physiological changes to carve the phenotypic plasticity [78]. However, one disadvantage of the analysis using whole brain is that brain tissue consists of a complex mixture of neuronal cell types, each of which has different biological signals, and that some of meaningful signals may be undetected because of their relative low abundance [71,78]. It was found that none of the candidate genes for aggressive behavior in the previous study [28], such as tryptophan hydroxylase (tph1 and tph2), were identified in studies using whole brain such as those in zebrafish [6] and cichlid [27], as well as in the present study. It is possible that these aggression-associated genes were up- or down-regulated in a specific part of the brain but could not be detected by using whole brain, possibly because any changes were diluted. It has been reported that dominant and subordinate zebrafish showed different number of DEGs in different brain regions after social interaction, such that some genes are overexpressed only in the hypothalamus and hindbrain of the dominant, but other genes are overexpressed in the telencephalon and optic tectum of both dominant and subordinate fish [28]. This suggests that difficulty of analysis for the expression of genes with low abundance in neuronal cells can be overcome by obtaining RNA from smaller area of the brain using microdissection or single cell sequencing method.

Brain-transcriptomic synchronization was observed in the analysis of the whole brain, suggesting that several parts of the brain could be involved in the synchronization events observed in this study. In recent years, an approach to map a specialized part of the brain by focusing on regions such as the hippocampus, which is involved in learning and memory [78, 79], has been developed. A comparison of gene expression profiles among brain regions (e.g., telencephalon, diencephalon, cerebellum) in fish responding to intruders [70] has been attempted, and analyses of changes in the conserved ‘social behavior network’—a collection of brain nuclei known to regulate social behaviors across vertebrates—based on detailed micro-dissected brain regions in cichlids [80], zebrafish [81], and swordtail fish [82] have been carried out. The limbic region, midbrain, and the cortical region are involved in a variety of motivations in mice [38, 39, 83]. Based on these studies, we propose that the telencephalon region could be a good candidate for further brain-transcriptomic synchronization analysis because it contains structures homologous to the hippocampus and amygdala in mammals and is associated with aggression in other vertebrates [84]. In particular, loci in this region including the dorsal medial (associated with the fear response), dorsal lateral (associated with spatial learning), lateral spectrum (associated with anxiety), pre-optic area (associated with reproduction), and cerebellum (associated with attention) are of interest [85]. A better understanding of the involvement of these specific brain loci in brain-transcriptomic synchronization is a promising area for future research.

Brain-transcriptomic synchronization reflects cooperative fighting aspect

In the present study, we propose that brain-transcriptomic synchronization might underlie the behavior synchronization between two interacting partners in vertebrates. Furthermore, we speculate that brain-transcriptomic synchronization may occur with collaborative behaviors because of the universality of behavioral synchronization under both competitive and collaborative circumstances, as discussed above.

As the paired fish had very similar behavioral repertoires overall, any transcriptomic commonalities could reflect transcriptomic responses to fighting effort or allocation of that effort rather than being related to behavioral synchrony per se. One possibility for the role of PIBS is that the brain-transcriptomic response in one individual might encode behaviors of its interacting opponent. Previous studies have identified neurons that are selectively activated during a partner’s action in songbirds, monkeys, and humans [8690]. Very recent studies reported that there is inter-brain correlation of neuron activity between two individuals of mice [38] and bats [39] during social interactions. Strikingly, in the case of the mice, correlated brain activity depends on both the encoding of one’s own behavior at the cellular level and the behavior of the interacting partner. Moreover, deeper similarities in neural responses were found among friends who view the same movie [91], or between infants and adults while they communicate and play with one another [92]. Also, highly sharing of emotional state within a judgment task was revealed in ravens (Corvus corax) [93]. Noticeably, it has been reported that different neuronal activity patterns induce different sets of activity-regulated genes [94]. Together, the existence of PIBS suggests that correlated brain gene expression may support synchronous neural activity in brains between two interacting partners, reflect cooperative aspects at the molecular level, and indicate the presence of basic forms of empathy in animals [86].

Brain-transcriptomic synchronization associated with fighting motivation

The neurotransmitter serotonin is involved in regulating various motivations in animals [95]. For example, serotonin plays an important role in prolonging fight duration in lobsters and crayfish [96]. Similarly, expression differences in the serotonin receptor subunits Htr1b and 5-Ht2c account for increased exercise motivation in mice, suggesting that higher motivation for exercise is modulated, at least in part, by the serotonin signaling mediated by the expression of Htr1b and its surrounding chromatin organization [97, 98]. Interestingly, we did find that two additional serotonin receptor subunits, htr7a and htr7c, which belong to the neuroactive ligand-receptor pathway were highly expressed and synchronized in four fighting pairs in our data (S14 Table), suggesting the possible involvement of these genes in the persistence of motivation in this system.

In summary, this study demonstrates that fighting behaviors and brain gene expression of the fish B. splendens within a fighting pair are highly synchronized. While the behavioral synchronization is observable at an early stage of fighting, the brain-transcriptomic synchronization is discernible after a 20-min fight and is strengthened after a 60-min fight, suggesting the importance of repeated signals for the brain to adjust the expression profiles of certain genes. We demonstrate that brain-transcriptomic synchronization is observed between individuals of a fighting pair in a pair-specific manner after sharing long-term interactions. This is a phenomenon that has been widely known at the behavioral level but has rarely been observed at the molecular level and could occur widely among animals. Altogether, this study provides a framework for an understanding of synchronization at the behavioral and molecular levels that occurs during a contest or collaboration in higher animals through social interactions.

Methods

Ethics statement

The animal experimentation procedures used in this study were approved by the Institutional Animal Care and Use Committee (IACUC) (Approval No.106171) of the National Cheng Kung University, Tainan, Taiwan.

Animals and maintenance

Male fighting fish B. splendens used in this experiment were obtained from a local fish shop in Thailand. During delivery, the fish were individually kept in 150-ml glass flasks. In the laboratory, they were transferred to 600-ml glass flasks and housed at 26 ± 2°C with a 12-h dark and 12-h light photoperiod for 2 days for adaptation. They were then transferred to individual tanks in the circulation aquarium system in the laboratory. Fish were fed twice a day with commercial food flakes. The average fish size was 5.2 ± 1.1 cm (standard length).

Behavioral assays

Fish were maintained in the aquarium system tanks for 1 week before a fighting experiment. To test for the interaction, a group of five non-fighting individuals (B1, B2, B3, B4, and B5) was used as the control group. Another 20 adult males, matched on the basis of standard length, were used for RNA sequencing. Three experimental sets were conducted in which 10 fish (five pairs: D20_11 vs. D20_12, D20_21 vs. D20_22, D20_31 vs. D20_32, D20_41 vs. D20_42, and D20_51 vs. D20_52) were paired for a 20-min fight, 10 fish (five pairs: D60_11 vs. D60_12, D60_21 vs. D60_22, D60_31 vs. D60_32, D60_41 vs. D60_42, and D60_51 vs. D60_52) were paired for a 60-min fight, and 34 fish were paired until the winner chased the loser; this last group was referred to as typical fighting pairs (S1 Table). Fish were individually recognizable based on their coloring, especially the color of their caudal, anal, and dorsal fins. Each pair of fish was introduced into a 1.7-L PVC tank (18 × 12.5 × 7.5 cm). An opaque PVC partition was temporally used to divide the tank into two parts when the two fish were first introduced into the tank to help them become stabilized. After 1 min, the opaque divider was removed. In the first experimental set, the fish were allowed to fight for 20 min; in the second set, the fighting duration was set for 60 min; and in the third set (referred to as the typical fighting pairs), the fish were allowed to fight until the winner chased the loser. Behavioral interactions were videotaped digitally, and the recordings were subsequently used for detailed behavioral analysis. The fish from each experimental set were immediately killed after the fighting interactions with a lethal dose of MS-222 (Syndel USA; 1,000–1,500 mg/L, anesthetic time was ~1 min at 600 ppm concentration), then stored at -80°C for subsequent analysis. The non-fighting group was killed once the fish were taken out of the tanks, and they were removed from the tanks after all fighting group fish were killed (S1C Fig).

The conflicts in all of the D20 and D60 fighting pairs in this study were not resolved. Data obtained from detailed scrutiny of behavioral patterns exhibited during interactions from the 17 typical fighting pairs revealed that a switching point at which the resolution of the fight took place was observed after a 60-min fight. During the 60-min fights, the fighting strength between two opponents within a fighting pair showed almost the same tendency, and the fighting sequences were consistent among all of the fighting pairs. Thus, it is likely that all pairs characterized as “D20” and “D60” were actually in the same stage with respect to fighting sequences.

All 25 brain samples, consisting of five non-fighting samples (B), 10 samples from the 20-min fights (D20), and 10 samples from the 60-min fights (D60), were used to generate RNA-seq data for transcriptomic analysis. The 17 typical fighting pairs were used only for behavioral analysis. With these typical pairs, we referred to the fish by a letter and number to indicate different sets of fighting experiments (e.g., M5 vs. M8 indicates “in set M, fish number 5 fights against fish number 8”) and to clarify between the winner and the loser.

Video analysis

Recorded videos (Nikon Cool Pix E5400) were analyzed using the video editing software Windows Movie Maker (Microsoft). Activities were recorded with respect to both the frequency and timing of biting/striking, surface-breathing, and mouth-locking. Biting/striking was defined as occurring when one fish made open-mouth contact with the other fish. Surface-breathing was defined as occurring when a fish ceased displaying and began gulping air at the surface, and mouth-locking was defined as occurring when one fish took hold of the other fish’s upper or lower lip.

Fighting behavioral structure

To examine the fighting behavioral structure, biting/striking attacks, surface-breathing events, and mouth-locking events were recorded in terms of their frequency and duration. An initial 60-min fight duration was analyzed in which the biting/striking, surface-breathing, and mouth-locking behaviors were broken down into 2-min windows with 1-min overlap for analysis. To consider the entire 60-min fighting process, the time series was broken up into short, overlapping segments (or time windows). Calculation for frequency and duration values are then determined for each of these time windows (1 min in this case). The resultant time sequence of all frequency and duration values from each time window can indicate behavioral changes over time. The length of the time window was determined based on the type of movement considered, along with the observed time scale for behavioral changes [99]. Pearson correlation coefficients for bite/strike and surface-breathing (S2 Table and S3 Table) were separately calculated between paired fish using the cor function in R/Bioconductor. The figure panels showing the combination of all behaviors (Fig 1C and 1D, S2 Fig) were generated using Excel.

Tissue processing and RNA extraction

Brains were rapidly dissected from the three groups of frozen fish (B, D20, and D60) in PBS buffer and were individually collected in 1.5-ml tubes containing 1,000 μl Trizol for further processing. Total RNA from each whole brain (0.01 g/brain) was extracted using Trizol and then purified with Quick-RNA MiniPrep (Zymo Research, USA). The concentration and purity of the RNA was determined by Qubit (Eugene, Oregon, USA) and a BioAnalyzer 2100 RNA Nano kit (Agilent, USA), respectively. The RNA quality integrity (RIN) values of the samples ranged from 6.3 to 8.8. Total extracted RNA was kept at –80°C until processing.

Stranded RNA sequencing

RNA-seq libraries were constructed using the TruSeq Stranded mRNA Library Prep Kit (Illumina, USA) with proper quality control, and the molar concentrations were normalized using a KAPA Library Quantification Kit (Kapa Biosystems, USA). The 25 RNA libraries were sequenced in two separate groups. Five RNA libraries (B1, B2, B3, D60_11, and D60_12) were sequenced on the Illumina HiSeq 2500 system at Yourgene Bioscience Co., Ltd. (Taipei, Taiwan) using paired-end sequencing. The other 20 samples (B4, B5, D20_11, D20_12, D20_21, D20_22, D20_31, D20_32, D20_41, D20_42, D20_51, D20_52, D20_61, D20_62, D60_21, D60_22, D60_41, D60_42, D60_51, and D60_52) were sequenced on the Illumina HiSeq 2500 system at the NGS High Throughput Genomics Core (Biodiversity Research Center, Academia Sinica, Taiwan) using single-end sequencing. The read lengths and total sequencing reads are shown in S4 Table.

Given that the samples had undergone two different sequencing methods, we examined whether the paired-end vs. single-end difference led to any biases in the data using a multidimensional scaling (MDS) plot. The MDS plot, which was color-coded based on the sequencing method, revealed that the two methods resulted in slight or no biases as all samples were clustered together (S7 Fig). We further evaluated the extent to which the two methods yielded equivalent results. We used only the forward-reads for the pair-end sequencing samples and analyzed all data in single-end mode throughout the pipeline to minimize any bias. The results obtained were comparable in terms of mapping rates and the total number of genes across all the samples (S4 Table) between the two methods, which is consistent with a previous study [100] and indicated sufficient quality for further downstream analyses.

Quality processing and read mapping

Adaptors were trimmed using the Cutadapt tool [101] with the following command line, -b GATCGGAAGAGCACACGTCTGAACTCCAGTCAC -b AGATCGGAAGAGCGTCGTGTAGGGAAAGAGTGT -o $1.output file name $1.inputfile name. The low-quality bases and reads were removed by using the fastq_quality_trimmer tool version 0.0.13 (http://hannonlab.cshl.edu/fastx_toolkit/) with the following parameters: fastq_quality_trimmer -t 20 -l 30 -Q 33 -i $1.input file name | fastq_quality_filter -q 20 -p 80 -Q 33 -o $l.output file name.

The B. splendens genome and its gene IDs were downloaded from http://gigadb.org [32]. The processed sequencing reads were mapped to the genome using Tophat version 2.1.1 [102], and its embedded aligner Bowtie2 version 2.1.0 with the default parameters [103]. The unique mapping reads (reads that matched the reference genome at only one position) were extracted using samtools with the following command line view, -q 4 $1.input file name >$1.output file name. The exon-mapped reads were counted by featureCounts [104]. The normalized expression levels of genes, represented by the trimmed mean of M-values (TMM) [105], were generated by the “edgeR” package (version 3.26.8) in R [106]. In total, 23,306 genes from each brain sample were generated. All raw data are available at Bioinformation and DDBJ Center (https://www.ddbj.nig.ac.jp/index-e.html), accession code number DRA009599.

Identification of DEGs

The DEGs were obtained from comparisons of gene expression between the non-fighting group (B) with the D20, and B with the D60. We included genes with at least one count per million in at least one sample. The p-values from all contrasts were adjusted via empirical FDR at once. Cut-off values of FDR ≤ 0.05 and |log FC| > 2 were used to select the DEGs (S6 Table; source code, edgeR).

Brain-transcriptomic synchronization analysis

To evaluate the transcriptomic synchronization, we first calculated Pearson correlation coefficients of log2-transformed TMM values (r) separately between two random fish from each group of the D20 and D60 using the pairs function in R. Combinations of all possibilities between two random fish among the 10 fish resulted in a total of 45 pairwise (five paired fish and 40 unpaired fish) for each fighting group (45 = 2-permutations of 10; S4 Fig). Next, to assess how similar the transcriptomes of paired fish were and to what extent the paired fish were synchronized, we compared the r values from paired fish to the r values from unpaired fish from the same treatment (e.g., D20_paired fish vs. D20_unpaired fish; D60_paired fish vs. D60_unpaired fish) using a permutation test, which is implemented by the "exactRankTests" package in R with all default settings. The hypergeometric test (23,306 total genes analyzed) was conducted by phyper function in the “dplyr package” in R.

Clustering

To differentiate transcriptomic responses among the three groups (B, D20, and D60), a principal component analysis (PCA) and a linear discriminant analysis (LDA) were done on the same data set of the top 50% most variable gene transcripts (11,653 genes sorted by standard deviation) from the entire gene set after undergoing DEG analysis in edgeR (n = 23,306 genes). Both the PCA and LDA analyses were implemented using the "MASS" package in R (S4E Fig; source code, LDA). To visualize transcriptomic synchronization, a grade of membership (GoM) analysis was performed using the entire gene set (n = 23,306 genes) with ‘7’ as the selected number of clusters based on the “CountClust” package in R [75] (source code, GoM). To cluster gene expression patterns of each fighting group and each individual, heatmaps were created using the pheatmap (Fig 3C) and ggplot2 (Fig 5B) functions in R. The consensus tree and bootstrap values for each tree node in Fig 3C and Fig 5B were obtained with the hclust function in the “pvclust” package in R/Bioconductor. A Venn diagram was constructed using webtools (http://bioinformatics.psb.ugent.be/webtools/Venn/).

WGCNA

Using normalized expression counts from all the genes that underwent differential expression analysis in edgeR (23,306 genes, described above), we characterized gene expression network dynamics. This analysis was performed using an open-source software, WGCNA package in R [74]. WGCNA shows patterns of co-expressed genes (modules, which are arbitrarily assigned a color name by the software package) based on the eigengene value (a variable derived from the first principal component of expression in a module that represents the gene expression profile in a module, see [74] for more details). Our goal was to identify modules associated with the fighting groups (B, D20, and D60). We calculated Pearson correlation coefficients between all gene pairs in all three fighting groups from the experiment: B, D20, and D60. Further, we used DAVID as described below to perform GO and KEGG pathway analyses on the chosen gene modules in WGCNA. The designed matrix for the WGCNA used to identify gene modules is shown in S15 Table. We used the default parameters for all WGCNA settings (source code, WGCNA).

Characterizing the synchronized gene set

We characterized the synchronized gene set with the following two-step method.

In the first step, the DEGs from the B vs. D20 and B vs. D60 comparisons (FDR < 0.05 and log FC > 0) were obtained as described above.

In the second step, we divided the procedure into three parts as follows.

(1) To determine whether a gene i is synchronized or not, the expression distance for each particular gene between two random fish among the total of 45 possible combinations (45 = 2-permutations of 10, five paired fish and 40 unpaired fish) was separately computed for each group (D20 and D60) using the following formula:

where Dij is the distance for gene i between a random pair j (A and B) and Aij and Bij are the TMM values for gene i of fish A and fish B from the pair j, respectively.

(2) A permutation test for each Di between the five paired fish and the 40 unpaired fish was separately conducted for each of the fighting groups, D20 (D20-paired fish vs. D20-unpaired fish) and D60 (D60-paired fish vs. D60-unpaired fish), to assess the statistical significance. The tests were implemented by the "exactRankTests" package in R with all default settings. A value of p < 0.05 was considered indicative of a synchronized gene.

(3) To determine in which fighting pair(s) gene i is synchronized, we first separately computed the expression distance for each particular gene i between two unpaired fish for each group (D20 and D60) using the following formula:

where Dup is the distance for a gene i between an unpaired fish k (C and D); Cik and Dik are the TMM values for a gene i of fish C and D from the pair k, respectively.

Then, a baseline for each synchronized gene (Bi) was identified by calculating the 25th percentile value for the distance of each particular gene (i) between the 40 unpaired fish for D20 and D60 separately (see Fig 6D, left). Finally, we compared the distance of synchronized gene i in pair j (Dij) with the baseline (Bi). If Dij < Bi, then the synchronized gene i is defined as synchronized in the pair j.

Obtain learning and memory genes

To find genes related to “memory” or “learning”, the sequences from the 1,522 synchronized genes in the D60 group were searched against the amino acid sequences related to “memory” or “learning” obtained from the UniProt database (https://www.uniprot.org/) using blastx with a cut-off E-value of 1e-20.

Gene annotation

The protein sequence table of B. splendens was obtained from http://gigadb.org/dataset/100433. These sequences were used for a BLAST analysis against protein databases to find the homologs of zebrafish (Danio rerio) using Blast2GO (https://www.blast2go.com/blast2go-pro). As a result, we obtained 15,754 well-annotated genes for the GO and KEGG pathways analyses.

Functional gene enrichment analysis

The significantly enriched Gene Ontology terms (biological process) and KEGG pathways were identified by the Database for Annotation Visualization and Integrated Discovery (DAVID) [107]. A p-value cut-off of < 0.05 was used to select for significantly enriched functional terms. The GO terms obtained from DAVID were summarized into larger and more general categories to get a general overview of the underlying biology using REVIGO [108]. Terms were grouped together if they were in a similar pathway and/or based on semantic similarity. GO enrichments along with their respective p-values are in S10S12 Tables

Supporting information

S1 Video [mp4]
Fighting between two males of . .

S1 Fig [a]
Dynamic fighting behaviors of two male . individuals.

S2 Fig [tif]
Behavioral differences among all fighting pairs.

S3 Fig [a]
Gene Ontology and KEGG pathways.

S4 Fig [a]
Correlation coefficient values (r), representative synchronized genes, WGCNA, and PCA.

S5 Fig [a]
Venn diagram for synchronized genes and heatmap for all gene set.

S6 Fig [d]
Brain-transcriptomic Synchronization Associated with Long-Term MemoryGenes.

S7 Fig [tif]
The MDS plot color coded by sequencing methods.

S1 Table [xlsx]
Seventeen typical fighting pairs and those parameters used to construct the representative schematic fighting structure.

S2 Table [xlsx]
Correlation coefficients of biting and striking attacks for eight typical fighting pairs used to visualize behavioral synchronization.

S3 Table [xlsx]
Correlation coefficients of surface-breathing events for eight typical fighting pairs used to visualize behavioral synchronization.

S4 Table [xlsx]
The mapping rate of 25 . brain samples from the non-fighting, D20, and D60 groups to the reference . genome.

S5 Table [xlsx]
Differential Expression Genes List (DEGs) (average TMM value).

S6 Table [xlsx]
Number of genes in each module from the WGCNA analysis.

S7 Table [xlsx]
Top 5 driving genes in each cluster extracted from the Grade of Membership (GoM) analysis.

S8 Table [xlsx]
List of DEGs generated from a comparison of the expression levels from the non-fighting group vs. the D60 group (FDR ≤ 0.05).

S9 Table [xlsx]
List of synchronized genes expressed in one to five pairs of fighting . .

S10 Table [xlsx]
GO terms among the synchronized genes and non-synchronized genes in the D60 group with respect to the biological process (BP) list (p < 0.05).

S11 Table [xlsx]
GO terms among the synchronized genes and non-synchronized genes in the D60 group with respect to the molecular function (MF) list (p < 0.05).

S12 Table [xlsx]
GO terms among the synchronized genes and non-synchronized genes in the D60 group with respect to the cellular component (CC) list (p < 0.05).

S13 Table [xlsx]
KEGG analysis using synchronized genes related to learning or memory and for which there were at least two genes in the pathway (p < 0.05).

S14 Table [xlsx]
Commonly enriched pathways and their constitutive components as represented among the synchronized genes in the D60 group (p < 0.05).

S15 Table [xlsx]
Behavioral matrix for the WGCNA analysis.


Zdroje

1. Darwin C. On the Origin of Species by Means of Natural Selection. Murray, London. 1859.

2. Arnott G, Elwood R. Information gathering and decision making about resource value in animal contests. Animal Behavior. 2008;76(3):529–42.

3. Arnott G, Elwood R. Assessment of fighting ability in animal contests. Anim Behaviour. 2009;77(5):991–1004.

4. Lindenfors P, Tullberg BS. Evolutionary aspects of aggression the importance of sexual selection. Adv Genet. 2011; 75:7–22. doi: 10.1016/B978-0-12-380858-5.00009-5 22078475

5. Taborsky B, Oliveira RF. Social competence: an evolutionary approach. Trends Ecol Evol. 2012;27(12):679–88. doi: 10.1016/j.tree.2012.09.003 23040461

6. Oliveira RF, Simoes JM, Teles MC, Oliveira CR, Becker JD, Lopes JS. Assessment of fight outcome is needed to activate socially driven transcriptional changes in the zebrafish brain. Proc Natl Acad Sci U S A. 2016;113(5): E654–61. doi: 10.1073/pnas.1514292113 26787876

7. Oliveira RF, McGregor PK, Latruffe C, PRSL B, Oliveira RF, McGregor PK, et al. Know thine enemy: fighting fish gather information from observing conspecific interactions. 1998:1045–9.

8. Lorenz K. On aggression. Routledge London & New York. 2005.

9. Kravitz EA, Huber R. Aggression in invertebrates. Curr Opin Neurobiol. 2003;13(6):736–43. doi: 10.1016/j.conb.2003.10.003 14662376

10. Tate M, McGoran RE, White CR. Life in a bubble: the role of the labyrinth organ in determining territory, mating and aggressive behaviours in anabantoids. 2017; 44:723–49.

11. Lin C-p, Lin H-c, Zoology B. Morphological and Biochemical Variations in the Gills of 12 Aquatic Air-Breathing Anabantoid Fish. 2014(May).

12. Monvises A, Nuangsaeng B, Sriwattanarothai N, Panijpan B. The Siamese fighting fish: well-known generally but little-known scientifically. SienceAsia. 2009; 35:8–16.

13. Simpson MJA. The Display of the Siamese Fighting Fish, Betta splendens. Animal Behavior Monographs. 1968;1(1): i-viii, 1–73.

14. Doutrelant C, McGregor PK, Oliveira RF. The effect of an audience on intrasexual communication in male Siamese fighting fish, Betta splendens. Behavioral Ecology. 2001;12(3):283–6.

15. Processes B. Type of intruder and reproductive phase influence male territorial defense in wild- caught Siamese fighting fish, Betta splendens. 2003.

16. Castro N, Ros A, Becker K, Oliveira R. Metabolic costs of aggressive behavior in the Siamese fighting fish, Betta splendens2006. 474–80 p.

17. Alton LA, Portugal SJ, White CR. Balancing the competing requirements of air-breathing and display behaviour during male-male interactions in Siamese fighting fish Betta splendens. Comp Biochem Physiol A Mol Integr Physiol. 2013;164(2):363–7. doi: 10.1016/j.cbpa.2012.11.012 23178457

18. Justus KT, Mendelson TC. Male preference for conspeci fi c mates is stronger than females’ in Betta. Behavioural Processes. 2018;151(February):6–10.

19. Schneider H. Measuring Agonistic Behavior in Zebrafish. Zebrafish Neurobehavioral Protocols. 2011; 51:125–34.

20. Matos RJ, McGregor PK. The effect of the sex of an audience on male-male displays of Siamese fighting fish (Betta splendens). Behaviour. 2002;138(9):1211–21.

21. Dzieweczynski TL, Leopard AK. The effects of stimulus type on consistency of responses to conflicting stimuli in Siamese fighting fish. Behav Processes. 2010;85(2):83–9. doi: 10.1016/j.beproc.2010.06.011 20600699

22. Verbeek P, Iwamoto T, Murakami N. Differences in aggression between wild-type and domesticated fighting fish are context dependent. Anim Behaviour. 2007;73(1):75–83.

23. Larson ET, O'Malley DM, Melloni RH Jr. Aggression and vasotocin are associated with dominant-subordinate relationships in zebrafish. Behav Brain Res. 2006;167(1):94–102. doi: 10.1016/j.bbr.2005.08.020 16213035

24. Clotfelter ED, O'Hare EP, McNitt MM, Carpenter RE, Summers CH. Serotonin decreases aggression via 5-HT1A receptors in the fighting fish Betta splendens. Pharmacol Biochem Behav. 2007;87(2):222–31. doi: 10.1016/j.pbb.2007.04.018 17553555

25. Clotfelter ED, McNitt MM, Carpenter RE, Summers CH. Modulation of monoamine neurotransmitters in fighting fish Betta splendens exposed to waterborne phytoestrogens. 2010:933–43.

26. Yang W, Wang Y, Zhu C, et al. De novo transcriptomic characterization of Betta splendens for identifying sex-biased genes potentially involved in aggressive behavior modulation and EST-SSR maker development. bioRxiv; 2018. doi: 10.1101/355354

27. Renn SC, Aubin-Horth N, Hofmann HA. Fish and chips: functional genomics of social plasticity in an African cichlid fish. J Exp Biol. 2008;211(Pt 18):3041–56. doi: 10.1242/jeb.018242 18775941

28. Filby AL, Paull GC, Hickmore TF, Tyler CR. Unravelling the neurophysiological basis of aggression in a fish model. BMC Genomics. 2010;11:498. doi: 10.1186/1471-2164-11-498 20846403

29. Sneddon LU, Margareto J, Cossins AR. The use of transcriptomics to address questions in behaviour: production of a suppression subtractive hybridisation library from dominance hierarchies of rainbow trout. Physiol Biochem Zool. 2005;78(5):695–705. doi: 10.1086/432141 16052453

30. Greenwood AK, Wark AR, Fernald RD, Hofmann HA. Expression of arginine vasotocin in distinct preoptic regions is associated with dominant and subordinate behaviour in an African cichlid fish. Proc Biol Sci. 2008;275(1649):2393–402. doi: 10.1098/rspb.2008.0622 18628117

31. Kroes RA, Panksepp J, Burgdorf J, Otto NJ, Moskal JR. Modeling depression: social dominance-submission gene expression patterns in rat neocortex. Neuroscience. 2006;137(1):37–49. doi: 10.1016/j.neuroscience.2005.08.076 16289586

32. Fan G, Chan J, Ma K, Yang B, Zhang H, Yang X, et al. Chromosome-level reference genome of the Siamese fighting fish Betta splendens, a model species for the study of aggression. Gigascience. 2018;7(11).

33. Axelrod R, Hamilton WD. The evolution of cooperation. science. 1981;211(4489):1390–6. doi: 10.1126/science.7466396 7466396

34. Bshary R, Grutter AS. Asymmetric cheating opportunities and partner control in a cleaner fish mutualism. Animal Behaviour. 2002;63(3):547–55.

35. Soares MC, Paula JR, Bshary R. Serotonin blockade delays learning performance in a cooperative fish. Animal cognition. 2016;19(5):1027–30. doi: 10.1007/s10071-016-0988-z 27107861

36. Paula JR, Messias JP, Grutter AS, Bshary R, Soares MC. The role of serotonin in the modulation of cooperative behavior. Behavioral Ecology. 2015;26(4):1005–12.

37. Cardoso SC, Paitio JR, Oliveira RF, Bshary R, Soares MC. Arginine vasotocin reduces levels of cooperative behaviour in a cleaner fish. Physiology & behavior. 2015;139:314–20.

38. Kingsbury L, Huang S, Wang J, Gu K, Golshani P, Wu YE. Correlated Neural Activity and Encoding of Behavior across Brains of Socially Interacting Animals Correlated Neural Activity and Encoding of Behavior across Brains of Socially Interacting Animals. Cell. 2019:1–18.

39. Zhang W, Yartsev MM. Article Correlated Neural Activity across the Brains of Socially Interacting Bats. Cell. 2019:1–16.

40. Kandel ER, Schwartz JH, Jessell TM. Principles of Neural Science, Fifth Edition. McGraw-Hill Companies, New York. 2000;4:1227–46.

41. Clayton DF. The genomic action potential. Neurobiol Learn Mem. 2000;74(3):185–216. doi: 10.1006/nlme.2000.3967 11031127

42. Cummings ME, Larkins-Ford J, Reilly CR, Wong RY, Ramsey M, Hofmann HA. Sexual and social stimuli elicit rapid and contrasting genomic responses. Proc Biol Sci. 2008;275(1633):393–402. doi: 10.1098/rspb.2007.1454 18055387

43. Maruska KP, Zhang A, Neboori A, Fernald RD. Social opportunity causes rapid transcriptional changes in the social behaviour network of the brain in an African cichlid fish. J Neuroendocrinol. 2013;25(2):145–57. doi: 10.1111/j.1365-2826.2012.02382.x 22958303

44. Oliveira RF, Silva JF, Simoes JM. Fighting zebrafish: characterization of aggressive behavior and winner-loser effects. Zebrafish. 2011;8(2):73–81. doi: 10.1089/zeb.2011.0690 21612540

45. Alaux C, Sinha S, Hasadsri L, Hunt GJ, Guzman-Novoa E, DeGrandi-Hoffman G, et al. Honey bee aggression supports a link between gene regulation and behavioral evolution. Proc Natl Acad Sci U S A. 2009;106(36):15400–5. doi: 10.1073/pnas.0907043106 19706434

46. Chandrasekaran S, Ament SA, Eddy JA, Rodriguez-Zas SL, Schatz BR, Price ND, et al. Behavior-specific changes in transcriptional modules lead to distinct and predictable neurogenomic states. Proc Natl Acad Sci U S A. 2011;108(44):18020–5. doi: 10.1073/pnas.1114093108 21960440

47. Burmeister SS, Jarvis ED, Fernald RD. Rapid behavioral and genomic responses to social opportunity. PLoS Biol. 2005;3(11): e363. doi: 10.1371/journal.pbio.0030363 16216088

48. Okuno H. Regulation and function of immediate-early genes in the brain: beyond neuronal activity markers. Neuroscience research. 2011;69(3):175–86. doi: 10.1016/j.neures.2010.12.007 21163309

49. Ghosh A, Greenberg ME. Calcium signaling in neurons: molecular mechanisms and cellular consequences. Science. 1995;268(5208):239–47. doi: 10.1126/science.7716515 7716515

50. Iii RJK, Govindarajan A, Jung H-y, Kang H, Tonegawa S. Translational Control by MAPK Signaling in Long-Term Synaptic Plasticity and Memory. 2004; 116:467–79.

51. Robinson GE, Ben-Shahar Y. Social behavior and comparative genomics: new genes or new gene regulation? Genes Brain Behav. 2002;1(4):197–203. doi: 10.1034/j.1601-183x.2002.10401.x 12882364

52. Howard MW, Rizzuto DS, Caplan JB, Madsen JR, Lisman J, Aschenbrenner-Scheibe R, et al. Gamma oscillations correlate with working memory load in humans. Cereb Cortex. 2003;13(12):1369–74. doi: 10.1093/cercor/bhg084 14615302

53. Chou MY, Amo R, Kinoshita M, Cherng BW, Shimazaki H, Agetsuma M, et al. Social conflict resolution regulated by two dorsal habenular subregions in zebrafish. Science. 2016;352(6281):87–90. doi: 10.1126/science.aac9508 27034372

54. Enquist M, Leimar O, Ljungberg T, Mallner Y, Segerdahl N. A test of the sequential assessment game: fighting in the cichlid fish Nannacara anomala. Animal Behavior. 1990;40(1):1–14.

55. Gotceitas V, Godin JGJ. Foraging under the Risk of Predation in Juvenile Atlantic Salmon (Salmo salar L.): Effects of Social Status and Hunger. Behavioral Ecology and Sociobiology. 1991;29(4):255–61.

56. Pellegrino G, Fadiga L, Fogassi L. Understanding motor events: a physiological. Experimental brain …. 1992:176–80.

57. Nakayama S, Johnstone RA, Manica A. Temperament and hunger interact to determine the emergence of leaders in pairs of foraging fish. PLoS One. 2012;7(8):e43747. doi: 10.1371/journal.pone.0043747 22952753

58. Nakayama S, Stumpe MC, Manica A, Johnstone RA. Experience overrides personality differences in the tendency to follow but not in the tendency to lead. Proc Biol Sci. 2013;280(1769):20131724. doi: 10.1098/rspb.2013.1724 23986110

59. Reale D, Reader SM, Sol D, McDougall PT, Dingemanse NJ. Integrating animal temperament within ecology and evolution. Biological reviews of the Cambridge Philosophical Society. 2007;82(2):291–318. doi: 10.1111/j.1469-185X.2007.00010.x 17437562

60. Croft DP, James R, Ward AJ, Botham MS, Mawdsley D, Krause J. Assortative interactions and social networks in fish. Oecologia. 2005;143(2):211–9. doi: 10.1007/s00442-004-1796-8 15682346

61. Firth JA, Cole EF, Ioannou CC, Quinn JL, Aplin LM, Culina A, et al. Personality shapes pair bonding in a wild bird social system. Nat Ecol Evol. 2018;2(11):1696–9. doi: 10.1038/s41559-018-0670-8 30275466

62. Bekinschtein P, Cammarota M, Igaz LM, Bevilaqua LR, Izquierdo I, Medina JH. Persistence of long-term memory storage requires a late protein synthesis- and BDNF- dependent phase in the hippocampus. Neuron. 2007;53(2):261–77. doi: 10.1016/j.neuron.2006.11.025 17224407

63. Safe S, Jin UH, Morpurgo B, Abudayyeh A, Singh M, Tjalkens RB. Nuclear receptor 4A (NR4A) family—orphans no more. J Steroid Biochem Mol Biol. 2016;157:48–60. doi: 10.1016/j.jsbmb.2015.04.016 25917081

64. Malenka RC, Nicoll RA. Long-term potentiation—a decade of progress? Science. 1999;285(5435):1870–4. doi: 10.1126/science.285.5435.1870 10489359

65. Cui WW, Low SE, Hirata H, Saint-Amant L, Geisler R, Hume RI, et al. The zebrafish shocked gene encodes a glycine transporter and is essential for the function of early neural circuits in the CNS. J Neurosci. 2005;25(28):6610–20. doi: 10.1523/JNEUROSCI.5009-04.2005 16014722

66. Mills BN, Albert GP, Halterman MW. Expression profiling of the MAP kinase phosphatase family reveals a role for DUSP1 in the glioblastoma stem cell niche. Cancer Microenvironment. 2017;10(1–3):57–68. doi: 10.1007/s12307-017-0197-6 28822081

67. Dolmetsch RE, Pajvani U, Fife K, Spotts JM, & Greenberg ME (2001) Signaling to the nucleus by an L-type calcium channel-calmodulin complex through the MAP kinase pathway. Science 294(5541):333–339. doi: 10.1126/science.1063395 11598293

68. Wu J, Huang KP, Huang FL. Participation of NMDA-mediated phosphorylation and oxidation of neurogranin in the regulation of Ca2+- and Ca2+/calmodulin-dependent neuronal signaling in the hippocampus. J Neurochem. 2003;86(6):1524–33. doi: 10.1046/j.1471-4159.2003.01963.x 12950461

69. Chi B, Wang Q, Wu G, Tan M, Wang L, Shi M, et al. Aly and THO are required for assembly of the human TREX complex and association of TREX components with the spliced mRNA. Nucleic acids research. 2013;41(2):1294–306. doi: 10.1093/nar/gks1188 23222130

70. Bukhari SA, Saul MC, Seward CH, Zhang H, Bensky M, James N, et al. Temporal dynamics of neurogenomic plasticity in response to social interactions in male threespined sticklebacks. PLoS genetics. 2017;13(7).

71. Saul MC, Seward CH, Troy JM, Zhang H, Sloofman LG, Lu X, et al. Transcriptional regulatory dynamics drive coordinated metabolic and neural response to social challenge in mice. Genome research. 2017;27(6):959–72. doi: 10.1101/gr.214221.116 28356321

72. Saul M.C., Blatti C., Yang W., Bukhari S.A., Shpigler H.Y., Troy J.M., Seward C.H., Sloofman L., Chandrasekaran S., Bell A.M. and Stubbs L., 2019. Cross‐species systems analysis of evolutionary toolkits of neurogenomic response to social challenge. Genes, Brain and Behavior, 18(1), p.e12502.

73. Wong RY, Hofmann HA. Behavioural genomics: an organismic perspective. e LS. 2001 May 30.

74. Langfelder P, Horvath S. WGCNA: an R package for weighted correlation network analysis. BMC bioinformatics. 2008;9:559. doi: 10.1186/1471-2105-9-559 19114008

75. Dey KK, Hsiao CJ, Stephens M. Visualizing the structure of RNA-seq expression data using grade of membership models. PLoS Genet. 2017;13(3):e1006599. doi: 10.1371/journal.pgen.1006599 28333934

76. Whitfield CW, Cziko A-M, Robinson GE. Gene expression profiles in the brain predict behavior in individual honey bees. Science. 2003;302(5643):296–9. doi: 10.1126/science.1086807 14551438

77. Robinson GE, Fernald RD, Clayton DF. Genes and social behavior. science. 2008;322(5903):896–900. doi: 10.1126/science.1159277 18988841

78. Hofmann HA. Functional genomics of neural and behavioral plasticity. Journal of neurobiology. 2003 Jan;54(1):272–82. doi: 10.1002/neu.10172 12486709

79. Kempadoo KA, Mosharov EV, Choi SJ, Sulzer D, Kandel ER. Dopamine release from the locus coeruleus to the dorsal hippocampus promotes spatial learning and memory. Proc Natl Acad Sci U S A. 2016;113(51):14835–40. doi: 10.1073/pnas.1616515114 27930324

80. Lee JL, Hynds RE. Divergent cellular pathways of hippocampal memory consolidation and reconsolidation. Hippocampus. 2013;23(3):233–44. doi: 10.1002/hipo.22083 23197404

81. Teles MC, Cardoso SD, Oliveira RF. Social plasticity relies on different neuroplasticity mechanisms across the brain social decision-making network in zebrafish. Frontiers in behavioral neuroscience. 2016;10:16. doi: 10.3389/fnbeh.2016.00016 26909029

82. Wong RY, Cummings ME. Expression patterns of neuroligin-3 and tyrosine hydroxylase across the brain in mate choice contexts in female swordtails. Brain, behavior and evolution. 2014;83(3):231–43. doi: 10.1159/000360071 24854097

83. Salamone JD, Yohn SE, Lopez-Cruz L, San Miguel N, Correa M. Activational and effort-related aspects of motivation: neural mechanisms and implications for psychopathology. Brain. 2016;139(Pt 5):1325–47. doi: 10.1093/brain/aww050 27189581

84. O’Connell L.A. and Hofmann H.A., 2012. Evolution of a vertebrate social decision-making network. Science, 336(6085), pp.1154–1157. doi: 10.1126/science.1218889 22654056

85. Desjardins JK, Fernald RD. What do fish make of mirror images? Biology Letters. 2010;6(6):744–7. doi: 10.1098/rsbl.2010.0247 20462889

86. Yoshida K, Saito N, Iriki A, Isoda M. Representation of others' action by neurons in monkey medial frontal cortex. Curr Biol. 2011;21(3):249–53. doi: 10.1016/j.cub.2011.01.004 21256015

87. Breveglieri R, Vaccari FE, Bosco A, Gamberini M, Fattori P, Galletti C. Neurons Modulated by Action Execution and Observation in the Macaque Medial Parietal Cortex. Curr Biol. 2019;29(7):1218–25 e3. doi: 10.1016/j.cub.2019.02.027 30880012

88. Cattaneo L, Fabbri-Destro M, Boria S, Pieraccini C, Monti A, Cossu G, et al. Impairment of actions chains in autism and its possible role in intention understanding. Proc Natl Acad Sci U S A. 2007;104(45):17825–30. doi: 10.1073/pnas.0706273104 17965234

89. Decety J., Bartal I.B.A., Uzefovsky F. and Knafo-Noam A., 2016. Empathy as a driver of prosocial behaviour: highly conserved neurobehavioural mechanisms across species. Philosophical Transactions of the Royal Society B: Biological Sciences, 371(1686), p.20150077.

90. Bailey MR, Simpson EH, Balsam PD. Neural substrates underlying effort, time, and risk-based decision making in motivated behavior. Neurobiol Learn Mem. 2016;133:233–56. doi: 10.1016/j.nlm.2016.07.015 27427327

91. Parkinson C, Kleinbaum AM, Wheatley T. Similar neural responses predict friendship. Nature Communications. 2018.

92. Piazza EA, Hasenfratz L, Hasson U, Lew-Williams C. Infant and adult brains are coupled to the dynamics of natural communication. Psychological Science. 2020;31(1):6–17.

93. Adriaense JEC, Martin JS, Schiestl M, Lamm C, Bugnyar T. Negative emotional contagion and cognitive bias in common ravens (Corvus corax). 2019:1–6.

94. Tyssowski K.M., DeStefino N.R., Cho J.H., Dunn C.J., Poston R.G., Carty C.E., Jones R.D., Chang S.M., Romeo P., Wurzelmann M.K. and Ward J.M., 2018. Different neuronal activity patterns induce different gene expression programs. Neuron, 98(3), pp.530–546. doi: 10.1016/j.neuron.2018.04.001 29681534

95. Goto Y, Grace AA. Dopaminergic modulation of limbic and cortical drive of nucleus accumbens in goal-directed behavior. Nat Neurosci. 2005;8(6):805–12. doi: 10.1038/nn1471 15908948

96. Huber R, Smith K, Delago A, Isaksson K, Kravitz EA. Serotonin and aggressive motivation in crustaceans: altering the decision to retreat. Proc Natl Acad Sci U S A. 1997;94(11):5939–42. doi: 10.1073/pnas.94.11.5939 9159179

97. Saul MC, Majdak P, Perez S, Reilly M, Garland T Jr., Rhodes JS. High motivation for exercise is associated with altered chromatin regulators of monoamine receptor gene expression in the striatum of selectively bred mice. Genes Brain Behav. 2017;16(3):328–41. doi: 10.1111/gbb.12347 27749013

98. Bailey MR, Goldman O, Bello EP, Chohan MO, Jeong N, Winiger V, et al. An Interaction between Serotonin Receptor Signaling and Dopamine Enhances Goal-Directed Vigor and Persistence in Mice. J Neurosci. 2018;38(9):2149–62. doi: 10.1523/JNEUROSCI.2088-17.2018 29367407

99. Dowd M, Joy R. Estimating behavioral parameters in animal movement models using a state‐augmented particle filter. Ecology. 2011;92(3):568–75. doi: 10.1890/10-0611.1 21608465

100. Sengupta S, Bolin JM, Ruotti V, Nguyen BK, Thomson JA, Elwell AL, et al. Single read and paired end mRNA-Seq Illumina libraries from 10 nanograms total RNA. JoVE (Journal of Visualized Experiments). 2011(56):e3340.

101. Embnet M. Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnetjournal. 1994;17(1):10–2.

102. Trapnell C, Pachter L, Salzberg SL. TopHat: discovering splice junctions with RNA-Seq. Bioinformatics. 2009;25(9):1105–11. doi: 10.1093/bioinformatics/btp120 19289445

103. Langmead B, Salzberg SL. Fast gapped-read alignment with Bowtie 2. Nat Methods. 2012;9(4):357–9. doi: 10.1038/nmeth.1923 22388286

104. Liao Y, Smyth GK, Shi W. featureCounts: an efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics. 2014;30(7):923–30. doi: 10.1093/bioinformatics/btt656 24227677

105. Robinson MD, Oshlack A. A scaling normalization method for differential expression analysis of RNA-seq data. Genome Biol. 2010;11(3):R25. doi: 10.1186/gb-2010-11-3-r25 20196867

106. Robinson MD, McCarthy DJ, Smyth GK. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2010;26(1):139–40. doi: 10.1093/bioinformatics/btp616 19910308

107. Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nature protocols. 2009;4(1):44. doi: 10.1038/nprot.2008.211 19131956

108. Supek F, Bošnjak M, Škunca N, Šmuc T. REVIGO summarizes and visualizes long lists of gene ontology terms. PloS one. 2011;6(7).


Článek vyšel v časopise

PLOS Genetics


2020 Číslo 6
Nejčtenější tento týden
Nejčtenější v tomto čísle
Kurzy

Zvyšte si kvalifikaci online z pohodlí domova

Důležitost adherence při depresivním onemocnění
nový kurz
Autoři: MUDr. Eliška Bartečková, Ph.D.

Koncepce osteologické péče pro gynekology a praktické lékaře
Autoři: MUDr. František Šenk

Sekvenční léčba schizofrenie
Autoři: MUDr. Jana Hořínková, Ph.D.

Hypertenze a hypercholesterolémie – synergický efekt léčby
Autoři: prof. MUDr. Hana Rosolová, DrSc.

Multidisciplinární zkušenosti u pacientů s diabetem
Autoři: Prof. MUDr. Martin Haluzík, DrSc., prof. MUDr. Vojtěch Melenovský, CSc., prof. MUDr. Vladimír Tesař, DrSc.

Všechny kurzy
Přihlášení
Zapomenuté heslo

Zadejte e-mailovou adresu, se kterou jste vytvářel(a) účet, budou Vám na ni zaslány informace k nastavení nového hesla.

Přihlášení

Nemáte účet?  Registrujte se

#ADS_BOTTOM_SCRIPTS#