DNA methylation in schizophrenia: progress and challenges of epigenetic studies.

Schizophrenia is a severe psychiatric disease affecting about 1% of the world's population, with significant effects on patients and society. Genetic studies have identified several candidate risk genes or genomic regions for schizophrenia, and epidemiological studies have revealed several environmental risk factors. However, the etiology of schizophrenia still remains largely unknown. Epigenetic mechanisms such as DNA methylation and histone modifications can explain the interaction between genetic and environmental factors at the molecular level, and accumulating evidence suggests that such epigenetic alterations are involved in the pathophysiology of schizophrenia. However, replication studies to validate previous findings and investigations of the causality of epigenetic alterations in schizophrenia are needed. Here, we review epigenetic studies of schizophrenia patients using postmortem brains or peripheral tissues, focusing mainly on DNA methylation. We also highlight the recent progress and challenges in characterizing the potentially complex and dynamic patterns of epigenomic variations. Such studies are expected to contribute to our understanding of schizophrenia etiology and should provide novel opportunities for the development of therapeutic drugs.

Schizophrenia: clinical course, epidemiology and pathogenesis

Schizophrenia is a severe and relatively common psychiatric disease, affecting approximately 1% of the population worldwide [1]. Patients with schizophrenia suffer from three categories of symptoms: positive symptoms (for example, auditory hallucinations and paranoia), negative symptoms (for example, anhedonia and social withdrawal) and cognitive impairment. The disease course typically starts with positive symptoms, and then progresses to a chronic stage with negative symptoms and cognitive impairment. Many patients experience the onset of schizophrenia around adolescence and usually need medication for the rest of their lives. The pattern of symptoms differs greatly between individuals, making the characterization of schizophrenia difficult. Consequently, the societal burden of schizophrenia is also significant. The World Health Organization estimated that schizophrenia contributed about 2% of worldwide disability-adjusted life-years [2]. Although various kinds of antipsychotic medications, such as clozapine and haloperidol, and other psychotropic medications, such as valproate (VPA), as well as social services are useful for the alleviation of symptoms, there are no curative therapies. In addition, antipsychotic drugs are supposed to be effective for the treatment of positive symptoms but not for the treatment of other symptoms. Therefore, it is imperative to decipher the etiology of schizophrenia to establish better therapeutic strategies [1].

Schizophrenia is partly a genetic disease, as evidenced by family, twin and adoption studies [3]. Traditionally, the risk of schizophrenia has been estimated to be up to 46% when both parents have schizophrenia. The concordance rate of monozygotic twins for schizophrenia is 48% and is always higher than that of dizygotic twins [3]. The fact that not all monozygotic twins are concordant for the onset of schizophrenia suggests that environmental factors are also important. For example, viral infection during the fetal period is a well-known high risk factor for schizophrenia [4]. Greater parental ages at the time of birth [5], reduced nutrition during the fetal period and after birth [6], perinatal abnormalities [7], social stress [8] and other factors [9-12] have also been reported as risk factors (Table 1). Those with a high genetic risk for schizophrenia are thought to be more vulnerable to the environmental factors, a phenomenon known as gene-environment interaction (G × E). For example, people with a family history of psychosis are more likely to suffer from schizophrenia if they have experienced viral infection during the fetal period [13].

Table 1

Examples of environmental factors implicated in schizophrenia

Environmental risk factors	Reference(s)
Viral infection during the fetal period	[4,13]
Greater parental ages at the time of birth	[5]
Reduced nutrition during the fetal period and after birth	[6]
Perinatal abnormalities	[7]
Social stress	[8]
Living in urban areas	[9]
Immigration	[10]
Winter/spring birth	[11]
Cannabis use	[12]

Numerous family, linkage and case-control genetic association studies have been performed. Genes such as COMT, DTNBP1, NRG1, RGS4, GRM3, G72, PPP3CC, CHRNA7 and PRODH have been relatively well replicated as genetic factors for schizophrenia [14-18]. Large-scale genome-wide association studies have also identified replicable genetic risk factors, including the major histocompatibility complex (MHC) region, and genes such as ZNF804A, TCF4, NRGN and MIR137 [19-23]. These main implicated genes and loci are summarized in Table 2.

Table 2

Examples of genetic factors implicated in schizophrenia

