Next Article in Journal
Comparison of High- and Low-LET Radiation-Induced DNA Double-Strand Break Processing in Living Cells
Next Article in Special Issue
Immunohistochemistry and Mutation Analysis of SDHx Genes in Carotid Paragangliomas
Previous Article in Journal
A Multi-Endpoint Approach to Base Excision Repair Incision Activity Augmented by PARylation and DNA Damage Levels in Mice: Impact of Sex and Age
Previous Article in Special Issue
Potential Associations Among Alteration of Salivary miRNAs, Saliva Microbiome Structure, and Cognitive Impairments in Autistic Children
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 21
Issue 18
10.3390/ijms21186599
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Gene Expression Changes in the Ventral Tegmental Area of Male Mice with Alternative Social Behavior Experience in Chronic Agonistic Interactions

by
Olga Redina
1,*,
Vladimir Babenko
1,2,
Dmitry Smagin
1,
Irina Kovalenko
1,
Anna Galyamina
1,
Vadim Efimov
1,2 and
Natalia Kudryavtseva
1
1
FRC Institute of Cytology and Genetics, Siberian Branch of Russian Academy of Sciences, 630090 Novosibirsk, Russia
2
Department of Natural Sciences, Novosibirsk State University, 630090 Novosibirsk, Russia
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2020, 21(18), 6599; https://doi.org/10.3390/ijms21186599
Submission received: 27 July 2020 / Revised: 3 September 2020 / Accepted: 4 September 2020 / Published: 9 September 2020
(This article belongs to the Special Issue Medical Genetics, Genomics and Bioinformatics – 2020)
Browse Figures
Review Reports Versions Notes

Abstract

:
Daily agonistic interactions of mice are an effective experimental approach to elucidate the molecular mechanisms underlying the excitation of the brain neurons and the formation of alternative social behavior patterns. An RNA-Seq analysis was used to compare the ventral tegmental area (VTA) transcriptome profiles for three groups of male C57BL/6J mice: winners, a group of chronically winning mice, losers, a group of chronically defeated mice, and controls. The data obtained show that both winners and defeated mice experience stress, which however, has a more drastic effect on defeated animals causing more significant changes in the levels of gene transcription. Four genes (Nrgn, Ercc2, Otx2, and Six3) changed their VTA expression profiles in opposite directions in winners and defeated mice. It was first shown that Nrgn (neurogranin) expression was highly correlated with the expression of the genes involved in dopamine synthesis and transport (Th, Ddc, Slc6a3, and Drd2) in the VTA of defeated mice but not in winners. The obtained network of 31 coregulated genes, encoding proteins associated with nervous system development (including 24 genes associated with the generation of neurons), may be potentially useful for studying their role in the VTA dopaminergic neurons maturation under the influence of social stress.

1. Introduction

The model of repeated agonistic interactions (model of chronic social conflicts) is known as an effective approach for studying the physiological and molecular mechanisms involved in the formation of distinctive behavior patterns in mice under aggressive social encounters [1,2,3,4]. Under experimental exposure to daily agonistic interactions, two groups of animals are formed: (1) with experience of social victories (winners) and (2) with chronic experience of defeats (defeated mice or losers). Mice with repeated winning experiences in everyday agonistic interactions develop marked aggressiveness, anxiety, impulsivity, and impaired motivated and cognitive behaviors [3]. Mice with chronic experience of social defeats develop mixed anxiety/depression-like states [5,6].
In previous publications of our group, it was shown that stress resulting from daily agonistic interactions causes crucial changes in the expression of many genes in various brain regions of adult mice [5,7,8,9,10,11,12,13]. Earlier, we described in detail changes in the transcription of the genes associated with the neurotransmitter systems functioning in the ventral tegmental area (VTA) of defeated male mice as compared to the controls that did not have an agonistic interaction experience [9].
The VTA is known as a brain region from which the mesocortical and mesolimbic dopaminergic pathways originate. More than half of the VTA neurons are dopaminergic, however, non-dopamine-releasing VTA neurons also contribute significantly to VTA functioning. Both dopamine-releasing and non-dopamine-releasing VTA neurons participate in the generation of output firing patterns by the integration of afferent signals with local inhibitory or excitatory inputs [14]. Ascending dopaminergic fibers of VTA neurons are projected onto multiple cortical and subcortical regions, mainly the nucleus accumbens and medial prefrontal cortex [15]. VTA neurons are thought to play a central role in reward circuits, attention, and motivational, emotional and addictive behaviors. It should be emphasized that due to functional and anatomical heterogeneity, VTA dopamine neurons are activated by opposite emotional stimuli, such as reward and aversion. The induction of the VTA dopaminergic activity is considered to play a critical role in the development of a number of psychiatric disorders under both acute and repeated stress [16,17].
The studying of the molecular mechanisms underlying the excitation of VTA neurons, in the context of the formation of alternative forms of social behavior under long-term negative and positive experiences in daily agonistic interactions, may help to identify the genes and the metabolic pathways contributing the most to the alternative behavioral phenotypes such as depression-like states and pathological aggression.
The current paper provides a comparative analysis of differentially expressed VTA genes (DEGs) in winners and defeated animals (losers) compared with controls in order to identify the molecular determinants underlying the formation of alternative behavior patterns in male mice under daily agonistic interactions.

2. Results

2.1. Differentially Expressed Genes (DEGs) in the VTA of Winners vs. Control Mice

A total of 13,750 expressed genes were identified. Among them, 44 genes were characterized as being differentially expressed in the VTA of aggressive mice (winners) compared to controls (Supplementary Materials Table S1). Of these, 28 genes (63.6%) had a lower level of transcription in the winners compared to controls. A functional annotation of DEGs revealed 14 genes associated with a behavior and neurological phenotype (Table 1). One of these genes (Tph2) is known as associated with abnormal aggression-related behavior and increased aggression towards mice. Three genes from the list (Slc17a7, Slc6a4 and Tph2) are associated with the terms abnormal fear/anxiety-related behavior and abnormal depression-related behavior. Increased expression levels of transcription of the Tph2 and Slc6a4 genes suggest an increase in the synthesis and transport of serotonin in the VTA of winning mice compared to controls.
A total of 9 out of 44 DEGs encode transcription factors (Table 2), five of which (Ercc2, Lbx1, Maf, Nkx6-1, Tcf7l2) are associated with the behavior/neurological phenotype and exhibit reduced expression levels in the VTA of winners. An analysis carried out in the Kyoto Encyclopedia of Genes and Genomes (KEGG) database did not reveal metabolic pathways significantly enriched by DEGs.

2.2. DEGs in the VTA of Losers vs. Controls

A total of 13,656 expressed genes were identified, among which 188 genes were characterized as differentially expressed in the VTA of defeated mice (losers) compared to controls (Supplementary Materials Table S2). Most DEGs, namely 112 genes (59.6%), showed a decreased level of transcription in defeated mice compared to controls. Altogether, 64 genes associated with the behavior/neurological phenotype were found among the DEGs (Table 1). Six of these genes (Camk2a, Egr1, En2, Fev, Gad2, Tph2) are known to be associated with abnormal aggression-related behavior and an increased aggression towards mice. Sixteen genes from the list in Table 1 are associated with abnormal fear/anxiety-related behavior, six of which (Camk2a, Gad2, Slc6a2, Slc6a3, Slc6a4, Tph2) are also known to be involved in abnormal depression-related behavior.
A total of 28 out of 188 DEGs encode transcription factors (Table 2), most of which (16 genes) are expressed at reduced levels in the VTA of defeated mice. Among the DEGs encoding transcription factors, 16 genes are associated with the behavior/neurological phenotype.
The analysis of DEGs detected by comparing the level of gene expression in defeated mice to controls revealed several metabolic pathways (Table 3). The expression characteristics of the DEGs linked with these metabolic pathways are presented in Supplementary Materials Table S3. According to the changes in the level of transcription of genes related to the metabolic pathways in defeated animals compared with control animals, the activation of GABAergic synapses should be assumed, which includes the activation of GABA synthesis (Gad1, Gad2) and GABA receptors (Gabra1, Gabrb2, Gabrg2); in addition, activation of the synthesis (Th, Ddc) and transport (Slc6a3) of dopamine can be assumed. These genes are also implicated in drug (nicotine, amphetamine, cocaine) dependence.

2.3. DEGs with a Changed Expression Both in Winners and Losers (Common DEGs)

A comparison of the corresponding DEG lists revealed 23 common genes which significantly changed their expression in the VTA of both winners and losers compared to controls (Table 4).
It can be seen that 22 of these genes either increased or decreased the level of transcription unidirectionally in both experimental groups, and only one gene (Nrgn, neurogranin) changed the level of transcription in winners and defeated animals in different directions with respect to the controls. In the VTA of winners, the Nrgn gene transcription level was decreased, and in the defeated mice the Nrgn gene transcription level was increased compared to control animals. Therefore, it can be concluded that, regardless of whether the animals win or are defeated, the level of transcription of many genes in the VTA of experimental mice changes unidirectionally, which is probably the result of the reaction of the animal organism to conditions of social stress. It can be assumed that genes whose expression level varies in different directions can make the most significant contribution to the formation of alternative (aggressive or depressive) behaviors. According to our results, Nrgn is one of such genes.
To identify more genes that changed the transcription levels in the VTA of the winners and losers in opposite directions, a comparative analysis of gene expression in these two experimental groups was performed.

2.4. DEGs in the VTA of Winners vs. Losers

When comparing the level of gene expression in the VTA of winners and losers, a total of 207 DEGs (Supplementary Materials Table S4) were identified, of which 89 genes (43.0%) had a decreased transcription level in defeated mice.
The significantly enriched metabolic pathways and expression characteristics of the related DEGs are shown in Table 5.
The changes in gene expression levels suggest an activation of the dopaminergic metabolic pathway in the VTA of defeated mice, differences in calcium signaling, and activation of several genes related to the axon guidance metabolic pathway, representing a key stage in the formation of neuronal networks, as well as a decrease in the expression of several extracellular matrix receptor interaction genes.
An analysis of DEG expressions in winners vs. losers, revealed six genes with an expression that level changed in opposite directions relative to controls (Table 6). Four of them encode transcription factors. Several genes (Ercc2, Nrgn, Otx2, and Tcf7l2) are annotated as being related to the behavior/neurological phenotype, and Ercc2 as involved in abnormal emotional/affective behaviors. Among these DEGs (listed in Table 6), there may be the genes whose expression patterns determine behavioral differences between winners and losers.

2.5. DEGs that Can Maximize the Differences between the Groups of Winners and Losers

