scMeFormer: a transformer-based deep learning model for imputing DNA methylation states in single cells enhances the detection of epigenetic alterations in schizophrenia.

DNA methylation (DNAm), a crucial epigenetic mark, plays a key role in gene regulation, mammalian development, and various human diseases. Single-cell technologies enable the profiling of DNAm states at cytosines within the DNA sequence of individual cells, but they often suffer from limited coverage of CpG sites. In this study, we introduce scMeFormer, a transformer-based deep learning model designed to impute DNAm states for each CpG site in single cells. Through comprehensive evaluations, we demonstrate the superior performance of scMeFormer compared to alternative models across four single-nucleus DNAm datasets generated by distinct technologies. Remarkably, scMeFormer exhibits high-fidelity imputation, even when dealing with significantly reduced coverage, as low as 10% of the original CpG sites. Furthermore, we applied scMeFormer to a single-nucleus DNAm dataset generated from the prefrontal cortex of four schizophrenia patients and four neurotypical controls. This enabled the identification of thousands of differentially methylated regions associated with schizophrenia that would have remained undetectable without imputation and added granularity to our understanding of epigenetic alterations in schizophrenia within specific cell types. Our study highlights the power of deep learning in imputing DNAm states in single cells, and we expect scMeFormer to be a valuable tool for single-cell DNAm studies.

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.

Electrophysiological measures from human iPSC-derived neurons are associated with schizophrenia clinical status and predict individual cognitive performance.

Neurons derived from human induced pluripotent stem cells (hiPSCs) have been used to model basic cellular aspects of neuropsychiatric disorders, but the relationship between the emergent phenotypes and the clinical characteristics of donor individuals has been unclear. We analyzed RNA expression and indices of cellular function in hiPSC-derived neural progenitors and cortical neurons generated from 13 individuals with high polygenic risk scores (PRSs) for schizophrenia (SCZ) and a clinical diagnosis of SCZ, along with 15 neurotypical individuals with low PRS. We identified electrophysiological measures in the patient-derived neurons that implicated altered Na+ channel function, action potential interspike interval, and gamma-aminobutyric acid-ergic neurotransmission. Importantly, electrophysiological measures predicted cardinal clinical and cognitive features found in these SCZ patients. The identification of basic neuronal physiological properties related to core clinical characteristics of illness is a potentially critical step in generating leads for novel therapeutics.

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.

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.

Dissecting transcriptomic signatures of neuronal differentiation and maturation using iPSCs.

Human induced pluripotent stem cells (hiPSCs) are a powerful model of neural differentiation and maturation. We present a hiPSC transcriptomics resource on corticogenesis from 5 iPSC donor and 13 subclonal lines across 9 time points over 5 broad conditions: self-renewal, early neuronal differentiation, neural precursor cells (NPCs), assembled rosettes, and differentiated neuronal cells. We identify widespread changes in the expression of both individual features and global patterns of transcription. We next demonstrate that co-culturing human NPCs with rodent astrocytes results in mutually synergistic maturation, and that cell type-specific expression data can be extracted using only sequencing read alignments without cell sorting. We lastly adapt a previously generated RNA deconvolution approach to single-cell expression data to estimate the relative neuronal maturity of iPSC-derived neuronal cultures and human brain tissue. Using many public datasets, we demonstrate neuronal cultures are maturationally heterogeneous but contain subsets of neurons more mature than previously observed.

iPSC-derived homogeneous populations of developing schizophrenia cortical interneurons have compromised mitochondrial function.

Schizophrenia (SCZ) is a neurodevelopmental disorder. Thus, studying pathogenetic mechanisms underlying SCZ requires studying the development of brain cells. Cortical interneurons (cINs) are consistently observed to be abnormal in SCZ postmortem brains. These abnormalities may explain altered gamma oscillation and cognitive function in patients with SCZ. Of note, currently used antipsychotic drugs ameliorate psychosis, but they are not very effective in reversing cognitive deficits. Characterizing mechanisms of SCZ pathogenesis, especially related to cognitive deficits, may lead to improved treatments. We generated homogeneous populations of developing cINs from 15 healthy control (HC) iPSC lines and 15 SCZ iPSC lines. SCZ cINs, but not SCZ glutamatergic neurons, show dysregulated Oxidative Phosphorylation (OxPhos) related gene expression, accompanied by compromised mitochondrial function. The OxPhos deficit in cINs could be reversed by Alpha Lipoic Acid/Acetyl-L-Carnitine (ALA/ALC) but not by other chemicals previously identified as increasing mitochondrial function. The restoration of mitochondrial function by ALA/ALC was accompanied by a reversal of arborization deficits in SCZ cINs. OxPhos abnormality, even in the absence of any circuit environment with other neuronal subtypes, appears to be an intrinsic deficit in SCZ cINs.

