Molecular landscapes of human hippocampal immature neurons across lifespan.

Immature dentate granule cells (imGCs) arising from adult hippocampal neurogenesis contribute to plasticity and unique brain functions in rodents1,2 and are dysregulated in multiple human neurological disorders3-5. Little is known about the molecular characteristics of adult human hippocampal imGCs, and even their existence is under debate1,6-8. Here we performed single-nucleus RNA sequencing aided by a validated machine learning-based analytic approach to identify imGCs and quantify their abundance in the human hippocampus at different stages across the lifespan. We identified common molecular hallmarks of human imGCs across the lifespan and observed age-dependent transcriptional dynamics in human imGCs that suggest changes in cellular functionality, niche interactions and disease relevance, that differ from those in mice9. We also found a decreased number of imGCs with altered gene expression in Alzheimer's disease. Finally, we demonstrated the capacity for neurogenesis in the adult human hippocampus with the presence of rare dentate granule cell fate-specific proliferating neural progenitors and with cultured surgical specimens. Together, our findings suggest the presence of a substantial number of imGCs in the adult human hippocampus via low-frequency de novo generation and protracted maturation, and our study reveals their molecular properties across the lifespan and in Alzheimer's disease.

Large-scale cellular-resolution gene profiling in human neocortex reveals species-specific molecular signatures.

Although there have been major advances in elucidating the functional biology of the human brain, relatively little is known of its cellular and molecular organization. Here we report a large-scale characterization of the expression of ∼1,000 genes important for neural functions by in situ hybridization at a cellular resolution in visual and temporal cortices of adult human brains. These data reveal diverse gene expression patterns and remarkable conservation of each individual gene's expression among individuals (95%), cortical areas (84%), and between human and mouse (79%). A small but substantial number of genes (21%) exhibited species-differential expression. Distinct molecular signatures, comprised of genes both common between species and unique to each, were identified for each major cortical cell type. The data suggest that gene expression profile changes may contribute to differential cortical function across species, and in particular, a shift from corticosubcortical to more predominant corticocortical communications in the human brain.

Investigating trait variability of gene co-expression network architecture in brain by controlling for genomic risk of schizophrenia.

The effect of schizophrenia (SCZ) genetic risk on gene expression in brain remains elusive. A popular approach to this problem has been the application of gene co-expression network algorithms (e.g., WGCNA). To improve reliability with this method it is critical to remove unwanted sources of variance while also preserving biological signals of interest. In this WCGNA study of RNA-Seq data from postmortem prefrontal cortex (78 neurotypical donors, EUR ancestry), we tested the effects of SCZ genetic risk on co-expression networks. Specifically, we implemented a novel design in which gene expression was adjusted by linear regression models to preserve or remove variance explained by biological signal of interest (GWAS genomic scores for SCZ risk-(GS-SCZ), and genomic scores- GS of height (GS-Ht) as a negative control), while removing variance explained by covariates of non-interest. We calculated co-expression networks from adjusted expression (GS-SCZ and GS-Ht preserved or removed), and consensus between them (representative of a "background" network free of genomic scores effects). We then tested the overlap between GS-SCZ preserved modules and background networks reasoning that modules with reduced overlap would be most affected by GS-SCZ biology. Additionally, we tested these modules for convergence of SCZ risk (i.e., enrichment in PGC3 SCZ GWAS priority genes, enrichment in SCZ risk heritability and relevant biological ontologies. Our results highlight key aspects of GS-SCZ effects on brain co-expression networks, specifically: 1) preserving/removing SCZ genetic risk alters the co-expression modules; 2) biological pathways enriched in modules affected by GS-SCZ implicate processes of transcription, translation and metabolism that converge to influence synaptic transmission; 3) priority PGC3 SCZ GWAS genes and SCZ risk heritability are enriched in modules associated with GS-SCZ effects. Overall, our results indicate that gene co-expression networks that selectively integrate information about genetic risk can reveal novel combinations of biological pathways involved in schizophrenia.