The genes with the most significant contribution to the formation of intergroup differences were sought using a partial-least squares discriminant analysis (PLS-DA). The location of the mice from the experimental groups (winners and losers) in the constructed PLS-DA axes is shown in Figure 1A. Next, a correlation analysis between the obtained coordinates of the mice and the level of gene expression in the winners and losers was performed. A histogram showing the distribution of expressed genes according to the obtained values of correlation coefficients is given in Figure 1B. Genes characterized by maximum correlation coefficients are believed to contribute most significantly to intergroup differences.
For further analysis in the KEGG database, 506 genes characterized by correlation coefficients greater than 0.90 or −0.90 (the list of these genes is presented in Supplementary Materials Table S5) were selected. Their analysis revealed the significance of several metabolic pathways (Table 7), among which melanogenesis and calcium signaling pathways are noted as the most significantly enriched. Among the DEGs assigned to the identified metabolic pathways, the most common gene is Camk2a, encoding calcium/calmodulin-dependent protein kinase II alpha.
Camk2a is one of the 63 DEGs characterized by correlation coefficients greater than 0.90 or −0.90. The detailed description of these DEGs is given in Supplementary Materials Table S6. The list of these genes includes 23 genes associated with the behavior/neurological phenotype, six of which (Ercc2, Bcl11b, Pitx2, Otx2, En1, Foxp2) encode transcription factors. Supplementary Materials Table S6 presents a total of 17 genes encoding transcription factors.
The list of genes presented in Supplementary Materials Table S6 includes four genes (Six3, Ercc2, Otx2, and Nrgn) whose expression levels in the winners and defeated mice changed in opposite directions. DEGs bidirectionally changing the level of transcription in winners and defeated mice and associated with the behavior/neurological phenotype (Ercc2, Otx2, Nrgn) can be considered as the key genes in the formation of alternative behaviors in the conditions of social stress. It should be noted that the most highly expressed of these genes are the Ercc2 and Nrgn genes.
According to the Gene Ontology, the Ercc2 gene is associated with biological processes such as neurogenesis (gliogenesis), the regulation of transcription from RNA polymerase II promoter (regulation of gene expression), programmed cell death, stem cell differentiation, response to stimulus, response to hypoxia, immune system development, and aging, and the Nrgn gene is associated with biological processes such as nervous system development, cell–cell signaling, trans-synaptic signaling, behavior, cognition, the regulation of signaling, learning or memory, response to a stimulus, signal transduction, regulation of synapse structure or activity, modulation of synaptic transmission, regulation of synaptic plasticity, associative learning, and intracellular signal transduction.
A list of 63 DEGs presented in Supplementary Materials Table S6 was used in STRING database utilities to identify possible associations between the encoded proteins. According to the results presented in Figure 2, the network includes 31 genes (shown in red) encoding proteins associated with nervous system development, 24 of which are associated with the generation of neurons (shown in blue).
Among them, several nodes were identified. One of the main nodes is the protein encoded by the Camk2a gene. An analysis revealed the interaction of Nrgn with Camk2a, which in its turn is linked to a group of proteins associated with a number of ionic (calcium and sodium) channels. In addition, Nrgn may be associated with glycine transporters (Slc6a5). In addition to Camk2a, proteins encoded by the Six3 and Otx2 genes can also be considered nodal ones. Since the expression of these genes, like the expression of the Nrgn gene, is changed in opposite directions in winners and defeated mice compared to controls, they can be considered as potential candidates for further studies.
Thus, we identified four genes (Nrgn, Ercc2, Otx2, and Six3) that change the level of transcription in opposite directions in winners and losers, thereby contributing to the formation of intergroup differences, including, possibly, the formation of alternative forms of behavior.

2.6. Correlation Between the Expression of Nrgn, Ercc2, Otx2, and Six3 and the Genes Involved in the Synthesis and Transport of Dopamine in the VTA Was Determined by a Cluster Analysis Using Expression Data for the Genes Differentially Expressed in the VTA of Winners and Losers

The results are presented in Supplementary Materials File S7. As can be seen, these four genes and Th, Ddc, and Slc6a3 genes encoding tyrosine hydroxylase, DOPA decarboxylase, and dopamine transporters, respectively, are in the same cluster. Therefore, it can be assumed that genes with oppositely changed expressions in winners and losers may be coregulated or functionally associated with the synthesis and transport of dopamine in the VTA.
According to our data, two genes encoding dopamine receptors, Drd2 and Drd5, are expressed in the VTA of mice. Based on the criteria of the Cufflinks/Cuffdiff packages used in our work, the Drd2 gene was not characterized as a DEG in all the comparisons; however, its transcription level increased 1.55 times in the VTA of defeated mice as compared to controls (p < 0.009) and only 1.02 times in the VTA of winners (p < 0.91) as compared to controls, which remained practically unchanged. Therefore, it can be assumed that a change in the level of Drd2 transcription may play a functionally important role in the VTA of defeated mice. The expression level of the Drd5 gene remained stable in all three experimental groups.
To understand whether changes in the level of expression of the genes characterized above as key ones can be related to alternative dopaminergic signaling cascades in winners and losers, the correlation between their level of expression and expression levels of the Drd2, Th, Ddc, and Slc6a3 genes was calculated in both experimental groups (Table 8). The results showed a strong correlation between the expression of Nrgn and the above four genes in defeated animals. In contrast to the group of winners, a statistically significant correlation between the expression of Nrgn and Ercc2 and dopaminergic genes was obtained in the VTA of defeated animals. Three nodal genes (Camk2a, Otx2, and Six3) were characterized by high correlation coefficients with the dopaminergic genes in the VTA of mice from both experimental groups.

3. Discussion

First of all, when discussing the results of this study, it is important to comment on the alterations in the expression levels of the genes involved in the synthesis or transport of neurotransmitters. The analysis showed that the Tph2 and Slc6a4 genes responsible for the synthesis and transport of serotonin, respectively, in the VTA are expressed at higher levels in both winners and losers compared to the controls. The role of serotonin in both aggression and anxiety/depression has been extensively studied and reviewed by many authors [18,19,20]. Although serotonin is believed to be a potent inhibitor of aggressive behavior, there are numerous data in favor and against this hypothesis (reviewed in [21]). The complexity of the mechanisms regulating serotonin-mediated behavior is noted in many studies [19,22]. It is reported that the pharmacological effects of drugs depend on age, environment or experimental conditions [18,21]. Therefore, when discussing the results of this study, one should bear in mind that they are mostly relevant to the development of mixed anxiety–depressive disorder and pathological aggressive behavior in mice within the framework of the used experimental model.
In the current study, in the experimental groups of mice with alternative types of behavior, unidirectional changes in the expression of serotonergic genes indicate that the synthesis and transport of serotonin are activated as a reaction of the serotonergic system to experimentally induced social stress experienced by animals and probably are not associated with the formation of alternative forms of behavior. This view is in good agreement with the opinion of other authors who believe that changes in serotonin levels associated with abnormal behavior are conditioned by different factors [19]. It should also be noted here that the defeated animals are likely to experience more severe stress as indicated by a significantly higher increase in the transcription of serotonergic genes in the VTA of defeated mice than in the winners compared to controls (Table 4).
As already noted in the introduction, most VTA neurons are dopaminergic. Statistically significant increases in the transcription levels of several genes (Ddc, Slc6a3, Th) controlling dopamine synthesis (Th, Ddc) and transport (Slc6a3) were detected in the VTA of defeated mice as compared to controls. In winning mice, the respective changes were less significant. Compared to controls, in the VTA of winning animals the increase in the level of transcription of the Ddc gene was characterized by a p value = 0.0028, whereas for Slc6a3 and Th it was statistically insignificant according to the annotations of the Cuffdiff program. At the same time, we previously studied the expression of these genes in larger groups of animals by real-time PCR and showed that in the VTA of the winning mice the level of transcription of the Slc6a3 and Th genes was significantly increased [23]. Therefore, it can be assumed that in the VTA of both winning (aggressive) and defeated mice, dopamine synthesis increases, however, as in the case of serotonin, the effect of stress on dopamine synthesis and transport in defeated mice is stronger than that in the winners.
GABAergic neurons also play an important role in VTA functioning [14]. Several genes (Gabra1, Gabrb2, Gabrg2, Gad1, Gad2) related to the functioning of GABAergic synapses demonstrated a significantly increased transcription in defeated mice compared to controls (Table 4 and Supplementary Materials Table S2). According to our RNA-Seq data, the level of transcription of all these five genes also increased in the group of aggressive animals, but the differences compared to the control did not reach statistically significant values. Accordingly, we can assume both an increase in the synthesis of GABA and an increase in the expression of its receptors (Gabra1, Gabrb2, Gabrg2) in both experimental groups of mice, but in aggressive animals, changes in the function of GABAergic synapses were less pronounced than in the defeated ones. The possible functional consequences of significant changes in these genes in defeated animals were discussed in [9]. Here we want to emphasize that the unidirectional activation of a number of genes that control the functioning of GABAergic neurons occurs in the VTA of mice of both experimental groups and that this activation is more pronounced in defeated animals. The fact that defeated mice are more stressed than winners also follows from almost four times the number of DEGs that were found in losers vs. controls than in winners vs. controls comparisons. These results are in very good agreement with our observation that both winners and defeated mice are characterized by increased anxiety [6,24,25,26].
It should be noted that in the current study the social stress experienced by animals caused a change in the expression level of multiple genes that encode transcription factors in both experimental groups of mice, suggesting a wide range of functional changes in the VTA and, with a high probability, in the areas of the brain onto which the VTA neurons project. As can be seen from Table 4, which represents common DEGs, three DEGs encoding transcription factors (Lbx1, Maf, Nkx6-1) reduced the level of transcription in both winners and defeated mice, therefore suggesting that these genes may play a key role in regulating the overall response to stress.
Almost all genes listed in Table 4 have their transcription levels changed in the same direction in both winners and defeated mice, and only one (Nrgn) has it changed in the opposite direction—its transcription level decreased in winners and increased in defeated mice. The Nrgn gene encodes neurogranin (Ng), which is enriched at dendritic spines and can enhance synaptic strength by targeting calmodulin (CaM) [27]. The bidirectional changes in the level of Nrgn transcription found in our experiment are in good agreement with the fact that neurogranin is associated with both long-term potentiation and long-term depression [28]. Thus, it can be assumed that this particular gene can be involved in the formation of alternative behaviors in the experimental groups of mice.
In physiological conditions, neurogranin forms a complex with CaM, and its CaM-binding affinity can be modulated by phosphorylation, oxidation and glutathiolation under the activation of protein kinase C or oxidant stress [29]. It was shown that mice with an Nrgn gene deletion displayed an apparently normal phenotype but were characterized by impairments in spatial and emotional learning and by changes in hippocampal short- and long-term plasticity (paired-pulse depression, synaptic fatigue, and long-term potentiation induction). These deficits were accompanied by a decreased basal level of the activated Ca(2+)/CaM-dependent kinase II (CaMKII) [30,31]. The behavioral testing revealed that neurogranin knockout mice were both hyperactive and socially withdrawn [32]. Neurogranin overexpression in neurons enhances synaptic strength as well as long-term potentiation (LTP) and learning by promoting calcium-mediated signaling [33]. Based on the above information, we can assume that the increase in Nrgn transcription in the VTA of defeated mice and its decrease in the VTA of winners observed in the current study are adaptive reactions of the system working to achieve a homeostatic balance.
It was shown that the restorative effects of neurogranin on synaptic depression and LTP deficits are dependent on the interaction of neurogranin and CaM and CaM-dependent activation of CaMKII [34]. The results of our experiment showed a coordinated increase in the expression of both Nrgn and Camk2a (encoding calcium/calmodulin-dependent protein kinase II alpha) genes in the VTA of defeated mice, while the level of Camk2a transcription was unchanged in the VTA of the winners exhibiting a decreased level of Nrgn transcription.
CaMKII is known as a mediator of Ca2+-linked signaling. It phosphorylates a wide range of substrates playing a multifunctional role in the coordination and regulation of Ca2+-mediated alterations in different intracellular events related to the neuronal functions: the synthesis and release of neurotransmitters, modulation of ion channel activity, cellular transport, cell morphology and neurite extension, synaptic plasticity, gene expression, etc. [35,36,37]. Camk2a encodes a calcium/calmodulin-dependent protein kinase subunit that plays a key role in synaptic plasticity, AMPA receptor transmission, LTP, and long-term memory formation, and its dysfunction underlies neuropsychiatric disorders such as drug addiction, schizophrenia, depression, epilepsy, and multiple neurodevelopmental disorders, perhaps through maladaptations in glutamate signaling and neuroplasticity [38,39]. The data obtained in the current study are in a good agreement with the findings that calcium signaling is an important link in the regulation of processes associated with neurogenesis in various neurological disorders [40].
In our work, when comparing the transcriptomes of the winning and defeated mice, the Camk2a gene was assigned to the DEGs with maximum contributions to the intergroup differences (Supplementary Materials Table S6). Figure 2 clearly demonstrates that the protein encoded by the Camk2a gene can be viewed as a nodal gene in the calcium modulating pathway, as well as in the regulation of sodium and potassium ion channels and a number of other processes responsible for functional differences in the VTA of winners and defeated animals.
The list of DEGs characterized as making a maximum contribution to intergroup differences contains three more genes (Ercc2, Otx2, Six3) that change the transcription level in winners and defeated mice in opposite directions. These three genes encode transcription factors.
The Ercc2 gene is associated with the behavior/neurological phenotype. In the VTA, a decrease in the level of transcription of this gene was detected in the winners, and an increase in its transcription was detected in the defeated mice (Table 6). Given that Ercc2 is involved in repairing the damage caused by the redox process [41], it can be assumed that oxidative stress processes in the VTA play a pivotal role in the formation of phenotypic features of defeated mice.
According to the results presented in Figure 2, the proteins encoded by the transcription factor genes Otx2 and Six3, the expressions of which are changed in opposite directions in winners and defeated mice compared to controls, can also be considered as key (nodal) genes. The investigations of the mechanisms involved in neural functions have shown that the expression of several regulatory genes, including Six3, cannot be initiated in the neural plate of Otx2−/− embryos [42,43], suggesting that the expressions of the Otx2 and Six3 genes are closely related.
Recently the Otx2 gene was characterized as an upstream mediator of increased susceptibility to social defeat stress in mice. Its transient knockdown in the VTA increased stress susceptibility, and its overexpression was associated with reverse effects [44]. The role for the OTX2 gene was confirmed in the study of stress-related depressive disorder pathophysiology in children [45]. Otx2 is involved in the regulation of the number of VTA neurons with efficient dopamine uptake [46] and is considered one of the genes with a key role in midbrain dopaminergic neuron development and the regulation of their survival and physiology [47].
To date, the Six3 gene has not been associated with the behavior/neurological phenotype, but its role in the control of the status of the neural progenitor cells has been reported. It was shown that the upregulation of Six3 expression kept the progenitor cells of the embryonic telencephalon in an undifferentiated state. The authors suggested that Six3 is involved in the control of the subtle equilibrium between the proliferation and differentiation of neural progenitor cells during mammalian neurogenesis [48]. A conditional deletion of the Six3 gene prevented the formation of most dopamine receptor DRD2-expressing medium spiny neurons being the principal projection neurons in the striatum [49]. To date, the role of the Six3 gene in the functioning of the VTA neurons has not been reported. However, based on the above-cited data on its role in the formation of the DRD2-expressing neurons of the striatum, it can be assumed that the Six3 gene can also participate in the regulation of D2R-expressing neurons functioning in the VTA. Our results from the cluster and correlation analyses confirm the plausibility of this hypothesis. Accordingly, a further study of the effects of Six3 gene expression in the VTA of the winners and defeated animals on dopaminergic signaling may be promising for understanding the formation of their behavioral characteristics.
The VTA is one of the major sites of dopamine synthesis. The VTA dopaminergic neurons are implicated in the regulation of emotional and motivational behaviors and are involved in the development of psychopathologies including depression [50], aggression [51], and addictions [52,53]. So, it was important to check if the prioritized genes could be implicated in VTA dopaminergic signaling. The results of the cluster analysis on expression data of the genes differentially expressed in the VTA of winners and defeated mice showed that the Nrgn, Ercc2, Otx2, Six3, and Camk2a genes were assigned to the same cluster as the dopaminergic (Th, Ddc, and Slc6a3) genes. According to the correlation analysis, the expression of Otx2, Six3, and Camk2a genes is coregulated with the dopaminergic genes in both experimental groups. The expression of two other genes (Nrgn and Ercc2) correlated with the expression of dopaminergic genes only in the VTA of losers indicating a possible coregulation or functional relationship of these genes in the VTA of defeated mice. It should be emphasized that the most significant correlations were found between the expression of dopaminergic genes and the expression of the Nrgn gene (see Table 8).
The development of technologies generating large amounts of gene expression data called for the use of mathematical methods capable of grouping the genes based on the similarity of their expression patterns. Such methods include a correlation analysis [54,55], which is a potent tool to identify coregulated genes or groups of genes associated with particular experimental or environmental conditions [56].
To date, the role of the Nrgn gene in dopaminergic signaling remains unclear, however, there is indirect evidence for its involvement [57]. Additionally, in favor of this hypothesis are the findings that neurogranin is enriched at dendritic spines [27] whereas D2R is localized to the VTA dendrites [58], which suggests their spatial colocalization.
According to our results, Nrgn is involved in the network of 31 coregulated genes, encoding proteins associated with the nervous system development, including 24 genes encoding proteins associated with the generation of neurons. Several of them (Otx2, En1) are known as regulators of terminal differentiation, survival and the maintenance of midbrain dopaminergic neurons [59]. We believe that other coregulated genes included in the network may also be essential for studying their role in VTA dopaminergic neuron maturation under the influence of social stress.
Of course, many genes with VTA expression levels significantly changed in one of experimental groups and not changed in the other group could be directly or indirectly related to the formation of behavioral characteristics in winners and defeated mice. However, we believe that the genes that change the level of transcription in opposite directions in winners and defeated mice deserve special research attention. As noted earlier, the Nrgn, Ercc2, Otx2, and Six3 genes that change their transcription levels in opposite directions in winners and defeated mice and make the maximum contribution to intergroup differences can be considered as key genes in the formation of alternative phenotypes associated with behavioral patterns in winning and defeated animals.

