Analysis of somatic mutations in 131 human brains reveals aging-associated hypermutability.

We analyzed 131 human brains (44 neurotypical, 19 with Tourette syndrome, 9 with schizophrenia, and 59 with autism) for somatic mutations after whole genome sequencing to a depth of more than 200×. Typically, brains had 20 to 60 detectable single-nucleotide mutations, but ~6% of brains harbored hundreds of somatic mutations. Hypermutability was associated with age and damaging mutations in genes implicated in cancers and, in some brains, reflected in vivo clonal expansions. Somatic duplications, likely arising during development, were found in ~5% of normal and diseased brains, reflecting background mutagenesis. Brains with autism were associated with mutations creating putative transcription factor binding motifs in enhancer-like regions in the developing brain. The top-ranked affected motifs corresponded to MEIS (myeloid ectopic viral integration site) transcription factors, suggesting a potential link between their involvement in gene regulation and autism.

Comprehensive identification of somatic nucleotide variants in human brain tissue.

BACKGROUND: Post-zygotic mutations incurred during DNA replication, DNA repair, and other cellular processes lead to somatic mosaicism. Somatic mosaicism is an established cause of various diseases, including cancers. However, detecting mosaic variants in DNA from non-cancerous somatic tissues poses significant challenges, particularly if the variants only are present in a small fraction of cells. RESULTS: Here, the Brain Somatic Mosaicism Network conducts a coordinated, multi-institutional study to examine the ability of existing methods to detect simulated somatic single-nucleotide variants (SNVs) in DNA mixing experiments, generate multiple replicates of whole-genome sequencing data from the dorsolateral prefrontal cortex, other brain regions, dura mater, and dural fibroblasts of a single neurotypical individual, devise strategies to discover somatic SNVs, and apply various approaches to validate somatic SNVs. These efforts lead to the identification of 43 bona fide somatic SNVs that range in variant allele fractions from ~ 0.005 to ~ 0.28. Guided by these results, we devise best practices for calling mosaic SNVs from 250× whole-genome sequencing data in the accessible portion of the human genome that achieve 90% specificity and sensitivity. Finally, we demonstrate that analysis of multiple bulk DNA samples from a single individual allows the reconstruction of early developmental cell lineage trees. CONCLUSIONS: This study provides a unified set of best practices to detect somatic SNVs in non-cancerous tissues. The data and methods are freely available to the scientific community and should serve as a guide to assess the contributions of somatic SNVs to neuropsychiatric diseases.

Genome-wide analyses of smoking behaviors in schizophrenia: Findings from the Psychiatric Genomics Consortium.

While 17% of US adults use tobacco regularly, smoking rates among persons with schizophrenia are upwards of 60%. Research supports a shared etiological basis for smoking and schizophrenia, including findings from genome-wide association studies (GWAS). However, few studies have directly tested whether the same or distinct genetic variants also influence smoking behavior among schizophrenia cases. Using data from the Psychiatric Genomics Consortium (PGC) study of schizophrenia (35476 cases, 46839 controls), we estimated genetic correlations between these traits and tested whether polygenic risk scores (PRS) constructed from the results of smoking behaviors GWAS were associated with schizophrenia risk or smoking behaviors among schizophrenia cases. Results indicated significant genetic correlations of schizophrenia with smoking initiation (rg = 0.159; P = 5.05 × 10-10), cigarettes-smoked-per-day (rg = 0.094; P = 0.006), and age-of-onset of smoking (rg = 0.10; P = 0.009). Comparing smoking behaviors among schizophrenia cases to the general population, we observe positive genetic correlations for smoking initiation (rg = 0.624, P = 0.002) and cigarettes-smoked-per-day (rg = 0.689, P = 0.120). Similarly, TAG-based PRS for smoking initiation and cigarettes-smoked-per-day were significantly associated with smoking initiation (P = 3.49 × 10-5) and cigarettes-smoked-per-day (P = 0.007) among schizophrenia cases. We performed the first GWAS of smoking behavior among schizophrenia cases and identified a novel association with cigarettes-smoked-per-day upstream of the TMEM106B gene on chromosome 7p21.3 (rs148253479, P = 3.18 × 10-8, n = 3520). Results provide evidence of a partially shared genetic basis for schizophrenia and smoking behaviors. Additionally, genetic risk factors for smoking behaviors were largely shared across schizophrenia and non-schizophrenia populations. Future research should address mechanisms underlying these associations to aid both schizophrenia and smoking treatment and prevention efforts.

