Shank3 Mice Carrying the Human Q321R Mutation Display Enhanced Self-Grooming, Abnormal Electroencephalogram Patterns, and Suppressed Neuronal Excitability and Seizure Susceptibility

Shank3, a postsynaptic scaffolding protein involved in regulating excitatory synapse assembly and function, has been implicated in several brain disorders, including autism spectrum disorders (ASD), Phelan-McDermid syndrome, schizophrenia, intellectual disability, and mania. Here we generated and characterized a Shank3 knock-in mouse line carrying the Q321R mutation (Shank3Q321R mice) identified in a human individual with ASD that affects the ankyrin repeat region (ARR) domain of the Shank3 protein. Homozygous Shank3Q321R/Q321R mice show a selective decrease in the level of Shank3a, an ARR-containing protein variant, but not other variants. CA1 pyramidal neurons in the Shank3Q321R/Q321R hippocampus show decreased neuronal excitability but normal excitatory and inhibitory synaptic transmission. Behaviorally, Shank3Q321R/Q321R mice show moderately enhanced self-grooming and anxiolytic-like behavior, but normal locomotion, social interaction, and object recognition and contextual fear memory. In addition, these mice show abnormal electroencephalogram (EEG) patterns and decreased susceptibility to induced seizures. These results indicate that the Q321R mutation alters Shank3 protein stability, neuronal excitability, repetitive and anxiety-like behavior, EEG patterns, and seizure susceptibility in mice.

The Shank3 Q321R mutation was identified as a de novo mutation in an individual with ASD who displayed symptoms including social and language deficits, repetitive behaviors (verbal repetitive behaviors, hair pulling, but no motor stereotypies), restricted interests, inattention and irritability (Moessner et al., 2007). This mutation has been shown to decrease excitatory synaptic targeting of Shank3 and Shank3-dependent dendritic spine development, decrease F-actin levels in spines, and suppress excitatory synaptic transmission in cultured hippocampal neurons (Durand et al., 2012). In a more recent study, this mutation was shown to enhance the interaction of Shank3 with Sharpin, but not with α-fodrin (Mameza et al., 2013), two known ligands of the ARR (ankyrin repeat region) domain of Shank3 (Bockers et al., 2001;Lim et al., 2001). In addition, the Q321R mutation has stronger influences on excitatory synapses, as compared with other Shank3 mutations such as R12C and R300C (Durand et al., 2012). These results indicate that the Q321R mutation exerts a significant influence on ASD-related behaviors and excitatory synapse development and function. However, in vivo functions of the Q321R mutation have not been explored.
In the present study, we generated and characterized a new Shank3-mutant mouse line carrying the Q321R mutation (Shank3 Q321R mice) and studied its in vivo effects. We found that this mutation leads to destabilization of Shank3 protein, decreased excitability in hippocampal CA1 pyramidal neurons, enhanced self-grooming and anxiolytic-like behavior, altered electroencephalogram (EEG) patterns, and decreased seizure susceptibility.

Structural Modeling of the Shank3 Protein Containing a Q321R Mutation
The structure of the SPN and ARR domains of the mouse Shank3 protein containing the p.Q321R missense mutation was modeled using the mutagenesis function in PyMOL software (version 1.3) (DeLano, 2009) based on the crystal structure of the SPN and ARR domains of the rat Shank3 protein (PDB ID: 5G4X). Energy minimization and loop flexible modeling were performed using Modeller software (Fiser et al., 2000). Electrostatic charge distribution surfaces were calculated and represented using PyMOL software (version 1.3) (DeLano, 2009). All structural figures were prepared using PyMOL software (version 1.3) (DeLano, 2009).

Stability Prediction of Mutant Shank3 Proteins
The stability of the SPN and ARR domains of Shank3 containing the ASD-risk missense mutations, p.R12C, p.L68P, p.A198G, p.R300C, or p.Q321R, were predicted using the in silico algorithm in I-Mutant 2.0 (version 2.0) 1 under conditions of pH 7.0 and 25 • C (Capriotti et al., 2005). I-Mutant 2.0 is a support vector machine (SVM)-based web server for automatic prediction of stability changes upon a single point mutation. We used the crystal structure of the SPN and ARR domains of Shank3 protein (PDB ID: 5G4X) as a template structure to calculate changes in the stability of the protein containing ASD-risk missense mutations. The G value was calculated by subtracting the unfolding Gibbs free energy of wild-type (WT) protein from that of the mutated protein (kcal/mol). A negative G value indicates a decrease in the stability of the mutated protein.

Animals
A mouse embryonic stem (ES) cell line harboring the Shank3 Q321R mutation was generated by Cyagen. ES cells were injected into C57BL/6N blastocysts to produce chimeric mice, which were crossed with WT C57BL/6J to produce heterozygous knock-in (KI) F1 mice. F1 mice were then crossed with protamine-Cre mice to remove the Neo cassette, followed by backcrossing into a C57BL/6J background for more than five generations. All mice used for experiments in the present study were obtained by mating heterozygotes (HT x HT). The average WT:HT:KI ratio was 0.94:1.94:1.12, consistent with the expected 1:2:1 Mendelian ratio. Breeding was successful in >95% of cases, and homozygous Shank3 Q321R mice showed survival rates and body weights comparable to those of WT mice. The mice were fed ad libitum under 12-h light-dark cycles, and 2-6 mice were grouped in each cage. Mice were bred and maintained according to the Animal Research Requirements of KAIST, and all procedures were approved by the Committee of Animal Research at KAIST (2016-30). For PCR genotyping, the following oligonucleotide primers were used to detect a band of 478 base pairs for WT allele and a band of 622 base pairs for the Q321R mutant allele: forward, 5 -CAT GAG GCA CCC TTT TCT GT-3 , reverse, 5 -TGT CCC TAA CCC CAA TGT GT-3 .

Behavioral Assays
Behavioral analyses were performed during light-off periods on male 2-4-month-old male or female mice in their home cages. WT and mutant littermates (partly pairs) from several mother mice were used to form a cohort. All mice were handled for 10 min per day for 3 days by experienced researchers prior to behavioral experiments. Mouse behavioral assays were performed in the following order: Laboras monitoring of 96-h movements, open-field test, novel-object-recognition test, repetitive behaviors, three-chamber social-interaction test, elevated plus-maze (EPM) test, light-dark test, courtship adult ultrasonic vocalization (USV), and contextual fear conditioning.

Laboras
Long-term locomotion and behaviors of mice were recorded and analyzed using the Laboratory Animal Behavior Observation Registration and Analysis System (LABORAS TM ) by Metris. Mice without habitation to the behavioral booth or Laboras cages were put into Laboras cages where recordings were conducted for 96 consecutive hours. Laboras results were not validated by manual analyses because similar Laboras validations have been reported (Van de Weerd et al., 2001;Quinn et al., 2003;Quinn et al., 2006;Dere et al., 2015).

