Decreased cortical muscarinic M1 receptors in schizophrenia are associated with changes in gene promoter methylation, mRNA and gene targeting microRNA.

Many studies have shown decreased cortical muscarinic M1 receptors (CHRM1) in schizophrenia (Sz), with one study showing Sz can be separated into two populations based on a marked loss of CHRM1 (-75%) in -25% of people (Def-Sz) with the disorder. To better understand the mechanism contributing to the loss of CHRM1 in Def-Sz, we measured specific markers of gene expression in the cortex of people with Sz as a whole, people differentiated into Def-Sz and people with Sz that do not have a deficit in cortical CHRM1 (Non-Def-Sz) and health controls. We now report that cortical CHRM1 gene promoter methylation and CHRM1 mRNA are decrease in Sz, Def-Sz and Non-Def-Sz but levels of the micro RNA (miR)-107, a CHRM1 targeting miR, are increased only in Def-Sz. We also report in vitro data strongly supporting the notion that miR-107 levels regulate CHRM1 expression. These data suggest there is a reversal of the expected inverse relationship between gene promoter methylation and CHRM1 mRNA in people with Sz and that a breakdown in gene promoter methylation control of CHRM1 expression is contributing to the global pathophysiology of the syndrome. In addition, our data argues that increased levels of at least one miR, miR-107, is contributing to the marked loss of cortical CHRM1 in Def-Sz and this may be a differentiating pathophysiology. These latter data continue to support the hypothesis that microRNAs (miRNA) have a role in the underlying neurobiology of Sz but argue they are differentially affected in subsets of people within that syndrome.

Introduction

It has been argued that being able to dissect the syndrome of schizophrenia (Sz) into its component disorders would be a significant advance towards understand the aetiologies of what should be more homogeneous component disorders.[1] Significantly, a growing body of data from genetic, neuroimaging and postmortem central nervous system (CNS) studies suggest that cortical muscarinic receptors (CHRMs) are involved in the pathophysiology of Sz.[2] More recently, we have identified a subgroup of people with Sz (∼25%) that can be separated into a discrete population because they have had a marked loss (∼75%) of cortical CHRM (Def-Sz).[3] Thus, for the first time, it has become possible to study Sz as a whole and then to subsequently interrogate data divided into those with Def-Sz and those people with the disorder who do not have the deficits in CHRM (Non-Def-Sz).

In perusing the discovery of Def-Sz it is important to acknowledge that the radioligand binding assay used to delineate the two groups of people with Sz is highly selective for the CHRM1.[4, 5] In addition, it has been reported that CHRM1,[6] but not CHRM2, CHRM3 or CHRM4,[6, 7] protein is decreased in the cortex of people with the disorder. From these data we hypothesised that people with Def-Sz have a marked loss of steady-state cortical CHRM1 and that such a loss of receptor must reflect a change in either gene expression or protein degradation. To begin to address this issue, we decided to determine if critical markers of gene expression (CHRM1 gene promoter methylation, CHRM1 mRNA, CHRM1 targeting microRNA (miRNA)) were altered in the cortex of people with Sz compared with age/sex matched controls and whether these changes were more marked in Def-Sz.

Materials and methods

Postmortem tissue collection and appraisal

Approval for this study was obtained from the Ethics Committee of the Victorian Institute of Forensic Medicine and the Mental Health Research and Ethics Committee of Melbourne Health. All tissue was collected following permission from the next of kin and, where a psychiatric history was identified, an extensive case history review was conducted using a structured diagnostic instrument (Diagnostic Instrument for Brain Studies[8]) allowing a consensus diagnosis according to DSM-IV criteria.[9] Comprehensive medication histories were obtained, the most recently prescribed doses of antipsychotic drugs were converted to a standardized drug dose,[10] those of anticholinergic drugs converted to benztropine equivalents and benzodiazepine doses converted to diazepam equivalents (Table 1). Duration of illness (time from first hospital admission to death) was also calculated.

Table 1

Donor demographics, treatment and CNS collection data for cases from which tissue was collected from BA 9 for this study

Cohort groups	Diagnoses	n	Sex (M/F)	Sui (Y/N)	Age (years)	pH	PMI (hour)	Brain weight (gms)	DOI (years)	APD	ACD	DZD
A	Control	63	47/16		43±2	6.33±0.02	41±1.8	1418±19
	Sz	69	51/18		44±2	6.31±0.02	40±1.5	1427±20	18±2	624±75	1.7±0.12	13±2.8
	Def-Sz	20	16/4	7/13	46±4	6.25±0.04	39±2.7	1461±45	20±4	593±124	1.8±0.09	3.5±0.6
	Non-Def-Sz	49	35/14	26/23	46±2	6.12±0.05	42±2.0	1413±23	18±2	638±94	1.5±0.18	19±4.0
B	Control	15	12/3		45±4	6.30±0.05	39±3.8	1355±38
	Sz	27	21/6		45±3	6.19±0.05	40±2.8	1393±34	18±3	602±109	1.8±0.10	5.7±1.8
	Def-Sz	14	11/3	5/9	44±5	6.25±0.06	40±3.5	1425±55	18±4	623±147	1.8±0.12	3.1±0.6
	Non-Def Sz	13	10/3	3/10	46±5	6.13±0.08	40±4.5	1355±36	18±4	568±169	1.7±0.20	12±4.4

