Next Article in Journal
Protective Effect of D-Panthenol in Rhabdomyolysis-Induced Acute Kidney Injury
Next Article in Special Issue
Mother–Fetus Immune Cross-Talk Coordinates “Extrinsic”/“Intrinsic” Embryo Gene Expression Noise and Growth Stability
Previous Article in Journal
Non-Coding RNAs Are Implicit in Chronic Myeloid Leukemia Therapy Resistance
Previous Article in Special Issue
Formation of Unique Placental Seed Capsules in the Maturation Process of the Tomato Fruit
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 23
Issue 20
10.3390/ijms232012269
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

Transcription Factors as Important Regulators of Changes in Behavior through Domestication of Gray Rats: Quantitative Data from RNA Sequencing

by
Dmitry Oshchepkov
*,
Irina Chadaeva
,
Rimma Kozhemyakina
,
Svetlana Shikhevich
,
Ekaterina Sharypova
,
Ludmila Savinkova
,
Natalya V. Klimova
,
Anton Tsukanov
,
Victor G. Levitsky
and
Arcady L. Markel
Institute of Cytology and Genetics, 630090 Novosibirsk, Russia
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2022, 23(20), 12269; https://doi.org/10.3390/ijms232012269
Submission received: 31 August 2022 / Revised: 28 September 2022 / Accepted: 12 October 2022 / Published: 14 October 2022
(This article belongs to the Special Issue Developmental Biology: Computational and Experimental Approaches)
Browse Figures
Versions Notes

Abstract

:
Studies on hereditary fixation of the tame-behavior phenotype during animal domestication remain relevant and important because they are of both basic research and applied significance. In model animals, gray rats Rattus norvegicus bred for either an enhancement or reduction in defensive response to humans, for the first time, we used high-throughput RNA sequencing to investigate differential expression of genes in tissue samples from the tegmental region of the midbrain in 2-month-old rats showing either tame or aggressive behavior. A total of 42 differentially expressed genes (DEGs; adjusted p-value  <  0.01 and fold-change  >  2) were identified, with 20 upregulated and 22 downregulated genes in the tissue samples from tame rats compared with aggressive rats. Among them, three genes encoding transcription factors (TFs) were detected: Ascl3 was upregulated, whereas Fos and Fosb were downregulated in tissue samples from the brains of tame rats brain. Other DEGs were annotated as associated with extracellular matrix components, transporter proteins, the neurotransmitter system, signaling molecules, and immune system proteins. We believe that these DEGs encode proteins that constitute a multifactorial system determining the behavior for which the rats have been artificially selected. We demonstrated that several structural subtypes of E-box motifs—known as binding sites for many developmental TFs of the bHLH class, including the ASCL subfamily of TFs—are enriched in the set of promoters of the DEGs downregulated in the tissue samples of tame rats’. Because ASCL3 may act as a repressor on target genes of other developmental TFs of the bHLH class, we hypothesize that the expression of TF gene Ascl3 in tame rats indicates longer neurogenesis (as compared to aggressive rats), which is a sign of neoteny and domestication. Thus, our domestication model shows a new function of TF ASCL3: it may play the most important role in behavioral changes in animals.

1. Introduction

It is well known that the development of the human civilization is inextricably linked with the domestication of animals. An important factor in the domestication of various species of wild animals has been the transformation of their aggressive behavior toward humans into tolerant and even friendly behavior. Despite a long and successful history of domestication, the genetics of tame and aggressive behavior of domesticated animals toward humans, as well as its involvement in the regulation of these genes, remains unclear.
Modern studies on the transcriptomes of various domestic animal species indicate an influence of a large number of genes on changes in the behavior, coloration, morphology, and physiology of domestic animals [1,2,3,4]. Some of these genes are common among the gene lists published in research articles—usually the genes that are involved in the control of behavior (for example, Bdnf, Grin1, Oxt, and dopamine receptor genes, among others) and various genes that are not related to each other [5,6]. Nonetheless, it is known that behavioral patterns are determined by multiple genes; therefore, one or several genes differentially expressed between tame and wild animals cannot explain the entire spectrum of so-called domestication syndrome [7], which includes alterations in skull morphology [8], coat/feather color [9], the reproductive cycle [10,11,12], and other features. However, the involvement of TF genes in transcriptome changes during the domestication of model animals has not yet been studied.
To determine the molecular mechanisms involved in changing of behavior toward humans in wild animals upon domestication, we used an experimental model: gray rats selected for 45 years for their reaction to humans (either tame or aggressive behavior) at the animal facility of the Institute of Cytology and Genetics, SB RAS (Novosibirsk, Russia). This animal model of domestication is a unique experimental system for analysis of the complex architecture of genetic determination of domestic behavior; in tame and aggressive rats, it is possible to correctly and adequately assess phenotypic manifestations of nonaggressive and aggressive reactions, respectively, upon exposure to a stimulus (for example, an experimenter’s gloved hand, i.e., the glove test) [13,14].
We previously performed genome-wide transcriptomic analysis of samples of the hypothalamus and hippocampus to identify genes differentially expressed between tame and aggressive rats [5,6]; these brain regions belong to the limbic system, which controls emotions [15,16]. The resultant lists of differentially expressed genes (DEGs) were found to contain several TF genes, including Ascl3 (achaete-scute complex-like 3), which is a member of the basic helix–loop–helix factor class (bHLH) [17]. Typically, members of the bHLH class can form homo- and heterodimers; ASCL3 belongs to the subfamily of achaete-scute-like factors of bHLH TFs, which also includes tissue-specific factors. Members of this subfamily can serve as either activators or repressors [17,18]. TFs of this subfamily participate in neurogenesis, tissue development, and cell differentiation and are critical for proper development of the nervous system [19,20]. Proteins encoded by the genes of this subfamily are directly involved in the segregation of neuroblasts from other epidermal cells; therefore, Ascl genes are classified as proneural; the expression of such genes enables ectodermal cells to differentiate into nerve cells [21].
According to the latest data [22], ASCL3 protein belongs to the set of candidate proteins involved in the development of mesencephalic dopaminergic neurons in the midbrain. This brain region contains clusters of serotonergic and dopaminergic neurons and is responsible for inhibitory control of aggressive behavior [23,24,25,26]. The midbrain also includes the periaqueductal gray and substantia nigra. The latter is a group of dopaminergic neurons that form a dopaminergic nucleus. These neurons project to structures of the so-called limbic system, including the previously studied hypothalamus and hippocampus. Therefore, the aim of our work was to determine whether a statistically significant change in gene expression in tame and aggressive rats is a consequence of the differential expression of transcription factor genes, in particular Ascl3, which can regulate these genes. Thus, in this work, to elucidate molecular mechanisms underlying the behavioral changes that occur in animals during domestication, we analyzed differences in gene expression profiles between gray rats (Rattus norvegicus) selected for either the absence or enhancement of a defensive response to humans (tame and aggressive strains, respectively). We applied two methods: (i) high-throughput RNA sequencing (RNA-seq) analysis of brain samples from the tegmental region of the midbrain (MT) region and (ii) motif enrichment analysis using a newly developed ESDEG tool; these methods helped to investigate coordinated changes in the transcriptome of the animal domestication model.

2. Results

2.1. Analysis of RNA-Seq Data from Brain Tissue Samples of Tame and Aggressive Rats

Using an Illumina NextSeq 550 system, we sequenced the MT transcriptome of three tame adult male gray rats (Rattus norvegicus) and three aggressive rats devoid of any family relationship among them (see Materials and Methods, Section 4.1). This procedure resulted in 182,197,974 reads (each 75 bp in length), which were deposited in the NCBI SRA database (PRJNA668014), as shown in Table 1. Among them, using the pipeline described in Materials and Methods (Section 4.2), we mapped 158,310,590 reads (86.9%) to rat reference genome Rn6 (assembly 2014) (see Table 1). Then, we identified 14,039 genes expressed within the MT of rats under our experimental conditions (Table 1). To minimize the false-positive error rates, we also applied Benjamini–Hochberg correction for multiple comparisons, resulting in 42 DEGs within the MT of tame rats versus aggressive rats (Table 1).
These DEGs are listed in Table 2. In particular, the Ascl3 gene encoding achaete-scute complex-like transcription factor 3 showed the highest significance of enrichment (PADJ = 3.16 × 10−4) and a log2 (fold change) of −2.82, ranking fifth among all upregulated DEGs in the tame rat samples (fold change or fold is a ratio of an expression level of a given gene in tame rats to that in aggressive rats, hereafter referred to as log2 fold; see Table 2).
Thus, we obtained a list of differentially expressed genes between tame and aggressive rats in the MT (42 genes in total; Table 2), among which genes encoding three TFs (ASCL3, FOS, and FOSB) were identified, in addition to ARC (known as the master regulator of synaptic function and of numerous neuronal signaling pathways), extracellular-matrix components (HSPA1A, HSPA1B, TAC3, and VIP), transporter protein (SLC4A5), neurotransmitter (HCRT), signaling molecules (serotonin receptors HTR3A and HTR5B, LILRB3L, and RETSAT), and immune system proteins (DEFB17 and VIP). The observed significant differences in the expression levels of these genes between tame and aggressive rats indicate major biochemical and physiological rearrangements in the activity of neurons at the metabolic level, along with changes in behavior. The proteins encoded by the DEGs differ both in terms of their functions and in the neurotransmitter systems with which they are affiliated; this state of affairs results in a multifactorial system determining the trait for which our rats were artificially selected.
Our list of DEGs was functionally annotated in the David Bioinformatics Resources database [27]. The only category statistically significantly (p = 0.0089) enriched in this gene set was a group of genes united by the term “Signal” (UniProt Category Domain, KW-0732): Htr3a, Fcgr2b, Lypd3, Cthrc1, Col15a1, Defb17, Hcrt, Lilrb3l, Mpeg1, Olfml1, Pcdhga1, Retsat, Spint1, Tac3, and Vip, in agreement with the current understanding of the relation of signaling pathways with behavioral response. The absence of significant enrichment with more specific terms can be explained by the multifactorial nature of behavioral responses and of the evolutionary process.

