Targeted knockout of a chemokine-like gene increases anxiety and fear responses.

Emotional responses, such as fear and anxiety, are fundamentally important behavioral phenomena with strong fitness components in most animal species. Anxiety-related disorders continue to represent a major unmet medical need in our society, mostly because we still do not fully understand the mechanisms of these diseases. Animal models may speed up discovery of these mechanisms. The zebrafish is a highly promising model organism in this field. Here, we report the identification of a chemokine-like gene family, samdori (sam), and present functional characterization of one of its members, sam2 We show exclusive mRNA expression of sam2 in the CNS, predominantly in the dorsal habenula, telencephalon, and hypothalamus. We found knockout (KO) zebrafish to exhibit altered anxiety-related responses in the tank, scototaxis and shoaling assays, and increased crh mRNA expression in their hypothalamus compared with wild-type fish. To investigate generalizability of our findings to mammals, we developed a Sam2 KO mouse and compared it to wild-type littermates. Consistent with zebrafish findings, homozygous KO mice exhibited signs of elevated anxiety. We also found bath application of purified SAM2 protein to increase inhibitory postsynaptic transmission onto CRH neurons of the paraventricular nucleus. Finally, we identified a human homolog of SAM2, and were able to refine a candidate gene region encompassing SAM2, among 21 annotated genes, which is associated with intellectual disability and autism spectrum disorder in the 12q14.1 deletion syndrome. Taken together, these results suggest a crucial and evolutionarily conserved role of sam2 in regulating mechanisms associated with anxiety.

Emotional responses, such as fear (responses to the appearance of well-defined and clearly present aversive stimuli) and anxiety (responses seen under aversive contexts in which the clear presence of fear inducing stimuli cannot be ascertained or these stimuli are continuously present and diffuse), are essential behavioral phenomena that have strong fitness components in all species (1). Anxiety-related disorders continue to represent a major unmet medical need in our society. Despite concerted efforts to develop anxiolytic and antidepressant pharmacotherapies, a large proportion of patients remain unresponsive to medication (2). This is due to the fact that we still do not fully understand the mechanisms of these diseases. To facilitate mechanistic analysis and to speed up the process of target identification and psychopharmacological characterization of compounds, animal models have been proposed.

Emotional responses, such as fear and anxiety, are regulated by various neuromodulators (3). Recently, studies have raised the possibility that chemokines may also possess neuromodulatory functions. Chemokines are small, secreted signaling proteins that commonly act as chemo-attractants to guide the migration of cells of the immune system. Importantly, chemokines are also expressed in neurons and glial cells of the central nervous system (CNS); while their role in the brain has mainly been associated with neuro-inflammatory responses, it is becoming more evident that chemokines in the brain may hold other functions (4–6).

Emotion dysregulation can lead to serious behavioral problems and impaired social interaction associated with a variety of disease conditions, including attention deficit hyperactivity disorder, bipolar disorder, and posttraumatic stress disorder (7, 8). Despite the paucity of proper treatment options for such disorders, the exploration of possible underlying genetic and biological mechanisms has been slow. The identification of novel genes associated with emotional behavior and understanding neural mechanisms of higher brain function remains a major goal of biology.

The use of animal models has greatly aided this research. The zebrafish (Danio rerio) has emerged as a promising model organism for the study of complex neuropsychiatric diseases because of the well-defined behavioral phenotypes of this species (9), and because of its evolutionarily conserved brain structures and functions. One example of such a brain structure is the habenula (Hb), which is connected with the subpallium and hypothalamus, and is implicated in fear, anxiety, addiction, and mood disorders (10, 11).

The zebrafish is useful for large-scale forward genetic screens to discover novel genes involved in CNS development, an approach we have utilized (12, 13). Over 10 y ago we discovered a chemokine-like gene family, called samdori (sam), using genetic screens (i.e., insertional mutagenesis with a unique dual reporter system) (). We present herein the body of work that has extended from initial gene discovery and encompasses various analyses of the physiological functions of this family.

A member of this family, sam2 exhibits exclusive brain-specific expression, predominantly in neurons of the dorsal Hb (dHb), as well as the telencephalic area and hypothalamus. Here, we report the development of a zebrafish null mutant (KO) in which this gene is silenced. We investigate the effect of this null mutation using behavioral assays developed to induce and measure fear and anxiety responses. To investigate the generalizability of our findings obtained with zebrafish, we also report herein the development and analysis of a mouse null mutant in which the mammalian homolog of the zebrafish sam2 gene is silenced. We found phenotypical alterations in the Sam2 KO mice that are consistent with those we observed in the zebrafish null mutants, suggesting that Sam2 is critical for regulating anxiety across such evolutionarily distant species as mammals and fish.

In addition, we also analyzed a variety of markers in the zebrafish sam2 KO (colocalization studies), which uncovered changes in mRNA expression of the gene encoding corticotropin-releasing hormone (CRH), a hormone involved in stress and anxiety responses in mammals (14). We followed up this investigation with whole-cell patch-clamp electrophysiological recordings from CRH neurons in the mouse hypothalamus, and found that Sam2 modulated GABAergic synaptic transmission on CRH neurons.

