Identification of a novel Shank2 transcriptional variant in Shank2 knockout mouse model of autism spectrum disorder
Molecular Brain volume 13, Article number: 54 (2020)
Autism spectrum disorder (ASD) is a group of neurodevelopmental disorders that are highly heterogeneous in clinical symptoms as well as etiologies. Mutations in SHANK2 are associated with ASD and accordingly, Shank2 knockout mouse shows ASD-like behavioral phenotypes, including social deficits. Intriguingly, two lines of Shank2 knockout (KO) mouse generated by deleting different exons (exon 6–7 or exon 7) showed distinct cellular phenotypes. Previously, we compared gene expressions between Shank2 KOs lacking exon 6–7 (e6–7 KO) and KOs lacking exon 7 (e7 KO) by performing RNA-seq. In this study, we expanded transcriptomic analyses to identify novel transcriptional variants in the KO mice. We found prominent expression of a novel exon (exon 4′ or e4’) between the existing exons 4 and 5 in the Shank2 e6–7 KO model. Expression of the transcriptional variant harboring this novel exon was confirmed by RT-PCR and western blotting. These findings suggest that the novel variant may function as a modifier gene, which contributes to the differences between the two Shank2 mutant lines. Furthermore, our result further represents an example of genetic compensation that may lead to phenotypic heterogeneity among ASD patients with mutations in the same gene.
Autism spectrum disorder (ASD) is a group of neurodevelopmental disorders that are diagnosed by impairments in social–communication and restricted, repetitive behaviors or interests . However, ASD is highly heterogeneous in terms of clinical symptoms and etiology. Although ASDs are highly heritable, and recent advances in next generation sequencing have revealed more than 1000 mutations associated with ASD, most of those identified alleles are rare variants that represent only a small fraction of ASD . Moreover, it is becoming clear that ASD involves multigenic interactions, as well as interactions between risk genes and various environmental factors, which further highlight its complexity [2, 3].
Mutations in SHANK genes are associated with ASD and other psychiatric conditions such as intellectual disabilities, schizophrenia, and Phelan-Mcdermid Syndrome [4, 5]. SHANK genes (SHANK1, 2, 3) encode the SH3 and multiple ankyrin repeat domain proteins, which are post-synaptic scaffold proteins of excitatory synapses [6, 7]. Among SHANK genes, over 80 and 30 ASD-associated mutations have been identified in SHANK3 and SHANK2, respectively (https://gene.sfari.org/database/human-gene/) [8, 9]. While Shank3 mutant mice have been most extensively characterized, previous studies demonstrated that Shank2 null mutant mice also exhibit ASD-like behavioral deficits and synaptic dysfunctions [7, 10,11,12]. Interestingly, Shank2 KO mice lacking exons 6 and 7 (Shank2 e6–7 KO) showed hypo-function of N-methyl-D-aspartate receptor (NMDAR) and impaired long-term potentiation (LTP), whereas the other Shank2 mutant lacking exon 7 of the Shank2 gene (e7 KO) showed hyper-function of NMDAR and enhanced LTP in the adult hippocampus [10, 11]. Recent studies showed that differences in gene expressions, experimental conditions, genetic background, and the developmental stages of the mice may, at least partially, explain the opposite cellular phenotypes in the two Shank2 mutant mouse models [13,14,15]. We have previously compared the gene expression profiles of Shank2 e6–7 and e7 KO mice by performing RNA-seq analyses, and found that the two KO mouse lines share only a few differentially expressed genes . In this study, we analyzed RNA-seq data from the hippocampi of two Shank2 mutant mouse models (e7 and e6–7 KO), and found a novel transcript that encodes a new exon in e6–7 and e7 KO mice. We confirmed the expression of this novel transcript at both mRNA and protein levels, which suggested that this novel variant may act as a modifier gene that contributes to differences between the two Shank2 mouse mutant lines.
Materials and methods
The RNA-Seq raw reads [GSE79824 (GSM2104342 - GSM2104353) and GSE47966] were processed and mapped against the mouse genome, as previously described . The uniquely mapped reads were used to reconstruct Shank2 gene transcripts (guided transcriptome reconstruction), as previously described . The transcribed regions based on the mapped RNA-Seq reads were analyzed to detect possible novel exons in the known gene. The spliced junction reads between a known exon and a novel exon were used to confirm the expression of the novel exon and to reconstruct the transcript isoforms. Conserved regions on the novel exon were evaluated using the comparative genomics tracks in the UCSC genome browser. Multiple sequence alignments among eight different species (human, rhesus, mouse, rat, dog, cow, chicken and frog) was performed using the MUSCLE algorithm implemented in EMBL-EBI bioinformatics tools . Based on the multiple sequence alignment result, a phylogenetic tree was constructed using the neighbor-joining method implemented in MEGA . Methylation information in the mouse genome was obtained from the “DNA Methylation” track hub in the UCSC genome browser provided by MethBase .
Quantitative real time polymerase chain reaction (PCR)
Total RNA samples were extracted from the prefrontal cortex, hippocampus, and cerebellar hemisphere region of 4-week and 8-week old male C57BL/6 N mice (Orient Bio). cDNA was extracted from total RNA by reverse-transcription using SuperScript III (Invitrogen) for RNA-seq. Primer sequences used were as follows: forward 5′-CGGTCTGAGTGAGATGGTCA-3′ (exon 4′) and reverse 5′-TAGGAGCCCACCGTGTAATG-3′ (exon 5). Quantitative RT-PCR was performed with SYBR Premix Ex Taq II (Tli RNase H Plus) (Takara) on the ABI7300 instrument. Initial activation was conducted at 95 °C for 30 s. This was followed by 40 cycles of 95 °C for 5 s and 60 °C for 30 s. Single peaks from dissociation cycles were confirmed, which were indicative of a single PCR product. The 2-ΔΔCt method was used for quantitative comparison, using GAPDH as the control. RT-PCR was performed using Kapa HiFi Hotstart Readymix (Kapa Biosystems). Initial activation at 94 °C was carried out for 30 s, and was followed by 34 cycles of 94 °C for 30 s and 58 °C for 30 s (Actb and ex4’f-ex5r), or 60 °C for 30 s (ex2f-ex4r and ex4f-ex6r) and 72 °C for 60 s. The RT-PCR products were visualized in agarose gels.
Western blot analyses were performed as previously described [14, 20]. Membranes were incubated with a rabbit anti-SHANK2 (ab171189, Abcam) or anti-tubulin (T4026, Sigma) antibodies overnight at 4 °C, followed by incubation with an appropriate HRP-conjugated secondary antibody for 2 h at room temperature. Chemiluminescence signals were acquired and analyzed using ChemiDoc (Bio-Rad).
Identification and validation of a novel transcript
We extended our comparative transcriptome analysis by measuring expression densities through the transcripts, including intronic regions. Previous studies have identified multiple splicing variants of Shank2, including short Shank2a and long Shank2b [21, 22]. In this study, we also identified a previously unknown intronic region via RNA sequencing mapping, which suggests that there may be an additional Shank2 splicing variant in the mouse hippocampus (Fig. 1a). To characterize the region in more detail and determine whether the identified expression region was a true transcribed region, we reconstructed KO and wild type (WT) Shank2 transcripts based on the RNA raw sequencing reads (Fig. 1a). The reads from the spliced junction between the known exons and newly identified transcribed regions were used to confirm the expression of the novel exon. We were able to reconstruct Shank2 transcripts in each sample (e6–7 KO, e6–7 WT littermate, e7 KO, and e7 WT littermate), and confirmed the previously known transcripts (Fig. 1a). We identified the transcribed regions, chr7:144238708–144,239,093 (e6–7 KO) and chr7:144238699–144,239,093 (e7 KO) (red box in Fig. 1a), which were located in exons 4 and 5 of the Shank2 transcript in the reconstructed e6–7 and e7 KO transcripts. This novel transcribed region has not been previously annotated in any gene annotation database, and we annotated the novel transcribed region as exon 4′. We found that the novel exon 4′ is connected with exon 5, but not with the upstream exons; therefore, exon 4′ is the first exon of this novel transcript. More interestingly, although expression of the novel transcript containing exon 4′ was detected in both e6–7 and e7 KO, but very weakly in e7 KO (Fig. 1a).
To verify the expression of the novel Shank2 transcript containing the newly discovered exon 4′, we designed PCR primers with the forward primer in exon 4′ (ex4’f) and the reverse primer in exon 5 (ex5r) (Table 1). Quantitative RT-PCR detected a single PCR product, confirming the existence of the exon 4′-5 containing transcript. Expression of the novel transcript was markedly increased by approximately 8 folds in the e6–7 KO compared with that of e6–7 WT littermate (Fig. 1b). However, consistent with the RNA-Seq results, expression of the novel transcript was considerably less increased in e7 KO (~ 1.71 fold compared with that of ex7 WT littermate, Fig. 1b). To further confirm that the novel transcript is translated into a protein of the expected size in the brains of e6–7 KO mice (Fig. 1c), we performed western blot analyses. A clear protein band was detected slightly below the band for Shank2a only from the cortex and hippocampus of e6–7 KO mice (Fig. 1d and Fig. S1). Both long (Shank2b) and short (Shank2a) Shank2 isoforms were not detected in e6–7 KO mice, which confirmed the specificity of the antibody (Fig. 1d and Fig. S1). Although expression of the novel transcript was also increased in e7 KO, it would be difficult to be detected in western blot due to its relatively lower expression level compared with that in e6–7 KO.
Characterization of the novel transcript
Unexpectedly, when we analyzed phylogenetic conservation of exon 4′ in the UCSC genome browser using comparative genomics tracks (conservation, Placentral Chain/Net, and Vertebrate Chain/Net), we found that it is well-conserved across vertebrate animals, except for fish species. As shown in Fig. 2a, there are two highly conserved regions. One of the conserved regions is located upstream of exon 4′ (left dotted black box in Fig. 2a), which could serve as the promoter for its expression. The other highly conserved region is located at the end of exon 4′ (right dotted black box in Fig. 2a), which could be the protein coding region in the novel transcript. We extracted DNA sequence that spans across the novel exon 4′ sequence, as well as 200 bp upstream of the exon 4′ sequence from eight different representative vertebrate species (except birds). Multiple sequence alignment and phylogenetic analyses were conducted. Expectedly, the overall sequence was well-aligned among the species, and the phylogenetic tree among the eight different species revealed that no evidence showing rapid evolutionary rate change was observed (Fig. S2).
We assessed the methylation profiles of regions near exon 4′ in several tissues using the DNA methylation hub track in the UCSC genome browser. Interestingly, there are overlaps between the exon 4′ region and hypomethylated regions in several tissues (Fig. 2b). Especially, the hypomethylated regions in brain tissues are associated with the tissue-specific epigenetic reconfigurations during mammalian brain development . Therefore, the novel exon 4′ is not only evolutionary conserved, but may also play functional roles in the brain.
Expression of novel transcripts in wild type mouse brain
To determine whether the novel transcript is constitutively expressed in wild type mouse brains, we investigated the public brain transcriptome data from the study by Lister et al. . They determined the temporal gene expression profiles of the mouse frontal cortex from fetal to 22-month-old mice using RNA sequencing. We collected and analyzed the RNA-Seq data for the following three data time points: fetal, 4-week-old, and 22-month-old. After constructing the Shank2 gene isoforms, as previously described, we confirmed that the novel transcripts were expressed in all time points (Fig. 3a). Furthermore, we performed RT-PCR using extracted RNAs from 4-week-old and 8-week-old wild type C57Bl/6 N mouse brain regions (prefrontal cortex, PFC; hippocampus, HPC; and cerebellum, CBL). Three different primer pairs were used to detect the known and novel transcripts (ex2f-ex4r, ex4f-ex6r: known transcript; ex4’f-ex5r: novel transcript) (Fig. 3a black arrows). As shown in Fig. 3b, we confirmed that both known and novel transcripts were expressed in each brain region.
We occasionally observe conflicting, heterogeneous phenotypes in the same genetic mutant lines associated with ASD. In addition to the two Shank2 KO lines studied here, mutant mouse lines of another Shank isoform, in which either exons 4–9 or exons 4–7 of Shank3 have been deleted, showing different synaptic and behavioral phenotypes when compared to each other [24,25,26]. Epac2 KO mice, another ASD mouse model that lacks either exon 1 or 3 also shows distinct behavioral and electrophysiological phenotypes [27,28,29]. However, the molecular mechanism underlying the different phenotypes in same gene mutant lines still remains unclear. In the previous study , we profiled whole transcriptome expressions for Shank2 e6–7 KO and e7 KO mouse models and found that expression of the Gabra2 gene was reduced in e6–7 KO mice, but not in e7 KO mice when compared to their own WT littermates; this may partially contribute to differences between Shank2 e6–7 KO and e7 KO mice.
In this study, we identified a novel Shank2 transcript, which has not been previously reported. It is rather surprising that the novel exon 4′ has not been annotated yet considering that we found that it is well-conserved across vertebrates. Relatively lower expression level of the transcript harboring the exon 4′ compared to the previously identified isoforms may have hindered the identification of this novel transcript. To validate this novel transcribed region, we computationally reconstructed transcript isoforms and experimentally confirmed the expression of the identified region using quantitative RT-PCR. In the reconstructed e6–7 KO and e7 KO transcripts, the novel transcribed region (novel exon 4′) was connected to the downstream exon 5. However, we were unable to detect the reconstructed transcripts containing the novel exon 4′ in WT samples because the RNA sequencing (RNA-Seq) depth was not enough. Furthermore, as shown in the RNA-Seq expression profile of the transcripts, expression level of the novel exon 4′ in e6–7 KO mice was greater than e7 KO mice, as determined by quantitative RT-PCR (Fig. 1b). This aberrant increase in exon 4′ expression in e6–7 KO mice may also contribute to the observed differences in the two Shank2 KO mutants, which remains to be examined. It would be interesting to examine whether ectopic expression of the cloned novel isoform in either WT or e7 KO mice can recapitulate, at least, some phenotypes in e6–7 KO.
Gene regulation is a complicated process and is finely tuned in different settings. Many studies have shown that epigenetic profiles such as methylation can affect differential isoform expressions . For example, intragenic DNA methylome map of the brain has shown that tissue-specific DNA methylation in the Shank3 gene regulates promoter activity, and that differential DNA methylation patterns are associated with alternative transcript expressions in tissue- and cell type-specific manners in the brain regions . Our study showed that regions near the novel exon 4′ were phylogenetically well-conserved across vertebrate species and were hypomethylated in several mouse tissues, including the brain tissues. We speculate that the deletion of exon 6 and 7 may have affected the methylation status of the region near the novel exon 4′, which subsequently increases the expression of the novel Shank2 transcript. Alternatively, other epigenetic modifications or chromatin reorganization, such as changes in nucleosome and positioning may have occurred in e6–7 KO. Further analyses would be required to test the suggested hypothesis.
Availability of data and materials
All data generated or analyzed during this study are included in this published article or are available from GEO DataSets (https://www.ncbi.nlm.nih.gov/gds).
Autism spectrum disorder
Real time-polymerase chain reaction
American Psychiatric Association., American Psychiatric Association. DSM-5 Task Force.: Diagnostic and statistical manual of mental disorders: DSM-5, 5th edn. Washington, D.C: American Psychiatric Association; 2013.
Devlin B, Scherer SW. Genetic architecture in autism spectrum disorder. Curr Opin Genet Dev. 2012;22(3):229–37.
Chaste P, Leboyer M. Autism risk factors: genes, environment, and gene-environment interactions. Dialogues Clin Neurosci. 2012;14(3):281–92.
Schmeisser MJ. Translational neurobiology in Shank mutant mice--model systems for neuropsychiatric disorders. Ann Anat = Anatomischer Anzeiger. 2015;200:115–7.
Choi SY, Pang K, Kim JY, Ryu JR, Kang H, Liu Z, Kim WK, Sun W, Kim H, Han K. Post-transcriptional regulation of SHANK3 expression by microRNAs related to multiple neuropsychiatric disorders. Mol Brain. 2015;8(1):74.
Monteiro P, Feng G. SHANK proteins: roles at the synapse and in autism spectrum disorder. Nat Rev Neurosci. 2017;18(3):147–57.
Yoo J, Bakes J, Bradley C, Collingridge GL, Kaang B-K. Shank mutant mice as an animal model of autism. Philos Trans R Soc B Biol Sci. 2014;369(1633):20130143.
Berkel S, Marshall CR, Weiss B, Howe J, Roeth R, Moog U, Endris V, Roberts W, Szatmari P, Pinto D, et al. Mutations in the SHANK2 synaptic scaffolding gene in autism spectrum disorder and mental retardation. Nat Genet. 2010;42(6):489–91.
Leblond CS, Heinrich J, Delorme R, Proepper C, Betancur C, Huguet G, Konyukh M, Chaste P, Ey E, Rastam M, et al. Genetic and functional analyses of SHANK2 mutations suggest a multiple hit model of autism spectrum disorders. PLoS Genet. 2012;8(2):e1002521.
Schmeisser MJ, Ey E, Wegener S, Bockmann J, Stempel AV, Kuebler A, Janssen AL, Udvardi PT, Shiban E, Spilker C, et al. Autistic-like behaviours and hyperactivity in mice lacking ProSAP1/Shank2. Nature. 2012;486(7402):256–60.
Won H, Lee HR, Gee HY, Mah W, Kim JI, Lee J, Ha S, Chung C, Jung ES, Cho YS, et al. Autistic-like social behaviour in Shank2-mutant mice improved by restoring NMDA receptor function. Nature. 2012;486(7402):261–5.
Jiang YH, Ehlers MD. Modeling autism by SHANK gene mutations in mice. Neuron. 2013;78(1):8–27.
Chung C, Ha S, Kang H, Lee J, Um SM, Yan H, Yoo YE, Yoo T, Jung H, Lee D, et al. Early correction of N-methyl-D-aspartate receptor function improves autistic-like social behaviors in adult Shank2(−/−) mice. Biol Psychiatry. 2019;85(7):534–43.
Lim CS, Kim H, Yu NK, Kang SJ, Kim T, Ko HG, Lee J, Yang JE, Ryu HH, Park T, et al. Enhancing inhibitory synaptic function reverses spatial memory deficits in Shank2 mutant mice. Neuropharmacology. 2017;112(Pt A):104–12.
Wegener S, Buschler A, Stempel AV, Kang SJ, Lim CS, Kaang BK, Shoichet SA, Manahan-Vaughan D, Schmitz D. Defective synapse maturation and enhanced synaptic plasticity in Shank2 Δex7-/- Mice. eNeuro. 2018;5(3):ENEURO.0398-17.2018.
Lee JH, Gao C, Peng G, Greer C, Ren S, Wang Y, Xiao X. Analysis of transcriptome complexity through RNA sequencing in normal and failing murine hearts. Circ Res. 2011;109(12):1332–41.
Chojnacki S, Cowley A, Lee J, Foix A, Lopez R. Programmatic access to bioinformatics tools from EMBL-EBI update: 2017. Nucleic Acids Res. 2017;45(W1):W550–3.
Kumar S, Stecher G, Li M, Knyaz C, Tamura K. MEGA X: molecular evolutionary genetics analysis across computing platforms. Mol Biol Evol. 2018;35(6):1547–9.
Song Q, Decato B, Hong EE, Zhou M, Fang F, Qu J, Garvin T, Kessler M, Zhou J, Smith AD. A reference methylome database and analysis pipeline to facilitate integrative and comparative epigenomics. PLoS One. 2013;8(12):e81148.
Wang H, Gong B, Vadakkan KI, Toyoda H, Kaang B-K, Zhuo M. Genetic evidence for adenylyl cyclase 1 as a target for preventing neuronal excitotoxicity mediated by N-methyl-D-aspartate receptors. J Biol Chem. 2007;282(2):1507–17.
Boeckers TM, Kreutz MR, Winter C, Zuschratter W, Smalla KH, Sanmarti-Vila L, Wex H, Langnaese K, Bockmann J, Garner CC, et al. Proline-rich synapse-associated protein-1/cortactin binding protein 1 (ProSAP1/CortBP1) is a PDZ-domain protein highly enriched in the postsynaptic density. J Neurosci. 1999;19(15):6506–18.
Lim S, Naisbitt S, Yoon J, Hwang JI, Suh PG, Sheng M, Kim E. Characterization of the Shank family of synaptic proteins. Multiple genes, alternative splicing, and differential expression in brain and development. J Biol Chem. 1999;274(41):29510–8.
Lister R, Mukamel EA, Nery JR, Urich M, Puddifoot CA, Johnson ND, Lucero J, Huang Y, Dwork AJ, Schultz MD, et al. Global epigenomic reconfiguration during mammalian brain development. Science (New York, NY). 2013;341(6146):1237905.
Bozdagi O, Sakurai T, Papapetrou D, Wang X, Dickstein DL, Takahashi N, Kajiwara Y, Yang M, Katz AM, Scattoni ML, et al. Haploinsufficiency of the autism-associated Shank3 gene leads to deficits in synaptic function, social interaction, and social communication. Mol Autism. 2010;1(1):15.
Peca J, Feliciano C, Ting JT, Wang W, Wells MF, Venkatraman TN, Lascola CD, Fu Z, Feng G. Shank3 mutant mice display autistic-like behaviours and striatal dysfunction. Nature. 2011;472(7344):437–42.
Wang X, McCoy PA, Rodriguiz RM, Pan Y, Je HS, Roberts AC, Kim CJ, Berrios J, Colvin JS, Bousquet-Moore D, et al. Synaptic dysfunction and abnormal behaviors in mice lacking major isoforms of Shank3. Hum Mol Genet. 2011;20(15):3093–108.
Lee K, Kobayashi Y, Seo H, Kwak JH, Masuda A, Lim CS, Lee HR, Kang SJ, Park P, Sim SE, et al. Involvement of cAMP-guanine nucleotide exchange factor II in hippocampal long-term depression and behavioral flexibility. Mol Brain. 2015;8:38.
Shibasaki T, Takahashi H, Miki T, Sunaga Y, Matsumura K, Yamanaka M, Zhang C, Tamamoto A, Satoh T, Miyazaki J, et al. Essential role of Epac2/Rap1 signaling in regulation of insulin granule dynamics by cAMP. Proc Natl Acad Sci U S A. 2007;104(49):19333–8.
Srivastava DP, Jones KA, Woolfrey KM, Burgdorf J, Russell TA, Kalmbach A, Lee H, Yang C, Bradberry MM, Wokosin D, et al. Social, communication, and cortical structural impairments in Epac2-deficient mice. J Neurosci. 2012;32(34):11864–78.
Yu N-K, Baek SH, Kaang B-K. DNA methylation-mediated control of learning and memory. Mol Brain. 2011;4(1):5.
Maunakea AK, Nagarajan RP, Bilenky M, Ballinger TJ, D'Souza C, Fouse SD, Johnson BE, Hong C, Nielsen C, Zhao Y, et al. Conserved role of intragenic DNA methylation in regulating alternative promoters. Nature. 2010;466(7303):253–7.
This work was supported by the National Honor Scientist Program (NRF-2012R1A3A1050385) through a grant to B.-K.K.; NRF-2017M3C7A1026959 to Y.-S.L.; NRF-2018H1A2A1061381 to G.P.; NRF-2017M3C9A6047623 to J.-H.L.
Ethics approval and consent to participate
Animal study was approved by the Institutional Animal Care and Use Committees at Seoul National University.
Consent for publication
The authors declare no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Western blot analysis of the Shank2 proteins. The antibody used for the analysis was same with the analysis presented in Fig. 1c. The amount of proteins loading was different and the running time was increased. Figure S2. Multiple sequence alignment and phylogenetic analysis of DNA sequences containing exon 4′ from eight vertebrate species. The novel exon 4′ region was highlighted by red boxes.
About this article
Cite this article
Lee, YS., Yu, NK., Chun, J. et al. Identification of a novel Shank2 transcriptional variant in Shank2 knockout mouse model of autism spectrum disorder. Mol Brain 13, 54 (2020). https://doi.org/10.1186/s13041-020-00595-4
- Autism spectrum disorder