Schizophrenia candidate gene ERBB4: covert routes of vulnerability to psychosis detected at the population level.

Prior genetic and functional evidence established ERBB4 as a probable schizophrenia susceptibility gene that may confer risk via modulating brain information processing dependent on the integrity of frontotemporal brain circuitry. Utilizing retrospective data drawn from the cross-sectional population-based Athens Study of Psychosis Proneness and Incidence of Schizophrenia (ASPIS) (n = 1127), we attempted to independently replicate and further extend previous findings by examining the effects of ERBB4 gene variants on 3 broad population-based psychosis-related phenotypes: verbal working memory (VWM), trait schizotypy, and stress-induced subclinical psychotic experiences (PE). Three common ERBB4 single nucleotide polymorphisms that were previously associated with schizophrenia and impaired frontotemporal-related information processing (rs7598440, rs839523, and rs707284), their haplotypes, and corresponding diplotypes were tested. VWM performance was significantly associated with rs839523 and rs707284 markers even after correction for multiple testing, thus validating reported findings that have implicated ERBB4 gene variation on working memory. No associations were detected between these ERBB4 variants and trait schizotypy. However, we were able to detect a significant effect of rs7598440 marker on PE expressed under stressful environmental conditions. Combined haplotype analysis of the above 3 markers, identified a "yin-yang" pattern of association, confirmed at the diplotype level. While GGG haplotype homozygotes were associated with "protective" effects on VWM performance and PE, AAA "risk" haplotype carriers were associated with worse VWM performance and simultaneously exhibited significantly elevated PE. This dual, possibly pleiotropic, impact on frontotemporal circuitry and increased sensitivity to psychosocial stress may represent subtle manifestations of ERBB4-related vulnerability to psychosis, expressed at the population level.

Neuregulin 1-ErbB4-PI3K signaling in schizophrenia and phosphoinositide 3-kinase-p110δ inhibition as a potential therapeutic strategy.

Neuregulin 1 (NRG1) and ErbB4, critical neurodevelopmental genes, are implicated in schizophrenia, but the mediating mechanisms are unknown. Here we identify a genetically regulated, pharmacologically targetable, risk pathway associated with schizophrenia and with ErbB4 genetic variation involving increased expression of a PI3K-linked ErbB4 receptor (CYT-1) and the phosphoinositide 3-kinase subunit, p110δ (PIK3CD). In human lymphoblasts, NRG1-mediated phosphatidyl-inositol,3,4,5 triphosphate [PI(3,4,5)P3] signaling is predicted by schizophrenia-associated ErbB4 genotype and PIK3CD levels and is impaired in patients with schizophrenia. In human brain, the same ErbB4 genotype again predicts increased PIK3CD expression. Pharmacological inhibition of p110δ using the small molecule inhibitor, IC87114, blocks the effects of amphetamine in a mouse pharmacological model of psychosis and reverses schizophrenia-related phenotypes in a rat neonatal ventral hippocampal lesion model. Consistent with these antipsychotic-like properties, IC87114 increases AKT phosphorylation in brains of treated mice, implicating a mechanism of action. Finally, in two family-based genetic studies, PIK3CD shows evidence of association with schizophrenia. Our data provide insight into a mechanism of ErbB4 association with schizophrenia; reveal a previously unidentified biological and disease link between NRG1-ErbB4, p110δ, and AKT; and suggest that p110δ is a previously undescribed therapeutic target for the treatment of psychiatric disorders.

Biological validation of increased schizophrenia risk with NRG1, ERBB4, and AKT1 epistasis via functional neuroimaging in healthy controls.