Finally, we investigated whether mutations in the human genome potentially affecting SAM2 function may have functional, behavioral consequences. We identified patients with microdeletions/microduplications in 12q14.1, including SAM2, who suffered from a variety of behavioral and intellectual problems. Using six copy-number variations (CNVs), we were able to refine the candidate gene region encompassing SAM2 among 21 annotated genes at 12q14.1.

In summary, we describe the discovery of a chemokine-like protein family and provide initial functional characterization of one of its members, sam2. Taken together, our multidisciplinary results demonstrate a potentially crucial role of this candidate gene in the regulation of emotion-related behaviors across evolutionarily distant species, and delineate a target that may allow the development of effective pharmacotherapy for human CNS disorders associated with altered fear, anxiety, and stress.

Results

Identification of samdori2.

Using insertional mutagenesis to identify novel functional genes, we isolated a chemokine-like gene, sam [also called TAFA (15) and FAM19A in the National Center for Biotechnology Information database] (). Subsequently, we also isolated other members of sam family based upon the sequences in the database. We have identified eight members in the zebrafish sam gene family, and found them highly conserved in mouse and human at the amino acid sequence level (). All Sam proteins exhibited 10 regularly spaced cysteine residues with a pattern of CX7CCX13CXCX14CX11CX4CX5CX10C, a sequence similar to that of CC-type chemokines (15).

Using whole-mount in situ hybridization in zebrafish, we found all sam genes to be exclusively expressed in the CNS at larval stages (). Among them, we focused on the functional analysis of the sam2 gene, because we found sam2 to be predominantly expressed in the dHb of adult fish (Fig. 1). To further characterize sam2-expressing cells, we performed whole-mount in situ hybridization in larval and adult fish brains using various neuronal markers. In addition to dHb, sam2-expressing cells were also detected in the telencephalon and otic neurons in larvae (). In situ hybridization of adult brain and serial sectioning revealed the detailed expressions of sam2 in the medial zone of the dorsal telencephalic area (Dm), the dorsal nucleus of the ventral telencephalon (Vd), preoptic area of hypothalamus (PPa), and other hypothalamic regions (ATN, anterior tuberal nucleus; Hc, caudal zone of periventricular hypothalamus; Hd, dorsal zone of periventricular hypothalamus; and PTN, posterior tuberal nucleus) (Fig. 1 and ). To assess the distribution of sam2 within excitatory and inhibitory circuits, we performed double in situ hybridization with GABAergic markers (16, 17) and glutamatergic markers (16, 17). We found sam2+ cells in GAD67/gad1b-expressing cells in the Vd, central nucleus of the ventral region (Vc), and PPa (Fig. 1 and ). In the diencephalon, sam2-expressing cells were observed adjacent to GAD65/gad2+ cells in the hypothalamic area, such as the ATN and PTN, and periventricular hypothalamus (). With the glutamatergic markers, sam2 mRNA was partially overlapped with vglut2b in the Dm but not in other forebrain regions (Fig. 1 ). In addition, sam2 was found to overlap with vglut2a in the dorsal region of the Hb (Fig. 1). Expression of sam2 in the Hb was further examined by costaining with dHb-specific kctd12.1 (lov) (18, 19) or ventral Hb (vHb)-specific marker, aoc1 or dHb/vHb marker fam84b (20) (). sam2 expression was found not to colocalize with known markers of dopaminergic (21), oxytocinergic (22), or hypocretinergic (23) neurons (Fig. 1 and ). The Hb has been implicated in the regulation of emotion, including anxiety, fear, depression, and reward (7, 24). The fish telencephalic regions are homologous to the mammalian amygdala (Dm) and striatum (Vd), structures that are known to play fundamental roles in emotion-related behaviors in mammals (25). Thus, expression of the sam2 gene detected in the dHb, telencephalic areas (Dm, Vd), and hypothalamic regions prompted us to investigate its function in the regulation of emotion.

Fig. 1.

Characterization of sam2-expressing cells in the adult zebrafish brain. (A and A′) Whole-mount two color in situ hybridization of a dissected zebrafish brain with sam2 and a dopaminergic neuron marker, th. Anterior is to the left; lateral (A) and dorsal view (A′). The sam2 expression region did not overlap with th+ neurons. Prominent sam2 expression is seen in the Vd and Hb as well as hypothalamic regions. (B–D) Section images of brain hybridized with sam2/vglut2b (B, sagittal section; C, cross-section) or sam2/gad1b (D). (E and F) Cross-section of Hb region stained with sam2 alone (E) or sam2/vglut2a (F) probes. AN, auditory nerve; D, area dorsalis telencephlali; Dc, central zone of area dorsalis telencephali; MO, medulla oblongata; OB, olfactory bulb; TeO, optic tectum; Vv, ventral nucleus of area ventral telencephali. (Scale bars, 100 µm.)

Normal Development of Neuronal Cells in sam2 KO Fish.

