Human Molecular Genetics Advance Access originally published online on January 9, 2008
Human Molecular Genetics 2008 17(8):1156-1168; doi:10.1093/hmg/ddn005
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Dysregulation of miRNA 181b in the temporal cortex in schizophrenia
1 Schizophrenia Research Institute, Sydney, NSW 2006, Australia 2 School of Biomedical Sciences, Faculty of Health, and Hunter Medical Research Institute, The University of Newcastle, University Drive, Callaghan, NSW 2308, Australia 3 Department of Infectious Diseases and Immunology 4 Department of Pathology, The University of Sydney, Blackburn Building Level 6, Sydney, NSW 2006, Australia
* To whom correspondence should be addressed. Tel: +61 249218670; Fax: +61 249217903; Email: murray.cairns{at}newcastle.edu.au
Received November 1, 2007; Accepted January 4, 2008
| ABSTRACT |
|---|
|
|
|---|
Analysis of global microRNA (miRNA) expression in postmortem cortical grey matter from the superior temporal gyrus, revealed significant up-regulation of miR-181b expression in schizophrenia. This finding was supported by quantitative real-time RT–PCR analysis of miRNA expression in a cohort of 21 matched pairs of schizophrenia and non-psychiatric controls. The implications of this finding are substantial, as this miRNA is predicted to regulate many target genes with potential significance to the development of schizophrenia. They include the calcium sensor gene visinin-like 1 (VSNL1) and the ionotropic AMPA glutamate receptor subunit (GRIA2), which were found to be down-regulated in the same cortical tissue from the schizophrenia group. Both of these genes were also suppressed in miR-181b transfected cells and shown to contain functional miR-181b miRNA recognition elements by reporter gene assay. This study suggests altered miRNA levels could be a significant factor in the dysregulation of cortical gene expression in schizophrenia.
| INTRODUCTION |
|---|
|
|
|---|
Schizophrenia is a severely debilitating psychiatric disorder characterized by a diverse range of symptoms. While extensive research has not determined the definitive cause(s), it is generally accepted that a number of influences including genetic, environmental and developmental factors are involved. Genetics in particular has been known for many years to play a key role in the development of schizophrenia. The pattern of inheritance, however, is far from clear with a complex array of associated loci and putative candidate genes each with a relatively small affect, disbursed across the entire genome (1–3). This picture is further complicated by environmental and epigenetic influences evident by the relatively low rate of concordance in monozygotic twins (
50%) (4,5). Recent research from a variety of quarters has lent support to the neurodevelopmental hypothesis, which suggests that genetic and environmental influences affect neurodevelopmental processes that cause abnormalities in brain plasticity and connectivity leading to schizophrenia in early adulthood. In this model, regulatory factors that control the stage and tissue expression of molecules involved in neurodevelopmental processes will be particularly important. While transcription factors and protein based signal transduction pathways are usually considered in this context, a new class of gene expression modulators consisting entirely of RNA is also emerging as a significant regulatory factor. These small non-coding RNA molecules known as microRNAs (miRNA) function as guide sequences for the cellular gene silencing pathway to bring about target gene suppression through a process known as post-transcriptional gene silencing (PTGS) (6). These riboregulators are now believed to play a major role in coordinating the regulation of gene expression during the differentiation and development of the brain (7). Indeed, a large proportion of human miRNAs display a brain-specific or brain-enriched expression pattern (8). There is even evidence to suggest they are involved in brain function and neural plasticity via the regulation of molecules involved in long-term potentiation (LTP) (9,10) and in the establishment and maintenance of dendrites (11).
Considering that each miRNA is potentially manipulating the expression of hundreds of target genes, the clinical implications of an abnormality or disturbance of this system are substantial, particularly if such abnormality occurs during development of the central nervous system. Neurological manifestations of abnormalities in miRNA-mediated regulation have in fact already been identified for fragile X mental retardation (9) and Tourettes syndrome (12). The miRNA pathway is also implicated in DiGeorge syndrome as one of the deleted genes, DGCR8 is associated with the microprocessor complex with Drosha and is essential for processing of the primary miRNA transcript (13). Interestingly, 25% of patients with DiGeorge syndrome also develop schizophrenia (14). In view of these findings, it is reasonable to suspect that alteration in miRNA gene silencing pathways may also be involved in the development of complex psychiatric disorders such as schizophrenia. This is supported by a recent study of miRNA expression in the prefrontal cortex, which identified a number of altered miRNAs (15). In the study reported here, we examined miRNA expression in cortical grey matter from the superior temporal gyrus (STG). The STG contains the primary and secondary auditory cortex thought to be involved in the generation of auditory hallucinations and appears to have reduced volume in schizophrenia (16,17). This analysis revealed significant up-regulation of miR-181b, which was found to correlate inversely with expression of schizophrenia-associated miR-181b target genes GRIA2 and visinin-like 1 (VSNL1) in STG from the same schizophrenia cohort.
| RESULTS |
|---|
|
|
|---|
miRNA expression in the temporal cortex
High throughput miRNA expression analysis in postmortem cortical grey matter from the STG was accomplished using a custom DNA microarray platform (18). Fluorescence signals observed in duplicate array features at a level 2-fold that of background were considered to be indicative of miRNA expression. This expression profile consisted of 76 miRNAs out of a total of 262 probes corresponding to miRBase version 7.0. Many of these miRNAs have previously been identified as either brain-specific, brain-enriched or constitutively expressed (Supplementary Material, Table S1). This is exemplified by strong expression of miR-125a/b, miR-124a, miR-9/9*, miR-128a/b as well as miR-138, miR-219 and miR-338 which are brain-specific or brain-enriched. The constitutively expressed let-7 group of miRNAs were also well represented.
Differential miRNA expression in schizophrenia
Differential expression of miRNA in postmortem STG from the schizophrenia cohort compared with their matched non-psychiatric controls was determined through two-class unpaired analysis using statistical analysis of microarrays software (Significance Analysis of Microarrays, SAM, version 2.23) (19). This analysis revealed significant up-regulation of two miRNAs, hsa-let-7g and hsa-miR-181b in the schizophrenia group, to levels 1.8-fold (P = 0.008) and 2.8-fold (P = 0.001), respectively (Fig. 1A).
|
To confirm the differential expression status of mature let-7g and miR-181b observed by the microarray analysis, RNA from the same tissue was subjected to northern blot analysis with 32P-labelled locked nucleic acid (LNA)-modified probes. As miR-181b expression in the STG was in the low to mid range compared with the entire profile and below that of the constitutively expressed let-7g, the naïve membrane was hybridized first with the miR-181b probe. After stripping, the membrane was re-probed for let-7g and then the U6 snRNA loading control, respectively, with the latter being used to normalize the expression value of each miRNA (Fig. 1B). Quantitative analysis with normalization to the U6 band intensity confirmed a significantly higher expression of miR-181b in the schizophrenia samples with an average increase of 1.5-fold (P = 0.041) compared with the controls. In contrast, the difference in let-7g expression between the two groups was not significant (P = 0.12; Fig. 1C).
Validation of miR-181b over-expression in schizophrenia
To further validate the apparent association between altered miR-181b expression in the STG and schizophrenia, tissue from an additional 14 matched pairs was obtained and the RNA extracted to generate a substantially larger cohort (schizophrenia n = 21, controls n = 21). As less tissue was available from these new cases yielding lower quantities of total RNA, a highly sensitive and specific quantitative real-time RT–PCR-based approach to determining miR-181b expression was established. This method utilized a miR-181b reverse primer consisting of a 9-nucleotide target-specific 3' segment and an adjoining 9-nucleotide arbitrary tag sequence at the 5' terminus. While the sequence-specific 9-mer component primed low-temperature cDNA extension on the miRNA target, the adjacent 9-mer tag provided the additional length required for primer recognition during the higher temperature PCR analysis phase. The relative expression of miR-181b from this RT–PCR analysis was consistent with the northern blot analysis confirming a significant up-regulation of this miRNA in the schizophrenia group, with an average increase of 1.24-fold (P = 0.049) that was observed in 14 out of 21 matched pairs. However, removal of one control sample, CTR1, that had a higher expression of miR-181b (differed by more than 2.7 SD units from the control mean), elevated the average fold change in miR-181b expression in schizophrenia across the remaining 20 pairs to 1.39-fold (P = 0.029).
To determine if the differential expression observed in miR-181b was due to alteration in 181b precursor RNA transcription, primers specific to the pre-181b-1 (transcribed from 1q31.3) and pre-181b-2 (transcribed from 9q33.3) were used for quantitative real-time RT–PCR. These were both normalized against U6 snRNA expression as described for the mature miRNA expression. The results of this analysis (Fig. 2C) indicated that the pre-181b-2 transcript is in much greater abundance in the STG than the pre-181b-1 transcript, with >10-fold higher relative expression. When the expression levels of these precursor miRNA transcripts in the schizophrenia samples were compared with the controls, there was no significant difference (Fig. 2C).
|
MiR-181b target prediction and pathway analysis
In order to gather some clues about the functional implication of elevated miR-181b concentrations in the cerebral cortex, an extensive list of putative target genes was assembled from a variety of publicly available miRNA data bases using the online TargetCombo web service (http://www.diana.pcbi.upenn.edu/cgi-bin/TargetCombo.cgi) (20). This provides the option of combining target predictions from the leading homology search algorithms miRanda, PicTar, TargetScanS and DIANA-microT (21–24). Each of these programmes screen for motifs with various miRNA homology attributes in the phylogenetically conserved segments of 3' UTR sequences taken from the genome database. Collectively, these searches identified >800 putative target genes (Supplementary Material, Table S2). A number of these genes are known to play an important role in human brain function and development, and could conceivably be involved in schizophrenia. This was exemplified by genes involved in synaptic transmission, including gamma amino butyric acid A receptor (GABRA1), glutamate receptors (GRM5, GRM7, GRIK2, GRID and GRIA2), serotonin receptors (HTR1B and HTR2C) and the cannabinoid receptor (CNR1). Some predicted target genes are also thought to be involved in brain development and neurodevelopmental disorders, for example ataxin 1 (ATXN1), zic family member 1 (ZIC1), SLIT and NTRK-like family, member 1 (SLTRK1) and fragile X mental retardation 1 (FMR1). The latter two involved in Tourettes syndrome and fragile X mental retardation, respectively, are already known to be associated with miRNA-related dysfunction.
To further understand the potential of miR-181b regulation, the predicted gene lists were subjected to functional annotation clustering and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis using the sdatabase for annotation, visualization, and integrated discovery (DAVID; http://david.abcc.ncifcrf.gov) (25). Interestingly, of the 144 conserved target genes predicted by the online version of the miRanda algorithm at the Memorial Sloan-Kettering Cancer Center Computational biology web site (http://www.microrna.org) (21), the most highly enriched gene ontology cluster related to development, nervous system development, neurogenesis and differentiation (Supplementary Material, Table S2). This was followed by a cluster of terms relating more broadly to signalling, signal transduction and transcription. When the combination gene list, consisting of 789 miR-181b target genes was submitted, a similar picture emerged except the ranking of the top two clusters was reversed. The KEGG pathways were also charted and ranked according to their P-value. The top five most significant terms for the large list included: LTP, MAPK signalling pathway, axon guidance and neurodegenerative disorders. In the smaller miRanda list developmentally related GO terms again featured highly with the top five consisting of development, neural development, cell differentiation, system development and organ development.
While miRNA target prediction algorithms provide a useful starting point for understanding the function of particular miRNA, they are also prone to a degree of over prediction leading to false positives (23,26). In order to relate the increased miR-181b activity in the context of schizophrenia-associated changes in the STG, the list of putative targets for this miRNA was cross matched with genes already shown in microarray experiments to be down regulated in the same tissue (Supplementary Material, Table S2e; Bowden et al., manuscript in preparation). This revealed a number of interesting target genes potentially regulated by miR-181b in the same tissue including the ionotropic glutamate receptor (GRIA2), fragile X- related 2 and VSNL1. Alteration of GRIA2 and VSNL1 expression in the dorsolateral prefrontal cortex and in the hippocampus, respectively, have also been shown to be associated with schizophrenia (27,28).
Validation of schizophrenia associated changes in VSNL1 and GRIA2 expression
To confirm the schizophrenia-associated down-regulation of VSNL1 and GRIA2 from the microarray analysis (Bowden et al., manuscript in preparation), STG RNAs (21 matched pairs) were subjected to quantitative real-time RT–PCR (Fig. 3). In the case of VSNL1, a significant decrease in expression was observed in the schizophrenia group averaging 1.91-fold with respect to the control cohort (P = 0.037) (Fig. 3A). Similarly, GRIA2 was also decreased by an average of 1.44-fold (P = 0.047) in the schizophrenia cohort (Fig. 3B).
|
Repression of miR-181b target genes in vitro
To gain further biological insight into the miRNA 181b target genes, in vitro cell cultures were transfected with synthetic miR-181b (siRNA-181b) and the mRNA expression analyzed using genome-wide Sentrix bead chip arrays (Illumina) (29). Interestingly, both VSNL1 and GRIA2 expression was suppressed by 67 and 41%, respectively, in response to elevated miR-181b in HEK293 cells and SH-SY5Y cells, respectively (Fig. 4).
|
Biological activity of miR-181b elements in VSNL1 and GRIA2
While circumstantial evidence supported an association between miR-181b expression and the repression of predicted target genes VSNL1 and GRIA2, both in the STG in schizophrenia and in vitro, a reporter gene assay was established to further investigate the biological validity of this link. For this purpose, miRNA recognition elements (MRE) from each of the target genes predicted to bind miR-181b were cloned into the 3'UTR of the firefly luciferase gene of the pMIR-REPORT vector (Fig. 5A). Mutant versions of these MREs containing single nucleotide polymorphisms in the seed-pairing region were also cloned into the pMIR-REPORT vector (Fig. 5B and C). Each of these reporter constructs and the control plasmid (pRNL-TK) encoding renilla luciferase were co-transfected (10 and 2.5 ng, respectively) with either siRNA-181b or LNA-modified antisense 181b (250 ng), into HEK-293 cells and their relative luciferase activity measured. In accordance with expectation, relative luciferase expression from reporter constructs containing the VSNL-1 MRE was found to be significantly reduced with respect to the control (average –20%, P = 0.007) when co-transfected with si181b, while being significantly increased in the presence of antisense 181b (average 28%, P = 0.012) (Fig. 5D). The constructs containing the mutant VSNL1 MRE displayed a smaller response to transfection, with suppression of 8% after co-transfection with si181b (P = 0.033), whereas it did not show a significant change in response to antisense 181b (Fig. 5D). The GRIA2 reporter gene displayed a similar response to siRNA-181b transfection, with an average suppression of 21% (P = 0.006), whereas antisense 181b transfection increased its expression 24% (P = 0.038). Surprisingly, luciferase activity from the construct containing the mutant GRIA2 MRE was also suppressed by si181b transfection (down 24%, P = 0.003), though the antisense 181b did not induce a significant change in its expression (Fig. 5E).
|
| DISCUSSION |
|---|
|
|
|---|
STG, miRNA expression and schizophrenia
In recent years, a number of genome-wide expression studies in schizophrenia including those from our own laboratory, have shown alterations in large numbers of genes (30,31). In view of the breadth of changes in the neural genomes activity in schizophrenia, it is reasonable to speculate that these could at least in part be due to changes in post-transcriptional gene silencing. To test this hypothesis, we initiated an investigation of miRNA expression patterns in the normal human cerebral cortex and in a matched cohort of tissue derived from subjects with schizophrenia. Postmortem grey matter from the STG of both normal and schizophrenia groups were dissected and the RNA extracted and subjected to high throughput analysis on custom miRNA-specific microarrays (18). While many postmortem expression studies of schizophrenia have focused on the prefrontal cortex, there is evidence to suggest the STG is also a participant in the pathophysiology of schizophrenia. Genome-wide analysis of schizophrenia-associated gene expression in the STG recently reported a large number of altered genes (Bowden et al., manuscript in preparation). In this study, we used the exact same tissue to enable a direct comparison between the schizophrenia-associated changes in gene and miRNA expression. The microarray analysis revealed up-regulation of let-7g and miR-181b, and this was also supported for miR-181b by northern blot hybridization and further validated by quantitative real-time RT–PCR in a larger cohort of 21 matched pairs.
miRNA 181b is encoded at two genomic loci, which give rise to primary transcripts for pre-miRNAs 181b-1 (MIRN181B1; 1q31.3) and pre-miRNA 181b-2 (MIRN181B2; 9q33.3). As the 181b-1 locus resides in the vicinity of a schizophrenia linkage region (1q23.3-1q31.1) (32) and the schizophrenia-associated genes for interleukin 10 (IL10; 1q31-32) and plexin A2 (PLXNA2; 1q32.2) (33–35), it is conceivable that its transcription could have been affected by schizophrenia-associated genetic alterations present in this region. To determine the relative contributions of each 181b transcript in the context of schizophrenia, we performed precursor-specific quantitative real-time RT–PCR. This analysis showed that the MIRN181B2 locus on chromosome 9 was actually far more active (
10-fold) than its counterpart on chromosome 1. While there was some increase in the average pre-miR-181b-2 expression in the schizophrenia group, this was not statistically significant. One interpretation of this finding is that the up-regulation of mature miR-181b expression was not due to changes at the level of transcription, but rather changes that affected its processing or maturation. This would be consistent with the hypothesis developed to explain miRNA expression changes observed in the dorsolateral prefrontal cortex (DLPFC) in schizophrenia (15). In this study, a subset of down-regulated miRNAs (which incidentally did not include miR-181b) shared an upstream consensus sequence, which was suggestive of a role in miRNA maturation. This concept is also supported by the loss of the DiGeorge critical region 8 (DGCR8; an accessory for the primary miRNA processing protein Drosha) (13) in the 21q11 micro deletion responsible for DiGeorge or Velo Cardio Facial syndrome, which is known to be the highest risk factor for schizophrenia (
25%) (14) apart from having a monozygotic twin with the disorder (
50%). However, a more specific effect could be operating in respect to the changes in miR-181b expression, observed here, perhaps through changes to pre-miRNA sequence elements involved in maturation of the miRNA. More general deficits in pri-miRNA processing suggested for the DLPFC were not supported here as pre-miR-181b was not altered in the STG and it does not contain the consensus sequence reported for miRNA altered in the DLPFC (15).
MiR-181b target genes
The implications of even a small increase in miR-181b expression in the cerebral cortex could be profound, particularly during development as it is capable of extending a repressive influence over hundreds of different target genes. While miRNA 181b is known to be expressed in a variety of mammalian tissues, a number of studies have shown brain enrichment (8,36,37) with expression in both neurons and neuroglia (38,39). It is also highly expressed in the thymus and lungs, and has been shown to be directly associated with B-cell lineage differentiation of hematopoietic progenitor cells (40). miR-181b was also found to be up-regulated in acute promyelocytic leukaemia cells in response to retinoic acid-induced differentiation (41). It was also up in thyroid papilliary carcinoma cells compared with those from normal thyroid tissue (42). Conversely, miR-181b was observed to be down-regulated in glioblastoma cells (39). Most of these findings, however, are consistent with a role for miR-181b in differentiation and development. This is broadly supported by miR-181b target gene predictions, which implicate an apparent bias towards genes involved in development, brain development and brain function.
While it is interesting to consider the potential of genome wide interactions, we were able to examine miRNA targets specifically in the context of genes shown previously to be down-regulated in schizophrenia and in the same tissue from the same cohort of controls and schizophrenia subjects (Supplementary Material, Table S2e; Bowden et al., manuscript in preparation). Among the 23 genes in this category, two targets including GluR2, the ionotropic glutamate/AMPA receptor gene (GRIA2) and the calcium sensor/trkB mRNA binding protein known as VSNL1 were particularly interesting in the context of schizophrenia. GRIA2 is a major ionotropic glutamate receptor subunit involved in fast excitatory neurotransmission. It has been shown to have an important function in the development of synaptic plasticity because of its involvement with NMDA receptors in the establishment of LTP, long-term depression (43) and by directly stimulating increased growth and density of dendritic spines (44). GRIA2 has also shown a consistent association with schizophrenia and is compatible with the glutamate hypofunction hypothesis, with an observed decrease in both mRNA and protein levels in postmortem samples from a number of brain regions including the medial temporal lobe, hippocampus and dorsolateral prefrontal cortex (27,45,46). The calcium sensor protein VSNL1 has also been shown to be differentially expressed in schizophrenia (28). Specific changes in VSNL1 expression have been reported in the rat brain in phencyclidine (PCP) and ketamine models of schizophrenia (47,48). VSNL1 is thought to be a calcium sensitive signal transduction molecule. Its location within hippocampal neurons has also been shown to be altered in response to stimulation by glutamate (49). Interestingly, VSNL1 is also a calcium-dependent double-stranded RNA-binding protein, that may provide activity-dependent trafficking of certain neuronal mRNAs to the dendrites and more specifically is known to bind the 3' UTR of mRNA for the neurotrophin receptor (trkB) (50).
To further support the plausibility of a relationship between GRIA2 and VSNL1 expression in the STG and miR-181b, we established a synthetic miRNA 181b transfection system and monitored the expression of theses genes in response to changes in miRNA 181b concentration in vitro. The expression of both of these genes responded in accordance with expectation by showing a substantial drop in expression after the addition of the miRNA (Fig. 4). However, in order to directly characterize the miR-181b target status of GRIA2 and VSNL1, the putative MRE from their respective 3' UTR were cloned downstream of the firefly luciferase gene to measure their response to miR-181b in vitro. The presence of either of these MRE and si181b was observed to significantly reduce expression of luciferase expression compared with the control siRNA. In both cases reporter gene expression was also elevated in response to antisense 181b transfection, presumably as a result of depleting the bioavailability of endogenous miR-181b. In most cases, these responses were absent or diminished substantially in the mutant MRE constructs, which contained point mutations in the miRNA seed regions. The one exception to this was the reporter gene carrying the mutant version of the GRIA2 MRE, which was also silenced significantly by si181b. This was not surprising as the mutant was predicted to be a viable miR-181b target gene (albeit a weaker one) and thus capable of responding to saturating levels of the cognate miRNA. In contrast, at lower endogenous concentrations it did not display any signs of the silencing observed with the wild-type reporter gene, evident by the absence of response to the antisense miR-181b. From these experiments, we were able to conclude that both VSNL1 and GRIA2 3'UTR segments are sufficient to support PTGS by miR-181b in vitro, enhancing their status as functioning targets of this miRNA.
In summary, we have considered the possibility that alteration of miRNA-mediated PTGS is associated with the pathophysiology of schizophrenia. In support of this hypothesis, miR-181b was shown to be up-regulated in grey matter from the STG. In a similar approach, miRNA expression was also shown to be altered in the DLPFC in schizophrenia, although changes in miRNA 181b expression was not reported in this tissue (15). Genetic analysis of polymorphisms in the vicinity of brain expressed miRNA genes hsa-miR-206 and hsa-miR-198, was recently determined to be weakly associated with schizophrenia (51). There has also been some interest in the analysis of hsa-miR-103b in the context of schizophrenia due to its location in the 22q11 locus, however, no association was observed (52). Time will tell if schizophrenia-associated changes in miRNA expression are due to direct genetic influence or some other upstream regulatory, epigenetic or mechanistic factors affecting miRNA maturation. In either case it has important implications for understanding the neurodevelopmental origins of the schizophrenia, particularly as changes in a given miRNA can affect the expression of hundreds of target genes. In the case of miR-181b there are up to 800 conserved targets predicted in the human genome that have the potential to interact with the miRNA. The true extent of its influence, however, will depend on the biological context and other so far unidentified factors. In respect to schizophrenia and changes in the STG, we have identified two important candidate target genes for miR-181b (GRIA2 and VSNL1) that are suppressed in the same tissue, and may be responding to changes in the local miRNA environment. If these specific effects are a manifestation of the observed changes in miRNA expression, they probably represent the tip of the iceberg in terms of their global regulatory influence. The full impact of this and other changes in miRNA expression will no doubt take some time to unravel and appreciate fully.
| MATERIALS AND METHODS |
|---|
|
|
|---|
Tissue collection
Fresh frozen postmortem STG grey matter tissue from 21 subjects with schizophrenia and 21 non-psychiatric controls was obtained through the NSW Tissue Resource Centre, The University of Sydney, Australia. The grey matter tissue was taken from the outer edge of blocks of STG tissue from the most caudal coronal brain slice containing the STG (Brodmanns Area 22). In all cases, a diagnosis of schizophrenia in accordance with DSM-IV criteria was confirmed by medical file review using the Item Group Checklist of the Schedules for Clinical Assessment in Neuropsychiatry and the Diagnostic Instrument for Brain Studies. Consent was obtained from the next of kin and subjects with a significant history of drug or alcohol abuse, or other condition or gross neuropathology that might have influenced agonal state were excluded. In addition, control subjects were excluded if there was a history of alcoholism or suicide. All subjects were of Caucasian descent. Subjects with schizophrenia were matched for gender, age, brain hemisphere, postmortem interval and pH (Table 1). This cohort of tissue contained 13 matched pairs that were previously analyzed for mRNA expression using microarray analysis (Bowden el al., manuscript in preparation).
|
Tissue dissection and RNA extraction
Postmortem cortical grey matter was dissected from the outer edge of frozen coronal sections (1 cm) using a fine diameter hole punch and scalpel. In each case,
50–60 mg grey matter was removed and immediately homogenized in 1 ml of Trizol reagent and the total RNA extracted according to the manufacturers instructions (Invitrogen). The RNA concentration and integrity was determined using an Experion bioanalyser (BioRad).
miRNA expression arrays
miRNAs were labelled directly using a ligation approach consisting of 9 µg of total RNA, in 50 mM HEPES, pH 7.8, 3.5 mM DTT, 20 mM MgCl2, 0.1 mM ATP, 10 µg/ml BSA, 10% DMSO, 500 ng 5'-phosphate-cytidyl-uridyl-Cy3-3' (Dharmacon) and 20 units T4 RNA ligase (Fermentas) (18,53). After incubating for 2 h on ice the labelled RNA was precipitated with 0.3 M sodium acetate, 2 volumes 100% ethanol and 20 µg glycogen at –20°C overnight. A synthetic reference library consisting of DNA oligonucleotides (representing the entirety of miRBase version 7.0) was labelled with Ulysis platinum conjugated AlexaFluor 647 (equivalent to Cy5) for detection in the control channel, using the labelling kit, according to the manufacturers instructions (Invitrogen). Unconjugated label was then removed by gel filtration through a Sephadex G-25 spin column (GE Healthcare). The labelled reference library was used at a 1/700 dilution, along side the Cy3 labelled miRNAs, in each array hybridization.
Microarrays were prepared using anti-sense DNA oligonucleotides corresponding to the miRBase Version 7.0 (Sanger Institute, UK) containing 261 human miRNAs sequences. The oligonucleotide probes were printed in duplicate onto GAPS-2 glass slides (Corning). The slides were then prepared and hybridized with the labelled miRNA and synthetic controls as described previously (18). Briefly, slides were pre-hybridized in 3x SSC, 0.1% SDS and 0.2% BSA for 1 h at 65°C and washed four times with RNAse-free water, once with 100% ethanol, and dried by centrifugation at 150g for 5 min. Hybridization chambers were created around each array using 17 mmx 28 mm disposable frame seals and cover slides (Bio-Rad). The labelled RNA sample was added to 100 µl hybridization buffer (400 mM Na2HPO4, pH 7.0, 0.8% BSA, 5% SDS, 12% formamide) and heated for 4 min at 95°C (in the dark). The mixture was injected into the chamber and hybridized for 2 h at 37°C in a rotary hybridization oven. The coverslips and frames were removed and the slides washed once in 2x SSC, 0.025% SDS at room temperature, three times in 0.8x SSC at room temperature and three times in ice-cold 0.4x SSC. Each slide was then dried by centrifugation for 10 min at 60g. Arrays were then scanned with a Genepix 4000B Scanner (Axon Instruments) and raw pixel intensities extracted with Genepix Pro 3.0 software (Axon Instruments).
A miRNA was considered expressed if its raw Cy3 pixel intensity was at least 200% above background. Raw Cy3 median pixel intensity values were background subtracted and normalized by median centring with respect to arrays using Cluster version 2.2 (Stanford University). Differential miRNA expression was analyzed using Significance Analysis of Microarrays SAM version 2.23 (Stanford University) (19) (available from http://www-stat.stanford.edu/~tibs/SAM/). The threshold for significance was set at 5% and a two-class comparison was performed using 5000 permutations of the data. A list of significantly altered miRNAs was compiled (false-discovery rate <5%).
Northern hybridization
Total RNA from each case (30 µg) was combined with equal volumes of loading dye (0.01% bromophenol blue, 10 mM EDTA and formamide), pre-heated at 95°C for 5 min then electrophoresed on a 16% denaturing (8.3 M urea) polyacrylamide sequencing gel. The RNA was then electro-transferred in a semi-dry blotter at 400 mA for 1 h to GeneScreen Plus nylon membrane (Perkin Elmer NEN) and immobilized by UV-cross linking, followed by baking at 80°C for 1 h. The membrane was then pre-hybridized and hybridized in PerfectHyb Plus hybridization buffer (Sigma) for 3 and 16h, respectively. The later was carried out with the addition of 50 pmol of 32P-labelled unmodified or LNA-modified antisense oligonucleotide prepared earlier using [
-32P] ATP (Perkin Elmer) and polynucleotide kinase (Fermentas) as described by previously (54). After low and high stringency washes, the radiolabelled membranes were imaged using a typhoon phosphorimager (Amersham) and analyzed using ImageQuant software (Amersham). The sequences of antisense probes for the detection of miR-181b, let-7g and U6 snRNA are presented in Table 2. Statistical analysis, consisting of a one-tailed t-test, was performed on the normalized intensity values to determine the significance of observed differences in average expression.
|
Quantitative real-time RT–PCR
Multiplex reverse transcription was performed on 500 ng of DNaseI-treated total RNA using either random hexamers (mRNA analysis), or a combination of reverse primers (miRNA analysis) specific for mature hsa-miR-181b, the U6 snRNA and β-actin, to a final concentration of 40 nM each (for sequences see Table 2). Reactions were performed using Superscript II reverse transcriptase in 1x first-strand buffer according to the manufacturers instructions (Invitrogen). Real-time PCR was performed essentially as previously described (31), in triplicate on diluted cDNA combined with Power SybrGreen master mix (Applied Biosystems) with 1 µM of the appropriate forward and reverse primers (Table 2), in a final volume of 25 µl using an ABI prism 7500 sequence detection system (PE Applied Biosystems). Relative miRNA expression was determined by the difference between their individual cycle threshold (Ct) value and that produced in the same sample for the U6 snRNA (
Ct). Similarly, relative mRNA expression ratio was normalized with respect to the geometric mean of β-actin and U6 snRNA expression. Differential expression of a given miRNA or mRNA was determined by the difference between the mean
Ct for the schizophrenia and control cohorts (
Ct) expressed as a ratio (2–
Ct) (55). To determine the significance of any difference in average expression in a given direction between the two cohorts, a paired one-tailed t-test was applied. No significant differences in the expression of normalizing genes (β-actin and U6 snRNA) were observed with respect to each other, between the schizophrenia and control cohorts.
Cell culture and siRNA transfection
HEK-293 and SH-SY5Y cell cultures were maintained as confluent monolayers at 37°C with 5% CO2 and 90% humidity in DMEM with 10% (vol/vol) fetal calf serum, 20 mM HEPES, 0.15% (wt/vol) sodium bicarbonate and 2 mM L-glutamine. HEK-293 cells were seeded into 10 cm petri dishes and transfected 24 h later using Lipofectamine 2000 (Invitrogen). SH-SY5Y cells were harvested and electroporated using the Nucleofector Kit V (Amaxa), before being seeded into 6-well plates at 1 x 106 cells/well. In each case transfections were performed according to manufacturers instructions with 100 nM siRNA (si181b) oligonucleotide (Table 2).
Target gene expression profiling
RNA was extracted directly from plates 24 h post-transfection using 2 ml of Trizol and a disposable cell scraper (Greiner Bio-One). miR-181b expression levels in transfected cells were analyzed using real-time PCR as described above. RNA was purified further using an RNAeasy MinElute kit (Qiagen) before amplification-labelling with a TotalPrep amplification kit (Ambion) and hybridization on HumanRef-8 whole-genome expression arrays (Illumina) according to the manufacturers instructions. Data were normalized and analyzed using Illumina Beadstudio 3.0 and GeneSpringGX 7.3.1 (Agilent Technologies, USA).
Target gene reporter assay
Validation of predicted miR-181b target genes VSNL1 and GRIA2 was accomplished by co-transfecting HEK293 cells with synthetic si181b or an LNA-modified antisense inhibitor and recombinant firefly luciferase reporter gene constructs containing 3' UTR sequences substituted from the target gene. Oligonucleotides encoding target gene MRE (or mutant controls) were annealed to form SpeI and HindIII restricted overhangs of a ligatable cassette compatible with SpeI and HindIII digested pMIR-REPORT vector (Ambion). Reporter gene silencing in response to miRNA co-transfection was monitored with respect to a spiked-in control plasmid expressing renilla luciferase using the dual luciferase reporter assay (Promega). To control for non-specific effects associated with siRNA transfection, the controls were co-transfected with siEGFP or miR-9* siRNA predicted to have little or no activity against the fusion transcripts tested here. HEK293 cells were cultured and transfected as described above (except in 24-well plates) using Lipofectamine2000 (Invitrogen).
| SUPPLEMENTARY MATERIAL |
|---|
|
|
|---|
Supplementary Material is available at HMG Online.
| FUNDING |
|---|
|
|
|---|
This study was supported by the Schizophrenia Research Institute, utilizing funding from NSW Health and the Henderson Foundation; a NARSAD Young Investigator Award (M.C.); a University of Newcastle pilot grant; and the M.C. Ainsworth Research Fellowship in Epigenetics (M.C.). Funding for the open access publication charge was provided by the Schizophrenia Research Institute.
| ACKNOWLEDGEMENTS |
|---|
Tissues were received from the Australian Brain Donor Programs NSW Tissue Resource Centre, which is supported by The University of Sydney, National Health and Medical Research Council of Australia, Schizophrenia Research Institute, National Institute of Alcohol Abuse and Alcoholism and NSW Department of Health. The authors would like to thank Professor Vaughan Carr for critical reading of the manuscript.
Conflict of Interest statement. None declared.
| REFERENCES |
|---|
|
|
|---|
- Badner J.A., Gershon E.S. Meta-analysis of whole-genome linkage scans of bipolar disorder and schizophrenia. Mol. Psychiatry (2002) 7:405–411.[CrossRef][ISI][Medline]
- Lewis C.M., Levinson D.F., Wise L.H., DeLisi L.E., Straub R.E., Hovatta I., Williams N.M., Schwab S.G., Pulver A.E., Faraone S.V., et al. Genome scan meta-analysis of schizophrenia and bipolar disorder, part II: Schizophrenia. Am. J. Hum. Genet. (2003) 73:34–48. Epub June 11, 2003.[CrossRef][ISI][Medline]
- Harrison P.J., Weinberger D.R. Schizophrenia genes, gene expression, and neuropathology: on the matter of their convergence. Mol. Psychiatry (2005) 10:40–68. image 5.[CrossRef][ISI][Medline]
- Tsuang M. Schizophrenia: genes and environment. Biol. Psychiatry (2000) 47:210–220.[CrossRef][ISI][Medline]
- Petronis A. The origin of schizophrenia: genetic thesis, epigenetic antithesis, and resolving synthesis. Biol. Psychiatry (2004) 55:965–970.[CrossRef][ISI][Medline]
- He L., Hannon G.J. MicroRNAs: small RNAs with a big role in gene regulation. Nat Rev Genet (2004) 5:522–531.[CrossRef][Medline]
-
Giraldez A.J., Cinalli R.M., Glasner M.E., Enright A.J., Thomson J.M., Baskerville S., Hammond S.M., Bartel D.P., Schier A.F. MicroRNAs regulate brain morphogenesis in zebrafish. Science (2005) 308:833–838. Epub March 17, 2005.
[Abstract/Free Full Text] - Sempere L.F., Freemantle S., Pitha-Rowe I., Moss E., Dmitrovsky E., Ambros V. Expression profiling of mammalian microRNAs uncovers a subset of brain-expressed microRNAs with possible roles in murine and human neuronal differentiation. Genome Biol (2004) 5:R13. Epub February 16, 2004.[CrossRef][Medline]
- Jin P., Zarnescu D.C., Ceman S., Nakamoto M., Mowrey J., Jongens T.A., Nelson D.L., Moses K., Warren S.T. Biochemical and genetic interaction between the fragile X mental retardation protein and the microRNA pathway. Nat. Neurosci. (2004) 7:113–117. Epub January 4, 2004.[CrossRef][ISI][Medline]
- Ashraf S.I., McLoon A.L., Sclarsic S.M., Kunes S. Synaptic protein synthesis associated with memory is regulated by the RISC pathway in Drosophila. Cell (2006) 124:191–205.[CrossRef][ISI][Medline]
- Schratt G.M., Tuebing F., Nigh E.A., Kane C.G., Sabatini M.E., Kiebler M., Greenberg M.E. A brain-specific microRNA regulates dendritic spine development. Nature (2006) 439:283–289.[CrossRef][Medline]
-
Abelson J.F., Kwan K.Y., ORoak B.J., Baek D.Y., Stillman A.A., Morgan T.M., Mathews C.A., Pauls D.L., Rasin M.R., Gunel M., et al. Sequence variants in SLITRK1 are associated with Tourettes syndrome. Science (2005) 310:317–320.
[Abstract/Free Full Text] - Gregory R.I., Yan K.P., Amuthan G., Chendrimada T., Doratotaj B., Cooch N., Shiekhattar R. The Microprocessor complex mediates the genesis of microRNAs. Nature (2004) 432:235–240. Epub November 7, 2004.[CrossRef][Medline]
-
Murphy K.C., Jones L.A., Owen M.J. High rates of schizophrenia in adults with velo-cardio-facial syndrome. Arch. Gen. Psychiatry (1999) 56:940–945.
[Abstract/Free Full Text] - Perkins D.O., Jeffries C.D., Jarskog L.F., Thomson J.M., Woods K., Newman M.A., Parker J.S., Jin J., Hammond S.M. microRNA expression in the prefrontal cortex of individuals with schizophrenia and schizoaffective disorder. Genome Biol (2007) 8:R27.[CrossRef][Medline]
- Rajarethinam R.P., DeQuardo J.R., Nalepa R., Tandon R. Superior temporal gyrus in schizophrenia: a volumetric magnetic resonance imaging study. Schizophr. Res. (2000) 41:303–312.[CrossRef][ISI][Medline]
-
Honea R., Crow T.J., Passingham D., Mackay C.E. Regional deficits in brain volume in schizophrenia: a meta-analysis of voxel-based morphometry studies. Am. J. Psychiatry (2005) 162:2233–2245.
[Abstract/Free Full Text] - Thomson J.M., Parker J., Perou C.M., Hammond S.M. A custom microarray platform for analysis of microRNA gene expression. Nat Methods (2004) 1:47–53. Epub September 29, 2004.[CrossRef][ISI][Medline]
-
Tusher V.G., Tibshirani R., Chu G. Significance analysis of microarrays applied to the ionizing radiation response. Proc. Natl. Acad. Sci. USA (2001) 98:5116–5121. Epub April 17, 2001.
[Abstract/Free Full Text] - Sethupathy P., Megraw M., Hatzigeorgiou A.G. A guide through present computational approaches for the identification of mammalian microRNA targets. Nat. Methods (2006) 3:881–886.[CrossRef][ISI][Medline]
- John B., Enright A.J., Aravin A., Tuschl T., Sander C., Marks D.S. Human MicroRNA targets. PLoS Biol (2004) 2:e363. Epub October 5, 2004.[CrossRef][Medline]
- Krek A., Grun D., Poy M.N., Wolf R., Rosenberg L., Epstein E.J., MacMenamin P., da Piedade I., Gunsalus K.C., Stoffel M., et al. Combinatorial microRNA target predictions. Nat. Genet. (2005) 37:495–500. Epub April 3, 2005.[CrossRef][ISI][Medline]
- Lewis B.P., Shih I.H., Jones-Rhoades M.W., Bartel D.P., Burge C.B. Prediction of mammalian microRNA targets. Cell (2003) 115:787–798.[CrossRef][ISI][Medline]
-
Kiriakidou M., Nelson P.T., Kouranov A., Fitziev P., Bouyioukos C., Mourelatos Z., Hatzigeorgiou A. A combined computational-experimental approach predicts human microRNA targets. Genes Dev. (2004) 18:1165–1178. Epub May 6, 2004.
[Abstract/Free Full Text] - Dennis G. Jr, Sherman B.T., Hosack D.A., Yang J., Gao W., Lane H.C., Lempicki R.A. DAVID: database for annotation, visualization, and integrated discovery. Genome Biol (2003) 4:3. Epub April 3, 2003.[CrossRef]
-
Yousef M., Nebozhyn M., Shatkay H., Kanterakis S., Showe L.C., Showe M.K. Combining multi-species genomic data for microRNA identification using a Naive Bayes classifier. Bioinformatics (2006) 22:1325–1334. Epub March 16, 2006.
[Abstract/Free Full Text] - Vawter M.P., Crook J.M., Hyde T.M., Kleinman J.E., Weinberger D.R., Becker K.G., Freed W.J. Microarray analysis of gene expression in the prefrontal cortex in schizophrenia: a preliminary study. Schizophr. Res. (2002) 58:11–20.[CrossRef][ISI][Medline]
- Bernstein H.G., Braunewell K.H., Spilker C., Danos P., Baumann B., Funke S., Diekmann S., Gundelfinger E.D., Bogerts B. Hippocampal expression of the calcium sensor protein visinin-like protein-1 in schizophrenia. Neuroreport (2002) 13:393–396.[CrossRef][ISI][Medline]
-
Kuhn K., Baker S.C., Chudin E., Lieu M.H., Oeser S., Bennett H., Rigault P., Barker D., McDaniel T.K., Chee M.S. A novel, high-performance random array platform for quantitative gene expression profiling. Genome Res. (2004) 14:2347–2356.
[Abstract/Free Full Text] - Mirnics K., Middleton F.A., Marquez A., Lewis D.A., Levitt P. Molecular characterization of schizophrenia viewed by microarray analysis of gene expression in prefrontal cortex. Neuron (2000) 28:53–67.[CrossRef][ISI][Medline]
- Weidenhofer J., Bowden N.A., Scott R.J., Tooney P.A. Altered gene expression in the amygdala in schizophrenia: up-regulation of genes located in the cytomatrix active zone. Mol. Cell. Neurosci. (2006) 31:243–250. Epub October 19, 2005.[CrossRef][ISI][Medline]
- Schwab S.G., Hallmayer J., Albus M., Lerer B., Eckstein G.N., Borrmann M., Segman R.H., Hanses C., Freymann J., Yakir A., et al. A genome-wide autosomal screen for schizophrenia susceptibility loci in 71 families with affected siblings: support for loci on chromosome 10p and 6. Mol. Psychiatry (2000) 5:638–649.[CrossRef][ISI][Medline]
- Yu L., Yang M.S., Zhao J., Shi Y.Y., Zhao X.Z., Yang J.D., Liu Z.J., Gu N.F., Feng G.Y., He L. An association between polymorphisms of the interleukin-10 gene promoter and schizophrenia in the Chinese population. Schizophr. Res. (2004) 71:179–183.[CrossRef][ISI][Medline]
- He G., Zhang J., Li X.W., Chen W.Y., Pan Y.X., Yang F.P., Gu N.F., Feng G.Y., Yang S.L., He J.Y., et al. Interleukin-10 -1082 promoter polymorphism is associated with schizophrenia in a Han Chinese sib-pair study. Neurosci. Lett. (2006) 394:1–4. Epub December 27, 2005.[CrossRef][ISI][Medline]
- Mah S., Nelson M.R., Delisi L.E., Reneland R.H., Markward N., James M.R., Nyholt D.R., Hayward N., Handoko H., Mowry B., et al. Identification of the semaphorin receptor PLXNA2 as a candidate for susceptibility to schizophrenia. Mol. Psychiatry (2006) 11:471–478.[CrossRef][ISI][Medline]
-
Krichevsky A.M., King K.S., Donahue C.P., Khrapko K., Kosik K.S. A microRNA array reveals extensive regulation of microRNAs during brain development. RNA (2003) 9:1274–1281.
[Abstract/Free Full Text] -
Kim J., Krichevsky A., Grad Y., Hayes G.D., Kosik K.S., Church G.M., Ruvkun G. Identification of many microRNAs that copurify with polyribosomes in mammalian neurons. Proc. Natl. Acad. Sci. USA (2004) 101:360–365. Epub December 22, 2003.
[Abstract/Free Full Text] -
Dostie J., Mourelatos Z., Yang M., Sharma A., Dreyfuss G. Numerous microRNPs in neuronal cells containing novel microRNAs. RNA (2003) 9:180–186.
[Abstract/Free Full Text] - Ciafre S.A., Galardi S., Mangiola A., Ferracin M., Liu C.G., Sabatino G., Negrini M., Maira G., Croce C.M., Farace M.G. Extensive modulation of a set of microRNAs in primary glioblastoma. Biochem. Biophys. Res. Commun. (2005) 334:1351–1358.[CrossRef][ISI][Medline]
-
Chen C.Z., Li L., Lodish H.F., Bartel D.P. MicroRNAs modulate hematopoietic lineage differentiation. Science (2004) 303:83–86. Epub December 4, 2003.
[Abstract/Free Full Text] - Garzon R., Pichiorri F., Palumbo T., Visentini M., Aqeilan R., Cimmino A., Wang H., Sun H., Volinia S., Alder H., et al. MicroRNA gene expression during retinoic acid-induced differentiation of human acute promyelocytic leukemia. Oncogene (2007) 26:4148–4157. Epub January 29, 2007.[CrossRef][ISI][Medline]
-
Pallante P., Visone R., Ferracin M., Ferraro A., Berlingieri M.T., Troncone G., Chiappetta G., Liu C.G., Santoro M., Negrini M., et al. MicroRNA deregulation in human thyroid papillary carcinomas. Endocr. Relat. Cancer (2006) 13:497–508.
[Abstract/Free Full Text] - Carroll R.C., Beattie E.C., von Zastrow M., Malenka R.C. Role of AMPA receptor endocytosis in synaptic plasticity. Nat. Rev. Neurosci. (2001) 2:315–324.[CrossRef][ISI][Medline]
- Passafaro M., Nakagawa T., Sala C., Sheng M. Induction of dendritic spines by an extracellular domain of AMPA receptor subunit GluR2. Nature (2003) 424:677–681.[CrossRef][Medline]
- Eastwood S.L., McDonald B., Burnet P.W., Beckwith J.P., Kerwin R.W., Harrison P.J. Decreased expression of mRNAs encoding non-NMDA glutamate receptors GluR1 and GluR2 in medial temporal lobe neurons in schizophrenia. Brain Res. Mol. Brain Res. (1995) 29:211–223.[Medline]
- Eastwood S.L., Kerwin R.W., Harrison P.J. Immunoautoradiographic evidence for a loss of alpha-amino-3-hydroxy-5-methyl-4-isoxazole propionate-preferring non-N-methyl-D-aspartate glutamate receptors within the medial temporal lobe in schizophrenia. Biol. Psychiatry (1997) 41:636–643.[CrossRef][ISI][Medline]
- Kajimoto Y., Shirakawa O., Kuno T., Nishino N., Nakai H. Delayed changes in neural visinin-like calcium-binding protein gene expression caused by acute phencyclidine administration. J. Neural Transm. Gen. Sect. (1995) 100:257–262.[CrossRef][ISI][Medline]
- Bernstein H.G., Becker A., Keilhoff G., Spilker C., Gorczyca W.A., Braunewell K.H., Grecksch G. Brain region-specific changes in the expression of calcium sensor proteins after repeated applications of ketamine to rats. Neurosci. Lett. (2003) 339:95–98.[CrossRef][ISI][Medline]
-
Spilker C., Dresbach T., Braunewell K.H. Reversible translocation and activity-dependent localization of the calcium-myristoyl switch protein VILIP-1 to different membrane compartments in living hippocampal neurons. J. Neurosci. (2002) 22:7331–7339.
[Abstract/Free Full Text] -
Mathisen P.M., Johnson J.M., Kawczak J.A., Tuohy V.K. Visinin-like protein (VILIP) is a neuron-specific calcium-dependent double-stranded RNA-binding protein. J. Biol. Chem. (1999) 274:31571–31576.
[Abstract/Free Full Text] - Hansen T., Olsen L., Lindow M., Jakobsen K.D., Ullum H., Jonsson E., Andreassen O.A., Djurovic S., Melle I., Agartz I., et al. Brain expressed micrornas implicated in schizophrenia etiology. PLoS ONE (2007) 2:e873.[CrossRef]
- Burmistrova O.A., Goltsov A.Y., Abramova L.I., Kaleda V.G., Orlova V.A., Rogaev E.I. MicroRNA in schizophrenia: genetic and expression analysis of miR-130b (22q11). Biochemistry (Mosc). (2007) 72:578–582.[CrossRef][Medline]
- Igloi G.L. Nonradioactive labeling of RNA. Anal. Biochem. (1996) 233:124–129.[CrossRef][ISI]




