- Open Access
Forebrain excitatory neuron-specific SENP2 knockout mouse displays hyperactivity, impaired learning and memory, and anxiolytic-like behavior
Molecular Brain volume 13, Article number: 59 (2020)
Sentrin/SUMO-specific protease 2 (SENP2) is a member of SENPs family involved in maturation of SUMO precursors and deSUMOylation of specific target, and is highly expressed in the central nervous system (CNS). Although SENP2 has been shown to modulate embryonic development, fatty acid metabolism, atherosclerosis and epilepsy, the function of SENP2 in the CNS remains poorly understood. To address the role of SENP2 in the CNS and its potential involvement in neuropathology, we generated SENP2 conditional knockout mice by crossing floxed SENP2 mice with CaMKIIα-Cre transgenic mice. Behavioral tests revealed that SENP2 ablation induced hyper-locomotor activity, anxiolytic-like behaviors, spatial working memory impairment and fear-associated learning defect. In line with these observations, our RNA sequencing (RNA-seq) data identified a variety of differential expression genes that are particularly enriched in locomotion, learning and memory related biologic process. Taken together, our results indicated that SENP2 plays a critical role in emotional and cognitive regulation. This SENP2 conditional knockout mice model may help reveal novel mechanisms that underlie a variety of neuropsychiatric disorders associated with anxiety and cognition.
SUMOylation is a dynamic and reversible post-translational modification that modulates diverse functions of target proteins, including protein stability, protein subcellular localization, protein-protein or protein-DNA interactions, and protein kinase activity . Cellular abundance of particular SUMO-conjugated substrates is regulated by a balance between SUMO conjugation and SUMO deconjugation. Sentrin-specific proteases (SENPs) catalyze the removal of SUMO from SUMO-conjugated target proteins as well as the cleavage of SUMO from its precursor proteins, thus playing a critical role in regulating the SUMOylation level of targets [2, 3]. In mammals, the SENP family consists of six members, which can be divide into three groups (SENP1 and SENP2; SENP3 and SENP5; SENP6 and SENP7) based on homology and function analysis . We recently showed that SENP1 participates in regulating nociceptive signaling in models of inflammatory pain and attenuates I/R injury induced cell death in a transient brain ischemia/reperfusion mouse model [4, 5]. Moreover, SENP2 has been reported to play a role in cardiac development , neuronal survive  and seizure . According to the Allen Brain Atlas, SENP2 mRNAs are highly expressed in the forebrain , but the function of SENP2 in the central nervous system (CNS) remains unclear. As SENP2 is required for expression of key developmental genes, global deletion of SENP2 is embryonically lethal [6,7,8, 10]. Thus, we developed a forebrain excitatory neuron-specific SENP2 knockout mouse model to examine SENP2 functions in the CNS. We found that these conditional knockout (cKO) animals display hyperactivity and reduced anxiety-like behavior, impaired learning and memory. Gene ontology (GO) analysis of the differential expression genes revealed enrichment for numerous cellular and molecular functional categories, including those related to “Cell death” and “Immune response”. Taken together, our results indicate that SENP2 plays an important role in the forebrain, and its absence leads to molecular and behavioral changes associated with locomotion, anxiety, learning and memory.
Generation of forebrain-specific SENP2 cKO mice
To determine the functional role of SENP2 in the CNS, we crossed floxed SENP2 (SENP2fl/fl) mice  with CaMKIIα-Cre transgenic mice  to generate SENP2 conditional knockout (cKO) mice, in which SENP2 was selectively removed from principal neurons of postnatal forebrain (Fig. 1a-b). Because the CaMKIIα-Cre transgene is expressed between postnatal days 14–21 in excitatory neurons in the forebrain , this allowed us to specifically assess SENP2 function in the postnatal forebrain without disrupting its contribution to early CNS development and/or causing embryonic lethality. In the cKO mice, SENP2 protein levels were strongly reduced in the forebrain (cortical and hippocampal) excitatory neurons (Fig. 1c). As determined by western blot analyses, loss of SENP2 protein expression occurred in the cortex (58.01 ± 4.90% of littermate controls) and hippocampus (42.96 ± 6.24% of littermate controls) of 6-week-old cKO mice, but not in the cerebellum (105.66 ± 9.99% of littermate controls), where the Cre recombinase is not expressed (Fig. 1d). Additionally, real-time quantitative PCR (qPCR) determined that SENP2 mRNA levels in the cKO mice were significantly reduced in cortex (45.82 ± 2.94% of littermate controls) and hippocampus (35.17 ± 6.07% of littermate controls) of 6-week-old cKO mice. However, no change was observed in the cKO mice at 3 weeks after birth (Fig. 1e).
Behavioral analyses of cKO mice
We performed a battery of behavioral tests to evaluate the behavioral phenotype of SENP2 cKO mice. In open field test, cKO mice travelled more distance (m) than their littermate controls during a 30 min recording period (Ctrl (Control): 94.66 ± 4.43, n = 19; cKO: 216.26 ± 15.81, n = 12; p < 0.0001, Welch’s t-test) (Fig. 2a-c). As abnormal exploratory behaviors may also be indicative of changes in anxiety, we assessed the number of entries, duration spent and distance traveled in center area. We observed that SENP2 cKO mice entered the center more frequently (Ctrl: 68.11 ± 6.00, n = 19; cKO: 94.75 ± 6.48, n = 12; p = 0.0068, Student’s t-test) (Fig. 2d-e), spent more time (s) in center area (91.64 ± 9.97, n = 19; cKO: 149.04 ± 11.42, n = 12; p = 0.0009, Student’s t-test) (Fig. 2f-g) and traveled more distance in center area than littermate controls (Ctrl: 9.93 ± 0.95, n = 19; cKO: 13.85 ± 1.00, n = 12; p < 0.0107, Student’s t-test) in a 30 min test session (Fig. 2h-i). Although a 5 min test session is often sufficient to assess the critical components of general exploratory locomotion, the most commonly used measure of overall exploratory/locomotor activity is the total distance traveled. SENP2 cKO mice spent longer time and entered more frequency in the center area than littermate control mice at 10 min, 15 min, 20 min(Fig. 2d, f), and traveled more distance at 10 min, 20 min (Fig. 2h), although no statistically significant difference between cKO and littermate control mice was observed at 5 min, 25 min and 30 min. These findings suggest that cKO mice showed great aspiration for exploring the center area, and thus displayed decreased anxiety-like behavior.
To further explore the anxiety-like behaviors, cKO mice were evaluated in two anxiety-related behavioral assays including elevated plus maze and novelty suppressed feeding. In the evaluated plus maze test, cKO mice gained more number of entries in the open arms (Ctrl: 6.84 ± 0.97, n = 19; cKO: 21.58 ± 3.91, n = 12; p = 0.0031, Welch’s t-test) (Fig. 2k) and less number of entries in the closed arms accordingly (Ctrl: 25.26 ± 1.32, n = 19; cKO: 19.00 ± 2.49, n = 12; p = 0.0213, Student’s t-test) (Fig. 2l). In addition, cKO mice spent more percentage of time in the open arms (Ctrl: 6.85 ± 1.54%, n = 19; cKO: 37.20 ± 6.90%, n = 12; p = 0.0010, Welch’s t-test) (Fig. 2m) and less percentage of time in the closed arms (Ctrl: 68.31 ± 2.94%, n = 19; cKO: 43.36 ± 7.81%, n = 12; p = 0.0096, Welch’s t-test) (Fig. 2n). Moreover, cKO mice significantly reduced the latency (s) to feed in the novel environment compared with littermate controls (Ctrl: 294.79 ± 36.24, n = 19; cKO: 89.42 ± 38.73, n = 12; p = 0.0008, Student’s t-test) (Fig. 2o). Taken together, these results demonstrated that cKO mice exhibited anxiolytic-like behavior.
Previous studies have suggested that SUMOylation plays an important role in learning and memory [12,13,14]. Given the critical function of SENP2 in regulating SUMOylation status and its high-level expression in the brain, we investigate whether the loss of SENP2 in forebrain excitatory neurons impairs learning and memory. As first, we assessed spatial working and reference memory of cKO mice using the Y-maze spontaneous alternation task. Compared with littermate controls, cKO mice displayed significantly reduced alternations (Ctrl: 72.18 ± 1.84%, n = 18; cKO: 60.35 ± 4.22%, n = 10; p = 0.0062, Student’s t-test) (Fig. 3a), suggesting an impairment of spatial working memory in cKO mice. Moreover, we examined associative learning and memory behaviors using a contextual fear conditioning protocol . As shown in Fig. 3b, littermate control mice responded well in training session and exhibited freezing behavior when reintroduced to the same context 24 h later. In contrast, cKO mice not only showed impaired learning ability during training, but also displayed extremely low freezing level in contextual test session, demonstrating an impairment of cKO mice in associative learning and memory. Furthermore, we found that the cKO mice had a substantially decreased nesting score in the nest building assay compared with littermate control mice (Ctrl: 4.30 ± 0.26, n = 10; cKO: 1.44 ± 0.24, n = 9; p < 0.0001, Student’s t-test) (Fig. 3c). In rodents, the nest building behavior represents a form of social behaviors, and that impaired nest building is considered to represent a negative phenotype of psychiatric diseases including schizophrenia . To determine whether cKO mice displayed other core characteristics of the neuropsychiatric disorders, we conducted prepulse inhibition (PPI) task to measure the sensorimotor gating of cKO mice . Compared with their littermates, cKO mice displayed normal startle reaction (Ctrl: 475.71 ± 78.32, n = 8; cKO: 522.48 ± 58.08, n = 8; p = 0.6389, Student’s t-test) (Fig. 3d). When assayed in a prepulse inhibition (PPI) task, there were no differences between cKO and littermate control mice in the extent of PPI at 3 increasing prepulse sound intensities (%PPI of 70 dB prepulse stimulus: Ctrl: 47.64 ± 5.23%, n = 8; cKO: 42.90 ± 3.01%, n = 8; p > 0.9999; %PPI of 74 dB prepulse stimulus: Ctrl: 47.86 ± 6.92%, n = 8; cKO: 41.00 ± 3.84%, n = 8; p = 0.9895; %PPI of 78 dB prepulse stimulus: Ctrl: 47.63 ± 6.44%, n = 8; cKO: 48.23 ± 2.21%, n = 8; p > 0.9999; two-way ANOVA analysis followed by Bonferroni’s post-hoc) (Fig. 3d). Thus, cKO mice had no obvious deficit in PPI.
Identification of SENP2-regulated transcripts in the cerebral cortex
To elucidate the molecular mechanism underlying the behavioral deficits exhibited in the cKO mice, we performed high-throughput RNA sequencing (RNA-seq) to identify genes with altered expressions when SENP2 is selectively removed in the forebrain excitatory neurons. RNAs were prepared from cerebral cortex tissues isolated from the brains of 6-week-old cKO and their littermate control mice. We sequenced RNA libraries from 3 biological replicates per genotype and evaluated the data by pearson correlation coefficient. Differential expression genes (DEGs) analysis revealed consistent changes between the two genotypes across all 3 replicates (Fig. 4a). We observed 863 up-regulated genes and 170 down-regulated genes in cerebral cortex of cKO mice using an adjusted p value < 0.05 and relative gene expression level > 2 fold change (Fig. 4b). A detailed comparative analysis of the gene expression profiles appears in Additional file 1: Table 1.
GO and KEGG enrichment analysis for DEGs
Gene ontology (GO) analysis of the DEGs revealed SENP2 ablation altered the expression of genes involved in several biological processes, including locomotion, learning and memory (Fig. 5a). Additionally, GO analysis of DEGs revealed enrichment for numerous cellular and molecular functional categories, including those related to “Cell death” (Fig. 5b) and “Immune response” (Fig. 5c). A detailed GO enrichment analysis of the related GO terms is included in the Additional file 1: Table 2. Together, these results suggest that DEGs between cKO and littermate control mice are enriched for locomotion, learning and memory, cell death risk related genes. To further explore the molecular signaling pathway related to behavioral phenotypes, we introduce KEGG (Kyoto Encyclopedia of Genes and Genomes) enrichment analysis for DEGs (Fig. 5d). A detailed KEGG enrichment analysis of the KEGG signaling pathways appears in the Additional file 1: Table 3. KEGG enrichment analysis showed that “Immune system” related pathways were significantly enriched such as “Inflammatory bowel disease (IBD)”, “Th17 cell differentiation”, “Hematopoietic cell lineage”, “Cytosolic DNA-sensing pathway” and so on. These results were in accordance with the results of the GO enrichment analysis of “Immune response related biologic process” (Fig. 5c). On the other hand, we observed that “Cell growth and death” related pathways were also enriched such as “Cellular senescence” and “P53 signaling pathway”. These pathways may involve in modulating the programmed cell death of cKO mice [18,19,20]. Previous studies revealed that “MAPK signaling pathway” involved in regulating anxiety and depression-like behavior in mice [21, 22]. Our observation also showed that SENP2 ablation decrease anxiety-like behavior (Fig. 2a-o). Based upon the enriched KEGG pathways, previous studies and our findings, the “MAPK signaling pathway” could be one of a possible mechanism in modulating anxiety-like behavior in cKO mice.
Multiple gene expression changes associated with anxiety-related behavior
As the SENP2 cKO mice exhibit decreased anxiety-like behavior (Fig. 2), we attempted to screen out the genes associated with anxiety-like behavior from DEGs. However, to the best of our knowledge, there are no available databases that record anxiety disorder related genes. Thus, we decided to retrieve the 1033 DEGs using the PubMed database. By searching for the literature, we identified that 17 genes associated with anxiolytic-like effect on behavior are up-regulated (Fig. 6a), and 9 genes associated with anxiety-like effect on behavior are down-regulated (Fig. 6b). In addition, the 26 genes of anxiety-like or anxiolytic-like effect on behavior have been further confirmed by transgenic mice. Taken together, these genes may have an important role in modulating the level of anxiety. The 17 genes associated with anxiolytic-like effect on behavior, including that APOE , TSPO , UCP2 , MT1 , MT2 , AIM2 , CNR2 , LCN2 , DLK1 , CNTF , MYD88 , HDC , ALOX5 , TLR4 , ROCK1 , MOV10 , DDC . The 9 genes associated with anxiety-like effect on behavior, including that LAMP5 , KALRN , CRHR1 , CCK , SNAP25 , OLFM2 , FAAH , IQSEC2 , EMX1 . These data suggest that either the anxiolytic-like effects gene up-regulated or anxiety-like effects gene down-regulated could be a reduced level of anxiety. However, to elucidate the detailed mechanism of decreased anxiety-like behavior, further studies are needed to screen and determine from the 26 candidates and confirm its role in regulating the level of anxiety in cKO mice.
Animal models are extremely useful in establishing causality between genetic mutations, synaptic changes, circuit dysfunctions and abnormal behaviors, thereby aiding us in understanding the pathogenesis of neurological and neuropsychiatric diseases [49, 50]. SENP2 is a member of the sentrin/SUMO-specific proteases (SENPs) family and implicated in embryonic development [6, 10], fatty acid metabolism , atherosclerosis , and epilepsy . However, the function of SENP2 in CNS and its potential contribution to the neuropathology remains unclear. In this study, we generated conditional knockout of SENP2 in excitatory neurons in the postnatal forebrain and determined that these cKO mice display comprehensive behavioral phenotypes including hyperactivity, reduced anxiety-like behavior, impaired learning and memory (Figs. 2 and 3). Consistently, RNA-seq results showed that the loss of SENP2 is associated with moderate changes in gene transcripts related to “locomotion”, “learning and memory” (Fig. 5a) and multiple gene expression changes associated with anxiety-related behavior (Fig. 6). In addition, GO enrichment analysis identified changes in genes related to “Cell death” and “Immune response” (Fig. 5b-c). Taken together, our results demonstrate that SENP2 plays important functional roles in the forebrain.
SUMOylation is a dynamic and reversible posttranslational protein modification that regulates the functions of target proteins . In the CNS, neuron-specific SUMO1–3 knockdown mice show less exploration of center area in open field test . Additionally, previous studies showed that hippocampus-dependent learning and memory is impaired by overexpression of a dominant negative Ubc9 peptide in the hippocampal CA1 area . These results suggested that the balance of SUMOylation/deSUMOylation is critical for forebrain mediated function, including hippocampus-dependent learning and memory, and anxiety-like behaviors. SENPs catalyze the removal of SUMO from SUMO-conjugated proteins, thus playing a critical role in regulating the SUMOylation level of targets [2, 3]. Previous study revealed that neuron-specific SENP2 knockout mice display hyperactivity and sudden death . Consistently, we report here that the forebrain excitatory neuron-specific SENP2 knockout mice displayed hyper-locomotor activity in the open field test. Moreover, using several different behavioral assays, we identified a reduced anxiety phenotype in the SENP2 cKO mice (Fig. 2j-o).
Previous study has reported that death of matured neurons in the forebrain increases the level of anxiety . To further explore the molecular mechanism related to behavioral phenotypes of cKO mice, we conducted RNA-seq on the cortex of cKO animals and observed that numerous DEGs related to programmed cell death biologic processes were enriched (Fig. 5b), suggesting that ablation of SENP2 may lead to neuronal death. However, we don’t yet have any direct evidence showing that conditional knockout of SENP2 led to programmed cell death in the forebrain excitatory neurons. Further study is certainly needed to investigate the potential linkage between neuronal death and anxiolytic-related behaviors in the SENP2 cKO mice. On the other hand, we identified 26 genes from 1033 DEGs that are implicated in anxiety-related behaviors. These include 17 up-regulated genes known to have anxiolytic-like effects and 9 down-regulated genes have anxiety-like effects (Fig. 6a-b). This finding thus should shed some light on the potential molecular mechanism for the anxiolytic behaviors exhibited by the SENP2 cKO mice.
The SENP2 cKO mice also exhibited deficits in associative learning and spatial working memory functions, which were assessed in both contextual fear conditioning test and Y-maze spontaneous alternation task respectively (Fig. 3a-b). Consistent with the observed cognitive impairments, the RNA-seq studies we conducted on the cortex of cKO animals detected 20 DEGs that are specifically related to learning and memory-related tasks (Fig. 5a). These findings suggest that SENP2-dependent protein modifications are important for learning and memory, and dysfunction of SENP2 contributes to impaired cognitive functions in mice. Moreover, GO analysis of the DEGs in the SENP2 cKO mice revealed enrichment of genes for several cellular and molecular functional categories, including those related to “Cell death” and “Immune response” (Fig. 5b-c). Given that SENP2 has been implicated in neural disorders such as epilepsy and neurodegeneration [7, 8], future studies should be directed to examine the role of these SENP2-regulated molecular pathways in the observed behaviors.
As the subcellular localization of SENP2 is predominantly nuclear, we performed the differential analysis of gene and transcript expression using high throughput RNA-seq. These analyses suggest that SENP2 knockout results in changes of gene transcription that are associated with various behaviors including locomotion, learning and memory, anxiety. One interpretation of these results is that SENP2 dependent post translational modification of proteins, such as transcription factors and co-factors, may lead to alterations in gene expression. However, it is unclear whether these changes in gene transcription are directly caused by the loss-of-function of SENP2. On the other hand, it is very likely that SENP2 can directly regulate neuronal functions by altering the status of SUMOylated proteins. We hope that, in future studies, different proteomics approaches will help further gain insight into the role of SENP2 plays in the CNS.
In summary, we generated a line of cKO mice to investigate the functions of SENP2 in the CNS. Our results demonstrate that dysfunction of SENP2 in forebrain excitatory neurons leads to behavioral changes associated with emotional and cognitive functions.
Materials and methods
All animal procedures were carried out in accordance with the guidelines for care and use of laboratory animals of Shanghai Jiao Tong University School of Medicine and approved by the Institutional Animal Care and Use Committee (IACUC). All mice used in this study were C57BL/6 background. SENP2fl/fl mice  were crossed with CaMKIIα-Cre mice  to generate SENP2 conditional knockout mice (cKO). PCR primers used for genotyping are the following: SENP2 loxP Forward: 5′-CTTCTGCTTCTCTTAGTGCT-3′ and SENP2 loxP Reverse: 5′-CTTCTGCTTCTCTTAGTGCT-3′, with the expected product sizes of 183 and 149 bps for SENP2fl/fl and WT mice respectively. The presence of CaMKIIα-Cre was identified by PCR using the following primers: Cre Forward: 5′-CGCTGGGCCTTGGGACTTCAC-3′ and Cre Reverse: 5′-CAGCATTGCTGTCACTTGGTC-3′, with the expected PCR product size of 403 bps.
Immunofluorescence staining was performed as described previously , with minor modifications. Mice were perfused intracardially with 4% paraformaldehyde. After an overnight postfixation in the same fixative at 4 °C, brain tissues were embedded in 2% agar and cut into 50 μm sections with a vibratome (Leica, VT1200S). Brain sections were blocked with 0.3% Triton X-100 and goat serum in PBS for 1 h at room temperature and then incubated with SENP2 (1:100, Abcam, ab58418) and CaMKIIα (1:400, ThermoFisher Scientific, MA1–048) primary antibody overnight at 4 °C. Next, brain sections were incubated with Alexa Fluor546-conjugated goat anti-rabbit (Invitrogen, 1:500, A-11035) or Alexa Fluor488-conjugated goat anti-mouse secondary antibodies (Invitrogen, 1:500, A-11029) for 1 h at room temperature and mounted on glass slides using a small brush. The fluorescence images were acquired using the confocal microscope (Leica, TCS SP8).
Western bolt was performed as previously described with minor modifications . Brain tissues were collected on ice and homogenized with lysis buffer I (50 mM Tris-HCl, pH 6.8, 2% SDS, 40 mM DTT, and 5% glycerol). After denaturation for 15 min at 95 °C, the samples were diluted 10-fold with lysis buffer II (50 mM Tris-HCl, pH 7.4, 150 mM NaCl and 1% NP-40) and ultra-sonicated for 5 s, followed by centrifugation for 10 min at 13,000×g (4 °C). The supernatant was transferred to a new tube and boiled with loading buffer for 15 min. The proteins were separated by SDS-PAGE, transferred to polyvinylidene fluoride (PVDF) membranes, blocked with 5% non-fat milk, and immunoblotted with SENP2 (1:1000, Santa Cruz, sc-376,731) and β-actin (1:10000, Santa Cruz, sc-130,065) antibodies.
Real-time quantitative polymerase chain reaction (qPCR)
qPCR was performed as previously described with minor modifications . Brain tissues (cortex, hippocampus and cerebellum) were isolated on ice from cKO and littermate control mice. The total RNAs were extracted from these brain tissues using TRIzol Reagent (Tiangen, Beijing, China). After removing genomic DNAs, the total RNAs were reverse transcribed to complementary DNA (cDNA) using the PrimeScript RT reagent Kit (Takara, Dalian, China) according to the manufacturer’s protocol. The cDNAs were then used as templates for the qPCR reactions, which were performed in a 10 μl volume with Power SYBR Green PCR Master Mix (CWBIO, Beijing, China) and 0.2 μM primers using the LightCycler 480 real-time PCR system (Roche, CA, USA). The qPCR primers used for cKO and littermate control mice were as follows: SENP2 forward: 5′-TTCTCGGCACCATTCTTCGCTTGT-3′, SENP2 reverse: 5′-TGCTGCAGGATCCAGAACTCATCA-3′. GAPDH forward: 5′-CATGGCCTTCCGTGTTCC-3′ and GAPDH reverse: 5′-GCCTGCTTCACCACCTTCTT-3′.
RNA sequencing (RNA-seq)
RNA extraction: Three pairs of cortical tissue were isolated on ice from 6-week-old cKO and littermate control mice and quickly frozen in liquid nitrogen. The samples were delivered to the company (Majorbio Bio-pharm Technology, Shanghai, Chain) to prepare RNA samples and conduct high-throughput RNA-seq. The RNA quality was determined by 2100 Bioanalyser (Agilent Technologies, Tokyo, Japan) and quantified using the ND-2000 (NanoDrop Technologies, DE, USA). The RNA integrity was also confirmed with agarose gel electrophoresis.
Library construction and sequencing: 5 μg high-quality RNA samples were used to construct sequencing library. Briefly, messenger RNA was isolated from total RNA samples and then fragmented by fragmentation buffer firstly. Secondly, double-stranded cDNA was synthesized using a SuperScript double-stranded cDNA synthesis kit (Invitrogen, CA, USA) with random primers (Illumina, CA, USA). Thirdly, the synthesized cDNA was subjected to end-repair, phosphorylation and ‘A’ base addition according to Illumina’s library construction protocol. Then, 200–300 bp cDNA target fragments were isolated using 2% agarose followed by PCR amplified, the isolated cDNA fragments were selected to construct sequencing library. Libraries quantified by TBS380, and sequenced with the Illumina HiSeq 4000 (Illumina, CA, USA).
Read mapping: The raw paired end reads were trimmed by “SeqPrep” software and quality controlled by “Sickle software. Then clean reads were separately aligned to reference genome with orientation mode using “TopHat” software. The criteria for mapping were that sequencing reads should be uniquely matched to the genome with less than 3 mismatches, without insertions or deletions. Then, the gene regions were expanded following depths of sites and the operon was obtained. In addition, the whole genome was segmented into multiple 15 kb fragments that share the same 5 kb fragments. If more than 2 consecutive fragments were without overlapped region and at least 2 reads mapped per fragment in the same orientation, we consider it was new transcribed region of gene .
Data analysis: To identify differentially expressed genes (DEGs) between cortical tissues of cKO and littermate control mice, each transcript expression level was calculated according to the fragments per kilobase of transcript per million mapped reads (FPKM) method. The gene abundances were quantified using the RSEM software, while the EdgeR software was used for differential expression analysis and quantified transcript read counts. In addition, the differential expression analysis was carried out on an online platform (www.majorbio.com) followed by multiple check calibration (BH). A differentially expressed gene is identified as relative read counts > 2 fold change and adjusted p value < 0.05. Similarly, Go enrichment analysis was conducted using the “Goatools” software followed by Fisher’s exact test. KEGG enrichment analysis was using the “KOBAS 2.0” software followed by Fisher’s exact test. If the adjusted p value < 0.05, we consider the GO terms or KEGG signing pathways were significantly enriched.
Experimental design for behavioral tests
Behavioral tests were performed on male mice of 8 to 12-weeks of age. Mice were housed in a room with 12 h light/dark circadian rhythm, suitable temperature (22–28 °C), adequate water and food. The behavioral tests were performed in the light-on phase of the cycle (10:00 a.m.-18:00 p.m.) except for otherwise noted. The experimenters were blinded to the genotype of each mouse during all tests and data analyses. The tests were performed in the following sequence: open field, elevated plus maze, novelty suppressed feeding, nest building, Y maze, and contextual fear condition. The startle response/prepulse inhibition test was performed using another set of male mice without any stressors. Tests were repeated using at least two different cohorts of mice. These tests were performed at intervals of 2–4 days.
Open field test
The open field test was performed as previously described with minor modifications . Mice were habituated in the testing room for 60 min and then introduced to the open field apparatus (40 × 40 × 30 cm) (MED Associates). The test mice were allowed to freely explore for 30 min in apparatus without interference. The distance traveled, number of entries in center area (20 × 20 cm), time spent in center area and distance traveled in center area were automated recorded by monitor system and software (EthoVision XT 12).
Elevated plus maze
Elevated plus maze test was performed as previously described . Mice were habituated in the testing room for 60 min and placed in the intersection of open and closed arms with the mouse head toward the open arm. The test mice were allowed to freely explore for 5 min in the apparatus (MED Associates) without interference. Monitor system and software (EthoVision XT 12) automated recorded the time that mice spent in open and closed arms as well as the number of entries in open and closed arms respectively.
Novelty suppressed feeding
Novelty suppressed feeding test was performed as previously described  with minor modifications. Briefly, mice were fasted for 24 h in the home cage before testing and then placed in a new feeding box in which food was fixed on one piece of round filter paper (10 cm in diameter) at the center of the plastic chamber (30 × 40 cm). The mice were allowed to freely explore and eat the food for 5 min. A monitor system (EthoVision XT 12) recorded the process of food eating including the latency to feeding.
The Y maze has three identical opaque arms (40 cm-long, 10 cm-wide, and 15 cm-high). We conducted the Y maze test as previously described . Briefly, mice were habituated the testing room for 60 min and placed in the distal end of one arm. The mice were allowed to freely explore for 5 min without any interference. The monitor system (EthoVision XT 12) recorded the locus of mice movement, while the sequence of entries was manually recorded and analyzed by the experimenter. The alternation ratio was calculated as follows:
Contextual fear condition
Contextual fear condition was conducted as previously described with minor modifications . Briefly, mice were habituated in the testing room for 60 min and then placed in a test chamber with a black and white plaid sticker on all sides for 10 min to habituate the apparatus (Ugo Basile) that constantly presents 100 lx bright light for 2 consecutive days. In the training session, mice were placed in the test chamber with a metal grid floor and received footshock (0.5 mA, 2 s) for 5 trials with 120 s interval. After the footshock, mice remained in the conditioning chamber for another 30 s and then were placed back to their home cage for 24 h. In the retrieval session, mice were placed into the same contextual chamber for 10 min test without footshock. The freezing values were recorded every 2-min block by the monitor system and software (ANY-maze).
Nest building test was performed as previously described with minor modifications . Briefly, mice were individually housed for 12 h (20:00 p.m.-8:00 a.m.) while providing a piece of cleansing paper into the cage. The next day, the nesting score was assessed based on the integrity of the paper by an experimenter blinded to mouse genotypes . Nests were given a score of 0–5 according to the following criteria: 1 = more than 90% of paper was intact; 2 = 50–90% of paper remained intact; 3 = more than 50% of paper was torn, but no identifiable nest site; 4 = more than 90% of paper was torn and a flat nest was built; 5 = more than 90% of paper was torn and the paper was transformed into a tridimensional nest.
Prepulse inhibition (PPI)
Prepulse inhibition test was carried out as previously described  with minor modifications. Mice were habituated to the testing room for 60 min, and placed into the test chamber for 5 min to acclimate the apparatus (MED Associates) that present the constant background white noise of 65 dB for 2 consecutive days. In stage I of the PPI session, we replaced the mice into the chamber for a 5 min acclimation. In stage II of PPI session, we presented 10 trials of high acoustic stimulus (120 dB) with an interval of 20 ms to make the mice accommodate the high acoustic stimulus. In stage III of PPI session, we random presented seven types of acoustic stimulus including: 1) the high acoustic stimulus (120 dB) only or 2) the low acoustic stimulus only (70 dB, 74 dB and 78 dB) or 3) the high acoustic stimulus paired with a low acoustic stimulus (70 dB paired with 120 dB or 74 dB paired with 120 dB or 78 dB paired with 120 dB) with a random interval of 30–100 ms. In our protocol, six blocks containing seven acoustic stimulus types were presented in a pseudorandom order and each acoustic stimulus type was presented once within a block. The acoustic reflex was measured by SR-Lab system (San Diego Instruments) and counted using the largest peaks of amplitude that the signal recorded in a 300 ms window after presented the acoustic stimulus. We respectively counted the startle amplitude of high acoustic stimulus (120 dB only) and the paired acoustic stimulus (70 dB paired with 120 dB or 74 dB paired with 120 dB or 78 dB paired with 120 dB). The startle amplitude of 120 dB stimulus represented the acoustic reflex of mice, and was used as the baseline value of PPI. The PPI was calculated as follows:
All data are presented as mean ± S.E.M and analyzed by GraphPad Prism8 software. Two-group comparison was processed with two-tailed, unpaired Student’s t test when the variance is equal (Fig. 2e, g, i, l, o, Fig. 3a, c) or Welch’s t-test when the variance is unequal (Fig. 2c, k, m, n). Multiple group comparison was performed with two-way ANOVA analysis followed by Bonferroni’s post-hoc to determine significance (Fig. 1d-e, Fig. 2b, d, f, h, Fig. 3b, d). *p < 0.05, **p < 0.01, ***p < 0.001 compared with littermate controls.
Data availability statement
The data that support the findings of this study are available from the corresponding author upon reasonable request.
Wilkinson KA, Nakamura Y, Henley JM. Targets and consequences of protein SUMOylation in neurons. Brain Res Rev. 2010;64(1):195–212.
Schorova L, Martin S. Sumoylation in synaptic function and dysfunction. Front Synaptic Neurosci. 2016;8:9.
Henley JM, Craig TJ, Wilkinson KA. Neuronal SUMOylation: mechanisms, physiology, and roles in neuronal dysfunction. Physiol Rev. 2014;94(4):1249–85.
Wang Y, Gao Y, Tian Q, Deng Q, Wang Y, Zhou T, et al. TRPV1 SUMOylation regulates nociceptive signaling in models of inflammatory pain. Nat Commun. 2018;9(1):1529.
Zhang H, Wang Y, Zhu A, Huang D, Deng S, Cheng J, et al. SUMO-specific protease 1 protects neurons from apoptotic death during transient brain ischemia/reperfusion. Cell Death Dis. 2016;7(11):e2484.
Kang X, Qi Y, Zuo Y, Wang Q, Zou Y, Schwartz RJ, et al. SUMO-specific protease 2 is essential for suppression of polycomb group protein-mediated gene silencing during embryonic development. Mol Cell. 2010;38(2):191–201.
Fu J, Yu HM, Chiu SY, Mirando AJ, Maruyama EO, Cheng JG, et al. Disruption of SUMO-specific protease 2 induces mitochondria mediated neurodegeneration. PLoS Genet. 2014;10(10):e1004579.
Qi Y, Wang J, Bomben VC, Li DP, Chen SR, Sun H, et al. Hyper-SUMOylation of the Kv7 potassium channel diminishes the M-current leading to seizures and sudden death. Neuron. 2014;83(5):1159–71.
Mouse Allen Brain Atlas. https://mouse.brain-map.org/experiment/show?id=276274.
Maruyama EO, Lin H, Chiu SY, Ivy Yu HM, Porter GA, Hsu W. Extraembryonic but not embryonic SUMO-specific protease 2 is required for heart development. Sci Rep. 2016;6:20999.
Minichiello L, Korte M, Wolfer D, Kühn R, Unsicker K, Cestari V, et al. Essential role for TrkB receptors in hippocampus-mediated learning. Neuron. 1999;24(2):401–14.
Wang L, Rodriguiz RM, Wetsel WC, Sheng H, Zhao S, Liu X, et al. Neuron-specific Sumo1-3 knockdown in mice impairs episodic and fear memories. J Psychiatry Neurosci. 2014;39(4):259–66.
Lee L, Dale E, Staniszewski A, Zhang H, Saeed F, Sakurai M, et al. Regulation of synaptic plasticity and cognition by SUMO in normal physiology and Alzheimer's disease. Sci Rep. 2014;4:7190.
Schwartz S, Truglio M, Scott MJ, Fitzsimons HL. Long-term memory in drosophila is influenced by histone deacetylase HDAC4 interacting with SUMO-conjugating enzyme Ubc9. Genetics. 2016;203(3):1249–64.
Phillips RG, LeDoux JE. Differential contribution of amygdala and hippocampus to cued and contextual fear conditioning. Behav Neurosci. 1992;106(2):274–85.
Jirkof P. Burrowing and nest building behavior as indicators of well-being in mice. J Neurosci Methods. 2014;234:139–46.
Mena A, Ruiz-Salas JC, Puentes A, Dorado I, Ruiz-Veguilla M, De la Casa LG. Reduced prepulse inhibition as a biomarker of schizophrenia. Front Behav Neurosci. 2016;10:202.
Childs BG, Baker DJ, Kirkland JL, Campisi J, van Deursen JM. Senescence and apoptosis: dueling or complementary cell fates? EMBO Rep. 2014;15(11):1139–53.
Haupt S, Berger M, Goldberg Z, Haupt Y. Apoptosis - the p53 network. J Cell Sci. 2003;116(20):4077–85.
Fricker M, Tolkovsky AM, Borutaite V, Coleman M, Brownet GC. Neuronal cell death. Physiol Rev. 2018;98(2):813–80.
Han Q, Huang H, Wang Y. Ghrelin exhibited antidepressant and anxiolytic effect via the p38-MAPK signaling pathway in hippocampus. Prog Neuro-Psychopharmacol Biol Psychiatry. 2019;93:11–20.
Wefers B, Hitz C, Hölter SM, Trümbach D, Hansen J, Weber P, et al. MAPK signaling determines anxiety in the juvenile mouse brain but depression-like behavior in adults. PLoS One. 2012;7(4):e35035.
Raber J. Role of apolipoprotein E in anxiety. Neural Plast. 2007;2007:91236.
Nothdurfter C, Rammes G, Baghai TC, Schüle C, Schumacher M, Papadopoulos V, et al. Translocator protein (18 kDa) as a target for novel anxiolytics with a favourable side-effect profile. J Neuroendocrinol. 2012;24(1):82–92.
Hermes G, Nagy D, Waterson M, Zsarnovszky A, Varela L, Hajos M, et al. Role of mitochondrial uncoupling protein-2 (UCP2) in higher brain functions, neuronal plasticity and network oscillation. Mol Metab. 2016;5(6):415–21.
Weil ZM, Hotchkiss AK, Gatien ML, Pieke-Dahl S, Nelson RJ. Melatonin receptor (MT1) knockout mice display depression-like behaviors and deficits in sensorimotor gating. Brain Res Bull. 2006;68(6):425–9.
Stefano C, Gabriella G. CCNP award paper: unveiling the role of melatonin MT2 receptors in sleep, anxiety and other neuropsychiatric diseases: a novel target in psychopharmacology. J Psychiatry Neurosci. 2014;39(1):6–21.
Wu PJ, Liu HY, Huang TN, Hsueh YP. AIM2 inflammasomes regulate neuronal morphology and influence anxiety and memory in mice. Sci Rep. 2016;6:32405.
Liu QR, Canseco-Alba A, Zhang HY, Tagliaferro P, Chung M, Dennis E, et al. Cannabinoid type 2 receptors in dopamine neurons inhibits psychomotor behaviors, alters anxiety, depression and alcohol preference. Sci Rep. 2017;7(1):17410.
Mucha M, Skrzypiec AE, Schiavon E, Attwood BK, Kucerova E, Pawlak R. Lipocalin-2 controls neuronal excitability and anxiety by regulating dendritic spine formation and maturation. Proc Natl Acad Sci. 2011;108(45):18436–41.
García-Gutiérrez MS, Navarrete F, Laborda J, Manzanares J. Deletion of Dlk1 increases the vulnerability to developing anxiety-like behaviors and ethanol consumption in mice. Biochem Pharmacol. 2018;158:37–44.
Peruga I, Hartwig S, Merkler D, Thöne J, Hovemann B, Juckel G, et al. Endogenous ciliary neurotrophic factor modulates anxiety and depressive-like behavior. Behav Brain Res. 2012;229(2):325–32.
Lim JE, Song M, Jin J, Kou J, Pattanayak A, Lalonde R, et al. The effects of MyD88 deficiency on exploratory activity, anxiety, motor coordination, and spatial learning in C57BL/6 and APPswe/PS1dE9 mice. Behav Brain Res. 2012;227(1):36–42.
Pittenger C. Histidine decarboxylase knockout mice as a model of the pathophysiology of tourette syndrome and related conditions. Handb Exp Pharmacol. 2017;241:189–215.
Joshi YB, Praticò D. Knockout of 5-lipoxygenase results in age-dependent anxiety-like behavior in female mice. PLoS One. 2011;6(12):e29448.
Femenia T, Qian Y, Arentsen T, Forssberg H, Diaz HR. Toll-like receptor-4 regulates anxiety-like behavior and DARPP-32 phosphorylation. Brain Behav Immun. 2018;69:273–82.
Greathouse KM, Henderson BW, Gentry EG, Herskowitz JH. Fasudil or genetic depletion of ROCK1 or ROCK2 induces anxiety-like behaviors. Behav Brain Res. 2019;373:11208.
Skariah G, Seimetz J, Norsworthy M, Lannom MC, Kenny PJ, Elrakhawy M, et al. Mov10 suppresses retroelements and regulates neuronal development and function in the developing brain. BMC Biol. 2017;15(1):54.
Lee NC, Chien YH, Hu MH, Liu WS, Chen PW, Wang WH, et al. Treatment of congenital neurotransmitter deficiencies by intracerebral ventricular injection of an adeno-associated virus serotype 9 vector. Hum Gene Ther. 2014;25(3):189–98.
Tiveron MC, Beurrier C, Céni C, Andriambao N, Combes A, Koehl M, et al. LAMP5 fine-tunes GABAergic synaptic transmission in defined circuits of the mouse brain. PLoS One. 2016;11(6):e0157052.
Mandela P, Yan Y, LaRese T, Eipper BA, Mains RE. Elimination of Kalrn expression in POMC cells reduces anxiety-like behavior and contextual fear learning. Horm Behav. 2014;66(2):430–8.
Marianne B, Stephan Z, Inge S, Hagemeyer TP, Deussing JM, Timpl P, et al. Limbic corticotropin-releasing hormone receptor 1 mediates anxiety-related behavior and hormonal adaptation to stress. Nat Neurosci. 2003;6:1100–7.
Del Boca C, Lutz PE, Le Merrer J, Koebel P, Kieffer BL. Cholecystokinin knock-down in the basolateral amygdala has anxiolytic and antidepressant-like effects in mice. Neuroscience. 2012;218:185–95.
Kataoka M, Yamamori S, Suzuki E, Watanabe S, Sato T, Miyaoka H, et al. A single amino acid mutation in SNAP-25 induces anxiety-related behavior in mouse. PLoS One. 2011;6(9):e25158.
Sultana A, Nakaya N, Dong L, Abu-Asab M, Qian H, Tomarev SI. Deletion of olfactomedin 2 induces changes in the AMPA receptor complex and impairs visual, olfactory, and motor functions in mice. Exp Neurol. 2014;261:802–11.
Moreira FA, Kaiser N, Monory K, Lutz B. Reduced anxiety-like behaviour induced by genetic and pharmacological inhibition of the endocannabinoid-degrading enzyme fatty acid amide hydrolase (FAAH) is mediated by CB1 receptors. Neuropharmacology. 2008;54(1):141–50.
Jackson MR, Loring KE, Homan CC, Thai MH, Määttänen L, Arvio M, et al. Heterozygous loss of function of IQSEC2/Iqsec2 leads to increased activated Arf6 and severe neurocognitive seizure phenotype in females. Life Sci Alliance. 2019;2(4):e201900386.
Cao BJ, Li Y. Reduced anxiety-- and depression-like behaviors in Emx1 homozygous mutant mice. Brain Res. 2002;937(1–2):32–40.
Kaiser T, Zhou Y, Feng G. Animal models for neuropsychiatric disorders: prospects for circuit intervention. Curr Opin Neurobiol. 2017;45:59–65.
Rothwell PE, Fuccillo MV, Maxeiner S, Hayton SJ, Gokce O, Lim BK, et al. Autism-associated neuroligin-3 mutations commonly impair striatal circuits to boost repetitive behaviors. Cell. 2014;158:198–212.
Koo YD, Choi JW, Kim M, Chae S, Ahn BY, Kim M, et al. SUMO-specific protease 2 (SENP2) is an important regulator of fatty acid metabolism in skeletal muscle. Diabetes. 2015;64(7):2420–31.
Heo KS, Chang E, Le NT, Cushman H, Yeh ET, Fujiwara K, et al. De-SUMOylation enzyme of sentrin/SUMO-specific protease 2 regulates disturbed flow-induced SUMOylation of ERK5 and p53 that leads to endothelial dysfunction and atherosclerosis. Circ Res. 2013;112(6):911–23.
Santos AL, Lindner AB. Protein posttranslational modifications: roles in aging and age-related disease. Oxidative Med Cell Longev. 2017;2017:5716409.
Wang AC, Jensen EH, Rexach JE, Vinters HV, Hsieh-Wilson LC. Loss of O-GlcNAc glycosylation in forebrain excitatory neurons induces neurodegeneration. Proc Natl Acad Sci. 2016;113(52):15120–5.
Gurusinghe CJ, Ehrlich D. Gelatin embedding of central nervous system tissue improves the quality of vibratome sections. Stain Technol. 1986;61(5):324–6.
Xie XL, Zhou WT, Zhang KK, Yuan Y, Qiu EM, Shen YW, et al. PCB52 induces hepatotoxicity in male offspring through aggravating loss of clearance capacity and activating the apoptosis: sex-biased effects on rats. Chemosphere. 2019;227:389–400.
Shoji H, Takao K, Hattori S, Miyakawa T. Age-related changes in behavior in C57BL/6J mice from young adulthood to middle age. Mol Brain. 2016;9:11.
Walf AA, Frye CA. The use of the elevated plus maze as an assay of anxiety-related behavior in rodents. Nat Protoc. 2007;2(2):322–8.
Cui W, Mizukami H, Yanagisawa M, Aida T, Nomura M, Isomura Y, et al. Glial dysfunction in the mouse habenula causes depressive-like behaviors and sleep disturbance. J Neurosci. 2014;34(49):16273–85.
Miedel CJ, Patton JM, Miedel AN, Miedel ES, Levenson JM. Assessment of spontaneous alternation, novel object recognition and limb clasping in transgenic mouse models of amyloid-β and tau neuropathology. J Vis Exp. 2017;123:55523.
Shoji H, Takao K, Hattori S, Miyakawa T. Contextual and cued fear conditioning test using a video analyzing system in mice. J Vis Exp. 2014;85:50871.
Ragot A, Pietropaolo S, Vincent J, Delage P, Zhang H, Allinquant B, et al. Genetic deletion of the histone Deacetylase 6 exacerbates selected behavioral deficits in the R6/1 mouse model for Huntington's disease. Brain Behav. 2015;5(9):e00361.
Deacon R. Assessing nest building in mice. Nat Protoc. 2006;1:1117–9.
Selcher JC, Nekrasova T, Paylor R, Landreth GE, Sweatt JD. Mice lacking the ERK1 isoform of MAP kinase are unimpaired in emotional learning. Learn Mem. 2001;8(1):11–9.
This study was supported by grants from the National Natural Science Foundation of China (31761163002, 31830031 and 31671053 to YL). The Science and Technology Commission of Shanghai Municipality (18JC1420302) and Innovative Research Team of High-Level Local Universities in Shanghai.
Ethics approval and consent to participate
All animal experimental procedures obey the care and use of laboratory animals of Shanghai Jiao Tong University School of Medicine and approved by the Institutional Animal Care and Use Committee (IACUC).
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional file 1: Table S1.
List of DEGs (Differentially Expressed Genes) from RNA-Seq analysis results. Table S2. List of GO term from RNA-Seq analysis results. Table S3. List of KEGG signing pathway from RNA-Seq analysis results.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Huang, D., Liu, H., Zhu, A. et al. Forebrain excitatory neuron-specific SENP2 knockout mouse displays hyperactivity, impaired learning and memory, and anxiolytic-like behavior. Mol Brain 13, 59 (2020). https://doi.org/10.1186/s13041-020-00591-8
- Conditional knockout mice
- Anxiety-like behavior
- Learning and memory