Deep learning predicts DNA methylation regulatory variants in the human brain and elucidates the genetics of psychiatric disorders.

There is growing evidence for the role of DNA methylation (DNAm) quantitative trait loci (mQTLs) in the genetics of complex traits, including psychiatric disorders. However, due to extensive linkage disequilibrium (LD) of the genome, it is challenging to identify causal genetic variations that drive DNAm levels by population-based genetic association studies. This limits the utility of mQTLs for fine-mapping risk loci underlying psychiatric disorders identified by genome-wide association studies (GWAS). Here we present INTERACT, a deep learning model that integrates convolutional neural networks with transformer, to predict effects of genetic variations on DNAm levels at CpG sites in the human brain. We show that INTERACT-derived DNAm regulatory variants are not confounded by LD, are concentrated in regulatory genomic regions in the human brain, and are convergent with mQTL evidence from genetic association analysis. We further demonstrate that predicted DNAm regulatory variants are enriched for heritability of brain-related traits and improve polygenic risk prediction for schizophrenia across diverse ancestry samples. Finally, we applied predicted DNAm regulatory variants for fine-mapping schizophrenia GWAS risk loci to identify potential novel risk genes. Our study shows the power of a deep learning approach to identify functional regulatory variants that may elucidate the genetic basis of complex traits.

Orthogonal approaches required to measure proteasome composition and activity in mammalian brain tissue.

Proteasomes are large macromolecular complexes with multiple distinct catalytic activities that are each vital to human brain health and disease. Despite their importance, standardized approaches to investigate proteasomes have not been universally adapted. Here, we describe pitfalls and define straightforward orthogonal biochemical approaches essential to measure and understand changes in proteasome composition and activity in the mammalian central nervous system. Through our experimentation in the mammalian brain, we determined an abundance of catalytically active proteasomes exist with and without a 19S cap(s), the regulatory particle essential for ubiquitin-dependent degradation. Moreover, we learned that in-cell measurements using activity-based probes (ABPs) are more sensitive in determining the available activity of the 20S proteasome without the 19S cap and in measuring individual catalytic subunit activities of each ß subunit within all neuronal proteasomes. Subsequently, applying these tools to human brain samples, we were surprised to find that post-mortem tissue retained little to no 19S-capped proteasome, regardless of age, sex, or disease state. In comparing brain tissues (parahippocampal gyrus) from patients with Alzheimer's disease (AD) and unaffected individuals, the available 20S proteasome activity was significantly elevated in severe cases of AD, an observation not previously noted. Taken together, our study establishes standardized approaches for the comprehensive investigation of proteasomes in mammalian brain tissue, and we reveal new insight into brain proteasome biology.

Deoxyhypusine hydroxylase: A novel therapeutic target differentially expressed in short-term vs long-term survivors of glioblastoma.

Glioblastoma (GB) is the most aggressive neoplasm of the brain. Poor prognosis is mainly attributed to tumor heterogeneity, invasiveness and drug resistance. Only a small fraction of GB patients survives longer than 24 months from the time of diagnosis (ie, long-term survivors [LTS]). In our study, we aimed to identify molecular markers associated with favorable GB prognosis as a basis to develop therapeutic applications to improve patients' outcome. We have recently assembled a proteogenomic dataset of 87 GB clinical samples of varying survival rates. Following RNA-seq and mass spectrometry (MS)-based proteomics analysis, we identified several differentially expressed genes and proteins, including some known cancer-related pathways and some less established that showed higher expression in short-term (<6 months) survivors (STS) compared to LTS. One such target found was deoxyhypusine hydroxylase (DOHH), which is known to be involved in the biosynthesis of hypusine, an unusual amino acid essential for the function of the eukaryotic translation initiation factor 5A (eIF5A), which promotes tumor growth. We consequently validated DOHH overexpression in STS samples by quantitative polymerase chain reaction (qPCR) and immunohistochemistry. We further showed robust inhibition of proliferation, migration and invasion of GB cells following silencing of DOHH with short hairpin RNA (shRNA) or inhibition of its activity with small molecules, ciclopirox and deferiprone. Moreover, DOHH silencing led to significant inhibition of tumor progression and prolonged survival in GB mouse models. Searching for a potential mechanism by which DOHH promotes tumor aggressiveness, we found that it supports the transition of GB cells to a more invasive phenotype via epithelial-mesenchymal transition (EMT)-related pathways.