Chromosomal location	Candidate gene	Encoded protein	Major function	Reference(s)
22q11	COMT	Catechol-O-methyltransferase	Monoamine metabolism	[14-16]
6p22	DTNBP1	Dystrobrevin binding protein 1	Organelle biogenesis	[15,16]
8p12	NRG1	Neuregulin 1	Signal transduction in cell-cell interactions	[14-18]
1q23.3	RGS4	Regulator of G protein signaling 4	Regulation of G proteins	[15,16]
7q21.1	GRM3	Metabotropic glutamate receptor 3	Excitatory glutamatargic neurotransmission	[15]
13q34	G72	D-Amino acid oxidase activator	Degradation of D-serine, involved in NMDA-type neurotransmission	[14,15]
8p21.3	PPP3CC	Catalytic subunit of protein phosphatase 3	Regulation of dopaminergic signal transduction	[14,16]
15q14	CHRNA7	Nicotinic acetylcholine receptor	Cholinergic synaptic transmission	[16]
22q11.21	PRODH	Proline dehydrogenase (oxidase) 1	Proline degradation	[14,16]
6p21.3-22.1	MHC region	MHC proteins	Crucial parts in immune system	[20,21,23]
2q32.1	ZNF804A	Zinc finger protein 804A	Not fully characterized	[19,22]
18q21.2	TCF4	Transcription factor 4	Involved in nervous system development	[21,23]
11q24.2	NRGN	Neurogranin	Postsynaptic protein kinase pathway	[23]
1p21.3	MIR137	MicroRNA 137	Not fully characterized	[21]
Copy number or structural variants
Chromosomal location	Representative genes in the region			Reference(s)
1q21.1	CHD1L, GPR89B, BCL9, PRKAB2 and others			[24-26]
15q11.2	CYFIP1, TUBGCP5, NIPA1, NIPA2			[24]
15q13.3	CHRNA7, KLF13, OTUD7A, TRPM1 and others			[24,26]
16p11.2	KCTD13, MAPK3, TAOK2, CDIPT and others			[29]
22q11.2	COMT, PRODH, SLC25A1, CLDN5, TBX1 and others			[26,28]

MHC, major histocompatibility complex; NMDA, N-methyl-D-aspartate.

In addition to these common variants, rare structural variants with large effect size have been identified (Table 2). For example, rare copy number variations (CNVs) in 1q21.1, 15q11.2 and 15q13.3 have been associated with the risk for schizophrenia [24-26]. Increased occurrence of rare CNVs in patients has also been reported [27]. Microdeletion in 22q11.2 [28] and microduplication in 16p11.2 [29] have also been reported as risk factors for schizophrenia. Intriguingly, genes disrupted by rare variants are enriched for neuronal functions such as N-methyl-D-aspartate (NMDA) receptor-related postsynaptic signaling [30]. There is no definitive consensus on the genetic risk factors for schizophrenia, but abnormalities in NMDA-receptor-related glutamatergic synaptic transmission have been suggested as a possible mechanism underlying the etiology of schizophrenia at the molecular level. Although many candidate genes and CNV alterations have been identified as genetic risk factors, the shortcomings of these findings include low reproducibility, low risk ratio, low penetrance, and lack of disease specificity. The heterogeneity of schizophrenia may be one reason for these limitations, but other mechanisms such as epigenetics might also underlie this phenomenon [31,32].

Epigenetic modifications are roughly divided into two categories: cytosine and chromatin/histone modifications. In mammals, DNA methylation of cytosine residues occurs on the cytosine of CpG dinucleotides, catalyzed by DNA methyltransferases (DNMTs). Methyl-CpG-binding protein 2 (MeCP2) and methyl-CpG-binding domain (MBD) proteins can bind to methylated CpG regions [33]. Chromatin conformation changes in response to histone modifications such as acetylation and methylation. Histone acetylation is a major histone modification, and the process is regulated by histone acetyltransferases and histone deacetylases (HDACs). DNA methylation and histone modifications affect gene expression by changing the profiles of proteins that bind to specific DNA regions; for example, the binding profiles of transcription factors to promoter regions. In addition, DNA methylation on the gene body is thought to affect mRNA processing such as alternative splicing [34,35].

Considering that the epigenome can change in response to environmental factors [36], epigenetic alterations could partly explain G × E, and therefore could help to explain the etiology of schizophrenia (Figure 1). In this review, we discuss recent studies of epigenetic alterations in schizophrenia using postmortem brains and peripheral tissue samples, focusing mainly on DNA methylation, a critical component of epigenetic regulation. Other forms of epigenetic modification have been reviewed recently [37-40]. We also consider important challenges and future directions in this field.

Figure 1

A schematic view of the possible relationship between schizophrenia symptoms, epigenetic alterations and other factors. Genetic factors and environmental factors can affect epigenetic status either independently or through their interactions, and some epigenetic alterations could lead to schizophrenia symptoms. Epigenetic alterations might also partly explain gene-environment interactions (G × E). Examples of genetic factors, environmental factors and epigenetic alterations are shown.

Epigenetic studies of schizophrenia

