- Open Access
A direct regulatory link between microRNA-137 and SHANK2: implications for neuropsychiatric disorders
Journal of Neurodevelopmental Disorders volume 10, Article number: 15 (2018)
Mutations in the SHANK genes, which encode postsynaptic scaffolding proteins, have been linked to a spectrum of neurodevelopmental disorders. The SHANK genes and the schizophrenia-associated microRNA-137 show convergence on several levels, as they are both expressed at the synapse, influence neuronal development, and have a strong link to neurodevelopmental and neuropsychiatric disorders like intellectual disability, autism, and schizophrenia. This compiled evidence raised the question if the SHANKs might be targets of miR-137.
In silico analysis revealed a putative binding site for microRNA-137 (miR-137) in the SHANK2 3′UTR, while this was not the case for SHANK1 and SHANK3. Luciferase reporter assays were performed by overexpressing wild type and mutated SHANK2-3′UTR and miR-137 in human neuroblastoma cells and mouse primary hippocampal neurons. miR-137 was also overexpressed or inhibited in hippocampal neurons, and Shank2 expression was analyzed by quantitative real-time PCR and Western blot. Additionally, expression levels of experimentally validated miR-137 target genes were analyzed in the dorsolateral prefrontal cortex (DLPFC) of schizophrenia and control individuals using the RNA-Seq data from the CommonMind Consortium.
miR-137 directly targets the 3′UTR of SHANK2 in a site-specific manner. Overexpression of miR-137 in mouse primary hippocampal neurons significantly lowered endogenous Shank2 protein levels without detectable influence on mRNA levels. Conversely, miR-137 inhibition increased Shank2 protein expression, indicating that miR-137 regulates SHANK2 expression by repressing protein translation rather than inducing mRNA degradation.
To find out if the miR-137 signaling network is altered in schizophrenia, we compared miR-137 precursor and miR-137 target gene expression in the DLPFC of schizophrenia and control individuals using the CommonMind Consortium RNA sequencing data. Differential expression of 23% (16/69) of known miR-137 target genes was detected in the DLPFC of schizophrenia individuals compared with controls. We propose that in further targets (e.g., SHANK2, as described in this paper) which are not regulated on RNA level, effects may only be detectable on protein level.
Our study provides evidence that a direct regulatory link exists between miR-137 and SHANK2 and supports the finding that miR-137 signaling might be altered in schizophrenia.
microRNAs are small non-coding, single-stranded RNA molecules that regulate gene expression by binding to the 3′UTRs of their target mRNAs through base pairing of their 6–8-nucleotide long seed region. They either initiate mRNA degradation or inhibit protein translation [1, 2]. Regulation of gene expression in the nervous system is highly complex and is fine-tuned by microRNAs . In neurons, mRNAs and microRNAs are compartmentalized to specific subcellular regions, like the synaptodendritic compartment, to regulate local protein synthesis in response to neuronal stimuli .
microRNA-137 (miR-137) is enriched in the human and mouse brain , especially in cortical regions and in the hippocampus . It plays a role in cell proliferation and differentiation and is present at the synapse [6, 7]. Furthermore, it regulates neuronal maturation, dendritic morphogenesis, and spine development . Homozygous knockout of Mir137 is embryonically lethal in mice, indicating that embryonic development is dependent on at least one functional allele . This emphasizes the importance of miR-137 in early developmental processes.
Genetic alterations of MIR137 have been associated with neurodevelopmental disorders. The rare chromosomal microdeletion 1p21.3 encompassing MIR137 and DPYD was identified in individuals with intellectual disability (ID), comorbid with autism spectrum disorder (ASD), and obesity [6, 10, 11]. In addition, 19/65 (29%) highest ranked ASD risk genes in the Simons Foundation Autism Research Initiative gene database [12, 13] (URL: https://sfari.org/resources/sfari-gene#refs, accessed May 2017) are predicted or confirmed miR-137 targets (Additional file 1: Table S1). miR-137 target genes are enriched among ASD risk genes (29%, CI 18.6–41.8%) compared with the random frequency (2.9%, predicted in the miRDB database ; CI 2.7–3.2%, 525/17860) (P ≤ 0.000001, Fisher’s exact test, two-sided). Furthermore, the MIR137 gene has been reported as a schizophrenia (SCZ) susceptibility locus because the common SNP rs1625579, located in an intron of MIR137, was associated with SCZ in several studies [15,16,17,18]. Putative miR-137 targets are enriched in SCZ risk loci , suggesting that convergent pathways connected by this microRNA contribute to SCZ etiology. A meta-analysis of SNP data across five different neuropsychiatric disorders identified a shared risk locus association with SCZ and ASD . Taken together, these findings indicate that deregulation of miR-137 is the key to various psychiatric and neurodevelopmental disorders.
SHANK proteins (SHANK1, SHANK2, and SHANK3) are postsynaptic scaffolding proteins with an important role in the formation, maintenance, and function of synapses in glutamatergic neurons of the brain . SHANK genes converge with miR-137 on several levels: (i) both are expressed at the synapse, (ii) both influence dendrite and spine formation in glutamatergic neurons [8, 21], and (iii) both have a strong link to neurodevelopmental and neuropsychiatric disorders like ID, ASD, and SCZ [18, 22, 23]. This suggests that SHANK genes might be targets of miR-137.
To validate a predicted miR-137 binding site in SHANK2, luciferase reporter assays were carried out in human SH-SY5Y cells and mouse primary hippocampal neurons. The wild type SHANK2-3′UTR was cloned into the psiCHECK™-2 vector (Promega, Mannheim, Germany) after PCR amplification on human genomic DNA from a healthy individual. PCR primers introduced the restriction sites (XhoI) which were used to introduce the PCR product into the psiCHECK™-2 vector via Gibson® Assembly (New England BioLabs, Frankfurt, Germany). Mutagenesis of the miR-137 binding site was performed on the wild type construct by PCR with the original cloning primers and mutagenesis primers, assembled using the SLiCE cloning system . The sequence of wild type and mutated SHANK2-3′UTR was validated by Sanger sequencing. Primer sequences are given in Additional file 1: Table S2.
The SH-SY5Y cells were seeded in antibiotic-free medium and transfected after 24 h with 200 ng of empty psiCHECK™-2 vector and wild type or mutated SHANK2 3′UTR constructs with 12 nM miR-137 or control mimics (mirVana® miRNA hsa-miR-137 mimic MIMAT0000429 and Pre-miR™ miRNA Precursor Molecules Negative Control #2, AM17111, Thermo Fisher Scientific, Darmstadt, Germany) in technical triplicates using Lipofectamine®2000 (Thermo Fisher Scientific, Darmstadt, Germany). miRNA mimics are small, double-stranded RNA molecules which mimic endogenous, mature microRNA molecules when introduced into the cells. Luciferase activity was measured 48 h posttransfection using the dual luciferase reporter assay system (Promega) and the luminometer Centro LB960 (Berthold Technologies, Bad Wildbad, Germany). Renilla luciferase activity was normalized against firefly luciferase activity, and the background was subtracted by normalization to the empty psiCHECK™-2 vector activity.
Primary hippocampal neurons were transfected on day in vitro (DIV) 5 using the same experimental setup as in the SH-SY5Y cells for the mimic experiments with 25 nM mimics. The amount of miR-137 mimics had to be optimized for each cell type separately. For the inhibitor experiments, the primary mouse hippocampal neurons were treated with 25 nM miR-137 or power inhibitor control (miRCURY LNA™ Power microRNA hsa-miR-137 4101446-111 or negative control B 199007-111 inhibitor 5′-Fluorescein-labeled, Exiqon, Vedbaek, Denmark) added drop-wise to the cells on DIV5. The miR-137 power inhibitor sequesters both the human hsa-miR-137 and mouse mmu-miR-137, thereby preventing it from binding to its targets. The control inhibitor does not influence miR-137 and was used to control for unspecific effects of the treatment. Experimental conditions were identical between both groups regarding cell number and cell types (Additional file 1: Figure S1). The cells were transfected with the luciferase reporter constructs with Lipofectamine®2000 on DIV6. Luciferase activity was measured 48 h posttransfection as described above.
The SH-SY5Y cells were obtained from the DSMZ (Leibniz Institute German collection of Microorganisms and Cell Cultures, Braunschweig, Germany, ACC209). The cells were kept in supplemented Dulbecco’s Modified Eagle’s Medium (DMEM) (15% fetal bovine serum (FBS), 1% penicillin/streptomycin (P/S), 1% non-essential amino acids). Hippocampal cultures were prepared from gestating CD-1® IGS mice (Charles River, Sulzfeld, Germany) at embryonic stage 15.5. The cells were seeded in plating media (DMEM + 10% FBS, 0.25% P/S, 2 mM L-glutamine) on poly-L-ornithine (Sigma Aldrich, Schnelldorf, Germany) coating. On DIV 1, a complete medium change to supplemented Neurobasal medium (Gibco® Neurobasal® media + 2% Gibco®B-27® supplement, 0.25% P/S, 0.5 mM L-glutamine) was performed. On DIV3, a half medium change was carried out with 5 μM cytosine β-D-arabinofuranoside in supplemented Neurobasal media, except for the cells that were harvested on DIV3. From DIV5 until the final culture day, a half medium change was performed with supplemented Neurobasal media every other day, except for the cells that were transfected or harvested.
Nucleofection of mouse hippocampal cells
Hippocampal cultures were transfected with the optimized hsa-miR-137 mimic concentration of 300 nM and negative miRNA control on DIV5 using the AD1 4D-Nucleofector™ Y Unit system (Lonza, Basel, Switzerland) according to manufacturer’s instructions. Cell types and cell numbers showed no difference after nucleofection between the two conditions (Additional file 1: Figure S1).
To determine the effect of miR-137 on endogenous Shank2 expression real-time quantitative PCR (RT-qPCR) was performed. For the quantification of Shank2 expression, RNA was isolated from mouse primary hippocampal neurons (DIV8) 72 h after transfection with using the Quick-RNA™ MicroPrep Kit (Zymo Research, Freiburg, Germany) and was transcribed with the SuperScript® VILO™ cDNA Synthesis Kit (Thermo Fisher Scientific, Darmstadt, Germany). RT-qPCR was performed to quantify Shank2 expression normalized against two reference genes (Gapdh and Hprt1). Primer sequences are provided in Additional file 1: Table S2.
The RNAqueous®-Micro Total RNA Isolation Kit (Thermo Fisher Scientific) was used for total RNA isolation. Ten nanograms of total RNA were transcribed using the miRCURY LNA™ microRNA PCR Universal cDNA synthesis Kit (Exiqon) with a spike-in loading control according to manufacturer’s instructions. RT-qPCR was performed with miRCURY LNA™ microRNA PCR ExiLENT SYBR® Green PCR sets for hsa-miR-137, U6, and spike-in (Exiqon), measuring the endogenous miR-137 expression profile in primary hippocampal neurons over a time period of 11 DIV. miR-137 expression levels were analyzed with the same method using 21 different human tissue RNA samples (specified in Additional file 1: Table S3).
Western blot analysis
To determine the effect of miR-137 on Shank2 protein expression, cellular protein was isolated using RIPA buffer from cell cultures on DIV10, and the lysates were run on Novex™WedgeWell™4–12% Tris-Glycine Gels (Thermo Fisher Scientific) and then blotted onto a PVDF membrane (Immobilon-FL, Millipore, Billerica, Massachusetts, USA) as recommended by the manufacturer. Primary antibodies anti-beta-III-tubulin (mouse, G7121 Promega) and anti-Shank2 (rabbit, ABIN656710 antibodies-online) were used followed by the secondary IRDye 680 donkey anti-mouse and IRDye 800CW donkey anti-rabbit antibodies (1706515 and 1706516, respectively, Li-COR). The membrane was scanned with the Odyssey® Infrared Imaging System (Li-COR) and quantified using ImageJ software39. Shank2 was normalized to beta-III-tubulin. All experiments were performed sequentially over several weeks. Hippocampal neurons were isolated and differentiated at different time points.
The tissue used for RNA sequencing was obtained from the following institutes: Mount Sinai School of Medicine, University of Pennsylvania, University of Pittsburgh, National Institute of Mental Health, F. Hoffman-La Roche Ltd., Takeda Pharmaceuticals Company Limited, Sage Bionetworks, Duke University, and University of North Carolina. Generation and analysis of the RNA sequencing data has been previously described in detail by the CommonMind Consortium . To analyze the expression levels of the validated miR-137 target genes and miR-137 precursor, we selected the average expression levels, the log2-fold change and the P values for the respective genes published on https://synapse.org in the “Results Explorer” on the “CommonMind Consortium Knowledge Portal” (data accessed June 2017). A weighted linear regression analysis had been performed for each gene considering covariates that influence expression levels. For each gene, the SCZ disease status coefficient had been statistically tested for being nonzero, implying an estimated effect for SCZ, above and beyond any other effect from covariates. This test produced a t statistic and a corresponding P value. Details about the statistical analysis are provided in the CommonMind Consortium paper .
Selection of experimentally validated miR-137 targets and network analysis
Validated targets of miR-137 were collected from miRTarBase (http://mirtarbase.mbc.nctu.edu.tw/, accessed January 2017) and PubMed (https://www.ncbi.nlm.nih.gov/pubmed/, accessed January 2017) database searches. miRTarBase targets were included if they were categorized as “strong evidence,” i.e., confirmed by experimental methods of validation including luciferase reporter assay, RT-qPCR, and/or Western blot analysis. Targets from the PubMed database search were considered as validated if a direct link was established by luciferase reporter assay, RT-qPCR, and/or Western blot analysis. Duplicates between databases were removed. A total of 71 validated targets were obtained; most of them have been investigated in cancer cells, whereas only 12 targets have been validated in neurons (details are shown in Additional file 1: Table S4). RNA sequencing data for the miR-137 target genes HCRT, SLC6A3, and SNAI1 was not available in the CommonMind dataset. The target gene network was analyzed through the use of QIAGEN’s Ingenuity®Pathway Analysis (IPA®, QIAGEN Redwood City, www.qiagen.com/ingenuity).
Statistical analyses were performed using the SPSS software (IBM Corp.; IBM SPSS Statistics for Windows, Version 22.0.; Armonk, NY, USA). Two-way ANOVAs were performed considering the influence of the experiment and the respective condition. The data is shown as mean ± SEM. The validated miR-137 target genes (including SHANK2) were ranked according to their point-wise P values from the CommonMind analysis, and the Benjamini-Hochberg method was used to correct for multiple testing with a false discovery rate of 10% (Additional file 1: Table S5). A two-sided Χ2 test with Yates correction was used to compare the expression of validated target genes of miR-137 and five control microRNAs. The differentially expressed targets of the five control microRNAs were pooled together. Genes that were targets of more than one control microRNA were only counted once. Eight target genes that were not differentially expressed in SCZ and control individuals overlapped between miR137 and the five controls and were excluded from analysis.
miR-137 directly targets the SHANK2 3′UTR
In silico analysis (TargetScanHuman Release 7.1 ) of all three SHANK genes identified a single, highly conserved binding site for miR-137 (MIMAT0000429) in the 3′UTR of SHANK2 (ENST00000449833.2, Additional file 1: Figure S2A), suggesting that this microRNA has the potential to bind SHANK2 mRNA. The analysis of miR-137 expression in 21 different human tissue samples revealed high expression in the central nervous system and marginal expression in other tested organs (Additional file 1: Figure S3). miR-137 was highly expressed in the fetal brain, and expression remained strong in the adult hippocampus, thalamus, and striatum. These are brain regions where SHANK2 is also expressed [27, 28].
To validate the predicted miR-137 binding site in SHANK2, we carried out luciferase reporter assays. Human neuroblastoma cells (SH-SY5Y) were transfected with a dual luciferase reporter plasmid (psiCHECK™-2) containing either the wild type (wt) SHANK2 3′UTR sequence or a miR-137 binding site-mutated 3′UTR (mut) together with either miR-137 or negative control miRNA mimics (Fig. 1a). The SH-SY5Y cells co-transfected with SHANK2 wild type 3′UTR reporter and miR-137 showed a 38% decrease in luciferase activity compared with the negative control (Fig. 1b; ***P = 5.2 × 10−9, two-way ANOVA). In contrast, no difference was observed between miR-137 and negative control miRNA for the mutated 3′UTR. This indicated that miR-137 specifically binds the target site, as mutation of the seed region in the miR-137 binding site eliminated the robust downregulatory effect of miR-137.
Alterations in the microRNA machinery have been described in various cancers ; therefore, we further analyzed the regulatory effect of miR-137 on Shank2 expression in mouse primary hippocampal neurons. These neurons express Shank2 and Mir137, and the miR-137 binding site is conserved in the mouse Shank2 3′UTR (Additional file 1: Figure S2B). We co-transfected SHANK2 wild type 3′UTR reporter and miR-137 and observed a 28% decrease in luciferase activity compared with the negative control (*P = 0.022, two-way ANOVA). In contrast, no effect was seen with the mutated 3′UTR (Fig. 1c). We concluded that miR-137 directly and specifically targets the predicted binding site in the 3′UTR of SHANK2 in the SH-SY5Y cells and primary hippocampal neurons.
miR-137 regulates Shank2 protein levels
We overexpressed miR-137 and negative control miRNA mimics in mouse hippocampal neurons to determine the effect on endogenous Shank2 expression. miR-137 overexpression did not alter endogenous Shank2 mRNA levels in hippocampal neurons (Fig. 1d). However, Western blotting revealed a 29% reduction of Shank2 protein (***P = 0.00003, two-way ANOVA) after miR-137 overexpression in hippocampal neurons compared with the negative control miRNA mimics (Fig. 1e).
To determine whether endogenous miR-137 is able to regulate Shank2 expression, we carried out luciferase assays and Western blot experiments in mouse primary hippocampal neurons after inhibition of miR-137 with a targeted miRCURY LNA™ Power microRNA inhibitor or treatment with a microRNA inhibitor control. First, we analyzed the endogenous miR-137 expression profile in primary hippocampal neurons. Expression levels remained constant until 5 days in vitro (DIV5), increased six-fold between DIV5 and DIV8, and remained high until DIV11 (Fig. 2a). Next, we treated cells with a miR-137 inhibitor or control inhibitor on DIV5, shortly before miR-137 reaches its highest endogenous expression levels. We performed luciferase reporter assays on DIV6 using SHANK2 3′UTR luciferase constructs. miR-137 inhibition increased the relative luciferase activity for the SHANK2 wt reporter construct by 70% (*P = 0.044, two-way ANOVA) compared with the control (Fig. 2b). Endogenous miR-137 was sequestered by the inhibitor; therefore, the regulatory effect on the SHANK2 3′UTR was lost and luciferase activity was increased. This was not observed with the mutated SHANK2 3′UTR reporter construct. Finally, we analyzed the consequences of endogenous miR-137 inhibition on Shank2 protein levels by Western blot analysis. Inhibition of miR-137 significantly increased endogenous Shank2 protein expression by 24% compared with control inhibitor (*P = 0.016, two-way ANOVA) (Fig. 2c).
Expression of miR-137 precursor and validated miR-137 target genes in the DLPFC of SCZ and control individuals
Previous studies have pointed to cellular, neurochemical and functional abnormalities in the dorsolateral prefrontal cortex (DLPFC) of SCZ individuals, implicating a role of this particular brain region in the pathophysiology of SCZ . In the DLPFC, SHANK2 and miR-137 are both expressed [27, 31] and we aimed to identify whether miR-137, SHANK2, and other known miR-137 targets are differentially expressed in the DLPFC of SCZ individuals compared to controls. We used a gene expression resource that provides RNA sequencing data in the DLPFC of deceased individuals with SCZ (n = 258) and control (n = 279) subjects (CommonMind Consortium) . Out of 16,423 genes analyzed, 693 genes were differentially expressed between SCZ and control individuals; 332 were upregulated and 361 were downregulated . The data for miR-137 precursor but not mature miR-137 expression levels were available in this data set. We examined the data for miR-137 precursor expression and found that it did not differ between SCZ and control individuals (MIR137HG, ENSG00000231269, P = 0.699). Next, we focused on the miR-137 signaling network and analyzed the expression of 69 known miR-137 target genes (68 validated targets from previously published studies summarized in Additional file 1: Table S4 and SHANK2). In total, 16 out of 69 genes were differentially expressed using the Benjamini-Hochberg method to correct for multiple testing with a false discovery rate of 10% (Table 1, Additional file 1: Table S5). This analysis was based on differential expression at the gene level. Four target genes (4/16; RORA, CPLX1, TCF4, SIRT1) showed significant expression differences on the whole-transcriptome level . All study-wide significant miR-137 target genes had modest fold changes with a mean of 1.09 and a range of 1.049–1.135 (inverting downregulated expression ratios). The majority of differentially expressed miR-137 target genes (12/16) showed elevated mRNA levels in SCZ individuals compared with controls. Five of the 12 upregulated genes (TCF4, SIRT1, XIAP, GRIA1, ZNF804A) had already been associated with SCZ. In addition, genes linked to ID (TCF4, CPLX1, CDK6, KDM5B), ASD (RORA, KDM5B), and bipolar disorder (ZNF804A) were also among the 16 differentially expressed target genes. Elevated mRNA expression of the SCZ-associated gene GRIA1 has been described in layer II/III and V pyramidal cells of the DLPFC in SCZ individuals .
In total, 23% of miR-137 targets were differentially expressed between SCZ and control individuals. We compared the expression of these targets to other microRNA targets. We selected the following five microRNAs that are expressed in the DLPFC but are not differentially expressed in SCZ individuals: let-7a, miR-21-5p, miR-93-5p, miR-451a, and miR-675-5p. Expression data for these microRNA precursors was not available in the CommonMind data set. We analyzed the number of validated target genes that were differentially expressed in SCZ individuals and revealed a lower frequency compared with the respective miR-137 targets (Additional file 1: Table S6, differentially expressed targets: let-7a 14%, miR-21-5p 11%, miR-93-5p 14%, miR-451a 0%, and miR-675-5p 0%). The numbers of differentially expressed target genes were significantly different between miR-137 and the five pooled controls (P = 0.031, Χ2 test, two-sided, Yates correction, Additional file 1: Table S6).
miR-137 acts cooperatively and synergistically with miR-124 and miR-128 [7, 33, 34]; therefore, changes in the expression of these two microRNAs may interfere with the expression of miR-137 target genes. We analyzed the precursor expression of both microRNAs in the CommonMind data set and found no evidence of differential expression. MIR124-2HG (ENSG00000254377) expression was not different between SCZ and control individuals (P = 0.8487, expression level 5.948). No data was available for MIR-128 (ENSG00000207654, ENSG00000207625) or for mature microRNAs. To investigate possible co-regulation with miR-137, we looked for additional binding sites in the 3′UTR of differentially expressed miR-137 target genes (Table 1). Several targets listed in Table 1 had putative binding sites for miR-124 (5 out of 16) and miR-128 (1 out of 16) (Additional file 1: Table S7).
Finally, we performed a functional network analysis of validated miR-137 targets (Table 2). Genes involved in cell death and survival, cellular movement, gene expression, and molecular transport were enriched. Many target genes were linked to developmental disorders, neurological diseases, and psychological disorders among which SCZ turned out to be the top hit.
In this study, we identified SHANK2 as a novel direct target of miR-137. We discovered that physiological levels of miR-137 regulate SHANK2 expression, most likely by repressing translation of SHANK2 protein rather than inducing mRNA degradation. Inhibiting effective protein translation is a known mechanism for local fine tuning of gene expression at postsynaptic sites and is in line with previously reported direct miR-137 targets, including Ephrin B2 (EFNB2) and the AMPA receptor subunit GluA1 (GRIA1) [35, 36]. miR-137 regulates the expression of several proteins that function at glutamatergic synapses, e.g., EFNB2, GluA1, and Mib1 and thereby influences neuronal maturation and signal transduction [8, 35, 36]. AMPA receptors are anchored to the postsynaptic Shank scaffold via PSD-95/Stargazin proteins . SHANK2 and AMPA receptors are both regulated by miR-137 and are important for synaptic maturation and plasticity; therefore, this control by miR-137 may have synergistic effects on synaptic regulation. Beside the regulation of postsynaptic genes, miR-137 has also been shown to regulate presynaptic genes and presynaptic neurotransmitter release . This is supported by the finding that altered miR-137 levels impact synaptic function and neuronal network formation in the mouse hippocampus .
Morphological effects caused by altered Shank2 levels in hippocampal neurons do not correlate to the morphological changes found when miR-137 levels are altered. Overexpression of SHANK2 increased spine volume while its downregulation led to reduced spine volume and increased dendritic arborization . In contrast, miR-137 overexpression or inhibition did not affect dendritic spine morphology, whereas miR-137 inhibition only led to a reduced spine density . This difference might be due to a more subtle regulatory effect of miR-137 on Shank2 protein expression compared to an RNAi-mediated knockdown as microRNAs are “fine-tuners” of protein expression. miR-137 regulates the expression of multiple glutamatergic synapse proteins; therefore, subtle regulation of SHANK2 expression by miR-137 is likely to be physiologically relevant. miR-137 also has subtle effects on the protein expression of its other targets . Different postsynaptic proteins are probably regulated by miR-137 at the same time, which may increase the intensity of the neuronal response.
We showed a regulatory influence of miR-137 on SHANK2 expression in mouse hippocampal neurons, whereas the relevance in other brain regions and in human neurons warrants future investigation.
miR-137 has been linked to various disorders including ID, ASD, and SCZ [6, 10, 11, 18]. Impaired synaptic plasticity and glutamatergic neurotransmission have been postulated as underlying pathological mechanisms . A heterozygous microdeletion on 1q21.3 encompassing the genes MIR137 and DPYD was previously described in ID and ASD patients [6, 10, 11]. Reduced levels of precursor and mature miR-137 concurrently with significantly increased levels of downstream target genes (MITF, EZH2, and KLF4) were determined in lymphoblastoid cells isolated from two patients with this 1q21.3 microdeletion . Based on our experimental results, we speculate that reduced miR-137 expression may also increase SHANK2 levels in these patients, which may contribute to the ID and ASD phenotype seen in these patients.
To further investigate the link between miR-137 and SCZ, we analyzed the expression of miR-137 precursor and known miR-137 target genes in postmortem DLPFC samples of SCZ individuals using the CommonMind gene expression data resource . No difference in miR-137 precursor expression was found between SCZ and control individuals. Our study was limited to miR-137 precursor expression analysis as no data of mature miR-137 expression was available. Previous studies have shown no difference in the levels of precursor and mature miR-137 between SCZ and control individuals [31, 41, 42]. miR-137 expression has been investigated in fibroblasts and fibroblast-derived neurons isolated from individuals homozygous for four schizophrenia-associated SNPs at the MIR137 locus. Increased endogenous miR-137 levels were identified in the minor compared to the major allele SNP group in the neurons, but not in fibroblasts . This indicates that regulation of miR-137 expression varies in different cell types. Analyzing miR-137 expression in specific cortical layers, cell types, or even cell compartments might reveal distinct local alterations in the DLPFC of SCZ patients.
Our network analysis revealed that many validated miR-137 target genes are linked to developmental, neurological, and psychiatric disorders, particularly SCZ. Therefore, we analyzed the expression of these genes in the DLPFC, a brain region that is affected in SCZ individuals. Expression of 23% (16/69) of the analyzed target genes was altered in the DLPFC of SCZ compared with control individuals. This change in expression was above the level of target enrichment in a pooled sample of five control microRNAs. This indicated that the miR-137 signaling network might be altered in the DLPFC of SCZ individuals. Five of these differentially expressed target genes have been associated with SCZ (TCF4, SIRT1, XIAP, GRIA1, ZNF804A) and two with ASD (RORA, KDM5B), suggesting some overlap between these two disorders. Four miR-137 target genes (RORA, CPLX1, TCF4, SIRT1) even show significant differences on the whole-transcriptome level, according to data provided by the CommonMind Consortium . Furthermore, analysis of the CommonMind RNA sequencing data confirmed elevated GRIA1 mRNA levels in the DLPFC of SCZ individuals , indicating impaired synaptic signaling in the DLPFC. The majority of differentially expressed miR-137 target genes (12/16) were only modestly elevated in SCZ individuals, suggesting that small effects of single genes accumulate in the DLPFC, presumably leading to a general impairment of miR-137 signaling. It is important to note that expression of miR-137 target genes is not regulated by miR-137 alone, but by multiple factors, e.g., other microRNAs, epigenetic mechanisms, subcellular localization, synaptic activity, medication, or other environmental factors. We analyzed miR-124 and miR-128, which act cooperatively with miR-137, and obtained no evidence that these two microRNAs influence the differential expression of miR-137 targets. miR-124 precursor expression was not different between SCZ and control individuals and the differentially expressed miR-137 target genes may only slightly be co-regulated by miR-124 and miR-128. Only 23% of the analyzed miR-137 target genes were differentially expressed on the mRNA level between SCZ and control individuals. Most miR-137 targets have been investigated in cancer cell lines and only 13 miR-137 target genes (including SHANK2) have been confirmed in neuronal cells; therefore, some of the confirmed targets may not be miR-137 targets in the brain. Expression of SHANK2 mRNA was not different in the DLPFC of SCZ and control individuals, which may mean that SHANK2 is not a relevant miR-137 target in the context of SCZ. However, we and others have identified neuronal miR-137 target genes which were regulated on protein [35, 36] and not on mRNA level. Therefore, we conclude that changes in mRNA expression do not reveal the full regulatory potential of miR-137. In this regard, proteomic data from postmortem human brains will be of high value for future SCZ studies.
We identified a direct regulatory link between microRNA-137 and SHANK2, which is of importance for a spectrum of disorders including ID, ASD, and SCZ. Furthermore, evidence was obtained that miR-137 target genes are differentially expressed between healthy controls and SCZ patients offering additional support for the involvement of miR-137 and its target genes in the pathogenesis of neuropsychiatric disorders. Further studies are warranted to address the mutual interaction in conditional knockout mouse models or human-induced pluripotent stem cell-derived neurons.
Autism spectrum disorders
Day in vitro
Dorsolateral prefrontal cortex
Dulbecco’s Modified Eagle’s Medium
Fetal bovine serum
miR-137 binding site-mutated 3′UTR
Real-time quantitative PCR
SHANK2 3′UTR wild type sequence
Bartel DP. MicroRNAs: target recognition and regulatory functions. Cell. 2009;136:215–33.
Hausser J, Zavolan M. Identification and consequences of miRNA-target interactions––beyond repression of gene expression. Nat Rev Genet. 2014;15:599–612.
Schratt G. Fine-tuning neural gene expression with microRNAs. Curr Opin Neurobiol. 2009;19:213–9.
Schratt G. microRNAs at the synapse. Nat Rev Neurosci. 2009;10:842–9.
Lagos-Quintana M, Rauhut R, Yalcin A, Meyer J, Lendeckel W, Tuschl T. Identification of tissue-specific microRNAs from mouse. Curr Biol. 2002;12:735–9.
Willemsen MH, Valles A, Kirkels LA, Mastebroek M, Olde Loohuis N, Kos A, Wissink-Lindhout WM, de Brouwer AP, Nillesen WM, Pfundt R, et al. Chromosome 1p21.3 microdeletions comprising DPYD and MIR137 are associated with intellectual disability. J Med Genet. 2011;48:810–8.
Silber J, Lim DA, Petritsch C, Persson AI, Maunakea AK, Yu M, Vandenberg SR, Ginzinger DG, James CD, Costello JF, et al. miR-124 and miR-137 inhibit proliferation of glioblastoma multiforme cells and induce differentiation of brain tumor stem cells. BMC Med. 2008;6:14.
Smrt RD, Szulwach KE, Pfeiffer RL, Li X, Guo W, Pathania M, Teng ZQ, Luo Y, Peng J, Bordey A, et al. MicroRNA miR-137 regulates neuronal maturation by targeting ubiquitin ligase mind bomb-1. Stem Cells. 2010;28:1060–70.
Crowley JJ, Collins AL, Lee RJ, Nonneman RJ, Farrell MS, Ancalade N, Mugford JW, Agster KL, Nikolova VD, Moy SS, Sullivan PF. Disruption of the microRNA 137 primary transcript results in early embryonic lethality in mice. Biol Psychiatry. 2015;77:e5–7.
Carter MT, Nikkel SM, Fernandez BA, Marshall CR, Noor A, Lionel AC, Prasad A, Pinto D, Joseph-George AM, Noakes C, et al. Hemizygous deletions on chromosome 1p21.3 involving the DPYD gene in individuals with autism spectrum disorder. Clin Genet. 2011;80:435–43.
Tucci A, Ciaccio C, Scuvera G, Esposito S, Milani D. MIR137 is the key gene mediator of the syndromic obesity phenotype of patients with 1p21.3 microdeletions. Mol Cytogenet. 2016;9:80.
Abrahams BS, Arking DE, Campbell DB, Mefford HC, Morrow EM, Weiss LA, Menashe I, Wadkins T, Banerjee-Basu S, Packer A. SFARI gene 2.0: a community-driven knowledgebase for the autism spectrum disorders (ASDs). Mol Autism. 2013;4:36.
Banerjee-Basu S, Packer A. SFARI gene: an evolving database for the autism research community. Dis Model Mech. 2010;3:133–5.
Wong N, Wang X. miRDB: an online resource for microRNA target prediction and functional annotations. Nucleic Acids Res. 2015;43:D146–52.
Guan F, Zhang B, Yan T, Li L, Liu F, Li T, Feng Z, Zhang B, Liu X, Li S. MIR137 gene and target gene CACNA1C of miR-137 contribute to schizophrenia susceptibility in Han Chinese. Schizophr Res. 2014;152:97–104.
Ripke S, O'Dushlaine C, Chambert K, Moran JL, Kahler AK, Akterin S, Bergen SE, Collins AL, Crowley JJ, Fromer M, et al. Genome-wide association analysis identifies 13 new risk loci for schizophrenia. Nat Genet. 2013;45:1150–9.
Schizophrenia Psychiatric Genome-Wide Association Study C. Genome-wide association study identifies five new schizophrenia loci. Nat Genet. 2011;43:969–76.
Schizophrenia Working Group of the Psychiatric Genomics C. Biological insights from 108 schizophrenia-associated genetic loci. Nature. 2014;511:421–7.
Cross-Disorder Group of the Psychiatric Genomics C, Genetic Risk Outcome of Psychosis C. Identification of risk loci with shared effects on five major psychiatric disorders: a genome-wide analysis. Lancet. 2013;381:1371–9.
Kreienkamp HJ. Scaffolding proteins at the postsynaptic density: shank as the architectural framework. Handb Exp Pharmacol. 2008:365–80.
Roussignol G, Ango F, Romorini S, Tu JC, Sala C, Worley PF, Bockaert J, Fagni L. Shank expression is sufficient to induce functional dendritic spine synapses in aspiny neurons. J Neurosci. 2005;25:3560–70.
Mahmoudi E, Cairns MJ. MiR-137: an important player in neural development and neoplastic transformation. Mol Psychiatry. 2017;22:44–55.
Grabrucker AM, Knight MJ, Proepper C, Bockmann J, Joubert M, Rowan M, Nienhaus GU, Garner CC, Bowie JU, Kreutz MR, et al. Concerted action of zinc and ProSAP/Shank in synaptogenesis and synapse maturation. EMBO J. 2011;30:569–81.
Zhang Y, Werling U, Edelmann W. SLiCE: a novel bacterial cell extract-based DNA cloning method. Nucleic Acids Res. 2012;40:e55.
Fromer M, Roussos P, Sieberts SK, Johnson JS, Kavanagh DH, Perumal TM, Ruderfer DM, Oh EC, Topol A, Shah HR, et al. Gene expression elucidates functional impact of polygenic risk for schizophrenia. Nat Neurosci. 2016;19:1442–53.
Agarwal V, Bell GW, Nam JW, Bartel DP. Predicting effective microRNA target sites in mammalian mRNAs. elife. 2015;4 https://doi.org/10.7554/eLife.05005.
Uhlen M, Fagerberg L, Hallstrom BM, Lindskog C, Oksvold P, Mardinoglu A, Sivertsson A, Kampf C, Sjostedt E, Asplund A, et al. Proteomics. Tissue-based map of the human proteome. Science. 2015;347:1260419.
Boeckers TM, Kreutz MR, Winter C, Zuschratter W, Smalla KH, Sanmarti-Vila L, Wex H, Langnaese K, Bockmann J, Garner CC, Gundelfinger ED. 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:6506–18.
Huang JT, Wang J, Srivastava V, Sen S, Liu SM. MicroRNA machinery genes as novel biomarkers for cancer. Front Oncol. 2014;4:113.
Bunney WE, Bunney BG. Evidence for a compromised dorsolateral prefrontal cortical parallel circuit in schizophrenia. Brain Res Brain Res Rev. 2000;31:138–46.
Guella I, Sequeira A, Rollins B, Morgan L, Torri F, van Erp TG, Myers RM, Barchas JD, Schatzberg AF, Watson SJ, et al. Analysis of miR-137 expression and rs1625579 in dorsolateral prefrontal cortex. J Psychiatr Res. 2013;47:1215–21.
O'Connor JA, Hemby SE. Elevated GRIA1 mRNA expression in layer II/III and V pyramidal cells of the DLPFC in schizophrenia. Schizophr Res. 2007;97:277–88.
Santos MC, Tegge AN, Correa BR, Mahesula S, Kohnke LQ, Qiao M, Ferreira MA, Kokovay E, Penalva LO. miR-124, -128, and -137 orchestrate neural differentiation by acting on overlapping gene sets containing a highly connected transcription factor network. Stem Cells. 2016;34:220–32.
Schouten M, Fratantoni SA, Hubens CJ, Piersma SR, Pham TV, Bielefeld P, Voskuyl RA, Lucassen PJ, Jimenez CR, Fitzsimons CP. MicroRNA-124 and -137 cooperativity controls caspase-3 activity through BCL2L13 in hippocampal neural stem cells. Sci Rep. 2015;5:12448.
Wu S, Zhang R, Nie F, Wang X, Jiang C, Liu M, Valenzuela RK, Liu W, Shi Y, Ma J. MicroRNA-137 inhibits EFNB2 expression affected by a genetic variant and is expressed aberrantly in peripheral blood of schizophrenia patients. EBioMedicine. 2016;12:133–42.
Olde Loohuis NF, Ba W, Stoerchel PH, Kos A, Jager A, Schratt G, Martens GJ, van Bokhoven H, Nadif Kasri N, Aschrafi A. MicroRNA-137 controls AMPA-receptor-mediated transmission and mGluR-dependent LTD. Cell Rep. 2015;11:1876–84.
Chen L, Chetkovich DM, Petralia RS, Sweeney NT, Kawasaki Y, Wenthold RJ, Bredt DS, Nicoll RA. Stargazin regulates synaptic targeting of AMPA receptors by two distinct mechanisms. Nature. 2000;408:936–43.
Siegert S, Seo J, Kwon EJ, Rudenko A, Cho S, Wang W, Flood Z, Martorell AJ, Ericsson M, Mungenast AE, Tsai LH. The schizophrenia risk gene product miR-137 alters presynaptic plasticity. Nat Neurosci. 2015;18:1008–16.
Berkel S, Tang W, Trevino M, Vogt M, Obenhaus HA, Gass P, Scherer SW, Sprengel R, Schratt G, Rappold GA. Inherited and de novo SHANK2 variants associated with autism spectrum disorder impair neuronal morphogenesis and physiology. Hum Mol Genet. 2012;21:344–57.
Coyle JT, Balu DT, Puhl MD, Konopaske GT. History of the concept of disconnectivity in schizophrenia. Harv Rev Psychiatry. 2016;24:80–6.
Beveridge NJ, Gardiner E, Carroll AP, Tooney PA, Cairns MJ. Schizophrenia is associated with an increase in cortical microRNA biogenesis. Mol Psychiatry. 2010;15:1176–89.
Santarelli DM, Beveridge NJ, Tooney PA, Cairns MJ. Upregulation of dicer and microRNA expression in the dorsolateral prefrontal cortex Brodmann area 46 in schizophrenia. Biol Psychiatry. 2011;69:180–7.
Nguyen A, Rauch TA, Pfeifer GP, Hu VW. Global methylation profiling of lymphoblastoid cell lines reveals epigenetic contributions to autism spectrum disorders and a novel autism candidate gene, RORA, whose protein product is reduced in autistic brain. FASEB J. 2010;24:3036–51.
Bi W, Cheung SW, Breman AM, Bacino CA. 4p16.3 microdeletions and microduplications detected by chromosomal microarray analysis: new insights into mechanisms and critical regions. Am J Med Genet A. 2016;170:2540–50.
Zweier C, Peippo MM, Hoyer J, Sousa S, Bottani A, Clayton-Smith J, Reardon W, Saraiva J, Cabral A, Gohring I, et al. Haploinsufficiency of TCF4 causes syndromal mental retardation with intermittent hyperventilation (Pitt-Hopkins syndrome). Am J Hum Genet. 2007;80:994–1001.
Kishi T, Fukuo Y, Kitajima T, Okochi T, Yamanouchi Y, Kinoshita Y, Kawashima K, Inada T, Kunugi H, Kato T, et al. SIRT1 gene, schizophrenia and bipolar disorder in the Japanese population: an association study. Genes Brain Behav. 2011;10:257–63.
Hussain MS, Baig SM, Neumann S, Peche VS, Szczepanski S, Nurnberg G, Tariq M, Jameel M, Khan TN, Fatima A, et al. CDK6 associates with the centrosome during mitosis and is mutated in a large Pakistani family with primary microcephaly. Hum Mol Genet. 2013;22:5199–214.
Rukova B, Staneva R, Hadjidekova S, Stamenov G, Milanova T, Toncheva D. Genome-wide methylation profiling of schizophrenia. Balkan J Med Genet. 2014;17:15–23.
Baek JH, Ha K, Kim Y, Yang SY, Cho EY, Choi Y, Ryu S, Lee YS, Park T, Hong KS. Association between the zinc finger protein 804A (ZNF804A) gene and the risk of schizophrenia and bipolar I disorder across diagnostic boundaries. Bipolar Disord. 2017;
Kang WS, Park JK, Kim SK, Park HJ, Lee SM, Song JY, Chung JH, Kim JW. Genetic variants of GRIA1 are associated with susceptibility to schizophrenia in Korean population. Mol Biol Rep. 2012;39:10697–703.
Vallianatos CN, Iwase S. Disrupted intricacy of histone H3K4 methylation in neurodevelopmental disorders. Epigenomics. 2015;7:503–19.
A.dS.C., F.-B.C., and S.B. are members of the Hartmut Hoffmann-Berling International Graduate School of Molecular and Cellular Biology (HBIGS). S.B. is indebted to the Baden-Württemberg Stiftung for the financial support of this research project by the Eliteprogram for Postdocs. G.A.R. is a member of CellNetworks Cluster of Excellence of the University of Heidelberg, Germany. We thank Birgit Weiss and Diana Porras Gonzalez for technical assistance. The study was supported by the EcTop3 Project of EXC81 CellNetworks (BMBF).
CommonMind Consortium––The data were generated as part of the CommonMind Consortium supported by funding from Takeda Pharmaceuticals Company Limited, F. Hoffman-La Roche Ltd., and NIH grants R01MH085542, R01MH093725, P50MH080405, R01MH097276, R01MH075916, P50MH096891, P50MH084053S1, R37MH057881 and R37MH057881S1, HHSN271201300031C, AG02219, AG05138, and MH06692. The brain tissue for the study was obtained from the following brain bank collections: the Mount Sinai NIH Brain and Tissue Repository, the University of Pennsylvania Alzheimer’s Disease Core Center, the University of Pittsburgh NeuroBioBank and Brain and Tissue Repositories, and the NIMH Human Brain Collection Core. CMC Leadership: Pamela Sklar, Joseph Buxbaum (Icahn School of Medicine at Mount Sinai), Bernie Devlin, David Lewis (University of Pittsburgh), Raquel Gur, Chang-Gyu Hahn (University of Pennsylvania), Keisuke Hirai, Hiroyoshi Toyoshiba (Takeda Pharmaceuticals Company Limited), Enrico Domenici, Laurent Essioux (F. Hoffman-La Roche Ltd.), Lara Mangravite, Mette Peters (Sage Bionetworks), Thomas Lehner, Barbara Lipska (NIMH).
This study was funded through the EcTop3 Project of EXC81 CellNetworks (BMBF).
Availability of data and materials
All data generated or analyzed during this study are included in this published article (and in Additional file 1).
Ethics approval and consent to participate
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Figure S1. Images of primary neuronal cultures pre‐ and posttreatment. Figure S2. High conservation of the miR‐137 binding site in the SHANK2‐3’UTR and of miR‐137 between different species. Figure S3. Relative expression levels of miR‐137 in different human tissues. Figure S4. Uncropped Western blot pictures. Table S1. ASD risk genes which are predicted or validated miR‐137 targets. Table S2. Primers used for cloning, mutagenesis, and screening. Primer sequences are all shown in 5′→3′ orientation. Table S3. Origin of total RNA samples used to measure hsa‐miR‐137 relative expression across different tissues (see Additional file 1: Figure S2 for results). Table S4. Experimentally validated miR‐137 targets. Table S5. Gene expression analysis of 69 validated miR‐137 target genes (including SHANK2) in the CommonMind RNA sequencing data. Table S6a. Gene expression analysis of validated targets from five different control microRNAs in the CommonMind RNA sequencing data. Genes labeled in gray withstand correction for multiple testing using the Benjamini-Hochberg method and a FDR of 10%. Table S6b. Comparison of the number of differentially expressed target genes of different microRNAs between SCZ and control individuals in the CommonMind RNASeq data. Table S7. Analysis of the 3′UTR of the differentially expressed miR-137 genes in the DLPFC between SCZ and control individuals for additional putative miR-124 and miR-128 binding sites. (PDF 1417 kb)