Cohort definitions: Cohort Group A: used for the determination of CHRM1 mRNA and methylation of the promoter region.

Cohort Group B: a subgroup selected at random from cohort A used for studies on miRNAs.

Abbreviations: ACD, last recorded anticholinergic drug doses in benztropine equivalents; APD, last recorded antipsychotic drug doses in chlorpromazine equivalent; Def-Sz, [3H]pirenzepine binding deficient schizophrenia; DOI, duration of illness; DZD, last recorded benzodiazepine doses in diazepam equivalents; F, female; M, male; Non-Def-Sz, subjects with schizophrenia without deficits in [3H]pirenzepine binding; N, no; PMI, postmortem interval; Sui, suicide; Sz, schizophrenia; Y, yes.

Cadavers were refrigerated within 5 h of being found and CNS tissue was rapidly processed according to a standardised protocol,[11] ensuring all CNS tissue was frozen to −70 °C within 30 min of being removed at autopsy. Consistent refrigeration of cadavers within 5 h was important because data suggests such rapid refrigeration would significantly inhibit autolysis.[12] To determine the quality of tissue preservation, the pH of CNS tissue was determined as described previously.[13] For witnessed deaths, postmortem interval (PMI) was the time between death and autopsy but when death was not witnessed PMI was taken as the interval mid-way between the donor last being seen alive and being found dead.

Tissue processing

For all components of this study tissue was taken from Brodmann's Area 9 (BA 9), which was defined as the lateral surface of the frontal lobe, including the middle frontal gyrus superior to the inferior frontal sulcus, of the left CNS hemisphere.

CHRM1 gene methylation and mRNA levels

Measurement of DNA methylation was completed using tissue from 132 donors; 69 people meeting the DSM-IV criteria for the diagnoses of Sz and 63 people with no history of psychiatric illness (controls) (Table 1, cohort A), these cases were all included in our previous study that identified a subset of subjects with Def-Sz.[3] The 69 people with Sz consisted of 49 people with Sz who did not have a deficit in [3H]pirenzepine binding (Non-Def-Sz) and 20 people with Sz who had Def-Sz.[3]

Genomic DNA extraction and bisulphite conversion

To measure CHRM1 promoter methylation, BA 9 was homogenised in DNA extraction buffer (100 mℳ NaCl, 10 mℳ Tris-HCl (pH 8.0), 25 mℳ EDTA, 0.5% SDS), containing 200 μg ml−1 Proteinase K and incubated for 3 h after which DNA was isolated using a standard phenol: chloroform method. The DNA was precipitated with absolute ethanol containing 10% sodium acetate (3 ℳ, pH 5.2) and the DNA pellets were washed with 70% ethanol and resuspended in 100 μl TE buffer (10 mℳ Tris-HCl (pH7.5), 1 mℳ EDTA). DNA was quantitated and its quality assessed spectroscopically using a Nanodrop (ThermoFisher Scientific, Walthan, MA, USA).

For each case, 2 μg of genomic DNA was bisulphite converted using the MethylEasyXceed Rapid bisulphite modification kit (Human Genetic Signatures, Sydney, NSW, Australia) following the manufacturer's instructions. Using this approach sodium bisulphite selectively converts unmethylated cytosine residues to uracil nucleotides, leaving methylated cytosine unchanged. All samples were processed in parallel to ensure consistency.

PCR amplification

Primers for CHRM1 were designed using PrimerExpress 3.0 following sequence bisulphite conversion using Methyl Primer Express v1.0 software. (Supplementary Table 1) for primer sequences. Three amplicons were designed, subsequent to a bioinformatic analysis, to characterise DNA methylation level and distribution at the CHRM1 gene promoter (See Supplementary Figure 1). To reduce technical variability, three replicate amplifications were performed in 15 μl reactions containing 10 ng bisulphite-converted DNA, 1X FastStart PCR master mix (Roche Diagnostics Australia Pty Castle Hill, NSW, Australia) and 3 nℳ each of forward and reverse primers. The amplification conditions were: initial denaturation at 95 °C for 2 min, followed by five cycles of 95 °C for 30 s, 56 °C for 30 s, 72 °C for 90 s, then 35 cycles of 95 °C for 30 s, 60 °C for 30 s, 72 °C for 90 s, and a final extension at 72 °C for 7 min.

DNA methylation

DNA methylation levels were determined as described previously[14] using MassARRAY EpiTYPER (Sequenom, QLD, San Diego, CA, USA), based on MALDI-TOF mass spectrometry. In short, bisulphite converted DNA was amplified with primers containing a T7-promoter tag and the amplification products cleaned with shrimp alkaline phosphatase (SAP, Sequenom, San Diego, CA, USA). After purification, in vitro transcription and T specific cleavage was performed. Each reaction was then spotted onto a Maldimatrix-containing SpectroCHIPS (Sequenom) and subjected to MALDI-TOF MS. The mass spectra was collected by MassARRAY Spectrometer and analysed by EpiTYPER v.1.0 software (Sequenom). Data was obtained for each CpG-containing fragment (unit), with most units containing single CpGs, and two units (0304_CpG5_6 and 06_CpG3_4) containing two CpGs. For simplicity, all CpG units are referred to as ‘sites' in the Results.

