Abstract
Shank2 is an abundant postsynaptic scaffolding protein that is known to regulate excitatory synapse assembly and synaptic transmission and has been implicated in various neurodevelopmental disorders, including autism spectrum disorders (ASD). Previous studies on Shank2-mutant mice provided mechanistic insights into their autistic-like phenotypes, but it remains unclear how transcriptomic patterns are changed in brain regions of the mutant mice in age- and gene dosage-dependent manners. To this end, we performed RNA-Seq analyses of the transcripts from the prefrontal cortex (PFC) of heterozygous and homozygous Shank2-mutant mice lacking exons 6 and 7 at juvenile (week 3) and adult (week 12) stages. Juvenile heterozygous Shank2-mutant mice showed upregulation of glutamate synapse-related genes, downregulation of ribosomal and mitochondrial genes, and transcriptomic changes that are opposite to those observed in ASD (anti-ASD) such as upregulation of ASD_down (downregulated in ASD), GABA neuron-related, and oligodendrocyte-related genes. Juvenile homozygous Shank2 mice showed upregulation of chromatin-related genes and transcriptomic changes that are in line with those occurring in ASD (pro-ASD) such as downregulation of ASD_down, GABA neuron-related, and oligodendrocyte-related genes. Adult heterozygous and homozygous Shank2-mutant mice both exhibited downregulation of ribosomal and mitochondrial genes and pro-ASD transcriptomic changes. Therefore, the gene dosage- and age-dependent effects of Shank2 deletions in mice include differential transcriptomic changes across distinct functional contexts, including synapses, chromatin, ribosomes, mitochondria, GABA neurons, and oligodendrocytes.
Introduction
Shank2 (also known as ProSAP1) is an abundant scaffolding protein known to regulate excitatory synapse assembly and synaptic transmission, plasticity, and signaling (Du et al., 1998; Boeckers et al., 1999; Lim et al., 1999; Naisbitt et al., 1999); reviewed in Boeckers et al. (2002); Grabrucker et al. (2011); Sala et al. (2015); Mossa et al. (2017); and Sheng and Kim (2000, 2011). Shank2 has also been implicated in various neurodevelopmental and neuropsychiatric disorders, including autism spectrum disorders (ASD), intellectual disability, developmental delay, and schizophrenia (Berkel et al., 2010, 2011; Pinto et al., 2010; Wischmeijer et al., 2011; Leblond et al., 2012, 2014; Prasad et al., 2012; Rauch et al., 2012; Sanders et al., 2012; Chilian et al., 2013; Liu et al., 2013; Schluth-Bolard et al., 2013; Guilmatre et al., 2014; Costas, 2015; Peykov et al., 2015a,b; Homann et al., 2016; Mossa et al., 2017; Yuen et al., 2017; Bai et al., 2018; Lu et al., 2018; Aspromonte et al., 2019; Zhou et al., 2019; Satterstrom et al., 2020; Wang et al., 2020). Mice carrying Shank2 mutations have enabled researchers to gain important mechanistic insights into how Shank2 deletion leads to various disease-related deficits at the molecular, synaptic, circuit, and behavioral levels (Schmeisser et al., 2012; Won et al., 2012; Ey et al., 2013, 2018; Ferhat et al., 2016; Ha et al., 2016; Ko et al., 2016; Peter et al., 2016; Lim et al., 2017; Pappas et al., 2017; Sato et al., 2017, 2020; Yoon et al., 2017; Heise et al., 2018; Kim et al., 2018; Lee S. et al., 2018; Wegener et al., 2018; Chung et al., 2019; de Chaumont et al., 2019; Han et al., 2020; Lee et al., 2020; Grabrucker et al., 2021; Horner et al., 2021); reviewed in Boeckers et al. (2002); Grabrucker et al. (2011); Jiang and Ehlers (2013); Guilmatre et al. (2014); Sala et al. (2015); Schmeisser (2015); Monteiro and Feng (2017); Mossa et al. (2017); Eltokhi et al. (2018); and Yoo et al. (2014). Many of the mechanistic studies have focused on how synaptic changes lead to behavioral abnormalities, partly because of the known functions of Shank2 at excitatory synapses. However, Shank2-mutant phenotypes likely involve various molecular and cellular changes, which should be explored comprehensively in various brain regions and cell types across different developmental stages and Shank2 gene dosages.
To begin such work, we herein analyzed the transcriptomic patterns in prefrontal cortex (PFC) regions of heterozygous (HT; Shank2+/–) and homozygous (HM; Shank2–/–) Shank2-mutant mice at juvenile (P21 or 3 weeks) and adult (P84 or 12 weeks) stages (W3-HT, W3-HM, W12-HT, and W12-HM). We chose to analyze the PFC due to its reported implications in cognitive functions, social functions, and ASD (Wang et al., 2011; Yizhar et al., 2011; Lee et al., 2016; Selimbeyoglu et al., 2017; Cao et al., 2018; Lazaro et al., 2019; Levy et al., 2019). We identified differential transcriptomic changes involving multiple biological functions and ASD-related genes, and found that they were distinct between W3-HT and W3-HM mice, but largely similar between W12-HT and W12-HM mice. These analyses also identified the potential involvement of genes associated with various functional contexts, including synapses, chromatin, ribosomes, mitochondria, GABA neurons, and oligodendrocytes.
Materials and Methods
Animals
Shank2 mice lacking exons 6 and 7 have been described (Won et al., 2012) (B6N.129S4-Shank2tm1Mgle/CsbdJ; Jackson 033667). Animals were housed under a 12-h (13:00–01:00) dark/light cycle and were fed ad libitum. The animal study was reviewed and approved by the Committee of Animal Research at KAIST (KA2020-99).
Tissue Preparation for RNA-Seq Analysis
Brains isolated from Shank2± and Shank2–/– mice (n = 5 brains from five mice for W3-HT, five for W3-HM, four for W12-HT, and four for W12-HM) were chilled in ice-cold phosphate-buffered saline solution. Mouse brains were not pooled, and samples from individual mouse brains were handled independently. Coronal brain sections (1 mm thickness) were used to dissect out PFC regions (AP + 2.8 to + 1.1 mm; ML −1 to +1 mm). The dissected tissues were kept in RNAlater solution (Ambion) at room temperature for 30 min and then kept at –70°C. Total RNA concentration was calculated by Quant-IT RiboGreen (Invitrogen, R11490). To determine the integrity of total RNAs, samples were run on the TapeStation RNA screentape (Agilent). Only high-quality RNA preparations, with RIN greater than 7.0, were used for RNA library construction. A library was prepared with 1 μg of total RNA for each sample by Illumina TruSeq mRNA Sample Prep kit (Illumina). The first step in the workflow involved purifying the poly-A-containing mRNA molecules using poly-T oligo-attached magnetic beads. Following purification, the mRNA is fragmented into small pieces using divalent cations under elevated temperature. The cleaved RNA fragments are copied into first-strand cDNAs using SuperScript II reverse transcriptase (Invitrogen) and random primers, which was followed by second-strand cDNA synthesis using DNA Polymerase I and RNase H. These cDNA fragments then went through an end repair process, the addition of a single “A” base, and then ligation of the indexing adapters. The products were then purified and enriched with PCR to create the final cDNA library. The libraries were quantified using qPCR according to the qPCR Quantification Protocol Guide (KAPA Library Quantification kits for Illumina Sequencing platforms) and qualified using the TapeStation D1000 ScreenTape (Agilent Technologies). Indexed libraries were then submitted to an Illumina NovaSeq (Illumina), and the paired-end (2 × 100 bp) sequencing was performed by Macrogen.
RNA-Seq Analysis
Transcript abundance was estimated with Salmon (v1.1.0) (Patro et al., 2017) in Quasi-mapping-based mode onto the Mus musculus genome (GRCm38) with GC bias correction (–gcBias). Quantified gene-level abundance data was imported to R (v.3.5.3) with the tximport (Soneson et al., 2015) package, and differential gene expression analysis was carried out using R/Bioconductor DEseq2 (v1.22.2) (Love et al., 2014). Some genes were filtered out when their adjusted p values were null value, failing to pass the independent filtering step in DESeq2, and when an HGNC (HUGO Gene Nomenclature Committee) human gene did not exist for a given mouse gene. For independent quality control, the percentage of mapped reads (y-axis) was plotted against the number of mapped reads (x-axis). The scatter plots were drawn using MultiQC tools and the statistics data from Salmon. Principle component analysis (PCA) was performed for the regularized log transform (rlog) of the normalized counts using plotPCA tools implemented in DEseq2. Normalized read counts were computed by dividing the raw read counts by size factors and fitted to a negative binomial distribution. The P values were adjusted for multiple testing for each developmental stage and comparisons with the Benjamini–Hochberg correction. Genes with an adjusted P value of less than 0.05 were considered as differentially expressed. Volcano plots were generated using the R ggplot2 (v.3.1.1) package. The Gene Ontology (GO) enrichment analyses were performed using DAVID software (version 6.8) (Huang da et al., 2009) on the human background. Mouse gene names were converted to human homologs using the Mouse Genome Informatics (MGI) database1.
Gene Set Enrichment Analysis
Gene Set Enrichment Analysis (GSEA)2 (Subramanian et al., 2005) was used to determine whether a priori-defined gene sets would show statistically significant differences in expression between WT and Shank2-mutant mice. Enrichment analysis was performed using GSEAPreranked (gsea-3.0.jar) module on gene set collections downloaded from Molecular Signature Database (MSigDB) v7.03. GSEAPreranked was applied using the list of all genes expressed, ranked by the fold change and multiplied by the inverse of the P value with recommended default settings (1,000 permutations and a classic scoring scheme). The False Discovery Rate (FDR), which adjusts different gene set sizes and multiple testings, was estimated to control the false-positive finding of a given Normalized Enrichment Score (NES) by comparing the tails of the observed and null distributions derived from 1,000 gene set permutations. The gene sets with an FDR of less than 0.05 were considered significantly enriched.
GSEA for biological functions were performed using the entire list of genes with transcriptional changes ranked by p values from the four transcript groups (W3-HT, W3-HM, W12-HT, and W12-HM) as inputs (Supplementary Table 1) and the Gene Ontology gene sets in the C5 subdomains at the GSEA site (2cellular component, biological process, molecular function; 10,185 gene sets as of now). For ASD-related/risk GSEA, the same inputs and the following gene sets were used (specific gene names in these gene sets are listed in Supplementary Table 9); DEG_Up_Voineagu, Co-Exp_Up_M16_Voineagu, DEG_Down_Voineagu, and Co-Exp_Down_M12_Voineagu (Voineagu et al., 2011; Werling et al., 2016), SFARI genes (all) (Abrahams et al., 2013), SFARI (high confidence, belonging to category 1) (Abrahams et al., 2013)4, FMRP targets (Darnell et al., 2011; Werling et al., 2016), DeNovoMis (protein-disrupting or missense rare de novo variants) (Iossifov et al., 2014; Werling et al., 2016), DeNovoVariants (protein-disrupting rare de novo variants) (Iossifov et al., 2014; Werling et al., 2016), and AutismKB (Autism KnowledgeBase) (Xu et al., 2012; Yang et al., 2018), and cell-type-specific gene sets (Kang et al., 2011).
Results
Differentially Expressed Gene Analyses of Heterozygous and Homozygous Shank2-Mutant Transcriptomes at Weeks 3 and 12
To determine transcriptomic changes in the PFC of heterozygous (HT; Shank2±) and homozygous (HM; Shank2–/–) Shank2-mutant mice at juvenile (3 weeks) and adult (12 weeks) stages (W3-HT, W3-HM, W12-HT, and W12-HM mice), relative to wild-type mice (W3-WT and W12-WT; five mice for W3-WT/HT/HM mice and four mice for W12-WT/HT/HM mice), we performed RNA-Seq and analyzed the obtained data. W3 and W12 transcriptomes showed an age-dependent increase in clustering, as shown by PCA (principal component analysis) plots (Supplementary Figure 1).
At week 3, eight and 50 DEGs were identified in the PFC of W3-HT and W3-HM mice, respectively, compared to WT mice at 3 weeks (Figures 1A–C and Supplementary Tables 1, 2; see W3-HT/HM tabs). Only two of the identified DEGs (Shank2 and NFIB; 3.6%) overlapped between 3W-HT and 3W-HM. At week 12, in contrast, we identified many more DEGs in 12W-HT (425; 191 upregulated and 234 downregulated) and 12W-HM (232; 83 upregulated and 149 downregulated) mice compared with W12-WT mice (Figures 1D–F and Supplementary Tables 1, 2; see W12-HT/HM tabs). Of them, 111 DEGs (20.3%) overlapped between 12W-HT and 12W-HM (see W12-HT-HM overlap tab).
FIGURE 1
DEG analyses based on RNA-Seq results from W3-HT, W3-HM, W12-HT, and W12-HM Shank2-mutant mice. (A–C) Summary table, volcano plot, and Venn diagram for DEGs from the prefrontal cortex of W3-HT and W3-HM Shank2-mutant mice, relative to WT mice. Red, upregulation; blue, downregulation (n = 5 mice for W3-WT/HT/HM; dotted lines indicate adjusted p value < 0.05, |FC| > 1.5, of which only adjusted p values but not FCs were used to define DEGs). (D–F) Summary table, volcano plot, and Venn diagram for the DEGs from the prefrontal cortex of W12-HT and W12-HM Shank2-mutant mice, relative to WT mice. Red, upregulation; blue, downregulation (n = 4 mice for W12-WT/HT/HM; dotted lines indicate adjusted p-value < 0.05, |FC| > 1.5). (G,H) Table summarizing all DEGs and a Venn diagram showing the overlaps in DEGs from W3-HT, W3-HM, W12-HT, and W12-HM Shank2-mutant mice. Red, upregulation; blue, downregulation.
A numerical comparison of the four DEG sets (Figure 1G) highlighted that there were age-dependent increases in the total numbers of DEGs in both heterozygous and homozygous Shank2 mice, with a greater increase seen in heterozygous mice. A Venn diagram plotting the DEG numbers showed that there were small overlaps between W3 and W12 DEGs within the HT and HM groups [0.2% (HT) and 5.2% (HM), respectively], compared to the greater overlaps of HT and HM DEGs at W3 (3.6%) and W12 (20.3%) (Figure 1H).
DAVID analysis of W12-HT transcripts (up- and downregulated together) revealed ribosome- and translation-related GO terms in the cellular component (CC), biological process (BP), molecular function (MF), and KEGG (Kyoto Encyclopedia of Genes and Genomes) domains (Figure 2A and Supplementary Table 3). GO terms related to the neuronal cell body, axon, and exosome were overrepresented in the CC domain but not in BP or MF domain.
FIGURE 2
DAVID gene ontology analyses of DEGs from W12-HT and W12-HM Shank2-mutant mice. (A,B) DAVID gene ontology (GO) analysis of DEGs from W12-HT and W12-HM Shank2-mutant mice (adjusted p-value < 0.05). Ribosome/translation-related GO terms are indicated in red.
DAVID analysis of W12-HM transcripts (up and downregulated together) highlighted similar GO terms (ribosome/translation-related and neuronal cell body/axon/ exosome-related) (Figure 2B). The up- and downregulated terms were pooled for these analyses to increase the strength of our analytic results.
GSEA of Heterozygous and Homozygous Shank2-Mutant Transcriptomes at 3 and 12 Weeks for Biological Functions
We additionally performed GSEA, which uses the entire list of genes with transcriptional changes (ranked by p values) rather than the relatively small number of DEGs above a particular p value or fold change, in order to identify altered biological pathways contributed by modest but coordinate changes in a more reproducible manner (Mootha et al., 2003; Subramanian et al., 2005) (see text footnote 4).
Gene Set Enrichment Analysis of W3-HT transcripts indicated positive enrichment of the ranked genes (meaning that the upregulated genes are enriched) for synapse- or excitatory postsynapse-related gene sets in the CC, BP, and MF domains of the C5 gene sets (Gene Ontology gene sets; 10,185 gene sets as of now), as shown by the top five gene sets enriched [Figure 3A; Supplementary Figures 2, 3; and Supplementary Table 4 (CC/BP/MF-Up tabs)]. Integration of all of the GSEA results (not just the top five or 10) and the use of EnrichmentMap Cytoscape App (Merico et al., 2010) to visualize clusters of biological functions further highlighted the apparent importance of synapse-related functions and (to a lesser extent) chromatin- and cell body/axon-related functions (Figure 3B). W3-HT transcripts were negatively enriched for ribosome- and mitochondria-related gene sets in the C5-CC, C5-BP, and C5-MF domains, and these findings were supported by EnrichmentMap visualization [Figures 3C,D; Supplementary Figures 2, 3; and Supplementary Table 4 (CC/BP/MF-Down tabs)]. These functional associations were generally stronger in CC and BP domains relative to the MF domain.
FIGURE 3
GSEA of transcriptomes from W3-HT and W3-HM Shank2-mutant mice for biological functions. (A–D) GSEA results for W3-HT transcripts showing the list of top five positively (red) and negatively (blue) enriched gene sets (A,C) and their integrated visualization generated using EnrichmentMap Cytoscape App (B,D). Only the results for C5-CC (cellular component) are shown here, and those for C5-BP (biological process) and C5-MF (molecular function) are shown in Supplementary Figures 2, 3. Note that, although the top five gene sets are shown in the main figure tables for simplicity (see the full lists in Supplementary Tables 4–7), EnrichmentMap visualization utilized all of the GSEA results (n = 5 mice for WT, HT, and HM, FDR < 0.05). (E–H) GSEA results for W3-HM transcripts (E,G) and their integrated visualization generated using EnrichmentMap (F,H) (n = 5 mice for WT, HT, and HM, FDR < 0.05). See Supplementary Figures 4, 5 for GSEA results for C5-BP and C5-MF.
Gene Set Enrichment Analysis of W3-HM transcripts indicated positive enrichment for chromatin- and RNA-related gene sets in the C5-CC, C5-BP, and C5-MF domains [Figures 3E,F; Supplementary Figures 4, 5; and Supplementary Table 5 (CC/BP/MF-Up tabs)]. In addition, W3-HM transcripts were, to a modest degree, negatively enriched for lipid/vacuole- and cation/anion transport-related gene sets [Figures 3G,H; Supplementary Figures 4, 5; and Supplementary Table 5 (CC/BP/MF-Down tabs)]. These results suggest that W3-HT (up and down) and W3-HM (up and down) transcripts are enriched in four distinct groups of genes.
Gene Set Enrichment Analysis of W12-HT transcripts indicated positive but weak enrichment for synapse- and chromatin-related gene sets in the C5-CC, C5-BP, and C5-MF domains [Figures 4A,B; Supplementary Figures 6, 7; and Supplementary Table 6 (CC/BP/MF-Up tabs)], which contrasts with the strong synapse-related enrichments of W3-HT transcripts. W12-HT transcripts showed strong negative enrichments for ribosome- and mitochondria-related gene sets [Figures 4C,D; Supplementary Figures 6, 7; and Supplementary Table 6 (CC/BP/MF-Down tabs)], similar to the results from W3-HT transcripts.
FIGURE 4
GSEA of transcriptomes from W12-HT and W12-HM Shank2-mutant mice for biological functions. (A–D) GSEA results for W12-HT transcripts showing the list of top five positively (red) and negatively (blue) enriched gene sets (A,C) and their integrated visualization using EnrichmentMap Cytoscape App (B,D). See Supplementary Figures 6, 7 for GSEA results for C5-BP and C5-MF (n = 4 mice for WT, HT, and HM, FDR < 0.05). (E–H) GSEA results for W12-HM transcripts (E,G) and their integrated visualization generated using EnrichmentMap (F,H) (n = 4 mice for WT, HT, and HM, FDR < 0.05). See Supplementary Figures 8, 9 for GSEA results for C5-BP and C5-MF.
W12-HM transcripts showed weakly positive enrichments for chromatin- and synapse-related gene sets [Figures 4E,F; Supplementary Figures 8, 9; and Supplementary Table 7 (CC/BP/MF-Up tabs)], similar to the results from W12-HT transcripts. In addition, W12-HM transcripts were negatively enriched for ribosome-related gene sets [Figures 4G,H; Supplementary Figures 8, 9; and Supplementary Table 7 (CC/BP/MF-Down tabs)], similar to the results from W3-HT (but not W3-HM) and W12-HT transcripts.
The abovementioned associations of the four transcript sets (W3-HT, W3-HM, W12-HT, and W12-HM) with key biological functions, including synapse, ribosome, mitochondria, and chromosome, were further highlighted by the list of specific genes that played key roles in the enrichments [Supplementary Table 8; a summary of key genes in the top gene set in each category (CC-up/down, BP-up/down, and MF-up/down domains of the four transcript sets) in the tables of Figures 3, 4]. For example, key upregulated synapse-related genes from W3-HT transcripts (CC domain) included Dmd (dystrophin), Slc8A1 (solute carrier family 8 member A1, a Na/Ca exchanger), Kcnb1 (potassium voltage-gated channel subfamily B member 1), Atp2b2 (ATPase plasma membrane Ca2+ transporting 2), and Grin2b (glutamate ionotropic receptor NMDA type subunit 2B). Key downregulated genes associated with ribosome/mitochondria-related functions from W3-HT transcripts (CC domain) included Rps27 (ribosomal protein S27), Malsu1 (mitochondrial assembly Of ribosomal large subunit 1), Larp4 (La ribonucleoprotein 4, an RNA-binding protein), Eif2d (eukaryotic translation initiation factor 2D), Nsun4 (NOP2/Sun RNA methyltransferase 4), and multiple Mrp (mitochondrial ribosomal protein) genes (i.e., MRPS6/MRPS7/MRPS5/MRPL22/MRPL4/MRPL1). Key upregulated chromatin-related genes from W3-HM transcripts included Nabp1 (nucleic acid binding protein 1), Cbx5 (Chromobox 5), Smarca1 (SWI/SNF related, matrix associated, actin dependent regulator of chromatin, subfamily A, member 1), Ago3 (Argonaute RISC catalytic component 3, with roles in RNA interference), Mef2c [myocyte enhancer factor 2C in the MADS box transcription enhancer factor 2 (MEF2) family].
GSEA of Shank2-Mutant Transcriptomes for ASD-Related/Risk Gene Sets
Previous studies have reported transcriptomic alterations in ASD (Garbett et al., 2008; Voineagu et al., 2011; Gupta et al., 2014; Parikshak et al., 2016; Velmeshev et al., 2019) and gene sets that are distinctly up- or downregulated in ASD such as DEG_Up_Voineagu, Co-Exp_Up_M16_Voineagu, DEG_Down_Voineagu, and Co-Exp_Down_M12_Voineagu (cortical samples with the age range of 2–56) (Voineagu et al., 2011; Werling et al., 2016) (see Supplementary Table 9 for the whole list of genes). In addition, multiple genes sets for ASD-risk genes have been reported; SFARI genes (all) (Abrahams et al., 2013), SFARI (high confidence, belonging to category 1) (Abrahams et al., 2013) (see text footnote 5), FMRP targets (Darnell et al., 2011; Werling et al., 2016), DeNovoMis (protein-disrupting or missense rare de novo variants) (Iossifov et al., 2014; Werling et al., 2016), DeNovoVariants (protein-disrupting rare de novo variants) (Iossifov et al., 2014; Werling et al., 2016), and AutismKB (Autism KnowledgeBase) (Xu et al., 2012; Yang et al., 2018; Supplementary Table 9). These ASD-related/risk gene sets were also found to be enriched in the transcriptomes of mice carrying ASD-risk mutations such as Chd8- and Tbr1-mutant mice in ways similar to those observed in ASD (Jung et al., 2018; Yook et al., 2019). We thus tested whether the four transcript groups (W3-HT, W3-HM, W12-HT, and W12-HM; whole list of genes ranked by p values) from Shank2-mutant mice show enrichments for these ASD-related/risk gene sets using the same GSEA approach.
In W3-HT, the transcript pattern of ASD-related genes was strongly opposite that seen in ASD patients and model mice: W3-HT transcripts were negatively and moderately enriched for gene sets that are upregulated in ASD (DEG_Up_Voineagu and Co-Exp_Up_M16_Voineagu) and positively enriched for gene sets that are downregulated in ASD (DEG_Down_Voineagu and Co-Exp_Down_M12_Voineagu) (Figure 5A). Moreover, W3-HT transcripts were positively enriched for ASD-risk gene sets that are usually downregulated in ASD likely through missense, non-sense, splice-site, frame-shift, and deletion mutations; these included SFARI genes (all), SFARI (high confidence), FMRP targets, DeNovoMis, DeNovoVariants, and AutismKB (Figure 5A). Therefore, W3-HT transcripts show patterns that are strongly opposite to those observed in ASD (here, called “anti-ASD”).
FIGURE 5
GSEA of transcriptomes from W3-HT, W3-HM, W12-HT, and W12-HM Shank2-mutant mice for ASD- and cell type-related gene sets. (A,B) GSEA results for all four transcript groups (W3-HT, W3-HM, W12-HT, and W12-HM) showing enrichment patterns for ASD- and brain disorder-related gene sets. The entire lists of genes with transcriptional changes ranked by p values were used as inputs. Note that the ASD-related gene sets include those that are upregulated in ASD (DEG_Up_ Voineagu and Co-Exp_Up_M16_ Voineagu), those that are downregulated in ASD (DEG_Down_Voineagu and Co-Exp_Down_M12_ Voineagu), and that ASD-risk genes are likely to be downregulated in ASD [SFARI (all), SFARI (high confidence), FMRP targets, DeNovoMis, DeNovoVariants, and AutismKB] (see Supplementary Table 9 for gene-set details) (n = 5 mice for W3-WT/HT/HM and four mice for W12-WT/HT/HM, FDR > 0.05). (C) GSEA results for all four transcript groups (W3-HT, W3-HM, W12-HT, and W12-HM) showing enrichment patterns for specific cell type-related gene sets (cell-type-specific; see Supplementary Table 9 for gene-set details). Note that ASD has been associated with downregulation of neuron (glutamate and GABA)- and oligodendrocyte-related genes and upregulation of glia (astrocyte and microglia)-related genes (n = 5 mice for W3-WT/HT/HM and four mice for W12-WT/HT/HM). (D) GSEA results for all four transcript groups (W3-HT, W3-HM, W12-HT, and W12-HM) showing enrichment patterns for cell-type-specific genes validated by single-cell RNA-Seq results of human cortical tissues (single-cell-type-specific; see Supplementary Table 9 for gene-set details) (n = 5 mice for W3-WT/HT/HM and four mice for W12-WT/HT/HM).
W3-HM transcripts were negatively enriched for gene sets that are upregulated in ASD (DEG_Up_ Voineagu and Co-Exp_Up_M16_ Voineagu) (Figure 5A), similar to the pattern in W3-HT transcripts (anti-ASD). However, W3-HM transcripts were negatively and strongly enriched for gene sets downregulated in ASD (DEG_Down_Voineagu and Co-Exp_Down_M12_ Voineagu) (Figure 5A), and thus also exhibited some “pro-ASD” transcription. W3-HM transcripts were positively enriched for SFARI genes (all), SFARI genes (high confidence), DeNovoMis, DeNovoVariants, and AutismKB but negatively enriched for FMRP targets (Figure 5A). These results from W3-HM mice partly deviate from the W3-HT pattern and suggest that an increase in the dosage of Shank2 deletion dampens the anti-ASD transcriptomic changes at W3.
W12-HT and W12-HM transcripts showed patterns in line with anti-ASD changes. The W12-HT/HM transcripts were negatively enriched for gene sets upregulated in ASD (DEG_Up_ Voineagu and Co-Exp_Up_M16_ Voineagu), although they were negatively enriched for some gene sets downregulated in ASD (DEG_Down_Voineagu and Co-Exp_Down_M12_ Voineagu). W12-HT/HM transcripts were also positively enriched for SFARI genes (all), SFARI (high confidence), FMRP targets, DeNovoMis, DeNovoVariants, and AutismKB (Figure 5A). These patterns are more similar to that of W3-HM than W3-HT.
The four transcript groups (W3-HT, W3-HM, W12-HT, and W12-HM) were minimally enriched for other brain disorder-related gene sets, compared with ASD-related/risk gene sets, in DisGeNet5 (Pinero et al., 2020; Figure 5B). One exception was the relatively strong enrichment of W3-HT, but not the other three transcripts, for genes related to epilepsy, in line with that (1) epilepsy is one of the comorbidities of ASD, (2) epilepsy symptoms were observed in individuals with SHANK2-related genetic variations (although rare) (Leblond et al., 2012), and (3) Shank2 expression is altered in epileptic humans and rat models (Fu et al., 2020). These results collectively suggest that W3-HT transcripts show unique and strong anti-ASD transcriptomic changes, whereas the other three transcript sets (W3-HM, W12-HT, and W12-HM) show both anti-ASD and pro-ASD changes.
GSEA of Shank2-Mutant Transcriptomes for Cell Type-Specific Gene Sets
Autism spectrum disorders has been associated with cell-type-specific changes in gene expression, including decreased neuronal gene expression and increased/decreased glial cell expression (increased astrocyte/microglia-related genes and decreased oligodendrocyte-related genes) (Voineagu et al., 2011; Werling et al., 2016). We thus tested if Shank2 transcripts from W3/12-HT/HM mice show any cell-type-specific enrichments, using cell-type-specific gene sets (Kang et al., 2011) and single-nucleus cell-type-specific gene sets (Velmeshev et al., 2019) (see Supplementary Table 9 for further details).
W3-HT transcripts were positively enriched for cortical glutamate neuron-related gene sets identified across different cortical layers (Figure 5C; Kang et al., 2011). Similarly, W3-HT transcripts were largely positively enriched for GABA neuron-related gene sets, including those for parvalbumin-positive neurons (Figure 5C; Kang et al., 2011). For glia-related gene sets, W3-HT transcripts were positively and strongly enriched for oligodendrocyte (mature)-, but not for astrocyte- or microglia-related gene sets.
W3-HM transcripts were also positively enriched for cortical glutamate neuron-related gene sets, similar to W3-HT transcripts; however, they were largely negatively enriched for GABA neuron-related gene sets, including those associated with parvalbumin-positive neurons (Figure 5C). W3-HM transcripts were negatively enriched for all three glia-related gene sets (astrocyte, microglia, and oligodendrocytes).
W12-HT and W12-HM transcripts were positively enriched for cortical glutamate neuron-related gene sets and largely negatively enriched for GABA neuron-related gene sets (Figure 5C). In addition, W12-HT and W12-HM transcripts were negatively enriched for glia-related gene sets (astrocytes, microglia, and oligodendrocytes).
Therefore, W3-HT transcripts are strongly anti-ASD based on upregulated neuronal (both glutamate and GABA) and oligodendrocyte genes, whereas the other three transcripts (W3-HM, W12-HT, and W12-HM) are partly anti-ASD (upregulated glutamate neuronal genes and downregulated astrocyte/microglia genes) and pro-ASD (downregulated oligodendrocyte genes).
Lastly, the four transcript groups (W3-HT, W3-HM, W12-HT, and W12-HM) were tested cell-type-specific genes validated by single-cell RNA-Seq results of human cortical tissues (Velmeshev et al., 2019). The overall enrichment patterns were similar to those described above (Figure 5C) in that W3-HT transcripts were uniquely and strongly anti-ASD, with upregulated neuronal and oligodendrocyte genes, and that the other three transcripts were less strongly anti-ASD, as shown by moderately upregulated neuronal and oligodendrocyte genes and largely downregulated astrocyte/microglia genes (Figure 5D).
Discussion
In the present study, we analyzed and compared the transcriptomic changes occurring in the PFC of four different Shank2-mutant mouse groups (W3-HT, W3-HM, W12-HT, and W12-HM) using DEG and GSEA analyses. Our goal was to provide insight into the molecular changes induced in the PFC by Shank2 heterozygous and homozygous deletions across juvenile and adult stages in mice.
Our analyses identified differential transcriptomic changes in heterozygous and homozygous Shank2-mutant mice across juvenile and adult stages and interesting biological functions, including those associated with synapse, ribosome, mitochondria, chromatin, GABA neuron, and oligodendrocyte. Before interpreting these results, the following limitations of the study should be pointed out. First, the fact that some biological functions are missing in a subset of the four transcript groups does not necessarily mean that they do not exist. They could be present, but the statistical powers may be low; we used only 4–5 mice per group. Second, the biological functions identified in the present study involving up- or downregulations of the genes are those that work only at the transcript level and are not supported by any functional results. In addition, the transcriptional changes may not reflect functional up- or downregulations but merely reflect compensatory transcriptional changes, which often move toward opposite directions, although additional hints could come from whether the changes are anti- or pro-ASD (see below).
One interesting transcriptomic change was the increase in synapse (particularly excitatory postsynapse)-related gene expression in W3-HT mice, which was not observed in W3-HM, W12-HT, or W12-HM mice (Figures 3, 4). Given that Shank2 is a core component of excitatory postsynaptic protein complexes (Kim and Sheng, 2004; Grabrucker et al., 2011; Sala et al., 2015; Mossa et al., 2017), and the loss of Shank2 exons 6 and 7 suppresses excitatory synaptic transmission mediated by NMDARs in the mPFC, hippocampus, and amygdala (Won et al., 2012; Lee et al., 2015; Chung et al., 2019), and that the W3-HT transcriptome has anti-ASD nature, the upregulation of excitatory synapse-related genes in W3-HT mice is likely to represent indirect changes that compensate for reduced excitatory transmission. However, although W3-HM mice show decreased NMDAR function in the mPFC, the previous studies did not directly test if W3-HT mice show reduced excitatory synaptic transmission in the PFC or hippocampus (Won et al., 2012; Chung et al., 2019).
Although W3-HM mice did not display altered synapse-related gene expression, they did show increased chromatin-related gene expression. This suggests that a stronger Shank2 deletion (HM relative to HT) converts the “color” of the transcriptomic changes from synapse- to chromatin-related gene expression changes. It is tempting to speculate that broader changes in the expression levels of chromatin remodeling-related genes are required to compensate for the deficits induced by a stronger Shank2 deletion. In line with this, human cortical neurons carrying SHANK2 mutations with excitatory synaptic hyperconnectivity show chromatin-related changes in gene expression (Zaslavsky et al., 2019). In addition, chromatin-related changes have frequently been associated with ASD (De Rubeis et al., 2014) and observed in mice lacking Shank3 (a relative of Shank2) (Zhu et al., 2014; Ma et al., 2018; Qin et al., 2018; Wang et al., 2019; Zhang et al., 2021).
The increased synaptic gene expression observed in W3-HT mice is markedly weakened in W12-HT mice, where only moderate increases in synapse- and chromatin-related gene expression were observed. It is possible that the continued absence of Shank2 expression across juvenile and adult stages dampens the compensatory gene expression changes. In addition, it could be that the period of active synapse development and flexible synapse-related transcriptomic changes comes to an end roughly around the third postnatal week and does not persist into adulthood (Sheng and Kim, 2011). Moreover, the prolonged heterozygous Shank2 deletion might have changed the nature of the compensatory responses from synapse-only changes to those involving both synapse- and chromatin-related gene expressions. Notably, W12-HM mice show transcript patterns that are similar to those of W12-HT mice (moderate increases in synapse- and chromatin-related gene expressions), suggesting that prolonged Shank2 deletion weakens the dosage-dependent effects of Shank2 loss at W12, which were strong at W3 between HT and HM mice. This could be attributable to age-dependent changes or some compensatory changes such as altered expression of other members of the Shank family (i.e., Shank3), as reported previously (Schmeisser et al., 2012).
Another notable transcriptomic change observed in W3-HT mice was the downregulations of ribosome- and mitochondria-related genes. This change was not observed in W3-HM mice, where only moderate decreases in lipid- and cation transport-related gene expression levels were observed. This suggests that Shank2 deletions of different dosages induce distinct changes in gene downregulations at W3, similar to the distinct upregulatory patterns seen at W3 (synapse- and chromatin-related genes in W3-HT and W3-HM mice, respectively). Notably, equally strong decreases in ribosome/mitochondria-related gene expression levels were observed in W12-HT mice, suggesting that, unlike the substantial dampening of synapse gene upregulations across W3 and W12, there is no age-dependent dampening of transcriptomic downregulations. In addition, W12-HM mice also showed strong decreases in ribosome-related gene expressions. It thus seems that decreases in ribosome/mitochondria-related gene expression levels represent the strongest and most widespread transcriptomic changes observed in three out of the four transcript sets (W3-HT, W12-HT, and W12-HM, but not W3-HM). Although we named some of the genes associated with ribosome/mitochondrial functions in Results based on the relative contributions to the enrichments, caution should be taken in narrowing down to specific genes because there are multiple gene sets associated with a single biological function (i.e., ribosome), and there are multiple genes in a single gene set that contributed to the enrichments. It is possible that moderate but coordinated transcriptional changes in many genes can lead to strong gene-set enrichments, which might be occurring in ribosome/mitochondria genes in Shank2-mutant mice.
It is unclear how ribosome/translation-related genes are downregulated in Shank2-mutant mice. Protein synthesis is thought to be reciprocally connected with synaptic proteins and functions under the context of ASD-risk gene functions (Santini and Klann, 2014). In addition, Shank2 is known to coordinate excitatory synaptic signaling pathways (Grabrucker et al., 2011; Sheng and Kim, 2011; Sala et al., 2015; Mossa et al., 2017). It is, therefore, possible that suppressed NMDAR or excitatory synaptic functions in Shank2-mutant mice are associated with the decreased ribosome/translation-related transcriptomic activities, perhaps through mTOR signaling, a central regulator of translational initiation and synaptic and neuronal functions implicated in various brain disorders, including ASD (Hoeffer and Klann, 2010; Costa-Mattioli and Monteggia, 2013; Borrie et al., 2017; Saxton and Sabatini, 2017; Switon et al., 2017; Winden et al., 2018; Bagni and Zukin, 2019; Sossin and Costa-Mattioli, 2019) and known to be regulated by synaptic receptors, including NMDARs and metabotropic glutamate receptors (mGluRs) (Hoeffer and Klann, 2010). However, our GSEA on Shank2-mutant transcriptomes did not detect any significant changes in mTOR signaling-related gene expressions, although altered mTOR signaling should be directly tested at the protein level using protein phosphorylation/activity as readouts. Intriguingly, bi-allelic but not mono-allelic deletion of Shank2 has been shown to increase Akt phosphorylation without changes in total protein levels during early neuronal differentiation in SH-SY5Y cells, which involved decreased apoptosis and increased cell proliferation (Unsicker et al., 2021), suggesting that Shank2 regulates mTOR signaling downstream of receptor tyrosine kinase activation.
Intriguingly, mitochondria-related gene expressions are changed largely in parallel with those of ribosome/translation-related genes in W3-HT, W12-HT, and W12-HM mice. Mitochondrial dysfunctions have been extensively implicated in ASD (Lombard, 1998; Oliveira et al., 2005; Weissman et al., 2008; Giulivi et al., 2010; Frye and Rossignol, 2011; Anitha et al., 2012; Rossignol and Frye, 2012; Frye et al., 2013; Legido et al., 2013; Goh et al., 2014; Siddiqui et al., 2016; Hollis et al., 2017; Khemakhem et al., 2017; El-Ansary et al., 2018; Rose et al., 2018; Frye, 2020; Singh et al., 2020). In addition, mitochondria critically regulate synapse formation and synaptic transmission and plasticity (Li et al., 2004, 2020; Vos et al., 2010; Sheng and Cai, 2012; Lee A. et al., 2018; Rossi and Pekkurnaz, 2019). It is possible that excitatory synaptic defects induced by Shank2 deletion, likely accompanying altered mTOR and ribosome/translation functions, may affect mitochondrial functions. Indeed, mTOR signaling regulates mitochondrial functions, in addition to protein synthesis, in mammalian cells to coordinate the energy consumption involved in protein translation with the energy production by mitochondria (Morita et al., 2015). Alternatively, mitochondria carry a few own genes required for oxidative phosphorylation, and mitochondrial ribosomes (mitoribosomes) are required for their translation (Amunts et al., 2015; Greber et al., 2015). Notably, downregulated ribosome-related genes from W3-HM transcripts included those for multiple MRPs (nuclear-encoded mitochondrial ribosomal, or mitoribosomal, proteins), known to regulate mitoribosome biogenesis and mitochondrial translation implicated in human diseases (Smits et al., 2010; Boczonadi and Horvath, 2014), and Malsu1, a regulator of mitoribosome biogenesis and mitochondrial translation (Rorbach et al., 2012; Wanschers et al., 2012; Fung et al., 2013). It is possible that mitochondrial dysfunctions in Shank2-mutant mice might induce downregulations of genes associated with mitoribosome biogenesis and mitochondrial translation.
GSEA for ASD-related gene sets indicate that W3-HT transcripts show enrichment patterns that are strongly opposite to those observed in ASD (Figure 5). The gene sets that are downregulated in ASD are negatively enriched (pro-ASD) in W3-HM, W12-HT, and W12-HM, but positively enriched in W3-HT (anti-ASD). In addition, the positive enrichment of W3-HT transcripts for synapse-related gene sets is opposite to the finding that synapse-related gene functions are downregulated in ASD post-mortem cortex (Voineagu et al., 2011). In these contexts, therefore, W3-HT transcripts again show the strongest and most unique anti-ASD patterns.
The four transcript groups are also similarly enriched for the ASD-related cell-type gene sets. W3-HT transcripts show positive enrichment for glutamate neuron-related gene sets (anti-ASD), and similar results are seen in the other three groups (W3-HM, W12-HT, and W12-HM). However, W3-HT transcripts are positively enriched for GABA neuron-related gene sets (anti-ASD), whereas the other three groups show negative enrichments, or less strong positive enrichments, for GABA neuronal genes (pro-ASD). Given the decreased neuronal gene expression in ASD (Voineagu et al., 2011) and mouse models of ASD (Jung et al., 2018; Yook et al., 2019), these results, together with the abovementioned results, indicate that W3-HT transcripts are uniquely and strongly anti-ASD. These early and strong transcriptomic changes might underlie the lack of strong autistic-like behavioral phenotypes in heterozygous Shank2-mutant mice, although hyperactivity was observed in these mice (Won et al., 2012). In addition, these results further suggest that decreased GABA neuron-related gene expression and consequent GABA neuronal dysfunction might promote ASD-like phenotypes in Shank2-mutant mice. Consistent with this possibility, parvalbumin-positive neuron-specific deletion of Shank2 (exon 6 + 7) leads to moderate hyperactivity, enhanced self-grooming, and suppressed seizure susceptibility in mice (Lee S. et al., 2018). In addition, promoting GABAergic synaptic transmission rescues cognitive deficits observed in adult homozygous Shank2-mutant mice (Lim et al., 2017).
W3-HT transcripts show strong positive enrichment for oligodendrocyte-related gene sets, whereas the other three groups show negative, or only moderate, enrichments. Given that downregulation of oligodendrocyte-related gene sets, but not those related to astrocytes or microglia, has been implicated in ASD (Voineagu et al., 2011), our present results suggest that W3-HT transcripts but not the other three are anti-ASD patterns. Whether this difference differentially affects the density and function of oligodendrocytes in Shank2-mutant mice remains to be determined, although abnormalities in oligodendrocytes and axonal myelination have been associated with ASD (Graciarena et al., 2018; Galvez-Contreras et al., 2020; Kawamura et al., 2020a,b).
In summary, our study demonstrates that heterozygous and homozygous Shank2 (exons 6 + 7) deletions in mice at W3 and W12 induce differential transcriptomic changes, and that these are distinct at W3 but become more similar at W12. Our study also identifies distinct functional groups of genes whose expression changes are associated with Shank2 deletions in mice, such as those related to synapses, chromatin, ribosomes, mitochondria, GABA neurons, and oligodendrocytes.
Statements
Data availability statement
The RNA-Seq datasets presented in this study can be found in the GEO (Gene Expression Omnibus) repository under Accession Number (GSE168211).
Ethics statement
The animal study was reviewed and approved by the committee of animal research at KAIST (KA2020-99).
Author contributions
SL, EL, HJ, and HK performed RNA-Seq analysis. EK and EL designed research and wrote the manuscript. All authors contributed to the article and approved the submitted version.
Funding
This work was supported by the Korea Institute of Science and Technology Information (K-19-L02-C07-S01 to HK), Samjung-Dalim Faculty Research Grant of Yonsei University College of Medicine (6-2019-0127 to EL), Korea National Research Foundation (NRF-2019R1A2C1084812 to EL), and the Institute for Basic Science (IBS-R002-D1 to EK).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnmol.2021.683196/full#supplementary-material
Supplementary Figure 1Principle component analysis (PCA) and scatter plots for W3-WT, W3-HT, W3-HM, W12-WT, W12-HT, and W12-HM transcriptomes. (A,B) PCA plots for WT, HT, and HM transcriptomes at W3 (A) and W12 (B). (C,D) Scatter plots for WT, HT, and HM transcriptomes combined at W3 (C) and W12 (D), for independent quality control of the sequencing results, showing the percentage of mapped reads plotted against the number of mapped reads, which ranged 25–36 million and ∼80%.
Supplementary Figure 2GSEA of transcriptomes from W3-HT Shank2-mutant mice for biological functions in the C5-BP (biological process) domain. (A,B) GSEA results for W3-HT transcripts showing the list of top five positively (red) and negatively (blue) enriched gene sets (A) and their integrated visualization generated using Cytoscape EnrichmentMap App (B) (n = 5 mice for WT, HT, and HM, FDR < 0.05).
Supplementary Figure 3GSEA of transcriptomes from W3-HT Shank2-mutant mice for biological functions in the C5-MF (molecular function) domain. (A,B) GSEA results for W3-HT transcripts showing the list of top five positively (red) and negatively (blue) enriched gene sets (A) and their integrated visualization generated using Cytoscape EnrichmentMap App (B) (n = 5 mice for WT, HT, and HM, FDR < 0.05).
Supplementary Figure 4GSEA of transcriptomes from W3-HM Shank2-mutant mice for biological functions in the C5-BP (biological process) domain. (A,B) GSEA results for W3-HM transcripts (A) and their integrated visualization generated using Cytoscape EnrichmentMap App (B) (n = 5 mice for WT, HT, and HM, FDR < 0.05).
Supplementary Figure 5GSEA of transcriptomes from W3-HM Shank2-mutant mice for biological functions in the C5-MF (molecular function) domain. (A,B) GSEA results for W3-HM transcripts (A) and their integrated visualization generated using Cytoscape EnrichmentMap App (B) (n = 5 mice for WT, HT, and HM, FDR < 0.05).
Supplementary Figure 6GSEA of transcriptomes from W12-HT Shank2-mutant mice for biological functions in the C5-BP (biological process) domain. (A,B) GSEA results for W12-HT transcripts showing the list of top five positively (red) and negatively (blue) enriched gene sets (A) and their integrated visualization generated using Cytoscape EnrichmentMap App (B) (n = 4 mice for WT, HT, and HM, FDR < 0.05).
Supplementary Figure 7GSEA of transcriptomes from W12-HT Shank2-mutant mice for biological functions in the C5-MF (molecular function) domain. (A,B) GSEA results for W12-HT transcripts showing the list of top five positively (red) and negatively (blue) enriched gene sets (A) and their integrated visualization generated using Cytoscape EnrichmentMap App (B) (n = 4 mice for WT, HT, and HM, FDR < 0.05).
Supplementary Figure 8GSEA of transcriptomes from W12-HM Shank2-mutant mice for biological functions in the C5-BP (biological process) domain. (A,B) GSEA results for W12-HM transcripts (A) and their integrated visualization generated using Cytoscape EnrichmentMap App (B) (n = 4 mice for WT, HT, and HM, FDR < 0.05).
Supplementary Figure 9GSEA of transcriptomes from W12-HM Shank2-mutant mice for biological functions in the C5-MF (molecular function) domain. (A,B) GSEA results for W12-HM transcripts (A) and their integrated visualization generated using Cytoscape EnrichmentMap App (B) (n = 4 mice for WT, HT, and HM, FDR < 0.05).
Supplementary Table 1Whole RNA-Seq results for the prefrontal cortex of W3-HT, W3-HM, W12-HT, and W12-HM Shank2-mutant mice.
Supplementary Table 2List of upregulated and downregulated DEGs in W3-HT, W3-HM, W12-HT, and W12-HM Shank2-mutant mice.
Supplementary Table 3DAVID analysis results (CC, BP, MF, and KEGG domains) for DEGs in W12-HT and W12-HM Shank2-mutant transcripts.
Supplementary Table 4GSEA results (C5-CC, C5-BP, and C5-MF domains) for W3-HT Shank2-mutant transcripts.
Supplementary Table 5GSEA results (C5-CC, C5-BP, and C5-MF domains) for W3-HM Shank2-mutant transcripts.
Supplementary Table 6GSEA results (C5-CC, C5-BP, and C5-MF domains) for W12-HT Shank2-mutant transcripts.
Supplementary Table 7GSEA results (C5-CC, C5-BP, and C5-MF domains) for W12-HM Shank2-mutant transcripts.
Supplementary Table 8Details on the top enriched gene sets in C5 subdomains (CC, BP, and MF) for W3-HT, W3-HM, W12-HT, and W12-HM transcript groups, showing information on individual genes contributed to the enrichments such as original rank and running enrichment score.
Supplementary Table 9ASD-related/risk and cell-type-specific gene sets used for GSEA of Shank2-mutant transcript groups (W3-HT, W3-HM, W12-HT, and W12-HM).
Abbreviations
- ASD
autism spectrum disorders
- DEG
differentially expressed genes
- GO
gene ontology
- GSEA
gene set enrichment analysis
- CC
cellular component
- BP
biological process
- MF
molecular function
- SFARI
Simons Foundation Autism Research Initiative
- AutismKB
Autism KnowledgeBase
- KEGG
Kyoto Encyclopedia of Genes and Genomes
- GABA
gamma amino butyric acid
Footnotes
1.^ http://www.informatics.jax.org/homology.shtml
2.^ http://software.broadinstitute.org/gsea
References
1
Abrahams B. S. Arking D. E. Campbell D. B. Mefford H. C. Morrow E. M. Weiss L. A. et al (2013). SFARI Gene 2.0: a community-driven knowledgebase for the autism spectrum disorders (ASDs).Mol. Autism4:36. 10.1186/2040-2392-4-36
2
Amunts A. Brown A. Toots J. Scheres S. H. W. Ramakrishnan V. (2015). Ribosome. The structure of the human mitochondrial ribosome.Science34895–98. 10.1126/science.aaa1193
3
Anitha A. Nakamura K. Thanseem I. Yamada K. Iwayama Y. Toyota T. et al (2012). Brain region-specific altered expression and association of mitochondria-related genes in autism.Mol. Autism3:12. 10.1186/2040-2392-3-12
4
Aspromonte M. C. Bellini M. Gasparini A. Carraro M. Bettella E. Polli R. et al (2019). Characterization of intellectual disability and autism comorbidity through gene panel sequencing.Hum. Mutat.401346–1363. 10.1002/humu.23822
5
Bagni C. Zukin R. S. (2019). A synaptic perspective of fragile X syndrome and autism spectrum disorders.Neuron1011070–1088. 10.1016/j.neuron.2019.02.041
6
Bai Y. Qiu S. Li Y. Li Y. Zhong W. Shi M. et al (2018). Genetic association between SHANK2 polymorphisms and susceptibility to autism spectrum disorder.IUBMB life70763–776. 10.1002/iub.1876
7
Berkel S. Marshall C. R. Weiss B. Howe J. Roeth R. Moog U. et al (2010). Mutations in the SHANK2 synaptic scaffolding gene in autism spectrum disorder and mental retardation.Nat. Genet.42489–491. 10.1038/ng.589
8
Berkel S. Tang W. Trevino M. Vogt M. Obenhaus H. A. Gass P. et al (2011). Inherited and de novo SHANK2 variants associated with autism spectrum disorder impair neuronal morphogenesis and physiology.Hum. Mol. Genet.21344–357. 10.1093/hmg/ddr470
9
Boczonadi V. Horvath R. (2014). Mitochondria: impaired mitochondrial translation in human disease.Int. J. Biochem. Cell Biol.4877–84. 10.1016/j.biocel.2013.12.011
10
Boeckers T. M. Bockmann J. Kreutz M. R. Gundelfinger E. D. (2002). ProSAP/Shank proteins - a family of higher order organizing molecules of the postsynaptic density with an emerging role in human neurological disease.J. Neurochem.81903–910. 10.1046/j.1471-4159.2002.00931.x
11
Boeckers T. M. Kreutz M. R. Winter C. Zuschratter W. Smalla K. H. Sanmarti-Vila L. et al (1999). 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.196506–6518. 10.1523/JNEUROSCI.19-15-06506.1999
12
Borrie S. C. Brems H. Legius E. Bagni C. (2017). Cognitive dysfunctions in intellectual disabilities: the contributions of the Ras-MAPK and PI3K-AKT-mTOR pathways.Annu. Rev. Genomics Hum. Genet.18115–142. 10.1146/annurev-genom-091416-035332
13
Cao W. Lin S. Xia Q. Q. Du Y. L. Yang Q. Zhang M. Y. et al (2018). Gamma oscillation dysfunction in mPFC leads to social deficits in neuroligin 3 R451C knockin mice.Neuron971253.e7–1260.e7. 10.1016/j.neuron.2018.02.001
14
Chilian B. Abdollahpour H. Bierhals T. Haltrich I. Fekete G. Nagel I. et al (2013). Dysfunction of SHANK2 and CHRNA7 in a patient with intellectual disability and language impairment supports genetic epistasis of the two loci.Clin. Genet.84560–565. 10.1111/cge.12105
15
Chung C. Ha S. Kang H. Lee J. Um S. M. Yan H. et al (2019). Early correction of N-Methyl-D-aspartate receptor function improves autistic-like social behaviors in adult Shank2(-/-) Mice.Biol. Psychiatry85534–543. 10.1016/j.biopsych.2018.09.025
16
Costa-Mattioli M. Monteggia L. M. (2013). mTOR complexes in neurodevelopmental and neuropsychiatric disorders.Nat. Neurosci.161537–1543. 10.1038/nn.3546
17
Costas J. (2015). The role of SHANK2 rare variants in schizophrenia susceptibility.Mol. Psychiatry20:1486. 10.1038/mp.2015.119
18
Darnell J. C. Van Driesche S. J. Zhang C. Hung K. Y. Mele A. Fraser C. E. et al (2011). FMRP stalls ribosomal translocation on mRNAs linked to synaptic function and autism.Cell146247–261. 10.1016/j.cell.2011.06.013
19
de Chaumont F. Ey E. Torquet N. Lagache T. Dallongeville S. Imbert A. et al (2019). Real-time analysis of the behaviour of groups of mice via a depth-sensing camera and machine learning.Nat. Biomed. Eng.3930–942. 10.1038/s41551-019-0396-1
20
De Rubeis S. He X. Goldberg A. P. Poultney C. S. Samocha K. Cicek A. E. et al (2014). Synaptic, transcriptional and chromatin genes disrupted in autism.Nature515209–215. 10.1038/nature13772
21
Du Y. Weed S. A. Xiong W. C. Marshall T. D. Parsons J. T. (1998). Identification of a novel cortactin SH3 domain-binding protein and its localization to growth cones of cultured neurons.Mol. Cell Biol.185838–5851. 10.1128/MCB.18.10.5838
22
El-Ansary A. Bjorklund G. Khemakhem A. M. Al-Ayadhi L. Chirumbolo S. Ben Bacha A. (2018). Metabolism-associated markers and childhood autism rating scales (CARS) as a measure of autism severity.J. Mol. Neurosci.65265–276. 10.1007/s12031-018-1091-5
23
Eltokhi A. Rappold G. Sprengel R. (2018). Distinct phenotypes of Shank2 mouse models reflect neuropsychiatric spectrum disorders of human patients with SHANK2 variants.Front. Mol. Neurosci.11:240. 10.3389/fnmol.2018.00240
24
Ey E. Torquet N. de Chaumont F. Levi-Strauss J. Ferhat A. T. Le Sourd A. M. et al (2018). Shank2 mutant mice display hyperactivity insensitive to methylphenidate and reduced flexibility in social motivation, but normal social recognition.Front. Mol. Neurosci.11:365. 10.3389/fnmol.2018.00365
25
Ey E. Torquet N. Le Sourd A. M. Leblond C. S. Boeckers T. M. Faure P. et al (2013). The Autism ProSAP1/Shank2 mouse model displays quantitative and structural abnormalities in ultrasonic vocalisations.Behav. Brain Res.256677–689. 10.1016/j.bbr.2013.08.031
26
Ferhat A. T. Torquet N. Le Sourd A. M. de Chaumont F. Olivo-Marin J. C. Faure P. et al (2016). Recording mouse ultrasonic vocalizations to evaluate social communication.J. Visual. Exp.112:53871. 10.3791/53871
27
Frye R. E. (2020). Mitochondrial dysfunction in autism spectrum disorder: unique abnormalities and targeted treatments.Semin. Pediatr. Neurol.35:100829. 10.1016/j.spen.2020.100829
28
Frye R. E. Melnyk S. Macfabe D. F. (2013). Unique acyl-carnitine profiles are potential biomarkers for acquired mitochondrial disease in autism spectrum disorder.Transl. Psychiatry3:e220. 10.1038/tp.2012.143
29
Frye R. E. Rossignol D. A. (2011). Mitochondrial dysfunction can connect the diverse medical symptoms associated with autism spectrum disorders.Pediatr. Res.6941R–47R. 10.1203/PDR.0b013e318212f16b
30
Fu Y. Liu D. Guo J. Long H. Xiao W. Xiao W. et al (2020). Dynamic change of shanks gene mRNA expression and DNA methylation in epileptic rat model and human patients.Mol. Neurobiol.573712–3726. 10.1007/s12035-020-01968-5
31
Fung S. Nishimura T. Sasarman F. Shoubridge E. A. (2013). The conserved interaction of C7orf30 with MRPL14 promotes biogenesis of the mitochondrial large ribosomal subunit and mitochondrial translation.Mol. Biol. Cell24184–193. 10.1091/mbc.e12-09-0651
32
Galvez-Contreras A. Y. Zarate-Lopez D. Torres-Chavez A. L. Gonzalez-Perez O. (2020). Role of oligodendrocytes and myelin in the pathophysiology of autism spectrum disorder.Brain Sci10:251. 10.3390/brainsci10120951
33
Garbett K. Ebert P. J. Mitchell A. Lintas C. Manzi B. Mirnics K. et al (2008). Immune transcriptome alterations in the temporal cortex of subjects with autism.Neurobiol. Dis.30303–311. 10.1016/j.nbd.2008.01.012
34
Giulivi C. Zhang Y. F. Omanska-Klusek A. Ross-Inta C. Wong S. Hertz-Picciotto I. et al (2010). Mitochondrial dysfunction in autism.JAMA3042389–2396. 10.1001/jama.2010.1706
35
Goh S. Dong Z. Zhang Y. DiMauro S. Peterson B. S. (2014). Mitochondrial dysfunction as a neurobiological subtype of autism spectrum disorder: evidence from brain imaging.JAMA Psychiatry71665–671. 10.1001/jamapsychiatry.2014.179
36
Grabrucker A. M. Schmeisser M. J. Schoen M. Boeckers T. M. (2011). Postsynaptic ProSAP/Shank scaffolds in the cross-hair of synaptopathies.Trends Cell Biol.21594–603. 10.1016/j.tcb.2011.07.003
37
Grabrucker S. Pagano J. Schweizer J. Urrutia-Ruiz C. Schon M. Thome K. et al (2021). Activation of the medial preoptic area (MPOA) ameliorates loss of maternal behavior in a Shank2 mouse model for autism.EMBO J.40:e104267. 10.15252/embj.2019104267
38
Graciarena M. Seiffe A. Nait-Oumesmar B. Depino A. M. (2018). Hypomyelination and oligodendroglial alterations in a mouse model of autism spectrum disorder.Front. Cell. Neurosci.12:517. 10.3389/fncel.2018.00517
39
Greber B. J. Bieri P. Leibundgut M. Leitner A. Aebersold R. Boehringer D. et al (2015). Ribosome. The complete structure of the 55S mammalian mitochondrial ribosome.Science348303–308. 10.1126/science.aaa3872
40
Guilmatre A. Huguet G. Delorme R. Bourgeron T. (2014). The emerging role of SHANK genes in neuropsychiatric disorders.Dev. Neurobiol.74113–122. 10.1002/dneu.22128
41
Gupta S. Ellis S. E. Ashar F. N. Moes A. Bader J. S. Zhan J. et al (2014). Transcriptome analysis reveals dysregulation of innate immune response genes and neuronal activity-dependent genes in autism.Nat. Commun.5:5748. 10.1038/ncomms6748
42
Ha S. Lee D. Cho Y. S. Chung C. Yoo Y. E. Kim J. et al (2016). Cerebellar Shank2 regulates excitatory synapse density, motor coordination, and specific repetitive and anxiety-like behaviors.J. Neurosci.3612129–12143. 10.1523/JNEUROSCI.1849-16.2016
43
Han K. A. Yoon T. H. Shin J. Um J. W. Ko J. (2020). Differentially altered social dominance- and cooperative-like behaviors in Shank2- and Shank3-mutant mice.Mol. Autism11:87. 10.1186/s13229-020-00392-9
44
Heise C. Preuss J. M. Schroeder J. C. Battaglia C. R. Kolibius J. Schmid R. et al (2018). Heterogeneity of cell surface glutamate and GABA receptor expression in shank and CNTN4 autism mouse models.Front. Mol. Neurosci.11:212. 10.3389/fnmol.2018.00212
45
Hoeffer C. A. Klann E. (2010). mTOR signaling: at the crossroads of plasticity, memory and disease.Trends Neurosci.3367–75. 10.1016/j.tins.2009.11.003
46
Hollis F. Kanellopoulos A. K. Bagni C. (2017). Mitochondrial dysfunction in autism spectrum disorder: clinical features and perspectives.Curr. Opin. Neurobiol.45178–187. 10.1016/j.conb.2017.05.018
47
Homann O. R. Misura K. Lamas E. Sandrock R. W. Nelson P. McDonough S. I. et al (2016). Whole-genome sequencing in multiplex families with psychoses reveals mutations in the SHANK2 and SMARCA1 genes segregating with illness.Mol. Psychiatry.211690–1695. 10.1038/mp.2016.24
48
Horner A. E. Norris R. H. McLaren-Jones R. Alexander L. Komiyama N. H. Grant S. G. N. et al (2021). Learning and reaction times in mouse touchscreen tests are differentially impacted by mutations in genes encoding postsynaptic interacting proteins SYNGAP1, NLGN3, DLGAP1, DLGAP2 and SHANK2.Genes Brain Behav.20:e12723. 10.1111/gbb.12723
49
Huang, da W. Sherman B. T. Lempicki R. A. (2009). Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources.Nat. Protoc.444–57. 10.1038/nprot.2008.211
50
Iossifov I. O’Roak B. J. Sanders S. J. Ronemus M. Krumm N. Levy D. et al (2014). The contribution of de novo coding mutations to autism spectrum disorder.Nature515216–221. 10.1038/nature13908
51
Jiang Y. H. Ehlers M. D. (2013). Modeling autism by SHANK gene mutations in mice.Neuron788–27. 10.1016/j.neuron.2013.03.016
52
Jung H. Park H. Choi Y. Kang H. Lee E. Kweon H. et al (2018). Sexually dimorphic behavior, neuronal activity, and gene expression in Chd8-mutant mice.Nat. Neurosci.211218–1228. 10.1038/s41593-018-0208-z
53
Kang H. J. Kawasawa Y. I. Cheng F. Zhu Y. Xu X. Li M. et al (2011). Spatio-temporal transcriptome of the human brain.Nature478483–489. 10.1038/nature10523
54
Kawamura A. Abe Y. Seki F. Katayama Y. Nishiyama M. Takata N. et al (2020a). Chd8 mutation in oligodendrocytes alters microstructure and functional connectivity in the mouse brain.Mol. Brain13:160. 10.1186/s13041-020-00699-x
55
Kawamura A. Katayama Y. Nishiyama M. Shoji H. Tokuoka K. Ueta Y. et al (2020b). Oligodendrocyte dysfunction due to Chd8 mutation gives rise to behavioral deficits in mice.Hum. Mol. Genet.291274–1291. 10.1093/hmg/ddaa036
56
Khemakhem A. M. Frye R. E. El-Ansary A. Al-Ayadhi L. Bacha A. B. (2017). Novel biomarkers of metabolic dysfunction is autism spectrum disorder: potential for biological diagnostic markers.Metab. Brain Dis.321983–1997. 10.1007/s11011-017-0085-2
57
Kim E. Sheng M. (2004). PDZ domain proteins of synapses.Nat. Rev. Neurosci.5771–781. 10.1038/nrn1517
58
Kim R. Kim J. Chung C. Ha S. Lee S. Lee E. et al (2018). Cell-type-specific Shank2 deletion in mice leads to differential synaptic and behavioral phenotypes.J. Neurosci.384076–4092. 10.1523/JNEUROSCI.2684-17.2018
59
Ko H. G. Oh S. B. Zhuo M. Kaang B. K. (2016). Reduced acute nociception and chronic pain in Shank2-/- mice.Mol Pain12:1744806916647056. 10.1177/1744806916647056
60
Lazaro M. T. Taxidis J. Shuman T. Bachmutsky I. Ikrar T. Santos R. et al (2019). Reduced prefrontal synaptic connectivity and disturbed oscillatory population dynamics in the CNTNAP2 model of autism.Cell Rep.272567–2578e6. 10.1016/j.celrep.2019.05.006
61
Leblond C. S. Heinrich J. Delorme R. Proepper C. Betancur C. Huguet G. et al (2012). Genetic and functional analyses of SHANK2 mutations suggest a multiple hit model of autism spectrum disorders.PLoS Genet.8:e1002521.
62
Leblond C. S. Nava C. Polge A. Gauthier J. Huguet G. Lumbroso S. et al (2014). Meta-analysis of SHANK mutations in autism spectrum disorders: a gradient of severity in cognitive impairments.PLoS Genet.10:e1004580. 10.1371/journal.pgen.1004580
63
Lee A. Hirabayashi Y. Kwon S. K. Lewis T. L. Jr. Polleux F. (2018). Emerging roles of mitochondria in synaptic transmission and neurodegeneration.Curr. Opin. Physiol.382–93. 10.1016/j.cophys.2018.03.009
64
Lee E. Rhim I. Lee J. W. Ghim J. W. Lee S. Kim E. et al (2016). Enhanced neuronal activity in the medial prefrontal cortex during social approach behavior.J. Neurosci.366926–6936. 10.1523/JNEUROSCI.0307-16.2016
65
Lee E. J. Lee H. Huang T. N. Chung C. Shin W. Kim K. et al (2015). Trans-synaptic zinc mobilization improves social interaction in two mouse models of autism through NMDAR activation.Nat. Commun.6:7168. 10.1038/ncomms8168
66
Lee S. Lee E. Kim R. Kim J. Lee S. Park H. et al (2018). Shank2 deletion in parvalbumin neurons leads to moderate hyperactivity, enhanced self-grooming and suppressed seizure susceptibility in mice.Front. Mol. Neurosci.11:209. 10.3389/fnmol.2018.00209
67
Lee Y. S. Yu N. K. Chun J. Yang J. E. Lim C. S. Kim H. et al (2020). Identification of a novel Shank2 transcriptional variant in Shank2 knockout mouse model of autism spectrum disorder.Mol. Brain13:54. 10.1186/s13041-020-00595-4
68
Legido A. Jethva R. Goldenthal M. J. (2013). Mitochondrial dysfunction in autism.Semin. Pediatr. Neurol.20163–175. 10.1016/j.spen.2013.10.008
69
Levy D. R. Tamir T. Kaufman M. Parabucki A. Weissbrod A. Schneidman E. et al (2019). Dynamics of social representation in the mouse prefrontal cortex.Nat. Neurosci.222013–2022. 10.1038/s41593-019-0531-z
70
Li S. Xiong G. J. Huang N. Sheng Z. H. (2020). The cross-talk of energy sensing and mitochondrial anchoring sustains synaptic efficacy by maintaining presynaptic metabolism.Nat. Metab.21077–1095. 10.1038/s42255-020-00289-0
71
Li Z. Okamoto K. Hayashi Y. Sheng M. (2004). The importance of dendritic mitochondria in the morphogenesis and plasticity of spines and synapses.Cell119873–887. 10.1016/j.cell.2004.11.003
72
Lim C. S. Kim H. Yu N. K. Kang S. J. Kim T. Ko H. G. et al (2017). Enhancing inhibitory synaptic function reverses spatial memory deficits in Shank2 mutant mice.Neuropharmacology112104–112. 10.1016/j.neuropharm.2016.08.016
73
Lim S. Naisbitt S. Yoon J. Hwang J. I. Suh P. G. Sheng M. et al (1999). Characterization of the shank family of synaptic proteins. Multiple genes, alternative splicing, and differential expression in brain and development.J. Biol. Chem.27429510–29518. 10.1074/jbc.274.41.29510
74
Liu Y. Du Y. Liu W. Yang C. Liu Y. Wang H. et al (2013). Lack of association between NLGN3, NLGN4, SHANK2 and SHANK3 gene variants and autism spectrum disorder in a Chinese population.PLoS One8:e56639. 10.1371/journal.pone.0056639
75
Lombard J. (1998). Autism: a mitochondrial disorder?Med. Hypothes.50497–500. 10.1016/S0306-9877(98)90270-5
76
Love M. I. Huber W. Anders S. (2014). Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2.Genome Biol.15:550. 10.1186/s13059-014-0550-8
77
Lu Z. A. Mu W. Osborne L. M. Cordner Z. A. (2018). Eighteen-year-old man with autism, obsessive compulsive disorder and a SHANK2 variant presents with severe anorexia that responds to high-dose fluoxetine.BMJ Case Rep.2018:bcr2018225119. 10.1136/bcr-2018-225119
78
Ma K. Qin L. Matas E. Duffney L. J. Liu A. Yan Z. (2018). Histone deacetylase inhibitor MS-275 restores social and synaptic function in a Shank3-deficient mouse model of autism.Neuropsychopharmacology431779–1788. 10.1038/s41386-018-0073-1
79
Merico D. Isserlin R. Stueker O. Emili A. Bader G. D. (2010). Enrichment map: a network-based method for gene-set enrichment visualization and interpretation.PLoS One5:e13984. 10.1371/journal.pone.0013984
80
Monteiro P. Feng G. (2017). SHANK proteins: roles at the synapse and in autism spectrum disorder.Nat. Rev. Neurosci.18147–157. 10.1038/nrn.2016.183
81
Mootha V. K. Lindgren C. M. Eriksson K. F. Subramanian A. Sihag S. Lehar J. et al (2003). PGC-1alpha-responsive genes involved in oxidative phosphorylation are coordinately downregulated in human diabetes.Nat. Genet.34267–273. 10.1038/ng1180
82
Morita M. Gravel S. P. Hulea L. Larsson O. Pollak M. St-Pierre J. et al (2015). mTOR coordinates protein synthesis, mitochondrial activity and proliferation.Cell Cycle14473–480. 10.4161/15384101.2014.991572
83
Mossa A. Giona F. Pagano J. Sala C. Verpelli C. (2017). SHANK genes in autism: defining therapeutic targets.Prog. Neuropsychopharmacol. Biol. Psychiatry84416–423. 10.1016/j.pnpbp.2017.11.019
84
Naisbitt S. Kim E. Tu J. C. Xiao B. Sala C. Valtschanoff J. et al (1999). Shank, a novel family of postsynaptic density proteins that binds to the NMDA receptor/PSD-95/GKAP complex and cortactin.Neuron23569–582. 10.1016/S0896-6273(00)80809-0
85
Oliveira G. Diogo L. Grazina M. Garcia P. Ataide A. Marques C. et al (2005). Mitochondrial dysfunction in autism spectrum disorders: a population-based study.Dev. Med. Child Neurol.47185–189. 10.1017/S0012162205000332
86
Pappas A. L. Bey A. L. Wang X. Rossi M. Kim Y. H. Yan H. et al (2017). Deficiency of Shank2 causes mania-like behavior that responds to mood stabilizers.JCI Insight2:e92052. 10.1172/jci.insight.92052
87
Parikshak N. N. Swarup V. Belgard T. G. Irimia M. Ramaswami G. Gandal M. J. et al (2016). Genome-wide changes in lncRNA, splicing, and regional gene expression patterns in autism.Nature540423–427. 10.1038/nature20612
88
Patro R. Duggal G. Love M. I. Irizarry R. A. Kingsford C. (2017). Salmon provides fast and bias-aware quantification of transcript expression.Nat. Methods14417–419. 10.1038/nmeth.4197
89
Peter S. Ten Brinke M. M. Stedehouder J. Reinelt C. M. Wu B. Zhou H. et al (2016). Dysfunctional cerebellar Purkinje cells contribute to autism-like behaviour in Shank2-deficient mice.Nat. Commun.7:12627. 10.1038/ncomms12627
90
Peykov S. Berkel S. Degenhardt F. Rietschel M. Nothen M. M. Rappold G. A. (2015a). Rare SHANK2 variants in schizophrenia.Mol. Psychiatry201487–1488. 10.1038/mp.2015.122
91
Peykov S. Berkel S. Schoen M. Weiss K. Degenhardt F. Strohmaier J. et al (2015b). Identification and functional characterization of rare SHANK2 variants in schizophrenia.Mol. Psychiatry201489–1498. 10.1038/mp.2014.172
92
Pinero J. Ramirez-Anguita J. M. Sauch-Pitarch J. Ronzano F. Centeno E. Sanz F. et al (2020). The DisGeNET knowledge platform for disease genomics: 2019 update.Nucleic Acids Res.48D845–D855. 10.1093/nar/gkz1021
93
Pinto D. Pagnamenta A. T. Klei L. Anney R. Merico D. Regan R. et al (2010). Functional impact of global rare copy number variation in autism spectrum disorders.Nature466368–372. 10.1038/nature09146
94
Prasad A. Merico D. Thiruvahindrapuram B. Wei J. Lionel A. C. Sato D. et al (2012). A discovery resource of rare copy number variations in individuals with autism spectrum disorder.G321665–1685. 10.1534/g3.112.004689
95
Qin L. Ma K. Wang Z. J. Hu Z. Matas E. Wei J. et al (2018). Social deficits in Shank3-deficient mouse models of autism are rescued by histone deacetylase (HDAC) inhibition.Nat. Neurosci.21564–575. 10.1038/s41593-018-0110-8
96
Rauch A. Wieczorek D. Graf E. Wieland T. Endele S. Schwarzmayr T. et al (2012). Range of genetic mutations associated with severe non-syndromic sporadic intellectual disability: an exome sequencing study.Lancet3801674–1682. 10.1016/S0140-6736(12)61480-9
97
Rorbach J. Gammage P. A. Minczuk M. (2012). C7orf30 is necessary for biogenesis of the large subunit of the mitochondrial ribosome.Nucleic Acids Res.404097–4109. 10.1093/nar/gkr1282
98
Rose S. Niyazov D. M. Rossignol D. A. Goldenthal M. Kahler S. G. Frye R. E. (2018). Clinical and molecular characteristics of mitochondrial dysfunction in autism spectrum disorder.Mol. Diagn. Ther.22571–593. 10.1007/s40291-018-0352-x
99
Rossi M. J. Pekkurnaz G. (2019). Powerhouse of the mind: mitochondrial plasticity at the synapse.Curr. Opin. Neurobiol.57149–155. 10.1016/j.conb.2019.02.001
100
Rossignol D. A. Frye R. E. (2012). Mitochondrial dysfunction in autism spectrum disorders: a systematic review and meta-analysis.Mol. Psychiatry17290–314. 10.1038/mp.2010.136
101
Sala C. Vicidomini C. Bigi I. Mossa A. Verpelli C. (2015). Shank synaptic scaffold proteins: keys to understanding the pathogenesis of autism and other synaptic disorders.J. Neurochem.135849–858. 10.1111/jnc.13232
102
Sanders S. J. Murtha M. T. Gupta A. R. Murdoch J. D. Raubeson M. J. Willsey A. J. et al (2012). De novo mutations revealed by whole-exome sequencing are strongly associated with autism.Nature485237–241. 10.1038/nature10945
103
Santini E. Klann E. (2014). Reciprocal signaling between translational control pathways and synaptic proteins in autism spectrum disorders.Sci. Signal.7:re10. 10.1126/scisignal.2005832
104
Sato M. Kawano M. Mizuta K. Islam T. Lee M. G. Hayashi Y. (2017). Hippocampus-dependent goal localization by head-fixed mice in virtual reality.eNeuro4:ENEURO.0369-16.2017. 10.1523/ENEURO.0369-16.2017
105
Sato M. Mizuta K. Islam T. Kawano M. Sekine Y. Takekawa T. et al (2020). Distinct mechanisms of over-representation of landmarks and rewards in the hippocampus.Cell Rep.32:107864. 10.1016/j.celrep.2020.107864
106
Satterstrom F. K. Kosmicki J. A. Wang J. Breen M. S. De Rubeis S. An J. Y. et al (2020). Large-scale exome sequencing study implicates both developmental and functional changes in the neurobiology of autism.Cell180568.e23–584.e23. 10.1016/j.cell.2019.12.036
107
Saxton R. A. Sabatini D. M. (2017). mTOR signaling in growth, metabolism, and disease.Cell169361–371. 10.1016/j.cell.2017.03.035
108
Schluth-Bolard C. Labalme A. Cordier M. P. Till M. Nadeau G. Tevissen H. et al (2013). Breakpoint mapping by next generation sequencing reveals causative gene disruption in patients carrying apparently balanced chromosome rearrangements with intellectual deficiency and/or congenital malformations.J. Med. Genet.50144–150. 10.1136/jmedgenet-2012-101351
109
Schmeisser M. J. (2015). Translational neurobiology in Shank mutant mice–model systems for neuropsychiatric disorders.Ann. Anat.200115–117. 10.1016/j.aanat.2015.03.006
110
Schmeisser M. J. Ey E. Wegener S. Bockmann J. Stempel V. Kuebler A. et al (2012). Autistic-like behaviours and hyperactivity in mice lacking ProSAP1/Shank2.Nature486256–260. 10.1038/nature11015
111
Selimbeyoglu A. Kim C. K. Inoue M. Lee S. Y. Hong A. S. O. Kauvar I. et al (2017). Modulation of prefrontal cortex excitation/inhibition balance rescues social behavior in CNTNAP2-deficient mice.Sci. Transl. Med.9:eaah6733. 10.1126/scitranslmed.aah6733
112
Sheng M. Kim E. (2000). The Shank family of scaffold proteins.J. Cell Sci.113(Pt 11), 1851–1856. 10.1242/jcs.113.11.1851
113
Sheng M. Kim E. (2011). The postsynaptic organization of synapses.Cold Spring Harbor perspectives in biology3:a005678. 10.1101/cshperspect.a005678
114
Sheng Z. H. Cai Q. (2012). Mitochondrial transport in neurons: impact on synaptic homeostasis and neurodegeneration.Nat. Rev. Neurosci.1377–93. 10.1038/nrn3156
115
Siddiqui M. F. Elwell C. Johnson M. H. (2016). Mitochondrial dysfunction in autism spectrum disorders.Autism Open Access61000190. 10.4172/2165-7890.1000190
116
Singh K. Singh I. N. Diggins E. Connors S. L. Karim M. A. Lee D. et al (2020). Developmental regression and mitochondrial function in children with autism.Ann. Clin. Transl. Neurol.7683–694. 10.1002/acn3.51034
117
Smits P. Smeitink J. van den Heuvel L. (2010). Mitochondrial translation and beyond: processes implicated in combined oxidative phosphorylation deficiencies.J. Biomed. Biotechnol.2010:737385. 10.1155/2010/737385
118
Soneson C. Love M. I. Robinson M. D. (2015). Differential analyses for RNA-seq: transcript-level estimates improve gene-level inferences.F1000Res4:1521. 10.12688/f1000research.7563.1
119
Sossin W. S. Costa-Mattioli M. (2019). Translational control in the brain in health and disease.Cold Spring Harb. Perspect. Biol.11:a032912. 10.1101/cshperspect.a032912
120
Subramanian A. Tamayo P. Mootha V. K. Mukherjee S. Ebert B. L. Gillette M. A. et al (2005). Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles.Proc. Natl. Acad. Sci. U.S.A.10215545–15550. 10.1073/pnas.0506580102
121
Switon K. Kotulska K. Janusz-Kaminska A. Zmorzynska J. Jaworski J. (2017). Molecular neurobiology of mTOR.Neuroscience341112–153. 10.1016/j.neuroscience.2016.11.017
122
Unsicker C. Cristian F. B. von Hahn M. Eckstein V. Rappold G. A. Berkel S. (2021). SHANK2 mutations impair apoptosis, proliferation and neurite outgrowth during early neuronal differentiation in SH-SY5Y cells.Sci. Rep.11:2128. 10.1038/s41598-021-81241-4
123
Velmeshev D. Schirmer L. Jung D. Haeussler M. Perez Y. Mayer S. et al (2019). Single-cell genomics identifies cell type-specific molecular changes in autism.Science364685–689. 10.1126/science.aav8130
124
Voineagu I. Wang X. Johnston P. Lowe J. K. Tian Y. Horvath S. et al (2011). Transcriptomic analysis of autistic brain reveals convergent molecular pathology.Nature474380–384. 10.1038/nature10110
125
Vos M. Lauwers E. Verstreken P. (2010). Synaptic mitochondria in synaptic transmission and organization of vesicle pools in health and disease.Front. Synapt. Neurosci.2:139. 10.3389/fnsyn.2010.00139
126
Wang F. Zhu J. Zhu H. Zhang Q. Lin Z. Hu H. (2011). Bidirectional control of social hierarchy by synaptic efficacy in medial prefrontal cortex.Science334693–697. 10.1126/science.1209951
127
Wang T. Hoekzema K. Vecchio D. Wu H. Sulovari A. Coe B. P. et al (2020). Large-scale targeted sequencing identifies risk genes for neurodevelopmental disorders.Nat. Commun.11:4932.
128
Wang Z. J. Zhong P. Ma K. Seo J. S. Yang F. Hu Z. et al (2019). Amelioration of autism-like social deficits by targeting histone methyltransferases EHMT1/2 in Shank3-deficient mice.Mol. Psychiatry252517–2533. 10.1038/s41380-019-0351-2
129
Wanschers B. F. Szklarczyk R. Pajak A. van den Brand M. A. Gloerich J. Rodenburg R. J. et al (2012). C7orf30 specifically associates with the large subunit of the mitochondrial ribosome and is involved in translation.Nucleic Acids Res.404040–4051. 10.1093/nar/gkr1271
130
Wegener S. Buschler A. Stempel A. V. Kang S. J. Lim C. S. Kaang B. K. et al (2018). Defective synapse maturation and enhanced synaptic plasticity in Shank2 Deltaex7(-/-) mice.eNeuro5:ENEURO.0398-17.2018. 10.1523/ENEURO.0398-17.2018
131
Weissman J. R. Kelley R. I. Bauman M. L. Cohen B. H. Murray K. F. Mitchell R. L. et al (2008). Mitochondrial disease in autism spectrum disorder patients: a cohort analysis.PLoS One3:e3815. 10.1371/journal.pone.0003815
132
Werling D. M. Parikshak N. N. Geschwind D. H. (2016). Gene expression in human brain implicates sexually dimorphic pathways in autism spectrum disorders.Nat. Commun.7:10717. 10.1038/ncomms10717
133
Winden K. D. Ebrahimi-Fakhari D. Sahin M. (2018). Abnormal mTOR activation in autism.Annu. Rev. Neurosci.411–23. 10.1146/annurev-neuro-080317-061747
134
Wischmeijer A. Magini P. Giorda R. Gnoli M. Ciccone R. Cecconi L. et al (2011). Olfactory receptor-related duplicons mediate a microdeletion at 11q13.2q13.4 associated with a syndromic phenotype.Mol. Syndromol.1176–184. 10.1159/000322054
135
Won H. Lee H. R. Gee H. Y. Mah W. Kim J. I. Lee J. et al (2012). Autistic-like social behaviour in Shank2-mutant mice improved by restoring NMDA receptor function.Nature486261–265. 10.1038/nature11208
136
Xu L. M. Li J. R. Huang Y. Zhao M. Tang X. Wei L. (2012). AutismKB: an evidence-based knowledgebase of autism genetics.Nucleic Acids Res.40D1016–D1022. 10.1093/nar/gkr1145
137
Yang C. Li J. Wu Q. Yang X. Huang A. Y. Zhang J. et al (2018). AutismKB 2.0: a knowledgebase for the genetic evidence of autism spectrum disorder.Database2018:bay106. 10.1093/database/bay106
138
Yizhar O. Fenno L. E. Prigge M. Schneider F. Davidson T. J. O’Shea D. J. et al (2011). Neocortical excitation/inhibition balance in information processing and social dysfunction.Nature477171–178. 10.1038/nature10360
139
Yoo J. Bakes J. Bradley C. Collingridge G. L. Kaang B. K. (2014). Shank mutant mice as an animal model of autism.Philos. Trans. R. Soc. Lond. Ser. B Biol. Sci.369:20130143. 10.1098/rstb.2013.0143
140
Yook C. Kim K. Kim D. Kang H. Kim S. G. Kim E. et al (2019). A TBR1-K228E mutation induces Tbr1 upregulation, altered cortical distribution of interneurons, increased inhibitory synaptic transmission, and autistic-like behavioral deficits in mice.Front. Mol. Neurosci.12:241. 10.3389/fnmol.2019.00241
141
Yoon S. Y. Kwon S. G. Kim Y. H. Yeo J. H. Ko H. G. Roh D. H. et al (2017). A critical role of spinal Shank2 proteins in NMDA-induced pain hypersensitivity.Mol. Pain13:1744806916688902. 10.1177/1744806916688902
142
Yuen R. K. C. Merico D. Bookman M. Thiruvahindrapuram B. Patel R. V. Whitney J. et al (2017). Whole genome sequencing resource identifies 18 new candidate genes for autism spectrum disorder.Nat. Neurosci.20602–611. 10.1038/nn.4524
143
Zaslavsky K. Zhang W. B. McCready F. P. Rodrigues D. C. Deneault E. Loo C. et al (2019). SHANK2 mutations associated with autism spectrum disorder cause hyperconnectivity of human neurons.Nat. Neurosci.22556–564. 10.1038/s41593-019-0365-8
144
Zhang F. Rein B. Zhong P. Shwani T. Conrow-Graham M. Wang Z. J. et al (2021). Synergistic inhibition of histone modifiers produces therapeutic effects in adult Shank3-deficient mice.Transl. Psychiatry11:99. 10.1038/s41398-021-01233-w
145
Zhou W. Z. Zhang J. Li Z. Lin X. Li J. Wang S. et al (2019). Targeted resequencing of 358 candidate genes for autism spectrum disorder in a Chinese cohort reveals diagnostic potential and genotype-phenotype correlations.Hum. Mutat.40801–815. 10.1002/humu.23724
146
Zhu L. Wang X. Li X. L. Towers A. Cao X. Wang P. et al (2014). Epigenetic dysregulation of SHANK3 in brain tissues from individuals with autism spectrum disorders.Hum. Mol. Genet.231563–1578. 10.1093/hmg/ddt547
Summary
Keywords
autism spectrum disorder, Shank2, prefrontal cortex, transcript, RNA-seq, dosage-dependence, age-dependence
Citation
Lee S, Kang H, Jung H, Kim E and Lee E (2021) Gene Dosage- and Age-Dependent Differential Transcriptomic Changes in the Prefrontal Cortex of Shank2-Mutant Mice. Front. Mol. Neurosci. 14:683196. doi: 10.3389/fnmol.2021.683196
Received
20 March 2021
Accepted
18 May 2021
Published
11 June 2021
Volume
14 - 2021
Edited by
Andreas Martin Grabrucker, University of Limerick, Ireland
Reviewed by
Simone Berkel, University Hospital Heidelberg, Germany; Rolf Sprengel, Max Planck Institute for Medical Research (MPIMF), Germany; Zhanyan Fu, Broad Institute, United States
Updates
Copyright
© 2021 Lee, Kang, Jung, Kim and Lee.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Eunjoon Kim, kime@kaist.ac.krEunee Lee, eunee@yuhs.ac
This article was submitted to Neuroplasticity and Development, a section of the journal Frontiers in Molecular Neuroscience
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.