To investigate the function of sam2, we created two sam2 KO zebrafish lines, taking advantage of targeted mutagenesis utilizing zinc-finger nucleases (ZFNs) (26, 27). The 5-bp deletion (sam2) and 17-bp deletion (sam2) of the two alleles induced frame-shift mutations (Fig. 2 ), which lead to alteration of the protein-coding sequence (). Both sam2 KO lines showed normal morphology during embryonic development (Fig. 2 ), survived to adulthood, and were fertile. To test whether sam2 is necessary for the development of specific neuronal cell-types, we examined the expression of various molecular markers in sam2 KO zebrafish (). Compared with wild-type fish, we could not detect any significant differences in the expression patterns of dopaminergic neuronal markers [tyrosine hydroxylase (th) (21), dopamine transporter (dat) (28), and nuclear receptor related 1 protein (nurr-1) (28)], serotonergic neuronal markers [tryptophan hydroxylase raphe (tphR) and 5-hydroxytryptophan (5-HT)] (29), oxytocinergic neuronal marker oxytocin (oxt) (22), and interpeduncular nucleus (IPN) marker somatostatin1.1 (sst1.1) (30). Moreover, sam2 KO embryos showed normal expression of neurohormonal genes (), including hypocretin/orexin (hcrt) (23), melanin concentrating hormone (mch) (31), agouti related protein homolog (agrp) (32), neuropeptide Y (npy) (33), and pet-1 (fev; ETS oncogene family) (29). In addition, we observed normal hypothalamus–locus coeruleus projections in progenies of the Tg(hcrt:mEGFP) (23) transgenic line crossed to the sam2 mutant (). Moreover, the targeting of Hb efferent axons innervating the IPN appeared to be largely unaffected in the Et(-1.0otpa:mmGFP)hd1 (34) transgenic embryos crossed to sam2 KO fish (Fig. 2 ). Therefore, we conclude that sam2 is not directly involved in the generation, migration, and axonal growth of neurons.

Fig. 2.

Generation of sam2 KO fish using ZFNs. (A) DNA sequencing analysis of two sam2 KO alleles (sam2 and sam2) with 5-bp and 17-bp deletion, respectively. Yellow mark, Left and Right ZFN-binding regions; red, spacer. (B and B′) Allele-specific genotyping of sam2 and sam2 KO fish by genomic PCR and digestion for BtgI site (CCGTGG) newly created in sam2. (C and C′) Normal morphology of sam2 KO embryo at 3 d. (D–E′) Projections of Hb efferent axons targeting the IPN in the Et(-1.0otpa:mmGFP)hd1 transgenic sam2 KO fish were not visibly affected. Lateral views (D and E) and dorsal views (D′ and E′). dIPN, dorsal IPN; FR, fasciculus retroflexus; lFR, Left FR; lHb, Left Hb; MR, median raphe; rFR, Right FR; rHb, Right HbvIPN, ventral IPN. (Scale bar, 100 µm.)

Increase of Anxiety-Like Responses in sam2 Mutant Fish.

We next examined whether the loss of sam2 affects behaviors associated with anxiety using the open (novel) tank (35, 36) and dark/light preference (scototaxis) (37) tests. Anxiety and fear are related phenomena, but can be distinguished. Responses induced by a present and clearly detectable aversive stimulus are regarded as fear responses (behavior aimed at the avoidance of the stimulus), whereas anxiety-like behavior is found under diffuse aversive conditions: that is, without the appearance or clear presence of such stimuli. Novelty is an aversive context without the presence of specific fear-inducing stimuli, and thus novel test situations are often used to quantify anxiety. In the novel tank, when placed individually, sam2 KO fish showed normal locomotor activity, as indicated by no changes in total distance traveled (Fig. 3) (sam2, 4,359.29 ± 218.68 cm, n = 18; sam2, 4,088.71 ± 249.83 cm, n = 28) and in average velocity (Fig. 3) (sam2, 7.63 ± 0.37 cm/s, n = 18; sam2, 7.02 ± 0.41 cm/s, n = 28), compared with control sibling fish. Moreover, there were no significant differences in the frequency of transition to the upper half of the test tank between sam2 KO and control fish (Fig. 3) (sam2, 30 ± 6.93, n = 18; sam2, 34.86 ± 4.19, n = 28), a measure of vertical exploration. However, sam2 KO fish exhibited behavioral characteristics indicative of elevated levels of anxiety, including an increased number of erratic movements (1, 38) (Fig. 3) (sam2, 52.67 ± 9.19, n = 18; sam2, 142.14 ± 12.50, n = 28) and number of freezing events (35, 36) (Fig. 3) (sam2, 0.17 ± 0.09, n = 18; sam2, 0.44 ± 0.10, n = 27). Freezing is a typical response seen in aversive contexts and can be used to quantify the level of anxiety in a variety of species, including zebrafish and rodents (36). Erratic movement has also been shown to be a characteristic sign of anxiety or fear in zebrafish (36, 39). Thus, our results indicate that sam2 may be involved in the regulation of anxiety-like responses in zebrafish. To further investigate anxiety-like behavior in sam2 KO fish, we also measured thigmotaxis (avoidance of open areas) in the novel tank. Increased thigmotaxis in novel contexts has been interpreted as a sign of anxiety (36). sam2 KO fish showed higher preference for the corner over the center of the test tank () (sam2 at the corner, 21.22 ± 0.56 s, at the center 34.15 ± 0.70 s, n = 12; sam2 at the corner, 26.83 ± 0.95 s, at the center, 29.14 ± 0.89 s, n = 12). The light/dark paradigm is another frequently employed test of fear and anxiety in zebrafish and rodents (37), and zebrafish have been shown to exhibit scototaxis (i.e., dark preference). Consistent with the results obtained in the novel tank task, in the light/dark paradigm, adult sam2 KO fish traveled a shorter distance (Fig. 3) (sam2, 1,533.59 ± 191.18 cm, n = 18; sam2, 788.88 ± 116.09 cm, n = 24), spent less time in the white arena () (sam2 in the black, 733.125 ± 80.66 s, in the white 168.38 ± 80.42 s, n = 18; sam2 in the black, 811.5 ± 64.81 s, in the white 104.08 ± 87.67 s, n = 24), and exhibited reduced frequency of black/white transitions compared with wild-type (Fig. 3) (sam2, 81.67 ± 8.74, n = 18; sam2, 38.69 ± 4.36, n = 24), alterations that further corroborate the notion that the null mutant fish suffered from elevated levels of anxiety. These results suggest that sam2 may regulate emotion-related behaviors.