Characterizing the nuclear and cytoplasmic transcriptomes in developing and mature human cortex uncovers new insight into psychiatric disease gene regulation.

Transcriptome compartmentalization by the nuclear membrane provides both stochastic and functional buffering of transcript activity in the cytoplasm, and has recently been implicated in neurodegenerative disease processes. Although many mechanisms regulating transcript compartmentalization are also prevalent in brain development, the extent to which subcellular localization differs as the brain matures has yet to be addressed. To characterize the nuclear and cytoplasmic transcriptomes during brain development, we sequenced both RNA fractions from homogenate prenatal and adult human postmortem cortex using poly(A)+ and Ribo-Zero library preparation methods. We find that while many genes are differentially expressed by fraction and developmental expression changes are similarly detectable in nuclear and cytoplasmic RNA, the compartmented transcriptomes become more distinct as the brain matures, perhaps reflecting increased utilization of nuclear retention as a regulatory strategy in adult brain. We examined potential mechanisms of this developmental divergence including alternative splicing, RNA editing, nuclear pore composition, RNA-binding protein motif enrichment, and RNA secondary structure. Intron retention is associated with greater nuclear abundance in a subset of transcripts, as is enrichment for several splicing factor binding motifs. Finally, we examined disease association with fraction-regulated gene sets and found nuclear-enriched genes were also preferentially enriched in gene sets associated with neurodevelopmental psychiatric disorders. These results suggest that although gene-level expression is globally comparable between fractions, nuclear retention of transcripts may play an underappreciated role in developmental regulation of gene expression in brain, particularly in genes whose dysregulation is related to neuropsychiatric disorders.

SUFI: an automated approach to spectral unmixing of fluorescent multiplex images captured in mouse and post-mortem human brain tissues.

BACKGROUND: Multispectral fluorescence imaging coupled with linear unmixing is a form of image data collection and analysis that allows for measuring multiple molecular signals in a single biological sample. Multiple fluorescent dyes, each measuring a unique molecule, are simultaneously measured and subsequently "unmixed" to provide a read-out for each molecular signal. This strategy allows for measuring highly multiplexed signals in a single data capture session, such as multiple proteins or RNAs in tissue slices or cultured cells, but can often result in mixed signals and bleed-through problems across dyes. Existing spectral unmixing algorithms are not optimized for challenging biological specimens such as post-mortem human brain tissue, and often require manual intervention to extract spectral signatures. We therefore developed an intuitive, automated, and flexible package called SUFI: spectral unmixing of fluorescent images. RESULTS: This package unmixes multispectral fluorescence images by automating the extraction of spectral signatures using vertex component analysis, and then performs one of three unmixing algorithms derived from remote sensing. We evaluate these remote sensing algorithms' performances on four unique biological datasets and compare the results to unmixing results obtained using ZEN Black software (Zeiss). We lastly integrate our unmixing pipeline into the computational tool dotdotdot, which is used to quantify individual RNA transcripts at single cell resolution in intact tissues and perform differential expression analysis, and thereby provide an end-to-end solution for multispectral fluorescence image analysis and quantification. CONCLUSIONS: In summary, we provide a robust, automated pipeline to assist biologists with improved spectral unmixing of multispectral fluorescence images.