The Sequenom technology, like all sequencing technologies,[15] gives results on single nucleotides that is of varying certainty and there are issues relating to the difficulties in measuring the status of every potential nucleotide methylation status.[16] Therefore to ensure our data analyses only included high quality data we applied a stringent quality control process to remove potentially unreliable measurements before analysis.[14] Firstly, all data from CpG-containing fragments flagged by EpiTYPER as having low or high mass (outside the MS analytical window), or overlapping MS peaks were discarded. Secondly, CpG analytic units that failed to produce data for >30% of samples were also discarded (unreliable CpG unit), and samples with more than 30% missing data points within an amplicon (unreliable sample) had all methylation values for that sample set to ‘missing'. Finally, technical replicates showing ⩾5% absolute difference from the median value of a set of technical replicates were set to ‘missing' and only samples with at least two successful technical replicates were analysed.

Quantifying CHRM1 mRNA

In situ hybridisation

In situ hybridisation was performed as described previously.[17] Briefly, 6 × 10 μm sections were cut from the BA 9 of each subject, fixed in 4% paraformaldehyde, placed in fresh 0.25% acetic anhydride in 0.1 ℳ triethanolamine HCl (pH 8.0), dehydrated, delipidated and equilibrated in 95% ethanol before being dried.

To increase specificity, the oligonucleotide probe contained a mix of three antisense sequences, the probes were complementary to bases 4–51, 721–768 and 811–853 of human CHRM1 mRNA (accession no.: X52068).[6] The oligonucleotides were labelled with [35S]ATP using the NEN Oligonucleotide 3′ end labelling system and purified using a NENSORB 20 purification cartridge (Perkin Elmer, Glen Waveryly, VIC, Australia).