Fig. 3.

Increase of anxiety-related behaviors in sam2 KO zebrafish. (A–E) Novel tank tests of individually placed fish. The distance traveled (A) (sam2, n = 18; sam2, n = 28; Cohen’s d = 0.24; Mann–Whitney U = 181, P = 0.11), average velocity (B) (sam2, n = 18; sam2, n = 28; Cohen’s d = 0.32; Mann–Whitney U = 171, P = 0.07), and transition to upper half (C) (sam2, n = 18; sam2, n = 28; Cohen’s d = 0.19; Mann–Whitney U = 202, P = 0.27) were not significantly changed in sam2 KO fish. However, the frequency of erratic behavior (D) (sam2, n = 18; sam2, n = 28; Cohen’s d = 1.65; Mann–Whitney U = 41, P < 0.00001) and the number of freezing fish (E) (sam2, n = 18; sam2, n = 27; Cohen’s d = 0.62; P = 0.042, Student’s t test) were increased in sam2 KO fish. (F and G) Black/white preference (scototaxis) test. Distance traveled in white arena (F) (sam2, n = 18; sam2, n = 24; Cohen’s d = 1.06; Mann–Whitney U = 85, P = 0.0032) and frequency of transitions (G) (sam2, n = 18; sam2, n = 24; Cohen’s d = 1.42; Mann–Whitney U = 61.5, P = 0.00016) were measured. *P < 0.05; **P < 0.01; ***P < 0.001; ns, not significant (P > 0.05).

Increase of Social Cohesion in sam2 KO Fish.

Zebrafish form groups, a behavior called shoaling. One of the adaptive functions of this social behavior has been shown to be predator-avoidance. Measures of shoal cohesion have been successfully utilized in the analysis and modeling of anxiety and fear in vertebrates, including zebrafish (1, 36). To test whether sam2 KO fish exhibit altered shoaling behavior, a group of five sam2 KO fish were placed in a novel tank at a time and monitored by 3D video tracking. We recorded the swim paths of the fish in the group for 30 min (Fig. 4 and ) (40). Because wild-type zebrafish typically take 10 min to habituate to a novel tank (41), we analyzed the behavior of fish at two time points: early phase (before habituation, 5–8 min) and late phase (after habituation, 13–16 min) of the session. During the early phase, both sam2 KO and control fish showed typical anxiety-like responses by staying close to the bottom of the tank. However, at the late phase control sibling fish started to swim in the middle and the top layers of the water (Fig. 4 ), a response interpreted as reduction of anxiety. In sharp contrast, sam2 KO fish remained near the bottom of the tank (Fig. 4 and Movie S1), suggesting impaired habituation to the novel environment (i.e., maintenance of elevated levels of anxiety) in sam2 KO fish. We next examined social cohesion in sam2 KO zebrafish by measuring the distance between a focal fish and each of its shoal members (interindividual distance) (41–43). Control fish significantly increased their mean interindividual distance from the early phase to the late phase, suggesting the reduction of novelty-induced anxiety (Fig. 4 ). However, sam2 KO fish did not show such changes in mean individual distance, and maintained the same strong shoal cohesion (Fig. 4 and ). The effect size (Cohen’s d) for cohesion of the control and sam2 KO fish was d = −0.063 at the early phase (5–8 min) and d = −0.80 at the late phase (13–16 min). These results confirm that sam2 KO mutants exhibit elevated levels of anxiety compared with wild-type fish.

Fig. 4.