CONTEXT: NRG1 is a schizophrenia candidate gene and plays an important role in brain development and neural function. Schizophrenia is a complex disorder, with etiology likely due to epistasis. OBJECTIVE: To examine epistasis between NRG1 and selected N-methyl-d-aspartate-glutamate pathway partners implicated in its effects, including ERBB4, AKT1, DLG4, NOS1, and NOS1AP. DESIGN: Schizophrenia case-control sample analyzed using machine learning algorithms and logistic regression with follow-up using neuroimaging on an independent sample of healthy controls. PARTICIPANTS: A referred sample of schizophrenic patients (n = 296) meeting DSM-IV criteria for schizophrenia spectrum disorder and a volunteer sample of controls for case-control comparison (n = 365) and a separate volunteer sample of controls for neuroimaging (n = 172). MAIN OUTCOME MEASURES: Epistatic association between single-nucleotide polymorphisms (SNPs) and case-control status; epistatic association between SNPs and the blood oxygen level-dependent physiological response during working memory measured by functional magnetic resonance imaging. RESULTS: We observed interaction between NRG1 5' and 3' SNPs rs4560751 and rs3802160 (likelihood ratio test P = .00020) and schizophrenia, which was validated using functional magnetic resonance imaging of working memory in healthy controls; carriers of risk-associated genotypes showed inefficient processing in the dorsolateral prefrontal cortex (P = .015, familywise error corrected). We observed epistasis between NRG1 (rs10503929; Thr286/289/294Met) and its receptor ERBB4 (rs1026882; likelihood ratio test P = .035); a 3-way interaction with these 2 SNPs and AKT1 (rs2494734) was also observed (odds ratio, 27.13; 95% confidence interval, 3.30-223.03; likelihood ratio test P = .042). These same 2- and 3-way interactions were further biologically validated via functional magnetic resonance imaging: healthy individuals carrying risk genotypes for NRG1 and ERBB4, or these 2 together with AKT1, were disproportionately less efficient in dorsolateral prefrontal cortex processing. Lower-level interactions were not observed between NRG1 /ERBB4 and AKT1 in association or neuroimaging, consistent with biological evidence that NRG1 × ERBB4 interaction modulates downstream AKT1 signaling. CONCLUSION: Our data suggest complex epistatic effects implicating an NRG1 molecular pathway in cognitive brain function and the pathogenesis of schizophrenia.

Genetic variation in AKT1 is linked to dopamine-associated prefrontal cortical structure and function in humans.

AKT1-dependent molecular pathways control diverse aspects of cellular development and adaptation, including interactions with neuronal dopaminergic signaling. If AKT1 has an impact on dopaminergic signaling, then genetic variation in AKT1 would be associated with brain phenotypes related to cortical dopaminergic function. Here, we provide evidence that a coding variation in AKT1 that affects protein expression in human B lymphoblasts influenced several brain measures related to dopaminergic function. Cognitive performance linked to frontostriatal circuitry, prefrontal physiology during executive function, and frontostriatal gray-matter volume on MRI were altered in subjects with the AKT1 variation. Moreover, on neuroimaging measures with a main effect of the AKT1 genotype, there was significant epistasis with a functional polymorphism (Val158Met) in catechol-O-methyltransferase [COMT], a gene that indexes cortical synaptic dopamine. This genetic interaction was consistent with the putative role of AKT1 in dopaminergic signaling. Supportive of an earlier tentative association of AKT1 with schizophrenia, we also found that this AKT1 variant was associated with risk for schizophrenia. These data implicate AKT1 in modulating human prefrontal-striatal structure and function and suggest that the mechanism of this effect may be coupled to dopaminergic signaling and relevant to the expression of psychosis.

Variants in the estrogen receptor alpha gene and its mRNA contribute to risk for schizophrenia.

