Systematic investigation of allelic regulatory activity of schizophrenia-associated common variants

Genome-wide association studies (GWASs) have successfully identified 145 genomic regions that contribute to schizophrenia risk, but linkage disequilibrium makes it challenging to discern causal variants. We performed a massively parallel reporter assay (MPRA) on 5,173 fine-mapped schizophrenia GWAS variants in primary human neural progenitors and identified 439 variants with allelic regulatory effects (MPRA-positive variants). Transcription factor binding had modest predictive power, while fine-map posterior probability, enhancer overlap, and evolutionary conservation failed to predict MPRA-positive variants. Furthermore, 64% of MPRA-positive variants did not exhibit expressive quantitative trait loci signature, suggesting that MPRA could identify yet unexplored variants with regulatory potentials. To predict the combinatorial effect of MPRA-positive variants on gene regulation, we propose an accessibility-by-contact model that combines MPRA-measured allelic activity with neuronal chromatin architecture.

Location matters: distinct DNA methylation patterns in GABAergic interneuronal populations from separate microcircuits within the human hippocampus

Recent studies describe distinct DNA methylomes among phenotypic subclasses of neurons in the human brain, but variation in DNA methylation between common neuronal phenotypes distinguished by their function within distinct neural circuits remains an unexplored concept. Studies able to resolve epigenetic profiles at the level of microcircuits are needed to illuminate chromatin dynamics in the regulation of specific neuronal populations and circuits mediating normal and abnormal behaviors. The Illumina HumanMethylation450 BeadChip was used to assess genome-wide DNA methylation in stratum oriens GABAergic interneurons sampled by laser-microdissection from two discrete microcircuits along the trisynaptic pathway in postmortem human hippocampus from eight control, eight schizophrenia, and eight bipolar disorder subjects. Data were analysed using the minfi Bioconductor package in R software version 3.3.2. We identified 11 highly significant differentially methylated regions associated with a group of genes with high construct-validity, including multiple zinc finger of the cerebellum gene family members and WNT signaling factors. Genomic locations of differentially methylated regions were highly similar between diagnostic categories, with a greater number of differentially methylated individual cytosine residues between circuit locations in bipolar disorder cases than in schizophrenia or control (42, 7, and 7 differentially methylated positions, respectively). These findings identify distinct DNA methylomes among phenotypically similar populations of GABAergic interneurons functioning within separate hippocampal subfields. These data compliment recent studies describing diverse epigenotypes among separate neuronal subclasses, extending this concept to distinct epigenotypes within similar neuronal phenotypes from separate microcircuits within the human brain.

Epigenetic Regulation of Glutamic Acid Decarboxylase 67 in a Hippocampal Circuit

GABAergic dysfunction in hippocampus, a key feature of schizophrenia (SZ), may contribute to cognitive impairment in this disorder. In stratum oriens (SO) of sector CA3/2 of the human hippocampus, a network of genes involved in the regulation of glutamic acid decarboxylase GAD67 has been identified. Several of the genes in this network including epigenetic factors histone deacetylase 1 (HDAC1) and death-associated protein 6 (DAXX), the GABAergic enzyme GAD65 as well as the kainate receptor (KAR) subunits GluR6 and 7 show significant changes in expression in this area in SZ. We have tested whether HDAC1 and DAXX regulate GAD67, GAD65, or GluR in the intact rodent hippocampus. Stereotaxic injections of lentiviral vectors bearing shRNAi sequences for HDAC1 and DAXX were delivered into the SO of CA3/2, followed by laser microdissection of individual transduced GABA neurons. Quantitative PCR (QPCR) analyses demonstrated that inhibition of HDAC1 and DAXX increased expression of GAD67, GAD65, and GluR6 mRNA. Inhibition of DAXX, but not HDAC1 resulted in a significant increase in GluR7 mRNA. Our data support the hypothesis that HDAC1 and DAXX play a central role in coordinating the expression of genes in the GAD67 regulatory pathway in the SO of CA3/2.

Toward dissecting the etiology of schizophrenia: HDAC1 and DAXX regulate GAD67 expression in an in vitro hippocampal GABA neuron model