Anxiety-like behavior and increased social cohesion in sam2 KO fish. Five fish (3-mo-old male siblings) were placed as a group in a novel tank. (A and A′) A snapshot of video tracking after 20-min novel tank test. (B and B′) Temporal 3D reconstructions of video tracking before habituation (early, 5–8 min, blue; Cohen’s d = −0.063) and after habituation (late, 13–16 min, red; Cohen’s d = −0.80) to the novel environment. (C and C′) Measurement of social cohesion as the mean individual gap. We measured the distance between the focal fish and one of its shoal members (lines a, b, c, and d; interindividual distances). P = 0.02 for control and P = 0.88 for sam2 KO fish (sam2, n = 5; sam2, n = 5). *P < 0.05; ns, not significant (P > 0.05).

Excessive Anxiety/Fear Behavior in sam2 KO Fish in Response to Alarm Substance.

The above behavioral paradigms are considered to induce anxiety-like responses, because they do not use specific clearly present aversive stimuli. To test fear responses, we conducted another test in which an aversive stimulus was specifically provided and was clearly present during the test. A clearly present aversive stimulus that may be delivered in a well-defined manner to the zebrafish is the alarm substance (38). Among other behavioral responses, the alarm substance has been shown to induce a reliable alteration in shoaling that may be used to measure fear in zebrafish (38). In our study, we examined changes in shoal cohesion in response to administration of the alarm substance (). Both control and sam2 KO fish exhibited increased shoal cohesion (reduced distance among shoal members) in response to the delivery of alarm substance. However, the relative decrease of interindividual distance in sam2 KO fish was more robust compared with that of wild-type control siblings. The effect size between before and after alarm treatment was d = −1.09 for control and d = −1.22 for sam2 KO fish, respectively. The effect size between control and sam2 KO fish was d = −0.96 before and d = −1.57 after alarm treatment. This was despite that the former had already displayed small interindividual distance even before the delivery of the alarm substance (). Thus, taken together these results demonstrate elevated anxiety and fear in sam2 KO fish.

Increase of mRNA Expression of Stress-Related crhb in sam2 KO Fish.

Aversive stimuli that induce fear and anxiety are expected to lead to physiological changes: for example, neurohormonal stress responses mediated by the hypothalamic–pituitary–adrenal (HPA) axis (14, 44). A previous study demonstrated that disruption of glucocorticoid-negative feedback of the HPA axis resulted in elevated cortisol levels and anxiety-like behavior in zebrafish (14). We hypothesized that the increased level of anxiety in sam2 KO fish may also be associated with dysregulation of the HPA axis. Using in situ hybridization, we found sam2 KO fish to exhibit elevated crhb expression, an important HPA axis marker, in the PPa compared with wild-type (see, for example, Fig. 6 ). Increased crhb mRNA levels in the PPa, a region homologous to the mammalian paraventricular nucleus (PVN), is consistent with previous reports of dysregulated HPA axis function in zebrafish and mice (14, 45).

Fig. 6.

Increase of mRNA expression of stress-related crhb in sam2 KO fish and spontaneous inhibitory postsynaptic currents onto CRH neurons by SAM2. (A–D) Increase of crhb mRNA expression in sam2 KO fish (sam2+/+, n = 6; sam2−/−, n = 6). Ventral views of the whole brain of control sam2+/+ (A and C) and sam2−/− KO fish (B and D). (A′ and B′) Higher magnifications of the boxed regions in A and B are the PPa in zebrafish, homologous to the mouse PVN. (Scale bars, 200 µm.) (E) Representative photomicrograph of the Cre-dependent TdTomato in the PVN CRH neurons. (Magnification: E, 10×.) (F) Representative voltage-clamp traces of sIPSCs in response to SAM2 application in PVN CRH neurons. (G) SAM2 application significantly increased sIPSC frequency. (H) SAM2 did not affect sIPSC amplitude. (I) Representative voltage-clamp trace of eIPSC onto CRH neurons in response to SAM2 application. (J) SAM2 did not affect the amplitude of eIPSCs. (K) SAM2 did not change the paired-pulse ratio. *P < 0.05; **P < 0.01.

Sam2 KO Mice also Exhibit Anxiety-Like and Fear-Related Behaviors.

To investigate whether the role of sam2 in anxiety and fear is evolutionarily conserved or specific only to the zebrafish, we generated a Sam2 KO mouse line () and compared the behavior of these mutants to wild-type littermates. Sam2 showed brain-specific expression in mouse similar to that of zebrafish, having a higher expression level in the hippocampus and Hb but not in the PVN (). We found Sam2 KO mice to spend significantly less time in the open arms on an elevated plus maze (Fig. 5) (Sam2 mice, 38.3 ± 9.16%, n = 11; Sam2 mice, 19.66 ± 6.66%, n = 8) and to show increased freezing in both the contextual (Fig. 5) (Sam2 mice, 42.86 ± 8.08%, n = 11; Sam2 mice, 69.74 ± 6.12%, n = 7) and cued fear-conditioning tests (Fig. 5) (Sam2 mice, 40.33 ± 6.32%, n = 11; Sam2 mice, 70.95 ± 4.97%, n = 7) compared with wild-type littermates. The elevated freezing in the null mutant mice was not due to a motor function-related change, as both mutant and wild-type animals exhibited similar levels of freezing on training day (Fig. 5) [two-way repeated-measures ANOVA, Time × Group interaction: F(3, 48) = 0.48, P = 0.69, Sam2 mice, n = 11; Sam2 mice, n = 7]. Body weight (Fig. 5) (Sam2 mice, 24.86 ± 0.59 g, n = 25; Sam2 mice, 24.07 ± 0.39 g, n = 10) and total locomotion on an open-field test (Fig. 5) [two-way repeated-measures ANOVA, Time × Group interaction: F(5, 75) = 0.32, P = 0.89, Sam2 mice, n =11; Sam2 mice, n = 6] were not significantly different between Sam2 KO and wild-type control mice. Collectively, the Sam2 KO mice exhibited phenotypical alterations reminiscent of the KO zebrafish, suggesting that Sam2 plays a potentially critical and evolutionally conserved role in regulating anxiety- and fear-related behavioral responses.