2.2. Gene Expression Analysis in the MT of Tame and Aggressive Rats by qPCR

The differential gene expression results obtained in the RNA-seq experiment (Table 1) were experimentally validated by qPCR (in the MT samples of eight tame and eight aggressive rats from the 90th generation of selection) on two selected genes—Ascl3 (achaete-scute complex-like transcription factor 3) and Defb17 (defensin beta 17)—as well as two reference genes: Ppia and Rpl30. The results of the RNA-seq analysis were successfully confirmed; expression levels of genes Ascl3 and Defb17 were significantly higher in the tame rats than in aggressive rats (t-test, p < 5 × 10−6, Figure 1). Expression results for the reference genes did not differ significantly between the tame and aggressive rats according to qPCR analysis.

2.3. Analysis of Enrichment with Binding Sites for TFs

Above, we revealed that TF genes Ascl3 and Fos/Fosb were significantly upregulated and downregulated, respectively, in the brain tissue samples of tame rats (Table 2). These TFs are included in the standard classification analysis of murine TFs [17]. Fos-related factors {1.1.2} (hereafter, numbers in curly brackets correspond to the classification of murine TFs based on their DNA-binding domains [17]) and MyoD/ASC-related factors {1.2.2}—respective TF families for these factors—possess well-known motifs in the publicly available Hocomoco library of motifs, which is based on a massive analysis of ChIP-seq data [28]. Therefore, we were able to test the overrepresentation of motifs of known TFs in the promoters of the DEGs. The standard methods for motif enrichment analysis utilizing the exact Fisher’s test appeared to be ineffective due to the small sizes of our promoter datasets (20 and 22 DEGs; see Table 2). Accordingly, we developed an original method called ESDEG to test the overrepresentation of motifs independent of dataset size (see Section 4.4, Materials and Methods). We applied ESDEG to the promoter sets of the genes upregulated and downregulated in the brain tissue samples of tame rats (Table 2), employing the library of motifs for known murine TFs (see Section 4.4, Materials and Methods). The names of the respective TF families were compiled based on their DNA-binding domains [17] for enriched motifs to simplify the interpretation of the results. Additionally, we performed a permutation test [29] to estimate the significance of pairwise similarity between the motifs that were found to be enriched in each set. The resulting tables of overrepresented motifs for known TFs detected in promoters of the genes upregulated and downregulated in the brain tissue samples of tame rats are shown in Figure 2A,B, respectively. The list of 25 motifs corresponds to the set of promoters of the genes downregulated in the brain tissue samples of tame rats, whereas the other list, corresponding to the set of promoters of the genes upregulated in the brain tissue samples of tame rats, contains only eight motifs. The IRF3 motif ranks first in the motif enrichment list for the genes downregulated in the brain tissue samples of tame rats. Because the Irf6 gene showed an upregulation trend in the RNA-seq experiment (although it does not fulfill the DEG criteria: padj = 0.051 < 0.1, log2 fold = (0.87, see Table S1), we wondered whether it could be responsible for the enrichment of the IRF3 motif. We found that the IRF6 motif does not possess a significant similarity to the rest of the TFs that are members of the interferon-regulatory factor family {3.5.3} (p > 0.05, [29]); thus, the relation of slightly upregulated IRF6 TF and IRF3 motif enrichment seems insignificant. A large proportion of the motif enrichment list for genes downregulated in the brain tissue samples of tame rats (9 of 25) are motifs for TFs of closely related families: E2A-related factors {1.2.1}, MyoD/ASC-related factors {1.2.2}, and Tal-related factors {1.2.3}. All the mentioned TF families belong to the bHLH class. The motifs of TFs from these families are very similar and are known as E boxes. These motifs showed a highly significant mutual similarity (Figure 2A, 0.01 < (log10(p-value) < 1 × 10−9), coinciding with the significant upregulation of the Ascl3 gene in the brain tissue samples of tame rats (PADJ = 3.16 × 10−4 and log2 fold = 2.82, see Table 2). A recent large all-against-all benchmarking of positional weight matrix models [30] was used for recognition of motifs. It was demonstrated here that the best-performing positional weight matrix for a given TF is often related to another TF, usually from the same family. Thus, we propose that the overrepresented bHLH motifs found in promoters of the genes downregulated in the brain tissue samples of tame rats match potential binding sites for ASCL3 TF with related partner TFs from the same bHLH class. Thus, we compiled a sizeable list of ASCL3 targets among genes downregulated in the brain tissue samples of tame rats; this is consistent with the fact that the ASCL3 TF can act as a repressor toward positively acting bHLH TFs [31]. We did not find any other TF genes in the lists of DEGs (Table 2) corresponding to other enriched motifs (Figure 2).
We noticed very significant enrichment with motifs for TFs from the POU domain factor family {3.1.10} for promoters of the genes upregulated in the brain tissue samples of tame rats. Although we did not detect any differential expression for any gene from the POU domain factor family, this finding is worth a mention because the TFs from this family regulate a wide range of developmental processes—from specification of the early embryo to terminal differentiation. Moreover, the most enriched motif in this group belongs to the POU2F2 TF, which is known to be a regulator of neuronal development and differentiation [32]. We suggest that this observation deserves further investigation. Furthermore, we noted that the milder threshold for the fold ratio (PADJ < 0.1, fold change > 1.5 or log2 fold > 0.585) abrogates the enrichment of motifs for bHLH TFs in the promoters of this extended list of DEGs (see Table S1) regarded as downregulated in the brain tissue samples of tame rats at this threshold, although it revealed the FOSL2_MOUSE.H11MO.0.A motif from the Hocomoco database in the motif enrichment list for the promoter set of the genes regarded as upregulated in the brain tissue samples of tame rats at this threshold (see Figure S1). Consequently, we believe that this motif enrichment may be consistent with the finding that FOS/FOSB TFs are downregulated in the brain tissue samples of tame rats (Table 2).
In conclusion, we identified a correspondence between the overexpression of TF gene Ascl3 in tame rats and changes in the expression of genes in the list of DEGs; TF gene Ascl3 is one of the top-scoring genes in the list of DEGs upregulated in the brain tissue samples of tame rats (Table 2). Additionally, we detected that the motifs matching potential binding sites for TF ASCL3 and significantly similar motifs for its putative partner TFs from the same bHLH class constitute a large proportion (9 of 25; Figure 2) of all motifs significantly enriched in the set of the genes downregulated in tame rats. Furthermore, Yoshida et al. [31] previously proved not only the possibility of an interaction between ASCL3 and another bHLH TF, MYOD1, but also that in vitro, ASCL3 acts as a repressor for target genes of other bHLH TFs that enhance transcription. Considering these facts together, we can conclude that ASCL3 is involved in the regulation of the set of the genes downregulated in tame rats, acting as a repressor. Moreover, the large number and variety of bHLH TF binding motifs enriched in the set of genes downregulated in tame rats (Figure 2B) suggest that ASCL3 may repress target genes of various TFs of the bHLH class. The proposed molecular mechanism for changes in gene expression between the tame and aggressive rats due to differences in ASCL3 expression is shown in Figure 3.

3. Discussion

3.1. Genes Encoding Transcription Factors Ascl3, Fos, and Fosb and Gene of Transcription Regulator Arc

Below we consider TF genes Ascl3, Fos, and Fosb and the gene of a transcription regulator, Arc, in detail. The presence of these four genes among the DEGs in the rats selected for their reaction to humans is of considerable importance because the products of these genes, by definition, regulate the expression of a considerable number of genes. Significantly upregulated TFs in aggressive rats include two genes, Fos and Fosb, which are fairly well studied, as well as a gene of an important transcription regulator, Arc. They belong to the group of immediate early genes and are well-known markers of neuron activation [33,34,35,36]. Tight regulation of the expression of immediate early genes is crucial for synaptic plasticity, learning, and memory [37,38]. At baseline, these genes are expressed weakly. Biological activities within neurons lead to a rapid increase in their expression and a subsequent return to baseline levels within a few hours [39]. Nevertheless, what exactly leads to the rapid upregulation of these genes is still unclear; immediate early genes, in particular Fos, Fosb, and Arc, show different temporal patterns of activation after exposure to stress [40,41]. FOSB is a key TF causing changes in gene expression in medium spiny neurons of the nucleus accumbens; this group of neurons in the ventral striatum is an important component of the mesolimbic pathway involved in the reward system and in the formation of pleasure, laughter, addiction, aggression, fear, and the placebo effect [42]. The Arc gene does not encode a classic TF but can localize to the nucleus, thereby regulating gene transcription [43], is involved in various neuronal signaling pathways [44,45], and controls network stability [46]. Arc is a conserved gene in vertebrates [47] and is predominantly expressed in glutamatergic neurons of the cortex and hippocampus. ARC is the main regulator of the expression of a substantial number of genes, thus controlling mRNA levels of more than 1900 genes involved in neuronal plasticity and intrinsic excitability, including many TFs [36].
We investigated the next level of possible regulation of TFs owing to Arc downregulation in the brain tissue samples of tame rats (Table 2) and its relation to the results of motif enrichment analysis (Figure 2). For instance, Leung et al. [36] have provided a list of 40 transcriptional regulators or TFs, the mRNA levels of which are altered when activity-dependent Arc expression is prevented. These proteins are of neuronal relevance and were designated as direct targets for the ARC protein in [36]. The Hocomoco database [28] indicates that only 11 of these 40 transcriptional regulators are well-studied TFs because they possess the motif models in the murine core collection. We next performed a detailed analysis of 6 of these 11 TFs found to be direct targets of ARC: Atf3, Fosb, Jun, Neurod2 (Ndf2), Nfil3, and Pou2f2. Notably, all of them, except Pou2f2, proved to be activated by ARC, whereas Pou2f2 is repressed [36]. The POU2F2 TF belongs to the POU2 (Oct-1/2-like factor) family {3.1.10.2}; this TF is underexpressed and has an intermediate expression pattern between tissue-specific and ubiquitous [48]. POU2F2 regulates the distribution of neurons and inhibits gene expression in neuronal cells [32]. Motifs for POU2F2 and for its closest homolog, TF POU2F1, are similar (p-value < 1 × 10−3, [29]); these are the top-two ranking motifs in the list of enriched motifs for the genes upregulated in the brain tissue samples of tame rats (Figure 2A, ESDEG tool, (log10(padj) > 72 and (log10(padj) > 6).
Although we failed to detect differential expression of the Pou2f2 gene, the known fact that the ARC protein [36], the gene of which is included in the list of repressed genes in tame rats (Table 2), can repress the Pou2f2 gene may explain the significant enrichment of Pou2f1/Pou2f2 motifs in upregulated genes in tame rats. We suggest that this observation deserves further investigation.
The FOSB TF belongs to the Fos-related factor family {1.1.2}. TFs from this family can dimerize with TFs from the closely related family of Jun-related factors {1.1.1}. TFs from the JUN and FOS families form the TF complex known as AP-1, and FOSB regulates cell proliferation, differentiation, and transformation [49]. Two members of the JUN and FOS families, TFs JUN and ATF3, were also detected in the list of 11 direct targets of ARC [36] and are annotated in the Hocomoco database. FOSB, JUN, and ATF3 are known to be expressed ubiquitously [48]. Among these three TF genes under study, only Fosb was significantly downregulated in the brain tissue samples of tame rats (Table 2), and we noted that the FOSL2_MOUSE.H11MO.0.A motif, which shares a high similarity with murine motifs for FOSB, JUN, and ATF3 (p-value < 1 × 10−9, [29]), is enriched in the promoter set of the genes regarded as upregulated in the brain tissue samples of tame rats at the mild threshold (fold change > 1.5, Figure S1). Neurod2 (NDF2) is a highly tissue-specific TF that is expressed only brain tissue [48] and belongs to the Tal-related factor family {1.2.3} of the bHLH TF class. This TF is implicated in neuronal determination and is essential for the repression of the genetic program for neuronal differentiation. Although NDF2 is not differentially expressed here, motifs specific to this TF are enriched in the list of promoters of the genes downregulated in the MT tissue samples of tame rats (see Figure 2B). NFIL3 is a ubiquitously expressed TF [48] from the C/EBP-related family {1.1.8}; it almost satisfied the criteria for upregulated DEGs in tame rats (PADJ = 0.085 < 0.1 and log2 fold = 0.71 < 1, i.e., fold change = 1.64, see Table S1). In conclusion, we documented significant enrichment with DNA motifs specific to a number of TFs, the genes of which are direct targets of ARC, but failed to detect their differential expression. The latter result may be due to the highly dynamic nature of immediate early genes and may become a subject of our further research.

3.2. The Effect of the SNP in the TATA Box of the Ascl3 Gene

In tame rats compared to aggressive rats, the Ascl3 gene is statistically significantly upregulated not only in the MT but also in the hypothalamus [5] and hippocampus [6], as well as in the periaqueductal gray (our unpublished data). We analyzed the initial stage of human ASCL3 gene transcription, which is initiated by the interaction of TATA-binding protein (TBP) with the gene promoter [50]. In version 151 of the dbSNP database [51], 18 single-nucleotide polymorphisms (SNPs) were found in the promoter region within the potential TBP-binding site. Using the SNP_TATA_Z-tester Web service [52], we made a prediction about the change in gene expression (as presented in refs. [53,54,55,56,57,58,59,60,61,62,63,64,65,66,67]; see Supplementary Materials, Figure S2 and Table S2). Among the 18 SNPs, two candidates were identified that statistically significantly change the affinity of TBP for this promoter, as described in detail in our previous work [6]. We tested this prediction for one SNP marker via an electrophoretic mobility shift assay (Supplementary Materials, Figure S3). An in vitro experiment regarding the influence of this SNP on the TBP–DNA interaction revealed that for the human ASCL3 gene, the affinity of TBP weakens by 40% (KD = 12 and 17 nM, respectively) as a result of the T/c substitution (rs1049743008, where T and c are the major and minor alleles, respectively) in the potential TBP-binding site. The Supplementary Materials include ASCL3 electropherograms that illustrate this finding. We can conclude that normal expression of the ASCL3 gene for humans corresponds to overexpression of this gene in tame rats, whereas the minor allele corresponds to aggressive rats, which is consistent with the mainstream understanding of the role of self-domestication in human evolution [68].

3.3. ASCL3 as an Important Regulator of Changes in Behavior through Domestication of Gray Rats

As mentioned above, the bHLH class is widely involved in cell differentiation and tissue development, including neural cell development and differentiation [18]. Our results show that ASCL3 has a major impact on these processes by reducing the expression of these genes in tame rats compared to aggressive rats.
The most important representative of the bHLH class of TFs regulating neural development is the achaete-scute complex-like (ASCL) transcription factor. This observation makes the latter proteins prime candidates for TFs, the target genes of which can be repressed by TF ASCL3. ASCL TFs are directly involved in neural development and neurogenesis [19,69]. In particular, TF ASCL1 is transiently expressed during nervous system development (including the development of the sense of smell and the autonomic nervous system [20]) and controls early and late phases of neurogenesis, the division of radial glia progenitor cells, and the migration of postmitotic neurons [70,71]. Under certain conditions, Ascl1 overexpression can cause dedifferentiation of adult hippocampal neurons into oligodendrocytes [72] and effectively induce the conversion of postnatal astrocytes into functional synapse-forming neurons in the dorsal midbrain of mice [73]. TF ASCL2 is specifically expressed in neuronal progenitors [69,70], whereas ASCL4 is thought to determine the differentiation pathway of nerve and epidermal cells by acting as a proneural protein [74].
Neurogenesis, as a process of differentiation of a neuronal progenitor and its integration into a neuronal network, is most active in the prenatal period; however, it continues in the adult organism, albeit less intensively [75]. In adult mammals, new neurons constantly form in two brain regions: the subgranular zone of the hippocampus and the subependymal zone of lateral ventricles. Therefore, the overexpression of the Ascl3 gene in the hippocampus of tame rats [6] could be expected based on our assumptions, although its significant underexpression in aggressive animals is difficult to interpret. In this work, we detected Ascl3 upregulation in the MT, raising the possibility of neurogenesis in the midbrain. Data reported in the literature on this subject are rather contradictory [76,77,78] and do not provide a definitive answer. Neurogenesis in the MT region, just as in other non-neurogenic areas of the brain, probably indicates immaturity of the organism and ongoing formation of neurons.
Therefore, we can theorize that in tame rats, neurogenesis is active for a longer period, continuing into the second month of life. On the other hand, TF ASCL3, by repressing the expression of genes that are responsible for neurogenesis, affects the rate of maturation. Accordingly, the genes whose expression levels change under the influence of TF ASCL3 become differentially expressed and determine the phenotype of tame behavior.
If our theory about prolonged activity of neurogenesis in 2-month-old tame rats owing to Ascl3 overexpression is correct, then it is possible that the puberty period in these animals is also delayed compared to that in aggressive rats because the neuroendocrine system directly controls the reproductive function. A study on the effect of long-term selection for behavior on the reproductive system activity of tame and aggressive rats revealed that the testosterone level and relative weights of the testes and testicular appendages are significantly lower in male tame rats than in aggressive rats at 2 months of age [79]. In this regard, the positive correlation between aggressiveness and testosterone levels is relevant, the latter being significantly higher in 2-month-old aggressive male gray rats in comparison with not only tame males but also unselected males [79].
It is believed that delayed neurogenesis is an indicator of neoteny [80]. In tame animals, the maturation period is longer than in wild animals; for example, in a research article on another model animal, the silver fox, it was demonstrated that prolonged artificial selection leads to increased neurogenesis in the hippocampus of tame foxes and that this process is inversely proportional to increased aggressiveness toward humans [81]. Taken together, the results with respect to prolonged neurogenesis and delayed puberty in tame rats (which are often used as an experimental model of behavior of domesticated animals) suggest that the manifestation of neoteny is a characteristic trait of animals during domestication [82,83,84,85].
Therefore, in tame and aggressive rats, a new function of TF ASCL3 was demonstrated; by changing the expression of gene networks, ASCL3 may play an important role in changes in the behavior of animals during domestication.

4. Materials and Methods

4.1. Experimental Animals

This study was conducted on adult male gray rats (R. norvegicus) artificially bred for more than 90 generations for either aggressive or tame behavior (as two outbred strains). The rats were kept under standard conditions at the Conventional Animal Facility at the ICG SB RAS (Novosibirsk, Russia) as described elsewhere [14,86,87]. A total of 22 rats (11 aggressive and 11 tame) were included in this study, each 2 months old and weighing 250–270 g, all from different unrelated litters. All the rats were decapitated. Using a handbook technique [88], we excised samples of the MT, which were then flash-frozen in liquid nitrogen and stored at −70 °C until use. Every effort was made to minimize the number of animals under study and to prevent their suffering. This work was conducted in accordance with the guidelines of the Declaration of Helsinki, Directive 2010/63/EU of the European Parliament, and of the European Council resolution of 22 September 2010.
The research protocol was approved by the Interinstitutional Commission on Bioethics at the ICG SB RAS, Novosibirsk, Russia (approval documentation No. 97 dated October 2021).

4.2. Total RNA Extraction, Library Construction, and RNA-Seq Data Analysis

Library preparation, sequencing, and bioinformatics manipulations can be found elsewhere [5,6]. We used the following criteria to determine whether a gene was differentially expressed: statistical significance at a p-value of <0.1 adjusted for multiple comparisons and a fold change of gene expression of >2.0 (or a fold change of >1.5).

4.3. qPCR

To selectively and independently verify the tame-versus-aggressive rat hippocampal DEGs found in the present study (Table 2), we performed a qPCR control assay on the total RNA extracted only from the remaining samples of the hypothalamus of tame (n = 8) and aggressive (n = 8) rats. First, with the help of TRIzol™, we isolated total RNA, purified it on Agencourt RNAClean XP Kit magnetic beads (Beckman, #A63987), and quantified it by means of a Qubit™ 2.0 fluorometer (Invitrogen/Life Technologies), along with a high-sensitivity RNA kit (Invitrogen, cat. #Q32852). Then, we synthesized cDNA using a reverse transcription kit (Syntol, #OT-1). Next, using the PrimerBLAST web service [89], we designed oligonucleotide primers for qPCR (Table 3). Then, we carried out qPCR on a LightCycler® 96 (Roche, Basel, Basel-Stadt, Switzerland) with an EVA Green I kit in three technical replicates. We determined qPCR efficiency by means of serial cDNA dilutions (standards). In line with the commonly accepted recommendations [90], we simultaneously analyzed two reference genes, namely Ppia (peptidylprolyl isomerase A) [91] and Rpl30 (ribosomal protein L30) [92].

4.4. Motif Enrichment Analysis

We performed this analysis on promoters of the genes upregulated and downregulated in the brain tissue samples of tame rats. For this purpose, we developed the ESDEG tool, which is available at https://github.com/ubercomrade/esdeg (accessed on 12 October 2022). The full set of 18,442 promoter regions (−2000; +1) relative to the transcription start sites for protein-coding genes annotated in an R. norvegicus genome (RGSC version Rnor_6.0, UCSC, Santa Cruz, USA, version Rn6, July 2014 assembly) were employed for the analysis; alternative transcripts were disregarded. The murine core collection including 356 motifs (nucleotide frequency matrices) from the HOCOMOCO database [28] served as a list of motifs, of which only 351 motifs were included in the analysis according to a filtering approach proposed in [29]. We decided to use the murine motif collection because this species is the best-studied model organism closest to the rat with the most complete motif collection accumulated to date.
To identify overrepresented motifs, we first performed motif recognition in the whole set of 18,442 promoter regions mentioned above, utilizing each from a collection of 351 nucleotide frequency matrices. Next, a table of recognition thresholds was compiled as described in [29], where expected recognition rates (ERRs) for each threshold were calculated as probabilities of site prediction in the whole-genome dataset of promoters. Given that the enrichment of sites depends on a chosen threshold, subsequent calculations were performed for 30 recognition thresholds (equidistant on a logarithmic scale of ERRs) in the range of 5 × 10−4 to 1 × 10−4. Next, all promoters were categorized into either a foreground or background group. The foreground group contained the promoters of the significantly upregulated or downregulated genes. We utilized the following criteria to determine whether a gene was differentially expressed: an adjusted p-value below 0.1 and a fold change above 2.0 (alternatively: above 1.5). The background group included the same number of promoters as the foreground group, but they were randomly selected according to the following criteria: an adjusted p-value above 0.1 and a fold change within the range of 0.80 to 1.25. Next, we calculated the frequency of a motif in each of the gene promoter groups. To estimate the significance of the motif enrichment, we applied the Monte Carlo approach. For a given motif and a recognition threshold, site frequency was calculated for the foreground group (AVFOR), whereas a background group was generated many times to estimate the site frequency distribution. Next, the average (AVBACK) and standard deviation (SDBACK) for the frequencies in the background group set were calculated. Then, a Z-score was calculated according to the formula (AVFOR − AVBACK)/SDBACK. This Z-score allowed us to compute the significance of enrichment (p-value) by fitting a site’s frequency distribution to a normal distribution. To confirm enrichment for a given motif, we combined multiple p-values, referring all recognition thresholds into a single unified p-value according to the Hartung method [93]. Thus, a unified p-value was calculated for each motif. Finally, we calculated an adjusted p-value (padj) from the set of unified p-values for all motifs with the Benjamini–Hochberg correction for multiple comparisons; in this way, our analysis involved multiple simultaneous statistical tests for various motifs.

4.5. Statistical Analysis

We performed the Mann–Whitney U test and Fisher’s Z-test using options in the STATISTICA standard toolbox (StatsoftTM, Hamburg, Germany).

5. Conclusions

We performed RNA-seq to investigate differential expression of genes in tissue samples from the tegmental region of the midbrain of rats showing either tame or aggressive behavior. The list of DEGs includes genes related to extracellular matrix components, transporter proteins, neurotransmitter system genes, signaling molecules, and immune-system proteins, as well as genes encoding transcription factors. A new function of the ASCL3 TF, which is upregulated in the brain tissue samples of tame rats, is proposed. We demonstrated that several structural subtypes of E-box motifs—known as binding sites for many developmental TFs of the bHLH class, including the ASCL subfamily of TFs—are enriched in the set of promoters of the DEGs downregulated in the tissue samples of tame rats. We showed that the increased expression of the Ascl3 gene in tame rats is the most likely cause of the statistically significant downregulation of the corresponding genes, as its product, ASCL3 TF, presumably acts as a repressor. In aggressive rats, the Ascl3 gene is expressed at a low level; therefore, ASCL3 TF does not inhibit the expression of the same genes, so they are upregulated compared to tame rats. Thus, we hypothesize that the upregulation of Ascl3 in tame rats plays an important role in changes of the behavior of animals during domestication, probably promoting longer neurogenesis (as compared to aggressive rats), which is a sign of neoteny and domestication.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijms232012269/s1.

Author Contributions

Conceptualization, D.O.; supervision, A.L.M.; writing of the manuscript, I.C., V.G.L. and D.O.; investigation, I.C., N.V.K., R.K. and S.S.; development of the motif enrichment analysis methodology, V.G.L. and A.T.; validation, E.S. and L.S.; data curation, I.C. and D.O. All authors have read and agreed to the published version of the manuscript.

Funding

The study was supported by the Russian Science Foundation (project 21-14-00240).

Institutional Review Board Statement

This work was conducted in accordance with the guidelines of the Declaration of Helsinki, Directive 2010/63/EU of the European Parliament, and of the European Council resolution of 22 September 2010. The research protocol was approved by the Interinstitutional Commission on Bioethics at the ICG SB RAS, Novosibirsk, Russia (approval documentation No. 97 dated October 2021).

Informed Consent Statement

Not applicable.

Data Availability Statement

The primary RNA-Seq data obtained in this work were deposited in the NCBI SRA database (ID = PRJNA668014).

Acknowledgments

The authors are grateful to Nikolai Shevchuk for fruitful discussions of the work and assistance with adapting it to the requirements on the English language for scientific publications. The authors are also thankful to the Multi-Access Center “Bioinformatics” for the use of computational resources as supported by Russian government project FWNR-2022-0020 and the Russian Federal Science and Technology Program for the Development of Genetic Technologies.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

ASCLachaete-scute complex like
bHLHbasic helix–loop–helix
DEGdifferentially expressed gene
MTtegmental region of the midbrain
PAGperiaqueductal gray
qPCRquantitative PCR
SNPsingle-nucleotide polymorphism
TFtranscription factor
TBPTATA-binding protein

References

  1. Rubin, C.J.; Zody, M.; Eriksson, J.; Meadows, J.R.S.; Sherwood, E.; Webster, M.T.; Jiang, L.; Ingman, M.; Sharpe, T.; Ka, S.; et al. Whole-genome resequencing reveals loci under selection during chicken domestication. Nature 2010, 464, 587–591. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Axelsson, E.; Ratnakumar, A.; Arendt, M.-L.; Maqbool, K.; Webster, M.T.; Perloski, M.; Liberg, O.; Arnemo, J.M.; Hedhammar, A.; Lindblad-Toh, K. The genomic signature of dog domestication reveals adaptation to a starch-rich diet. Nature 2013, 495, 360–364. [Google Scholar] [CrossRef] [PubMed]
  3. Carneiro, M.; Rubin, C.-J.; Di Palma, F.; Albert, F.W.; Alföldi, J.; Barrio, A.M.; Pielberg, G.; Rafati, N.; Sayyab, S.; Turner-Maier, J.; et al. Rabbit genome analysis reveals a polygenic basis for phenotypic change during domestication. Science 2014, 345, 1074–1079. [Google Scholar] [CrossRef] [Green Version]
  4. Rosenfeld, C.S.; Hekman, J.P.; Johnson, J.L.; Lyu, Z.; Ortega, M.T.; Joshi, T.; Mao, J.; Vladimirova, A.V.; Gulevich, R.G.; Kharlamova, A.V.; et al. Hypothalamic transcriptome of tame and aggressive silver foxes (Vulpes vulpes) identifies gene expression differences shared across brain regions. Genes Brain Behav. 2020, 19, e12614. [Google Scholar] [CrossRef] [PubMed]
  5. Chadaeva, I.; Ponomarenko, P.; Kozhemyakina, R.; Suslov, V.; Bogomolov, A.; Klimova, N.; Shikhevich, S.; Savinkova, L.; Oshchepkov, D.; Kolchanov, N.A.; et al. Domestication explains two-thirds of differential-gene-expression variance between domestic and wild animals, the remaining one-third reflects intraspecific and interspecific variation. Animals 2021, 11, 2667. [Google Scholar] [CrossRef]
  6. Oshchepkov, D.; Chadaeva, I.; Kozhemyakina, R.; Zolotareva, K.; Khandaev, B.; Sharypova, E.; Ponomarenko, P.; Bogomolov, A.; Klimova, N.V.; Shikhevich, S.; et al. Stress reactivity, susceptibility to hypertension, and differential expression of genes in hypertensive compared to normotensive patients. Int. J. Mol. Sci. 2022, 23, 2835. [Google Scholar] [CrossRef]
  7. Wilkins, A.S.; Wrangham, R.W.; Fitch, W.T. The “domestication syndrome” in mammals: A unified explanation based on neural crest cell behavior and genetics. Genetics 2014, 197, 795–808. [Google Scholar] [CrossRef] [Green Version]
  8. Kruska, D. The effect of domestication of brain size and composition in the mink (Mustela vison). J. Zool. 1996, 239, 645–661. [Google Scholar] [CrossRef]
  9. Jensen, P. Domestication—From behaviour to genes and back again. Appl. Anim. Behav. Sci. 2006, 97, 3–15. [Google Scholar] [CrossRef]
  10. Setchell, B.P. Domestication and reproduction. Anim. Reprod. Sci. 1992, 28, 195–202. [Google Scholar] [CrossRef]
  11. Belyaev, D.K.; Plyusnina, I.Z.; Trut, L.N. Domestication in the silver fox (Vulpes fulvus desm): Changes in physiological boundaries of the sensitive period of primary socialization. Appl. Anim. Behav. Sci. 1984, 13, 359–370. [Google Scholar] [CrossRef]
  12. Trut, L.; Oskina, I.; Kharlamova, A. Animal evolution during domestication: The domesticated fox as a model. Bioessays 2009, 31, 349–360. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Belyaev, D.K.; Borodin, P.M. The influence of stress on variation and its role in evolution. Biol. Zent. 1982, 101, 705–714. [Google Scholar]
  14. Plyusnina, I.; Oskina, I. Behavioral and adrenocortical responses to open-field test in rats selected for reduced aggressiveness toward humans. Physiol. Behav. 1997, 61, 381–385. [Google Scholar] [CrossRef]
  15. Hrabovszky, E.; Halász, J.; Meelis, W.; Kruk, M.R.; Liposits, Z.; Haller, J. Neurochemical characterization of hypothalamic neurons involved in attack behavior: Glutamatergic dominance and co-expression of thyrotropin-releasing hormone in a subset of glutamatergic neurons. Neuroscience 2005, 133, 657–666. [Google Scholar] [CrossRef] [PubMed]
  16. Ito, R.; Lee, A.C.H. The role of the hippocampus in approach-avoidance conflict decision-making: Evidence from rodent and human studies. Behav. Brain Res. 2016, 313, 345–357. [Google Scholar] [CrossRef] [PubMed]
  17. Wingender, E.; Schoeps, T.; Dönitz, J. TFClass: An expandable hierarchical classification of human transcription factors. Nucleic Acids Res. 2013, 41, D165–D170. [Google Scholar] [CrossRef] [PubMed]
  18. Torres-Machorro, A.L. Homodimeric and heterodimeric interactions among vertebrate basic helix-loop-helix transcription factors. Int. J. Mol. Sci. 2021, 22, 12855. [Google Scholar] [CrossRef]
  19. Sanes, D.H. Development of the Nervous System, 3rd ed.; Sanes, D., Reh, T., Harris, W., Eds.; Academic Press: Cambridge, MA, USA, 2012. [Google Scholar] [CrossRef]
  20. Guillemot, F.; Lo, L.-C.; Johnson, J.E.; Auerbach, A.; Anderson, D.J.; Joyner, A.L. Mammalian achaete-scute homolog 1 is required for the early development of olfactory and autonomic neurons. Cell 1993, 75, 463–476. [Google Scholar] [CrossRef]
  21. Gyoja, F. Basic helix-loop-helix transcription factors in evolution: Roles in development of mesoderm and neural tissues. Genesis 2017, 55, e23051. [Google Scholar] [CrossRef] [PubMed]
  22. Niu, Y.; Moghimyfiroozabad, S.; Moghimyfiroozabad, A.; Tierney, T.S.; Alavian, K.N. The factors for the early and late development of midbrain dopaminergic neurons segregate into two distinct evolutionary clusters. Brain Disord. 2021, 1, 100002. [Google Scholar] [CrossRef]
  23. Shaikh, M.B.; Brutus, M.; Siegel, H.E.; Siegel, A. Differential control of aggression by the midbrain. Exp. Neurol. 1984, 83, 436–442. [Google Scholar] [CrossRef]
  24. Coimbra, B.; Soares-Cunha, C.; Vasconcelos, N.; Domingues, A.V.; Borges, S.; Sousa, N.; Rodrigues, A.J. Role of laterodorsal tegmentum projections to nucleus accumbens in reward-related behaviors. Nat. Commun. 2019, 10, 4138. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. 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]
  26. Mahadevia, D.; Saha, R.; Manganaro, A.; Chuhma, N.; Ziolkowski-Blake, A.; Morgan, A.A.; Dumitriu, D.; Rayport, S.; Ansorge, M.S. Dopamine promotes aggression in mice via ventral tegmental area to lateral septum projections. Nat. Commun. 2021, 12, 6796. [Google Scholar] [CrossRef]
  27. Sherman, B.T.; Hao, M.; Qiu, J.; Jiao, X.; Baseler, M.W.; Lane, H.C.; Imamichi, T.; Chang, W. DAVID: A web server for functional enrichment analysis and functional annotation of gene lists (2021 update). Nucleic Acids Res. 2022, 50, W216–W221. [Google Scholar] [CrossRef] [PubMed]
  28. Kulakovskiy, I.V.; Vorontsov, I.E.; Yevshin, I.S.; Sharipov, R.N.; Fedorova, A.D.; Rumynskiy, E.I.; Medvedeva, Y.A.; Magana-Mora, A.; Bajic, V.B.; Papatsenko, D.A.; et al. HOCOMOCO: Towards a complete collection of transcription factor binding models for human and mouse via large-scale ChIP-Seq analysis. Nucleic Acids Res. 2018, 46, D252–D259. [Google Scholar] [CrossRef] [PubMed]
  29. Levitsky, V.; Zemlyanskaya, E.; Oshchepkov, D.; Podkolodnaya, O.; Ignatieva, E.; Grosse, I.; Mironova, V.; Merkulova, T. A single ChIP-seq dataset is sufficient for comprehensive analysis of motifs co-occurrence with MCOT package. Nucleic Acids Res. 2019, 47, e139. [Google Scholar] [CrossRef] [PubMed]
  30. Ambrosini, G.; Vorontsov, I.; Penzar, D.; Groux, R.; Fornes, O.; Nikolaeva, D.D.; Ballester, B.; Grau, J.; Grosse, I.; Makeev, V.; et al. Insights gained from a comprehensive all-against-all transcription factor binding motif benchmarking study. Genome Biol. 2020, 21, 114. [Google Scholar] [CrossRef] [PubMed]
  31. Yoshida, S.; Ohbo, K.; Takakura, A.; Takebayashi, H.; Okada, T.; Abe, K.; Nabeshima, Y. Sgn1, a basic helix-loop-helix transcription factor delineates the salivary gland duct cell lineage in mice. Dev. Biol. 2001, 240, 517–530. [Google Scholar] [CrossRef] [PubMed]
  32. Harris, A.; Masgutova, G.; Collin, A.; Toch, M.; Hidalgo-Figueroa, M.; Jacob, B.; Corcoran, L.M.; Francius, C.; Clotman, F. Onecut factors and Pou2f2 regulate the distribution of V2 interneurons in the mouse developing spinal cord. Front. Cell. Neurosci. 2019, 13, 184. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Dragunow, M.; Faull, R. The use of c-fos as a metabolic marker in neuronal pathway tracing. J. Neurosci. Methods 1989, 29, 261–265. [Google Scholar] [CrossRef]
  34. Kovács, K.J. Measurement of immediate-early gene activation-c-fos and beyond. J. Neuroendocrinol. 2008, 20, 665–672. [Google Scholar] [CrossRef]
  35. Hamann, C.S.; Bankmann, J.; Mora Maza, H.; Kornhuber, J.; Zoicas, I.; Schmitt-Böhrer, A. Social fear affects limbic system neuronal activity and gene expression. Int. J. Mol. Sci. 2022, 23, 8228. [Google Scholar] [CrossRef]
  36. Leung, H.-W.; Foo, G.; VanDongen, A. Arc regulates transcription of genes for plasticity; excitability and alzheimer’s disease. Biomedicines 2022, 10, 1946. [Google Scholar] [CrossRef] [PubMed]
  37. Guzowski, J.F.; Timlin, J.A.; Roysam, B.; McNaughton, B.L.; Worley, P.F.; Barnes, C.A. Mapping behaviorally relevant neural circuits with immediate-early gene expression. Curr. Opin. Neurobiol. 2005, 15, 599–606. [Google Scholar] [CrossRef]
  38. Boutros, S.W.; Kessler, K.; Unni, V.K.; Raber, J. Infusion of etoposide in the CA1 disrupts hippocampal immediate early gene expression and hippocampus-dependent learning. Sci. Rep. 2022, 12, 12834. [Google Scholar] [CrossRef]
  39. Lamprecht, R.; Dudai, Y. Transient expression of c-Fos in rat amygdala during training is required for encoding conditioned taste aversion memory. Learn. Mem. 1996, 3, 31–41. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  40. Honkaniemi, J.; Kononen, J.; Kainu, T.; Pyykönen, I.; Pelto-Huikko, M. Induction of multiple immediate early genes in rat hypothalamic paraventricular nucleus after stress. Brain Res. Mol. Brain Res. 1994, 25, 234–241. [Google Scholar] [CrossRef]
  41. Ons, S.; Martí, O.; Armario, A. Stress-induced activation of the immediate early gene Arc (activity-regulated cytoskeleton-associated protein) is restricted to telencephalic areas in the rat brain: Relationship to c-fos mRNA. J. Neurochem. 2004, 89, 1111–1118. [Google Scholar] [CrossRef] [PubMed]
  42. Grueter, B.A.; Robison, A.J.; Neve, R.L.; Nestler, E.J.; Malenka, R.C. ∆FosB differentially modulates nucleus accumbens direct and indirect pathway function. Proc. Natl. Acad. Sci. USA 2013, 110, 1923–1928. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Korb, E.; Wilkinson, C.L.; Delgado, R.N.; Lovero, K.L.; Finkbeiner, S. Arc in the nucleus regulates PML-dependent GluA1 transcription and homeostatic plasticity. Nat. Neurosci. 2013, 16, 874–883. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Alberi, L.; Liu, S.; Wang, Y.; Badie, R.; Smith-Hicks, C.; Wu, J.; Pierfelice, T.J.; Abazyan, B.; Mattson, M.P.; Kuhl, D.; et al. Activity-induced Notch signaling in neurons requires Arc/Arg3.1 and is essential for synaptic plasticity in hippocampal networks. Neuron 2011, 69, 437–444. [Google Scholar] [CrossRef] [Green Version]
  45. Messaoudi, E.; Kanhema, T.; Soulé, J.; Tiron, A.; Dagyte, G.; da Silva, B.; Bramham, C.R. Sustained Arc/Arg3.1 synthesis controls long-term potentiation consolidation through regulation of local actin polymerization in the dentate gyrus in vivo. J. Neurosci. 2007, 27, 10445–10455. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Peebles, C.L.; Yoo, J.; Thwin, M.T.; Palop, J.J.; Noebels, J.L.; Finkbeiner, S. Arc regulates spine morphology and maintains network stability in vivo. Proc. Natl. Acad. Sci. USA 2010, 107, 18173–18178. [Google Scholar] [CrossRef] [Green Version]
  47. Epstein, I.; Finkbeiner, S. The Arc of cognition: Signaling cascades regulating Arc and implications for cognitive function and disease. Semin. Cell Dev. Biol. 2018, 77, 63–72. [Google Scholar] [CrossRef]
  48. Fagerberg, L.; Hallström, B.M.; Oksvold, P.; Kampf, C.; Djureinovic, D.; Odeberg, J.; Habuka, M.; Tahmasebpoor, S.; Danielsson, A.; Edlund, K.; et al. Analysis of the human tissue-specific expression by genome-wide integration of transcriptomics and antibody-based proteomics. Mol. Cell. Proteom. 2014, 13, 397–406. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  49. Hess, J.; Angel, P.; Schorpp-Kistner, M. AP-1 subunits: Quarrel and harmony among siblings. J. Cell Sci. 2004, 117 Pt 25, 5965–5973. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  50. Sharypova, E.B.; Drachkova, I.A.; Chadaeva, I.V.; Ponomarenko, M.P.; Savinkova, M.P. An experimental study of the effects of SNPs in the TATA boxes of the GRIN1; ASCL3 and NOS1 genes on interactions with the TATA-binding protein. Vavilovskii Zhurnal Genet. Sel. 2022, 26, 227–233. [Google Scholar] [CrossRef]
  51. Day, I.N. dbSNP in the detail and copy number complexities. Hum. Mutat. 2010, 31, 2–4. [Google Scholar] [CrossRef] [PubMed]
  52. Ponomarenko, M.; Rasskazov, D.; Arkova, O.; Ponomarenko, P.; Suslov, V.; Savinkova, L.; Kolchanov, N. How to use SNP_TATA_Comparator to find a significant change in gene expression caused by the regulatory SNP of this gene’s promoter via a change in affinity of the TATA-binding protein for this promoter. Biomed. Res. Int. 2015, 2015, 359835. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Ponomarenko, M.; Kleshchev, M.; Ponomarenko, P.; Chadaeva, I.; Sharypova, E.; Rasskazov, D.; Kolmykov, S.; Drachkova, I.; Vasiliev, G.; Gutorova, N.; et al. Disruptive natural selection by male reproductive potential prevents underexpression of protein-coding genes on the human Y chromosome as a self-domestication syndrome. BMC Genet. 2020, 21 (Suppl. S1), 89. [Google Scholar] [CrossRef] [PubMed]
  54. Haeussler, M.; Raney, B.; Hinrichs, A.; Clawson, H.; Zweig, A.; Karolchik, D.; Casper, J.; Speir, M.; Haussler, D.; Kent, W. Navigating protected genomics data with UCSC Genome Browser in a box. Bioinformatics 2015, 31, 764–766. [Google Scholar] [CrossRef] [Green Version]
  55. Ponomarenko, P.; Savinkova, L.; Drachkova, I.; Lysova, M.; Arshinova, T.; Ponomarenko, M.; Kolchanov, N. A step-by-step model of TBP/TATA box binding allows predicting human hereditary diseases by single nucleotide polymorphism. Dokl. Biochem. Biophys. 2008, 419, 88–92. [Google Scholar] [CrossRef] [PubMed]
  56. Delgadillo, R.; Whittington, J.; Parkhurst, L.; Parkhurst, L. The TBP core domain in solution variably bends TATA sequences via a three-step binding mechanism. Biochemistry 2009, 48, 1801–1809. [Google Scholar] [CrossRef] [PubMed]
  57. Coleman, R.A.; Pugh, B.F. Evidence for functional binding and stable sliding of the TATA binding protein on nonspecific DNA. J. Biol. Chem. 1995, 270, 13850–13859. [Google Scholar] [CrossRef] [Green Version]
  58. Berg, O.G.; von Hippel, P.H. Selection of DNA binding sites by regulatory proteins: Statistical-mechanical theory and application to operators and promoters. J. Mol. Biol. 1987, 193, 723–750. [Google Scholar] [CrossRef]
  59. Flatters, D.; Lavery, R. Sequence-dependent dynamics of TATA-Box binding sites. Biophys. J. 1998, 75, 372–381. [Google Scholar] [CrossRef] [Green Version]
  60. Hahn, S.; Buratowski, S.; Sharp, P.; Guarente, L. Yeast TATA-binding protein TFIID binds to TATA elements with both consensus and nonconsensus DNA sequences. Proc. Natl. Acad. Sci. USA 1989, 86, 5718–5722. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  61. Bucher, P. Weight matrix descriptions of four eukaryotic RNA polymerase II promoter elements derived from 502 unrelated promoter sequences. J. Mol. Biol. 1990, 212, 563–578. [Google Scholar] [CrossRef]
  62. Karas, H.; Knuppel, R.; Schulz, W.; Sklenar, H.; Wingender, E. Combining structural analysis of DNA with search routines for the detection of transcription regulatory elements. Comput. Applic. Biosci. 1996, 12, 441–446. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  63. Ponomarenko, M.P.; Ponomarenko, J.V.; Frolov, A.S.; Podkolodny, N.L.; Savinkova, L.K.; Kolchanov, N.A.; Overton, G.C. Identification of sequence-dependent features correlating to activity of DNA sites interacting with proteins. Bioinformatics 1999, 15, 687–703. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  64. Waardenberg, A.; Basset, S.; Bouveret, R.; Harvey, R. CompGO: An R package for comparing and visualizing Gene Ontology enrichment differences between DNA binding experiments. BMC Bioinform. 2015, 16, 275. [Google Scholar] [CrossRef] [Green Version]
  65. Pugh, B. Purification of the human TATA-binding protein, TBP. Methods Mol. Biol. 1995, 37, 359–367. [Google Scholar] [CrossRef]
  66. Savinkova, L.; Drachkova, I.; Arshinova, T.; Ponomarenko, P.; Ponomarenko, M.; Kolchanov, N. An experimental verification of the predicted effects of promoter TATA-box polymorphisms associated with human diseases on interactions between the TATA boxes and TATA-binding protein. PLoS ONE 2013, 8, e54626. [Google Scholar] [CrossRef] [Green Version]
  67. Drachkova, I.; Savinkova, L.; Arshinova, T.; Ponomarenko, M.; Peltek, S.; Kolchanov, N. The mechanism by which TATA-box polymorphisms associated with human hereditary diseases influence interactions with the TATA-binding protein. Hum. Mutat. 2014, 35, 601–608. [Google Scholar] [CrossRef]
  68. Herbeck, Y.E.; Zakharov, I.K.; Trapezov, O.V.; Shumny, V.K. Evolution compressed in time. Philos. Sci. 2013, 1, 115–139. [Google Scholar]
  69. Liu, Z.; Wang, X.; Jiang, K.; Ji, X.; Zhang, Y.A.; Chen, Z. TNFα-induced up-regulation of Ascl2 affects the differentiation and proliferation of neural stem cells. Aging Dis. 2019, 10, 1207–1220. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  70. Johnson, J.E.; Birren, S.J.; Anderson, D.J. Two rat homologues of Drosophila achaete-scute specifically expressed in neuronal precursors. Nature 1990, 346, 858–861. [Google Scholar] [CrossRef]
  71. Pacary, E.; Heng, J.L.; Azzarelli, R.; Riou, P.; Castro, D.; Lebel-Potter, M.; Parras, C.; Bell, D.M.; Ridley, A.J.; Parsons, M.; et al. Proneural transcription factors regulate different steps of cortical neuron migration through Rnd-mediated inhibition of RhoA signaling. Neuron 2011, 69, 1069–1084. [Google Scholar] [CrossRef] [Green Version]
  72. Jessberger, S.; Gage, F.H. Fate plasticity of adult hippocampal progenitors: Biological relevance and therapeutic use. Trends Pharmacol. Sci. 2009, 30, 61–65. [Google Scholar] [CrossRef] [PubMed]
  73. Liu, Y.; Miao, Q.; Yuan, J.; Han, S.; Zhang, P.; Li, S.; Rao, Z.; Zhao, W.; Ye, Q.; Geng, J.; et al. Ascl1 converts dorsal midbrain astrocytes into functional neurons in vivo. J. Neurosci. 2015, 35, 9336–9355. [Google Scholar] [CrossRef] [Green Version]
  74. Henrique, D.; Tyler, D.; Kintner, C.; Heath, J.K.; Lewis, J.H.; Ish-Horowicz, D.; Storey, K.G. Cash4, a novel achaete-scute homolog induced by Hensen’s node during generation of the posterior nervous system. Genes Dev. 1997, 11, 603–615. [Google Scholar] [CrossRef] [Green Version]
  75. Kempermann, G. The pessimist’s and optimist’s views of adult neurogenesis. Cell 2011, 145, 1009–1011. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  76. Farzanehfar, P.; Lu, S.S.; Dey, A.; Musiienko, D.; Baagil, H.; Horne, M.K.; Aumann, T.D. Evidence of functional duplicity of Nestin expression in the adult mouse midbrain. Stem Cell Res. 2017, 19, 82–93. [Google Scholar] [CrossRef] [PubMed]
  77. Farzanehfar, P.; Horne, M.K.; Aumann, T.D. Can valproic acid regulate neurogenesis from Nestin+ cells in the adult midbrain? Neurochem. Res. 2017, 42, 2127–2134. [Google Scholar] [CrossRef]
  78. Fauser, M.; Pan-Montojo, F.; Richter, C.; Kahle, P.J.; Schwarz, S.C.; Schwarz, J.; Storch, A.; Hermann, A. Chronic-progressive dopaminergic deficiency does not induce midbrain neurogenesis. Cells 2021, 10, 775. [Google Scholar] [CrossRef]
  79. Prasolova, L.A.; Gerbek, Y.E.; Gulevich, R.G.; Shikhevich, S.G.; Konoshenko, M.Y.; Kozhemyakina, R.V.; Oskina, I.N.; Plyusnina, I.Z. The effects of prolonged selection for behavior on the stress response and activity of the reproductive system of male grey mice (Rattus norvegicus) Russ. J. Genet. 2014, 50, 959–966. [Google Scholar] [CrossRef]
  80. Penz, O.K.; Fuzik, J.; Kurek, A.B.; Romanov, R.; Larson, J.; Park, T.J.; Harkany, T.; Keimpema, E. Protracted brain development in a rodent model of extreme longevity. Sci. Rep. 2015, 5, 11592. [Google Scholar] [CrossRef] [Green Version]
  81. Huang, S.; Slomianka, L.; Farmer, A.J.; Kharlamova, A.V.; Gulevich, R.G.; Herbeck, Y.E.; Trut, L.N.; Wolfer, D.P.; Amrein, I. Selection for tameness, a key behavioral trait of domestication, increases adult hippocampal neurogenesis in foxes. Hippocampus 2015, 25, 963–975. [Google Scholar] [CrossRef]
  82. Konoshenko, M.Y.; Plyusnina, I.Z. Play fighting in rats, selected for the elimination and the enhancement of aggressiveness towards human. Zhurnal Vyss. Nervn. Deiatelnosti Im. I P Pavlov. 2017, 67, 93–105. [Google Scholar]
  83. Trut, L.N.; Plyusnina, I.Z.; Os’kina, I.N. An experiment on fox domestication and debatable issues of evolution of the dog. Genetika 2004, 40, 794–807. [Google Scholar] [CrossRef] [PubMed]
  84. Driscoll, C.A.; Macdonald, D.W.; O’Brien, S.J. From wild animals to domestic pets; an evolutionary view of domestication. PNAS 2009, 106 (Suppl. S1), 9971–9978. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  85. Price, E.O. Behavioral development in animals undergoing domestication. Appl. Anim. Behav. Sci. 1999, 65, 245–271. [Google Scholar] [CrossRef]
  86. Naumenko, E.; Popova, N.; Nikulina, E.; Dygalo, N.; Shishkina, G.; Borodin, P.; Markel, A. Behavior, adrenocortical activity; and brain monoamines in Norway rats selected for reduced aggressiveness towards man. Pharmacol. Biochem. Behav. 1989, 33, 85–91. [Google Scholar] [CrossRef]
  87. Plyusnina, I.Z.; Solov’eva, M.Y.; Oskina, I.N. Effect of domestication on aggression in gray Norway rats. Behav. Genet. 2011, 41, 583–592. [Google Scholar] [CrossRef] [PubMed]
  88. Paxinos, G.; Watson, C.R. The Rat Brain in Stereotaxic Coordinates, 7th ed.; Academic Press: London, UK, 2013; p. 472. [Google Scholar]
  89. Ye, J.; Coulouris, G.; Zaretskaya, I.; Cutcutache, I.; Rozen, S.; Madden, T.L. Primer-BLAST: A tool to design target-specific primers for polymerase chain reaction. BMC Bioinform. 2012, 13, 134. [Google Scholar] [CrossRef] [Green Version]
  90. Bustin, S.A.; Benes, V.; Garson, J.A.; Hellemans, J.; Huggett, J.; Kubista, M.; Mueller, R.; Nolan, T.; Pfaffl, M.W.; Shipley, G.L.; et al. The MIQE guidelines: Minimum information for publication of quantitative real-time PCR experiments. Clin. Chem. 2009, 55, 611–622. [Google Scholar] [CrossRef] [Green Version]
  91. Gholami, K.; Loh, S.Y.; Salleh, N.; Lam, S.K.; Hoe, S.Z. Selection of suitable endogenous reference genes for qPCR in kidney and hypothalamus of rats under testosterone influence. PLoS ONE 2017, 12, e0176368. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  92. Penning, L.C.; Vrieling, H.E.; Brinkhof, B.; Riemers, F.M.; Rothuizen, J.; Rutteman, G.R.; Hazewinkel, H.A. A validation of 10 feline reference genes for gene expression measurements in snap-frozen tissues. Vet. Immunol. Immunopathol. 2007, 120, 212–222. [Google Scholar] [CrossRef] [PubMed]
  93. Hartung, J. A note on combining dependent tests of significance. Biom. J. 1999, 41, 849–855. [Google Scholar] [CrossRef]
Ijms 23 12269 g001 550
Figure 1. Normalized mRNA gene levels in the midbrain tissues of tame rats (white bars) versus aggressive rats (gray bars) for two assays (Ascl3 and Defb17) and two reference (Ppia and Rpl30) genes. Genes Ascl3 and Defb17 are statistically significantly overexpressed in the MT. The bar heights denote mean values, error bars represent standard error, and asterisks denote statistical significance (p < 0.05) according to both the Mann–Whitney U test and Fisher’s Z-test.
Figure 1. Normalized mRNA gene levels in the midbrain tissues of tame rats (white bars) versus aggressive rats (gray bars) for two assays (Ascl3 and Defb17) and two reference (Ppia and Rpl30) genes. Genes Ascl3 and Defb17 are statistically significantly overexpressed in the MT. The bar heights denote mean values, error bars represent standard error, and asterisks denote statistical significance (p < 0.05) according to both the Mann–Whitney U test and Fisher’s Z-test.
Ijms 23 12269 g001
Ijms 23 12269 g002 550
Figure 2. Ranking and classification of the overrepresented motifs for known TFs in the promoters of DEGs selected according to criteria PADJ  <  0.01 and log2 fold > 1. Panels (A,B) show 8/25 top-scoring motifs detected for the promoters of genes upregulated and downregulated in the brain tissue samples of tame rats, respectively (see Table 2). Motif overrepresentation was estimated by the ESDEG tool (see Materials and Methods, Section 4.4). Motifs of known TFs were compiled from the murine core collection of the Hocomoco database (see Materials and Methods). On the left, each panel shows four columns, indicating the ranks of motifs, the significance of their enrichment computed by ESDEG [(log10(padj)], the respective TF families according to the classification of murine TFs based on their DNA-binding domains [17], and the names of TFs (Hocomoco ID). The yellow color in panel B indicates motifs for TFs from the bHLH class belonging to three closely related families (E2A-related factors {1.2.1}, MyoD/ASC-related factors {1.2.2}, and Tal-related factors {1.2.3}) possessing a similar E-box binding motif. On the right, each panel presents the bottom left corner of the symmetrical matrix, indicating pairwise similarities of motifs as estimated by the permutation test described earlier [29]. Empty cells imply significantly distinct motifs (p-value > 0.05), the shades from light orange to dark orange denote significantly similar motifs (p-value < 0.05), and the numbers indicate the logarithm of pairwise significance: (log10(p-value).
Figure 2. Ranking and classification of the overrepresented motifs for known TFs in the promoters of DEGs selected according to criteria PADJ  <  0.01 and log2 fold > 1. Panels (A,B) show 8/25 top-scoring motifs detected for the promoters of genes upregulated and downregulated in the brain tissue samples of tame rats, respectively (see Table 2). Motif overrepresentation was estimated by the ESDEG tool (see Materials and Methods, Section 4.4). Motifs of known TFs were compiled from the murine core collection of the Hocomoco database (see Materials and Methods). On the left, each panel shows four columns, indicating the ranks of motifs, the significance of their enrichment computed by ESDEG [(log10(padj)], the respective TF families according to the classification of murine TFs based on their DNA-binding domains [17], and the names of TFs (Hocomoco ID). The yellow color in panel B indicates motifs for TFs from the bHLH class belonging to three closely related families (E2A-related factors {1.2.1}, MyoD/ASC-related factors {1.2.2}, and Tal-related factors {1.2.3}) possessing a similar E-box binding motif. On the right, each panel presents the bottom left corner of the symmetrical matrix, indicating pairwise similarities of motifs as estimated by the permutation test described earlier [29]. Empty cells imply significantly distinct motifs (p-value > 0.05), the shades from light orange to dark orange denote significantly similar motifs (p-value < 0.05), and the numbers indicate the logarithm of pairwise significance: (log10(p-value).
Ijms 23 12269 g002
Ijms 23 12269 g003 550
Figure 3. Proposed molecular mechanism for changes in gene expression between tame and aggressive rats due to differences in ASCL3 expression and its involvement in the regulation of these genes. TFs of the bHLH class acting as transcription activators are shown in red/pink/orange. Blue color shows TF ASCL3 acting as a repressor for target genes of various TFs of the bHLH class.
Figure 3. Proposed molecular mechanism for changes in gene expression between tame and aggressive rats due to differences in ASCL3 expression and its involvement in the regulation of these genes. TFs of the bHLH class acting as transcription activators are shown in red/pink/orange. Blue color shows TF ASCL3 acting as a repressor for target genes of various TFs of the bHLH class.
Ijms 23 12269 g003
Table
Table 1. Summary statistics for differentially expressed genes (DEGs) in the tegmental region of the midbrain (MT) transcriptomes of three tame adult male rats (Rattus norvegicus) and three aggressive rats.
Table 1. Summary statistics for differentially expressed genes (DEGs) in the tegmental region of the midbrain (MT) transcriptomes of three tame adult male rats (Rattus norvegicus) and three aggressive rats.
GroupCount
Total number of sequence reads (NCBI SRA ID: PRJNA668014)182,197,974
Reads mapped to reference rat genome RGSC Rnor_6.0, UCSC Rn6, July 2014 (%)158,310,590 (86.9%)
Expressed genes identified14,039
Statistically significant DEGs (|log2 fold| > 1, PADJ < 0.1, Fisher’s Z-test with Benjamini–Hochberg correction)42
During the processing of the RNA-seq data, significant differences were revealed between the aggressive and tame strains of rats in the transcriptome profile of the MT.
Table
Table 2. DEGs in the MT of tame vs. aggressive rats.
Table 2. DEGs in the MT of tame vs. aggressive rats.
#No.Rat Gene, NameSymbollog2 foldPADJ
1defensin beta 17Defb177.311.86 × 10−3
2leukocyte immunoglobulin-like receptor, subfamily B member 3-likeLilrb3l5.719.15 × 10−2
3lipoxygenase homology PLAT domains 1Loxhd14.027.91 × 10−2
4sucrase-isomaltaseSi3.308.66 × 10−2
5achaete-scute family bHLH transcription factor 3Ascl3 *2.823.16 × 10−4
6protocadherin gamma subfamily A, 1Pcdhga12.802.21 × 10−2
7 LOC1009108022.256.64 × 10−2
8Fc gamma receptor IIbFcgr2b2.012.46 × 10−2
9cysteine-rich secretory protein LCCL domain containing 1Crispld11.497.40 × 10−2
10collagen triple helix repeat containing 1Cthrc11.446.13 × 10−2
11solute carrier family 7 member 11Slc7a111.379.84 × 10−2
12nicotinamide nucleotide adenylyltransferase 1Nmnat11.311.76 × 10−11
135-hydroxytryptamine receptor 3AHtr3a1.304.33 × 10−2
14macrophage expressed 1Mpeg11.288.30 × 10−4
15schlafen family member 13Slfn131.266.30 × 10−2
16syncoilin, intermediate filament proteinSync1.195.31 × 10−2
17olfactomedin-like 1Olfml11.146.73 × 10−2
18phosphotriesterase relatedPter1.137.16 × 10−2
19collagen type XV alpha 1 chainCol15a11.125.45 × 10−2
20MORN repeat containing 1Morn11.101.43 × 10−4
21somatostatin receptor 2Sstr2−1.012.46 × 10−2
22retinol saturaseRetsat−1.031.33 × 10−2
23glycerol-3-phosphate dehydrogenase 1Gpd1−1.052.46 × 10−2
24aquaporin 9Aqp9−1.229.89 × 10−2
25serine peptidase inhibitor, Kunitz type 1Spint1−1.291.43 × 10−4
26activity-regulated cytoskeleton-associated proteinArc−1.356.13 × 10−2
27t-complex-associated testis expressed 1Tcte1−1.416.30 × 10−2
28tetratricopeptide repeat domain 22Ttc22−1.445.08 × 10−2
29keratin2Krt2−1.548.63 × 10−5
30Fos proto-oncogene, AP-1 transcription factor subunitFos*−1.636.87 × 10−2
31synaptophysin-like 2Sypl2−1.658.72 × 10−2
32Ly6/Plaur domain containing 3Lypd3−1.865.08 × 10−2
335-hydroxytryptamine (serotonin) receptor 5BHtr5b−1.949.17 × 10−2
34FosB proto-oncogene, AP-1 transcription factor subunitFosb*−1.952.46 × 10−2
35heat shock protein family A (Hsp70) member 1AHspa1a−2.104.25 × 10−2
36heat shock protein family A (Hsp70) member 1BHspa1b−2.193.42 × 10−8
37vasoactive intestinal peptideVip−2.358.93 × 10−3
38 RGD1565611−2.946.69 × 10−2
39tachykinin precursor 3Tac3−3.642.46 × 10−2
40protein tyrosine phosphatase, non-receptor type 20Ptpn20−4.499.17 × 10−2
41hypocretin neuropeptide precursorHcrt−5.156.56 × 10−2
42hemoglobin, beta adult major chainHbb-b1−7.781.09 × 10−5
Note. log2 fold: log2-transformed fold change (i.e., ratio of an expression level of a given gene in tame rats to that in aggressive rats); PADJ: statistical significance according to Fisher’s Z-test with Benjamini–Hochberg correction for multiple comparisons. Results are ordered by log2 fold change. Positive log2 fold change values indicate genes upregulated in tame rat samples; negative log2 fold change values indicate genes upregulated in aggressive rat samples. Transcription factor genes are marked with an asterisk.
Table
Table 3. qPCR primers selected using the publicly available PrimerBLAST web service [89].
Table 3. qPCR primers selected using the publicly available PrimerBLAST web service [89].
No.GeneForward, 5′→3′Reverse, 5′→3′
DEGs identified in the hippocampus of tame versus aggressive adult male rats [This Work]
1Ascl3CCTCTGCTGCCCTTTTCCAGACTTGACTCGCTGCCTCTCT
2Defb17TGGTAGCTTGGACTTGAGGAAAGAATGCAGCAGTGTGTTCCAGGTC
Reference Genes
3PpiaTTCCAGGATTCATGTGCCAGCTTGCCATCCAGCCACTC
4Rpl30CATCTTGGCGTCTGATCTTGTCAGAGTCTGTTTGTACCCC
Notes. Regarding the DEGs subjected to this qPCR verification, see Table 2; reference rat genes: Ppia (peptidylprolyl isomerase A [91]) and Rpl30 (ribosomal protein L30 [92]).
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Oshchepkov, D.; Chadaeva, I.; Kozhemyakina, R.; Shikhevich, S.; Sharypova, E.; Savinkova, L.; Klimova, N.V.; Tsukanov, A.; Levitsky, V.G.; Markel, A.L. Transcription Factors as Important Regulators of Changes in Behavior through Domestication of Gray Rats: Quantitative Data from RNA Sequencing. Int. J. Mol. Sci. 2022, 23, 12269. https://doi.org/10.3390/ijms232012269

AMA Style

Oshchepkov D, Chadaeva I, Kozhemyakina R, Shikhevich S, Sharypova E, Savinkova L, Klimova NV, Tsukanov A, Levitsky VG, Markel AL. Transcription Factors as Important Regulators of Changes in Behavior through Domestication of Gray Rats: Quantitative Data from RNA Sequencing. International Journal of Molecular Sciences. 2022; 23(20):12269. https://doi.org/10.3390/ijms232012269

Chicago/Turabian Style

Oshchepkov, Dmitry, Irina Chadaeva, Rimma Kozhemyakina, Svetlana Shikhevich, Ekaterina Sharypova, Ludmila Savinkova, Natalya V. Klimova, Anton Tsukanov, Victor G. Levitsky, and Arcady L. Markel. 2022. "Transcription Factors as Important Regulators of Changes in Behavior through Domestication of Gray Rats: Quantitative Data from RNA Sequencing" International Journal of Molecular Sciences 23, no. 20: 12269. https://doi.org/10.3390/ijms232012269

APA Style

Oshchepkov, D., Chadaeva, I., Kozhemyakina, R., Shikhevich, S., Sharypova, E., Savinkova, L., Klimova, N. V., Tsukanov, A., Levitsky, V. G., & Markel, A. L. (2022). Transcription Factors as Important Regulators of Changes in Behavior through Domestication of Gray Rats: Quantitative Data from RNA Sequencing. International Journal of Molecular Sciences, 23(20), 12269. https://doi.org/10.3390/ijms232012269

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