The main findings of epigenetic studies on schizophrenia using human samples are summarized in Table 3. In analyzing DNA methylation, one of the traditional approaches is bisulfite sequencing. Sodium bisulfite treatment converts unmethylated cytosine to uracil, but does not convert methylated cytosine. Thus, sequencing DNA after bisulfite modification can allow discrimination between methylated and unmethylated cytosine. Other major approaches include utilization of methylation-sensitive restriction enzymes, antibodies against methylated cytosine, and methyl-cytosine binding proteins such as MBD1 and MBD2.

Table 3

Epigenetic alterations in schizophrenia reported in studies using human samples

Sample tissue	Approach	Number of samples	Main findings	Implicated genes	Reference
Postmortem occipital cortex and prefrontal cortex	Bisulfite sequencing*	15 Sz, 15 CT	Hypermethylation of RELN in Sz	RELN encodes a protein related to cell positioning and neuronal migration during brain development	[46]
Postmortem prefrontal cortex	Bisulfite sequencing*	11 Sz, 12 CT	Hypermethylation of SOX10 in Sz; downregulation of SOX10 transcription in Sz	SOX10 encodes an oligodendrocyte-specific transcription factor	[53]
Postmortem frontal cortex	Bisulfite sequencing*; methylation-specific PCR	35 Sz, 35 CT	Hypomethylation of MB-COMT promoter in Sz	MB-COMT encodes the membrane-bound form of catechol-O-methyltransferase, which is involved in dopamine metabolism	[51]
Postmortem cerebellum	Pyrosequencing after bisulfite-PCR*	15 Sz, 15 CT	No significant difference in COMT methylation		[52]
Postmortem frontal cortex	CpG island microarray; pyrosequencing after bisulfite-PCR*	35 Sz, 35 CT	Methylation changes in genes related to glutamatergic/GABAergic neurotransmission, brain development, mitochondrial function and stress response in Sz	Examples include: NR3B and GRIA2, both of which are involved in glutamatergic neurotransmission, and implicated in the etiology of schizophrenia	[49]
Postmortem prefrontal cortex	Pyrosequencing after bisulfite-PCR*	15 Sz, 15 CT	No significant difference in RELN methylation		[50]
Postmortem parahippocampus gyrus	Bisulfite sequencing*	6 right and 7 left hemisphere of Sz; 5 right and 6 left hemisphere of CT	Hypermethylation of FOXP2 in left parahippocampus gyrus in Sz	FOXP2 encodes a transcriptional factor required for the development of language	[56]
Postmortem frontal lobe	Bisulfite sequencing*Methylation specific PCR	35 Sz, 35 CT	Hypermethylation at and around -1438A/G SNP, hypomethylation at and around T102C SNP of HTR2A promoter in Sz	HTR2A encodes serotonin receptor 2A, one of the main target molecules of antipsychotic drugs	[57]
Peripheral blood	Bisulfite sequencing*	1 MZ pair concordant for Sz; 1 MZ pair discordant for Sz	SZ twin in a pair discordant for SZ had more similar DRD2 methylation profiles to the affected concordant twin pair than to its unaffected cotwin	DRD2 encodes the D2 subtype of the dopamine receptor	[74]
Peripheral blood	High-performance liquid chromatography	210 Sz (124 male, 86 female); 237 CT (108 male, 129 female)	Global hypomethylation in male Sz		[66]
Peripheral blood	Radiolabeled [3H]dCTP-extension assay	28 Sz, 26 CT	No difference in global DNA methylation		[68]
Peripheral blood	Bisulfite sequencing*	30 Sz, 30 CT	Hypermethylation in the neighborhood of disease-associated SNP rs1816071; hypomethylation at three CpG sites; hypermethylation at four CpG sites, in GABRB2 in Sz	GABRB2 encodes beta 2 subunit of GABA(A) receptor	[71]
Peripheral blood	High-resolution melt assay*	40 Sz, 67 CT	Hypermethylation 5HTR1A promoter in Sz	5HTR1A encodes serotonin receptor 1A	[70]
Saliva	Bisulfite sequencing*; methylation specific PCR	63 Sz, 76 CT	Hypomethylation of T102C SNP in 5HTR2A in Sz	5HTR2A encodes serotonin receptor 2A	[72]
Saliva	Bisulfite sequencing*; methylation specific PCR	63 Sz, 76 CT	Hypomethylation of MB-COMT promoter in Sz	MB-COMT encodes membrane-bound form of catechol-O-methyltransferase, which is involved in dopamine metabolism	[73]
Peripheral blood	Bead array (Illumina Infinium Human-Methylation27)*	11 MZ pairs discordant for Sz	Hypomethylation of ST6GALNAC1 promoter in psychosis; PUS3 is a top-ranked differentially methylated gene in Sz; differential methylation in psychological-disorder-related genes such as ADAMTS3 and SLC6A3 in Sz	ST6GALNAC1 encodes an enzyme that transfers sialic acid to O-linked N-acetylgalactosamine residues, and PUS3 encodes pseudouridylate synthase 3, both of which require further investigation	[62]
Peripheral blood	Global assay (restriction enzyme and pyrosequencing); pyrosequencing after bisulfite-PCR*	177 Sz, 171 CT; (S-COMT: 47 Sz, 47 CT; SLC6A4: 67 Sz, 81 CT)	Global hypomethylation in SzEarly-onset-associated lower global methylation; haloperidol-associated higher global methylation; hypermethylation of S-COMT in Sz; no significant difference in SLC6A4 methylation	S-COMT encodes the soluble form of catechol-O-methyltransferase, which is involved in dopamine metabolism; SLC6A4 encodes a serotonin transporter	[67]