Fig. 5.

Increase of anxiety-related behaviors in Sam2 KO mice. (A) Sam2 KO mice have normal body weight (Sam2, n = 25; Sam2, n = 10; Cohen’s d = 0.31; unpaired two-tailed, Mann–Whitney U = 124, P = 0.98). (B) Sam2 KO mice show normal total locomotion in an open-field test [Sam2, n = 11; Sam2, n = 6; two-way repeated-measures ANOVA, Time × Group interaction: F(5, 75) = 0.32, P = 0.89]. (C) Sam2 KO mice spent significantly less time in the open arms on an elevated plus maze (Sam2, n = 11; Sam2, n = 8; Cohen’s d = 0.74; unpaired one-tailed, Mann–Whitney U = 23, P = 0.045). (D) Sam2 KO mice show similar levels of freezing response during training days [Sam2, n = 11; Sam2, n = 7; two-way repeated-measures ANOVA, Time × Group interaction: F(3, 48) = 0.48, P = 0.69]. (E) Sam2 KO mice show higher freezing behavior in contextual test 24 h after training (Sam2, n = 11; Sam2, n = 7; Cohen’s d = 1.21; unpaired two-tailed, Mann–Whitney U = 15, P = 0.035). (F) Sam2 KO mice show higher freezing behavior in cued test 24 h after training (Sam2, n = 11; Sam2, n = 7; Cohen’s d = 1.75; unpaired two-tailed, Mann–Whitney U = 7, P = 0.0028). *P < 0.05; **P < 0.01; ns, not significant (P > 0.05).

Increase in the Frequency of Spontaneous Inhibitory Postsynaptic Currents by SAM2.

Increased activity in the stress axis is one of the most enduring findings in both animal models and human patients with anxiety. Elevated fear and anxiety responses in mice are accompanied by both Crh mRNA overexpression and increased excitability of the PVN neurons (14, 45). Similarly, selective knockout of Crh in the PVN results in reduced anxiety (46). Therefore, we hypothesized that Sam2 could regulate CRH neuron excitability in the PVN. To investigate this hypothesis, we exposed PVN CRH neurons to a 5-min-long bath application of SAM2 (1 µM) (). In response to this treatment, we found a significant increase in the frequency of spontaneous inhibitory postsynaptic currents (sIPSC); that is, 162.5 ± 39.76% of baseline (n = 11, P < 0.01, Wilcoxon signed-rank matched-pairs test) (Fig. 6 ). We found no significant change in the sIPSC amplitude, however (100.9 ± 5.3% of baseline, n = 11, P = 0.15, Wilcoxon signed-rank matched-pairs test) (Fig. 6). We also observed no significant differences in evoked IPSC (eIPSC) amplitude and paired-pulse ratio (Fig. 6 ). Thus, we conclude that SAM2 is able to significantly increase the frequency, but not the amplitude, of spontaneous GABAergic currents onto CRH neurons, rather than the amplitude of evoked GABAergic currents. These observations indicate that SAM2 may control tonic GABAergic inputs onto CRH neurons.

Comparative Genomic Mapping Using Six CNVs at 12q14.1.