Estrogen modifies human emotion and cognition and impacts symptoms of schizophrenia. We hypothesized that the variation in the estrogen receptor alpha (ESR1) gene and cortical ESR1 mRNA is associated with schizophrenia. In a small case-control genetic association analysis of postmortem brain tissue, genotype CC (rs2234693) and haplotypes containing the C allele of a single-nucleotide polymorphism (SNP) in intron1 (PvuII) were more frequent in African American schizophrenics (P = 0.01-0.001). In a follow-up family-based association analysis, we found overtransmission of PvuII allele C and a PvuII C-containing haplotype (P = 0.01-0.03) to African American and Caucasian patients with schizophrenia. Schizophrenics with the 'at risk' PvuII genotype had lower ESR1 mRNA levels in the frontal cortex. Eighteen ESR1 splice variants and decreased frequencies of the wild-type ESR1 mRNA were detected in schizophrenia. In one patient, a unique ESR1 transcript with a genomic insert encoding a premature stop codon and a truncated ESR1 protein lacking most of the estrogen binding domain was the only transcript detected. Using a luciferase assay, we found that mRNA encoding a truncated ESR1 significantly attenuates gene expression at estrogen-response elements demonstrating a dominant negative function. An intron 6 SNP [rs2273207(G)] was associated with an ESR1 splice variant missing exon seven. The T allele of another intron 6 SNP was part of a 3' haplotype less common in schizophrenia [rs2273206(T), rs2273207(G), rs2228480(G)]. Thus, the variation in the ESR1 gene is associated with schizophrenia and the mechanism of this association may involve alternative gene regulation and transcript processing.

Evidence that the BLOC-1 protein dysbindin modulates dopamine D2 receptor internalization and signaling but not D1 internalization.

The schizophrenia susceptibility gene dystrobrevin-binding protein 1 (DTNBP1) encodes dysbindin, which along with its binding partner Muted is an essential component of the biogenesis of lysosome-related organelles complex 1 (BLOC-1). Dysbindin expression is reduced in schizophrenic brain tissue, but the molecular mechanisms by which this contributes to pathogenesis and symptomatology are unknown. We studied the effects of transfection of DTNBP1 siRNA on cell surface levels of dopamine D2 receptor (DRD2) in human SH-SY5Y neuroblastoma cells and in rat primary cortical neurons. DTNBP1 siRNA decreased dysbindin protein, increased cell surface DRD2 and blocked dopamine-induced DRD2 internalization. MUTED siRNA produced similar effects. In contrast, decreased dysbindin did not change dopamine D1 receptor (DRD1) levels, or its basal or dopamine-induced internalization. The DRD2 agonist quinpirole reduced phosphorylation of CREB (cAMP response element-binding protein) in dysbindin downregulated cells, demonstrating enhanced intracellular signaling caused by the upregulation of DRD2. This is the first demonstration of a schizophrenia susceptibility gene exerting a functional effect on DRD2 signaling, a pathway that has long been implicated in the illness. We propose a molecular mechanism for pathogenesis in which risk alleles in DTNBP1, or other factors that also downregulate dysbindin, compromise the ability of BLOC-1 to traffic DRD2 toward degradation, but has little effect on DRD1 trafficking. Impaired trafficking of DRD2 decreases dopamine-induced internalization, and with more receptors retained on the cell surface, dopamine stimulation produces excess intracellular signaling. Such an increase in DRD2 signaling relative to DRD1 would contribute to the imbalances in dopaminergic neurotransmission characteristic of schizophrenia.

Alpha7 nicotinic acetylcholine receptor mRNA expression and binding in postmortem human brain are associated with genetic variation in neuregulin 1.

Studies in cell culture and in animals suggest that neuregulin 1 (NRG1), a probable schizophrenia susceptibility gene, regulates the expression of the alpha7 nicotinic acetylcholine receptors (nAChRs). We hypothesized that schizophrenia-associated allelic variations within the NRG1 gene, via their effects on NRG1 isoform expression, would be associated with alterations in nAChR alpha7 receptor levels. We examined the effects of four disease-associated single-nucleotide polymorphisms (SNPs) in the 5' region of the NRG1 gene on nAChR alpha7 mRNA transcript expression in both the dorsolateral prefrontal cortex (DLPFC) and hippocampus of normal controls and patients with schizophrenia using quantitative real-time PCR. NRG1 risk alleles at SNPs SNP8NRG221132 and rs6994992 predicted significantly lower nAChR alpha7 mRNA expression in the DLPFC. Haplotypes containing the risk alleles at the above SNPs were also associated with lower expression of nAChR alpha7 in the DLPFC. The genotype effect for rs6994992 and the haplotype effect were more pronounced within the schizophrenic patient group. To determine whether receptor levels follow that of mRNA expression, we performed receptor binding and autoradiography using [(125)I] alpha-bungarotoxin in the DLPFC. Consistent with the mRNA findings, we found a decrease in binding in risk allele carriers of SNP8NRG221132 as compared with heterozygous individuals. Together, these results suggest that the molecular mechanism of the association between NRG1 risk alleles and schizophrenia may include down-regulation of nAChR alpha7 expression.