Associations of psychiatric disease and ageing with FKBP5 expression converge on superficial layer neurons of the neocortex.

Identification and characterisation of novel targets for treatment is a priority in the field of psychiatry. FKBP5 is a gene with decades of evidence suggesting its pathogenic role in a subset of psychiatric patients, with potential to be leveraged as a therapeutic target for these individuals. While it is widely reported that FKBP5/FKBP51 mRNA/protein (FKBP5/1) expression is impacted by psychiatric disease state, risk genotype and age, it is not known in which cell types and sub-anatomical areas of the human brain this occurs. This knowledge is critical to propel FKBP5/1-targeted treatment development. Here, we performed an extensive, large-scale postmortem study (n = 1024) of FKBP5/1, examining neocortical areas (BA9, BA11 and ventral BA24/BA24a) derived from subjects that lived with schizophrenia, major depression or bipolar disorder. With an extensive battery of RNA (bulk RNA sequencing, single-nucleus RNA sequencing, microarray, qPCR, RNAscope) and protein (immunoblot, immunohistochemistry) analysis approaches, we thoroughly investigated the effects of disease state, ageing and genotype on cortical FKBP5/1 expression including in a cell type-specific manner. We identified consistently heightened FKBP5/1 levels in psychopathology and with age, but not genotype, with these effects strongest in schizophrenia. Using single-nucleus RNA sequencing (snRNAseq; BA9 and BA11) and targeted histology (BA9, BA24a), we established that these disease and ageing effects on FKBP5/1 expression were most pronounced in excitatory superficial layer neurons of the neocortex, and this effect appeared to be consistent in both the granular and agranular areas examined. We then found that this increase in FKBP5 levels may impact on synaptic plasticity, as FKBP5 gex levels strongly and inversely correlated with dendritic mushroom spine density and brain-derived neurotrophic factor (BDNF) levels in superficial layer neurons in BA11. These findings pinpoint a novel cellular and molecular mechanism that has potential to open a new avenue of FKBP51 drug development to treat cognitive symptoms in psychiatric disorders.

Molecular phenotypes associated with antipsychotic drugs in the human caudate nucleus.

Antipsychotic drugs are the current first-line of treatment for schizophrenia and other psychotic conditions. However, their molecular effects on the human brain are poorly studied, due to difficulty of tissue access and confounders associated with disease status. Here we examine differences in gene expression and DNA methylation associated with positive antipsychotic drug toxicology status in the human caudate nucleus. We find no genome-wide significant differences in DNA methylation, but abundant differences in gene expression. These gene expression differences are overall quite similar to gene expression differences between schizophrenia cases and controls. Interestingly, gene expression differences based on antipsychotic toxicology are different between brain regions, potentially due to affected cell type differences. We finally assess similarities with effects in a mouse model, which finds some overlapping effects but many differences as well. As a first look at the molecular effects of antipsychotics in the human brain, the lack of epigenetic effects is unexpected, possibly because long term treatment effects may be relatively stable for extended periods.

Single molecule in situ hybridization reveals distinct localizations of schizophrenia risk-related transcripts SNX19 and AS3MT in human brain.