Mutations (i.e., microdeletions or duplications) arising at the 12q14.1 region are characterized by a previously unreported feature set including craniofacial anomalies, cryptorchidism, intellectual disability, autism, and behavioral problems. Using six informative CNVs, we were able to refine the candidate gene region encompassing SAM2 at 12q14.1 (Fig. 7). We found deletions in two females, GC42855 and GC48823, sized 6.2 Mb (chr12: 62,097,144 to 68,266,543/hg19) and 5.1 Mb (chr12: 58,011,515 to 63,115,073/hg19), respectively (Fig. 7). GC48823 is 4-y-old with autism and seizures, whereas GC42855 is 6-y-old with developmental delay. A 1-Mb minimal overlapping region of these two deletions at 12q14.1 contains four genes including SAM2. Furthermore, we identified two more small CNVs in the DECIPHER (database of chromosomal imbalance and phenotype in humans using Ensembl resources) human genome variants database (47) overlapping the gene. Among them, patient 288660 with an 887-kb duplication (chr12: 61,271,591 to 62,158,216/hg19) has attention deficit hyperactivity disorder, autism spectrum disorder, and generalized joint laxity. However, this duplication overlaps only the 3′ end of the gene and may not disrupt its function. More importantly, patient 290951 has inherited an 880-kb deletion (chr12: 62,038,772 to 62,918,916/hg19) encompassing SAM2 and USP15 in a three-generation pedigree from the maternal grandfather through his mother with significant intellectual disability. Patient 290951 had a clinical diagnosis of autism spectrum disorder, behavioral difficulties, and severe speech and language delay. Overall he performed in the mild range of intellectual disability, as did his mother and maternal grandfather, who all had the same deletion. His unaffected brother, sister, and uncle did not inherit the deletion (Fig. 7 ). In mice, Usp15 is ubiquitously expressed in various tissues. Usp15 KO mice had normal survival rate and did not show abnormalities in development, although its deficiency enhanced T cell response against bacterial infection and tumor (48). A 61-kb minimal overlapping region of these four CNVs flanked by two blue vertical dotted lines encompasses the 56-kb 3′ end of SAM2. These data suggest that SAM2 contributes to the neurological phenotype of these patients. DECIPHER also has two additional small CNVs in this region, not including SAM2: child 251128 with behavioral problems, low-set ears, intellectual disability, and speech delay, has inherited an 886-kb deletion (chr12: 62,799,199 to 63,684,842/hg19) from a parent with similar phenotypes, and patient 287965 with a 224-kb duplication (chr12: 61,544,595 to 61,768,510/hg19) has global developmental delay.

Fig. 7.

CNV mapping at 12q14.1 with six human patients, implicating SAM2 as a candidate gene. Deletions and duplications are depicted in red and blue, respectively. (A) Alignment of three microdeletions and one microduplication at 12q14.1 has refined the candidate gene region to a 61 kb (arrows), which contains 56-kb 3′ end of SAM2. Notably, the microdeletion in a 6-y-old boy 290951 with autism spectrum disorder and intellectual disability contains SAM2, as well as USP15 and part of MON2. (B) Pedigree showing three family members affected by intellectual disability in three generations. We have confirmed the deletion in each family member; II-2 (del/+), III-2 (del/+), III-3 (+/+), IV-1 (+/+), IV-2 (+/+), and IV-3 (del/+). (C) SAM2 copy number in the family as determined by qPCR. The proband IV-3 has inherited the deletion from the maternal grandfather II-2 through his mother III-2 and affected members have only one SAM2 copy. The patient’s siblings are normal having two SAM2 copies. (D) Transcript levels of SAM2 in patient, his mother, sister, and white male control as revealed by qRT-qPCR. SAM2 transcripts are reduced approximately by half in both the patient and his mother relative to control.

Discussion

The functions of chemokines are well documented in the immune system. However, the potential role of CNS-expressed chemokines and chemokine-like proteins as neuromodulators is still largely unknown. In this study, we identified a chemokine-like sam gene family, showing CNS-specific expression. We provide empirical evidence for the involvement of a chemokine-like protein, SAM2, in vertebrate brain function and behavior associated with the regulation of anxiety and fear responses.

The targeted knockout of a single gene, sam2, in zebrafish produced an anxiogenic phenotype reminiscent of changes induced by Hb ablation. For example, ablation of the dHb in zebrafish has been shown to increase freezing and other behavioral signs of anxiety (25, 49–51). The anxiety-like and fear-related behavioral changes we detected in our zebrafish null mutants were also observed in our mouse KO model, highlighting the evolutionarily conserved role of Sam2 in the regulation of anxiety and fear in vertebrates.

Our in situ hybridization results showed the presence of sam2 transcripts in the telencephalic areas (Dm, Vd) and hypothalamic regions of the zebrafish brain. The Dm and Vd are known as the homologous regions for the mammalian amygdala and striatum, respectively. The involvement of these regions in anxiety has been demonstrated by numerous behavioral and molecular studies using both mammalian species and zebrafish (25, 52). In addition to traditional measures of anxiety and fear—such as freezing, thigmotaxis, or erratic movement—we also found evidence of increased anxiety in our zebrafish null mutants by measuring their shoaling responses. Shoaling is an adaptive and natural antipredatory behavior utilized in the analysis and modeling of vertebrate anxiety and fear (1). Furthermore, sam2 KO fish, when exposed to the alarm substance, exhibited significantly increased shoaling (reduction of interindividual distance) compared with wild-type. Also notably, novelty-induced increase of shoaling remained persistently elevated in the null mutants, whereas it showed the expected normal habituation (time-dependent reduction) in wild-type control fish (Fig. 4). Consistent with these findings, our Sam2 KO mice also exhibited increased anxiety-like responses, for example, on the elevated plus maze (Fig. 5).