Open-Field Test
The open-field test was performed as described previously (Gould et al., 2009). A subject mouse was gently introduced into the center zone of the white acryl chamber (40 × 40 × 40 cm) and recorded of movements for 1 h. Locomotor activities and time spent in the center region (20 × 20 cm) of the open field arena were analyzed using EthoVision XT 10 (Noldus).

Novel Object Recognition Test
Novel object recognition memory (Antunes and Biala, 2012) was measured as follows. Briefly, after 1-h habituation in the white acryl chamber 1 day before the experimental day, the subject mouse was allowed to explore two objects of the same shape and material placed symmetrically in the center of the chamber for 10 min. Twenty-four hour later, one object was replaced by a novel object, and the subject mouse was allowed to explore the two objects (one novel and the other familiar) for 10 min. Exploration was defined as events in which the mouse sniffs the object within 1 cm distance or climbs up the object. Time spent in sniffing each object was analyzed by EthoVision XT 10 (Noldus) using the movies from the first 5 min to minimize saturation effects. Genotype comparisons were made using the discrimination index [time spent exploring novel object -time spent exploring familiar object)/(time spent exploring novel + familiar objects) * 100], as previously described (Stefanko et al., 2009;Vogel-Ciernia and Wood, 2014).

Repetitive Behaviors
Autism-relevant repetitive behaviors (Silverman et al., 2010) was measured as follows. The subject mice were allowed to freely move in a clean new home cage with bedding (∼60 lux). Time spent in digging and self-grooming during the 10-min session from a single recording was measured manually in a blind manner. Digging was defined as a mouse uses its head or forelimbs to dig out beddings. Self-grooming was defined as a mouse stroking or scratching its face or body area, or licking its body.

Three-Chamber Social Interaction Test
The three-chamber test, known to measure social approach and social novelty-recognition behavior Nadler et al., 2004;Silverman et al., 2010), was performed as previously described . Briefly, the subject mice were isolated for 4 days before this assay. Three-chambered apparatus (60 × 40 × 20 cm) consisted of the left, center, and right chamber with two entrances to the center chamber. Two empty containers were located in the corner of left and right chamber but not the center chamber. This assay consisted of three sessions. First, the subjects were allowed to freely explore all three chambers for 10 min. Second, the subject mouse was allowed to explore the containers with a stranger 1 (S1) or a novel object (O) for 10 min. Finally, the novel object was replaced with a stranger 2 (S2), and the subject mouse was allowed to explore S1 and S2 for 10 min. During the interval between each session, the subject was gently guided to the center chamber and the two entrances were blocked. Analysis of time spent in sniffing stranger/object was performed using EthoVision XT 10 (Noldus)

Elevated Plus-Maze Test
The EPM test (Pellow et al., 1985) was performed as follows. Briefly, the EPM apparatus, located 75 cm above the floor, consisted of two open arms (30 × 5 × 0.5 cm), two closed arms (30 × 5 × 30 cm), and a center zone with access to both arms. Light conditions were 180 lux for open arms and 20 lux for closed arms. A subject mouse was gently introduced into the center zone and allowed to move freely the open and closed arms for 8 min. Time spent in open/closed arms was measured using EthoVision XT 10 (Noldus).

Light-Dark Test
The light-dark test (Bourin and Hascoet, 2003) was performed as previously described (Ikeda et al., 1995) with minor modifications. Animals were placed in the light chamber with their heads toward the opposite wall from the dark chamber, and allowed to explore the light-dark apparatus (20 × 30 × 20 cm, 670 lux for light chamber, 20 × 13 × 20 cm, 5 lux for dark chamber), which has a 5-cm wide entrance between the two chambers. Time spent in the light chamber was analyzed using EthoVision XT 10 (Noldus).

Direct Social Interaction Test
Direct social interaction test was performed as described previously . Mice, which were isolated in their home cage for 4 days, were used in this assay. On the day of habituation, mice were individually habituated in gray acryl box (33 × 33 × 22 cm) for 10 min. On the day for experiments, two mice in the same genotype that never met each other previously were put into the habituated box and allowed to freely interact each other for 10 min while all behaviors were recorded. Time spent in interaction such as nose-to-nose sniffing, following, nose-to-tail sniffing, and other interactions such as body contacts and huddling were measured manually in a blind manner.

Courtship Adult Ultrasonic Vocalization
Courtship ultrasonic vocalizations (Portfors, 2007) was measured as follows. Male adult subject mice were socially isolated in their home cages for 4 days to allow them to recognize the cage as their own territory. Female adult stranger mice were groupcaged prior to these experiments on the assumption that group housing might synchronize female cycles. A subject male mouse was habituated to a novel test cage for 5 min while recording its basal vocalizations in the absence of a female stranger. Next, a female stranger mouse was introduced into the subject's cage, and the two mice were allowed to interact freely with each other while recording courtship USVs of the subject mouse for 5 min. Because courtship USVs in the context of male-female encounters are mainly produced by males (Maggio and Whitney, 1985;Egnor and Seagraves, 2016), USVs were not categorized as arising from male or female mice. Avisoft SASLab Pro was used to analyze USVs.

Contextual Fear Conditioning and Extinction Test
Contextual fear conditioning and extinction (Phillips and LeDoux, 1992;Milad and Quirk, 2002) were measured as described previously (Lee et al., 2012). The fear conditioning test consisted of two sessions performed on consecutive days. On day 1, subject mice were introduced to the conditioning box and allowed to freely explore the environment for 120 s, and then received a series of 0.8 mA electric foot shocks (1 s; without associated auditory cue) for five times at intervals of 120 s. After the last shock the mice were left in the box for an additional 120 s, making the total experimental time 12 min. On day 2, the mice were placed in the same conditioning box and allowed to freely explore for 10 min without any stimuli, and the freezing levels of the mice on the test day (day 2) were quantified for the first 3 min. These mice, which were tested for 24-h fear memory, were exposed to the same context for 5 min to measure fear extinction every day during the following 7 days (days 3-9). All freezing behaviors were recorded and analyzed using FreezeFrame software (Coulbourn Instruments).

Hot Plate Test
Hot plate test was performed as previously described (Ankier, 1974) with minor modifications. The subject mice were placed in the clear acryl chamber (10 × 10 × 15 cm), which was located on an aluminum heat controller (Sun Electronics Co.) maintained at 50-51 • C. Times taken for the animals to lick their fore-/hindpaws or to jump were measured.