Genome-wide association studies have identified single nucleotide polymorphisms (SNPs) associated with schizophrenia risk. Integration of RNA-sequencing data from postmortem human brains with these risk SNPs identified transcripts associated with increased schizophrenia susceptibility, including a class of exon 9-spliced isoforms of Sorting nexin-19 (SNX19d9) and an isoform of Arsenic methyltransferase (AS3MT) splicing out exons 2 and 3 (AS3MTd2d3). However, the biological function of these transcript variants is unclear. Defining the cell types where these risk transcripts are dominantly expressed is an important step to understand function, in prioritizing specific cell types and/or neural pathways in subsequent studies. To identify the cell type-specific localization of SNX19 and AS3MT in the human dorsolateral prefrontal cortex (DLPFC), we used single-molecule in situ hybridization techniques combined with automated quantification and machine learning approaches to analyze 10 postmortem brains of neurotypical individuals. These analyses revealed that both pan-SNX19 and pan-AS3MT were more highly expressed in neurons than non-neurons in layers II/III and VI of DLPFC. Furthermore, pan-SNX19 was preferentially expressed in glutamatergic neurons, while pan-AS3MT was preferentially expressed in GABAergic neurons. Finally, we utilized duplex BaseScope technology, to delineate the localization of SNX19d9 and AS3MTd2d3 splice variants, revealing consistent trends in spatial gene expression among pan-transcripts and schizophrenia risk-related transcript variants. These findings demonstrate that schizophrenia risk transcripts have distinct localization patterns in the healthy human brains, and suggest that SNX19 transcripts might disrupt the normal function of glutamatergic neurons, while AS3MT may lead to disturbances in the GABAergic system in the pathophysiology of schizophrenia.

dotdotdot: an automated approach to quantify multiplex single molecule fluorescent in situ hybridization (smFISH) images in complex tissues.

