Divergent landscapes of A-to-I editing in postmortem and living human brain.

Adenosine-to-inosine (A-to-I) editing is a prevalent post-transcriptional RNA modification within the brain. Yet, most research has relied on postmortem samples, assuming it is an accurate representation of RNA biology in the living brain. We challenge this assumption by comparing A-to-I editing between postmortem and living prefrontal cortical tissues. Major differences were found, with over 70,000 A-to-I sites showing higher editing levels in postmortem tissues. Increased A-to-I editing in postmortem tissues is linked to higher ADAR and ADARB1 expression, is more pronounced in non-neuronal cells, and indicative of postmortem activation of inflammation and hypoxia. Higher A-to-I editing in living tissues marks sites that are evolutionarily preserved, synaptic, developmentally timed, and disrupted in neurological conditions. Common genetic variants were also found to differentially affect A-to-I editing levels in living versus postmortem tissues. Collectively, these discoveries offer more nuanced and accurate insights into the regulatory mechanisms of RNA editing in the human brain.

Divergent landscapes of A-to-I editing in postmortem and living human brain.

Adenosine-to-inosine (A-to-I) editing is a prevalent post-transcriptional RNA modification within the brain. Yet, most research has relied on postmortem samples, assuming it is an accurate representation of RNA biology in the living brain. We challenge this assumption by comparing A-to-I editing between postmortem and living prefrontal cortical tissues. Major differences were found, with over 70,000 A-to-I sites showing higher editing levels in postmortem tissues. Increased A-to-I editing in postmortem tissues is linked to higher ADAR1 and ADARB1 expression, is more pronounced in non-neuronal cells, and indicative of postmortem activation of inflammation and hypoxia. Higher A-to-I editing in living tissues marks sites that are evolutionarily preserved, synaptic, developmentally timed, and disrupted in neurological conditions. Common genetic variants were also found to differentially affect A-to-I editing levels in living versus postmortem tissues. Collectively, these discoveries illuminate the nuanced functions and intricate regulatory mechanisms of RNA editing within the human brain.

Analysis of Genetically Regulated Gene Expression Identifies a Prefrontal PTSD Gene, SNRNP35, Specific to Military Cohorts.

To reveal post-traumatic stress disorder (PTSD) genetic risk influences on tissue-specific gene expression, we use brain and non-brain transcriptomic imputation. We impute genetically regulated gene expression (GReX) in 29,539 PTSD cases and 166,145 controls from 70 ancestry-specific cohorts and identify 18 significant GReX-PTSD associations corresponding to specific tissue-gene pairs. The results suggest substantial genetic heterogeneity based on ancestry, cohort type (military versus civilian), and sex. Two study-wide significant PTSD associations are identified in European and military European cohorts; ZNF140 is predicted to be upregulated in whole blood, and SNRNP35 is predicted to be downregulated in dorsolateral prefrontal cortex, respectively. In peripheral leukocytes from 175 marines, the observed PTSD differential gene expression correlates with the predicted differences for these individuals, and deployment stress produces glucocorticoid-regulated expression changes that include downregulation of both ZNF140 and SNRNP35. SNRNP35 knockdown in cells validates its functional role in U12-intron splicing. Finally, exogenous glucocorticoids in mice downregulate prefrontal Snrnp35 expression.

A molecular approach to treating cognition in schizophrenia by calcium channel blockade: An open-label pilot study of the calcium-channel antagonist isradipine.

Cognitive impairment is a prominent and difficult to treat symptom in schizophrenia (SZ), which is directly related to functional disability. A variant in the gene coding for the alpha 1C subunit of L-type voltage gated calcium channel (CACNA1C) has been shown to negatively affect several neurocognitive domains. We conducted a 4-week, open label, pilot study of isradipine, a calcium channel blocker, to determine its feasibility, safety, and efficacy in improving cognition in SZ patients. Ten adults with stable SZ were started on a flexible dose of isradipine 5 mg/day (up to 10 mg/day) for 4 weeks. Weekly in-person visits tracked side effects and symptoms while neurocognition and functional capacity were assessed at baseline and week 4. There were no serious adverse events reported. Newly emergent side effects were dizziness (1 new incidence at week 4); difficulty sleeping (2 new incidences at week 4); and decreased energy (3 new incidences at week 4). 1 patient discontinued medication and was withdrawn. Treatment did not exacerbate clinical symptoms. Although power is limited, results indicate no clear benefit on neurocognition but a positive effect (baseline mean = 6.8 ± 1.3 to week 4 mean = 7.9 ± 1.1; t = 2.91, p = 0.017) on functional capacity was noted. This open label, pilot study provides preliminary evidence that isradipine is a relatively safe medication when used adjunctively in SZ patients. This study suggests that isradipine offers no clear cognitive and only minimal functional benefit; however, additional studies may be warranted in symptomatic patients, or those with specific CACNA1C genotypes.