Von Frey Test
The von Frey test was performed according to the standard operating procedures of the Jackson Laboratory Mouse Neurobehavioral Phenotyping Facility (JAX-MNBF). The von Frey test apparatus consisted of a wire floor grid and a clear observation arena with aerated lids (Stoelting Co.). A subject mouse was gently introduced in the observation arena and habituated for 1 h until the activity levels become low. Beginning with the left hind paw, stimuli of von Frey filaments (Stoelting Co.) were presented perpendicularly to the plantar surface of the hind paws of the subject mice for two trials. The initial stimulus started with a von Frey filament with 0.4 g bending force. When there is a withdrawal response, the strength of the stimulus was gradually decreased to lower bending forces (0.16 to 0.02 g), while the stimulus was increased to a higher bending force (0.6 g) when there was no response. The minimum threshold required to induce a withdrawal response for left and right paws was recorded. A withdrawal response is defined as the mouse lifting up their hind paws when a filament is pushed with increasing pressure until it bends. All experimental procedures and analysis were performed in a double-blinded manner.

Electroencephalography
Electroencephalography (EEG) recording was performed as previously described with minor modifications ). An EEG driver was composed of a small header pin connector socket (Hirose Electric, H2021-ND) connected with six stainless steel screws (1 mm × 3 mm) through soldered coated stainless steel wire. All screws were implanted on the skull (two for bilateral frontal, +1.8 mm AP, ±1.0 mm ML; two for bilateral parietal, −2.0 mm AP and ±1.8 mm ML; two for animal ground and reference, −1.0 mm AP and ±1.0 mm ML from lambda). After recovery 1 week after the surgery, EEG recordings were performed using a Cheetah Data Acquisition System (Neuralynx). The subject mice were allowed to freely move around in a white acryl box (25 × 25 × 35 cm) with synchronized video recording. After 20-min habituation in the box, EEGs were recorded for 40 min. EEG data were analyzed using a customized MATLAB code. Baseline EEGs were analyzed using 5-min serial samplings, and the total power spectral density (PSD) was averaged per mouse by applying Fast Fourier Transform (FFT). EEG frequency ranges were defined as follows; Delta, 0-4 Hz; theta, 4-12 Hz; alpha, 12-30 Hz; low gamma, 30-80 Hz; high gamma, 80-130 Hz.

PTZ-Induced Seizure Susceptibility
To determine the susceptibility of mice to PTZ-induced seizures, we performed the experiments using mice that did not receive EEG driver surgery for cleaner results. After intraperitoneal injection of PTZ (Sigma; 50 mg/kg), subject mice were placed in a clean new home cage. Video recordings for 20 min were used to analyze seizure stages defined as follows; stage 1, behavioral arrest; stage 2, myoclonic (jerk) seizures; stage 3: general tonic-clonic seizures, as previously described (Naydenov et al., 2014). The seizure susceptibility score was defined as follows; 0.2 × (latency to stage 1) + 0.3 × (latency to stage 2) + 0.5 × (latency to stage 3).

Statistical Analysis
All behavioral, biochemical, and electrophysiological experiments were performed by researchers blinded to the experimental conditions. Statistical analyses were performed using GraphPad Prism 7 and MATLAB. Data sets, after removal of outliers, were analyzed for significance using unpaired two-tailed Student's t-test, Mann-Whitney U-test, one-sample t-test, log-rank test, and two-way analysis of variance (ANOVA) with post hoc multiple comparisons test when appropriate. The n values and statistical significance values are indicated in figure legends. Statistical significance values are indicated as follows in the figures: * P < 0.05; * * P < 0.01; * * * P < 0.001; ns, not significant. Additional statistical details and results, as well as information on sex, age and numbers of mice used, are presented in Supplementary Table S1.

Molecular Modeling of the Shank3 Q321R Mutation
The Q321R mutation in the SHANK3 gene (Moessner et al., 2007), which strongly influences the synaptic targeting of Shank3 and the development and function of dendritic spines and excitatory synapses (Durand et al., 2012), is located in the ARR domain of the protein. The ARR domain of Shank3 has been shown to interact with the N-terminal SPN domain of Shank3 in an intramolecular manner, leading to suppression of the binding of the ARR domain to Sharpin and α-fodrin (Mameza et al., 2013). These results suggest the possibility that these intra-and intermolecular interactions of the SPN and ARR domains of Shank3 might underlie the effects of the Q321R mutation.
Using the recently reported X-ray crystal structure of the N-terminal SPN and ARR domains of Shank3 (Lilja et al., 2017), we attempted molecular modeling of the SPN and ARR domains harboring the Q321R mutation to determine whether the mutation could interfere with the interaction between these two domains ( Figure 1A). We also included the following additional ASD-risk Shank3 mutations in SPN and ARR domains in our modeling (Durand et al., 2007;Moessner et al., 2007;Gauthier et al., 2010): L68P, which disrupts the interaction between the SPN and ARR domains and enhances the interactions of the ARR with Sharpin (Mameza et al., 2013); R12C, which affects excitatory synapse structure and function (Durand et al., 2012) and disrupts SPN binding to GTP-bound Ras and Rap small GTPases to suppress integrin signaling (Lilja et al., 2017); and R300C, which affects excitatory synapse structure and function (Durand et al., 2012).
We found that the Q321R mutation, unlike the L68P mutation, is located in the C-terminal region of the ARR domain, away from the interface of the SPN and ARR domains, and thus is unlikely to affect the SPN-ARR interaction ( Figure 1A). We next tested whether these mutations affect the stability of SPN and ARR domains by calculating changes in free energy. All five mutations, including the Shank3 Q321R mutation, induced significant decreases in the stability of the protein, with the A198G mutation exerting the strongest effect ( Figure 1B), although the impact of the A198G mutation has not been explored in previous in vitro studies (Durand et al., 2012;Mameza et al., 2013). In addition, an analysis of surface charge distribution indicated that the Q321R mutation induces a substantial increase in the local positive charge density ( Figure 1C). These results suggest that the Q321R mutation is less likely to affect the intramolecular SPN-ARR interaction, but more likely to affect the protein stability or intermolecular interactions of Shank3. (F) The levels of Shank3a, but not Shank3c/d or Shank3e, protein variants, are substantially decreased in whole-brain lysates from WT and Shank3 Q321R/Q321R and Shank3 +/Q321R mice (2 months; female), as determined by immunoblot analysis using Shank3-specific antibodies targeting amino acids 1431-1590 in the proline-rich region of the protein that can detect all three protein splice variants. n = 3 mice (WT, HT, and HM), * P < 0.05, * * P < 0.01, ns, not significant, One sample t-test. (G) Normal gross morphology of the brains of WT and homozygous Shank3 Q321R/Q321R mice (2 months; female), as shown by immunostaining of coronal sections for the neuronal marker, NeuN. Scale bar, 1 mm.
The mutated codon corresponding to the Shank3 Q321R mutation is located in exon 8 of the Shank3 gene, which encodes the C-terminal region of the ARR domain of the Shank3 protein. Therefore, among the multiple variants of the Shank3 protein, Shank3a and Shank3b, containing the ARR domain (Figure 2E), would be most strongly affected compared with ARR-lacking Shank3c/d/e variants (Lim et al., 1999;Maunakea et al., 2010;Wang et al., 2011Wang et al., , 2014Waga et al., 2014;Zhu et al., 2014). Indeed, the levels of Shank3a, the longest (∼240 kDa) variant, which is readily detected by an antibody targeting the proline-rich region (aa 1431-1590), were significantly decreased (by ∼18%) in homozygous Shank3 Q321R/Q321R brains (Figure 2F), in line with the predicted destabilization of the Shank3 protein noted above. Detection of Shank3b protein, which lacks a large portion of the proline-rich region, was not possible because the antibody used targets the proline-rich region. In heterozygous Shank3 +/Q321R brains, Shank3a levels were comparable to those in WT mice, indicative of less severe protein degradation of the mutant Shank3 protein in the heterozygous Shank3-mutant brain.
Shank3 Q321R/Q321R and Shank3 +/Q321R mice were born in the expected Mendelian ratios. In addition, Shank3 Q321R/Q321R mice displayed no detectable abnormalities in the size or gross morphology of the brain, as shown by staining for the neuronal marker, NeuN ( Figure 2G).
When female mice were tested for mEPSCs and mIPSCs, the frequency, but not amplitude, of mIPSCs was increased in CA1 pyramidal neurons (P22-26), whereas mEPSC frequency and amplitude in these neurons (P21-27) were normal (Supplementary Figures S1A,B). Given that Shank proteins are mainly present at excitatory, but not inhibitory, synapses (Boeckers et al., 1999;Naisbitt et al., 1999;Tu et al., 1999;Yao et al., 1999;Bockmann et al., 2002;Heise et al., 2016), these results likely represent a secondary change rather than the direct consequence of Shank3 deletion in pyramidal neurons. We next measured the neuronal excitability of CA1 pyramidal neurons from Shank3 Q321R/Q321R mice (P22-26; male). Neuronal excitability was decreased, as evidenced by the relationship between the amount of injected current and the number of action potentials fired, but the current-voltage relationship, input resistance, and Sag ratio were normal (Figures 3C-F). These results collectively suggest that the Shank3 Q321R mutation suppresses neuronal excitability in male CA1 pyramidal neurons, which would suppress the output function of these neurons.
In addition to the hippocampus, we measured mEPSCs and mIPSCs in the striatum, a brain region strongly implicated in the pathophysiology of Shank3-related autistic-like phenotypes in mice (Peca et al., 2011;Peixoto et al., 2016;Lee et al., 2017). We found that the frequency and amplitude of mEPSCs in dorsolateral striatal neurons in Shank3 Q321R/Q321R mice (P28-43; male) were normal ( Figure 3G). In addition, mIPSC frequency and amplitude were normal in these neurons (P28-43; male) ( Figure 3H). These results suggest that the Shank3 Q321R homozygous mutation does not affect spontaneous excitatory or inhibitory synaptic transmission in the dorsolateral striatum in male mice.