Identification and prioritization of gene sets associated with schizophrenia risk by co-expression network analysis in human brain.

Schizophrenia polygenic risk is plausibly manifested by complex transcriptional dysregulation in the brain, involving networks of co-expressed and functionally related genes. The main purpose of this study was to identify and prioritize co-expressed gene sets in a hierarchical manner, based on the strength of the relationships with clinical diagnosis and with polygenic risk for schizophrenia. Weighted Gene Co-expression Network Analysis (WGCNA) was applied to RNA-quality-adjusted DLPFC RNA-Seq data from the LIBD Postmortem Human Brain Repository (90 controls, 74 schizophrenia cases; all Caucasians) to construct co-expression networks and detect "modules" of co-expressed genes. After multiple internal and external validation procedures, modules of selected interest were tested for enrichment in biological ontologies, for association with schizophrenia polygenic risk scores (PRSs) and with diagnosis, and also for enrichment in genes within the significant GWAS loci reported by the Psychiatric Genomic Consortium (PGC2). The association between schizophrenia genetic signals and modules of co-expression converged on one module showing not only a significant association with both diagnosis and PRS but also significant overlap with 36 PGC2 loci genes, deemed the strongest candidates for drug targets. This module contained many genes involved in synaptic signaling and neuroplasticity. Fifty-three PGC2 genes were in modules associated only with diagnosis and 59 in modules unrelated to diagnosis or PRS. Our study highlights complex relationships between gene co-expression networks in the brain and clinical state and polygenic risk for SCZ and provides a strategy for using this information in selecting and prioritizing potentially targetable gene sets for therapeutic drug development.

Correction: iPSC-derived homogeneous populations of developing schizophrenia cortical interneurons have compromised mitochondrial function.

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

Schizophrenia risk variants influence multiple classes of transcripts of sorting nexin 19 (SNX19).

Genome-wide association studies (GWAS) have identified many genomic loci associated with risk for schizophrenia, but unambiguous identification of the relationship between disease-associated variants and specific genes, and in particular their effect on risk conferring transcripts, has proven difficult. To better understand the specific molecular mechanism(s) at the schizophrenia locus in 11q25, we undertook cis expression quantitative trait loci (cis-eQTL) mapping for this 2 megabase genomic region using postmortem human brain samples. To comprehensively assess the effects of genetic risk upon local expression, we evaluated multiple transcript features: genes, exons, and exon-exon junctions in multiple brain regions-dorsolateral prefrontal cortex (DLPFC), hippocampus, and caudate. Genetic risk variants strongly associated with expression of SNX19 transcript features that tag multiple rare classes of SNX19 transcripts, whereas they only weakly affected expression of an exon-exon junction that tags the majority of abundant transcripts. The most prominent class of SNX19 risk-associated transcripts is predicted to be overexpressed, defined by an exon-exon splice junction between exons 8 and 10 (junc8.10) and that is predicted to encode proteins that lack the characteristic nexin C terminal domain. Risk alleles were also associated with either increased or decreased expression of multiple additional classes of transcripts. With RACE, molecular cloning, and long read sequencing, we found a number of novel SNX19 transcripts that further define the set of potential etiological transcripts. We explored epigenetic regulation of SNX19 expression and found that DNA methylation at CpG sites near the primary transcription start site and within exon 2 partially mediate the effects of risk variants on risk-associated expression. ATAC sequencing revealed that some of the most strongly risk-associated SNPs are located within a region of open chromatin, suggesting a nearby regulatory element is involved. These findings indicate a potentially complex molecular etiology, in which risk alleles for schizophrenia generate epigenetic alterations and dysregulation of multiple classes of SNX19 transcripts.