Characterization of Single Gene Copy Number Variants in Schizophrenia.

BACKGROUND: Genetic studies of schizophrenia have implicated numerous risk loci including several copy number variants (CNVs) of large effect and hundreds of loci of small effect. In only a few cases has a specific gene been clearly identified. Rare CNVs affecting a single gene offer a potential avenue to discovering schizophrenia risk genes. METHODS: CNVs were generated from exome sequencing of 4913 schizophrenia cases and 6188 control subjects from Sweden. We integrated two CNV calling methods (XHMM and ExomeDepth) to expand our set of single-gene CNVs and leveraged two different approaches for validating these variants (quantitative polymerase chain reaction and NanoString). RESULTS: We found a significant excess of all rare CNVs (deletions: p = .0004, duplications: p = .0006) and single-gene CNVs (deletions: p = .04, duplications: p = .03) in schizophrenia cases compared with control subjects. An expanded set of CNVs generated from integrating multiple approaches showed a significant burden of deletions in 11 of 21 gene sets previously implicated in schizophrenia and across all genes in those sets (p = .008), although no tests survived correction. We performed an extensive validation of all deletions in the significant set of voltage-gated calcium channels among CNVs called from both exome sequencing and genotyping arrays. In total, 4 exonic, single-gene deletions were validated in schizophrenia cases and none in control subjects (p = .039), of which all were identified by exome sequencing. CONCLUSIONS: These results point to the potential contribution of single-gene CNVs to schizophrenia, indicate that the utility of exome sequencing for CNV calling has yet to be maximized, and note that single-gene CNVs should be included in gene-focused studies using other classes of variation.

Contributions of common genetic variants to risk of schizophrenia among individuals of African and Latino ancestry.