Normal Locomotion and Moderate
Anxiolytic-Like Behavior in Heterozygous Shank3 +/Q321R Mice To determine the behavioral impacts of the Shank3 Q321R mutation, we next subjected Shank3 +/Q321R and Shank3 Q321R/Q321R mice to a battery of behavioral tests. Previous studies on Shank3-mutant mice have shown that Shank3 heterozygosity can lead to ASD-related behavioral abnormalities (Bozdagi et al., 2010;Yang et al., 2012;Duffney et al., 2015;Jaramillo et al., 2016Jaramillo et al., , 2017. We thus first analyzed behavioral abnormalities in heterozygous Shank3 +/Q321R mice. Continuous monitoring of mouse movements for four consecutive days in Laboras cages, representing a familiar environment, revealed that Shank3 +/Q321R mice have normal locomotor activity compared with WT mice (Figures 4A,B). Similarly, Shank3 +/Q321R mice displayed normal locomotor activity in the open-field test, a novel environment (Figures 4C,D). These results suggest that locomotor activities in both familiar and novel environments are not affected by the heterozygous Shank3 Q321R mutation.
Shank3 +/Q321R and WT mice spent comparable amounts of time in the center region of the open-field arena ( Figure 4E). In addition, Shank3 +/Q321R mice spent normal amounts of time in the open arms of the EPM (Figure 4F). In contrast, these mice spent more time in the light chamber of the lightdark apparatus ( Figure 4G). Together, these results suggest that the heterozygous Shank3 Q321R mutation leads to moderate anxiolytic-like behavior but does not affect locomotor activity.

Normal Social Interaction, Moderately
Enhanced Social Communication and Self-Grooming, and Suppressed Digging in Heterozygous Shank3 +/Q321R Mice Shank3 +/Q321R and WT mice showed comparable levels of social approach in the three-chamber test, as shown by time spent sniffing social and object targets ( Figure 5A). Social novelty recognition could not be assessed because WT mice failed to recognize a novel stranger mouse. In tests for USVs in adult male mice encountering a novel female mouse (courtship USVs), Shank3 +/Q321R mice emitted normal numbers of USVs, but the mean duration of each call was increased (Figure 5B), indicative of a moderate and abnormal increase in social communication.
In tests for repetitive behaviors, Shank3 +/Q321R mice displayed enhanced self-grooming and suppressed digging in novel home cages with bedding (10 min; ∼60 lux) (Figures 5C,D), suggesting that the suppressed digging might result from the enhanced selfgrooming. However, long-term (96-h) monitoring of behavior in Laboras cages showed that Shank3 +/Q321R mice exhibit normal self-grooming during the first 10 min, the first 1 or 12 h, the entire session (96 h), and during light-off and light-on periods (48 h each) (Figures 5E-I). These results collectively suggest that the heterozygous Shank3 Q321R mutation does not affect social approach, moderately enhances social communication and self-grooming, and suppresses digging in mice.
Homozygous Shank3 Q321R/Q321R Mice Show Behaviors That Are Largely Similar to Those Observed in Heterozygous Shank3 +/Q321R Mice To test whether there are any dose-dependent effects of the Shank3 Q321R mutation on mouse behaviors, we subjected homozygous Shank3 Q321R/Q321R mice to a battery of behavioral tests that were used for heterozygous Shank3 +/Q321R mice. Shank3 Q321R/Q321R mice showed normal levels of locomotor activity in Laboras cages and in the open-field test (Figures 6A-D). In addition, these mice showed moderately increased anxiolytic-like behaviors, as shown by normal levels of center time in the open-field test but increased open-arm time in the EPM and increased light-box time in the light-dark test (Figures 6E-G), similar to heterozygous Shank3 +/Q321R mice.
Social behavior was normal in homozygous Shank3 Q321R/Q321R mice, as shown by the three-chamber test (Figure 7A), similar to heterozygous Shank3 +/Q321R mice. Again, social novelty recognition could not be assessed because WT mice failed to recognize a novel stranger mouse. Shank3 Q321R/Q321R and WT mice showed comparable levels of social interaction in the direct social interaction test, as shown by the total time spent interacting with an age-and sex-matched male stranger mouse ( Figure 7B).
Homozygous Shank3 Q321R/Q321R mice displayed normal courtship USVs, as shown by the number of USV calls as well as the duration of each USV calls (Figure 7C), partly similar to heterozygous Shank3 +/Q321R mice that displayed a normal number of courtship USV calls but increased duration of each USV calls. Lastly, homozygous Shank3 Q321R/Q321R mice displayed increased self-grooming and decreased digging in novel home cages with bedding (10 min) (Figures 7D,E), but normal Laboras-cage self-grooming (Figures 7F-J), largely similar to heterozygous Shank3 +/Q321R mice.
Because the differences between novel home-cage and Laboras-cage environments, where both Shank3 +/Q321R and Shank3 Q321R/Q321R mice showed positive and negative repetitive self-grooming, respectively, include the novelty of the space (less novel in new home cages and more novel in Laboras cages) and the light intensity (bright light in home cages and complete darkness in Laboras cages). To differentiate these factors, we first tested a new condition, novel homecage environment without light, and, intriguingly, could not observe enhanced self-grooming in Shank3 Q321R/Q321R mice ( Supplementary Figure 2A), suggesting that the presence of light is important. However, the presence of light in Laboras cages did not enhance self-grooming in Shank3 Q321R/Q321R mice during the first 10 min (Supplementary Figure 2B), similar to the results from Laboras cages without light. Furthermore, the presence of light in a novel open-field arena did not enhance self-grooming in Shank3 Q321R/Q321R mice (Supplementary Figure 2C). These results indicate that self-grooming in Shank3 Q321R/Q321R mice is enhanced selectively in novel home cages in the presence of light. In addition, these results suggest that novelty of space/environment per se is not important, but, rather, an increase in stress associated with the light in novel home cages (but not in other environments) strongly enhances self-grooming in Shank3 Q321R/Q321R mice.
Together, these results suggest that homozygous Shank3 Q321R/Q321R mice show behavioral abnormalities that are largely similar to those observed in heterozygous Shank3 +/Q321R mice such as moderately enhanced anxiolytic-like behavior and self-grooming, although there were minor differences in a sub-parameter of USVs (duration of each USV calls). In addition, the largely similar behaviors of Shank3 +/Q321R and Shank3 Q321R/Q321R mice suggest that there is no strong dose-dependent effect of the Shank3 Q321R mutation on mouse behaviors.

Normal Novel Object-Recognition and Contextual Fear Memory in Homozygous Shank3 Q321R/Q321R Mice
In the novel object-recognition test, in which a mouse is exposed to two identical objects on day 1 and then to one original object and a new object on day 2, Shank3 Q321R/Q321R and WT mice showed similar discrimination index scores for the novel object ( Figure 8A), indicative of normal object-recognition memory in the mutant mice.
In the contextual fear-conditioning test, Shank3 Q321R/Q321R mice showed normal acquisition of fear memory on the training day (day 1) (Figure 8B). Twenty-four hours later (day 2), Shank3 Q321R/Q321R and WT mice showed comparable levels of freezing in the same context, suggestive of normal retrieval of contextual fear memory. When these mice were subsequently exposed to the same context every day for seven consecutive FIGURE 5 | Normal social interaction, moderately enhanced social communication and self-grooming, and suppressed digging in heterozygous Shank3 +/Q321R mice. (A) Normal social approach in Shank3 +/Q321R mice (2-3 months; male) in the three-chamber test, as shown by time spent sniffing. S1, a stranger; O, object; S2, new stranger. Social novelty recognition, measured by the preference for new stranger (S2) over old stranger (S1), could not be determined due the lack of normal social novelty recognition in WT mice. n = 10 mice (WT), 13 mice (HT), * * * P < 0.001, ns, not significant, Welch's t-test, Mann-Whitney U-test, and Student's t-test. (B) Moderately increased courtship USVs emitted by Shank3 +/Q321R mice (2-3 months; male) upon encounter with a novel female stranger, as shown by the normal number of USVs but increased mean duration of each USV calls. n = 12 mice (WT), 12 mice (HT), * P < 0.05, ns, not significant, Student's t-test. (C) Enhanced self-grooming in Shank3 +/Q321R mice (2-3 months; male) in home cages with bedding (10 min), as shown by total self-grooming time. n = 12 mice (WT), 13 mice (HT), * P < 0.05, Student's t-test. (D) Suppressed digging in Shank3 +/Q321R mice (2-3 months; male) in home cages with bedding (10 min), as shown by total digging time. n = 12 mice (WT), 13 mice (HT), * P < 0.05, Welch's t-test. (E-I) Normal self-grooming in Shank3 +/Q321R mice (2-3 months; male) in Laboras cages, as shown by total self-grooming duration. n = 13 mice (WT), 13 mice (HT), ns, not significant, two-way ANOVA (genotype main effect p-value = 0.4087 in E), Mann-Whitney test (F,H,I), Student's t-test (G). FIGURE 6 | Normal locomotor activity and moderate anxiolytic-like behavior in homozygous Shank3 Q321R/Q321R mice. (A,B) Normal locomotor activity in Shank3 Q321R/Q321R mice (2-4 months; male), in Laboras cages, where locomotor activity was measured together with other movements for four consecutive days in the absence of habituation. Note that the WT data in this panel and other panels in this figure are identical to those shown in Figure 4 because WT, heterozygous, and homozygous mice were tested together in the same behavioral tests. OFF/ON, light-off/on. n = 13 mice (WT) and 14 mice (Q321R), ns, not significant, two-way ANOVA (A, genotype main effect p-value = 0.9356) and Student's t-test (B). (C,D) Normal locomotor activity in Shank3 Q321R/Q321R mice (2 months; male) in the open-field test, as shown by the distance moved. Note that Shank3 Q321R mice spent a normal amount of time in the center region of the open-field arena, suggestive of normal anxiety-like behavior. n = 13 mice (WT) and 14 mice (Q321R), * * * P < 0.001, ns, not significant, repeated measures two-way ANOVA (C) and Student's t-test days (day 3-9) for fear extinction, they displayed comparable decreases in freezing levels, suggestive of normal fear extinction in Shank3 Q321R/Q321R mice. In addition, Shank3 Q321R/Q321R mice showed normal levels of somatosensory functions, as determined by the hot plate and von Frey tests (Figures 8C,D). These results collectively suggest that the homozygous Shank3 Q321R mutation does not affect objectrecognition memory, acquisition, or extinction of contextual fear memory, or somatosensory functions.

Female Homozygous Shank3 Q321R/Q321R Mice Show Normal Anxiety-Like Behavior and Reduced Self-Grooming
To test whether there is a male-female difference in the impacts of the Shank3 Q321R mutation on behaviors, we subjected female Shank3 Q321R/Q321R mice to the behavioral tests in which male Shank3 Q321R/Q321R mice showed abnormal behaviors. To our surprise, female Shank3 Q321R/Q321R mice FIGURE 7 | Continued Frontiers in Molecular Neuroscience | www.frontiersin.org FIGURE 7 | Normal social interaction and communication, moderately enhanced stress-induced self-grooming, and suppressed digging in homozygous Shank3 Q321R/Q321R mice. (A) Normal social approach in Shank3 Q321R/Q321R mice (3 months; male) in the three-chamber test, as shown by time spent sniffing social (S1) or object (O) target. Social novelty recognition, measured by the preference for new stranger (S2) over old stranger (S1), could not be determined due the lack of normal social novelty recognition in WT mice. Note that the WT data in this panel and other panels in this figure are identical to those shown in Figure 5 because WT, heterozygous, and homozygous mice were tested together in the same behavioral tests. n = 10 mice (WT), 14 mice (Q321R), * * * P < 0.001, ns, not significant, Welch's t-test (for S1 vs. O) and Student's t-test (for S1 vs. S2). (B) Normal social interaction in Shank3 Q321R/Q321R mice (3 months; male) in bidirectional direct social-interaction tests, as shown by total time spent in social interaction. n = 10 pairs (WT), 13 pairs (Q321R), ns, not significant, Student's t-test. (C) Normal USVs emitted by Shank3 Q321R/Q321R mice (4 months; male) upon encounter with a novel female stranger (courtship USVs), as shown by the total number of USVs and the mean duration of each USV calls. n = 12 mice (WT), 14 mice (Q321R), ns, not significant, Student's t-test. (D) Enhanced self-grooming in Shank3 Q321R/Q321R mice (2 months; male) in home cages with bedding (10 min), as shown by total self-grooming time. n = 12 mice (WT), 13 mice (Q321R), * * * P < 0.001, Student's t-test. (E) Suppressed digging in Shank3 Q321R/Q321R mice (2 months; male) in home cages with bedding (10 min), as shown by total digging time. n = 12 mice (WT), 13 mice (Q321R), * P < 0.05, Welch's t-test. (F-J) Normal self-grooming in Shank3 Q321R/Q321R mice (2 months; male) in Laboras cages, where self-grooming was measured together with other movements for four consecutive days in the absence of habituation. OFF/ON, light-off/on. Note that there are no genotype differences during the first 10 min, 1 or 12 h, the entire session (96 h), or during light-off and light-on periods (48 h each). n = 13 mice (WT) and 14 mice (Q321R), ns, not significant, two-way ANOVA (genotype main effect p-value = 0.9754 in F), Student's t-test (H) and Mann-Whitney U-test (G,I,J).
showed normal levels of anxiety-like behaviors in both EPM and light-dark tests (Supplementary Figures 1C,D), dissimilar to male Shank3 Q321R/Q321R mice that showed anxiolytic-like behavior in both tests (Figures 6F,G).
Furthermore, female Shank3 Q321R/Q321R mice showed normal levels of self-grooming and digging in novel home cages with bedding (10 min) (Supplementary Figures 1E,F), again dissimilar to male Shank3 Q321R/Q321R mice that showed enhanced self-grooming and reduced digging in novel home cages (Figures 7D,E). These results collectively suggest that female Shank3 Q321R/Q321R mice do not show anxiolytic-like behavior or enhanced self-grooming, dissimilar to male Shank3 Q321R/Q321R mice.

Abnormal EEG Patterns and Decreased
Susceptibility to Induced Seizures in Homozygous Male Shank3 Q321R/Q321R Mice Abnormal electroencephalogram (EEG) patterns have been observed in ASD (Wang et al., 2013), ASD and PMS associated with SHANK3 mutations (Moessner et al., 2007;Soorya et al., 2013;Figura et al., 2014;Holder and Quach, 2016), as well as in Shank3-mutant mouse models (Han et al., 2013;Dhamne et al., 2017). In particular, the individual with ASD carrying the Q321R mutation was reported to display abnormal EEG, bilateral epileptiform discharges without seizures, and severe sleep disorders (Moessner et al., 2007). We thus attempted bilateral EEG recordings in the frontal and parietal lobes of the male Shank3 Q321R/Q321R brain using EEG drivers implanted on the skull. The power of EEG in the delta range, but not other frequency ranges (0-4 Hz), was decreased in the frontal lobes, whereas EEG power in the alpha frequency range (12-30 Hz) was increased in the parietal lobes (Figures 9A,B).
Abnormalities in the left-right hemispheric asymmetry of EEGs have also been observed in ASD (Wang et al., 2013). We thus analyzed the changes in EEGs in left and right hemispheres of the Shank3 Q321R/Q321R brain separately. In the frontal lobe, left and right hemispheres showed similar decreases in normalized EEG power in the delta range, but not other ranges (Supplementary Figure 3A). Intriguingly, in the parietal lobe, the right but not left hemisphere showed decreased EEG power in the delta range and increased EEG power in theta and alpha ranges (Supplementary Figure 3B). These results suggest that the decreased delta-band EEG power in the frontal lobe likely involves both hemispheres, whereas the increase alpha-band EEG power in the parietal lobe likely involves the right parietal lobe.
These abnormal EEG patterns in the Shank3 Q321R/Q321R brain and reduced neuronal excitability in Shank3 Q321R/Q321R hippocampal CA1 neurons (Figure 3) suggest the possibility that the Shank3 Q321R/Q321R mutation disrupts the balance between neuronal excitation and inhibition in the brain, a mechanism suggested to underlie ASD (Rubenstein and Merzenich, 2003). To test this, we assessed the susceptibility of Shank3 Q321R/Q321R mice to seizures induced by the GABA A receptor antagonist pentylenetetrazole (PTZ).
We found that male Shank3 Q321R/Q321R mice are less susceptible to PTZ-induced seizures compared with WT mice, as shown by the latency to seizure stage 1 (although not stage 2), latency-based seizure susceptibility index, and final seizure stages reached (Figure 10). Spontaneous behavioral seizures were not observed in the absence of PTZ injection. These results collectively suggest that the homozygous Shank3 Q321R mutation induces abnormal EEG patterns and decreased PTZ-induced seizure susceptibility in male mice.

DISCUSSION
Our results indicate that the Shank3 Q321R mutation, located in exon 8 of the Shank3 gene encoding the C-terminal part of the ARR domain of the protein, leads to substantial destabilization of the ARR-containing Shank3a protein variant and synaptic, behavioral, and EEG/seizure abnormalities in mice.
Our in vivo data from Shank3 Q321R/Q321R mice indicate that the Q321R mutation induces a substantial decrease (∼18% of WT levels) in the total level of the ARR-containing Shank3a protein variant in the Shank3 Q321R/Q321R brain ( Figure 2F). In addition, our molecular modeling analyses suggest that the Q321R mutation is less likely to affect the intramolecular SPN-ARR interaction but more likely to affect the protein stability or intermolecular interactions of Shank3. The ARR domain of Shank3 is known to interact with Sharpin and α-fodrin (Bockers et al., 2001;Lim et al., 2001), and the Q321R mutation has FIGURE 8 | Normal object-recognition and fear memory in homozygous male Shank3 Q321R/Q321R mice. (A) Normal levels of object recognition memory in Shank3 Q321R/Q321R mice (2-4 months; male) in the novel object recognition test, as shown by the discrimination index over a familiar object and a novel object (see Materials and Methods for details) presented 24 h after exploring two identical objects on the first day. n = 11 mice (WT), 13 mice (Q321R), ns, not significant, Student's t-test. (B) Normal acquisition, retrieval, and extinction of contextual fear memory Shank3 Q321R/Q321R mice (3-6 months; male), as shown by freezing levels. Mice were given 5 foot shocks (2-min intervals) during the 12-min fear memory acquisition phase (day 1), and were consecutively exposed to the same context 24 h after the training (day 2) and also during days 3-9 for fear extinction. n = 13 mice (WT), 17 mice (Q321R), repeated measures two-way ANOVA. (C) Normal somatosensory function in Shank3 Q321R/Q321R mice (6 months; male) in the hot plate test, as shown by latency to first licking/jumping. n = 17 mice (WT) and 20 mice (Q321R), ns, not significant, Student's t-test (for latency to first licking), and Mann-Whitney U-test (for latency to first jumping). (D) Normal somatosensory function in Shank3 Q321R/Q321R mice (2 months; male) in the von Frey test, as shown by threshold for response to stimulation. n = 10 mice (WT) and 11 mice (Q321R), ns, not significant, Mann-Whitney U-test. been shown to enhance Shank3 binding to Sharpin but not α-fodrin (Mameza et al., 2013). Given that Shank proteins can be ubiquitinated and deubiquitinated in an activity-dependent manner to regulate excitatory synapse structure and function (Ehlers, 2003;Kerrisk Campbell and Sheng, 2018) and that Sharpin is a component of the E3 ligase complex termed LUBAC (linear ubiquitin chain assembly complex) (Ikeda et al., 2011;Tokunaga et al., 2011;Iwai et al., 2014;Hrdinka and Gyrd-Hansen, 2017), it is tempting to speculate that the strong decrease in the levels of the mutant Shank3 protein (Q321R) might involve enhanced protein ubiquitination. Dissimilar to the results from Shank3 Q321R/Q321R mice, however, Shank3 +/Q321R mice show normal levels of Shank3 proteins, as compared with WT mice (Figure 2F). Therefore, distinct pathophysiological mechanisms of the Shank3 Q321R mutation at the protein level (i.e., decreased protein levels vs. abnormal protein function) might lead to similar behavioral abnormalities in Shank3 Q321R/Q321R and Shank3 +/Q321R mice, likely by acting on similar sets of neurons in the brain.
in Table 1). These results suggest that homozygous deletion of ARR-containing Shank3 protein variants does not affect mEPSCs in hippocampal or striatal neurons. Potential reasons for these results could be that ARR-containing Shank3 protein variants only minimally affect excitatory synapse development and function in hippocampal and striatal neurons or that some compensatory changes have occurred. Similar to our results, mEPSC frequency and amplitude are unchanged in Shank3 mouse lines carrying similar ARR deletions, including Shank3 9 mice (lacking exon 9)  and Shank3 4−9 mice (lacking exons 4-9) Jaramillo et al., 2016). Notably, however, heterozygous deletion of Shank3 exons 4-9 in mice induces an increase in mEPSC frequency and a decrease in mEPSC amplitude (Bozdagi et al., 2010), suggesting different effects of homozygous and heterozygous deletions.
Whereas mEPSCs and mIPSCs were unchanged in male Shank3 Q321R/Q321R CA1 pyramidal neurons, and mEPSCs were unchanged in female Shank3 Q321R/Q321R CA1 pyramidal neurons, the frequency but not amplitude of mIPSCs was increased in female CA1 neurons (Figure 3 and Supplementary  Figures 1A,B). These changes may involve some presynaptic changes given that Shank3 proteins are mainly localized at excitatory postsynaptic sites. Alternatively, the Shank3 protein with Q321R mutation might be localized at subcellular sites other than excitatory postsynaptic sites or even at inhibitory postsynaptic sites in CA1 pyramidal neurons, although our Shank3 antibodies were not good enough for immunohistochemical analyses at the electron microscopic level. Intriguingly, a similar increase in mIPSC frequency has been reported in Shank3 9 CA1 pyramidal neurons , although mIPSCs were not affected in Shank3 4−9 CA1 pyramidal neurons .
(G) Final seizure stages that animals reached, expressed as a percentage of animals. * P < 0.05, Chi-square test.
Shank3 mice carrying heterozygous and homozygous deletion of exons 13-16 encoding the PDZ domain of Shank3 (Peca et al., 2011) display markedly increased neuronal excitability in association with increased input resistance and decreased hyperpolarization-activated currents (I h ). Moreover, human neurons with a SHANK3 exon 13 deletion display similar changes (Yi et al., 2016). This difference may be attributable to differences in the specific exons of the Shank3 gene affected: exon 8 in our study of the Q321R mutation and exons 13-16 in these previous studies. The former encompasses mainly ARR-containing Shank3 protein variants, whereas the latter encompasses almost all Shank3 protein variants. Although further details remain to be elucidated, our electrophysiology results suggest that, at minimum, the decreased neuronal excitability and increased mIPSC frequency in Shank3 Q321R/Q321R CA1 pyramidal neurons would substantially suppress the output functions of these neurons. Behaviorally, heterozygous Shank3 +/Q321R mice show moderately enhanced social communication and self-grooming and moderate anxiolytic-like behaviors in the light-dark tests (although open-field center activity and EPM performance were normal), while showing normal levels of locomotor activity and three-chamber social interaction (Figures 4, 5). Homozygous Shank3 Q321R/Q321R mice show behaviors that are largely similar to those in heterozygous Shank3 +/Q321R mice, displaying moderately enhanced self-grooming, moderate anxiolytic-like behaviors in EPM and light-dark tests (although open-field center activity was normal), while showing normal levels of locomotor activity, social interaction (three-chamber and direct social interaction), social communication (courtship USV), and object-recognition and contextual fear memory (Figures 6-8; summarized in Table 1).
Notably, Shank3 9 mice, which carry a similar alteration in the ARR domain, show behaviors that are largely similar to those of Shank3 Q321R/Q321R mice, displaying normal levels of locomotor activity, anxiety-like behavior (open-field center time), social interaction and communication (three-chamber and pup USV), and object-recognition memory, although self-grooming was also normal in these mice . Behavioral impacts in other ARR-deletion-carrying mouse lines are more severe compared with those reported here, a difference that is likely attributable to larger deletions in the ARR domain of the Shank3 protein. Specifically, self-grooming is enhanced in Shank3 4−9 mice, similar to our mice, but these former mice additionally show altered social interaction and communication and impaired novel object-recognition memory (Bozdagi et al., 2010;Wang et al., 2011;Yang et al., 2012;Jaramillo et al., 2016;Dhamne et al., 2017). In addition, Shank3 4−7 mice show abnormal social novelty recognition, but normal self-grooming (Peca et al., 2011), again different from our mice.
Identifying synaptic mechanisms and neural circuits associated with the behavioral phenotypes in Shank3 Q321R mice will require additional investigation. Notably, however, while male Shank3 Q321R/Q321R mice show normal mEPSCs and mIPSCs in the hippocampus and anxiolytic-like behavior, enhanced self-grooming, and suppressed digging, female Shank3 Q321R/Q321R mice show increased mIPSC frequency in the hippocampus and normal anxiety-like behavior, self-grooming, and digging (Figures 3, 6, 7 and Supplementary Figures 1C-F). Therefore, the increased mIPSC frequency in female mice might represent a part of the compensatory changes that might have normalized the anxiety-like behaviors and repetitive behaviors in female Shank3 Q321R/Q321R mice. In addition, the most strongly affected Shank3 protein variant in our mice was Shank3a, which is abundantly expressed in the striatum among many brain regions . In addition, the striatopallidal pathway has been shown to regulate self-grooming in Shank3 13−16 mice (Wang et al., 2017). In addition, the striatum has emerged as one of the key brain regions involved in the regulation of anxiety-like behaviors, in addition to the amygdala, bed nucleus of the stria terminalis, hippocampus, and ventromedial prefrontal cortex (Lago et al., 2017). Therefore, it is tempting to speculate that striatal dysfunctions might contribute to the self-grooming and anxiolytic-like behavior in Shank3 Q321R mice. Although male Shank3 Q321R/Q321R mice displayed normal levels of excitatory and inhibitory spontaneous synaptic transmission (mEPSCs and mIPSCs) in the dorsolateral striatum (Figures 3G,H), additional investigations might reveal some other striatal dysfunctions in synaptic transmission and plasticity. For instance, previous studies have reported decreased NMDA/AMPA ratio but normal mEPSCs and mIPSCs in the dorsal striatum of homozygous and heterozygous Shank3 4−9 mice (Jaramillo et al., 2016).
Shank3 Q321R/Q321R mice (males) display abnormal baseline EEG patterns, with a decreased delta band in the frontal lobe and increased alpha band in the parietal lobe (Figure 9). In addition, we found left-right asymmetry in the parietal, but not frontal, lobe of the Shank3 Q321R brain. Understanding the biological significances of these results would have to involve, for instance, additional analyses of sleep behaviors and rhythms in Shank3 Q321R/Q321R mice and acquisition of EEG-and sleeprelated details from the individual carrying the Q321R mutation (Moessner et al., 2007). Nonetheless, our data are in line with the previous studies in ASD that have reported abnormalities in the power, coherence, and asymmetry of EEGs in the delta and alpha ranges in addition to other frequency ranges (Wang et al., 2013). Moreover, the individual carrying the SHANK3 Q321R mutation was shown to display abnormal EEG patterns (Moessner et al., 2007), although it is unclear whether the abnormal EEGs were observed in similar frequency ranges. Notably, a previous study on Shank3 13−16 mice reported abnormally enhanced EEGs in the low gamma range (30-80 Hz), known to be associated with parvalbumin-positive GABAergic neurons (Sohal et al., 2009), but normal EEGs in other frequency ranges (Dhamne et al., 2017). This difference could again reflect differences in the exons deleted in the two mutant mouse lines, with Shank3 13−16 mice lacking a larger number of Shank3 splice variants.
Our Shank3 Q321R/Q321R mice (males) are more resistant to PTZ-induced seizures and do not display spontaneous behavioral seizures (Figure 10). These results are in line with the suppressed excitability of CA1 pyramidal neurons in these animals. Similar to our results, a previous study on Shank3 13−16 mice revealed a decrease in susceptibility to PTZ-induced seizures (Dhamne et al., 2017) and lack of spontaneous behavioral seizures (Peca et al., 2011). In addition, Shank3 4−22 mice show tonic hyperactivity in the corticostriatal-thalamic axis, revealed by multi-site in vivo recordings, but do not show spontaneous behavioral seizures . In contrast, a transgenic mouse line carrying duplicated Shank3 shows enhanced epileptiform spikes in the dentate gyrus and electrographic seizures (Han et al., 2013). These results collectively suggest that Shank3 is an important regulator of excitatory drive in the brain. Perhaps more importantly, the individual with the Q321R mutation was shown to display bilateral epileptiform discharges without seizures, in addition to abnormal EEG patterns (Moessner et al., 2007). Therefore, the EEG and seizure phenotypes of our Shank3 Q321R mice may serve as potential biomarkers for future studies.