Regional Heterogeneity in Gene Expression, Regulation, and Coherence in the Frontal Cortex and Hippocampus across Development and Schizophrenia.

The hippocampus formation, although prominently implicated in schizophrenia pathogenesis, has been overlooked in large-scale genomics efforts in the schizophrenic brain. We performed RNA-seq in hippocampi and dorsolateral prefrontal cortices (DLPFCs) from 551 individuals (286 with schizophrenia). We identified substantial regional differences in gene expression and found widespread developmental differences that were independent of cellular composition. We identified 48 and 245 differentially expressed genes (DEGs) associated with schizophrenia within the hippocampus and DLPFC, with little overlap between the brain regions. 124 of 163 (76.6%) of schizophrenia GWAS risk loci contained eQTLs in any region. Transcriptome-wide association studies in each region identified many novel schizophrenia risk features that were brain region-specific. Last, we identified potential molecular correlates of in vivo evidence of altered prefrontal-hippocampal functional coherence in schizophrenia. These results underscore the complexity and regional heterogeneity of the transcriptional correlates of schizophrenia and offer new insights into potentially causative biology.

Dysregulated protocadherin-pathway activity as an intrinsic defect in induced pluripotent stem cell-derived cortical interneurons from subjects with schizophrenia.

We generated cortical interneurons (cINs) from induced pluripotent stem cells derived from 14 healthy controls and 14 subjects with schizophrenia. Both healthy control cINs and schizophrenia cINs were authentic, fired spontaneously, received functional excitatory inputs from host neurons, and induced GABA-mediated inhibition in host neurons in vivo. However, schizophrenia cINs had dysregulated expression of protocadherin genes, which lie within documented schizophrenia loci. Mice lacking protocadherin-α showed defective arborization and synaptic density of prefrontal cortex cINs and behavioral abnormalities. Schizophrenia cINs similarly showed defects in synaptic density and arborization that were reversed by inhibitors of protein kinase C, a downstream kinase in the protocadherin pathway. These findings reveal an intrinsic abnormality in schizophrenia cINs in the absence of any circuit-driven pathology. They also demonstrate the utility of homogenous and functional populations of a relevant neuronal subtype for probing pathogenesis mechanisms during development.

Polygenic risk score increases schizophrenia liability through cognition-relevant pathways.

Cognitive deficit is thought to represent, at least in part, genetic mechanisms of risk for schizophrenia, with recent evidence from statistical modelling of twin data suggesting direct causality from the former to the latter. However, earlier evidence was based on inferences from twin not molecular genetic data and it is unclear how much genetic influence 'passes through' cognition on the way to diagnosis. Thus, we included direct measurements of genetic risk (e.g. schizophrenia polygenic risk scores) in causation models to assess the extent to which cognitive deficit mediates some of the effect of polygenic risk scores on the disorder. Causal models of family data tested relationships among key variables and allowed parsing of genetic variance components. Polygenic risk scores were calculated from summary statistics from the current largest genome-wide association study of schizophrenia and were represented as a latent trait. Cognition was also modelled as a latent trait. Participants were 1313 members of 1078 families: 416 patients with schizophrenia, 290 unaffected siblings, and 607 controls. Modelling supported earlier findings that cognitive deficit has a putatively causal role in schizophrenia. In total, polygenic risk score explained 8.07% [confidence interval (CI) 5.45-10.74%] of schizophrenia risk in our sample. Of this, more than a third (2.71%, CI 2.41-3.85%) of the polygenic risk score influence was mediated through cognition paths, exceeding the direct influence of polygenic risk score on schizophrenia risk (1.43%, CI 0.46-3.08%). The remainder of the polygenic risk score influence (3.93%, CI 2.37-4.48%) reflected reciprocal causation between schizophrenia liability and cognition (e.g. mutual influences in a cyclical manner). Analysis of genetic variance components of schizophrenia liability indicated that 26.87% (CI 21.45-32.57%) was associated with cognition-related pathways not captured by polygenic risk score. The remaining variance in schizophrenia was through pathways other than cognition-related and polygenic risk score. Although our results are based on inference through statistical modelling and do not provide an absolute proof of causality, we find that cognition pathways mediate a significant part of the influence of cumulative genetic risk on schizophrenia. We estimate from our model that 33.51% (CI 27.34-43.82%) of overall genetic risk is mediated through influences on cognition, but this requires further studies and analyses as the genetics of schizophrenia becomes better characterized.