Schizophrenia (SZ) is associated with GABA neuron dysfunction in the hippocampus, particularly the stratum oriens of sector CA3/2. A gene expression profile analysis of human postmortem hippocampal tissue followed by a network association analysis had shown a number of genes differentially regulated in SZ, including the epigenetic factors HDAC1 and DAXX. To characterize the contribution of these factors to the developmental perturbation hypothesized to underlie SZ, lentiviral vectors carrying short hairpin RNA interference (shRNAi) for HDAC1 and DAXX were used. In the hippocampal GABA neuron culture model, HiB5, transduction with HDAC1 shRNAi showed a 40% inhibition of HDAC1 mRNA and a 60% inhibition of HDAC1 protein. GAD67, a enzyme associated with GABA synthesis, was increased twofold (mRNA); the protein showed a 35% increase. The expression of DAXX, a co-repressor of HDAC1, was not influenced by HDAC1 inhibition. Transduction of HiB5 cells with DAXX shRNAi resulted in a 30% inhibition of DAXX mRNA that translated into a 90% inhibition of DAXX protein. GAD1 mRNA was upregulated fourfold, while its protein increased by ~30%. HDAC1 expression was not altered by inhibition of DAXX. However, a physical interaction between HDAC1 and DAXX was demonstrated by co-immunoprecipitation. Inhibition of HDAC1 or DAXX increased expression of egr-1, transcription factor that had previously been shown to regulate the GAD67 promoter. Our in vitro results point to a key role of both HDAC1 and DAXX in the regulation of GAD67 in GABAergic HiB5 cells, strongly suggesting that these epigenetic/transcription factors contribute to mechanisms underlying GABA cell dysfunction in SZ.

Circuit- and Diagnosis-Specific DNA Methylation Changes at γ-Aminobutyric Acid-Related Genes in Postmortem Human Hippocampus in Schizophrenia and Bipolar Disorder

IMPORTANCE

Dysfunction related to γ-aminobutyric acid (GABA)-ergic neurotransmission in the pathophysiology of major psychosis has been well established by the work of multiple groups across several decades, including the widely replicated downregulation of GAD1. Prior gene expression and network analyses within the human hippocampus implicate a broader network of genes, termed the GAD1 regulatory network, in regulation of GAD1 expression. Several genes within this GAD1 regulatory network show diagnosis- and sector-specific expression changes within the circuitry of the hippocampus, influencing abnormal GAD1 expression in schizophrenia and bipolar disorder.

OBJECTIVE

To investigate the hypothesis that aberrant DNA methylation contributes to circuit- and diagnosis-specific abnormal expression of GAD1 regulatory network genes in psychotic illness.

DESIGN, SETTING, AND PARTICIPANTS

This epigenetic association study targeting GAD1 regulatory network genes was conducted between July 1, 2012, and June 30, 2014. Postmortem human hippocampus tissue samples were obtained from 8 patients with schizophrenia, 8 patients with bipolar disorder, and 8 healthy control participants matched for age, sex, postmortem interval, and other potential confounds from the Harvard Brain Tissue Resource Center, McLean Hospital, Belmont, Massachusetts. We extracted DNA from laser-microdissected stratum oriens tissue of cornu ammonis 2/3 (CA2/3) and CA1 postmortem human hippocampus, bisulfite modified it, and assessed it with the Infinium HumanMethylation450 BeadChip (Illumina, Inc). The subset of CpG loci associated with GAD1 regulatory network genes was analyzed in R version 3.1.0 software (R Foundation) using the minfi package. Findings were validated using bisulfite pyrosequencing.

MAIN OUTCOMES AND MEASURES

Methylation levels at 1308 GAD1 regulatory network-associated CpG loci were assessed both as individual sites to identify differentially methylated positions and by sharing information among colocalized probes to identify differentially methylated regions.

RESULTS

A total of 146 differentially methylated positions with a false detection rate lower than 0.05 were identified across all 6 groups (2 circuit locations in each of 3 diagnostic categories), and 54 differentially methylated regions with P < .01 were identified in single-group comparisons. Methylation changes were enriched in MSX1, CCND2, and DAXX at specific loci within the hippocampus of patients with schizophrenia and bipolar disorder.