*The limitation of bisulfite treatment (this approach does not allow precise discrimination between various forms of cytosine modifications).

CT, normal controls; DZ, dizygotic twin; GABA, γ-aminobutyric acid, MZ, monozygotic twin; SNP, single nucleotide polymorphism, Sz, schizophrenia (patients).

Studies using postmortem brain samples

Postmortem brains are the most important samples for epigenetic research on schizophrenia, but there are many confounding factors in such studies, including terminal condition (that is, the condition of the samples at death, or their 'agonal state'), cause of death and postmortem condition. In addition, demographic factors such as age, gender and medical history may also affect experimental data. The pH of brain samples, which reflects the agonal state of the subject, is a major factor that affects transcriptome analysis [41-43]. However, Ernst et al. [44] reported that the pH of brain samples has no significant effect on DNA methylation status. Although DNA methylation is thought to be chemically stable, systematic evaluation of the effects of other factors on epigenetic alterations present in postmortem brain samples is necessary. Nevertheless, studies of postmortem brain samples of schizophrenia patients have revealed epigenetic alterations in several genes in specific brain regions, although some findings have not been confirmed by other groups (Table 3).

The GABAergic-neuron-expressed genes GAD67 and RELN have been shown to be underexpressed at both the mRNA and the protein level in postmortem brains from schizophrenia patients [45]. GAD67 encodes an enzyme that catalyzes production of γ-aminobutyric acid (GABA), which is a major neurotransmitter of inhibitory neurons. Altered GABAergic neurotransmission has been described in schizophrenia, and GABA receptor agonists, mainly categorized as benzodiazepines, are occasionally prescribed to schizophrenia patients to alleviate anxiety or sleep disturbance, as an adjuvant therapy to antipsychotic drugs. RELN encodes an extracellular matrix protein involved in cell positioning and neuronal migration during brain development. GAD67 and RELN underexpression in the prefrontal cortex appears to be correlated with their promoter hypermethylation [46,47]. In addition, the gene encoding DNA methyltransferase 1 (DNMT1) was found to be overexpressed in GABAergic interneurons in postmortem brain samples from schizophrenia patients [48]. However, other groups found no correlation between schizophrenia and RELN promoter methylation levels in (pre)frontal cortices [49,50].

COMT has also been identified as a candidate gene in schizophrenia, and encodes catechol-O-methyltransferase, which is involved in dopamine metabolism. Dopamine is thought to be involved in the molecular pathways underpinning positive symptoms, including hallucinations and delusions. This gene is located on 22q11, and microdeletion of this region has shown strong correlation with schizophrenia [28]. The methylation status of COMT in schizophrenia is controversial. COMT promoter hypomethylation and corresponding COMT mRNA overexpression in left frontal lobes have been reported [51]. However, another group did not confirm the association between COMT methylation levels and schizophrenia in cerebellum, attributing differences in COMT expression to variation in SNPs [52].

We previously reported that SOX10, which encodes an oligodendrocyte-specific transcription factor, tended to be highly methylated in the prefrontal cortices of schizophrenia patients, and this correlated with the reduced expression of SOX10 and other oligodendrocyte-related genes [53]. SOX10 underexpression can result in oligodendrocyte dysfunction with the downregulation of key oligodendrocyte and myelination genes, which has been observed in studies utilizing postmortem brain samples from schizophrenia patients [43,54,55].

DNA methylation alterations of other genes, such as FOXP2 (which encodes a transcription factor required for language development) [56] and HTR2A (which encodes serotonin receptor 2A, one of the main target molecules of antipsychotic drugs) [57], have also been reported in postmortem brain samples from schizophrenia patients. Mill et al. [49] first performed comprehensive DNA methylation analysis in postmortem frontal cortices from schizophrenia patients, investigating more than 10,000 CpG islands using CpG-island microarrays. They reported DNA methylation changes in genes related to glutamatergic/GABAergic neurotransmission, brain development, mitochondrial function and stress response.