4. Materials and Methods

4.1. Animals

The work was carried out using 10–12-week old C57BL/6J male mice obtained from the Animal Breeding Facility, Branch of Institute of Bioorganic Chemistry of the Russian Academy of Sciences (Pushchino, Moscow region, Russia). Mice were kept under standard conditions at 22+/−2 °C with a 12/12 h light–dark cycle (light period started from 8:00 AM) and dry laboratory food and water given ad libitum. All procedures were in compliance with the European Communities Council Directive 210/63/EU on 22 September 2010. The study was approved by Scientific Council N 9 of the Institute of Cytology and Genetics SB RAS of 24 March 2010, N 613 (Novosibirsk, Russia).

4.2. Generation of Alternative Forms of Social Behavior in Male Mice under Agonistic Interactions

Alternative forms of social behavior were generated as a result of the daily agonistic interactions (intermale confrontations) of male mice as described previously [1,3]. Pairs of weight-matched mice were each placed in a cage (14 × 28 × 10 cm) bisected by a perforated transparent partition allowing the animals to see, hear and smell each other, but preventing physical contact. The animals were left undisturbed for two or three days to adapt to unfamiliar housing conditions before they were exposed to encounters. Every day at 14:00–17:00 p.m. (Russian local time), the cage lid was replaced by a transparent one, and after 5 min (the period necessary to activate agonistic interactions), the partition was removed for 10 min for intermale confrontation. The mouse that attacked, bit, and chased the opponent was considered the winner. The superiority of the winner was established by the results of two or three encounters with the same opponent. A mouse that showed only defensive behavior (sideways postures, upright postures, withdrawal, lying on the back or freezing) was defined as the loser. If the aggressive attacks were very active and long, then the interactions between the males were stopped after 3 min (or even earlier) by restoring the partition, to prevent damage to the defeated male. That means that the painful effects of agonistic interactions are absent in this model for defeated mice.
Each defeated mouse (loser) was exposed to the same winner for three days, and then, to continue the agonistic interactions, the defeated mouse was placed in an unfamiliar cage with an unfamiliar winner behind the partition. Each winner remained in its original cage. The intermale confrontations procedure was performed once a day for 21 days and yielded an equal number of winners and losers.
Three groups of animals were used: (1) controls—mice without a consecutive experience of agonistic interactions; (2) winners—group of aggressive mice chronically winning during 21 days in daily agonistic interactions (intermale confrontations); (3) losers (defeated mice)—mice with chronic experience of defeats during 21 days in daily agonistic interactions. Animals with alternative types of social experience demonstrated the development of various pathological behaviors. Mice with a long experience of aggression and victories (winners) showed increased aggressiveness, hyperactivity, stereotyped behaviors, anxiety, impaired social recognition, irritability, autistic spectrum symptoms, a condition similar to drug addiction, etc. (reviewed in [60]). The losers developed mixed anxiety/depression-like behaviors accompanied by full immobility, strong anxiety, avoidance of any social interactions, helplessness, indifference, etc. [6,61].
A group of 21 pairs of male C57BL/6J mice was used to generate alternative forms of social behavior. Winners and losers with the most expressed behavioral phenotypes were selected for the transcriptome analysis. The following criteria were used. For the losers, during the activation period (5 min before a fight) the chronically defeated mice demonstrated all symptoms of depressive behavior: they did not approach the partition, sit in the corner of the cage opposite to the partition, or in a corner or into litter with their nose; they were characterized by immobility, freezing under winners’ attacking, and a demonstrated indifference in all experimental situations (without behavioral reactions); there were no inversions of behavior to the opposite one after a change of aggressors, and they showed avoidance and passive defense when attacked by the aggressor. Aggressive winners, during the activation period demonstrated a strong aggressive motivation, and every day immediately attacked the losers after partition removal, stopping only for rest, displaying a manic motivation to bite losers in spite of full submission.
Control animals and all experimental mice were simultaneously decapitated. Experimental mice were decapitated 24 h after the last agonistic interaction. The brain regions were dissected by the same experimenter according to the map presented in the Allen Mouse Brain Atlas [62]. All biological samples were placed in RNAlater solution (Thermo Fisher Scientific, Waltham, MA, USA) and stored at −70 °C until sequencing.

4.3. RNA-Seq Analysis

The frozen VTA samples from male mice: winners (n = 3), losers (n = 3) and controls (n = 3), were sent to JSC Genoanalytica (Moscow, Russia) which specializes in RNA-Seq analysis. For transcriptome profiling, mRNA was extracted using the Dynabeads mRNA Purification Kit (Ambion, Thermo Fisher Scientific, Waltham, MA, USA). cDNA libraries were constructed using the NEBNext mRNA Library Prep Reagent Set for Illumina (NEB, Ipswich, MA USA). All experimental procedures were performed in accordance with the manufacturer’s protocols. A single-end sequencing of cDNA libraries was performed on the Illumina Hiseq 1500 platform (Illumina Sequencing, San Diego, CA, USA) with a read length of 50 bases. All samples were analyzed as biological replicates. The quality metrics of the mapped data (Supplementary Materials Table S8) were collected using Spliced Transcripts Alignment to a Reference (STAR) software (Cold Spring Harbor, NY, USA) [63]. The sequencing data were mapped to the mouse reference genome sequence (GRCm38.p3) available in GenBank using the TopHat2 aligner (Center for Bioinformatics and Computational Biology, University of Maryland, College Park, MD, USA) [64]. The RNA-Seq data sets are available at European Nucleotide Archive (Accession PRJEB39532)
Cufflinks/Cuffdiff programs were run to estimate the gene expression levels in FPKM (fragments per kilobase of transcript per million mapped reads) and to identify the differentially expressed genes (DEGs) in the experimental and control groups. Genes were considered differentially expressed at a false discovery rate (q value) <5% [65]. Only annotated gene sequences were used in the analysis.
DEGs obtained by comparing winners with losers were considered as bidirectionally changing the level of expression if their transcription level changed in different directions by at least 1.5 times relative to the control in the control/winner and control/loser comparisons.