Publisher Correction: Variations in Dysbindin-1 are associated with cognitive response to antipsychotic drug treatment.

In the original version of this Article, references in the Methods section incorrectly referred to references in the Supplementary References section. The relevant references (now numbered 20, 27, 42, 47, 69-80) have been removed from the Supplementary References section of the Supplementary Information file and added to the References section of the main manuscript, in both the PDF and HTML versions of the Article.

Developmental and genetic regulation of the human cortex transcriptome illuminate schizophrenia pathogenesis.

Genome-wide association studies have identified 108 schizophrenia risk loci, but biological mechanisms for individual loci are largely unknown. Using developmental, genetic and illness-based RNA sequencing expression analysis in human brain, we characterized the human brain transcriptome around these loci and found enrichment for developmentally regulated genes with novel examples of shifting isoform usage across pre- and postnatal life. We found widespread expression quantitative trait loci (eQTLs), including many with transcript specificity and previously unannotated sequence that were independently replicated. We leveraged this general eQTL database to show that 48.1% of risk variants for schizophrenia associate with nearby expression. We lastly found 237 genes significantly differentially expressed between patients and controls, which replicated in an independent dataset, implicated synaptic processes, and were strongly regulated in early development. These findings together offer genetics- and diagnosis-related targets for better modeling of schizophrenia risk. This resource is publicly available at http://eqtl.brainseq.org/phase1 .

Multi-Trait Analysis of GWAS and Biological Insights Into Cognition: A Response to Hill (2018).

Hill (Twin Research and Human Genetics, Vol. 21, 2018, 84-88) presented a critique of our recently published paper in Cell Reports entitled 'Large-Scale Cognitive GWAS Meta-Analysis Reveals Tissue-Specific Neural Expression and Potential Nootropic Drug Targets' (Lam et al., Cell Reports, Vol. 21, 2017, 2597-2613). Specifically, Hill offered several interrelated comments suggesting potential problems with our use of a new analytic method called Multi-Trait Analysis of GWAS (MTAG) (Turley et al., Nature Genetics, Vol. 50, 2018, 229-237). In this brief article, we respond to each of these concerns. Using empirical data, we conclude that our MTAG results do not suffer from 'inflation in the FDR [false discovery rate]', as suggested by Hill (Twin Research and Human Genetics, Vol. 21, 2018, 84-88), and are not 'more relevant to the genetic contributions to education than they are to the genetic contributions to intelligence'.

Variations in Dysbindin-1 are associated with cognitive response to antipsychotic drug treatment.

Antipsychotics are the most widely used medications for the treatment of schizophrenia spectrum disorders. While such drugs generally ameliorate positive symptoms, clinical responses are highly variable in terms of negative symptoms and cognitive impairments. However, predictors of individual responses have been elusive. Here, we report a pharmacogenetic interaction related to a core cognitive dysfunction in patients with schizophrenia. We show that genetic variations reducing dysbindin-1 expression can identify individuals whose executive functions respond better to antipsychotic drugs, both in humans and in mice. Multilevel ex vivo and in vivo analyses in postmortem human brains and genetically modified mice demonstrate that such interaction between antipsychotics and dysbindin-1 is mediated by an imbalance between the short and long isoforms of dopamine D2 receptors, leading to enhanced presynaptic D2 function within the prefrontal cortex. These findings reveal one of the pharmacodynamic mechanisms underlying individual cognitive response to treatment in patients with schizophrenia, suggesting a potential approach for improving the use of antipsychotic drugs.

Convergence of placenta biology and genetic risk for schizophrenia.