Altered histone modifications have also been reported in relation to schizophrenia. Akbarian et al. [58] reported that high levels of H3-(methyl)arginine 17 were associated with downregulated expression of metabolic genes in the prefrontal cortex of schizophrenia patients. Huang et al. [59] reported reduced GAD67 expression associated with decreased levels of trimethylated lysine 4 of histone H3 (H3K4me3), a hallmark of gene activation, in prefrontal cortex samples derived from schizophrenia patients, and suggested that chromatin remodeling mechanisms at GABAergic gene promoters might be involved in the pathophysiology of schizophrenia. Tang et al. [60] reported the underexpression of eight genes, including GAD67 and HTR2C (which encodes serotonin receptor 2C), associated with hypoacetylation at lysine 9/14 of histone H3 (H3K9K14) in postmortem prefrontal cortex samples from young schizophrenia patients, indicating the transcriptional effect of histone alterations on the investigated genes.

A number of interesting epigenetic alterations of several important genes, possibly related to schizophrenia, have been identified in postmortem brain samples so far. However, most of these studies involved candidate gene approaches, and more comprehensive analyses of epigenetic changes in brains of patients with schizophrenia as well as the effects of confounding factors such as medication remain to be explored.

Studies using blood and other peripheral samples

In contrast to brain samples, peripheral tissue samples can be obtained from living patients with minimally invasive techniques, and this also enables larger and multiple sample collection. Peripheral tissue samples can also be used in longitudinal studies, in order to investigate the course of disease. Although epigenetic status greatly differs between tissue types, brain samples and peripheral tissue samples could have common epigenetic status in certain genomic regions, and thus could be used to identify epigenetic markers of disease. One study reported that COMT promoter methylation status was nearly identical in brain and blood samples in schizophrenia patients [61]. Dempster et al. [62] reported that methylation alterations in ST6GALNAC1 were observed both in peripheral blood samples and in postmortem brain samples obtained from schizophrenia patients. ST6GALNAC1 encodes α-N-acetylgalactosaminide α-2,6-sialyltransferase 1, an enzyme that transfers sialic acid to O-linked N-acetylgalactosamine residues, but its function has not been fully elucidated. Similarly, studies of bipolar disorder also demonstrated that methylation status at specific CpG sites in SLC6A4 (which encodes a serotonin transporter, one of the main targets of antidepressants) and human leukocyte antigen complex group 9 gene (HCG9, which is located within the MHC class I region) showed a common pattern in brain and blood samples [63,64]. More recently, Davies et al. [65] also reported that DNA methylation differences in blood samples and those in postmortem brain samples were well correlated in the investigated genomic regions, indicating the utility of peripheral tissue samples in studying mental disorders. Several studies have reported epigenetic alterations in peripheral samples (including blood or saliva) from schizophrenia patients; however in some cases the findings have not been supported by work from other groups (Table 3).

Leukocyte DNA from schizophrenia patients has been reported to be globally hypomethylated [66,67]. However, another group reported no global methylation differences between schizophrenia patients and controls [68], and others have reported overexpression of DNMT1 and DNMT3A mRNA in peripheral blood samples from schizophrenia patients [69]. Given that DNMT proteins promote DNA methylation, these results seem to be contradictory. Interestingly, Melas et al. [67] reported global hypomethylation, especially in patients with early onset schizophrenia (a severe form of disease) compared with controls, indicating a correlation between disease severity and global methylation status. They also reported hypermethylation of S-COMT (which encodes the soluble form of COMT) and no differential methylation of SLC6A4 in peripheral blood cells derived from schizophrenia patients. Other groups have reported hypermethylation of the promoter region of 5HTR1A (which encodes serotonin receptor 1A) in leukocytes [70], and hypermethylation of the region around the schizophrenia-associated SNP (rs1816071) of GABRB2 (which encodes the β2 subunit of GABA(A) receptor) in leukocytes [71]. Further studies will be needed to understand the patterns of DNA methylation in leukocytes from schizophrenia patients.

Altered DNA methylation has also been observed in saliva samples, including hypomethylation of 5HTR2A [72] (in agreement with a study of postmortem brain samples [57]), and hypomethylation of COMT [73] (in agreement with a study of postmortem brain samples [51]) in saliva from schizophrenia patients.