4.4. Functional Annotation

The functional annotation of DEGs was performed using DAVID (The Database for Annotation, Visualization and Integrated Discovery) gene annotation tool [66]. The Mus musculus genome was utilized as the background list for the over-representation analysis. The Gene Ontology (GO) option in DAVID and Kyoto Encyclopedia of Genes and Genomes (KEGG) Pathway Database [67] were used to identify the significantly (p < 0.05) enriched biological processes and metabolic pathways. To construct the functional protein association network the STRING database was used [68].
To define the association of DEGs with behavior/neurological phenotype the Neurological Disease Portal (Phenotypes, Mouse) in Rat Genome Database (RGD [69] was employed.
An atlas of combinatorial transcriptional regulation in mice and men [70] was used to reveal the DEGs encoding the transcription factor genes.

4.5. Statistical Methods

The acquired RNA-Seq data (in FPKM values) were log transformed, centered, and normalized. Then, the data sets were scaled by the principal coordinates (PCOs) method based on Euclidean metric distances and the partial-least squares discriminant analysis (PLS-DA) was conducted using the pattern of covariation for linear combinations between two blocks of variables [71]. The constructed PLS-DA Axis maximized the distance between the analyzed groups of animals. After that, the Pearson correlation between gene expression levels and the coordinates of animals along the first functionally meaningful synthetic axis (PLS-DA Axis 1) was employed to define a set of expressed genes that are characterized by the highest correlation coefficients. Such genes are thought to contribute the most to intergroup differences. Normalized RNA-Seq data were also used for correlation analysis. The software packages Statistica 6.0 (StatSoft, Tulsa, OK, USA) and JACOBI4 [72] were used for the multivariate analysis and representation of the data.
Agglomerative Hierarchical Clustering (AHC) analysis was carried out by XLStat software -(Version 2016.02, Addinsoft Inc., Paris France). We used the unweighted pair-group average agglomerative method, noncentered mode, with the Pearson pairwise correlation matrix as an input.

5. Conclusions

The analysis of gene expression profiles in the VTA of mice with alternative social behaviors showed that during daily social confrontations, the transcription of many genes in both chronically defeated and aggressive mice changes unidirectionally, which indicates that mice in both experimental groups are under stress. However, the results unambiguously indicated that the consequences of long-term social stress of agonistic interactions are much more pronounced in defeated mice, which develop a mixed anxiety/depression-like state under conditions of chronic unavoidable psychoemotional social stress.
The comprehensive bioinformatics analysis carried out in this work made it possible to identify the potentially important candidate genes (Nrgn, Ercc2, Otx2, and Six3) that may be involved in the formation of the behavioral differences in losers and winners. It was first shown that Nrgn (neurogranin) expression was highly correlated with the expression of the genes involved in dopamine synthesis and transport (Th, Ddc, Slc6a3, and Drd2) in the VTA of defeated mice but not in winners. This study may be useful for understanding specific and nonspecific molecular events occurring in VTA neurons during the formation of human disorders associated with aggressive or depressive-like states.
The limitations of this study include the following:
  • It is necessary to remember the complexity of the mechanisms regulating behavior. Accordingly, the results obtained in this study may be valid for mixed anxiety/depressive disorder and pathological aggressive behavior in mice within the framework of the used experimental model. However, we tried to describe the methodology of the experiment in detail in order to increase the replicability of the study.
  • The sequencing and bioinformatic analysis carried out in this work is only the beginning of the process of studying the role of the identified candidate genes (Nrgn and other highlighted genes) in the formation of behavioral features under the influence of social confrontations. To establish the relationship between the transcription of the candidate genes and the external (behavioral) phenotype, additional studies at the proteome and metabolome levels are required.

Supplementary Materials

The following are available online at https://www.mdpi.com/1422-0067/21/18/6599/s1, Supplementary Table S1. Genes differentially expressed in VTA of winners (aggressive) mice as compared to controls; Supplementary Table S2. Genes differentially expressed in VTA of losers as compared to controls; Supplementary Table S3. Metabolic pathways enriched by genes differentially expressed in VTA of losers as compared to control mice; Supplementary Table S4. Genes differentially expressed in VTA of winners (aggressive mice) as compared to losers; Supplementary Table S5. Genes making the maximum contribution to intergroup (winners vs. losers) differences (correlation more than 0.90 or −0.90); Supplementary Table S6. DEGs characterized by correlation more than 0.90 or −0.90; Supplementary File S7. Agglomerative hierarchical clustering (AHC) of VTA DEGs in winners vs. losers comparison; Supplementary Table S8. The summary statistics for the sequenced libraries.

Author Contributions

Conceptualization, N.K.; Data curation, V.B., D.S., I.K. and A.G.; Formal analysis, O.R., V.B., D.S., I.K., A.G. and V.E.; Funding acquisition, N.K.; Investigation, O.R., V.B., D.S., I.K. and A.G.; Methodology, V.E. and N.K.; Software, V.E.; Supervision, N.K.; Writing—original draft, O.R.; Writing—review and editing, O.R., V.B., D.S., I.K., A.G., V.E. and N.K. All authors have read and agreed to the published version of the manuscript.

Funding

The current study was supported by Russian Science Foundation (grant No 19-15-00026). Preparation and maintenance of experimental animals was carried out in the Conventional vivarium of the Institute of Cytology and Genetics SB RAS and was funded under Budjet project No. 0324-2019-0041-C-01.

Acknowledgments

The authors are grateful to JSC Genoanalytica (Moscow, Russia) for conducting the technological part of the experiment and the primary statistical analysis.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Kudryavtseva, N.N. The sensory contact model for the study of aggressive and submissive behaviors in male mice. Aggress. Behav. 1991, 17, 285–291. [Google Scholar] [CrossRef]
  2. Golden, S.A.; Covington, H.E., III; Berton, O.; Russo, S.J. A standardized protocol for repeated social defeat stress in mice. Nat. Protoc. 2011, 6, 1183–1191. [Google Scholar] [CrossRef]
  3. Kudryavtseva, N.N.; Smagin, D.A.; Kovalenko, I.L.; Vishnivetskaya, G.B. Repeated positive fighting experience in male inbred mice. Nat. Protoc. 2014, 9, 2705–2717. [Google Scholar] [CrossRef]
  4. Shimamoto, A. Social Defeat Stress, Sex, and Addiction-Like Behaviors. Int. Rev. Neurobiol. 2018, 140, 271–313. [Google Scholar] [CrossRef]
  5. Smagin, D.A.; Kovalenko, I.L.; Galyamina, A.G.; Bragin, A.O.; Orlov, Y.L.; Kudryavtseva, N.N. Dysfunction in Ribosomal Gene Expression in the Hypothalamus and Hippocampus following Chronic Social Defeat Stress in Male Mice as Revealed by RNA-Seq. Neural Plast. 2016, 2016, 3289187. [Google Scholar] [CrossRef] [Green Version]
  6. Galyamina, A.G.; Kovalenko, I.L.; Smagin, D.A.; Kudryavtseva, N.N. Interaction of depression and anxiety in the development of mixed anxiety/depression disorder: Experimental studies of the mechanisms of comorbidity (review). Neurosci. Behav. Physiol. 2017, 47, 699–713. [Google Scholar] [CrossRef]
  7. Kovalenko, I.L.; Smagin, D.A.; Galyamina, A.G.; Orlov, Y.L.; Kudryavtseva, N.N. Changes in the Expression of Dopaminergic Genes in Brain Structures of Male Mice Exposed to Chronic Social Defeat Stress: An RNA-seq Study. Mol. Biol. 2016, 50, 184–187. [Google Scholar] [CrossRef]
  8. Kudryavtseva, N.N.; Smagin, D.A.; Kovalenko, I.L.; Galyamina, A.G.; Vishnivetskaya, G.B.; Babenko, V.N.; Orlov, Y.L. Serotonergic genes in the development of anxiety/depression-like state and pathology of aggressive behavior in male mice: RNA-seq data. Mol. Biol. 2017, 51, 288–300. [Google Scholar] [CrossRef]
  9. Galyamina, A.G.; Kovalenko, I.L.; Smagin, D.A.; Kudryavtseva, N.N. Altered Expression of Neurotransmitters Systems’ Genes in the Ventral Tegmental Area of Depressive Male Mice: Data of RNA-Seq. Zhurnal Vysshey Nervnoy Deyatel’nosti IP Pavlova 2017, 67, 113–128. [Google Scholar]
  10. Smagin, D.A.; Kovalenko, I.L.; Galyamina, A.G.; Orlov, Y.L.; Babenko, V.N.; Kudryavtseva, N.N. Heterogeneity of Brain Ribosomal Genes Expression Following Positive Fighting Experience in Male Mice as Revealed by RNA-Seq. Mol. Neurobiol. 2018, 55, 390–401. [Google Scholar] [CrossRef]
  11. Babenko, V.N.; Smagin, D.A.; Kudryavtseva, N.N. RNA-Seq Mouse Brain Regions Expression Data Analysis: Focus on ApoE Functional Network. J. Integr. Bioinform. 2017, 14. [Google Scholar] [CrossRef]
  12. Babenko, V.N.; Smagin, D.A.; Galyamina, A.G.; Kovalenko, I.L.; Kudryavtseva, N.N. Altered Slc25 family gene expression as markers of mitochondrial dysfunction in brain regions under experimental mixed anxiety/depression-like disorder. BMC Neurosci. 2018, 19, 79. [Google Scholar] [CrossRef]
  13. Smagin, D.A.; Galyamina, A.G.; Kovalenko, I.L.; Babenko, V.N.; Kudryavtseva, N.N. Aberrant Expression of Collagen Gene Family in the Brain Regions of Male Mice with Behavioral Psychopathologies Induced by Chronic Agonistic Interactions. BioMed Res. Int. 2019, 2019, 7276389. [Google Scholar] [CrossRef] [Green Version]
  14. Morales, M.; Margolis, E.B. Ventral tegmental area: Cellular heterogeneity, connectivity and behaviour. Nat. Rev. Neurosci. 2017, 18, 73–85. [Google Scholar] [CrossRef]
  15. Han, X.; Jing, M.Y.; Zhao, T.Y.; Wu, N.; Song, R.; Li, J. Role of dopamine projections from ventral tegmental area to nucleus accumbens and medial prefrontal cortex in reinforcement behaviors assessed using optogenetic manipulation. Metab. Brain Dis. 2017, 32, 1491–1502. [Google Scholar] [CrossRef]
  16. Holly, E.N.; Miczek, K.A. Ventral tegmental area dopamine revisited: Effects of acute and repeated stress. Psychopharmacology (Berl.) 2016, 233, 163–186. [Google Scholar] [CrossRef] [Green Version]
  17. Barker, D.J.; Root, D.H.; Zhang, S.; Morales, M. Multiplexed neurochemical signaling by neurons of the ventral tegmental area. J. Chem. Neuroanat. 2016, 73, 33–42. [Google Scholar] [CrossRef] [Green Version]
  18. Glick, A.R. The role of serotonin in impulsive aggression, suicide, and homicide in adolescents and adults: A literature review. Int. J. Adolesc. Med. Health 2015, 27, 143–150. [Google Scholar] [CrossRef]
  19. Dell’Osso, L.; Carmassi, C.; Mucci, F.; Marazziti, D. Depression, Serotonin and Tryptophan. Curr. Pharm. Des. 2016, 22, 949–954. [Google Scholar] [CrossRef]
  20. Kraus, C.; Castren, E.; Kasper, S.; Lanzenberger, R. Serotonin and neuroplasticity-Links between molecular, functional and structural pathophysiology in depression. Neurosci. Biobehav. Rev. 2017, 77, 317–326. [Google Scholar] [CrossRef] [Green Version]
  21. Kudryavtseva, N.N. Serotonergic Control of Aggressive Behavior: Novel Approaches--New Interpretations (Review). Zhurnal Vysshey Nervnoy Deyatel’nosti IP Pavlova 2015, 65, 546–563. [Google Scholar]
  22. Coccaro, E.F.; Fanning, J.R.; Phan, K.L.; Lee, R. Serotonin and impulsive aggression. CNS Spectr. 2015, 20, 295–302. [Google Scholar] [CrossRef]
  23. Bondar, N.P.; Boyarskikh, U.A.; Kovalenko, I.L.; Filipenko, M.L.; Kudryavtseva, N.N. Molecular implications of repeated aggression: Th, Dat1, Snca and Bdnf gene expression in the VTA of victorious male mice. PLoS ONE 2009, 4, e4190. [Google Scholar] [CrossRef]
  24. Avgustinovich, D.F.; Gorbach, O.V.; Kudryavtseva, N.N. Comparative analysis of anxiety-like behavior in partition and plus-maze tests after agonistic interactions in mice. Physiol. Behav. 1997, 61, 37–43. [Google Scholar] [CrossRef]
  25. Kudryavtseva, N.N.; Bondar, N.P.; Avgustinovich, D.F. Association between experience of aggression and anxiety in male mice. Behav. Brain Res. 2002, 133, 83–93. [Google Scholar] [CrossRef]
  26. Smagin, D.A.; Park, J.H.; Michurina, T.V.; Peunova, N.; Glass, Z.; Sayed, K.; Bondar, N.P.; Kovalenko, I.N.; Kudryavtseva, N.N.; Enikolopov, G. Altered Hippocampal Neurogenesis and Amygdalar Neuronal Activity in Adult Mice with Repeated Experience of Aggression. Front. Neurosci. 2015, 9, 443. [Google Scholar] [CrossRef]
  27. Petersen, A.; Gerges, N.Z. Neurogranin regulates CaM dynamics at dendritic spines. Sci. Rep. 2015, 5, 11135. [Google Scholar] [CrossRef] [Green Version]
  28. Gerendasy, D.D.; Sutcliffe, J.G. RC3/neurogranin, a postsynaptic calpacitin for setting the response threshold to calcium influxes. Mol. Neurobiol. 1997, 15, 131–163. [Google Scholar] [CrossRef]
  29. Li, H.Y.; Li, J.F.; Lu, G.W. Neurogranin: A brain-specific protein. Sheng Li Ke Xue Jin Zhan 2003, 34, 111–115. [Google Scholar]
  30. Pak, J.H.; Huang, F.L.; Li, J.; Balschun, D.; Reymann, K.G.; Chiang, C.; Westphal, H.; Huang, K.P. Involvement of neurogranin in the modulation of calcium/calmodulin-dependent protein kinase II, synaptic plasticity, and spatial learning: A study with knockout mice. Proc. Natl. Acad. Sci. USA 2000, 97, 11232–11237. [Google Scholar] [CrossRef] [Green Version]
  31. Diez-Guerra, F.J. Neurogranin, a link between calcium/calmodulin and protein kinase C signaling in synaptic plasticity. IUBMB Life 2010, 62, 597–606. [Google Scholar] [CrossRef]
  32. Huang, F.L.; Huang, K.P. Methylphenidate improves the behavioral and cognitive deficits of neurogranin knockout mice. Genes Brain Behav. 2012, 11, 794–805. [Google Scholar] [CrossRef] [Green Version]
  33. Huang, K.P.; Huang, F.L.; Jager, T.; Li, J.; Reymann, K.G.; Balschun, D. Neurogranin/RC3 enhances long-term potentiation and learning by promoting calcium-mediated signaling. J. Neurosci. 2004, 24, 10660–10669. [Google Scholar] [CrossRef] [Green Version]
  34. Kaleka, K.S.; Gerges, N.Z. Neurogranin restores amyloid beta-mediated synaptic transmission and long-term potentiation deficits. Exp. Neurol. 2016, 277, 115–123. [Google Scholar] [CrossRef]
  35. Braun, A.P.; Schulman, H. The multifunctional calcium/calmodulin-dependent protein kinase: From form to function. Annu. Rev. Physiol. 1995, 57, 417–445. [Google Scholar] [CrossRef]
  36. Hudmon, A.; Schulman, H. Structure-function of the multifunctional Ca2+/calmodulin-dependent protein kinase II. Biochem. J. 2002, 364, 593–611. [Google Scholar] [CrossRef]
  37. Yamauchi, T. Neuronal Ca2+/calmodulin-dependent protein kinase II--discovery, progress in a quarter of a century, and perspective: Implication for learning and memory. Biol. Pharm. Bull. 2005, 28, 1342–1354. [Google Scholar] [CrossRef] [Green Version]
  38. Robison, A.J. Emerging role of CaMKII in neuropsychiatric disease. Trends Neurosci. 2014, 37, 653–662. [Google Scholar] [CrossRef]
  39. Zalcman, G.; Federman, N.; Romano, A. CaMKII Isoforms in Learning and Memory: Localization and Function. Front. Mol. Neurosci. 2018, 11, 445. [Google Scholar] [CrossRef]
  40. Glaser, T.; Arnaud Sampaio, V.F.; Lameu, C.; Ulrich, H. Calcium signalling: A common target in neurological disorders and neurogenesis. Semin. Cell Dev. Biol. 2019, 95, 25–33. [Google Scholar] [CrossRef]
  41. Lerner, L.K.; Moreno, N.C.; Rocha, C.R.R.; Munford, V.; Santos, V.; Soltys, D.T.; Garcia, C.C.M.; Sarasin, A.; Menck, C.F.M. XPD/ERCC2 mutations interfere in cellular responses to oxidative stress. Mutagenesis 2019, 34, 341–354. [Google Scholar] [CrossRef] [PubMed]
  42. Rhinn, M.; Dierich, A.; Shawlot, W.; Behringer, R.R.; Le Meur, M.; Ang, S.L. Sequential roles for Otx2 in visceral endoderm and neuroectoderm for forebrain and midbrain induction and specification. Development 1998, 125, 845–856. [Google Scholar] [PubMed]
  43. Acampora, D.; Avantaggiato, V.; Tuorto, F.; Briata, P.; Corte, G.; Simeone, A. Visceral endoderm-restricted translation of Otx1 mediates recovery of Otx2 requirements for specification of anterior neural plate and normal gastrulation. Development 1998, 125, 5091–5104. [Google Scholar]
  44. Pena, C.J.; Kronman, H.G.; Walker, D.M.; Cates, H.M.; Bagot, R.C.; Purushothaman, I.; Issler, O.; Loh, Y.E.; Leong, T.; Kiraly, D.D.; et al. Early life stress confers lifelong stress susceptibility in mice via ventral tegmental area OTX2. Science 2017, 356, 1185–1188. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Kaufman, J.; Wymbs, N.F.; Montalvo-Ortiz, J.L.; Orr, C.; Albaugh, M.D.; Althoff, R.; O’Loughlin, K.; Holbrook, H.; Garavan, H.; Kearney, C.; et al. Methylation in OTX2 and related genes, maltreatment, and depression in children. Neuropsychopharmacology 2018, 43, 2204–2211. [Google Scholar] [CrossRef] [Green Version]
  46. Di Salvio, M.; Di Giovannantonio, L.G.; Acampora, D.; Prosperi, R.; Omodei, D.; Prakash, N.; Wurst, W.; Simeone, A. Otx2 controls neuron subtype identity in ventral tegmental area and antagonizes vulnerability to MPTP. Nat. Neurosci. 2010, 13, 1481–1488. [Google Scholar] [CrossRef] [Green Version]
  47. Blaudin de Thé, F.-X.; Rekaik, H.; Prochiantz, A.; Fuchs, J.; Joshi, R.L. Neuroprotective Transcription Factors in Animal Models of Parkinson Disease. Neural Plast. 2016, 2016, 6097107. [Google Scholar] [CrossRef] [Green Version]
  48. Appolloni, I.; Calzolari, F.; Corte, G.; Perris, R.; Malatesta, P. Six3 controls the neural progenitor status in the murine CNS. Cereb. Cortex 2008, 18, 553–562. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  49. Xu, Z.; Liang, Q.; Song, X.; Zhang, Z.; Lindtner, S.; Li, Z.; Wen, Y.; Liu, G.; Guo, T.; Qi, D.; et al. SP8 and SP9 coordinately promote D2-type medium spiny neuron production by activating Six3 expression. Development 2018, 145. [Google Scholar] [CrossRef] [Green Version]
  50. Kaufling, J. Alterations and adaptation of ventral tegmental area dopaminergic neurons in animal models of depression. Cell Tissue Res. 2019, 377, 59–71. [Google Scholar] [CrossRef]
  51. Tadano, T.; Abe, Y.; Morikawa, Y.; Asao, T.; Hozumi, M.; Takahashi, N.; Tan-no, K.; Kisara, K. Involvement of dopaminergic neurons in mouse-killing aggression in rats. Methods Find. Exp. Clin. Pharmacol. 1997, 19, 527–531. [Google Scholar] [PubMed]
  52. Solecki, W.; Wilczkowski, M.; Pradel, K.; Karwowska, K.; Kielbinski, M.; Drwiega, G.; Zajda, K.; Blasiak, T.; Soltys, Z.; Rajfur, Z.; et al. Effects of brief inhibition of the ventral tegmental area dopamine neurons on the cocaine seeking during abstinence. Addict. Biol. 2019, e12826. [Google Scholar] [CrossRef] [PubMed]
  53. Simmons, S.C.; Wheeler, K.; Mazei-Robison, M.S. Determination of circuit-specific morphological adaptations in ventral tegmental area dopamine neurons by chronic morphine. Mol. Brain 2019, 12, 10. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Ma, Y.; Luo, T.; Dong, D.; Wu, X.; Wang, Y. Characterization of long non-coding RNAs to reveal potential prognostic biomarkers in hepatocellular carcinoma. Gene 2018, 663, 148–156. [Google Scholar] [CrossRef]
  55. Ma, L.; Deng, C. Identification of a novel four-lncRNA signature as a prognostic indicator in cirrhotic hepatocellular carcinoma. PeerJ 2019, 7, e7413. [Google Scholar] [CrossRef] [Green Version]
  56. Grant, G.R.; Manduchi, E.; Stoeckert, C.J., Jr. Analysis and Management of Microarray Gene Expression Data. Curr. Protoc. Mol. Biol. 2007, 77. [Google Scholar] [CrossRef]
  57. Caetano-Anolles, K.; Rhodes, J.S.; Garland, T., Jr.; Perez, S.D.; Hernandez, A.G.; Southey, B.R.; Rodriguez-Zas, S.L. Cerebellum Transcriptome of Mice Bred for High Voluntary Activity Offers Insights into Locomotor Control and Reward-Dependent Behaviors. PLoS ONE 2016, 11, e0167095. [Google Scholar] [CrossRef] [Green Version]
  58. Garzon, M.; Duffy, A.M.; Chan, J.; Lynch, M.K.; Mackie, K.; Pickel, V.M. Dopamine D(2) and acetylcholine alpha7 nicotinic receptors have subcellular distributions favoring mediation of convergent signaling in the mouse ventral tegmental area. Neuroscience 2013, 252, 126–143. [Google Scholar] [CrossRef] [Green Version]
  59. Arenas, E.; Denham, M.; Villaescusa, J.C. How to make a midbrain dopaminergic neuron. Development 2015, 142, 1918–1936. [Google Scholar] [CrossRef] [Green Version]
  60. Kudryavtseva, N.N. The psychopathology of repeated aggression: A neurobiological aspect. In Perspectives on the Psychology of Aggression; Morgan, J.P., Ed.; Nova Science Publishers Inc.: New York, NY, USA, 2006; pp. 35–64. [Google Scholar]
  61. Kudryavtseva, N.N.; Bakshtanovskaya, I.V.; Koryakina, L.A. Social model of depression in mice of C57BL/6J strain. Pharmacol. Biochem. Behav. 1991, 38, 315–320. [Google Scholar] [CrossRef]
  62. Allen Institute for Brain Science. Allen Mouse Brain Atlas 2004. Available online: http://mouse.brain-map.org/static/atlas (accessed on 9 September 2020).
  63. Dobin, A.; Davis, C.A.; Schlesinger, F.; Drenkow, J.; Zaleski, C.; Jha, S.; Batut, P.; Chaisson, M.; Gingeras, T.R. STAR: Ultrafast universal RNA-seq aligner. Bioinformatics 2013, 29, 15–21. [Google Scholar] [CrossRef] [PubMed]
  64. Kim, D.; Pertea, G.; Trapnell, C.; Pimentel, H.; Kelley, R.; Salzberg, S.L. TopHat2: Accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions. Genome Biol. 2013, 14, R36. [Google Scholar] [CrossRef] [Green Version]
  65. Trapnell, C.; Hendrickson, D.G.; Sauvageau, M.; Goff, L.; Rinn, J.L.; Pachter, L. Differential analysis of gene regulation at transcript resolution with RNA-seq. Nat. Biotechnol. 2013, 31, 46–53. [Google Scholar] [CrossRef] [PubMed]
  66. Huang, D.W.; Sherman, B.T.; Lempicki, R.A. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat. Protoc. 2009, 4, 44–57. [Google Scholar] [CrossRef]
  67. Kanehisa, M. Post-genome Informatics; Kyoto Encyclopedia of Genes and Genomes; Oxford University Press: Oxford, UK, 2000; 148p, Available online: http://www.genome.jp/kegg (accessed on 26 March 2020).
  68. Szklarczyk, D.; Gable, A.L.; Lyon, D.; Junge, A.; Wyder, S.; Huerta-Cepas, J.; Simonovic, M.; Doncheva, N.T.; Morris, J.H.; Bork, P.; et al. STRING v11: Protein-protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets. Nucleic Acids Res. 2019, 47, D607–D613. [Google Scholar] [CrossRef] [Green Version]
  69. Neurological Disease Portal, Rat Genome Database Web Site, Medical College of Wisconsin, Milwaukee, Wisconsin, USA. Available online: http://rgd.mcw.edu/ (accessed on 18 October 2019).
  70. Ravasi, T.; Suzuki, H.; Cannistraci, C.V.; Katayama, S.; Bajic, V.B.; Tan, K.; Akalin, A.; Schmeier, S.; Kanamori-Katayama, M.; Bertin, N.; et al. An atlas of combinatorial transcriptional regulation in mouse and man. Cell 2010, 140, 744–752. [Google Scholar] [CrossRef] [Green Version]
  71. Barker, M.; Rayens, W. Partial least squares for discrimination. J. Chemom. 2003, 17, 166–173. [Google Scholar] [CrossRef]
  72. Polunin, D.; Shtaiger, I.; Efimov, V. JACOBI4 software for multivariate analysis of biological data. bioRxiv 2019, 803684. [Google Scholar] [CrossRef] [Green Version]
Ijms 21 06599 g001 550
Figure 1. (A). Axis maximizing the distances between winning and defeated mice; (B). The distribution of expressed genes along the axis representing the correlation between gene expression and mice coordinates along the partial-least squares discriminant analysis (PLS-DA) Axis 1.
Figure 1. (A). Axis maximizing the distances between winning and defeated mice; (B). The distribution of expressed genes along the axis representing the correlation between gene expression and mice coordinates along the partial-least squares discriminant analysis (PLS-DA) Axis 1.
Ijms 21 06599 g001
Ijms 21 06599 g002 550
Figure 2. The relationships between DEGs that make the most significant contribution to the intergroup (winners vs. losers) differences according to the PLS-DA analysis. STRING database utilities were used to identify possible associations between encoded proteins. The genes encoding proteins associated with the nervous system development are shown in red; the genes associated with generation of neurons are shown in blue.
Figure 2. The relationships between DEGs that make the most significant contribution to the intergroup (winners vs. losers) differences according to the PLS-DA analysis. STRING database utilities were used to identify possible associations between encoded proteins. The genes encoding proteins associated with the nervous system development are shown in red; the genes associated with generation of neurons are shown in blue.
Ijms 21 06599 g002
Table
Table 1. Differentially Expressed Genes (DEGs) associated with the behavior/neurological phenotype.
Table 1. Differentially Expressed Genes (DEGs) associated with the behavior/neurological phenotype.
Gene SymbolLocusExpression in VTA (FPKM)Log2 (Fold Change) Winners/Controlsq_ValueDefinition
Control MiceWinners
Winners vs. Control Mice
Dnaja1chr4:40720153-4075788541.90104.771.322.05 × 10−2DnaJ heat shock protein family (Hsp40) member A1
Ercc2#chr7:19382024-1940410453.428.72−2.622.05 × 10−2excision repair cross-complementing rodent repair deficiency, complementation group 2
Lbx1chr19:45232619-452353777.183.19−1.172.05 × 10−2ladybird homeobox homolog 1 (Drosophila)
Mafchr8:115682935-1157079409.665.98−0.693.58 × 10−2avian musculoaponeurotic fibrosarcoma oncogene homolog
Myh6chr14:54936101-549666071.4759.905.352.05 × 10−2myosin, heavy polypeptide 6, cardiac muscle, alpha
Ngfrchr11:95568820-955876983.552.01−0.822.05 × 10−2nerve growth factor receptor (TNFR superfamily, member 16)
Nkx6-1chr5:101658055-1016653615.012.79−0.853.58 × 10−2NK6 homeobox 1
Nptx2chr5:144545886-1445574783.281.59−1.042.05 × 10−2neuronal pentraxin 2
Nrgnchr9:37544492-3755274522.3514.15−0.663.58 × 10−2neurogranin
Slc17a7#,§,©chr7:45163920-4517632810.7420.590.942.05 × 10−2solute carrier family 17 (sodium-dependent inorganic phosphate cotransporter), member 7
Slc6a4#,§,©chr11:76998596-770323435.039.990.992.05 × 10−2solute carrier family 6 (neurotransmitter transporter, serotonin), member 4
Sox18chr2:181669836-1816716409.605.15−0.902.05 × 10−2SRY (sex determining region Y)-box 18
Tcf7l2chr19:55741714-559336936.563.44−0.932.05 × 10−2transcription factor 7 like 2, T cell specific, HMG box
Tph2#,∆,¶,§,©chr10:115078553-1151850227.5213.790.874.87 × 10−2tryptophan hydroxylase 2
Losers vs. Control Mice
Gene SymbolLocusControl MiceLosersLog2 (Fold Change) Losers/Controlsq_ValueDefinition
Akap12chr10:4266328-435960530.7419.84−0.630.0159A kinase (PRKA) anchor protein (gravin) 12
Anxa4chr6:86736652-867936455.402.37−1.190.0063annexin A4
Bag3chr7:128523582-12854697921.999.86−1.160.0063BCL2-associated athanogene 3
Cacna1e#,§chr1:154390344-1548844127.2611.610.680.0275calcium channel, voltage-dependent, R type, alpha 1E subunit
Calcachr7:114631477-11463635723.975.14−2.220.0063calcitonin/calcitonin-related polypeptide, alpha
Camk2a#,∆,¶,§,©chr18:60925325-6098899321.0432.700.640.0275calcium/calmodulin-dependent protein kinase II alpha
Cdkn1achr17:29090985-291007225.823.43−0.760.0275cyclin-dependent kinase inhibitor 1A (P21)
Cdkn1cchr7:143458338-1434610745.593.08−0.860.0063cyclin-dependent kinase inhibitor 1C (P57)
Col25a1chr3:130180797-1305998863.255.240.690.0205collagen, type XXV, alpha 1
Colqchr14:31523083-315773832.140.78−1.450.0063collagen-like tail subunit (single strand of homotrimer) of asymmetric acetylcholinesterase
Dao#,§chr5:114003734-1140256765.172.55−1.020.0063D-amino acid oxidase
Ddc#,§chr11:11814099-1189814414.2935.121.300.0063dopa decarboxylase
Dnaja1chr4:40720153-4075788541.15155.221.920.0063DnaJ heat shock protein family (Hsp40) member A1
E2f1chr2:154559399-1545698924.752.38−0.990.0063E2F transcription factor 1
Egr1#,∆chr18:34861206-348649569.6215.750.710.0115early growth response 1
En1#,§chr1:120602486-1206079915.3210.270.950.0063engrailed 1
En2#,∆chr5:28165695-281721663.095.950.940.0063engrailed 2
Fev#,∆,¶,§chr1:74881508-748854081.253.951.660.0063FEV (ETS oncogene family)
Foxp2#chr6:14901169-154424501.593.150.980.0063forkhead box P2
Gabra1#,§chr11:42130920-4218306645.2677.940.780.0063gamma-aminobutyric acid (GABA) A receptor, subunit alpha 1
Gabrb2chr11:42419746-4263269419.5529.080.570.0446gamma-aminobutyric acid (GABA) A receptor, subunit beta 2
Gabrg2#,§chr11:41909958-4200087425.1438.210.600.0159gamma-aminobutyric acid (GABA) A receptor, subunit gamma 2
Gad1chr2:70489939-7060202948.9379.730.700.0063glutamate decarboxylase 1
Gad2#,∆,§,©chr2:22622057-2269407834.0454.290.670.0243glutamic acid decarboxylase 2
Gjc2chr11:59175476-5918321341.3224.81−0.740.0115gap junction protein, gamma 2
Inhbbchr1:119415464-1194222483.442.15−0.680.0374inhibin beta-B
Irs4chrX:141710996-1417252541.763.330.920.0063insulin receptor substrate 4
Isl1chr13:116298269-1163096931.610.49−1.720.0063ISL1 transcription factor, LIM/homeodomain
Lbx1chr19:45232619-452353777.040.71−3.320.0063ladybird homeobox homolog 1 (Drosophila)
Mafchr8:115682935-1157079409.485.55−0.770.0063avian musculoaponeurotic fibrosarcoma oncogene homolog
Mafbchr2:160363676-1603670654.472.51−0.840.0063v-maf musculoaponeurotic fibrosarcoma oncogene family, protein B (avian)
Mecomchr3:29951264-300132041.230.60−1.040.0063MDS1 and EVI1 complex locus
Mmp9chr2:164948218-1649558492.260.89−1.340.0063matrix metallopeptidase 9
Myh6chr14:54936101-549666071.4420.923.860.0243myosin, heavy polypeptide 6, cardiac muscle, alpha
Nefhchr11:4938562-5042794287.73119.70−1.270.0063neurofilament, heavy polypeptide
Neflchr14:68082365-68092351914.72422.75−1.110.0063neurofilament, light polypeptide
Nefmchr14:68119465-68131373690.57269.10−1.360.0063neurofilament, medium polypeptide
Neurod2chr11:98325416-983296452.921.59−0.880.0243neurogenic differentiation 2
Ngbchr12:87097530-871025399.2520.211.130.0063neuroglobin
Ngfrchr11:95568820-955876983.481.91−0.860.0159nerve growth factor receptor (TNFR superfamily, member 16)
Nkx6-1chr5:101658055-1016653614.911.66−1.560.0063NK6 homeobox 1
Nrgnchr9:37544492-3755274521.9542.780.960.0063neurogranin
Ntng2chr2:29193828-2925370725.6114.93−0.780.0275netrin G2
Otx2chr14:48657676-486676440.732.031.490.0115orthodenticle homeobox 2
Pax7chr4:139737093-1398329680.451.641.850.0063paired box 7
Phox2bchr5:67094327-670991267.031.56−2.170.0063paired-like homeobox 2b
Prss12#,§chr3:123446912-1235066024.231.89−1.160.0063protease, serine 12 neurotrypsin (motopsin)
Pvalb#,§chr15:78191044-78206351119.0970.07−0.770.0063parvalbumin
Scn3achr2:65456990-655675193.866.640.780.0346sodium channel, voltage-gated, type III, alpha
Slc17a8chr10:89574019-896212490.902.121.240.0063solute carrier family 17 (sodium-dependent inorganic phosphate cotransporter), member 8
Slc18a3chr14:32407354-3246600433.8313.30−1.350.0063solute carrier family 18 (vesicular monoamine), member 3
Slc6a2#,§,©chr8:92960641-930016677.524.87−0.630.0346solute carrier family 6 (neurotransmitter transporter, noradrenalin), member 2
Slc6a3#,§,©chr13:73536746-735786721.039.513.210.0063solute carrier family 6 (neurotransmitter transporter, dopamine), member 3
Slc6a4#,§,©chr11:76998596-770323434.9422.712.200.0063solute carrier family 6 (neurotransmitter transporter, serotonin), member 4
Slc6a5chr7:49909938-4996386145.0613.79−1.710.0063solute carrier family 6 (neurotransmitter transporter, glycine), member 5
Spp1chr5:104435110-10444105354.8023.57−1.220.0063secreted phosphoprotein 1
Syt1chr10:108497647-10901098351.3286.930.760.0063synaptotagmin I
Tgfb3chr12:86056580-860791596.994.46−0.650.0374transforming growth factor, beta 3
Th#,§chr7:142892670-14290196010.9321.961.010.0063tyrosine hydroxylase
Tph2#,∆,¶,§,©chr10:115078553-1151850227.4026.971.870.0063tryptophan hydroxylase 2
Vdr#chr15:97854426-979082960.760.29−1.410.0115vitamin D receptor
Vgfchr5:137025191-137033857119.2478.61−0.600.0473VGF nerve growth factor inducible
Vipr2chr12:116077725-1161462611.830.39−2.240.0063vasoactive intestinal peptide receptor 2
Whrnchr4:63414854-6349613021.7213.85−0.650.0275whirlin
Genes associated with: #—abnormal emotion/affect behavior; ∆—abnormal aggression-related behavior; —increased aggression towards mice; §-abnormal fear/anxiety-related behavior; ©—abnormal depression-related behavior.
Table
Table 2. DEGs encoding transcription factors.
Table 2. DEGs encoding transcription factors.
Gene SymbolLocusExpression in VTA (FPKM)log2 (Fold Change) Winners/Controlsq_ValueDefinition
Control MiceWinners
Winners vs. Control Mice
Ddx5chr11:106780156-106818861121.77381.381.652.05 × 10−2DEAD (Asp-Glu-Ala-Asp) box polypeptide 5
Ercc2 *chr7:19382024-1940410453.428.72−2.622.05 × 10−2excision repair cross-complementing rodent repair deficiency,
complementation group 2
Fubp3chr2:31572650-3161759011.2160.262.432.05 × 10−2far upstream element (FUSE) binding protein 3
Lbx1 *chr19:45232619-452353777.183.19−1.172.05 × 10−2ladybird homeobox homolog 1 (Drosophila)
Maf *chr8:115682935-1157079409.665.98−0.693.58 × 10−2avian musculoaponeurotic fibrosarcoma oncogene homolog
Nkx6-1 *chr5:101658055-1016653615.012.79−0.853.58 × 10−2NK6 homeobox 1
Parp1chr1:180568937-18060138912.93102.042.982.05 × 10−2poly (ADP-ribose) polymerase family, member 1
Tcf24chr1:9960162-99674850.390.891.194.87 × 10−2transcription factor 24
Tcf7l2 *chr19:55741714-559336936.563.44−0.932.05 × 10−2transcription factor 7 like 2, T cell specific, HMG box
Losers vs. Control Mice
Gene SymbolLocusControl MiceLoserslog2 (Fold Change) Losers/Controlsq_ValueDefinition
Ctdsplchr9:118921134-1190443535.1426.562.370.0063CTD (carboxy-terminal domain, RNA polymerase II,
polypeptide A) small phosphatase-like
E2f1 *chr2:154559399-1545698924.752.38−0.990.0063E2F transcription factor 1
Egr1 *,#,∆chr18:34861206-348649569.6215.750.710.0115early growth response 1
En1 *,#,§chr1:120602486-1206079915.3210.270.950.0063engrailed 1
En2 *,#,∆chr5:28165695-281721663.095.950.940.0063engrailed 2
Fev *,#,∆,¶,§chr1:74881508-748854081.253.951.660.0063FEV (ETS oncogene family)
Foxp2 *,#chr6:14901169-154424501.593.150.980.0063forkhead box P2
Hoxb2chr11:96323125-963541743.720.64−2.540.0063homeobox B2
Hoxb3chr11:96323125-963541742.980.31−3.260.0063homeobox B3
Hoxd3chr2:74710043-747651421.820.08−4.530.0115homeobox D3
Isl1 *chr13:116298269-1163096931.610.49−1.720.0063ISL1 transcription factor, LIM/homeodomain
Lbx1 *chr19:45232619-452353777.040.71−3.320.0063ladybird homeobox homolog 1 (Drosophila)
Maf *chr8:115682935-1157079409.485.55−0.770.0063avian musculoaponeurotic fibrosarcoma oncogene homolog
Mafb *chr2:160363676-1603670654.472.51−0.840.0063v-maf musculoaponeurotic fibrosarcoma oncogene family, protein B (avian)
Neurod2 *chr11:98325416-983296452.921.59−0.880.0243neurogenic differentiation 2
Nkx6-1 *chr5:101658055-1016653614.911.66−1.560.0063NK6 homeobox 1
Otx1chr11:21994763-220016510.622.502.010.0063orthodenticle homeobox 1
Otx2 *chr14:48657676-486676440.732.031.490.01154orthodenticle homeobox 2
Pax2chr19:44746709-448382668.343.02−1.470.0063paired box 2
Pax7 *chr4:139737093-1398329680.451.641.850.0063paired box 7
Phox2b *chr5:67094327-670991267.031.56−2.170.0063paired-like homeobox 2b
Rorcchr3:94372793-943982742.971.40−1.080.0063RAR-related orphan receptor gamma
Scrt2chr2:152081528-15209580213.307.04−0.920.0063scratch family zinc finger 2
Sox14chr9:99874105-998761700.742.211.570.0063SRY (sex determining region Y)-box 14
Tal1chr4:115056425-1150717582.454.250.800.0205T cell acute lymphocytic leukemia 1
Tbx20chr9:24720811-247743030.860.25−1.770.0063T-box 20
Vdr *,#chr15:97854426-979082960.760.29−1.410.0115vitamin D receptor
Zkscan16chr4:58943559-589627057.7212.670.710.0063zinc finger with KRAB and SCAN domains 16
Genes associated with: *—behavior/neurological phenotype; #—abnormal emotion/affect behavior; ∆—abnormal aggression-related behavior; —increased aggression towards mice; §—abnormal fear/anxiety-related behavior.
Table
Table 3. Metabolic pathways most significantly altered in defeated mice compared to controls.
Table 3. Metabolic pathways most significantly altered in defeated mice compared to controls.
TermCountp ValueGenes
Nicotine addiction54.25 × 10−4Slc17a8, Gabrg2, Gabra1, Gabrb2, Chrna6
Axon guidance65.86 × 10−3Plxnc1, Nrp1, Sema7a, Sema3f, Ntng2, Slit1
GABAergic synapse57.55 × 10−3Gabrg2, Gad2, Gabra1, Gabrb2, Gad1
Amyotrophic lateral sclerosis (ALS)41.06 × 10−2Prph, Nefh, Nefl, Nefm
Amphetamine addiction42.20 × 10−2Ddc, Slc6a3, Th, Camk2a
Chronic myeloid leukemia42.65 × 10−2E2f1, Cdkn1a, Tgfb3, Mecom
Bladder cancer35.20 × 10−2E2f1, Cdkn1a, Mmp9
Transcriptional misregulation in cancer56.03 × 10−2Maf, Cdkn1a, Mmp9, Ngfr, Fev
Retrograde endocannabinoid signaling46.48 × 10−2Slc17a8, Gabrg2, Gabra1, Gabrb2
Cocaine addiction37.12 × 10−2Ddc, Slc6a3, Th
Taurine and hypotaurine metabolism29.44 × 10−2Gad2, Gad1
Table
Table 4. DEGs with changed expression both for the winners and losers (common DEGs).
Table 4. DEGs with changed expression both for the winners and losers (common DEGs).
Gene SymbolLog2 (Fold_Change) Winners/Controlsq_ValueLog2 (Fold_Change) Losers/Controlsq_ValueDefinition
Apba31.142.05 × 10−21.386.29 × 10−3amyloid beta (A4) precursor protein-binding, family A, member 3
Brd3−1.912.05 × 10−2−1.696.29 × 10−3bromodomain containing 3
Cdk12−5.452.05 × 10−2−5.436.29 × 10−3cyclin-dependent kinase 12
Dnaja11.322.05 × 10−21.926.29 × 10−3DnaJ heat shock protein family (Hsp40) member A1
Egfl7−0.952.05 × 10−2−0.976.29 × 10−3EGF-like domain 7
Gcn1l1−4.172.05 × 10−2−3.966.29 × 10−3GCN1 general control of amino-acid synthesis 1-like 1 (yeast)
Gfap−0.902.05 × 10−2−1.436.29 × 10−3glial fibrillary acidic protein
Gm13889−0.733.58 × 10−2−0.826.29 × 10−3predicted gene 13889
Lbx1*−1.172.05 × 10−2−3.326.29 × 10−3ladybird homeobox homolog 1 (Drosophila)
Maf *−0.693.58 × 10−2−0.776.29 × 10−3avian musculoaponeurotic fibrosarcoma oncogene homolog
Myh14−0.644.87 × 10−2−0.796.29 × 10−3myosin, heavy polypeptide 14
Myh65.352.05 × 10−23.862.43 × 10−2myosin, heavy polypeptide 6, cardiac muscle, alpha
Myoc−1.892.05 × 10−2−2.086.29 × 10−3myocilin
Ngfr−0.822.05 × 10−2−0.861.59 × 10−2nerve growth factor receptor (TNFR superfamily, member 16)
Nkx6-1 *−0.853.58 × 10−2−1.566.29 × 10−3NK6 homeobox 1
Nrgn−0.663.58 × 10−20.966.29 × 10−3neurogranin
Prph−0.842.05 × 10−2−2.226.29 × 10−3peripherin
Rln3−1.702.05 × 10−2−3.546.29 × 10−3relaxin 3
Slc6a40.992.05 × 10−22.206.29 × 10−3solute carrier family 6 (neurotransmitter transporter, serotonin), member 4
Smo−0.852.05 × 10−2−0.823.74 × 10−2smoothened, frizzled class receptor
Taf1d4.612.05 × 10−24.586.29 × 10−3TATA-box binding protein associated factor, RNA polymerase I, D
Tph20.874.87 × 10−21.876.29 × 10−3tryptophan hydroxylase 2
Vezt1.734.87 × 10−22.106.29 × 10−3vezatin, adherens junctions transmembrane protein
* Genes, encoding transcription factors.
Table
Table 5. Metabolic pathways most significantly altered in winners compared to defeated mice.
Table 5. Metabolic pathways most significantly altered in winners compared to defeated mice.
Gene SymbolMiceLog2 (Fold_Change) Losers/Winnersq_ValueDefinition
Winners FPKMLosers FPKM
Amphetamine addiction, p Value = 3.30 × 10−3
Camk2a22.6933.700.576.12 × 10−3calcium/calmodulin-dependent protein kinase II alpha
Ddc21.2136.210.776.12 × 10−3dopa decarboxylase
Fos3.976.680.756.12 × 10−3FBJ osteosarcoma oncogene
Slc6a31.539.772.686.12 × 10−3solute carrier family 6 (neurotransmitter transporter, dopamine), member 3
Th12.7822.600.826.12 × 10−3tyrosine hydroxylase
Dopaminergic synapse, p Value = 7.74 × 10−3
Caly79.62118.760.586.12 × 10−3calcyon neuron-specific vesicular protein
Camk2a22.6933.700.576.12 × 10−3calcium/calmodulin-dependent protein kinase II alpha
Ddc21.2136.210.776.12 × 10−3dopa decarboxylase
Fos3.976.680.756.12 × 10−3FBJ osteosarcoma oncogene
Slc6a31.539.772.686.12 × 10−3solute carrier family 6 (neurotransmitter transporter, dopamine), member 3
Th12.7822.600.826.12 × 10−3tyrosine hydroxylase
Adrenergic signaling in cardiomyocytes, p Value = 9.81 × 10−3
Cacna2d314.229.93−0.522.68 × 10−2calcium channel, voltage-dependent, alpha2/delta subunit 3
Cacnb314.4921.100.546.12 × 10−3calcium channel, voltage-dependent, beta 3 subunit
Camk2a22.6933.700.576.12 × 10−3calcium/calmodulin-dependent protein kinase II alpha
Scn4b73.3550.16−0.551.54 × 10−2sodium channel, type IV, beta
Scn5a0.430.840.996.12 × 10−3sodium channel, voltage-gated, type V, alpha
Tnnt22.354.390.901.09 × 10−2troponin T2, cardiac
Amyotrophic lateral sclerosis, p Value = 1.14 × 10−2
Nefh253.37123.42−1.046.12 × 10−3neurofilament, heavy polypeptide
Nefl710.19436.14−0.706.12 × 10−3neurofilament, light polypeptide
Nefm557.04277.66−1.006.12 × 10−3neurofilament, medium polypeptide
Prph26.4410.15−1.386.12 × 10−3peripherin
Axon guidance, p Value = 3.10 × 10−2
Nrp15.3841.272.946.12 × 10−3neuropilin 1
Plxnc13.034.520.581.54 × 10−2plexin C1
Sema3f5.268.440.684.25 × 10−2sema domain, immunoglobulin domain (Ig), short basic domain, secreted, (semaphorin) 3F
Slit18.8913.380.596.12 × 10−3slit homolog 1 (Drosophila)
Unc5d2.723.780.484.25 × 10−2unc-5 netrin receptor D
ECM-receptor interaction, p Value = 4.74 × 10−2
Col27a13.922.29−0.786.12 × 10−3collagen, type XXVII, alpha 1
Spp159.6424.33−1.296.12 × 10−3secreted phosphoprotein 1
Sv2b20.7214.58−0.511.09 × 10−2synaptic vesicle glycoprotein 2 b
Sv2c36.9226.11−0.501.09 × 10−2synaptic vesicle glycoprotein 2c
Tyrosine metabolism, p Value = 5.00 × 10−2
Dbh21.5911.27−0.946.12 × 10−3dopamine beta hydroxylase
Ddc21.2136.210.776.12 × 10−3dopa decarboxylase
Th12.7822.600.826.12 × 10−3tyrosine hydroxylase
Cocaine addiction, p Value = 7.49 × 10−2
Ddc21.2136.210.776.12 × 10−3dopa decarboxylase
Slc6a31.539.772.686.12 × 10−3solute carrier family 6 (neurotransmitter transporter, dopamine), member 3
Th12.7822.600.826.12 × 10−3tyrosine hydroxylase
MAPK signaling pathway, p Value = 8.15 × 10−2
Cacna1h3.635.580.626.12 × 10−3calcium channel, voltage-dependent, T type, alpha 1H subunit
Cacna2d314.229.93−0.522.68 × 10−2calcium channel, voltage-dependent, alpha2/delta subunit 3
Cacnb314.4921.100.546.12 × 10−3calcium channel, voltage-dependent, beta 3 subunit
Fos3.976.680.756.12 × 10−3FBJ osteosarcoma oncogene
Hspb116.739.18−0.876.12 × 10−3heat shock protein 1
Mecom1.060.62−0.792.32 × 10−2MDS1 and EVI1 complex locus
Table
Table 6. Expression of genes that changed the level of transcription in comparisons of Losers vs. Control and Winners vs. Control in different directions and characterized as DEGs when comparing Losers vs. Winners.
Table 6. Expression of genes that changed the level of transcription in comparisons of Losers vs. Control and Winners vs. Control in different directions and characterized as DEGs when comparing Losers vs. Winners.
Gene SymbolLosers/ControlWinners/ControlLosers/WinnersDefinition
Controls FPKMLosers FPKMLog2 (Fold_Change) Losers/Controlsq ValueControls FPKMWinners FPKMLog2 (Fold_Change) Winners/Controlsq ValueLog2 (Fold_Change) Losers/winnersq Value
Ercc2 *,#,∆52.51143.521.455.07 × 10−253.428.72−2.622.05 × 10−24.066.12 × 10−3excision repair cross-complementing rodent repair deficiency, complementation group 2
Loxl21.040.69−0.599.99 × 10−11.0646.255.452.05 × 10−2−6.046.12 × 10−3lysyl oxidase-like 2
Nrgn#21.9542.780.966.29 × 10−322.3514.15−0.663.58 × 10−21.626.12 × 10−3neurogranin
Otx2 *,#0.732.031.491.15 × 10−20.740.32−1.196.96 × 10−12.676.12 × 10−3orthodenticle homeobox 2
Six3 *1.936.561.777.09 × 10−11.960.68−1.529.99 × 10−13.286.12 × 10−3sine oculis-related homeobox 3
Tcf7l2 *,#6.449.720.591.07 × 10−16.563.44−0.932.05 × 10−21.526.12 × 10−3transcription factor 7 like 2, T cell specific, HMG box
*—DEGs encoding transcription factors; Genes associated with: #—behavior/neurological phenotype; ∆—abnormal emotion/affect behavior.
Table
Table 7. Metabolic pathways that include genes that maximize intergroup differences when comparing winners with defeated mice.
Table 7. Metabolic pathways that include genes that maximize intergroup differences when comparing winners with defeated mice.
Term IDPathwayCountp ValueGenes
mmu04916Melanogenesis88.60 × 10−3Wnt3, Adcy9, Adcy6, Mitf, Edn1, Crebbp, Camk2a *, Fzd7
mmu04020Calcium signaling pathway119.76 × 10−3Gnal, P2rx6, Slc8a2 *, Adcy9, Ryr1, Cacna1g, Ppp3cc, Plcd4, Camk2a*, Htr5a, Cacna1b
mmu04713Circadian entrainment72.77 × 10−2Adcy9, Adcy6, Per2, Ryr1, Cacna1g, Per3, Camk2a *
mmu04141Protein processing in endoplasmic reticulum94.39 × 10−2Hsph1, Eif2ak1, Dnajb11, Man1a, Ubqln2, Ubqln1, Eif2ak2, Hspa8, Bcap31
mmu04350TGF-beta signaling pathway64.95 × 10−2Acvr2b, Nog, Zfyve9, Crebbp, Chrd, Pitx2 *
mmu04261Adrenergic signaling in cardiomyocytes85.00 × 10−2Adcy9, Adcy6, Scn4b *, Cacnb3 *, Myh6, Cacna2d3 *, Scn5a *, Camk2a *
*—Genes, differentially expressed in VTA of winners and defeated mice (according to the Cufflinks/Cuffdiff programs).
Table
Table 8. Correlation between expression of the key genes (Nrgn, Camk2a, Ercc2, Otx2, and Six3) and genes related to the dopamine signaling pathway.
Table 8. Correlation between expression of the key genes (Nrgn, Camk2a, Ercc2, Otx2, and Six3) and genes related to the dopamine signaling pathway.
Correlation in VTA of Winners
Gene SymbolNrgnCamk2aErcc2Otx2Six3
Drd20.7220.920 **0.6050.895 *0.902 *
Slc6a30.5850.881 *0.5180.885 *0.893 *
Th0.5170.6960.2020.5270.537
Ddc0.4750.865 *0.3240.7310.749
Correlation in VTA of Losers
Gene SymbolNrgnCamk2aErcc2Otx2Six3
Drd20.991 ***0.970 **0.909 *0.924 **0.944 **
Slc6a30.953 **0.965 **0.881 *0.917 **0.941 **
Th0.910 *0.7890.7470.6510.703
Ddc0.865 *0.942 **0.860 *0.951 **0.954 **
Level of significance (two-tailed test): * <0.05; ** <0.01; *** <0.001.

Share and Cite

MDPI and ACS Style

Redina, O.; Babenko, V.; Smagin, D.; Kovalenko, I.; Galyamina, A.; Efimov, V.; Kudryavtseva, N. Gene Expression Changes in the Ventral Tegmental Area of Male Mice with Alternative Social Behavior Experience in Chronic Agonistic Interactions. Int. J. Mol. Sci. 2020, 21, 6599. https://doi.org/10.3390/ijms21186599

AMA Style

Redina O, Babenko V, Smagin D, Kovalenko I, Galyamina A, Efimov V, Kudryavtseva N. Gene Expression Changes in the Ventral Tegmental Area of Male Mice with Alternative Social Behavior Experience in Chronic Agonistic Interactions. International Journal of Molecular Sciences. 2020; 21(18):6599. https://doi.org/10.3390/ijms21186599

Chicago/Turabian Style

Redina, Olga, Vladimir Babenko, Dmitry Smagin, Irina Kovalenko, Anna Galyamina, Vadim Efimov, and Natalia Kudryavtseva. 2020. "Gene Expression Changes in the Ventral Tegmental Area of Male Mice with Alternative Social Behavior Experience in Chronic Agonistic Interactions" International Journal of Molecular Sciences 21, no. 18: 6599. https://doi.org/10.3390/ijms21186599

APA Style

Redina, O., Babenko, V., Smagin, D., Kovalenko, I., Galyamina, A., Efimov, V., & Kudryavtseva, N. (2020). Gene Expression Changes in the Ventral Tegmental Area of Male Mice with Alternative Social Behavior Experience in Chronic Agonistic Interactions. International Journal of Molecular Sciences, 21(18), 6599. https://doi.org/10.3390/ijms21186599

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop