|Accueil | Aperçu | Revues | Soumettre | Nous Contacter | English|
Author contributions: R.B., J.L., P.D.G., K.F., and S.A. designed research; R.B., Y.J., W.M., C.W., J.S.T., and J.L. performed research; W.K., K.G., C.W., A.S., M.P.V., P.C., T.R., W.E.B., and K.M. contributed unpublished reagents/analytic tools; R.B., Y.J., W.M., M.J., and A.D. analyzed data; R.B. and S.A. wrote the paper.
Little is known about chromosomal loopings involving proximal promoter and distal enhancer elements regulating GABAergic gene expression, including changes in schizophrenia and other psychiatric conditions linked to altered inhibition. Here, we map in human chromosome 2q31 the 3D configuration of 200 kb of linear sequence encompassing the GAD1 GABA synthesis enzyme gene locus, and we describe a loop formation involving the GAD1 transcription start site and intergenic noncoding DNA elements facilitating reporter gene expression. The GAD1-TSS-50kbLoop was enriched with nucleosomes epigenetically decorated with the transcriptional mark, histone H3 trimethylated at lysine 4, and was weak or absent in skin fibroblasts and pluripotent stem cells compared with neuronal cultures differentiated from them. In the prefrontal cortex of subjects with schizophrenia, GAD1-TSS-50kbLoop was decreased compared with controls, in conjunction with downregulated GAD1 expression. We generated transgenic mice expressing Gad2 promoter-driven green fluorescent protein-conjugated histone H2B and confirmed that Gad1-TSS-55kbLoop, the murine homolog to GAD1-TSS-50kbLoop, is a chromosomal conformation specific for GABAergic neurons. In primary neuronal culture, Gad1-TSS-55kbLoop and Gad1 expression became upregulated when neuronal activity was increased. We conclude that 3D genome architectures, including chromosomal loopings for promoter-enhancer interactions involved in the regulation of GABAergic gene expression, are conserved between the rodent and primate brain, and subject to developmental and activity-dependent regulation, and disordered in some cases with schizophrenia. More broadly, the findings presented here draw a connection between noncoding DNA, spatial genome architecture, and neuronal plasticity in development and disease.
Psychosis and defective cognition in some cases with schizophrenia and bipolar disorder are thought to be rooted in molecular and cellular defects of the cortical inhibitory (GABAergic) circuitry, resulting in altered synchronization of electrical activity across widespread brain regions (Fatemi et al., 2005; Zhang et al., 2008; Lewis et al., 2012; Uhlhaas and Singer, 2012). The empirical framework for this hypothesis is based, among others, on a large body of postmortem literature reporting downregulated RNA and protein levels for the rate-limiting GABA synthesis enzyme, 67 kDa glutamic acid decarboxylase (GAD67), in multiple subtypes of cortical interneurons, including those commonly defined as fast-spiking, parvalbumin+ (PV+), and low-treshold spiking somatostatin+ cells (Akbarian and Huang, 2006; Fung et al., 2010; Konradi et al., 2011a; Schmidt and Mirnics, 2012). There is evidence that 30% or 40% of subjects with schizophrenia are affected by a more robust deficit in cortical gene expression, affecting GAD67, the transcription factor LHX6, and various neuropeptides selectively expressed by interneurons (Volk et al., 2012).
Of note, single nucleotide polymorphisms in the proximal GAD1 promoter confer genetic risk for schizophrenia, impaired working memory performance, and accelerated loss of gray matter (Addington et al., 2005; Straub et al., 2007), possibly in conjunction with altered expression of the cation chloride cotransporters NKCC1/KCC2, two key regulators of postsynaptic GABAA receptor-mediated currents (Hyde et al., 2011). These findings, in conjunction with altered chromatin structures at the proximal GAD1 promoter, including a shift from open toward repressive chromatin-associated histone modifications and abnormal DNA methylation signatures (Huang et al., 2007; Tang et al., 2011; Grayson and Guidotti, 2012), firmly point to transcriptional dysregulation as a key factor for GABAergic dysfunction in schizophrenia. However, the regulation of gene expression in a vertebrate cell goes far beyond the genetic and epigenetic architectures of proximal promoters and transcription start sites (TSSs). Instead, chromosomes and gene expression units inside the cell nucleus are organized as highly complex dynamic 3D structures that includes chromosomal loopings via direct physical interactions of enhancer and promoter elements with TSSs (Deng and Blobel, 2010; Ribeiro de Almeida et al., 2011; Singh et al., 2012). To date, however, the regulation of supranucleosomal and higher-order chromatin structures, beyond the level of DNA methylation and post-translational modifications of the core histones (which together with 146 bp of DNA wrapped around them comprise the nucleosome as the fundamental unit of chromatin), essentially remains unexplored in normal or diseased human brain.
Here, we explore 3D chromosomal architectures encompassing 200 kb of chromosome 2q31 in prefrontal cortex (PFC) of subjects with schizophrenia and controls, cultured human and mouse (chromosome 2qC2) neurons and their pluripotent precursors, and in transgenic animals tagged for GABAergic chromatin. We identify, for the first time, conserved regulatory DNA elements in rodent and primate genomes that are engaged in a loop formation with the GAD1/Gad1 TSS and present multiple lines of evidence that dynamic regulation of 3D genome structures plays a critical role for orderly development and function of GABAergic systems in the human brain.
Fibroblasts from a 43-year-old (hDF6, ATCC) healthy female were reprogrammed using a published protocol (Maherali et al., 2008) with minor modifications. Briefly, passage 6 cells were reprogrammed by infection of 3 × 105 fibroblasts with a mixture of lentivirus particles (Openbiosystems) containing genomes encoding OCT4, SOX2, KLF4, and human c-MYC at a multiplicity of infection of 5:5:5:2.5 (hDF6). Infected cells were grown for 4 d in fibroblast growth medium (DMEM/10% FCS), but 3 d after infection, each 100,000 fibroblasts were passaged into 10 cm tissue culture dishes and grown on feeder cells until candidate colonies were harvested. On day 4, the fibroblast growth medium was replaced with human ES cell medium (80% DMEM/F12) (Invitrogen), 20% knock-out serum replacement (Invitrogen), 1× nonessential amino acids (Invitrogen), 1 mm glutamine (Invitrogen), 100 μm β-mercaptoethanol (Invitrogen), and 5 ng/ml fibroblast growth factor 2 (FGF2) (R&D Systems). The culture medium was replaced daily until all candidate colonies were harvested. This approach yielded colonies with embryonic stem cell morphology at a frequency of 8 × 10−5 (data not shown). From day 30 onwards, candidate iPS cell colonies were harvested using glass knives (Vitrolife), transferred into culture wells of 6-well plate containing fresh mouse embryonic fibroblast feeder cells and expanded to colonies with a diameter of 1–2 mm. Candidate iPS cells with a modal chromosome number of 46 were then maintained further feeder free conditions on Matrigel in the presence of mTESR1.
Pluripotency was determined initially by immunofluorescent detection of endogenous expression of Nanog and Tra1-60. Surface pluripotency marker expression for Tra1–60 and SSEA4 was confirmed using fluorescence activated cell sorting. We next confirmed by qRT-PCR the expression of OCT4 and REX1 mRNAs. We assessed pluripotency further through targeted differentiation into all three germlayers and detection of Sox17 (endoderm), Otx2 (ectoderm), or Brachyury (mesoderm) protein expression by immunofluorescence. Finally, we also produced embryoid bodies from our iPSCs and assessed expression of Nestin and Neuropilin 2 (neuroectoderm), SOX17 and Alpha-Feto Protein (endoderm), and KDR, Matrilin 1, Myf5, and Brachyury (mesoderm). Line W(K)6 iPS was grown feeder free on Matrigel (BD Biosciences; 356231) and collected with dispase, and embroid bodies were formed using Aggrewell 400 plates (Stem Cell Technologies; 27845). After growing the iPS cells for 1–2 d, embroid bodies were collected and moved into nonadherent flasks, with HES media changed daily for 4 d and on the fifth day changed to neural induction (NI) media containing FGF2. On day 6, the aggregates were plated on coated tissue culture plates (poly-l-ornithine/laminin) in NI media + FGF2 and processed similar as described above for H9ESC. For chromosome conformation capture (3C) assays, 8–15 × 106 stem cells (or differentiated neurons, see below) were used as input.
Depending on the strength of the chromosomal interactions and sensitivity of the PCR assays, conventional chromosome conformation capture (3C) often requires 106 to 107 nuclei as input (Miele and Dekker, 2009). Therefore, we modified existing protocols to induce neural differentiation in a larger number of cells. Briefly, the W6 iPS were grown feeder free on Matrigel (BD Biosciences; 356231) and collected with dispase, and embroid bodies were formed using Aggrewell 400 plates (Stem Cell Technologies; 27845). After growing the iPS cells for 1–2 d, embroid bodies were collected and moved into nonadherent flasks, with HES media changed daily for 4 d and on the fifth day changed to NI media containing FGF2. On day 6, the aggregates were plated on coated tissue culture plates (poly-l-ornithine/laminin) in NI media + FGF2, changing the media every second day. By day 11–14, neural rosettes were forming. At approximately day 14, these primary neural rosettes were mechanically removed from the plate and grown in suspension for 2 d in a nonadherent flask, then replated in NI + FGF2. After 2 d, the media is changed to a modified (by the addition of N2 supplement) neural proliferation media (NP) + FGF2 + leukemia inhibitory factor (LIF) (Dhara et al., 2008). The cells are grown for a further 6–10 d in changing the media every second day. At approximately day 24–28, the secondary neural rosettes were collected and grown for 2 d in suspension and then replated on (poly-l-ornithine/laminin) in NP media + FGF2 + LIF with the addition of 200 μm Sonic hedgehog for 24 h, grown for another 7–12 d. Then, tertiary neural rosettes were collected in 1 × HBSS, incubated for 10–30 min, and then titrated gently to break up the rosettes, and the cells replated in NP media + FGF2 + LIF with the addition of 200 μm Sonic hedgehog for 24 h. After 1–2 d, the cells are grown in neural differentiation media (Sara et al., 2005) for 6–12 d and then harvested.
Fibroblasts (from two unrelated donors, including hDF6) were grown in batches of four T175 flasks in DMEM, 15% Hyclone serum + 0.5 ml glutamine and 1 × nonessential amino acids. After confluency (day 6), cells were harvested and processed for 3C as above.
For the clinical studies, specimens from the rostral third of the dorsolateral prefrontal cortex of 10 subjects diagnosed with schizophrenia (4 females/5 males, age range 40–87 years; tissue pH, 6.1–7.2) and 7 controls (2 females, 6 males; age range 41–90 years; tissue pH, 6.1–7.2) were included in this study. An independent set of 6 controls (4 females, 2 males; age range 16–81 years), and RNA integrity number from 4.1 to 9.0 was used for RNAseq experiments. Specimens were from a brain bank at the University of California at Irvine and Davis. Procedures for tissue collection, neuropathologic examination (to rule out degenerative and neurologic disease), diagnosing schizophrenia using DSM-IV-based diagnostic criteria were described previously (Akbarian et al., 1995; Huang et al., 2007).
Transgenic mice were generated by injection of circular modified bacterial artificial chromosome (BAC) into fertilized C57BL/6 mouse oocytes by the University of California Irvine Transgenic Mouse Facility (http://www.research.uci.edu/tmf/), using BAC RP23 27D24, containing the mouse glutamate decarboxylase 2 (mGad2) locus (chromosome 2: 22,453,784–22,682,644, NCBI Build 37.1) (http://bacpac.chori.org/), including the mGad2 gene (chromosome 2: 22,478,000–22,550,009, + strand). BAC modifications were performed to insert H2BGFP-FRT-neo-FRT at the Gad2 start codon, and then the selection marker (neomycin) cassette was removed via FRT directed recombination, the BAC modifications further confirmed with restriction mapping and sequencing.
All animal experiments were approved by the Animal Use and Care Committees of the participating institutions. Primary hippocampal cultures were prepared from early postnatal (P0-P1) mouse as described previously (Brewer et al., 1993) with some modifications (Futai et al., 2013). Briefly, hippocampi were dissected and trypsinized, after which neurons were dissociated and plated onto coverslips coated with poly-d-lysine (40 μg/ml) and laminin (4 μg/ml) at a density of 500 cell/mm2. Neurons were maintained in B27 containing medium for 14 d in vitro and treated with picrotoxin (100 μm), tetrodotoxin (1 μm), or equal volume of DMSO for 15 h. Total RNA was extracted using RNAqueous Micro Kit (Ambion) and reverse transcribed using High Capacity RNA-to-cDNA Kit (Applied Biosystems). Real-time quantitative PCR was then performed on an ABI 7500 Fast sequence detection system (Applied Biosystems) using TaqMan assay (Applied Biosystems). Gene expression levels in drug-treated groups or control groups (DMSO-treated) were analyzed in 3 independent cultures in PCR duplicates. Glyceraldehyde-3-phosphate-dehydrogenase (Gapdh) was used as reference gene, and fold change of Gad1 expression level in drug-treated groups relative to control group was determined using the 2−ΔΔCT method [ΔΔCT values = (CTgad1 − CTgapdh)treated − (CTgad1 − CTgapdh)control](Livak and Schmittgen, 2001).
For antipsychotic drug studies, adult male C57BL/6J mice, 10–15 weeks of age, were treated acutely or for 21 d with once daily intraperitoneal injections of saline or haloperidol (0.5 mg/kg) or clozapine (5 mg/kg) (Sigma). Tissue (PFC) was harvested 60 min after the last treatment.
Gad2-H2BGFP mice were anesthetized and perfused transcardially with cold 4% PFA in 1× PBS. Brains were harvested and postfixed in 4% PFA overnight and then cryoprotected in 30% sucrose in 0.1 m PBS. Brains were sectioned 30 μm thick on a freezing microtome. Immunostaining was done with free-floating sections using Tyramide-Plus signal amplification system. Briefly, sections were incubated with primary antibodies against calretinin (1:1000, SWANT) or parvalbumin (1:1000, SWANT), followed by biotin-tagged secondary antibody incubation (1:200; Jackson ImmunoResearch Laboratories), ABC (Vector Laboratories) treatment, and Tyramide-Plus signal amplification (1:50, PerkinElmer Life and Analytical Science). Sections were mounted to slides using VECTASHIELD mounting medium with DAPI.
Cerebral cortex (300–500 mg of gray matter/specimen) was homogenized in 1 × PBS-buffered, 1.5% formaldehyde to cross-link higher-order chromatin (of note, based on flow cytometry of purified nuclei from postmortem specimens, there are at least 3–4 × 106 physically intact nuclei in a 300 mg piece of cortical gray matter) (Jiang et al., 2008). After 15 min, cross-linking was stopped by the addition of glycine to 0.125 m final concentration) and samples spun in a tabletop centrifuge c at 4000 rpm for 10 min at 4°C; then the pellet was resuspended in 2 ml of cellular lysis buffer (10 mm Tris-HCl, pH 8.0/10 mm NaCl/0.2% Ige cal CA630) and incubated at room temperature for 45 min. Next, samples were homogenized by constant pipetting and spun at 5000 rpm for 5 min, the supernatant removed and the pelleted nuclei stored at −80°C.
Media was aspirated, followed by a quick wash with 37°C 1 × PBS. Then 5 ml of 1 × PBS/1.5% formaldehyde was added to the cultures, the cells scrapped, and kept in suspension for 10 min. Then cross-linking was stopped by the addition of glycine to 0.125 m and cells harvested by centrifugation at 1200 rpm/5 min. The pellet is treated with 1 × cellular lysis buffer after adding one-tenth volume of protease inhibitor mixture (Sigma), resuspended, and left at room temperature for 15 min, followed by homogenization and nuclei pelleting (5000 rpm for 5 min) and two washes with 1 × New England Biolabs2 [NEB2] restriction enzyme buffer).
Nuclei were resuspended in 250 μl of 1 × NEB2, divided into 5 × 50 μl aliquots, to each of which 312 μl of 1 × NEB2 and 38 μl of 1% SDS were added, followed by gentle mixing and incubation at 65°C for 10 min to remove proteins that are not directly associated with DNA and facilitate more efficient digestion. Then, 44 μl of Triton X-100 was added to each sample to quench the SDS by gentle pipetting, and samples were digested with 400 units of HindIII (NEB) at 37°C overnight with gentle shaking.
To inactivate HindIII, 86 μl of 10% SDS (Bio-Rad) was added and samples were incubated at 65°C for 30 min. Then, each sample was added to 7.61 ml of ligation reaction mixture. Each mixture consisted of 745 μl of 10% Triton X-100, 745 μl of 10 × ligation buffer [1 m Tris HCl, pH 7.5, 1 m MgCl2, 1 m DTT dithiothreitol (Bio-Rad)], 80 μl of 10 mg/ml bovine serum albumin (NEB), 80 μl of 100 mm ATP (Sigma) and 5960 μl of H2O. To each mixture, 50 μl of T4 DNA ligase (1 U/μl; Invitrogen) was added and incubated at 16°C for 4 h followed by reverse cross-linking at 65°C overnight with simultaneous Proteinase K digestion [50 μl of 10 mg/ml Proteinase K (Sigma) reconstituted with TE buffer pH8.0] to digest chromatin-associated proteins. To improve ligated DNA recovery, another 50 μl of Proteinase K was added and incubated at 65°C for 2 h.
Ligated samples were phenol extracted (pH 8.0, Fisher) followed by phenol-chloroform (1:1) (pH 8.0, Omnipur) followed by ethanol precipitation. DNA was dissolved in 1 × TE buffer (pH 8.0).
For 3C studies in Gad2-H2BGFP transgenic mice, 0.2–0.3 × 106 GFP+ nuclei and 48 × 106 GFP− nuclei were separately sorted and collected from cerebral cortex using fluorescence-activated nuclei sorting as described previously (Jiang et al., 2008). Additional 3C studies were conducted in whole cortex of adult C57BL/6J wild-type mice or on the primary neuronal cultures described above.
We mapped physical interactions of noncontiguous DNA elements within a 200 kb region of chromosome chromosome 2 (171,573,000–171,797,000; HG19), encompassing GAD1 (Fig. 1) by PCR across the ligation junctions of interacting fragments. Primers were designed according to 3C guidelines (Miele and Dekker, 2009), with 30 bp in length and positioned within 60 bp of the restriction enzyme cut site, a melting temperature of 60°C, and with a GC% content of 45%. From a total of 55 HindIII fragments, 11 fragments yielding primer sequences in accordance with these primer selection criteria were probed by specific primers positioned in the 5′ to 3′ direction on the sense strand amplifying toward the end of the respective restriction fragment were used. The anchor primer for all PCR represented the HindIII fragment containing the TSS of GAD1. All primers were 30–32 nt in length, melting temperature = 60°C, and positioned < 200 bp from a HindIII restriction site. PCR products were resolved in agarose gels and the level of interaction between any two fragments in the 3C libraries measured semiquantitatively. To this end, amplicons generated by PCR across a single ligation junction run in triplicate on 2% agarose/1xTAE/130V/25 min and band intensities measured by UVP Bioimaging system/Labworks version 4.5 software, after background subtractions from blank gel regions. We used a BAC (clone ID-RP11–570C16, representing chromosome 2: 171552700–171745852; http://bacpac.chori.org) to control for primer efficiency (Dekker, 2006; Miele and Dekker, 2009). To control for ligation efficiencies across independent 3C experiments/samples, we mapped physical interactions of noncontiguous DNA elements within 75 kb of chromosome 16 (62,275,000–62,350,000; HG19) in a previously described “gene desert” region with the two nearest genes (CDH8 and CDH11) located 200 kb upstream and 2.6 Mb downstream (Ferraiuolo et al., 2010).
Furthermore, we mapped physical interactions of noncontiguous DNA elements within a 70 kb region of chromosome mouse chr2qC2 (chromosome 2: 70,340,000–70,410,069; mm9), to interrogate loop formations of DNA sequences positioned up to 100 kb upstream of Gad1 TSS.
Luciferase reporter gene assays were performed using the minimal promoter (TATA box) luciferase vector (pGL4.23 [luc2/minP] vector; Promega; E8411) in HEK-293 cell cultures. Specifically, 115 bp (chromosome 2: 171622200–171622315) was cloned into the multiple cloning site of the luciferase vector and sequence-verified clones transfected into HEK-293 cells with Lipofectamine-2000 (Invitrogen). Furthermore, control plasmids were generated with a 125 bp sequence (chromosome 2: 171645000–171645125) from an upstream region restriction fragment that was not surrounded by H3K4me1 or H3K27ac peaks in various Encyclopedia of DNA Elements Consortium (ENCODE) cell lines and did not harbor AP-1 motifs. As an additional control, Renilla luciferase vector (Prl-TK; Promega; E2241) was cotransfected into every well to correct for variations in transfection efficiency. Cells were harvested and lysed 48 h after transfection using lysis buffer from the Promega luciferase assay kit (catalog #E1910) and firefly luciferase units quantified (kit luciferase substrate) using the Bio-Rad luminometer followed by measurement of Renilla luciferase units by chemiluminescence (kit Stop& Glo reagent) according to the manufacturer's protocol. Relative luciferase units were calculated relative to minimal TATA box driven luciferase units.
Cross-linked nuclei were restriction digested as described in the 3C method using HindIII (NEB) enzyme. Digested restriction fragments were then washed with FSB (20 mm Tris-Cl, pH 7.5, 5 mm EDTA, 50 mm NaCl final concentration) followed by immunoprecipitation (anti-H3K4me3 1 μg/μl; 1:1000 dilution; ab8898, Abcam) at 4°C overnight in FSB with phenyl methyl sulfonyl fluoride (0.1 μm final concentration), dithiothreitol (3 μm final concentration), and benzamidine (0.1 μm final concentration). Antibody–chromatin complexes were then isolated using protein G agarose slurry followed by elution as described previously (Huang et al., 2007). The immunoprecipitated chromatin fragments were then ligated with T4 DNA ligase (Invitrogen) as in the 3C protocol step. Ligated fragments were then reverse cross-linked and digested with Proteinase K (Invitrogen, catalog #25530–031) and purified to generate the 3C-ChIP loop library using phenol/phenol-chloroform extraction followed by ethanol precipitation as in 3C. Quantification of interacting fragments from H3K4me3 ChIP-3C, and for control, IgG-3C and H3K4me3 ChIP-3C (no ligase) libraries was done exactly as for 3C quantification from conventional 3C assays described above. The H3K4me3 ChIP-3C, IgG-3C, and H3K4me3 ChIP-3C (no ligase) assays and libraries were processed in parallel using the same protocol and solutions.
Chromatin from cerebral cortex from one hemisphere of adult C57BL/6J mice (N = 3) was prepared by formalin crosslinking and sonication as described previously (Jiang et al., 2010), for ChIP with an antibody raised against the C terminus of the cAMP-response element binding protein (CREB)-binding protein (CBP/CREBBP) (Abs c-583, rabbit polyclonal, Santa Cruz Biotechnology). For PCR-based quantifications of ChIP assays, input DNA and immunoprecipitated DNA were amplified by SYBR Green-based real-time PCR. Primers were designed to interrogate four positions in the Gad1 locus, with the following chromosome 2 (mm9) coordinates: (1) 70330778–70330906, −77 kb upstream of the Gad1 TSS; (2) 70352784–70352899, −55 kb to Gad1-TSS; (3) 70399957–70400075 at Gad1-TSS; and (4) 70437331–70437444, +37 kb downstream from Gad1-TSS.
RNA was extracted from ~75–100 mg of gray matter of 6 adult control PFC specimens with Rneasy Lipid Tissue mini kit (catalog #74804, QIAGEN), treated with DNase I, purified, and diluted to 20 ng/μl. Sequencing libraries were prepared according to the NuGen Ovation RNASeq version 2 protocol, and run on the paired-end 50 bp module in Illumina HiSeq 2000 (Eurofins MWG; Operon). RNA-seq raw reads that passed the QC metric, which is called as the “chastity filter” by Illumina, were aligned to the UCSC Homo sapiens reference genome build 19 using TopHat version 1.4.0 default parameters, allowing a maximum of 2 mismatches when mapping reads to the reference and multiple alignments per read (up to 40 by default). Raw reads that passed Illumina's chastity filter varied from 22,470,801 to 60,734,424. TopHat uses the bowtie version 0.12.8 algorithm to perform the alignment in first phase and then identifies potential splice junctions by using these mapping results and aligns previously unmapped reads in second phase to generate alignment files. The resulting aligned bam files were processed by SAMtools (Li et al., 2009), and the read coverage of indexed alignment files for each sample was visualized with the Integrative Genomic Viewer. For intranuclear RNAseq, FACS-sorted neuronal NeuN+ nuclei (Jiang et al., 2008) (>3 × 106) from adult C57BL/6J cerebral cortex were subjected to total RNA isolation using trizol (Invitrogen) reagent. Ribsosomal and other RNAs <200 nt were removed using RNA clean and Concentrator kit (Zymo Research) followed by DNA digestion. A total of 20 ng/μl (10 μl total = 200 ng) RNA was fragmented at 95°C for 5 min using Invitrogen fragmentation buffer (from SuperScript III). First-strand synthesis was performed using Random primer (Invitrogen; catalog #48190011) and Superscript III. First-strand cDNA was cleaned using Ampure XP beads followed by second-strand synthesis using dUTP mixture (dUTP 20 mm, dATP, dCTP, dGTP 10 mm) and DNA polymerase I (10U/μl, NEB). Double-stranded cDNA was end-repaired using Klenow DNA polymerase END-repair kit (Epicenter) followed by A-tailing and adapter ligation. Ampure XP beads were again used to purify the adapter ligated double-stranded DNA followed by Uracil-DNA Glycosylase (NEB, M0280S) incubation to remove dUTP containing second strand. Strand-specific single-stranded DNA was again purified using Ampure XP beads followed by PCR using adapter specific barcoded PCR primers. Purified library was confirmed using bioanalyzer high-sensitive DNA ChiP analysis followed by deep-sequencing using Illumina multiplexing protocol and bioinformatical analyses as described above.
Data comparing two groups were analyzed with two-tailed t tests, paired or unpaired, where appropriate. Data on GAD1 RNA expression in PFC versus fibroblasts (FIB) was analyzed using the Wilcoxon signed rank test. Data with more than two groups were analyzed by one-way ANOVA having either “locus” or for the picrotoxin/tetrodoxin-based neuronal activity experiments “treatment” as between group factors, followed by Bonferroni′s multiple-comparison test. The 3C data from mouse cortex of animals subjected to systemic administration of saline versus haloperidol or clozapine were analyzed by two-way ANOVA having “locus” and “treatment” as between group factors, followed by Tukey′s honest significant difference test.
Cortical gray matter from the rostral frontal lobe of four adult postmortem specimens, ranging in (1) age from 58 to 81 years, (2) postmortem interval 7–11 h, and (3) tissue pH 6.1–7.2 was used to prepare 3C libraries from DNA ligase-treated HindIII digests, to map and quantify physical interactions of noncontiguous DNA elements by PCR. Together, we explored 200 kb on chromosome 2 (171,573,000–171,797,000; HG19) surrounding the GAD1 TSS (Fig. 1A).
Physical interaction frequencies from our PFC 3C assays are shown in Figure 1B. Because chromatin fibers are flexible, resulting in much higher interactions of neighboring restriction fragments within 20–25 kb range in 3C assays (Gheldof et al., 2006; Miele et al., 2006), we primarily interrogated physical interactions of GAD1 TSS with DNA sequences that on the linear genome are positioned >25 kB from the TSS. With the anchor primer on the 11.6 kb HindIII fragment harboring GAD1 TSS (Fig. 1B, green box), the 3C assays from the PFC consistently revealed in 4 of 4 specimens a robust interaction with a restriction fragment positioned −50 kb upstream of the TSS (red box in Fig. 1B). This interaction was highly significant in the PFC (F(9,30) = 3021, p < 0.001 with Bonferroni's post hoc for multiple comparisons) and tissue specific because no evidence for this chromosomal looping was found in 3C assays from cultured fibroblasts from two unrelated donors (Fig. 1B). In contrast to these brain/fibroblast differences in GAD1 higher-order chromatin, the physical interaction of noncontiguous DNA elements across a 75 kb intergenic portion on chromosome 16 were similar between PFC and fibroblasts (mean ± SEM: PFC, 0.30 ± 0.02; fibroblasts 0.32 ± 0.01). Importantly, no PCR product was observed in PFC 3C assays from the same samples when processed without the critical DNA ligase step (which conjugates noncontiguous DNA elements) before phenol extraction and DNA purification (Fig. 1B). These findings, together, strongly suggest that chromosome 2q31 in human PFC, but not fibroblasts, includes a loop formation (referred to as GAD1-TSS-50kbLoop hereafter) between the GAD1 TSS with intergenic sequences positioned ~50 kb further upstream.
We noticed that both the GAD1 TSS and its looping partner sequences positioned 50 kb upstream of GAD1 were epigenetically decorated with sharp peaks for histone H3 trimethylated at lysine 4 (H3K4me3) in neuronal chromatin extracted from prefrontal cortex (Cheung et al., 2010) (Fig. 1A). The H3K4me3 modification is an epigenetic mark sharply regulated around TSSs and other regulatory sequences associated with open chromatin and active gene expression (Guenther et al., 2005, 2007). In human cerebral cortex, the H3K4me3 mark shows a preference for CpG enriched sequences, specifically those not tagged by DNA cytosine methylation and repressive histone marks (Maunakea et al., 2010). The H3K4me3 peak at the GAD1 TSS (Fig. 1A) extended across 2.5 kb (HG19, chromosome 2: 171,672,200–171,674,500) and included all alternative gene proximal promoters previously associated with GAD1 expression (Chen et al., 2011). Furthermore, the −50 kb upstream sequence was defined by sharp peaks of histone H3 acetylated at lysine 27 (H3K27ac) in conjunction with a corresponding spike of histone H3 monomethylated at lysine 4 (H3K4me1) in ENCODE ChIP-seq tracks for H1 human embryonic stem cells (H1ESC), and the human umbilical vein endothelial and K562 erythroblastoid cell lines (Dunham et al., 2012) (Fig. 1A). Intriguingly, in tissue homogenates from frontal cortex and hippocampus and in cultured neurons, the 50 kb upstream fragment interacting with GAD1 TSS showed similar types of enrichment for H3K4me1 and H3K27 (Fig. 2A, browser tracks were loaded with data from Zhu et al., 2013). Of note, such sharp peaks of H3K27ac in conjunction with H3K4me1 and H3K4me3 often define DNA elements that function as enhancers and regulate gene expression at sites of promoters positioned further upstream or downstream (Barski et al., 2007; Maston et al., 2012). Therefore, we hypothesized that GAD1-TSS-50kbLoop is associated with the open chromatin fraction from the PFC, specifically with the H3K4me3-tagged GAD1 TSS interacting with nucleosomes positioned 50 kb further upstream that are epigenetically decorated with the same mark. In support of this hypothesis, GAD1-TSS-50kbLoop was readily detectable in PFC, a tissue with robust levels of GAD1 gene expression, in contrast to fibroblasts, which lack both the loop and detectable levels of GAD1 gene transcript (Fig. 1B,C). To test this hypothesis directly, we used3C-ChIP Loop (Simonis et al., 2007), which is a modified version of the 3C approach to study chromosomal loopings in selected chromatin fractions enriched by immunoprecipitation. To this end, we immunoprecipitated HindIII digested chromatin from human PFC with anti-H3K4me3 antibody followed by PCR to interrogate ligation products representing GAD1-TSS-50kbLoop (Fig. 1D). Indeed, pulldown with anti-H3K4me3 antibody resulted in much stronger enrichment for GAD1-TSS-50kbLoop compared with pulldown with IgG control; furthermore, no PCR products from noncontiguous DNA sequences were observed in 3C-ChIP Loop experiments without ligase (Fig. 1D). We conclude that GAD1-TSS-50kbLoop is a higher-order chromatin structure in which at least one or both sides of the loop are enriched with the H3K4me3 mark.
As discussed above, the loop formation in chromosome 2q31, resulting in a physical interaction between the sequences surrounding the GAD1 TSS with DNA positioned 50 kb further upstream, is associated with GAD1 gene expression and open chromatin-associated histone methylation and acetylation in human PFC. Of note, the interacting DNA upstream of GAD1 is associated with sharp histone modification peaks associated with enhancer functions, including H3K27ac and H3K4me1 in peripheral cell lines and H3K4me3 in PFC neurons, but in the absence of detectable (above background) cis-associated transcript at the peak site (Fig. 1A, red box; Fig. 2A) in our deeply sequenced total RNA extracts from PFC of six control subjects (Fig. 2B). These findings, together, would suggest that the upstream components of GAD1-TSS-50kbLoop function most likely as enhancers by facilitating transcription when brought in proximal position to a gene promoter. To explore this, we focused on 115 bp (chromosome 2: 171622200–171622315, HG19) which, based on TFSearch (http://www.cbrc.jp/research/db/TFSEARCH.html) (Heinemeyer et al., 1998), was highly enriched for a total of seven activating protein AP-1 early response transcription factor motifs (TGAG/CTCA), reflecting a much higher AP-1 motif density, or 60 AP-1 sites/kb, compared with the surrounding 15 kb with only 1.2 AP-1 sites/kb. The AP-1 motif is thought to mediate a significant portion of activity-dependent gene expression both in human and rodent PFC (Covington et al., 2010). We fused this 115 bp sequence 5′ from a minimal promoter sequence and explored changes in transcriptional activity with a luciferase assay in HEK-293 cells. Indeed, in 3 of 3 experiments, addition of the 100 bp element resulted in a robust, >10-fold increase in luciferase expression and activity, compared with a control sequence void of AP-1 motifs (Fig. 3). These findings support the hypothesis that GAD1-TSS-50kbLoop plays a role in the regulation of GAD1 transcription, by positioning upstream enhancer sequences in close proximity to GAD1 TSS.
The above experiments suggest that physical interactions of open chromatin-associated DNA elements at the GAD1 locus, including GAD1-TSS-50kbLoop, are tissue specific and readily detectable in PFC. Of note, PFC is defined by robust expression of GAD1 RNA, in contrast to skin fibroblasts, which lack both detectable levels of gene expression and a tractable 3D chromosomal architecture at this locus (Fig. 1B,C). To further test the association between GAD1-TSS-50kbLoop and active gene expression, we studied GAD1 higher-order chromatin and gene expression both in induced pluripotent stem cells (which were generated from one of the two fibroblast cultures used for 3C assays shown in Fig. 1B,C) and during their transition into differentiated neural cultures primarily comprised of microtubule-associated protein 2-immunoreactive neurons (Fig. 4A,B). Indeed, GAD1-TSS-50kbLoop, but no other DNA elements across a 220 kb segment of chromosome 2q31, showed a differentiation-induced, significant increase in interaction with the GAD1 TSS (Fig. 4B,C), in conjunction with increased GAD1 RNA expression (Fig. 4D). Levels of GAD1-TSS-50kbLoop in the pluripotent stem cells were ~50% lower compared with the neural cultures differentiated from them (Fig. 4D), which would be consistent with the notion that pluripotent cells maintain some degree of open chromatin, including histone acetylation and H3K4-trimethylation, at genes poised for expression during subsequent development and differentiation (Vastenhouw and Schier, 2012). The observed increase in GAD1-TSS-50kbLoop structures upon neural differentiation was highly specific for this type of higher-order chromatin because no loop structures were observed in 3C assays from cell cultures without the ligation step (Fig. 4C). Furthermore, the physical interactions of noncontiguous DNA elements across a 75 kb intergenic portion on chromosome 16, which is insensitive to cellular differentiation (Ferraiuolo et al., 2010), were very similar between stem cells (0.31 ± 0.01) and differentiated neural cultures (0.33 ± 0.01).
We noticed that, in mouse cerebral cortex, histone H3K4me3 landscapes around the chromosome 2qC2 Gad1 locus (Fig. 5A) show a striking resemblance to the epigenetic profile encountered in neuronal chromatin from human (Fig. 1A). Both mouse and primate cortex exhibits sharp H3K4me3 peak profiles both at the GAD1/Gad1 TSS and a conserved CpG-rich sequence positioned −55 kb (−50 kb) upstream of Gad1 (GAD1) gene. The CpG-rich sequences at −50 kb from the GAD1 TSS (HG19, chromosome 2: 171,627,226–171,628,062) include a 302 bp sequence that shows 88% sequence homology to mouse (mm9, chromosome 2: 70,347,243–70,347,545). This degree of conservation is much higher than what is expected, given that the 10 kb of surrounding sequence shows a much lower degree, or 47.2%, of sequence conservation (mm9 chromosome 2: 70342243–70352545 to hg19 chromosome 2: 171621775–171629528). Furthermore, as for human PFC (tissue homogenate) (Fig. 2B), RNAseq from neuronal nuclei of mouse cerebral cortex did not reveal evidence for cis-associated transcript at the H3K4me3 peak site positioned −55 kb upstream of Gad1 (mouse chromosome 2qC2) (Fig. 5A, red box). Furthermore, inspection of the H3K4me3 ChIp-seq datasets from PFC neurons and nonhuman primate prefrontal cortex, including macaque and chimpanzee (Shulha et al., 2012b), revealed epigenetic profiles at the Gad1 locus very similar to those observed in human and mouse, with H3K4me3 peaks limited to the Gad1 TSS and the conserved CpG-rich sequence positioned ~50 kb further upstream: Rhesus macaque genome RM3, chromosome 12: 34110698–34112216 and chimpanzee (Pan troglodytes) genome PT3 chromosome 2B: 175198622–175200113. Given the conserved H3K4me3 landscapes between mouse and primate brains, we hypothesized that 3D genome architectures in mouse cerebral cortex will include a homolog to the aforementioned GAD1-TSS-50kbLoop from the human PFC, with a conserved enhancer sequence positioned 55 kb upstream of Gad1 directly interacting with the Gad1 promoter. Indeed, 3C assays from mouse cerebral cortex consistently revealed highest levels of interaction between the Gad1 TSS and the HindIII restriction fragment harboring the H3K4me3-decorated, CpG-rich sequence positioned ~55 kb further upstream (Fig. 5A). We conclude that GAD1-TSS-50kbLoop is not unique to (human) chromosome 2q31 because very similar chromatin architectures, including Gad1-TSS-55kbLoop, are observed in mouse brain.
Of note, GAD1-TSS-50kbLoop in human is lacking in fibroblasts but highly regulated in brain where it is tightly linked to open chromatin and active expression of the GAD1 gene (Fig. 1B–D) and upregulated during the course of neuronal differentiation (Fig. 4). Therefore, we hypothesized that this chromosomal loop is a defining feature of GABAergic neurons, which in brain are the only cell type expressing high levels of GAD1/GAD67 GABA synthesis enzyme (and of a closely related gene, GAD2/GAD65). To address this question, we first generated BAC transgenic mice expressing GFP-conjugated histone H2B under control of the Gad2 gene promoter. We confirmed widespread expression of the transgene in cortical, hippocampal, and striatal GABAergic neurons (Fig. 5B) and then conducted 3C assays separately for GFP+ and GFP− nuclei collected by fluorescence-activated nuclei sorting from cerebral cortex and hippocampus of adult Gad2-H2BGFP mice. Strikingly, Gad1-TSS-55kbLoop was highly specific for the population of GABAergic (GFP+) nuclei and completely absent in chromatin of non-GABAergic cells (Fig. 5C). This effect was consistently observed in 3 of 3 experiments. Therefore, Gad1-TSS-55kbLoop, a higher-order chromatin structure associated with Gad1 gene expression, is part of the genome architecture specific for GABAergic neurons.
Activity-dependent regulation of Gad1/Gad67 gene expression is important for the orderly balance of inhibition and excitation in the cerebral cortex both in rodent and primates. For example, in monkeys, monocular deprivation of sensory input into visual cortex results after 5 d in a strong and prolonged downregulation of Gad1/Gad67 RNA in the ocular dominance columns innervated by the affected eye (Benson et al., 1994). In the rat, loss of hippocampal pyramidal neurons and hippocampal–prefrontal connectivity results in downregulated Gad1/Gad67 expression in hippocampus (Shetty and Turner, 2001) and prefrontal cortex (Lipska et al., 2003). Decreased Gad1/Gad67 expression, in turn, results in reduced levels of cytosolic GABA levels and lower amplitudes of miniature inhibitory synaptic postsynaptic currents (mIPSCs) (Lau and Murthy, 2012). Given that activity-regulated fine-tuning of Gad1/Gad67 expression is of pivotal importance for homeostatic plasticity both in primate and rodent cerebral cortex, we asked whether higher-order chromatin structures around the GAD1 locus, including GAD1-TSS-50kbLoop and its mouse homolog, Gad1-TSS-55kbLoop are involved in these mechanisms. We explored this question in primary cultures from mouse hippocampus, which were exposed for 15 h either to vehicle (DMSO) or the GABAA receptor antagonist picrotoxin (PTX) (100 μm) to increase general neuronal activity, or the sodium channel blocker TTX (1 μm) for silencing. As expected, Gad1/Gad67 RNA levels were consistently increased (3 of 3 experiments) by 70% in PTX-treated cultures but decrease by ~40% after TTX, compared with vehicle controls (Fig. 5D). This was associated with a consistent (3 of 3 experiments) 50% increase in physical interactions of the Gad1 TSS with its partnering sequences in Gad1-TSS-55kbLoop after treatment with PTX, compared with vehicle control (Fig. 5D). In contrast, Gad1-TSS-55kbLoop was not affected by TTX treatment. We conclude that Gad1 higher-order chromatin is dynamically regulated upon acute changes in neuronal activity, particularly in the context of increased excitation.
Our results show that DNA elements positioned 50 (human) or 55 (mouse) kb upstream of GAD1/Gad1 facilitate transcription and engage in a loop formation with the gene's TSS that is upregulated upon increased neuronal activity and transcription. Therefore, we hypothesized that both portions of the loop, including chromatin at site of the Gad1 TSS and its partnering sequences positioned −55 kb further upstream, show increased occupancy for proteins, such as the CREB-binding protein and histone acetyltransferase, CREBBP/CBP, which, like its close homolog EP300, acts as a broad facilitator of gene expression (Alarcon et al., 2004; Korzus et al., 2004; Wood et al., 2006; Kim et al., 2010). To explore this, we conducted ChIP assays on adult mouse cerebral cortex and measured CBP binding in Gad1-TSS-55kbLoop chromatin and additional positions in chromosome 2qC2. Indeed, CBP load was significantly higher at the Gad1-TSS and its partnering sequences −55 kb further upstream, compared with chromatin positioned −77 kb upstream, or +37 kb downstream from the TSS (Fig. 5E).
Subjects on the mood and psychosis spectrum, including many cases diagnosed with schizophrenia or bipolar disorder, are frequently affected by dysregulated gene expression in GABAergic interneurons in multiple brain regions (Guidotti et al., 2000; Konradi et al., 2011b). In the PFC (in particular, 30% or 40% of subjects with schizophrenia) are affected by a robust deficit in GAD1/GAD67 gene expression (Volk et al., 2012). Consistent with this observation, open chromatin-associated histone modifications, including trimethylation of H3K4 and acetylation of other lysine residues, are decreased in nucleosomes surrounding the GAD1 TSS in the PFC of a subset of cases with schizophrenia (Huang et al., 2007; Tang et al., 2009). Given that the murine homolog of GAD1-TSS-50kbLoop is specific for GABAergic neuron chromatin (Fig. 5C), and given the tight coregulation of this loop with GAD1/Gad1 expression both in and ex vivo in the human and the mouse (Figs. 1B–D, ,44D, and and55D), we hypothesized that this type of higher-order chromatin is altered in cases with schizophrenia that are affected by loss of GAD1 gene expression. To explore this, we conducted 3C assays on PFC tissue of 10 cases with schizophrenia and 7 controls, including 5 case-control pairs (matched for age/gender/autolysis time) previously found to show decreased GAD1/GAD67 RNA expression and H3K4 trimethylation in the disease case (Huang et al., 2007). Indeed, comparison of the general cohort of schizophrenics and controls, and separate analyses of the 5 case-control pairs revealed a significant (at least 35–45%) decrease in physical interaction frequencies between the GAD1-TSS and its partnering sequences positioned −50 kb further upstream (Fig. 6A,B). Consistent with the findings from the aforementioned postmortem studies that loss of transcription-associated epigenetic decorations of sequences GAD1 TSS are not caused by antipsychotic medication (Huang et al., 2007; Tang et al., 2009), no quantitative changes in the strength of Gad1-TSS-55kbLoop were observed in cerebral cortex of mice chronically treated with the conventional antipsychotic and dopamine D2-like receptor antagonist, or atypical antipsychotic clozapine (Fig. 5A).
Mapping physical interactions of noncontiguous DNA elements across 200 kb surrounding the GAD1/GAD67 GABA synthesis enzyme gene locus in human PFC, we discovered a loop formation GAD1-TSS-50kbLoop. The existence of GAD1-TSS-50kbLoop, and Gad1-TSS-55kbLooop, as its mouse homolog, provides testimony for actively regulated higher-order chromatin structures in the human brain because GAD1-TSS-50kbLoop was absent in 3C assays in which the critical DNA ligase step was omitted; furthermore, the loop was coregulated with GAD1 gene expression and therefore much stronger in cultured neurons compared with their undifferentiated precursors or the donors' skin fibroblasts. In mature cerebral cortex, the Gad genes are selectively transcribed in cortical interneurons; and indeed, our 3C assays on sorted nuclei revealed that Gad1-TSS-55kbLoop was specific for GABAergic neurons. These findings, together, suggest that GAD1-TSS-50kbLoop serves as a classical promoter/enhancer loop (Ong and Corces, 2011), by moving conserved DNA elements positioned 50 kb upstream of GAD1 into close proximity with the gene's TSS to facilitate its expression. In support of this model, the upstream interacting sequences are decorated with sharp peaks of histone H3 acetylated at lysine 27 (H3K27ac) in conjunction with monomethyl and trimethyl H3K4 (H3K4me1/3), which is of interest because co-occurrence of these histone modifications often defines enhancer elements (Barski et al., 2007; Zhou et al., 2011; Maston et al., 2012). In addition, we observed increased occupancy of the transcriptional facilitator, CBP/CREBBP, a protein that on a genome-wide scale preferentially locates both to gene proximal promoters and enhancers (Visel et al., 2009; Kim et al., 2010), at the two interacting sequences in Gad1-TSS-55kbLoop. Furthermore, we deeply sequenced in the present study total RNA from tissue homogenates of human PFC, and from mouse cortex neuronal nuclei, and observed that none of the upstream DNA elements in GAD1-TSS-50kbLoop and Gad1-TSS-55kbLoop was transcribed in the cortex. Instead, these sequences harbored AP-1 transcription factor motifs; and, as expected for an enhancer sequence, when fused to generic promoters, we observed increased expression of the luciferase reporter.
Chromosomal loopings, including the GAD1-TSS-50kbLoop described here, are among the most highly regulated supranucleosomal structures and pivotal for transcriptional regulation, by enabling distal regulatory enhancer or silencer elements positioned a few or many hundred kilobases apart from a gene, to interact directly with that specific promoter (Gaszner and Felsenfeld, 2006; Wood et al., 2010). Surprisingly, however, only very few studies explored loop formations in brain tissue (Horike et al., 2005; Dhar et al., 2009; Jiang et al., 2010; Yasui et al., 2011; Shen et al., 2012; Shulha et al., 2012b). Our study is the first to (1) monitor the emergence of higher-order chromatin at a GABAergic gene locus during the process of neural differentiation, and demonstrate that epigenetic architectures, including a specific promoter/enhancer loop, are (2) conserved between rodent and primate (including human) Gad1/GAD1, (3) sensitive to changes in neuronal activity, and (4) significantly dysregulated in PFC of subjects with schizophrenia affected by loss of GAD1 gene expression. It remains to be determined whether the weakening of prefrontal GAD1-TSS-50kbLoop in a specific disease case is the result of genetic factors (including the GAD1 risk haplotypes in selected populations) (Addington et al., 2005; Lundorf et al., 2005; Straub et al., 2007; Ayalew et al., 2012) or the outcome of a developmental defect (Benes, 2012) or altered glutamate NMDA receptor signaling (Behrens et al., 2007; Homayoun and Moghaddam, 2007; Lisman et al., 2008; Belforte et al., 2010) or simply an adaptive response to cortical dysfunction, or combination thereof. In either case, a prolonged decrease of GAD1 is likely to be detrimental for higher cognitive function and the inhibitory networks in the cortex, given that Gad1 is important for homeostatic synaptic scaling by controlling GABA release and mIPSCs (Lau and Murthy, 2012), and the orderly formation of perisomatic innervation of cortical pyramidal neurons and basket neuron innervation fields in the adolescent cortex (Chattopadhyaya et al., 2007). Furthermore, GAD1 shares genetic risk factors that lead to alterations in the expression of the cation chloride cotransporters NKCC1/KCC2, which are of pivotal importance for the fine-tuning the inhibitory/excitatory effects of postsynaptic GABAA receptor-mediated currents (Hyde et al., 2011).
In the present study, we first mapped physical interactions between noncontiguous DNA elements in postmortem brain specimens collected several hours after death, then explored these types of higher-order chromatin in cell culture and animal models. The implications of our study, when viewed as a proof of principle, are obvious: Only 1% of the human genome encodes protein. However, the vast majority of noncoding sequences are engaged in various open and repressive higher-order chromatin structures involved in genome organization and function (Lieberman-Aiden et al., 2009), but in the human brain these mechanisms until now remain completely unexplored. Because both the intergenic portions of the genome and its intronic sequences each harbor at least 40% of DNA polymorphisms and structural variants associated with complex traits and common medical (including psychiatric) conditions (Stranger et al., 2011), the findings presented here will pave the way for a more comprehensive mapping of 3D genome architectures in normal and diseased human brain. Such type of information, in conjunction with genome-wide DNA methylation and histone modification mappings (Mill et al., 2008; Kano et al., 2012; Labonte et al., 2012; Numata et al., 2012; Shulha et al., 2012a; Guintivano et al., 2013), will provide deep insights into the coordinated regulation of the epigenetic landscape of neuronal and glial genomes, including the role of DNA sequence variants, in specific disease cases diagnosed with schizophrenia and other psychiatric disease. Finally, the observation that the GAD1-associated chromosomal loop and the H3K4 trimethylation landscape associated with it are conserved between rodent and primate brain would suggest that higher-order chromatin, just as recently reported for DNA methylation, histone regulation, and transcription factor occupancies, shows robust homologies between species from different mammalian lineages (Xiao et al., 2012). Thus, the “genome in 3D” is fully open to study in human brain and preclinical models alike.
This work was supported by National Institutes of Health Grants R01 MH09332, R21NS076958, and P50MH096890, and by the Brain & Behavior Research Foundation. We thank Dr. Yin Guo for excellent technical assistance.
The authors declare no competing financial interests.