In addition to the elevated anxiety-like and fear responses in our null mutant zebrafish, we found an increase in crhb mRNA expression in the hypothalamus. CRH neurons are under constant GABAergic suppression in vivo (53), but importantly, our experimental manipulation (i.e., bath application of SAM2) increased sIPSC frequency. These data indicate that SAM2 plays a role in tonic GABAergic suppression of CRH excitability, which may, at least in part, be responsible for the elevated anxiety-like responses observed in both KO models. It is also worth considering that the downstream effects of CRH and the resulting release of glucocorticoids may also contribute to the anxiety phenotype (14, 54).

Identification of the specific neural circuits mediating the effects of SAM2 first requires the determination of neuronal cell types or neurotransmitter systems involved. In case of glucocorticoid receptors, their actions on fear and anxiety behaviors have recently been shown to be mediated via glutamatergic but not GABAergic neural circuits in the mouse forebrain (45). We performed in situ hybridization analysis with GABAergic, glutamatergic, and sam2 markers in the adult fish brain to determine colocalization. sam2 expression within the Vd, Vc, Hc, and Hd was found to be predominantly confined to GAD65/67+ GABAergic neurons, whereas expression in the Dm revealed colocalization with vglut2b+ glutamatergic neurons. Whether sam2-expressing neurons employ GABA or glutamate as a second transmitter remains to be determined (55). Also, whether sam2 regulates the neuronal activity of GABAergic or glutamatergic neurons by neuromodulatory action will need to be ascertained in the future. Nevertheless, the expression and colocalization patterns reported in our current study open the possibility of a neuromodulatory role of sam2 in both excitatory and inhibitory neurotransmitter systems comprising anxiety circuitry.

The specific molecular mechanisms underlying increased anxiety-like behavioral responses observed in sam2 KO fish will require further characterization. Nevertheless, the behavioral changes we found appeared to be unrelated to the serotonergic neurotransmitter system, as the expression of indicative markers, such as tphR and slc6a4a, were not affected by the loss of sam2 () (29, 56). However, after the novel tank assays, sam2 KO fish did show higher level of c-fos expression in the entire brain compared with wild-type (), confirming their abnormal response to novelty but also implying potential involvement of neurotransmitters other than the glutamatergic and GABAergic systems.

Finally, comparative genomic mapping with microdeletions and microduplications in human patients with intellectual disability, autism, and behavioral tantrums has identified SAM2 as a candidate gene for these phenotypes in individuals with 12q14.1 CNVs (Fig. 7). Although the latter results do not prove the role of SAM2 in these diseases, these results, together with our findings with zebrafish and mice, present an interesting possibility: evolutionarily distant species, such as fish, rodents, and humans, may possess similar mechanisms underlying anxiety and emotion, and these mechanisms may include sam2. We thus propose that the discovery of the samdori gene family and the functional information of one of its member sam2, may have important translational relevance and therapeutic potential.

Materials and Methods

Full methods and any associated references are available in .

Animals.

All experimental protocols and procedures were approved and conducted according to the approved guidelines and regulations of the Animal Ethics Committee of Chungnam National University (CNU-00866).

Isolation of samdori Gene Family.

We established a unique insertional mutagenesis system based on the Sleeping Beauty transposon, called the golden fish project (). In addition to the GFP gene, we used a melanin-concentrating hormone gene as a transgene reporter for visual screening of mutants. During the course of mutagenesis, we mapped the insertion site in a chemokine gene, named samdori (sam). The full-length cDNAs of zebrafish sam-1a, 1b, 2, 3a, 3b, 4, 5a, and 5b genes were isolated using the RACE technique. The sam2 KO zebrafish lines (sam2 and sam2) and Sam2 KO mice were generated by the targeted mutagenesis utilizing ZFNs (26, 27) and transcription activator-like endonucleases (57), respectively.

Behavior Tests in Zebrafish and Mice.

Male zebrafish siblings of identical age (3-mo-old) and size were used for the following behavioral tests: open novel tank (35, 36, 41), scototaxis (37), and shoaling behavior (1, 42). The open-field test, elevated plus maze, and fear conditioning were employed in mice. All values represent mean ± SEM.

Human Genetics Data.

Patients GC42855 and GC48823 were identified among a cohort of ∼32,700 patients undergoing clinical, oligonucleotide-based array comparative genomic hybridization (aCGH) (58, 59) at Signature Genomic Laboratories. The DECIPHER database is a publicly accessible web-based resource of human genome variants with associated phenotypes, which facilitates the identification and interpretation of pathogenic genetic variation in thousands of patients with rare disorders worldwide (47). The description of case 290951 was ascertained from the DECIPHER database. We tested the deletion in 290951 family members; II-2 (affected, has deletion), III-2 (mother, affected, has deletion), III-3 (uncle, unaffected, does not have deletion), IV-1 (elder sister, unaffected, no deletion), IV-2 (elder brother, unaffected, no deletion), IV-3 (the proband, affected, has the deletion). This study was approved by the Institutional Review Board of Augusta University, Augusta, Georgia (HAC 0904264) for collection of deidentified phenotypic and genotypic data of human subjects or, with informed consent, for additional studies completed on blood or DNA. Informed consent, blood, and cells, were obtained from the family of 290951, while deidentified data were used for the rest of the analyses.