CONCLUSIONS AND RELEVANCE

This work demonstrates diagnosis- and circuit-specific DNA methylation changes at a subset of GAD1 regulatory network genes in the human hippocampus in schizophrenia and bipolar disorder. These genes participate in chromatin regulation and cell cycle control, supporting the concept that the established GABAergic dysfunction in these disorders is related to disruption of GABAergic interneuron physiology at specific circuit locations within the human hippocampus.

Selective epigenetic alteration of layer I GABAergic neurons isolated from prefrontal cortex of schizophrenia patients using laser-assisted microdissection

Among the most consistent results of studies of post-mortem brain tissue from schizophrenia patients (SZP) is the finding that in this disease, several genes expressed by GABAergic neurons are downregulated. This downregulation may be caused by hypermethylation of the relevant promoters in affected neurons. Indeed, increased numbers of GABAergic interneurons expressing DNA methyltransferase 1 (DNMT1) mRNA have been demonstrated in the prefrontal cortex (PFC) of SZP using in situ hybridization. The present study expands upon these findings using nested competitive reverse transcription-polymerase chain reaction combined with laser-assisted microdissection to quantitate the extent of DNMT1 mRNA overexpression in distinct populations of GABAergic neurons obtained from either layer I or layer V of the PFC of SZP. In a cohort of eight SZP and eight non-psychiatric subject (NPS) post-mortem BA9 tissue samples, DNMT1 mRNA was found to be selectively expressed in GABAergic interneurons and virtually absent in pyramidal neurons. DNMT1 mRNA expression was approximately threefold higher in GABAergic interneurons microdissected from layer I of SZP relative to the same neurons microdissected from NPS. GABAergic interneurons obtained from layer V of the same samples displayed no difference in DNMT1 mRNA expression between groups. In the same samples, the GABAergic neuron-specific glutamic acid-decarboxylase(67) (GAD(67)) and reelin mRNAs were underexpressed twofold in GABAergic interneurons isolated from layer I of SZP relative to GABAergic interneurons microdissected from layer I of NPS, and unaltered in GABAergic interneurons of layer V. These findings implicate an epigenetically mediated layer I GABAergic dysfunction in the pathogenesis of schizophrenia, and suggest novel strategies for treatment of the disease.

An epigenetic mouse model for molecular and behavioral neuropathologies related to schizophrenia vulnerability

Reelin and glutamic acid decarboxylase (GAD)67 expressed by cortical gamma-aminobutyric acid-ergic interneurons are down-regulated in schizophrenia. Because epidemiological studies of schizophrenia fail to support candidate gene haploinsufficiency of Mendelian origin, we hypothesize that epigenetic mechanisms (i.e., cytosine hypermethylation of CpG islands present in the promoter of these genes) may be responsible for this down-regulation. Protracted l-methionine (6.6 mmolkg for 15 days, twice a day) treatment in mice elicited in brain an increase of S-adenosyl-homocysteine, the processing product of the methyl donor S-adenosyl-methionine, and a marked decrease of reelin and GAD67 mRNAs in both WT and heterozygous reeler mice. This effect of l-methionine was associated with an increase in the number of methylated cytosines in the CpG island of the reelin promoter region. This effect was not observed for GAD65 or neuronal-specific enolase and was not replicated by glycine doses 2-fold greater than those of l-methionine. Prepulse inhibition of startle declined at a faster rate as the prepulsestartle interval increased in mice receiving l-methionine. Valproic acid (2 mmolkg for 15 days, twice a day) reverted l-methionine-induced down-regulation of reelin and GAD67 in both WT and heterozygous reeler mice, suggesting an epigenetic action through the inhibition of histone deacetylases. The same dose of valproate increased acetylation of histone H3 in mouse brain nearly 4-fold. This epigenetic mouse model may be useful in evaluating drug efficacy on schizophrenia vulnerability. Hence the inhibition of histone deacetylases could represent a pharmacological intervention mitigating epigenetically induced vulnerability to schizophrenia in individuals at risk.