Studies of peripheral samples from monozygotic twins have also revealed interesting findings. Petronis et al. [74] analyzed DNA methylation in the regulatory region of DRD2 (which encodes the D2 subtype of the dopamine receptor) in lymphocytes derived from two sets of monozygotic twins, one concordant and one discordant for schizophrenia. They reported that the affected twin of the discordant pair was epigenetically more similar to the affected concordant twins than to the unaffected co-twin, by analyses of methylation differences at each CpG site [74]. More recently, Dempster et al. [62] analyzed genome-wide DNA methylation status in peripheral blood DNA samples obtained from monozygotic twins discordant for major psychosis, including schizophrenia and bipolar disorder. They found significant DNA hypomethylation in the promoter of ST6GALNAC1 in the affected twins with major psychosis, and this hypomethylation was also observed in independent postmortem brain samples [62].

Recently, Kano et al. [75] reported that olfactory cells of schizophrenia patients exhibited alterations in the expression of several genes, including MGST1 (which encodes microsomal glutathione S-transferase 1), and alterations in H3K4me3 and trimethylated lysine 27 of histone H3 (H3K27me3), which are hallmarks of gene activation and repression, respectively, in the same genes. The utilization of olfactory cells is advantageous in psychiatric research because olfactory cells in olfactory bulbs are thought to share cellular traits with neurons in the brain.

Although the precise mechanisms that affect epigenetic status in peripheral tissue samples and their pathophysiological significance in schizophrenia remain unclear, epigenetic alterations in peripheral tissues are useful for development of potential biomarkers for monitoring disease status. Performing more comprehensive epigenome-scale studies using larger sample sizes would be expected to provide further insights into these mechanisms and to facilitate biomarker identification.

Investigating the epigenetic actions of drugs

VPA, an HDAC inhibitor and the best-known psychotropic drug targeting epigenetic mechanisms, is often prescribed to schizophrenia patients, especially to affective schizophrenia patients, as an adjuvant to antipsychotic drugs. It has been actively investigated using animal models. Several groups have reported the therapeutic effects of VPA in animal models with schizophrenia-like behavior or with neglected upbringing, which were associated with modifications in DNA methylation or histone acetylation [76-78]. Here, we focus on studies investigating the epigenetic effects of VPA using peripheral samples from schizophrenia patients.

VPA has been reported to induce a significant increase in the levels of acetylated histone H3 (acH3) and acetylated histone H4 (acH4) in lymphocytes. The increase in acH3 levels positively correlated with VPA concentration, and the VPA-induced increases in acH3 and acH4 were found to be smaller in schizophrenia patients [79]. More recently, Gavin et al. [80,81] reported reduced baseline acH3 levels and elevated dimethylated lysine 9 of histone H3 levels in lymphocytes obtained from schizophrenia patients, which were less responsive to HDAC inhibitors than controls. They also reported that VPA induced a significant increase in GAD67 mRNA levels, acetylation of H3K9K14, and attachment of acetylated H3K9K14 to the promoter of GAD67 in lymphocytes from human subjects, including schizophrenia patients [82]. Taken together, these findings suggest that histone acetylation may be reduced in several genomic regions in schizophrenia patients, and that VPA treatment can increase the levels of histone acetylation in such genomic regions.

Further studies on VPA and other psychotropic drugs using human peripheral samples are needed to better understand the relationship between the effects of these drugs on epigenetic status and their mechanisms of action.

Variations in the epigenome

Most of the previous epigenetic studies on schizophrenia have focused on DNA methylation, but recent research has identified other forms of cytosine modification. The physiological and pathophysiological roles and distribution of such newly discovered cytosine modifications are yet to be understood, but it may be necessary to investigate their importance in schizophrenia. Recent research has also revealed the heterogeneity of DNA methylation patterns between tissue types or cell types, and other types of epigenomic variations, for example during development and aging, which represent additional challenges for epigenetic studies on schizophrenia.

Dynamic cytosine modifications in the brain

DNA methylation had been considered to be chemically stable, but recent research has revealed that cytosine modification is a dynamic process. Recent studies have shown that 5-hydroxymethylcytosine (5hmC) is abundant in mammalian brain cells [83] and is converted from 5-methylcytosine (5mC) by ten-eleven translocation (TET) proteins [84]. 5hmC might be an intermediate product of active demethylation, although its existence had been considered controversial in mammals [85,86]. In addition, several other modifications have been described, including 5-hydroxymethyluracil (5hmU), 5-formylcytosine (5fC) and 5-carboxylcytosine (5caC) [85,87-89]. 5fC and 5caC might also be intermediate products of demethylation [88,89]. The pathways that might be involved in these cytosine modifications are shown in Figure 2.

Figure 2