Defining the environmental context in which genes enhance disease susceptibility can provide insight into the pathogenesis of complex disorders. We report that the intra-uterine environment modulates the association of schizophrenia with genomic risk (in this study, genome-wide association study-derived polygenic risk scores (PRSs)). In independent samples from the United States, Italy, and Germany, the liability of schizophrenia explained by PRS is more than five times greater in the presence of early-life complications (ELCs) compared with their absence. Patients with ELC histories have significantly higher PRS than patients without ELC histories, which is confirmed in additional samples from Germany and Japan. The gene set composed of schizophrenia loci that interact with ELCs is highly expressed in placenta, is differentially expressed in placentae from complicated in comparison with normal pregnancies, and is differentially upregulated in placentae from male compared with female offspring. Pathway analyses reveal that genes driving the PRS-ELC interaction are involved in cellular stress response; genes that do not drive such interaction implicate orthogonal biological processes (for example, synaptic function). We conclude that a subset of the most significant genetic variants associated with schizophrenia converge on a developmental trajectory sensitive to events that affect the placental response to stress, which may offer insights into sex biases and primary prevention.

Schizophrenia polygenic risk score predicts mnemonic hippocampal activity.

The use of polygenic risk scores has become a practical translational approach to investigating the complex genetic architecture of schizophrenia, but the link between polygenic risk scores and pathophysiological components of this disorder has been the subject of limited research. We investigated in healthy volunteers whether schizophrenia polygenic risk score predicts hippocampal activity during simple memory encoding, which has been proposed as a risk-associated intermediate phenotype of schizophrenia. We analysed the relationship between polygenic risk scores and hippocampal activity in a discovery sample of 191 unrelated healthy volunteers from the USA and in two independent replication samples of 76 and 137 healthy unrelated participants from Europe and the USA, respectively. Polygenic risk scores for each individual were calculated as the sum of the imputation probability of reference alleles weighted by the natural log of odds ratio from the recent schizophrenia genome-wide association study. We examined hippocampal activity during simple memory encoding of novel visual stimuli assessed using blood oxygen level-dependent functional MRI. Polygenic risk scores were significantly associated with hippocampal activity in the discovery sample [P = 0.016, family-wise error (FWE) corrected within Anatomical Automatic Labeling (AAL) bilateral hippocampal-parahippocampal mask] and in both replication samples (P = 0.033, FWE corrected within AAL right posterior hippocampal-parahippocampal mask in Bari sample, and P = 0.002 uncorrected in the Duke Neurogenetics Study sample). The relationship between polygenic risk scores and hippocampal activity was consistently negative, i.e. lower hippocampal activity in individuals with higher polygenic risk scores, consistent with previous studies reporting decreased hippocampal-parahippocampal activity during declarative memory tasks in patients with schizophrenia and in their healthy siblings. Polygenic risk scores accounted for more than 8% of variance in hippocampal activity during memory encoding in discovery sample. We conclude that polygenic risk scores derived from the most recent schizophrenia genome-wide association study predict significant variability in hippocampal activity during memory encoding in healthy participants. Our findings validate mnemonic hippocampal activity as a genetic risk associated intermediate phenotype of schizophrenia, indicating that the aggregate neurobiological effect of schizophrenia risk alleles converges on this pattern of neural activity.awy004media15749593779001.

Large-Scale Cognitive GWAS Meta-Analysis Reveals Tissue-Specific Neural Expression and Potential Nootropic Drug Targets.

Here, we present a large (n = 107,207) genome-wide association study (GWAS) of general cognitive ability ("g"), further enhanced by combining results with a large-scale GWAS of educational attainment. We identified 70 independent genomic loci associated with general cognitive ability. Results showed significant enrichment for genes causing Mendelian disorders with an intellectual disability phenotype. Competitive pathway analysis implicated the biological processes of neurogenesis and synaptic regulation, as well as the gene targets of two pharmacologic agents: cinnarizine, a T-type calcium channel blocker, and LY97241, a potassium channel inhibitor. Transcriptome-wide and epigenome-wide analysis revealed that the implicated loci were enriched for genes expressed across all brain regions (most strongly in the cerebellum). Enrichment was exclusive to genes expressed in neurons but not oligodendrocytes or astrocytes. Finally, we report genetic correlations between cognitive ability and disparate phenotypes including psychiatric disorders, several autoimmune disorders, longevity, and maternal age at first birth.