For hybridisation, the probe was prepared at 1 × 106 dpm of [35S]ATP per 50μl of hybridisation buffer (formamide (50%), NaCl (600 mℳ, Tris-HCl (80 mℳ), ethylenediamine tetraacetate (4 mℳ, sodium pyrophosphate (0.1%), SDS (0.2%), sodium heparin (2 mg ml−1), dextran sulphate (10%) and dithiothreitol (100 mℳ)) and incubated with tissue sections for 20 h at 37 °C, in a humidified chamber. Two sections were incubated with the radioactive oligonucleotide probe (total binding) while two sections were incubated with radioactive probes in the presence of excess non-radioactive probes (non-specific binding). Hybridisation to mRNA was confirmed by incubating two sections, preincubated in 0.02% RNaseA, with the radioactive probe (negative control).

After hybridisation, all sections were washed[18] and dried before being apposed to BAS-SR2025 imaging plates (Fujifilm, Tokyo, Japan), with [14C]microscales (GE Healthcare, Rydalmere, NSW, Australia) for 4 weeks before being scanned in the BAS 5000 high resolution phosphoimager (Fujifilm).[19] The signal intensity of the phosphoimages were measured against those on the [14C]microscales using AIS image analysis software. Results were expressed as TB minus NSB in dpm × 103 per mg estimated wet weight tissue equivalents.

miRNA

This component of the study used Brodmann's Area 9 from 14 people with Def-Sz, 13 people with Non-Def Sz and 15 control people (Table 1B).

CHRM1 miRNA target site analysis

CHRM1 target site prediction was accomplished using the miRNA target database miRGen (http://diana.pcbi.upenn.edu/cgi-bin/miRGen/v3/Targets.cgi) (Supplementary Table 2). This service provides meta-analysis of established miRNA target prediction algorithms that search the 3′ UTR sequences for homology and putative affinity with known miRNA sequences.[20] These algorithms also have weightings that enable the appreciation of conservation between 3′ UTR motifs in vertebrate species that are more likely to have functional significance. An important feature of the miRGen database is that it provides output on the amalgamation of predictions from multiple algorithms to yield a broad low stringency appraisal of a given target gene or more stringent outputs that consider the intersection of miRNA target sites that satisfy multiple algorithms including miRanda, PictaR, TargetScanS and DIANA-microT.

RNA purification and first-strand complementary DNA synthesis

Total RNA was isolated from 100 mg frozen cortical samples with 1.0 ml TRIzol reagent (Invitrogen, Mulgrave, VIC, Australia) for human tissue and QIAGEN (Doncaster, VIC, Australia) RNA extraction kit for mouse tissue, according to the manufacturer's instructions. Following Trizol extraction, the RNA was treated with 0.1 U μl−1 DNAse in 1 × DNAse buffer (Invitrogen, Australia) to remove genomic DNA contamination, extracted using phenol (pH 4.2): chloroform extraction and ethanol precipitated. complementary DNA was synthesised as described previously.[21]

miRNA expression

miRNA expression analysis was performed using miRNA specific relative quantitative real-time PCR as described previously.[22] Primer sequences used in reverse transcription reaction and quantitative real-time PCR are provided in Supplementary Table 3. Relative expression was determined with respect to the geometric means of three constitutively expressed small nuclear RNA (U6, U44 and U49) by subtracting the geometric means of their cycle threshold value from the cycle threshold value derived for the miRNA.[23] All reactions were subjected to post amplification dissociation curve analysis to provide guidance on reaction specificity. Initial confirmation of amplicon length was also accomplished through gel electrophoresis.

Reporter gene assay

Luciferase reporter gene analysis of miRNA function was conducted as described previously.[24] Briefly, HEK-293 cell cultures were seeded into 24-well plates and maintained in DMEM with 10% (vol/vol) foetal calf serum, 20 mℳ HEPES, 0.15% (wt/vol) sodium bicarbonate and 2 mℳ L-glutamine at 37 °C with 5% CO2 and 90% humidity. Transfections were performed 24 h postseeding with recombinant CHRM1-Firefly luciferase reporter gene constructs, 100 nℳ synthetic miRNA or LNA modified anti-miR oligonucleotide in Lipofectamine 2000 according to manufacturer's instructions (Invitrogen). Reporter gene constructs containing 3′ UTR sequences encoding target gene miRNA recognition elements were produced through the ligation of a synthetic double stranded cassette compatible with SpeI and Hind III digested pMIR-REPORT vector (Invitrogen) (see Supplementary Figure 2). Reporter gene silencing in response to miRNA cotransfection was monitored with respect to a control plasmid expressing Renilla luciferase (pRL-TK) using the dual luciferase reporter assay (Promega, Alexandria, NSW, Australia). Non-specific effects associated with transfection were controlled by comparison to cells cotransfected with mutant miRNAs or mutant anti-miRs.

[3H]pirenzepine binding

The data on [3H]pirenzepine binding used in this study was originally generated in an earlier study.[3] The data pertaining to the relevant cohorts was collated and then reanalysed within the parameters of each portion of this study.

Statistics

The distribution of all data was determined using the D'Agostino and Pearson omnibus normality test as this is best for determining data distribution in small cohorts.[25] Demographic, pharmacological and CNS collection data were compared using a student's t-test at the levels of Sz and a one-way analysis of variance when comparing Def-Sz and Non-Def-Sz to controls. Where data were found to be mixed between parametric and non-parametric distributions the General Linear Model was used to identify significant variance whereas Mann–Whitney U or Kruskal–Wallis tests were used when all data to be analysed was shown to have a non-parametric distribution. Gender frequency across diagnostic cohorts was assessed using the χ test. Potential biological relationships and confounding factors were integrated using linear regression analysis. Due to small cohort sizes only strong relationships (r=0.49) where the regression deviated significantly were considered as being of interest.[26]

Most analyses were conducted using Prism 5.01 (Graphpad Software, La Jolla, CA, USA) but the General Linear Model was computed using Minitab 15 (Minitab, State College, PA, USA) while the methylation data was analysed using Stata 11 (StataCorp., 2010; College Station, TX, USA).

Results and discussion

Demographics, pharmacological history and CNS collection data

mRNA and methylation cohorts:

The cases used in this component of the study were part of the cohorts used in the study that identified Def-Sz[27] with the cohort sizes being limited by availability of tissue from BA 9. There were no significant differences in mean age (P=0.63; P=0.92), PMI (P=0.75; P=0.84), CNS pH (P=0.19; P=0.10), CNS weight (P=0.54, P=0.23) or gender frequency (P=1.00; P=0.11) when comparing Sz to control or Def-Sz, Non-Def-Sz and controls, respectively, for the larger cohorts (Table 1: Cohort Groups A). In comparing Def-Sz and Non-Def-Sz there were no significant differences in mean duration of illness (P=0.87), the last recorded dose of antipsychotic (P=0.73) or anticholinergic (P=0.15) drugs but the Def-Sz had lower final recorded doses of benzodiazepines (P=0.01). The frequency of benzodiazepine prescription did not differ between Def-Sz and Non-Def-Sz (P=0.60). This would suggest that the symptoms of schizophrenia targeted by benzodiazepines (distress, insomnia and behavioural disturbances secondary to psychosis[28]) are more responsive to benzodiazepine treatment in subjects with Def-Sz by a mechanism that has yet to be understood.

miRNA cohort:

The cohorts for this study were selected based on diagnoses, whether or not they were shown to be Def-Sz or Non-Def-Sz by measuring [3H]pirenzepine binding, the availability of tissue from BA 9 and limitations in sample processing capacity of the technique. Together these factors contributed to resulted in smaller cohorts being studied when measuring miRNA. There were no significant differences in mean age (P=0.55; P=0.83), PMI (P=0.71; P=0.93), CNS pH (P=0.11; P=0.16), CNS weight (P=0.44; P=0.29) or gender frequency (P=1.00; P=0.99)when comparing Sz to control or Def-Sz, Non-Def-Sz and controls, respectively (Table 1: Cohort Groups B). In comparing Def-Sz and Non-Def-Sz there were no significant differences in mean duration of illness (P=0.96), the last recorded dose of antipsychotic drugs (P=0.59) or anticholinergic drugs (P=0.50) but the Def-Sz had lower final recorded doses of benzodiazepines (P=0.04).

CHRM1 Gene methylation

Ten potential CpG sites were identified across the putative CHRM1 promoter region (Supplementary Figure 1). There was a significant reduction in levels of methylation in four of these sites in people with Sz (Table 2) with further analyses revealing a significant reduction in levels of DNA methylation in three sites in people with Def-Sz and Non-Def-Sz (Table 3), plus an addition site in the Non-Def-Sz. Examining the pooled methylation data showed that there was a significant decrease in methylation associated with the diagnoses of Sz but not with other factors such as suicide, gender, age, tissue pH or PMI (Table 4). In addition, there were very similar reductions in DNA methylation in all CpG sites in tissue from both Def-Sz and Non-Def-Sz compared with controls, which meant there was no difference in DNA methylation between the two forms of Sz.

Table 2

Difference in mean CHRM1 DNA methylation at individual CpG units for schizophrenia vs control groups (after adjustment for death by suicide, sex, age, pH and PMI)

CpG unit	Genome coordinate	Number		Methylation		Difference (Sz minus Controls)
		Sz	Control	Sz	Control	Δ Mean Methylation	P	95% CI
chrm1_0401_CpG2	62445428	56	54	16.0%	17.9%	−2.9%	0.17	−7.0%	1.2%
chrm1_0401_CpG3	62445440	64	60	7.3%	7.4%	−0.4%	0.55	−1.8%	1.0%
chrm1_0304_CpG2	62445521	55	49	15.8%	17.0%	−2.2%	0.072	−4.5%	0.2%
chrm1_0304_CpG4	62445617	48	44	25.6%	28.0%	−5.3%	0.006	−9.1%	−1.5%
chrm1_0304_CpG5_6	62445641	54	52	18.7%	20.8%	−3.2%	0.017	−5.9%	−0.6%
chrm1_06_CpG2	62445645	58	55	15.1%	16.4%	−2.7%	0.010	−4.8%	−0.6%
chrm1_06_CpG3_4	62445657	63	60	15.2%	15.7%	−0.9%	0.32	−2.5%	0.8%
chrm1_06_CpG7	62445700	61	56	26.7%	26.8%	−1.1%	0.34	−3.4%	1.2%
chrm1_06_CpG8	62445953	58	57	32.8%	34.5%	−2.4%	0.046	−4.9%	0.0%
chrm1_06_CpG9	62445992	59	55	41.7%	44.1%	−2.6%	0.075	−5.4%	0.3%

Abbreviations: CI, confidence interval; Sz, schizophrenia.

Bold entries indicate P<0.05.

Table 3

Difference in mean methylation for Def and Non-Def schizophrenia groups minus controls at individual CpG units (after adjustment for sex, age pH, PMI and death by suicide)

CpG unit	Genome coordinate	Number			Def-Sz vs controls				Non-Def-Sz vs controls
		Sz		Control	Δ Mean methylation	P	95% CI		Δ Mean methylation	P	95% CI
		Def	Non-Def
chrm1_0401_CpG2	62445428	18	38	54	−3.9%	0.12	−8.8%	1.0%	−2.2%	0.34	−6.8%	2.4%
chrm1_0401_CpG3	62445440	20	44	60	−0.7%	0.40	−2.4%	1.0%	−0.2%	0.78	−1.8%	1.3%
chrm1_0304_CpG2	62445521	16	39	49	−2.5%	0.08	−5.4%	0.3%	−1.9%	0.16	−4.5%	0.7%
chrm1_0304_CpG4	62445617	14	34	44	−6.8%	<0.01	−11.6%	−2.1%	−4.4%	0.04	−8.6%	−0.3%
chrm1_0304_CpG5_6	62445641	16	38	52	−3.8%	0.02	−7.1%	−0.6%	−2.8%	0.07	−5.7%	0.2%
chrm1_06_CpG2	62445645	17	41	55	−3.1%	0.01	−5.5%	−0.6%	−2.4%	0.04	−4.7%	−0.1%
chrm1_06_CpG3_4	62445657	20	43	60	−0.8%	0.42	−2.9%	1.2%	−0.9%	0.37	−2.8%	1.0%
chrm1_06_CpG7	62445700	20	41	56	−1.0%	0.46	−3.7%	1.7%	−1.2%	0.36	−3.8%	1.4%
chrm1_06_CpG8	62445953	17	41	57	−2.1%	0.16	−5.0%	0.8%	−2.7%	0.05	−5.4%	0.0%
chrm1_06_CpG9	62445992	20	39	55	−2.2%	0.19	−5.5%	1.1%	−2.9%	0.08	−6.1%	0.3%

Abbreviations: CI, confidence interval; chrm1, cortical muscarinic M1 receptors; Def-Sz, [3H]pirenzepine binding deficient schizophrenia; Non-Def-Sz, subjects with schizophrenia without deficits in [3H]pirenzepine binding; Sz, schizophrenia.

Bold entries indicate P<0.05.

Table 4

Predictors of CHRM1 DNA methylation

CHRM1 DNA methylation
Covariate	Coefficient	P	95% CI
Schizophrenia vs controls
Schizophrenia	−2.2%	<0.01	−3.8%	−0.6%
Suicide	1.7%	0.09	−0.3%	3.7%
Sex (males minus females)	−0.5%	0.53	−2.0%	1.0%
Age (z-score)	−0.4%	0.33	−1.1%	0.4%
pH (z-score)	−0.4%	0.36	−1.2%	0.4%
Post mortem interval (z-score)	−0.1%	0.86	−0.8%	0.7%
Def-Sz, non-Def-Sz and controls
non-Def_Sz minus controls	−2.0%	0.03	−3.9%	−0.2%
Def_Sz minus controls	−2.4%	0.02	−4.5%	−0.4%
Def-Sz minus non-Def-Sz	−0.4%	0.70	−2.5%	1.7%

Abbreviations: CI, confidence interval; Def-Sz, [3H]pirenzepine binding deficient schizophrenia; Non-Def-Sz, subjects with schizophrenia without deficits in [3H]pirenzepine binding; Sz, schizophrenia.

Bold entries indicate P<0.05.

In situ hybridisation and in situ radioligand binding

Levels of oligonucleotide hybridisation to CHRM1 mRNA and [3H]pirenzepine binding were homogenous throughout the cortical laminae and therefore a single integrated measure was taken across the cortex. Levels of CHRM1 mRNA (P=0.01) and [3H]pirenzepine binding (P<0.001) were significantly lower in the cortex of people with Sz compared with controls (Figure 1a).

Figure 1

Levels (Mean±s.e.m.) of [3H]pirenzepine binding and oligonucleotide CHRM1 probe hybridisation to Brodmann's area 9 (BA 9) from subjects with. (a): schizophrenia and controls and (b): Def-Sz, Non-Def-Sz and controls.

There was a significant variation in levels of oligonucleotide hybridisation to CHRM1 mRNA with group between people with Def-Sz, Non-Def-Sz and controls (F=9.09, d.f.=2148, P<0.0001) due to lower levels of mRNA in the cortex from people with Def-Sz (P=0.001) and Non-Def-Sz (P=0.004) compared with controls (Figure 1b) but levels of CHRM1 mRNA in Def-Sz were not significantly different to that in Non-Def-Sz (P=0.63). As would be predicted, there was also significant variation in levels of [3H]pirenzepine binding with group (F=95.6, d.f.=2148, P<0.0001) due to decreased levels of binding in Def-Sz (P<0.0001) and Non-Def-Sz (P=0.02) compared with controls.

Initially, to identify potential miRNA regulators of CHRM1 we performed a low stringency search of the miRGen miRNA target prediction database, using the ‘union' of their algorithms.[20] This suggested that CHRM1 could be influenced by up to 32 different miRNA (Supplementary Table 2). However, a more stringent appraisal of targets using the ‘intersection' of different algorithms identified a single miR-107 as the miRNA most likely to target CHRM1 expression. The miR107 target site in the CHRM1 3′ UTR is only one of two that display a highly conserved complete seven nucleotide seed-paring region (Supplementary Figure 2) and, additionally, has previously been shown to be increased in postmortem superior temporal cortex and dorsolateral prefrontal cortex from people with Sz.[23] Similar analyses suggested it was highly unlikely that miR-16 would target CHRM1 expression and so levels of this miR were measured to determine if any change in miR-107 was due to a generalised change in levels of miR in the cortex of people with Sz.

In the cohorts used in the study of miRNA, levels of cortical [3H]pirenzepine binding (P=0.0004) and CHRM1 mRNA (P=0.03) were decreased in BA 9 from people with Sz (Figure 2a) but neither levels of miR-16 (P=0.50) nor miR-107 (P=0.19) varied significantly with diagnosis (Figure 2b).

Figure 2

(a): Levels (Mean±s.e.m.) of [3H]pirenzepine binding and oligonucleotide CHRM1 probe hybridisation to BA 9 from the cohort of subjects with schizophrenia and control subjects used for the study of miRNA and. (b): Levels (Mean±s.e.m.) of miR-16 and miR-107 in the same region from the individuals. (c): Levels (Mean±s.e.m.) of [3H]pirenzepine binding and oligonucleotide CHRM1 probe hybridisation to BA 9 from the cohort of subjects with Def-Sz, Non-Def-Sz and control subjects used for the study of miRNA and (d): Levels (Mean±s.e.m.) of miR-16 and miR-107 in the same region from the individuals.

In comparing data from Def-Sz, Non-Def-Sz and controls, levels of miR-107 (F=4.25, d.f.=2,40, P=0.02) and [3H] pirenzepine binding (F=72.3, d.f.=2,40, P<0.0001) varied with group, levels of CHRM1 mRNA showed a trend to varying with group (F=2.68, d.f.=2,40, P=0.08) whereas levels of miR-16 did not vary with group (F=0.65, d.f.=2,40, P=0.52) (Figures 2c and d). The variance in miR-107 was due to increased levels of that miRNA in the cortex of people with Def-Sz compared with controls (P<0.05) but levels of miR-107 were not significantly different between the Non-Def-Sz and controls. The variation in [3H]pirenzepine binding with group was due to decreased levels of binding in cortex from people with Def-Sz (P<0.001), but not Non-Def-Sz, compared with controls. The trend to variance in levels of CHRM1 mRNA was predominantly due to lower levels in the Def-Sz group (t-test compared with controls P=0.03) rather from the Non-Def-Sz (t-test compared with controls P=0.67).

Potential biological relationships

For the cohorts of people used for the study of [3H]pirenzepine binding, CHRM1 mRNA and gene promoter methylation, there were no strong correlations between levels of [3H]pirenzepine binding and levels of mRNA (Supplementary Table 4). There were no correlations between any levels of CpG methylation and any demographic, brain collection of pharmacological variables (data not shown). Significantly, in the cohorts used for the study of miRNA there were correlations between levels of miR-16 and miR-107, but not any of the other measures in all diagnostic cohorts (Figure 3 and Supplementary Table 5).

Figure 3

The relationship between levels of miR-16 and miR-107 in BA 9 from subjects with schizophrenia, Def-Sz, Non-Def-Sz and controls.

CHRM1 reporter expression is modulated by intracellular miR-107

To further investigate the role of miR-107 in the regulation of CHRM1 we cloned a segment of its 3′ UTR downstream of the Firefly luciferase reporter gene in pMIR-REPORT and analysed reporter gene expression in HEK293 cells cotransfected with synthetic miR-107 or LNA modified antimiR-107 (miR-107 inhibitor). In each case, the Firefly luciferase gene dosage was normalised by Renilla luciferase expressed from cotransfected pRL-TK vector and compared with control transfections containing mutant versions of miR-107 (and the corresponding LNA modified antimiRs). Following transfection, analysis of luciferase activity demonstrated that the addition of miR-107 was capable of reducing CHRM1 reporter gene expression by 20% (P=0.004) (Figure 4). When the endogenous miR-107 was depleted or inhibited by its cognate anti-miRs, the reporter gene expression was increased by 62% (P=0.03). These observations, showing gene silencing in the presence of miR-107 over-expression and a reduction of this effect when miR-107 is inhibited, demonstrate that this region of CHMR1 3' UTR is potently regulated by miR-107, at least in vitro.

Figure 4

Effect (Mean±s.e.m.) of miR-107 and antimiR-107 on the activity of a CHRM1 reported gene construct. (a): P<0.05; (b): P<0.005.

Potential confounds

The only strong and significant relationships between experimental measures and potential numeric confounds was between [3H]pirenzepine binding and final recorded benzodiazepine drug dose (r=0.68, P<0.001) in people with Sz and between miR-107 and PMI (r=0.47, P<0.01) in Non-Def-Sz (Supplementary Table 4). Effects of these relationships on the original analyses using an analysis of covariance showed that final recorded benzodiazepine drug dose as a covariate did not affect the variation of [3H]pirenzepine binding with diagnoses (F=101, d.f.=12 129, P<0.0001). When PMI was included as a covariate in the analyses of miR-107 it is notable that there is an increase in the variance associated with group (F=4.94, d.f.=1,2,39, P=0.01).

There was no effect of either gender (P=0.63) or suicide (P=0.17) on levels of [3H]pirenzepine binding in either people with Sz as a whole or their age/ sex matched controls. Levels of CHRM1 mRNA did not vary with gender or suicide in these cohorts (P=0.41 and P=0.26, respectively). In the cohorts used for the study of miRNA; levels of [3H]pirenzepine binding (P=0.63), CHRM1 mRNA (P=0.87), miR-16 (P=0.78) or miR-107 (P=0.30) did not vary with gender when people with Sz were compared with controls. Similarly, levels of [3H]pirenzepine binding (P=0.46), CHRM1 mRNA (P=0.50), miR-16 (P=0.26) and miR-107 (P=0.13) did not vary with suicide.

This study explores the potential mechanisms that could contribute to widely replicated finding of reduced levels of cortical [3H]pirenzepine binding in people with Sz[3, 6, 29, 30, 31, 32] taking into account that, under the conditions used in most protocols, the radioligand will show a⩾80% selectivity for CHRM1.[4, 5] In addition, it explores whether changes in the status of markers of CHRM1 expression may differentiate people with Def-Sz, who have a ∼75% loss of cortical [3H] pirenzepine binding,[3] from Non-Def-Sz. Consistent with previous findings[6, 33] we show a decrease in levels of cortical CHRM1 mRNA in Sz that is equivocal across Def-Sz and Non-Def-Sz. Our data also shows a decrease in levels of methylation of the putative promoter region of the CHRM1 gene in Sz. There was a significant decrease in the levels of methylation in four CpG sites in Non-Def-Sz whereas only three of these sites had significant decrease in levels of gene methylation in Def-Sz. This being noted, our overall analyses suggested differences in CHRM1 gene promoter methylation did not permit separation of the two groups. The presence of low levels of CHRM1 binding and mRNA has a symmetry but it is generally argued that increased gene methylation is associated with gene silencing and lower levels of mRNA.[34] Thus, on the one hand our data is consistent with the suggestion that Sz is associated with a hypo-methylation state[35] but, on the other hand, our data also suggests such a state may not be associated with the expected increase in gene expression in the CNS of people with the disorder. Thus, our data leads us to postulate that the pathophysiology of Sz may interfere with the relationship between gene promoter methylation and levels of the gene expression.

Our study showed levels of miR-107, but not miR-16, were increased in the cortex from people with Def-Sz, but not Non-Def-Sz. Significantly, from sequence homology, miR-107 was predicted to target CHRM1 expression whereas miR-16 would not target CHRM1 mRNA. The notion that miR-107 targets CHRM1 expression is supported by our data showing changing miR-107 levels appropriately changes the activity of a reporter gene linked to levels of CHRM1 expression. Therefore our data suggests that CHRM1 gene expression in people with Def-Sz is affected by the increased activity of at least one CHRM1 targeting miRNA and this unique pathophysiology could be causing the large decrease in cortical CHRM1 in these people. Overall, our data suggests that the changes in the regulation of CHRM1 expression is not a straightforward progression whereby changes in levels of mRNA predict changes in levels of protein; rather our data is in keeping with the concept of a derangement of the highly complex system of checks and balances regulating expression in the CNS.[36]

This study, for the first time, to our knowledge, offers an explanation as to an underlying mechanism that may be contributing to the marked decrease in cortical CHRM1 levels in Def-Sz. Our data suggests this mechanism involves increased levels of at least one miR, miR-107, in the cortex of people with Def-Sz and our in vitro data shows this increase in miR-107 would act to decrease CHRM1 expression and/or translation. Such a functional link is supported by our data showing a negative correlation between [3H]pirenzepine binding and levels of miR-107 (r=−0.3671, P=0.001) in the cortex of people with Def-Sz that was not detectable in tissue from Non-Def-Sz and controls. We interpret these data as suggesting a stronger relationship between miR-107 and levels of CHRM1 in Def-Sz that would be expected if miR-107 was a significant factor in its pathophysiology. Importantly, levels of a miR-16, that does not target CHRM1 expression, is not altered in the cortex of Def-Sz and Non-Def-Sz showing the changes in miR-107 are not reflecting generalised changes in miR in the cortex of people with Def-Sz. Coalescing our data on CHRM1 binding, CHRM1 mRNA and miR-107 we now hypothesise that the marked decrease in CHRM1 levels in Def-Sz is due to the cumulative effects of decrease gene expression being amplified by the effects of increased CHRM1 targeting miRNA activity.

Our data also adds to a growing body of data implicating gene methylation in the pathophysiology of Sz.[37, 38, 39, 40, 41] However, our data differs in that it shows a decrease in CHRM1 expression (as measured by levels of mRNA) in the presence of decreased gene methylation, which would normally be expected to cause a decrease in gene silencing (that is, increased expression).[34] Significantly, the expected inverse relationship between methylation and levels of gene expression has been reported in the face of decreased gene copy number in transgenic animals.[42] Levels of gene methylation have long been known to be subject to feed-back loop regulation and it has been suggested that disruption of such processes could result in disease states.[43] Thus, our data may reflect a disruptive feed-back loop involving methylation of CHRM1 as part of the pathophysiology of Sz because CNS gene methylation is a dynamic process that varies during life[44] and the decreased methylation of the CHRM1 promoter region we report in Sz could be part of a feedback mechanism trying to turn on the expression of CHRM1 in an attempt to compensate for low levels of cortical mRNA and/or protein.

In conclusion, our study has identified changes in multiple markers (gene promoter methylation, levels of mRNA and protein) of CHRM1 gene expression in the dorsolateral prefrontal cortex from people with Sz. However, being able to analyse our data at the subsyndrome level has revealed an increase in a CHRM1 mRNA targeting miR (miR-107) only in the cortex of people with Def-Sz. We therefore hypothesise that this extra burden on CHRM1 gene expression is a significant contributor to the marked loss in cortical CHRM1 by which this subset of people is defined. Our data also supports that of others that have shown changes in levels of miRNA in the CNS of subjects with schizophrenia the suggest a role for miRNAs in the pathophysiology of Sz[23] but extends this hypothesis to suggest that changes in miRNA levels may differentiate some subsets of people within the syndrome. However, we acknowledge that our sample size for this portion of our study is small and therefore our results must be regarded with caution until further data is obtained.

Our findings on the mechanisms underpinning changes in CHRM1 in Sz comes at a critical time given the development of CHRM1 allosteric modulates which, for the first time, open up the possibility of specifically targeting that receptor to bring about therapeutic benefits.[45] Significantly, data on these new compounds support the hypothesis that they will be useful in reducing the symptoms of Sz[46] but our data may indicate they may be useful in treating selective subgroups of people with the disorder. Hence, a high priority should be given to developing more selective CHRM1 ligands for neuroimaging, rather than the pan-muscarinic receptor ligand used to show decreased CHRMs in living people with Sz,[47] to facilitate the testing of these new drugs in subjects who have a marked loss of CHRM1 compared to those who have not. Finally, our data strongly supports the argument that subdividing the syndrome of Sz will facilitate an improved understanding of its component disorders[1] and that this understanding may lead to more personalised medicines with which to treat these component disorders.[48]