Cytosine modification pathways. Recent research has revealed that cytosine modifications can be changed dynamically. DNA methyltransferases (DNMTs) catalyze the conversion of cytosine to 5-methylcytosine (5mC). Ten-eleven translocation (TET) proteins catalyze the synthesis of 5-hydroxymethylcytosine (5hmC), 5-formylcytosine (5fC) or 5-carboxylcytosine (5caC). Both 5fC and 5caC are thought to be replaced by cytosine by DNA repair processes such as base excision repair. Thus, 5hmC, 5fC and 5caC are thought to be intermediate products of demethylation. The synthesis of 5-hydroxymethyluracil (5hmU) or thymine is thought to be catalyzed by activation-induced deaminases/apolipoprotein B mRNA editing enzyme, catalytic polypeptide-like (AID/APOBEC) cytidine deaminases from 5hmC or 5mC, respectively. However, recent studies have questioned these reactions [131].

Previous studies of schizophrenia have not provided definitive results on the presence of 5mC or 5hmC, partly because conventional bisulfite sequencing cannot discriminate between these forms. Therefore, studies will be required to investigate their distribution and potential pathophysiological roles in schizophrenia using techniques that can distinguish between these cytosine modifications [90-94].

Other variations in the epigenome

Epigenetic diversity has been reported to accumulate in brain cells, germ cells, peripheral blood cells and other tissues over time [95-97], although the effects of aging on specific epigenetic alterations are unclear. Hernandez et al. [98] reported that DNA methylation tends to increase with chronological age in the human brain. On the other hand, Heyn et al. [99] reported that aged individuals had less genome-wide DNA methylation in peripheral blood leukocytes than newborns. Numata et al. [100] reported that dynamic changes in DNA methylation occur in human prefrontal cortex especially during the prenatal period, typically with prenatal decreases in methylation and postnatal increases in methylation.

Monozygotic twins have the same DNA sequences, except in rare cases such as individual de novo mutations [101]. The epigenomes of monozygotic twins are also more similar than those of dizygotic twins [102]. Although monozygotic twins are epigenetically indistinguishable early in life, they each experience different patterns of DNA methylation and chromatin modifications, thus affecting gene expression in adolescence and adulthood [103]. Recent studies have revealed that monozygotic twins can have significantly different epigenetic status in genes that have been implicated in psychiatric disorders, specifically DRD4 (which encodes D4 subtype of the dopamine receptor), SLC6A4, and MAOA (which encodes monoamine oxidase A, an enzyme that deaminates dopamine and serotonin) as early as ages 5 to 10 years [104]. Epigenetic differences could therefore be one reason for monozygotic twin discordance in schizophrenia [74,101].

Epigenetic variations could also be partly explained by differences in DNA sequence, such as SNPs and CNVs [100,105,106]. For example, Numata et al. [100] reported more DNA methylation in one allele of HLA-DQB1 (which encodes a part of human leukocyte antigen class II molecules, with important functions in the immune system) that contained a SNP (rs1063355) than in the other allele in human postmortem prefrontal cortex samples.

In addition to genetic and environmental factors, Petronis argued that epigenetic status, especially DNA methylation, could also be altered by chance [107]. This stochasticity of changes in DNA methylation might also contribute to epigenetic variation and might be an important factor in complex diseases such as schizophrenia.

Among human tissues, the brain has a particularly complex anatomical structure with functionally different subregions. Epigenetic differences have been reported across different brain subregions [108]. However, it remains unclear which specific subregions are involved in the pathophysiology of schizophrenia, and there is study bias for regions of the frontal cortex and hippocampus, because these regions have been well characterized in previous research. As well as epigenetic differences at the subregion level, each subregion contains various differentiated cells such as neurons and glia. Extensive transcriptome analyses have revealed that each of these cell types has unique expression profiles [109,110], implicating the presence of distinct gene regulation, a large part of which could be attributable to epigenetic mechanisms. Using postmortem brain tissues, neuronal and non-neuronal nuclei can be separated by cell sorting based on NeuN, a neuronal marker [111,112]. Recently, oligodendrocyte nuclei were also successfully separated from other non-neuronal nuclei, using an anti-Olig2-antibody [113]. DNA extracted from the separated fractions can be used for epigenetic analysis of histone modifications and DNA methylation [114,115]. Interestingly, studies of neurons have shown not only unique DNA methylation signatures, but also higher interindividual variation compared to non-neurons [114]. These results suggest that neurons might be more susceptible to epigenetic variations due to environmental factors.

Even analyses at the cell-population level may not fully resolve the epigenetic complexity in the brain, because neurons and glia can be further divided into numerous cell types, and each cell type might have a different epigenetic profile. Although laser-assisted microdissection can be useful when analyzing particular cells, other techniques will be required to analyze the epigenome at the single cell level. Recently, Navin et al. [116] quantified genomic CNVs in a single nucleus by whole-genome amplification followed by next-generation sequencing of flow-sorted nuclei, and Bendall et al. [117] reported the results of single-cell mass cytometry. Furthermore, Kantlehner et al. [118] presented a candidate method for single-cell DNA methylation analysis utilizing methylation sensitive restriction enzymes and PCR amplification.