Risk for premenstrual dysphoric disorder is associated with genetic variation in ESR1, the estrogen receptor alpha gene.

BACKGROUND: Premenstrual dysphoric disorder (PMDD) is a heritable mood disorder that is triggered by gonadal steroids during the luteal phase in susceptible women. METHODS: We performed haplotype analyses of estrogen receptors alpha and beta (ESR1 and ESR2) in 91 women with prospectively confirmed PMDD and 56 control subjects to investigate possible sources of the genetic susceptibility to affective dysregulation induced by normal levels of gonadal steroids. We also examined associations with the valine (Val)158methionine (Met) single nucleotide polymorphism (SNP) of the gene for catechol-O-methyltransferase (COMT), an enzyme involved in estrogen metabolism and prefrontal cortical activation. RESULTS: Four SNPs in intron 4 of ESR1 showed significantly different genotype and allele distributions between patients and control subjects. Significant case-control differences were seen in sliding-window analyses of two-, three-, and four-marker haplotypes but only in those haplotypes containing SNPs in intron 4 that were positive in the single-locus analysis. No significant associations were observed with ESR2 or with the COMT Val158Met polymorphism, although the significant associations with ESR1 were observed only in those with the Val/Val genotype. CONCLUSIONS: These are the first positive (albeit preliminary) genetic findings in this reproductive endocrine-related mood disorder and involve the receptor for a hormone that is pathogenically relevant.

Genetic evidence implicating DARPP-32 in human frontostriatal structure, function, and cognition.

Dopamine- and cAMP-regulated phosphoprotein of molecular weight 32 kDa (DARPP-32), encoded by PPP1R1B, is a pivotal integrator of information in dopaminoceptive neurons, regulating the response to neuroleptics, psychotomimetics, and drugs of abuse, and affecting striatal function and plasticity. Despite extensive preclinical work, there are almost no data on DARPP-32 function in humans. Here, we identify, through resequencing in 298 chromosomes, a frequent PPP1R1B haplotype predicting mRNA expression of PPP1R1B isoforms in postmortem human brain. This haplotype was associated with enhanced performance on several cognitive tests that depend on frontostriatal function. Multimodal imaging of healthy subjects revealed an impact of the haplotype on neostriatal volume, activation, and the functional connectivity of the prefrontal cortex. The haplotype was associated with the risk for schizophrenia in 1 family-based association analysis. Our convergent results identify a prefrontal-neostriatal system affected by variation in PPP1R1B and suggest that DARPP-32 plays a pivotal role in cognitive function and possibly in the pathogenesis of schizophrenia.

Expression of DISC1 binding partners is reduced in schizophrenia and associated with DISC1 SNPs.