Multiplex single-molecule fluorescent in situ hybridization (smFISH) is a powerful method for validating RNA sequencing and emerging spatial transcriptomic data, but quantification remains a computational challenge. We present a framework for generating and analyzing smFISH data in complex tissues while overcoming autofluorescence and increasing multiplexing capacity. We developed dotdotdot (https://github.com/LieberInstitute/dotdotdot) as a corresponding software package to quantify RNA transcripts in single nuclei and perform differential expression analysis. We first demonstrate robustness of our platform in single mouse neurons by quantifying differential expression of activity-regulated genes. We then quantify spatial gene expression in human dorsolateral prefrontal cortex (DLPFC) using spectral imaging and dotdotdot to mask lipofuscin autofluorescence. We lastly apply machine learning to predict cell types and perform downstream cell type-specific expression analysis. In summary, we provide experimental workflows, imaging acquisition and analytic strategies for quantification and biological interpretation of smFISH data in complex tissues.

Somatic LINE-1 retrotransposition in cortical neurons and non-brain tissues of Rett patients and healthy individuals.

Mounting evidence supports that LINE-1 (L1) retrotransposition can occur postzygotically in healthy and diseased human tissues, contributing to genomic mosaicism in the brain and other somatic tissues of an individual. However, the genomic distribution of somatic human-specific LINE-1 (L1Hs) insertions and their potential impact on carrier cells remain unclear. Here, using a PCR-based targeted bulk sequencing approach, we profiled 9,181 somatic insertions from 20 postmortem tissues from five Rett patients and their matched healthy controls. We identified and validated somatic L1Hs insertions in both cortical neurons and non-brain tissues. In Rett patients, somatic insertions were significantly depleted in exons-mainly contributed by long genes-than healthy controls, implying that cells carrying MECP2 mutations might be defenseless against a second exonic L1Hs insertion. We observed a significant increase of somatic L1Hs insertions in the brain compared with non-brain tissues from the same individual. Compared to germline insertions, somatic insertions were less sense-depleted to transcripts, indicating that they underwent weaker selective pressure on the orientation of insertion. Our observations demonstrate that somatic L1Hs insertions contribute to genomic diversity and MeCP2 dysfunction alters their genomic patterns in Rett patients.

TRIOBP-1 Protein Aggregation Exists in Both Major Depressive Disorder and Schizophrenia, and Can Occur through Two Distinct Regions of the Protein.

The presence of proteinopathy, the accumulation of specific proteins as aggregates in neurons, is an emerging aspect of the pathology of schizophrenia and other major mental illnesses. Among the initial proteins implicated in forming such aggregates in these conditions is Trio and F-actin Binding Protein isoform 1 (TRIOBP-1), a ubiquitously expressed protein involved in the stabilization of the actin cytoskeleton. Here we investigate the insolubility of TRIOBP-1, as an indicator of aggregation, in brain samples from 25 schizophrenia patients, 25 major depressive disorder patients and 50 control individuals (anterior cingulate cortex, BA23). Strikingly, insoluble TRIOBP-1 is considerably more prevalent in both of these conditions than in controls, further implicating TRIOBP-1 aggregation in schizophrenia and indicating a role in major depressive disorder. These results were only seen using a high stringency insolubility assay (previously used to study DISC1 and other proteins), but not a lower stringency assay that would be expected to also detect functional, actin-bound TRIOBP-1. Previously, we have also determined that a region of 25 amino acids in the center of this protein is critical for its ability to form aggregates. Here we attempt to refine this further, through the expression of various truncated mutant TRIOBP-1 vectors in neuroblastoma cells and examining their aggregation. In this way, it was possible to narrow down the aggregation-critical region of TRIOBP-1 to just 8 amino acids (333-340 of the 652 amino acid-long TRIOBP-1). Surprisingly our results suggested that a second section of TRIOBP-1 is also capable of independently inducing aggregation: the optionally expressed 59 amino acids at the extreme N-terminus of the protein. As a result, the 597 amino acid long version of TRIOBP-1 (also referred to as "Tara" or "TAP68") has reduced potential to form aggregates. The presence of insoluble TRIOBP-1 in brain samples from patients, combined with insight into the mechanism of aggregation of TRIOBP-1 and generation of an aggregation-resistant mutant TRIOBP-1 that lacks both these regions, will be of significant use in further investigating the mechanism and consequences of TRIOBP-1 aggregation in major mental illness.

NURR1 and ERR1 Modulate the Expression of Genes of a DRD2 Coexpression Network Enriched for Schizophrenia Risk.

Multiple schizophrenia (SCZ) risk loci may be involved in gene co-regulation mechanisms, and analysis of coexpressed gene networks may help to clarify SCZ molecular basis. We have previously identified a dopamine D2 receptor (DRD2) coexpression module enriched for SCZ risk genes and associated with cognitive and neuroimaging phenotypes of SCZ, as well as with response to treatment with antipsychotics. Here we aimed to identify regulatory factors modulating this coexpression module and their relevance to SCZ. We performed motif enrichment analysis to identify transcription factor (TF) binding sites in human promoters of genes coexpressed with DRD2. Then, we measured transcript levels of a group of these genes in primary mouse cortical neurons in basal conditions and upon overexpression and knockdown of predicted TFs. Finally, we analyzed expression levels of these TFs in dorsolateral prefrontal cortex (DLPFC) of SCZ patients. Our in silico analysis revealed enrichment for NURR1 and ERR1 binding sites. In neuronal cultures, the expression of genes either relevant to SCZ risk (Drd2, Gatad2a, Slc28a1, Cnr1) or indexing coexpression in our module (Btg4, Chit1, Osr1, Gpld1) was significantly modified by gain and loss of Nurr1 and Err1. Postmortem DLPFC expression data analysis showed decreased expression levels of NURR1 and ERR1 in patients with SCZ. For NURR1 such decreased expression is associated with treatment with antipsychotics. Our results show that NURR1 and ERR1 modulate the transcription of DRD2 coexpression partners and support the hypothesis that NURR1 is involved in the response to SCZ treatment.SIGNIFICANCE STATEMENT In the present study, we provide in silico and experimental evidence for a role of the TFs NURR1 and ERR1 in modulating the expression pattern of genes coexpressed with DRD2 in human DLPFC. Notably, genetic variations in these genes is associated with SCZ risk and behavioral and neuroimaging phenotypes of the disease, as well as with response to treatment. Furthermore, this study presents novel findings on a possible interplay between D2 receptor-mediated dopamine signaling involved in treatment with antipsychotics and the transcriptional regulation mechanisms exerted by NURR1. Our results suggest that coexpression and co-regulation mechanisms may help to explain some of the complex biology of genetic associations with SCZ.

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.

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.

Developmental effects of maternal smoking during pregnancy on the human frontal cortex transcriptome.

Cigarette smoking during pregnancy is a major public health concern. While there are well-described consequences in early child development, there is very little known about the effects of maternal smoking on human cortical biology during prenatal life. We therefore performed a genome-wide differential gene expression analysis using RNA sequencing (RNA-seq) on prenatal (N = 33; 16 smoking-exposed) as well as adult (N = 207; 57 active smokers) human postmortem prefrontal cortices. Smoking exposure during the prenatal period was directly associated with differential expression of 14 genes; in contrast, during adulthood, despite a much larger sample size, only two genes showed significant differential expression (FDR < 10%). Moreover, 1,315 genes showed significantly different exposure effects between maternal smoking during pregnancy and direct exposure in adulthood (FDR < 10%)-these differences were largely driven by prenatal differences that were enriched for pathways previously implicated in addiction and synaptic function. Furthermore, prenatal and age-dependent differentially expressed genes were enriched for genes implicated in non-syndromic autism spectrum disorder (ASD) and were differentially expressed as a set between patients with ASD and controls in postmortem cortical regions. These results underscore the enhanced sensitivity to the biological effect of smoking exposure in the developing brain and offer insight into how maternal smoking during pregnancy affects gene expression in the prenatal human cortex. They also begin to address the relationship between in utero exposure to smoking and the heightened risks for the subsequent development of neuropsychiatric disorders.

qSVA framework for RNA quality correction in differential expression analysis.

RNA sequencing (RNA-seq) is a powerful approach for measuring gene expression levels in cells and tissues, but it relies on high-quality RNA. We demonstrate here that statistical adjustment using existing quality measures largely fails to remove the effects of RNA degradation when RNA quality associates with the outcome of interest. Using RNA-seq data from molecular degradation experiments of human primary tissues, we introduce a method-quality surrogate variable analysis (qSVA)-as a framework for estimating and removing the confounding effect of RNA quality in differential expression analysis. We show that this approach results in greatly improved replication rates (>3×) across two large independent postmortem human brain studies of schizophrenia and also removes potential RNA quality biases in earlier published work that compared expression levels of different brain regions and other diagnostic groups. Our approach can therefore improve the interpretation of differential expression analysis of transcriptomic data from human tissue.

Integrated DNA methylation and gene expression profiling across multiple brain regions implicate novel genes in Alzheimer's disease.

Late-onset Alzheimer's disease (AD) is a complex age-related neurodegenerative disorder that likely involves epigenetic factors. To better understand the epigenetic state associated with AD, we surveyed 420,852 DNA methylation (DNAm) sites from neurotypical controls (N = 49) and late-onset AD patients (N = 24) across four brain regions (hippocampus, entorhinal cortex, dorsolateral prefrontal cortex and cerebellum). We identified 858 sites with robust differential methylation collectively annotated to 772 possible genes (FDR < 5%, within 10 kb). These sites were overrepresented in AD genetic risk loci (p = 0.00655) and were enriched for changes during normal aging (p < 2.2 × 10-16), and nearby genes were enriched for processes related to cell-adhesion, immunity, and calcium homeostasis (FDR < 5%). To functionally validate these associations, we generated and analyzed corresponding transcriptome data to prioritize 130 genes within 10 kb of the differentially methylated sites. These 130 genes were differentially expressed between AD cases and controls and their expression was associated with nearby DNAm (p < 0.05). This integrated analysis implicates novel genes in Alzheimer's disease, such as ANKRD30B. These results highlight DNAm differences in Alzheimer's disease that have gene expression correlates, further implicating DNAm as an epigenetic mechanism underlying pathological molecular changes associated with AD. Furthermore, our framework illustrates the value of integrating epigenetic and transcriptomic data for understanding complex disease.