Variation in the epigenome between tissues and cells represents an enormous challenge in epigenetic studies of schizophrenia. Considering that the human brain consists of a larger number of glial cells compared with neuronal cells, it will be necessary to separate distinct cell populations as much as possible when analyzing epigenetic status in brain samples.

Transgenerational epigenetic inheritance

DNA methylation can be faithfully transmitted from a mother cell to daughter cells during mitosis [119], but a long-standing question is whether transgenerational transmission occurs during the generation of offspring [120].

Mammalian germ-cell-derived DNA is thought to be fully demethylated in early embryos and in primordial germ cells [87]. However, several recent reports have argued that some epigenetic marks survive in early embryos and remain until adulthood. Two groups studying mice reported that epigenetic alterations induced by social stress could be transmitted to the next generation [121,122]. Another report demonstrated in mice that parental diet could affect metabolism-related epigenetic status of the offspring [123]. Other groups have reported the existence of transgenerational epigenetic inheritance in nematode and fly models [124,125]. However, further studies are needed to investigate epigenetic inheritance in humans. If epigenetic alterations were shown to be transmitted transgenerationally, they could partly explain the heritability of schizophrenia.

Conclusions and future perspectives

Given that epigenetic status, including DNA methylation and chromatin/histone modifications, affects gene expression, epigenetic alterations could result in functional deficits in brain cells in schizophrenia. Many studies have identified epigenetic alterations in samples obtained from schizophrenia patients, and here we have focused mainly on studies of DNA methylation. However, due to the small sample sizes and technical difficulties, most of the early results were not well replicated [126]. Therefore, the mechanisms underlying the epigenetic alterations in schizophrenia and how these might lead to the symptoms of schizophrenia remain to be understood.

To investigate the importance of previous findings, it will be necessary to perform independent replication studies using larger sample sizes. Most importantly, more comprehensive genome-scale analysis using array-based platforms or next-generation sequencing will be needed [127]. Among these technologies, bisulfite-modified whole-genome sequencing approaches have been utilized for methylome analysis at base pair resolution [35,128]. Recent studies have revealed complex and dynamic patterns of DNA methylation in the human brain. It will be important to investigate various cytosine and histone modifications in specific brain cell types to understand their potential roles in the pathophysiology of schizophrenia. Recent technologies that can discriminate between various cytosine modifications on a genome-wide scale [90-94] as well as cell-sorting-based analysis [113-115] will empower further characterization of cytosine modifications in brain samples of schizophrenia patients.

Previous research using human samples has shown only associations between specific epigenetic changes and schizophrenia, so it remains unclear whether the identified epigenetic alterations are the causes or the results of schizophrenia pathology. Several groups have reported changes in gene expression associated with epigenetic alterations, but further studies are needed to examine these relationships. After careful and extensive replication studies, the relationship between candidate epigenetic alterations and schizophrenia pathology must be validated using cellular and animal models. Methods for site-specific epigenetic engineering could be used to investigate causality [129,130], although these technologies are still immature. If causality between epigenetic alterations and schizophrenia was demonstrated, therapeutic drugs could be developed by searching for molecules affecting the target epigenetic alterations [40].

In conclusion, despite many challenges regarding the approaches for sampling and detection of epigenetic alterations, epigenetic studies are expected to contribute to our understanding of the pathophysiology of schizophrenia and should provide novel opportunities for the development of therapeutic drugs.

Abbreviations

acH3: acetylated histone H3; acH4: acetylated histone H4; CNV: copy number variation; DNMT: DNA methyltransferase; G × E: gene-environment interaction; GABA: γ-aminobutyric acid; H3K27me3: trimethylated lysine 27 of histone H3; H3K4me3: trimethylated lysine 4 of histone H3; H3K9K14: lysine 9/14 of histone H3; HDAC: histone deacetylase; MBD: methyl-CpG-binding domain; MeCP2: methyl-CpG-binding protein 2; MHC: major histocompatibility complex; NMDA: N-methyl-D-aspartate; SNP: single nucleotide polymorphism; TET: ten-eleven translocation; VPA: valproate; 5caC: 5-carboxylcytosine; 5fC: 5-formylcytosine; 5hmC: 5-hydroxymethylcytosine; 5hmU: 5-hydroxymethyluracil; 5mC: 5-methylcytosine.

Competing interests

The authors declare that they have no competing interests.