DISC1 has been identified as a schizophrenia susceptibility gene based on linkage and SNP association studies and clinical data suggesting that risk SNPs impact on hippocampal structure and function. In cell and animal models, C-terminus-truncated DISC1 disrupts intracellular transport, neural architecture and migration, perhaps because it fails to interact with binding partners involved in neuronal differentiation such as fasciculation and elongation protein zeta-1 (FEZ1), platelet-activating factor acetylhydrolase, isoform Ib, PAFAH1B1 or lissencephaly 1 protein (LIS1) and nuclear distribution element-like (NUDEL). We hypothesized that altered expression of DISC1 and/or its molecular partners may underlie its pathogenic role in schizophrenia and explain its genetic association. We examined the expression of DISC1 and these selected binding partners as well as reelin, a protein in a related signaling pathway, in the hippocampus and dorsolateral prefrontal cortex of postmortem human brain patients with schizophrenia and controls. We found no difference in the expression of DISC1 or reelin mRNA in schizophrenia and no association with previously identified risk DISC1 SNPs. However, the expression of NUDEL, FEZ1 and LIS1 was each significantly reduced in the brain tissue from patients with schizophrenia and expression of each showed association with high-risk DISC1 polymorphisms. Although, many other DISC1 binding partners still need to be investigated, these data implicate genetically linked abnormalities in the DISC1 molecular pathway in the pathophysiology of schizophrenia.

Evidence of novel neuronal functions of dysbindin, a susceptibility gene for schizophrenia.

Genetic variation in dysbindin (DTNBP1: dystrobrevin-binding protein 1) has recently been shown to be associated with schizophrenia. The dysbindin gene is located at chromosome 6p22.3, one of the most promising susceptibility loci in schizophrenia linkage studies. We attempted to replicate this association in a Japanese sample of 670 patients with schizophrenia and 588 controls. We found a nominally significant association with schizophrenia for four single nucleotide polymorphisms and stronger evidence for association in a multi-marker haplotype analysis (P = 0.00028). We then explored functions of dysbindin protein in primary cortical neuronal culture. Overexpression of dysbindin induced the expression of two pre-synaptic proteins, SNAP25 and synapsin I, and increased extracellular basal glutamate levels and release of glutamate evoked by high potassium. Conversely, knockdown of endogenous dysbindin protein by small interfering RNA (siRNA) resulted in the reduction of pre-synaptic protein expression and glutamate release, suggesting that dysbindin might influence exocytotic glutamate release via upregulation of the molecules in pre-synaptic machinery. The overexpression of dysbindin increased phosphorylation of Akt protein and protected cortical neurons against neuronal death due to serum deprivation and these effects were blocked by LY294002, a phosphatidylinositol 3-kinase (PI3-kinase) inhibitor. SiRNA-mediated silencing of dysbindin protein diminished Akt phosphorylation and facilitated neuronal death induced by serum deprivation, suggesting that dysbindin promotes neuronal viability through PI3-kinase-Akt signaling. Genetic variants associated with impairments of these functions of dysbindin could play an important role in the pathogenesis of schizophrenia.

Candidate genes for nicotine dependence via linkage, epistasis, and bioinformatics.

Many smoking-related phenotypes are substantially heritable. One genome scan of nicotine dependence (ND) has been published and several others are in progress and should be completed in the next 5 years. The goal of this hypothesis-generating study was two-fold. First, we present further analyses of our genome scan data for ND published by Straub et al. [1999: Mol Psychiatry 4:129-144] (PMID: 10208445). Second, we used the method described by Cox et al. [1999: Nat Genet 21:213-215] (PMID: 9988276) to search for epistatic loci across the markers used in the genome scan. The overall results of the genome scan nearly reached the rigorous Lander and Kruglyak [1995: Nat Genet 11:241-247] criteria for "significant" linkage with the best findings on chromosomes 10 and 2. We then looked for correspondence between genes located in the 10 regions implicated in affected sibling pair (ASP) and epistatic linkage analyses with a list of genes suggested by microarray studies of experimental nicotine exposure and candidate genes from the literature. We found correspondence between linkage and microarray/candidate gene studies for genes involved with the mitogen-activated protein kinase (MAPK) signaling system, nuclear factor kappa B (NFKB) complex, neuropeptide Y (NPY) neurotransmission, a nicotinic receptor subunit (CHRNA2), the vesicular monoamine transporter (SLC18A2), genes in pathways implicated in human anxiety (HTR7, TDO2, and the endozepine-related protein precursor, DKFZP434A2417), and the micro 1-opioid receptor (OPRM1). Although the hypotheses resulting from these linkage and bioinformatic analyses are plausible and intriguing, their ultimate worth depends on replication in additional linkage samples and in future experimental studies.