Temporal dynamics and genetic control of transcription in the human prefrontal cortex.

Previous investigations have combined transcriptional and genetic analyses in human cell lines, but few have applied these techniques to human neural tissue. To gain a global molecular perspective on the role of the human genome in cortical development, function and ageing, we explore the temporal dynamics and genetic control of transcription in human prefrontal cortex in an extensive series of post-mortem brains from fetal development through ageing. We discover a wave of gene expression changes occurring during fetal development which are reversed in early postnatal life. One half-century later in life, this pattern of reversals is mirrored in ageing and in neurodegeneration. Although we identify thousands of robust associations of individual genetic polymorphisms with gene expression, we also demonstrate that there is no association between the total extent of genetic differences between subjects and the global similarity of their transcriptional profiles. Hence, the human genome produces a consistent molecular architecture in the prefrontal cortex, despite millions of genetic differences across individuals and races. To enable further discovery, this entire data set is freely available (from Gene Expression Omnibus: accession GSE30272; and dbGaP: accession phs000417.v1.p1) and can also be interrogated via a biologist-friendly stand-alone application (http://www.libd.org/braincloud).

Characteristics of the cation cotransporter NKCC1 in human brain: alternate transcripts, expression in development, and potential relationships to brain function and schizophrenia.

Early in development, GABA, an inhibitory neurotransmitter in adults, is excitatory. NKCC1 (SLC12A2) encodes one of two cation chloride cotransporters mediating the conversion of GABA from excitatory to inhibitory. Using 3' and 5' RACE and PCR, we verified previously characterized alternative transcripts of NKCC1a (1-27) and NKCC1b (1-27(Δ21)), identified new NKCC1 transcripts, and explored their expression patterns during human prefrontal cortical development. A novel ultra-short transcript (1-2a) was expressed preferentially in the fetus. Expression of NKCC1b and 1-2a were decreased in schizophrenia compared with controls (NKCC1b: 0.8-fold decrease, p = 0.013; 1-2a: 0.8-fold decrease, p = 0.006). Furthermore, the expression of NKCC1b was associated with NKCC1 polymorphism rs3087889. The minor allele at rs3087889, associated with reduced NKCC1b expression (homozygous for major allele: N = 37; homozygous for minor allele: N = 15; 1.5-fold decrease; p < 0.01), was also associated with a modest increase in schizophrenia risk in a case-control sample (controls: N = 435; cases: N = 397, OR = 1.5). This same allele was then found associated with cognitive (n = 369) and fMRI (n = 313) intermediate phenotypes associated with schizophrenia-working memory (Cohen's d = 0.35), global cognition or g (d = 0.18), and prefrontal inefficiency (d = 0.36) as measured by BOLD fMRI during a working memory task. Together, these preclinical and clinical results suggest that variation in NKCC1 may increase risk for schizophrenia via alterations of mRNA expression at the molecular level and impairment of optimal prefrontal function at the macro or systems level.

DNA methylation signatures in development and aging of the human prefrontal cortex.

The human prefrontal cortex (PFC), a mastermind of the brain, is one of the last brain regions to mature. To investigate the role of epigenetics in the development of PFC, we examined DNA methylation in ∼14,500 genes at ∼27,000 CpG loci focused on 5' promoter regions in 108 subjects range in age from fetal to elderly. DNA methylation in the PFC shows unique temporal patterns across life. The fastest changes occur during the prenatal period, slow down markedly after birth and continue to slow further with aging. At the genome level, the transition from fetal to postnatal life is typified by a reversal of direction, from demethylation prenatally to increased methylation postnatally. DNA methylation is strongly associated with genotypic variants and correlates with expression of a subset of genes, including genes involved in brain development and in de novo DNA methylation. Our results indicate that promoter DNA methylation in the human PFC is a highly dynamic process modified by genetic variance and regulating gene transcription. Additional discovery is made possible with a stand-alone application, BrainCloudMethyl.

Transcript-specific associations of SLC12A5 (KCC2) in human prefrontal cortex with development, schizophrenia, and affective disorders.

The neuron-specific K(+)-Cl(-) cotransporter SLC12A5, also known as KCC2, helps mediate the electrophysiological effects of GABA. The pattern of KCC2 expression during early brain development suggests that its upregulation drives the postsynaptic switch of GABA from excitation to inhibition. We previously found decreased expression of full-length KCC2 in the postmortem hippocampus of patients with schizophrenia, but not in the dorsolateral prefrontal cortex (DLPFC). Using PCR and rapid amplification of cDNA ends, we discovered several previously unrecognized alternative KCC2 transcripts in both human adult and fetal brain in addition to the previously identified full-length (NM_020708.3) and truncated (AK098371) transcripts. We measured the expression levels of four relatively abundant truncated splice variants, including three novel transcripts (ΔEXON6, EXON2B, and EXON6B) and one previously described transcript (AK098371), in a large human cohort of nonpsychiatric controls across the lifespan, and in patients with schizophrenia and affective disorders. In SH-SY5Y cell lines, these transcripts were translated into proteins and expressed at their predicted sizes. Expression of the EXON6B transcript is increased in the DLPFC of patients with schizophrenia (p = 0.03) but decreased in patients with major depression (p = 0.04). The expression of AK098371 is associated with a GAD1 single nucleotide polymorphism (rs3749034) that previously has been associated with GAD67 expression and risk for schizophrenia. Our data confirm the developmental regulation of KCC2 expression, and provide evidence that KCC2 transcripts are differentially expressed in schizophrenia and affective disorders. Alternate transcripts from KCC2 may participate in the abnormal GABA signaling in the DLPFC associated with schizophrenia.

Expression of GABA signaling molecules KCC2, NKCC1, and GAD1 in cortical development and schizophrenia.

GABA signaling molecules are critical for both human brain development and the pathophysiology of schizophrenia. We examined the expression of transcripts derived from three genes related to GABA signaling [GAD1 (GAD67 and GAD25), SLC12A2 (NKCC1), and SLC12A5 (KCC2)] in the prefrontal cortex (PFC) and hippocampal formation of a large cohort of nonpsychiatric control human brains (n = 240) across the lifespan (from fetal week 14 to 80 years) and in patients with schizophrenia (n = 30-31), using quantitative RT-PCR. We also examined whether a schizophrenia risk-associated promoter SNP in GAD1 (rs3749034) is related to expression of these transcripts. Our studies revealed that development and maturation of both the PFC and hippocampal formation are characterized by progressive switches in expression from GAD25 to GAD67 and from NKCC1 to KCC2. Previous studies have demonstrated that the former leads to GABA synthesis, and the latter leads to switching from excitatory to inhibitory neurotransmission. In the hippocampal formation, GAD25/GAD67 and NKCC1/KCC2 ratios are increased in patients with schizophrenia, reflecting a potentially immature GABA physiology. Remarkably, GAD25/GAD67 and NKCC1/KCC2 expression ratios are associated with rs3749034 genotype, with risk alleles again predicting a relatively less mature pattern. These findings suggest that abnormalities in GABA signaling critical to brain development contribute to genetic risk for schizophrenia.

DISC1 splice variants are upregulated in schizophrenia and associated with risk polymorphisms.

Disrupted-In-Schizophrenia-1 (DISC1) is a promising susceptibility gene for major mental illness, but the mechanism of the clinical association is unknown. We searched for DISC1 transcripts in adult and fetal human brain and tested whether their expression is altered in patients with schizophrenia and is associated with genetic variation in DISC1. Many alternatively spliced transcripts were identified, including groups lacking exon 3 (Delta3), exons 7 and 8 (Delta7Delta8), an exon 3 insertion variant (extra short variant-1, Esv1), and intergenic splicing between TSNAX and DISC1. Isoforms Delta7Delta8, Esv1, and Delta3, which encode truncated DISC1 proteins, were expressed more abundantly during fetal development than during postnatal ages, and their expression was higher in the hippocampus of patients with schizophrenia. Schizophrenia risk-associated polymorphisms [non-synonymous SNPs rs821616 (Cys704Ser) and rs6675281 (Leu607Phe), and rs821597] were associated with the expression of Delta3 and Delta7Delta8. Moreover, the same allele at rs6675281, which predicted higher expression of these transcripts in the hippocampus, was associated with higher expression of DISC1Delta7Delta8 in lymphoblasts in an independent sample. Our results implicate a molecular mechanism of genetic risk associated with DISC1 involving specific alterations in gene processing.

Localization of cancer susceptibility genes by genome-wide single-nucleotide polymorphism linkage-disequilibrium mapping.

With the large numbers of single nucleotide polymorphisms (SNPs) available and new technologies that permit high throughput genotyping, we have investigated the possibility of the localization of disease genes with genome-wide panels of SNP markers and taking advantage of the linkage-disequilibrium (LD) between the disease gene and closely linked markers. For this purpose, we selected cases from the Ashkenazi Jewish population, in which the mutant alleles are expected to be identical by descent from a common founder and the regions of LD encompassing these mutant alleles are large. As a validation of this approach for localization, we performed two trials: one in autosomal recessive Bloom syndrome, in which a unique mutation of the BLM gene is present at elevated frequencies in cases, and the other in autosomal dominant hereditary nonpolyposis colorectal cancer (HNPCC), in which a unique mutation of MSH2 is present at elevated frequencies. In the Bloom syndrome trial, we genotyped 3,258 SNPs in 10 Jewish Bloom syndrome cases and 31 non-Bloom syndrome Jewish persons as a comparison group. In the HNPCC trial, we genotyped 8,549 SNPS in 13 Jewish HNPCC cases whose colon cancers exhibited microsatellite instability and in 63 healthy Jews as a comparison group. To identify significant associations, we performed (a) Fisher's exact test comparing genotypes at each locus in cases versus controls and (b) a haplotype analysis by estimating the frequency of haplotypes with the expectation-maximization algorithm and comparing haplotype frequencies in cases versus controls by logistic regression and a maximum likelihood ratio method. In the Bloom syndrome trial, by Fisher's exact test, statistically significant association was detected at a single locus, TSC0754862, which is a locus 1.7 million bp from BLM. Two-locus, three-locus, and four-locus haplotypes that included TSC0754862 and flanked BLM were also statistically more frequent in cases versus controls. In the HNPCC trial, although a significant P value was not obtained by the single SNP genotype analysis, significant associations were detected for several multilocus haplotypes in an 11-million-bp region that contained the MSH2 gene. This work demonstrates the power of the LD mapping approach in an isolated population and its general applicability to the identification of novel cancer-causing genes.

GAD2 Alternative Transcripts in the Human Prefrontal Cortex, and in Schizophrenia and Affective Disorders.

Genetic variation and early adverse environmental events work together to increase risk for schizophrenia. γ-aminobutyric acid (GABA), the major inhibitory neurotransmitter in adult mammalian brain, plays a major role in normal brain development, and has been strongly implicated in the pathobiology of schizophrenia. GABA synthesis is controlled by two glutamic acid decarboxylase (GAD) genes, GAD1 and GAD2, both of which produce a number of alternative transcripts. Genetic variants in the GAD1 gene are associated with increased risk for schizophrenia, and reduced expression of its major transcript in the human dorsolateral prefrontal cortex (DLPFC). No consistent changes in GAD2 expression have been found in brains from patients with schizophrenia. In this work, with the use of RNA sequencing and PCR technologies, we confirmed and tracked the expression of an alternative truncated transcript of GAD2 (ENST00000428517) in human control DLPFC homogenates across lifespan besides the well-known full length transcript of GAD2. In addition, using quantitative RT-PCR, expression of GAD2 full length and truncated transcripts were measured in the DLPFC of patients with schizophrenia, bipolar disorder and major depression. The expression of GAD2 full length transcript is decreased in the DLPFC of schizophrenia and bipolar disorder patients, while GAD2 truncated transcript is increased in bipolar disorder patients but decreased in schizophrenia patients. Moreover, the patients with schizophrenia with completed suicide or positive nicotine exposure showed significantly higher expression of GAD2 full length transcript. Alternative transcripts of GAD2 may be important in the growth and development of GABA-synthesizing neurons as well as abnormal GABA signaling in the DLPFC of patients with schizophrenia and affective disorders.

CHRNA7 and CHRFAM7A mRNAs: co-localized and their expression levels altered in the postmortem dorsolateral prefrontal cortex in major psychiatric disorders.

OBJECTIVE: CHRNA7, coding α-7 nicotinic acetylcholine receptor (α7 nAChR), is involved in cognition through interneuron modulation of dopamine and glutamate signaling. CHRNA7 and its partially duplicated chimeric gene CHRFAM7A have been implicated in schizophrenia through linkage and association studies. METHOD: Expression of CHRNA7 and CHRFAM7A mRNA was measured in the postmortem prefrontal cortex in more than 700 subjects, including patients with schizophrenia, bipolar disorder, major depression, and normal comparison subjects. The effects of antipsychotics and nicotine, as well as associations of CHRNA7 SNPs with gene expression, were explored. Fluorescent in-situ hybridization was used to examine coexpression of both transcripts in the human cortex. RESULTS: CHRFAM7A expression and CHRFAM7A/CHRNA7 ratios were higher in fetal compared with postnatal life, whereas CHRNA7 expression was relatively stable. CHRFAM7A expression was significantly elevated in all diagnostic groups, while CHRNA7 expression was reduced in the schizophrenia group and increased in the major depression group compared with the comparison group. CHRFAM7A/CHRNA7 ratios were significantly increased in the schizophrenia and bipolar disorder groups compared with the comparison group. There was no effect of nicotine or antipsychotics and no association of SNPs in CHRNA7 with expression. CHRNA7 and CHRFAM7A mRNAs were expressed in the same neuronal nuclei of the human neocortex. CONCLUSIONS: These data show preferential fetal CHRFAM7A expression in the human prefrontal cortex and suggest abnormalities in the CHRFAM7A/CHRNA7 ratios in schizophrenia and bipolar disorder, due mainly to overexpression of CHRFAM7A. Given that these transcripts are coexpressed in a subset of human cortical neurons and can interact to alter function of nAChRs, these results support the concept of aberrant function of nAChRs in mental illness.

DNA methylation changes in the postmortem dorsolateral prefrontal cortex of patients with schizophrenia.

BACKGROUND: Schizophrenia is a complex psychiatric disorder with a lifetime morbidity rate of 0.5-1.0%. The pathophysiology of schizophrenia still remains obscure. Accumulating evidence indicates that DNA methylation, which is the addition of a methyl group to the cytosine in a CpG dinucleotide, might play an important role in the pathogenesis of schizophrenia. METHODS: To gain further insight into the molecular mechanisms underlying schizophrenia, a genome-wide DNA methylation profiling (27,578 CpG dinucleotides spanning 14,495 genes) of the human dorsolateral prefrontal cortex (DLPFC) was conducted in a large cohort (n = 216) of well characterized specimens from individuals with schizophrenia and non-psychiatric controls, combined with an analysis of genetic variance at ~880,000 SNPs. RESULTS: Aberrant DNA methylation in schizophrenia was identified at 107 CpG sites at 5% Bonferroni correction (p < 1.99 × 10(-6)). Of these significantly altered sites, hyper-DNA methylation was observed at 79 sites (73.8%), mostly in the CpG islands (CGIs) and in the regions flanking CGIs (CGI: 31 sites; CGI shore: 35 sites; CGI shelf: 3 sites). Furthermore, a large number of cis-methylation quantitative trait loci (mQTL) were identified, including associations with risk SNPs implicated in schizophrenia. CONCLUSIONS: These results suggest that altered DNA methylation might be involved in the pathophysiology and/or treatment of schizophrenia, and that a combination of epigenetic and genetic approaches will be useful to understanding the molecular mechanism of this complex disorder.

Analysis of copy number variations in brain DNA from patients with schizophrenia and other psychiatric disorders.

BACKGROUND: Clinical studies have identified several regions of the genome with copy number variations (CNVs) associated with diverse neurodevelopmental behavioral disorders. METHODS: We analyzed 1 million (M) single nucleotide polymorphism genotype arrays for evidence of previously reported recurrent CNVs and enriched genome-wide CNV burden in DNA from 600 brains, including 441 individuals with various psychiatric diagnoses. We explored gene expression in the dorsolateral prefrontal cortex in selected cases with CNVs and in other subjects with Illumina BeadArrays (568 subjects in total) and additionally in 66-92 subjects with quantitative real-time polymerase chain reaction. RESULTS: The CNVs in previously reported genomic regions were identified in 4 of 193 patients with the diagnosis of schizophrenia (1q21.1, 11q25, 15q11.2, 22q11), 4 of 238 patients with mood disorders (11q25, 15q11.2, 22q11), and 1 of 10 patients with autism (2p16.3). No evidence of increased genome-wide CNV burden was observed in cases with schizophrenia or mood disorders, although the study is underpowered to observe rare events. Messenger RNA expression patterns suggested incomplete molecular penetrance of observed CNVs. CONCLUSIONS: Our data confirm in brain DNA the presence of certain recurrent CNVs in a small percentage of patients with psychiatric diagnoses.

Evidence of sex-modulated association of ZNF804A with schizophrenia.

BACKGROUND: The single nucleotide polymorphism (SNP) rs1344706 in ZNF804A (2q32.1) has been associated with schizophrenia in a genome-wide association study (GWAS). A recent candidate gene study, which replicated the positive association with rs1344706, identified another positive SNP (rs7597593) in ZNF804A associated with schizophrenia. METHODS: We performed an association study of rs7597593 in four GWAS cohorts of European ancestry. Postmortem human brain expression data of normal Caucasian individuals (n = 89) was also analyzed for examining the effect of rs7597593 on ZNF804A messenger RNA expression, using logistic regression and linear regression. RESULTS: We found that rs7597593 was significantly associated with schizophrenia in the combined GWAS datasets (n = 5023, odds ratio [OR](combined) = 1.15, p = .0011). Analysis of stratification by sex showed that the association was driven by the female subjects (OR = 1.29, p = .0002) and was not significant in male subjects (OR = 1.08, p = .148) in the combined sample of four cohorts. A sex by genotype interaction was near significant in both the Genetic Association Information Network sample (p = .0532) and the combined sample of four cohorts (p(combined) = .0531). Gene expression analysis showed no main effects but a significant female-specific association (p(female) = .047, p(male) = .335) and sex by genotype interaction (p = .0166) for rs7597593. CONCLUSIONS: Our data suggest a clinical and molecular modulation by sex of the association of ZNF804A SNP rs7597593 and risk of schizophrenia.

Characterization of STIP, a multi-domain nuclear protein, highly conserved in metazoans, and essential for embryogenesis in Caenorhabditis elegans.

We report here the identification and characterization of STIP, a multi-domain nuclear protein that contains a G-patch, a coiled-coil, and several short tryptophan-tryptophan repeats highly conserved in metazoan species. To analyze their functional role in vivo, we cloned nematode stip-1 genes and determined the spatiotemporal pattern of Caenorhabditis elegans STIP-1 protein. RNA analyses and Western blots revealed that stip-1 mRNA was produced via trans-splicing and translated as a 95-kDa protein. Using reporter constructs, we found STIP-1 to be expressed at all developmental stages and in many tissue/cell types including worm oocyte nuclei. We found that STIP-1 is targeted to the nucleus and forms large polymers with a rod-like shape when expressed in mammalian cells. Using deletion mutants, we mapped the regions of STIP-1 involved in nuclear import and polymer assembly. We further showed that knockdown of C. elegans stip-1 by RNA interference arrested development and resulted in morphologic abnormalities around the 16-cell stage followed by 100% lethality, suggesting its essential role in worm embryogenesis. Importantly, the embryonic lethal phenotype could be faithfully rescued with Drosophila and human genes via transgenic expression. Our data provide the first direct evidence that STIP have a conserved essential nuclear function across metazoans from worms to humans.

Localization of breast cancer susceptibility loci by genome-wide SNP linkage disequilibrium mapping.

We studied the feasibility of a novel approach to localize breast cancer susceptibility genes, using a low-density genome-wide panel of single-nucleotide polymorphisms and taking advantage of large regions of linkage disequilibrium (LD) flanking Jewish disease genes in high-risk cases. With Affymetrix GeneChip arrays, we genotyped 8,576 polymorphisms in three sets of Ashkenazi Jewish breast cancer cases: a "validation" set of 27 breast cancer cases, all of whom carried the BRCA2*6174delT founder mutation; a "field" set of 19 breast cancer cases from male breast cancer kindreds, which simulated conditions for finding new genes; and a "test" set of 57 probands from breast cancer kindreds (4 or more cases/kindred), in which mutations in BRCA1 and BRCA2 had been excluded. To identify associations, we compared the frequency of genotypes and haplotypes in cases vs. controls by the Fisher's exact test and a maximum likelihood ratio test. In the "validation" set, we demonstrated the presence of a region of linkage disequilibrium on BRCA2*6174delT chromosomes that spanned over 5 million bases. In the "field" set, we showed that this large region of linkage disequilibrium flanking BRCA2 was detectable despite the presence of heterogeneity in the sample set. Finally, in the "test" set, at least three regions of interest emerged that could contain novel breast cancer genes, one of which had been identified previously by linkage analysis. While these results demonstrate the feasibility of genome-wide association strategies, further application of this approach will critically depend on optimizing the density and distribution of SNPs and the size and type of study design.

Intra-nuclear trafficking of the BLM helicase to DNA damage-induced foci is regulated by SUMO modification.

The Bloom syndrome gene, BLM, encodes a RecQ DNA helicase that when absent from the cell results in genomic instability and cancer predisposition. We show here that BLM is a substrate for small ubiquitin-like modifier (SUMO) modification, with lysines at K317, K331, K334 and K347 being preferred sites of modification. Unlike normal BLM, a double mutant BLM protein with lysine to arginine substitutions at residues 317 and 331 was not modified by SUMO, and it failed to localize efficiently to the PML nuclear bodies. Rather, double mutant BLM protein induced the formation of DNA damage-induced foci (DDI) that contained BRCA1 protein and phosphorylated histone H2AX. Double mutant BLM only partially complemented the genomic instability phenotypes of Bloom syndrome cells as assessed by sister-chromatid exchange and micronuclei formation assays. These results constitute evidence that BLM is a DNA damage sensor that signals the formation of DDI, and they establish SUMO modification as a negative regulator of BLM's signaling function.