One of the domesticated species; the dog has been selectively bred for various aims by human. The dog has many breeds, which are artificially selected for specific behaviors and morphologies. Dogs contribute their life to human as working dogs for guide, rescue, detection or etc. Working dogs requires good personality, such as gentleness, robustness and patience for performing their special duty. Many studies have concentrated on finding genetic marker for selecting the high-quality working dog. In this study, we confirmed quantitative expression patterns of eight genes ( ABAT ; 4-Aminobutyrate Aminotransferase, PLCB1 ; Phospholipase C, Beta 1, SLC10A4 ; Solute Carrier Family 10, Member 4, WNT1 ; Wingless-Type MMTV Integration Site Family, Member 1, BARX2 ; BarH-Like Homeobox 2, NEUROD6 ; Neuronal Differentiation 6, SEPT9; Septin 9 and TBR1 ; T-Box, Brain, 1) among brains tissues from four dog breeds (Beagle, Sapsaree, Shepherd and Jindo), because these genes were expressed and have functions in brain mostly. Specially, BARX2 , SEPT9 , SLC10A4 , TBR1 and WNT1 genes were highly expressed in Beagle and Jindo, and Sapsaree and German Shepherd were vice versa. The biological significance of total genes was estimated by database for annotation, visualization and integrated discovery (DAVID) to determine a different gene ontology (GO) class. In these analyses, we suppose to these eight genes could provide influential information for brain development, and intelligence of organisms. Taken together, these results could provide clues to discover biomarker related to functional traits in brain, and beneficial for selecting superior working dogs. The dog ( Canis lupus familiaris ) had been domesticated from wild gray wolves ( Canis lupus ). Based on archeological data, the dogs had had first artificial selection that from 100,000 to 15,000 years ago in multiple locations, including Europe, the Middle East and East Asia [48 , 49] . Nowadays, the dog population is separated into more than 400 breeds exist in worldwide [36] . The dogs evolved through a mutually valuable relationship with human beings, and their abilities have been developed to perform an outstanding variety of working or special roles. These roles include military watch, security guards, shepherds, guides, rescue and pets [45] . Working dogs are required to good personality of gentleness, robustness and patience for performing their special duty. Several animal personality are influenced by the activity of specific genes [17 , 42] . Specially, the brain specific genes controls behavior, personality, or aggression, therefore it is needed to study for confirming the gene expression patterns in brain. Microarray studies have been performed to assess changes in behavior with gene expression patterns in the brain [46] . We tried to identify such expression changes of eight genes in brain of four domesticated dog breeds (Beagle, German Shepherd, Sapsaree and Jindo), then compared to their expression patterns. A major goal in the brain and behavioral sciences is to identify genes that influence social behavior and understand how their gene products influence the structure and function of the nervous system [17 , 42 , 46] . Thus, the aim of this study was to predict the relationships among personality traits (calmness, trainability, dog sociability, boldness, and etc.) of dogs. Gene expression profiles may reflect the complete personality of regulatory pathways; therefore we describe the gene expression patterns of eight genes ( ABAT, BARX2, NEUROD6, SEPT9, SLC10A4, PLCB1, TBR1, WNT1 ) in dog breeds. The domestic dogs display an extraordinary level of phenotypic diversity in personality and behavior, because dog breeding was introduced as various methods by human during the nineteenth century [48 , 50] . As the results of artificial selection, each dog breeds have fixed phenotypic traits [12 , 47] . In this point, it is important to select the suitable traits of special dog from many dog breeds. In this study, we presented each gene expression patterns in four dog breeds, which provides the clue to study of each dog breeds of specific traits. The Beagle is one breed of small to medium-sized dog, and developed primarily for tracking wild animals as detection dogs. Their great sense of smell and tracking ability is used to detect food items, drug, or explosive detection [6] . Due to their strength, intelligence, trainability and obedience, German Shepherds are generally used in search-and-rescue, police, or military dogs around the world. German Shepherds have abilities of tracking, patrolling, or detection by training [39] . Therefore, the German Shepherd is the widely used breeds for military, police or scent-work roles. As the Korean traditional breeds, Sapsaree is known as dauntlessness and loyalty aspects. Jindo is also one of the Korean traditional breeds, and well known for its unwavering loyalty and gentle nature. Because the Sapsaree and Jindo are an active and intelligent dog, it requires frequent interaction with people or another dog in the family [14 , 28] . Although they are frequently used as pet dog or guide dog, they are not frequently used as special dogs. In order to use these breeds as special dogs, gene expression patterns in brain could be good clues. One post-mortem brain tissue sample was extracted from four dog breeds (Beagle, German Shepherd, Sapsaree, and Jindo). Beagle and Sapsaree brain tissues were obtained by Chungnam National University with approval by the Animal Ethics Committee (CNU-00199). German Shepherd and Jindo brain tissues by Rural Development Administration (Jeonju, Korea), and the animals received care in accordance with the standard guidelines for the Care and Use of Laboratory Animals provided by the National Institute of Animal Science Animal Care Committee, and the experiment was executed with approval from the animal ethics committee under the operation rule of animal experiment ethics at the National Institute of Animal Science (approval number: 2014-085). Cellular RNA of dog brain tissues was isolated by TRIzol® Reagent (Invitrogen, Carlsbad, CA, USA) to purify total RNA according to the manufacturer’s guidelines. After RNA isolation, the quality and quantity of the resulting single-stranded RNAs were assessed using a ND-1000 spectrophotometer (NanoDrop, Wilmington, DE, USA). Total RNA was treated and reverse-transcribed using a Prime-Script RT reagent kit with genomic DNA Eraser (Takara Bio, Shiga, Japan) according to the instructions of the manufacturer. Eight pairs of primer were designed to detect the mRNA of each gene according to the open reading frames (ORFs) in the whole genome ( Table 1 ) Quantitative real-time RT-PCR was performed with the Rotor-Gene Q system (QIAGEN, Hilden, Germany) with a gene-specific primer set ( Table 1 ). Each of amplification reaction mixture (20 μl) contained 7 μl of H ₂ O, 10 μl of QuantiTech SYBR Green PCR Master Mix (QIAGEN, Hilden, NW, Germany), 1 μl each of forward and reverse primers at 10 nmol/ μl, and 1 μl of cDNA template. In addition, to confirm non-specific background amplification, we amplified a template control without cDNA. Real-time RT-PCR amplifications for target genes and housekeeping genes were conducted as follows: 30 cycles each of 95℃ for 15 s, annealing temperature for 15 s, and 72℃ for 15 s. Annealing temperature range is set from 54 to 57℃ depends on genes. Melting curve analysis was performed for 30 s at 65-99℃. To guarantee reproducibility, we amplified all samples in triplicate and the results were averaged. As a standard control, we used GAPDH (Glyceraldehyde-3-phosphate dehydrogenase) in gene expression, for normalization of real-time RT-PCR amplification. Functional annotation of eight genes was analyzed by the DAVID (Database for Annotation, Visualization and integrated Discovery) [16] . DAVID calculates P -values to demonstrate GO terms enrichment, where P-values less than 0.05 are considered to be strongly enriched in the annotation category after Benjamini multiple test correction. We then performed gene ontology (GO) term analyses and each genes were grouped as GO terms. The results of significant GO terms were queried to the REViGO program in order to construct a scatterplot and interactive graph [44] , and then we summarized GO terms on the 2D semantic space by semantic similarities. P -values were originated from the Benjamini and Hochberg false discovery rate (FDR), and the color of circles indicates enriched GO terms with FDR <0.05. Our gene expression data indicate that similar expression patterns were presented in four breeds of dog. Generally, Beagle and Jindo are dominant expressed patterns, whereas German shepherd and Sapsaree are lower expressed than Beagle and Jindo ( Fig. 1 ). ABAT has several psychiatry roles, and regulate specially gamma-aminobutyric acid (GABA) in neuronal cells by catalyzing the degradation of GABA. Thus low-expressed ABAT gene induces to enhanced amount of GABA in the synaptic junctions, then increases GABA-mediated signaling by means of the GABA receptors [26] . Specially, GABA has no means to penetrate the blood-brain barrier, so it must be synthesized in the brain. Therefore to elucidate ABAT expression patterns in the brain is important for neuronal processes and brain development [51] . In the brain of four dog breeds, ABAT is high-expressed in Beagle, and low-expressed in Shepherd. The two breeds of Sapsaree and Jindo, originated from Korea, have similar expression patterns. BARX2, SLC10A4, TBR1 and WNT1 genes have same expression patterns. BARX2 gene is known as transcription factor of cell adhesion [33] . The mechanisms of cell adhesion is very important for brain morphology, and functions such as learning, signal transduction, and memory [40] . When the early development of the nervous system, neurons maintain synapses by formation of cell-cell adhesions. Differently expressed patterns could provide the formation of brain cell, and the features of personality or intelligence. TBR1 has a role in glutamatergic projection in neuron differentiation. Glutamatergic neurons express receptors for the excitatory neurotransmitter glutamate as opposed to receptors for the inhibitory neurotransmitter GABA [23] . Solute carriers (SLCs) have a role for the transmembrane transport of various materials, such as amino acids, sugars, or inorganic ions. SLC10A4, one of the solute carrier family, has a role of carrying solution in blood or body fluid. SLC10A1 genes were involved in blood-brain barrier [19] , and SLC10A4 also has possible important roles in brain. Therefore, the problems of SLC in the brain could induce mental illness such as ADHD, depression or psychiatric disorders [38] . The study of this gene in the dog could help to elucidate the canine mental illness or behavior disorders. Proto-oncogene protein WNT1, has been originally considered as a candidate gene for Joubert syndrome, an autosomal recessive disorder with cerebellar hypoplasia as a leading feature. Joubert syndrome is a rare genetic disorder that affects the cerebellum, an area of the brain that controls balance and coordination [31 , 32] . Wnt signaling pathway controlled by WNT1, is stimulated by BDNF (Brain-derived neurotrophic factor), then induce the proliferation and differentiation of neural stem cells [13] . These genes same expression patterns have crucial roles in brain, and each gene pathway studies will be performed as further studies. Three genes, NEUROD6, SEPT9 and PLCB1 , were similar expressed patterns and highest expressed in Beagle. NEUROD6 is transcription factors, and associated to differentiation or development of the nervous system. It also regulates fasciculation and targeted axogenesis in the neocortex [8] . SEPT9 (septin 9) is involved in cytokinesis and cell cycle control, and known as a candidate for the ovarian tumor suppressor gene [ 20 ]. Mutations in this gene cause hereditary neuralgic amyotrophy [30] , and a chromosomal translocation involving this gene on chromosome 17 and the MLL gene on chromosome 11 results in acute myelomonocytic leukemia [29] . SEPT9 is highest expressed in Beagle, but we could not detect to specific expression patterns. SEPT9 is related to not only neural diseases but also diseases such as leukemia or tumor, therefore the specific expression patterns is need to be specified. PLCB1 catalyzes the formation of inositol 1,4,5-trisphosphate and diacylglycerol from phosphatidylinositol 4,5-bisphosphate by using calcium as a cofactor [21] . This reaction plays an important role in the intracellular transduction of many extracellular signals. The general features of total eight genes in the organisms were listed, and these genes have crucial roles for cellular metabolisms, DNA binding, or cell signaling ( Table 2 ). The variable gene expression patterns in four breeds of dog brain can provide clues for the studies of breed-specific traits. A GO analysis was performed in order to confirm the biological function of the eight genes. The gene set was confirmed as cellular GO terms, such as cell differentiation, multicellular organismal process, and cellular developmental process. The full list of statistically significant GO terms is enumerated ( Table 3 ). Total 12 functions were enriched, and the most significant function is cell differentiation ( P -value = 0.038) with four genes. Multicellular organismal process includes six of eight genes, except SLC10A4 and PLCB1 genes. BARX2, TBR1, and WNT1 genes were enriched in positive regulation of various biological processes, such as DNA-dependent transcription, RNA metabolic process, and nucleic acid metabolic process. Interestingly, these two genes were high-expressed in Beagle and Jindo, whereas low-expressed in Shepherd and Sapsaree. According to these facts, we summarized that positive regulation of biological process is enhanced in the brain of two breeds, Beagle and Jindo. We used the REViGO program to get the functional relationship of GO terms in the network structure. The scatter plots and integrated graph show the cell metabolism-related GO terms, including cell differentiation, multicellular organismal process, cellular developmental process, and positive regulation of transcription, DNA-template ( Fig. 2 ). WNT1 gene is included in all GO terms, and crucial factor of each network structure. WNT1 regulate the wnt signaling pathway, then induce neural stem cell differentiation [13] . Therefore, WNT1 gene can provide a clue for the study of wnt signaling brain cell, and their mental features. SEPT9 gene is also included in the scatter plots and integrated graph. SEPT9 gene is related to cell division, such as cell cytokinesis and cell cycle [20] . SEPT9 is also high-expressed in human lymphoid and malignant brain tumors [43] . We predict to these genes can provide crucial information for brain development, and intelligence of organisms. In conclusion, we selected eight genes related to brain functions, then confirmed their expression patterns, functions, and network. Total genes have a similar expression patterns, and their function is related to cell differentiation, cellular developmental process and multicellular organismal process. This data provides the fundamental clues for the studies of brain functions.