Schizophrenia is a common, chronic and debilitating neuropsychiatric syndrome affecting tens of millions of individuals worldwide. While rare genetic variants play a role in the etiology of schizophrenia, most of the currently explained liability is within common variation, suggesting that variation predating the human diaspora out of Africa harbors a large fraction of the common variant attributable heritability. However, common variant association studies in schizophrenia have concentrated mainly on cohorts of European descent. We describe genome-wide association studies of 6152 cases and 3918 controls of admixed African ancestry, and of 1234 cases and 3090 controls of Latino ancestry, representing the largest such study in these populations to date. Combining results from the samples with African ancestry with summary statistics from the Psychiatric Genomics Consortium (PGC) study of schizophrenia yielded seven newly genome-wide significant loci, and we identified an additional eight loci by incorporating the results from samples with Latino ancestry. Leveraging population differences in patterns of linkage disequilibrium, we achieve improved fine-mapping resolution at 22 previously reported and 4 newly significant loci. Polygenic risk score profiling revealed improved prediction based on trans-ancestry meta-analysis results for admixed African (Nagelkerke's R2 = 0.032; liability R2 = 0.017; P < 10-52), Latino (Nagelkerke's R2 = 0.089; liability R2 = 0.021; P < 10-58), and European individuals (Nagelkerke's R2 = 0.089; liability R2 = 0.037; P < 10-113), further highlighting the advantages of incorporating data from diverse human populations.

Comparative genetic architectures of schizophrenia in East Asian and European populations.

Schizophrenia is a debilitating psychiatric disorder with approximately 1% lifetime risk globally. Large-scale schizophrenia genetic studies have reported primarily on European ancestry samples, potentially missing important biological insights. Here, we report the largest study to date of East Asian participants (22,778 schizophrenia cases and 35,362 controls), identifying 21 genome-wide-significant associations in 19 genetic loci. Common genetic variants that confer risk for schizophrenia have highly similar effects between East Asian and European ancestries (genetic correlation = 0.98 ± 0.03), indicating that the genetic basis of schizophrenia and its biology are broadly shared across populations. A fixed-effect meta-analysis including individuals from East Asian and European ancestries identified 208 significant associations in 176 genetic loci (53 novel). Trans-ancestry fine-mapping reduced the sets of candidate causal variants in 44 loci. Polygenic risk scores had reduced performance when transferred across ancestries, highlighting the importance of including sufficient samples of major ancestral groups to ensure their generalizability across populations.

Synergistic effects of common schizophrenia risk variants.

The mechanisms by which common risk variants of small effect interact to contribute to complex genetic disorders are unclear. Here, we apply a genetic approach, using isogenic human induced pluripotent stem cells, to evaluate the effects of schizophrenia (SZ)-associated common variants predicted to function as SZ expression quantitative trait loci (eQTLs). By integrating CRISPR-mediated gene editing, activation and repression technologies to study one putative SZ eQTL (FURIN rs4702) and four top-ranked SZ eQTL genes (FURIN, SNAP91, TSNARE1 and CLCN3), our platform resolves pre- and postsynaptic neuronal deficits, recapitulates genotype-dependent gene expression differences and identifies convergence downstream of SZ eQTL gene perturbations. Our observations highlight the cell-type-specific effects of common variants and demonstrate a synergistic effect between SZ eQTL genes that converges on synaptic function. We propose that the links between rare and common variants implicated in psychiatric disease risk constitute a potentially generalizable phenomenon occurring more widely in complex genetic disorders.

Global landscape and genetic regulation of RNA editing in cortical samples from individuals with schizophrenia.

RNA editing critically regulates neurodevelopment and normal neuronal function. The global landscape of RNA editing was surveyed across 364 schizophrenia cases and 383 control postmortem brain samples from the CommonMind Consortium, comprising two regions: dorsolateral prefrontal cortex and anterior cingulate cortex. In schizophrenia, RNA editing sites in genes encoding AMPA-type glutamate receptors and postsynaptic density proteins were less edited, whereas those encoding translation initiation machinery were edited more. These sites replicate between brain regions, map to 3'-untranslated regions and intronic regions, share common sequence motifs and overlap with binding sites for RNA-binding proteins crucial for neurodevelopment. These findings cross-validate in hundreds of non-overlapping dorsolateral prefrontal cortex samples. Furthermore, ~30% of RNA editing sites associate with cis-regulatory variants (editing quantitative trait loci or edQTLs). Fine-mapping edQTLs with schizophrenia risk loci revealed co-localization of eleven edQTLs with six loci. The findings demonstrate widespread altered RNA editing in schizophrenia and its genetic regulation, and suggest a causal and mechanistic role of RNA editing in schizophrenia neuropathology.

Publisher Correction: Common schizophrenia alleles are enriched in mutation-intolerant genes and in regions under strong background selection.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Publisher Correction: Gene expression imputation across multiple brain regions provides insights into schizophrenia risk.

In the HTML version of the article originally published, the author group 'The Schizophrenia Working Group of the Psychiatric Genomics Consortium' was displayed incorrectly. The error has been corrected in the HTML version of the article.

Gene expression imputation across multiple brain regions provides insights into schizophrenia risk.

Transcriptomic imputation approaches combine eQTL reference panels with large-scale genotype data in order to test associations between disease and gene expression. These genic associations could elucidate signals in complex genome-wide association study (GWAS) loci and may disentangle the role of different tissues in disease development. We used the largest eQTL reference panel for the dorso-lateral prefrontal cortex (DLPFC) to create a set of gene expression predictors and demonstrate their utility. We applied DLPFC and 12 GTEx-brain predictors to 40,299 schizophrenia cases and 65,264 matched controls for a large transcriptomic imputation study of schizophrenia. We identified 413 genic associations across 13 brain regions. Stepwise conditioning identified 67 non-MHC genes, of which 14 did not fall within previous GWAS loci. We identified 36 significantly enriched pathways, including hexosaminidase-A deficiency, and multiple porphyric disorder pathways. We investigated developmental expression patterns among the 67 non-MHC genes and identified specific groups of pre- and postnatal expression.

Contribution of Rare Copy Number Variants to Bipolar Disorder Risk Is Limited to Schizoaffective Cases.

BACKGROUND: Genetic risk for bipolar disorder (BD) is conferred through many common alleles, while a role for rare copy number variants (CNVs) is less clear. Subtypes of BD including schizoaffective disorder bipolar type (SAB), bipolar I disorder (BD I), and bipolar II disorder (BD II) differ according to the prominence and timing of psychosis, mania, and depression. The genetic factors contributing to the combination of symptoms among these subtypes are poorly understood. METHODS: Rare large CNVs were analyzed in 6353 BD cases (3833 BD I [2676 with psychosis, 850 without psychosis, and 307 with unknown psychosis history], 1436 BD II, 579 SAB, and 505 BD not otherwise specified) and 8656 controls. CNV burden and a polygenic risk score (PRS) for schizophrenia were used to evaluate the relative contributions of rare and common variants to risk of BD, BD subtypes, and psychosis. RESULTS: CNV burden did not differ between BD and controls when treated as a single diagnostic entity. However, burden in SAB was increased relative to controls (p = .001), BD I (p = .0003), and BD II (p = .0007). Burden and schizophrenia PRSs were increased in SAB compared with BD I with psychosis (CNV p = .0007, PRS p = .004), and BD I without psychosis (CNV p = .0004, PRS p = 3.9 × 10-5). Within BD I, psychosis was associated with increased schizophrenia PRSs (p = .005) but not CNV burden. CONCLUSIONS: CNV burden in BD is limited to SAB. Rare and common genetic variants may contribute differently to risk for psychosis and perhaps other classes of psychiatric symptoms.

Differential activity of transcribed enhancers in the prefrontal cortex of 537 cases with schizophrenia and controls.

Transcription at enhancers is a widespread phenomenon which produces so-called enhancer RNA (eRNA) and occurs in an activity-dependent manner. However, the role of eRNA and its utility in exploring disease-associated changes in enhancer function, and the downstream coding transcripts that they regulate, is not well established. We used transcriptomic and epigenomic data to interrogate the relationship of eRNA transcription to disease status and how genetic variants alter enhancer transcriptional activity in the human brain. We combined RNA-seq data from 537 postmortem brain samples from the CommonMind Consortium with cap analysis of gene expression and enhancer identification, using the assay for transposase-accessible chromatin followed by sequencing (ATACseq). We find 118 differentially transcribed eRNAs in schizophrenia and identify schizophrenia-associated gene/eRNA co-expression modules. Perturbations of a key module are associated with the polygenic risk scores. Furthermore, we identify genetic variants affecting expression of 927 enhancers, which we refer to as enhancer expression quantitative loci or eeQTLs. Enhancer expression patterns are consistent across studies, including differentially expressed eRNAs and eeQTLs. Combining eeQTLs with a genome-wide association study of schizophrenia identifies a genetic variant that alters enhancer function and expression of its target gene, GOLPH3L. Our novel approach to analyzing enhancer transcription is adaptable to other large-scale, non-poly-A depleted, RNA-seq studies.

Author Correction: Expression-based drug screening of neural progenitor cells from individuals with schizophrenia.

In the originally published version of this Article, the affiliation details for Eric E. Schadt and Radoslav Savic incorrectly omitted 'Sema4, a Mount Sinai venture, Stamford, Connecticut, USA'. This has been corrected in both the PDF and HTML versions of the Article.

Integrative transcriptome analyses of the aging brain implicate altered splicing in Alzheimer's disease susceptibility.

Here we use deep sequencing to identify sources of variation in mRNA splicing in the dorsolateral prefrontal cortex (DLPFC) of 450 subjects from two aging cohorts. Hundreds of aberrant pre-mRNA splicing events are reproducibly associated with Alzheimer's disease. We also generate a catalog of splicing quantitative trait loci (sQTL) effects: splicing of 3,006 genes is influenced by genetic variation. We report that altered splicing is the mechanism for the effects of the PICALM, CLU and PTK2B susceptibility alleles. Furthermore, we performed a transcriptome-wide association study and identified 21 genes with significant associations with Alzheimer's disease, many of which are found in known loci, whereas 8 are in novel loci. These results highlight the convergence of old and new genes associated with Alzheimer's disease in autophagy-lysosomal-related pathways. Overall, this study of the transcriptome of the aging brain provides evidence that dysregulation of mRNA splicing is a feature of Alzheimer's disease and is, in some cases, genetically driven.

Expression-based drug screening of neural progenitor cells from individuals with schizophrenia.

A lack of biologically relevant screening models hinders the discovery of better treatments for schizophrenia (SZ) and other neuropsychiatric disorders. Here we compare the transcriptional responses of 8 commonly used cancer cell lines (CCLs) directly with that of human induced pluripotent stem cell (hiPSC)-derived neural progenitor cells (NPCs) from 12 individuals with SZ and 12 controls across 135 drugs, generating 4320 unique drug-response transcriptional signatures. We identify those drugs that reverse post-mortem SZ-associated transcriptomic signatures, several of which also differentially regulate neuropsychiatric disease-associated genes in a cell type (hiPSC NPC vs. CCL) and/or a diagnosis (SZ vs. control)-dependent manner. Overall, we describe a proof-of-concept application of transcriptomic drug screening to hiPSC-based models, demonstrating that the drug-induced gene expression differences observed with patient-derived hiPSC NPCs are enriched for SZ biology, thereby revealing a major advantage of incorporating cell type and patient-specific platforms in drug discovery.

Evaluation of chromatin accessibility in prefrontal cortex of individuals with schizophrenia.

Schizophrenia genome-wide association studies have identified >150 regions of the genome associated with disease risk, yet there is little evidence that coding mutations contribute to this disorder. To explore the mechanism of non-coding regulatory elements in schizophrenia, we performed ATAC-seq on adult prefrontal cortex brain samples from 135 individuals with schizophrenia and 137 controls, and identified 118,152 ATAC-seq peaks. These accessible chromatin regions in the brain are highly enriched for schizophrenia SNP heritability. Accessible chromatin regions that overlap evolutionarily conserved regions exhibit an even higher heritability enrichment, indicating that sequence conservation can further refine functional risk variants. We identify few differences in chromatin accessibility between cases and controls, in contrast to thousands of age-related differential accessible chromatin regions. Altogether, we characterize chromatin accessibility in the human prefrontal cortex, the effect of schizophrenia and age on chromatin accessibility, and provide evidence that our dataset will allow for fine mapping of risk variants.

Cell-specific histone modification maps in the human frontal lobe link schizophrenia risk to the neuronal epigenome.

Risk variants for schizophrenia affect more than 100 genomic loci, yet cell- and tissue-specific roles underlying disease liability remain poorly characterized. We have generated for two cortical areas implicated in psychosis, the dorsolateral prefrontal cortex and anterior cingulate cortex, 157 reference maps from neuronal, neuron-depleted and bulk tissue chromatin for two histone marks associated with active promoters and enhancers, H3-trimethyl-Lys4 (H3K4me3) and H3-acetyl-Lys27 (H3K27ac). Differences between neuronal and neuron-depleted chromatin states were the major axis of variation in histone modification profiles, followed by substantial variability across subjects and cortical areas. Thousands of significant histone quantitative trait loci were identified in neuronal and neuron-depleted samples. Risk variants for schizophrenia, depressive symptoms and neuroticism were significantly over-represented in neuronal H3K4me3 and H3K27ac landscapes. Our Resource, sponsored by PsychENCODE and CommonMind, highlights the critical role of cell-type-specific signatures at regulatory and disease-associated noncoding sequences in the human frontal lobe.

Landscape of Conditional eQTL in Dorsolateral Prefrontal Cortex and Co-localization with Schizophrenia GWAS.

Causal genes and variants within genome-wide association study (GWAS) loci can be identified by integrating GWAS statistics with expression quantitative trait loci (eQTL) and determining which variants underlie both GWAS and eQTL signals. Most analyses, however, consider only the marginal eQTL signal, rather than dissect this signal into multiple conditionally independent signals for each gene. Here we show that analyzing conditional eQTL signatures, which could be important under specific cellular or temporal contexts, leads to improved fine mapping of GWAS associations. Using genotypes and gene expression levels from post-mortem human brain samples (n = 467) reported by the CommonMind Consortium (CMC), we find that conditional eQTL are widespread; 63% of genes with primary eQTL also have conditional eQTL. In addition, genomic features associated with conditional eQTL are consistent with context-specific (e.g., tissue-, cell type-, or developmental time point-specific) regulation of gene expression. Integrating the 2014 Psychiatric Genomics Consortium schizophrenia (SCZ) GWAS and CMC primary and conditional eQTL data reveals 40 loci with strong evidence for co-localization (posterior probability > 0.8), including six loci with co-localization of conditional eQTL. Our co-localization analyses support previously reported genes, identify novel genes associated with schizophrenia risk, and provide specific hypotheses for their functional follow-up.