CONCLUSION
Our data suggest that the Shank3 Q321R mutation in mice has significant influences on Shank3 protein stability, hippocampal neuronal excitability, anxiety-like and repetitive behaviors, EEG, and seizure susceptibility.

ETHICS STATEMENT
The mice were bred and maintained according to the Animal Research Requirements of KAIST, and all procedures were approved by the Committee of Animal Research at KAIST (2016-30).

AUTHOR CONTRIBUTIONS
JL generated the mice. Y-EY and TY performed the electrophysiological and behavioral experiments. Y-EY performed the immunoblot experiments. SL performed the immunohistochemistry and behavioral experiments. Y-EY, SL, and TY performed the EEG experiments. DK performed the structural modeling. H-MH, Y-CB, and EK wrote the manuscript.

FUNDING
This study was supported by the National Research Foundation of Korea (NRF-2017R1A5A2015391 to Y-CB) and the Institute for Basic Science (IBS-R002-D1 to EK).

FIGURE S3 | Abnormal EEG patterns in the left and right hemispheres in
Shank3 Q321R/Q321R mice. (A,B) Altered left and right hemispheric EEG powers in the frontal and parietal lobes of the Shank3 Q321R/Q321R brain (3 months; male). Note that the frontal lobe displays decreased delta power in both left and right hemispheres, whereas the parietal lobe displays decreased delta power and increased theta and alpha powers only in the right hemisphere. n = 10 mice (WT) and 10 mice (Q321R) for frontal lobe and 9 mice (WT) and 10 mice (Q321R) for parietal lobe, * P < 0.05, * * P < 0.01, ns, not significant, Student's t-test, Mann-Whitney U-test, Welch's t-test.