Search

Most transcriptome-wide association studies (TWASs) so far focus on European ancestry and lack diversity. To overcome this limitation, we aggregated genome-wide association study (GWAS) summary statistics, whole-genome sequences and expression quantitative trait locus (eQTL) data from diverse ancestries. We developed a new approach, TESLA (multi-ancestry integrative study using an optimal linear combination of association statistics), to integrate an eQTL dataset with a multi-ancestry GWAS. By exploiting shared phenotypic effects between ancestries and accommodating potential effect heterogeneities, TESLA improves power over other TWAS methods. When applied to tobacco use phenotypes, TESLA identified 273 new genes, up to 55% more compared with alternative TWAS methods. These hits and subsequent fine mapping using TESLA point to target genes with biological relevance. In silico drug-repurposing analyses highlight several drugs with known efficacy, including dextromethorphan and galantamine, and new drugs such as muscle relaxants that may be repurposed for treating nicotine addiction.

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.

Ancestral differences in genomic variation are determining factors in gene regulation; however, most gene expression studies have been limited to European ancestry samples or adjusted for ancestry to identify ancestry-independent associations. We instead examined the impact of genetic ancestry on gene expression and DNA methylation (DNAm) in admixed African/Black American neurotypical individuals to untangle effects of genetic and environmental factors. Ancestry-associated differentially expressed genes (DEGs), transcripts, and gene networks, while notably not implicating neurons, are enriched for genes related to immune response and vascular tissue and explain up to 26% of heritability for ischemic stroke, 27% of heritability for Parkinson’s disease, and 30% of heritability for Alzhemier’s disease. Ancestry-associated DEGs also show general enrichment for heritability of diverse immune-related traits but depletion for psychiatric-related traits. The cell-type enrichments and direction of effects vary by brain region. These DEGs are less evolutionarily constrained and are largely explained by genetic variations; roughly 15% are predicted by DNAm variation implicating environmental exposures. We also compared Black and White Americans, confirming most of these ancestry-associated DEGs. Our results highlight how environment and genetic background affect genetic ancestry differences in gene expression in the human brain and affect risk for brain illness.SummaryWe examine the impact of genetic ancestry on gene expression and DNA methylation of admixed African/Black Americans, highlighting how genetic and environmental background affect risk for brain illness.

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.

Norepinephrine (NE) neurons in the locus coeruleus (LC) project widely throughout the central nervous system, playing critical roles in arousal and mood, as well as various components of cognition including attention, learning, and memory. The LC-NE system is also implicated in multiple neurological and neuropsychiatric disorders. Importantly, LC-NE neurons are highly sensitive to degeneration in both Alzheimer’s and Parkinson’s disease. Despite the clinical importance of the brain region and the prominent role of LC-NE neurons in a variety of brain and behavioral functions, a detailed molecular characterization of the LC is lacking. Here, we used a combination of spatially-resolved transcriptomics and single-nucleus RNA-sequencing to characterize the molecular landscape of the LC region and the transcriptomic profile of LC-NE neurons in the human brain. We provide a freely accessible resource of these data in web-accessible formats.

Schizophrenia is a complex neuropsychiatric disorder with sexually dimorphic features, including differential symptomatology, drug responsiveness, and male incidence rate. To date, only the prefrontal cortex has been examined in large-scale transcriptome analyses for sex differences in schizophrenia. Here, we examined the BrainSeq Consortium RNA-sequencing and genotypes for the caudate nucleus (n=399), dorsolateral prefrontal cortex (DLPFC; n=377), and hippocampus (n=394) to characterize sex differences in schizophrenia. We identified genomic features (genes, transcripts, exons, and exon-exon junctions) associated with sex, sex-specific expression in schizophrenia, and sex-interacting expression quantitative trait loci (si-eQTL) associated with schizophrenia risk. We found 878 unique genes with sex differences across brain regions, including ANK3, which shows male-biased expression in the caudate nucleus. X-chromosome dosage was significantly decreased in the hippocampus of female and male individuals with schizophrenia. Our sex interaction model revealed 15 novel junctions dysregulated for schizophrenia in a sex-specific manner. Sex-specific schizophrenia analysis identified dozens of expressed, sex-specific features with enrichment in the transcriptional response of cellular stress. Finally, our si-eQTL analysis revealed 974 unique genes, 14 of which are associated with schizophrenia risk. Overall, our results increased the number of annotated sex-biased features, identified sex-specific schizophrenia genes, and provided the first annotation of si-eQTL in the human DLPFC and hippocampus. Altogether, these results point to the importance of sex-informed analysis of sexually dimorphic traits and inform personalized therapeutic strategies in schizophrenia.

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.

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.

Posttraumatic stress disorder (PTSD) is a debilitating neuropsychiatric disease that is highly comorbid with major depressive disorder (MDD) and bipolar disorder. The overlap in symptoms is hypothesized to stem from partially shared genetics and underlying neurobiological mechanisms. To delineate conservation between transcriptional patterns across PTSD and MDD, the authors examined gene expression in the human cortex and amygdala in these disorders.RNA sequencing was performed in the postmortem brain of two prefrontal cortex regions and two amygdala regions from donors diagnosed with PTSD (N=107) or MDD (N=109) as well as from neurotypical donors (N=109).The authors identified a limited number of differentially expressed genes (DEGs) specific to PTSD, with nearly all mapping to cortical versus amygdala regions. PTSD-specific DEGs were enriched in gene sets associated with downregulated immune-related pathways and microglia as well as with subpopulations of GABAergic inhibitory neurons. While a greater number of DEGs associated with MDD were identified, most overlapped with PTSD, and only a few were MDD specific. The authors used weighted gene coexpression network analysis as an orthogonal approach to confirm the observed cellular and molecular associations.These findings provide supporting evidence for involvement of decreased immune signaling and neuroinflammation in MDD and PTSD pathophysiology, and extend evidence that GABAergic neurons have functional significance in PTSD.

Schizophrenia (SCZ) is a brain disorder originating during neurodevelopment with complex genetic and environmental etiologies. Despite decades of clinical evidence of altered striatal function in affected patients, its cellular and molecular underpinnings remain unclear. Here, to explore neurodevelopmental alterations in the striatum associated with SCZ, we established a method for the differentiation of iPS cells into ventral forebrain organoids. Given substantial genetic heterogeneity among individuals, which can obscure disease-associated phenotypes, we generated organoids from postmortem dural fibroblast-derived iPS cells of 3 patients and 4 healthy control individuals with nonoverlapping polygenic risk score (PRS) for SCZ and whose genotype and postmortem caudate transcriptomic data were profiled in the Brainseq neurogenomics consortium. Single cell RNA sequencing (scRNA-seq) analyses of the organoids revealed differences in developmental trajectory between SCZ cases and controls in which inhibitory neurons from patients exhibited accelerated maturation. Furthermore, we found a significant overlap of genes upregulated in the inhibitory neurons in SCZ organoids with upregulated genes in postmortem caudate tissues from patients with SCZ compared with control individuals, including the donors of our iPS cell cohort. Our findings suggest that striatal neurons in the patients with SCZ carry abnormalities that originated during early brain development and a ventral forebrain striatal organoid model can recapitulate those neurodevelopmental phenotypes in a dish.

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.

We used the 10x Genomics Visium platform to define the spatial topography of gene expression in the six-layered human dorsolateral prefrontal cortex (DLPFC). We identified extensive layer-enriched expression signatures, and refined associations to previous laminar markers. We overlaid our laminar expression signatures onto large-scale single nuclei RNA sequencing data, enhancing spatial annotation of expression-driven clusters. By integrating neuropsychiatric disorder gene sets, we showed differential layer-enriched expression of genes associated with schizophrenia and autism spectrum disorder, highlighting the clinical relevance of spatially-defined expression. We then developed a data-driven framework to define unsupervised clusters in spatial transcriptomics data, which can be applied to other tissues or brain regions where morphological architecture is not as well-defined as cortical laminae. We lastly created a web application for the scientific community to explore these raw and summarized data to augment ongoing neuroscience and spatial transcriptomics research (http://research.libd.org/spatialLIBD).

X-linked Dystonia-Parkinsonism (XDP) is an inherited, X-linked, adult-onset movement disorder characterized by degeneration in the neostriatum. No therapeutics alter disease progression. The mechanisms underlying regional differences in degeneration and age of onset are unknown. Developing therapeutics that target XDP-related mechanisms requires a deeper understanding of how XDP-relevant features vary in health and disease. XDP is due, in part, to either a partial loss of TAF1 function and/or a SVA-driven pathological gain of function. A disease-specific SINE-VNTR-Alu (SVA) retrotransposon insertion occurs within intron 32 of TAF1, a subunit of TFIID involved in transcription initiation. While all XDP males are usually clinically affected, females are heterozygous carriers generally not manifesting the full syndrome. As a resource for disease modeling, we characterized eight iPSC lines from XDP female carrier individuals, and identified isogenic lines where one clonal iPSC line expressed the wild-type X, and the two other clonal iPSC lines expressed the XDP haplotype. Furthermore, we characterized XDP-relevant transcript expression variation in humans, and found that SVA-F expression decreases slightly after 30 years of age in the neurotypical human brain and that TAF1 is modestly decreased in the majority of female samples. Uniquely in the caudate nucleus, TAF1 expression is not sexually dymorphic and decreased after 15 years of age. These findings indicate that regional-, age- and sex-specific mechanisms regulate TAF1, highlighting the importance of disease-relevant models and postmortem tissue analysis. We propose that the decreased TAF1 expression in the adult caudate may synergize with the XDP-specific partial loss of TAF1 function in patients, thereby passing a minimum threshold of TAF1 function, and triggering degeneration in the neostriatum.Significance StatementXDP is an inherited, X-linked, adult-onset movement disorder characterized by degeneration in the neostriatum. No therapeutics alter disease progression. Developing therapeutics requires a deeper understanding of how XDP-relevant features vary in health and disease. XDP is possibly due to a partial loss of TAF1 function. While all XDP males are usually affected, females are heterozygous carriers generally not manifesting the full syndrome. As a resource for disease modeling, we characterized eight stem cell lines from XDP female carrier individuals. Furthermore, we found that, uniquely in the caudate nucleus, TAF1 expression decreases after adolescence in healthy humans. We hypothesize that the decrease of TAF1 after adolescence in human caudate, in general, may underlie the vulnerability of the adult neostriatum in XDP.

OBJECTIVE: The authors sought to study the transcriptomic and genomic features of completed suicide by parsing the method chosen, to capture molecular correlates of the distinctive frame of mind of individuals who die by suicide while reducing heterogeneity. METHOD: The authors analyzed gene expression (RNA sequencing) from postmortem dorsolateral prefrontal cortex of patients who died by suicide with violent versus non-violent means and other non-suicide patients with the same psychiatric disorders, and also of neurotypicals (total n=329). They then examined genomic risk-scores (GRS) for each psychiatric disorder included, and GRS for cognition (IQ) and for suicide attempt, testing how they predict diagnosis or traits (total n=888). RESULTS: Patients who died by suicide by violent means (vS-pt) showed a transcriptomic pattern remarkably divergent from each of the other patient groups but less from neurotypicals; consistently, their genomic profile of risk was relatively low for their diagnosed illness as well as for suicide attempt and relatively high for IQ, i.e., more similar to neurotypicals than other patients. Differentially expressed genes (DEGs) associated with vS-pt pointed to purinergic signaling in microglia, showing similarities to GWAS of Drosophila aggression. WGCNA revealed that these DEGs were co-expressed in a context of mitochondrial metabolic activation unique to vS-pt. CONCLUSIONS: These findings suggest that vS-pt are in part biologically separable from other patients with the same diagnoses, and their behavioral outcome may be less dependent on genetic risk for conventional psychiatric disorders while associated with an alteration of purinergic signaling and mitochondrial metabolism.

DNA methylation (DNAm) is a key epigenetic regulator of gene expression across development. The developing prenatal brain is a highly dynamic tissue, but our understanding of key drivers of epigenetic variability across development is limited. We therefore assessed genomic methylation at over 39 million sites in the prenatal cortex using whole genome bisulfite sequencing and found loci and regions in which methylation levels are dynamic across development. We saw that DNAm at these loci was associated with nearby gene expression and enriched for enhancer chromatin states in prenatal brain tissue. Additionally, these loci were enriched for genes associated with psychiatric disorders and genes involved with neurogenesis. We also found autosomal differences in DNAm between the sexes during prenatal development, though these have less clear functional consequences. We lastly confirmed that the dynamic methylation at this critical period is specifically CpG methylation, with very low levels of CpH methylation. Our findings provide detailed insight into prenatal brain development as well as clues to the pathogenesis of psychiatric traits seen later in life.

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 withDRD2.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. Ourin silicoanalysis 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 ofDRD2coexpression partners and support the hypothesis that NURR1 is involved in the response to SCZ treatment.SIGNIFICANCE STATEMENTIn the present study, we providein silicoand experimental evidence for a role of the TFs NURR1 and ERR1 in modulating the expression pattern of genes coexpressed withDRD2in 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.

Introduction: It is reasonable to estimate that tens of millions of men in the United States played high school football. There is societal concern that participation in football confers risk for later-in-life mental health problems. The purpose of this study is to examine whether there is an association between a personal history of playing high school football and death by suicide.Methods: The subjects were obtained from the Lieber Institute for Brain Development (LIBD) brain donation program in collaboration with the Office of the Medical Examiner at Western Michigan University Homer Stryker MD School of Medicine. Donor history was documented via medical records, mental health records, and telephone interviews with the next-of-kin.Results: The sample included 198 men aged 50 or older (median = 65.0 years, interquartile range = 57–75). There were 34.8% who participated in contact sports during high school (including football), and 29.8% participated in high school football. Approximately one-third of the sample had suicide as their manner of death (34.8%). There was no statistically significant difference in the proportions of suicide as a manner of death among those men with a personal history of playing football compared to men who did not play football or who did not play sports (p = 0.070, Odds Ratio, OR = 0.537). Those who played football were significantly less likely to have a lifetime history of a suicide attempt (p = 0.012, OR = 0.352). Men with mood disorders (p < 0.001, OR = 10.712), substance use disorders (p < 0.020, OR = 2.075), and those with a history of suicide ideation (p < 0.001, OR = 8.038) or attempts (p < 0.001, OR = 40.634) were more likely to have suicide as a manner of death. Moreover, those men with a family history of suicide were more likely to have prior suicide attempts (p = 0.031, OR = 2.153) and to have completed suicide (p = 0.001, OR = 2.927).Discussion: Suicide was related to well-established risk factors such as a personal history of a mood disorder, substance abuse disorder, prior suicide ideation, suicide attempts, and a family history of suicide attempts. This study adds to a steadily growing body of evidence suggesting that playing high school football is not associated with increased risk for suicidality or suicide during adulthood.

The molecular diversity of glia in the human hippocampus and their temporal dynamics over the lifespan remain largely unknown. Here, we performed single-nucleus RNA sequencing to generate a transcriptome atlas of the human hippocampus across the postnatal lifespan. Detailed analyses of astrocytes, oligodendrocyte lineages, and microglia identified subpopulations with distinct molecular signatures and revealed their association with specific physiological functions, age-dependent changes in abundance, and disease relevance. We further characterized spatiotemporal heterogeneity of GFAP-enriched astrocyte subpopulations in the hippocampal formation using immunohistology. Leveraging glial subpopulation classifications as a reference map, we revealed the diversity of glia differentiated from human pluripotent stem cells and identified dysregulated genes and pathological processes in specific glial subpopulations in Alzheimer's disease (AD). Together, our study significantly extends our understanding of human glial diversity, population dynamics across the postnatal lifespan, and dysregulation in AD and provides a reference atlas for stem-cell-based glial differentiation.

Immature dentate granule cells (imGCs) arising from adult hippocampal neurogenesis contribute to plasticity and unique brain functions in rodents(1,2) and they are dysregulated in multiple human neurological disorders(3–5). Little is known about molecular characteristics of adult human hippocampal imGCs and even their existence is under debate(1,6–8). Here we performed single-nucleus RNA-sequencing (snRNA-seq) aided by a validated machine learning-based analytic approach to identify imGCs and quantify their abundance in the human hippocampus during infant, child, adolescent, adult, and aging stages. We identified common molecular hallmarks of human imGCs across the lifespan and discovered age-dependent transcriptional dynamics unique to human imGCs that suggest changes in cellular functionality, niche interactions, and disease relevance, in contrast to those in mice9. We also found a decreased number of imGCs with altered gene expression in Alzheimer’s disease (AD). 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 significant 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 AD.

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-cell gene expression technologies are powerful tools to study cell types in the human brain, but efforts have largely focused on cortical brain regions. We therefore created a single-nucleus RNA-sequencing resource of 70,615 high-quality nuclei to generate a molecular taxonomy of cell types across five human brain regions that serve as key nodes of the human brain reward circuitry: nucleus accumbens, amygdala, subgenual anterior cingulate cortex, hippocampus, and dorsolateral prefrontal cortex. We first identified novel subpopulations of interneurons and medium spiny neurons (MSNs) in the nucleus accumbens and further characterized robust GABAergic inhibitory cell populations in the amygdala. Joint analyses across the 107 reported cell classes revealed cell-type substructure and unique patterns of transcriptomic dynamics. We identified discrete subpopulations of D1- and D2-expressing MSNs in the nucleus accumbens to which we mapped cell-type-specific enrichment for genetic risk associated with both psychiatric disease and addiction.

Development of resistance to chemo- and immunotherapies often occurs following treatment of melanoma brain metastasis (MBM). The brain microenvironment (BME), particularly astrocytes, cooperate toward MBM progression by upregulating secreted factors, among which we found that monocyte chemoattractant protein-1 (MCP-1) and its receptors, CCR2 and CCR4, were overexpressed in MBM compared with primary lesions. Among other sources of MCP-1 in the brain, we show that melanoma cells altered astrocyte secretome and evoked MCP-1 expression and secretion, which in turn induced CCR2 expression in melanoma cells, enhancing in vitro tumorigenic properties, such as proliferation, migration, and invasion of melanoma cells. In vivo pharmacological blockade of MCP-1 or molecular knockout of CCR2/CCR4 increased the infiltration of cytotoxic CD8+ T cells and attenuated the immunosuppressive phenotype of the BME as shown by decreased infiltration of Tregs and tumor-associated macrophages/microglia in several models of intracranially injected MBM. These in vivo strategies led to decreased MBM outgrowth and prolonged the overall survival of the mice. Our findings highlight the therapeutic potential of inhibiting interactions between BME and melanoma cells for the treatment of this disease.

DNA methylation (DNAm) is an epigenetic regulator of gene expression and a hallmark of gene-environment interaction. Using whole-genome bisulfite sequencing, we have surveyed DNAm in 344 samples of human postmortem brain tissue from neurotypical subjects and individuals with schizophrenia. We identify genetic influence on local methylation levels throughout the genome, both at CpG sites and CpH sites, with 86% of SNPs and 55% of CpGs being part of methylation quantitative trait loci (meQTLs). These associations can further be clustered into regions that are differentially methylated by a given SNP, highlighting the genes and regions with which these loci are epigenetically associated. These findings can be used to better characterize schizophrenia GWAS-identified variants as epigenetic risk variants. Regions differentially methylated by schizophrenia risk-SNPs explain much of the heritability associated with risk loci, despite covering only a fraction of the genomic space. We provide a comprehensive, single base resolution view of association between genetic variation and genomic methylation, and implicate schizophrenia GWAS-associated variants as influencing the epigenetic plasticity of the brain., The authors provide a comprehensive, single base resolution view of association between genetic variation and DNA methylation in human brain. They also show that heritability attributed to schizophrenia GWAS-associated variants reflects the epigenetic plasticity of the brain.

While the role of genomic risk for schizophrenia on brain gene co-expression networks has been described, the patterns of its manifestations are varied and complex. To acquire a deeper understanding of this issue, we implemented a novel approach to network construction by manipulating the RNA-Seq expression input to “integrate” or “remove” the “modulatory” effects of genomic risk for schizophrenia. We created co-expression networks in DLPFC from the adjusted expression input and compared them in terms of gene overlap and connectivity. We used linear regression models to remove variance explained by RNA quality, cell type proportion, age, sex and genetic ancestry. We also created co-expression networks based on the genomic profile of a normative trait, height, as a “negative control”; we also applied the same analytical approach in two independent samples: LIBD Human Brain Repository (HBR) (N=78 brains, European ancestry) and Common Mind Consortium (CMC) (N=116 brains, European ancestry). In addition to direct comparisons, we explored the biological plausibility of the differential gene clusters between co-expression networks by testing them for enrichment in relevant gene ontologies and gene sets of interest (PGC2-CLOZUK GWAS significant loci genes, height GWAS significant loci genes, genes in synaptic ontologies-SynGO and genes of the “druggable genome”). We identify several key aspects of the role of genomic risk for schizophrenia in brain co-expression networks: 1) Variability of co-expression modules with “integration” or “removal” of genomic profiles of complex traits (normal or pathological); 2) Biological plausibility of gene sets represented in the differential co-expression contrasts and potential relevance for illness etiopathogenesis; 3) Non-preferential mapping of schizophrenia GWAS loci genes to network areas apparently influenced by the genomic risk score. Overall, our study supports the notion that genomic risk for schizophrenia has an extensive and non-linear effect on brain gene co-expression networks that possibly manifests as a molecular background for gene-gene, gene-environment interactions that affect various biological pathways.

Background DNA methylation (DNAm) is a critical regulator of both development and cellular identity and shows unique patterns in neurons. To better characterize maturational changes in DNAm patterns in these cells, we profile the DNAm landscape at single-base resolution across the first two decades of human neocortical development in NeuN+ neurons using whole-genome bisulfite sequencing and compare them to non-neurons (primarily glia) and prenatal homogenate cortex. Results We show that DNAm changes more dramatically during the first 5 years of postnatal life than during the entire remaining period. We further refine global patterns of increasingly divergent neuronal CpG and CpH methylation (mCpG and mCpH) into six developmental trajectories and find that in contrast to genome-wide patterns, neighboring mCpG and mCpH levels within these regions are highly correlated. We integrate paired RNA-seq data and identify putative regulation of hundreds of transcripts and their splicing events exclusively by mCpH levels, independently from mCpG levels, across this period. We finally explore the relationship between DNAm patterns and development of brain-related phenotypes and find enriched heritability for many phenotypes within identified DNAm features. Conclusions By profiling DNAm changes in NeuN-sorted neurons over the span of human cortical development, we identify novel, dynamic regions of DNAm that would be masked in homogenate DNAm data; expand on the relationship between CpG methylation, CpH methylation, and gene expression; and find enrichment particularly for neuropsychiatric diseases in genomic regions with cell type-specific, developmentally dynamic DNAm patterns.

Background Previous findings suggest that differences in brain expression of a human-specific long intergenic noncoding RNA (LINC01268; GRCh37/hg19: LOC285758) may be linked to suicide by violent methods. We sought to replicate and extend these findings in a new sample and translate the results to the behavioral level in living healthy subjects. Methods We examined RNA sequencing data in human brains to confirm the prior postmortem association of the long intergenic noncoding RNA specifically with suicide by violent means. In addition, we used a genetic variant associated with LINC01268 expression to detect association in healthy subjects with trait aggression and with in vivo prefrontal physiology related to behavioral control. Finally, we performed weighted gene coexpression network analysis and gene ontology analysis to identify biological processes associated with a LINC01268 coexpression network. Results In the replication sample, prefrontal expression of LINC01268 was again higher in suicides by violent means (n = 65) than in both nonsuicides (n = 78; p = 1.29 × 10−6) and suicides by nonviolent means (n = 46; p = 1.4 × 10−6). In the living cohort, carriers of the minor allele of a single nucleotide polymorphism associated with increased LINC01268 expression in brain scored higher on a lifetime aggression questionnaire and show diminished engagement of prefrontal cortex (Brodmann area 10) when viewing angry faces during functional magnetic resonance imaging. Weighted gene coexpression network analysis highlighted the immune response. Conclusions These results suggest that LINC01268 influences emotional regulation, aggressive behavior, and suicide by violent means; the underlying biological dynamics may include modulation of genes potentially engaged in the immune response.

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.

X-linked dystonia-parkinsonism (XDP) is an inherited, X-linked, adult-onset movement disorder characterized by degeneration in the neostriatum. No therapeutics alter disease progression. The mechanisms underlying regional differences in degeneration and adult onset are unknown. Developing therapeutics requires a deeper understanding of how XDP-relevant features vary in health and disease. XDP is possibly due, in part, to a partial loss ofTAF1function. A disease-specific SINE-VNTR-Alu (SVA) retrotransposon insertion occurs within intron 32 ofTAF1, a subunit of TFIID involved in transcription initiation. While all XDP males are usually clinically affected, females are heterozygous carriers generally not manifesting the full syndrome. As a resource for disease modeling, we characterized eight iPSC lines from three XDP female carrier individuals for X chromosome inactivation (XCI) status and identified clonal lines that express either the wild-type X or XDP haplotype. Furthermore, we characterized XDP-relevant transcript expression in neurotypical humans, and found that SVA-F expression decreases after 30 years of age in the brain and thatTAF1is decreased in most female samples. Uniquely in the caudate nucleus,TAF1expression is not sexually dimorphic and decreased after adolescence. These findings indicate that regional-specific, age-specific, and sex-specific mechanisms regulateTAF1, highlighting the importance of disease-relevant models and postmortem tissue. We propose that the decreasedTAF1expression in the adult caudate may synergize with the XDP-specific partial loss ofTAF1function in patients, thereby passing a minimum threshold ofTAF1function, and triggering degeneration in the neostriatum.

Variability between human pluripotent stem cell (hPSC) lines remains a challenge and opportunity in biomedicine. We identify differences in the spontaneous self-organization of individual hPSC lines during self-renewal that lie along a fundamental axis of in vivo development. Distinct stable and dynamic transcriptional elements were revealed by decomposition of RNA-seq data in pluripotency and early lineage emergence. Stable differences within pluripotency predicted regional bias in the dynamics of neural differentiation that were also observed in large collections of hPSC lines. Using replicate human induced PSC (hiPSC) lines and paired adult tissue, we demonstrate that cells from individual humans expressed unique transcriptional signatures that were maintained throughout life. In addition, replicate hiPSC lines from one donor showed divergent expression phenotypes driven by distinct chromatin states. These stable transcriptional states are under both genetic and epigenetic control and predict bias in subsequent forebrain and hindbrain neural trajectories. These results define mechanisms controlling transcription in pluripotency to first establish human neural diversity. Data availability: GSE164055; https://nemoanalytics.org/p?l=KimEtAL2021&g=GBX2. O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=200 SRC="FIGDIR/small/435870v1_ufig1.gif" ALT="Figure 1"> View larger version (45K): org.highwire.dtl.DTLVardef@1b5dc19org.highwire.dtl.DTLVardef@851d5forg.highwire.dtl.DTLVardef@14b703forg.highwire.dtl.DTLVardef@1b63ad0_HPS_FORMAT_FIGEXP M_FIG C_FIG

Most studies of gene expression in the brains of individuals with schizophrenia have focused on cortical regions, but subcortical nuclei such as the striatum are prominently implicated in the disease, and current antipsychotic drugs target the striatum's dense dopaminergic innervation. Here, we performed a comprehensive analysis of the genetic and transcriptional landscape of schizophrenia in the postmortem caudate nucleus of the striatum of 443 individuals (245 neurotypical individuals, 154 individuals with schizophrenia and 44 individuals with bipolar disorder), 210 from African and 233 from European ancestries. Integrating expression quantitative trait loci analysis, Mendelian randomization with the latest schizophrenia genome-wide association study, transcriptome-wide association study and differential expression analysis, we identified many genes associated with schizophrenia risk, including potentially the dopamine D2 receptor short isoform. We found that antipsychotic medication has an extensive influence on caudate gene expression. We constructed caudate nucleus gene expression networks that highlight interactions involving schizophrenia risk. These analyses provide a resource for the study of schizophrenia and insights into risk mechanisms and potential therapeutic targets.

Multispectral fluorescence imaging coupled with linear unmixing is a form of image data collection and analysis that uses multiple fluorescent dyes - each measuring a specific biological signal - that are simultaneously measured and subsequently “unmixed” to provide a read-out for each individual signal. This strategy allows for measuring multiple signals in a single data capture session - for example, 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 postmortem 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 (https://github.com/LieberInstitute/SUFI). 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 demonstrate these remote sensing algorithms’ performance 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 that is used to quantify individual RNA transcripts at single cell resolution in intact tissues and perform differential expression analysis of smFISH data, and thereby provide a one-stop solution for multispectral fluorescence image analysis and quantification. In summary, we provide a robust, automated pipeline to assist biologists with improved spectral unmixing of multispectral fluorescence images.

Identification and characterisation of novel targets for treatment is a priority in the field of psychiatry.FKBP5is 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 thatFKBP5/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 propelFKBP5/1-targeted treatment development. Here, we performed an extensive, large-scale postmortem study (n=1024) ofFKBP5/1 examining prefrontal cortex (BA9, BA11, BA24) 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, aging and genotype on corticalFKBP5/1 expression including in a cell-type specific manner. We identified consistently heightenedFKBP5/1 levels in psychopathology and with age, but not genotype, with these effects strongest in schizophrenia. Using single-nucleus RNA sequencing (snRNAseq) and targeted histology, we established that these disease- and aging-effects onFKBP5/1 expression were most pronounced in excitatory supragranular neurons. We then found that this increase inFKBP5levels likely impacts on synaptic plasticity, asFKBP5gex levels strongly and inversely correlated with dendritic mushroom spine density and brain-derived neurotrophic factor (BDNF) levels in supragranular neurons. These findings pinpoint a novel cellular and molecular mechanism that has significant potential to open a new avenue of FKBP51 drug development to treat cognitive symptoms in psychiatric disorders.

Post-traumatic stress disorder (PTSD) is a debilitating neuropsychiatric disease with a projected lifetime risk of 8.7%. PTSD is highly comorbid with depressive disorders including major depressive disorder (MDD) and bipolar disorder (BD). It is hypothesized that the overlap in symptoms stems from partially shared underlying neurobiological mechanisms. To better understand shared and unique transcriptional patterns of PTSD and MDD we performed RNA-sequencing in the postmortem brain of two prefrontal cortex (PFC) regions and two amygdala (AMY) regions, from neurotypical donors (N=109) as well as donors with diagnoses of PTSD (N=107) or MDD (N=109) across 1285 RNA-seq samples. We identified a small number of differentially expressed genes (DEGs) specific to PTSD, mostly in the cortex compared to amygdala. PTSD-specific DEGs were preferentially enriched in cortistatin-expressing cells, a subpopulation of somatostatin interneurons. These PTSD DEGs also showed strong enrichment for gene sets associated with immune-related pathways and microglia, largely driven by decreased expression of these genes in PTSD donors. While we identified a greater number of DEGs for MDD, there were only a few that were specific to MDD as they showed high overlap with PTSD DEGs. Finally, we used weighted gene co-expression network analysis (WGCNA) as an orthogonal approach to confirm the observed cellular and molecular associations. These findings highlight the sub-population of cortistatin-expressing interneurons as having potential functional significance in PTSD and provide supporting evidence for dysregulated neuroinflammation and immune signaling in MDD and PTSD pathophysiology.

Recent genetic studies have identified variants associated with bipolar disorder (BD), but it remains unclear how brain gene expression is altered in BD and how genetic risk for BD may contribute to these alterations. Here, we obtained transcriptomes from subgenual anterior cingulate cortex and amygdala samples from post-mortem brains of individuals with BD and neurotypical controls, including 511 total samples from 295 unique donors. We examined differential gene expression between cases and controls and the transcriptional effects of BD-associated genetic variants. We found two coexpressed modules that were associated with transcriptional changes in BD: one enriched for immune and inflammatory genes and the other with genes related to the postsynaptic membrane. Over 50% of BD genome-wide significant loci contained significant expression quantitative trait loci (QTL) (eQTL), and these data converged on several individual genes, including SCN2A and GRIN2A. Thus, these data implicate specific genes and pathways that may contribute to the pathology of BP.

Background and AimsThe physical, emotional, and social impacts of opioid abuse are well known; although preclinical models reveal the neurobiological pathways altered through opioid abuse, comprehensive assessments of gene expression in human brain samples are lacking. The goals of the present study were to compare gene expression in the prefrontal cortex between brain samples of individuals who died of acute opioid intoxication and group-matched controls, and to test if differential gene expression was enriched in gene sets related to opioid use.DesignCross-sectional study using human brains donated to the Lieber Institute for Brain Development. Study groups included 72 brain samples from individuals who died of acute opioid intoxication, 53 group-matched psychiatric control samples, and 28 group-matched normal control samples.SettingMaryland, USA.ParticipantsPostmortem tissue samples of the dorsolateral prefrontal cortex from 153 deceased individuals (Mage = 35.42, SD = 9.43 years; 62% male; 77% White).MeasurementsWhole transcriptome RNA sequencing was used to generate exon counts, and differential expression was tested using limma-voom. Analyses controlled for relevant sociodemographic characteristics, technical covariates, and cryptic relatedness and batch effects using quality surrogate variable analysis. Gene set enrichment analyses (GSEA) also were conducted.FindingsSixteen genes were differentially expressed (i.e., FDR-corrected p < .10) in opioid samples compared to control samples. The top differentially expressed gene, NPAS4 (FDR adjusted p = .005), was downregulated in opioid samples and has previously been implicated in cocaine use. Enrichment analyses did not provide evidence for enrichment in pathways obviously related to opioid use.ConclusionsNPAS4 is differentially expressed in the prefrontal cortex of subjects that died of an opioid overdose, providing evidence for another gene with functional relevance to opioid use and overdose.

Altering the metabotropic glutamate receptor 3 (mGluR3) by pharmacology or genetics is associated with differences in learning and memory in animals and humans. GRM3 (the gene coding for mGluR3) is also genome-wide associated with risk for schizophrenia. The neurotransmitterThe effects of genetic variation in FOLH1 on mRNA expression in human brain and NAAG levels using 7-T magnetic resonance spectroscopy (MRS) were measured. NAAG levels and FOLH1 genetic variation were correlated with measures of cognition in subjects with psychosis and unaffected subjects. Additionally, FOLH1 genetic variation was correlated with neural activity during working memory, as measured by functional MRI (fMRI).A missense mutation in FOLH1 (rs202676 G allele) was associated with increased FOLH1 mRNA in the dorsolateral prefrontal cortex of brains from unaffected subjects and schizophrenia patients. This FOLH1 variant was associated with decreased NAAG levels in unaffected subjects and patients with psychosis. NAAG levels were positively correlated with visual memory performance. Carriers of the FOLH1 variant associated with lower NAAG levels had lower IQ scores. Carriers of this FOLH1 variant had less efficient cortical activity during working memory.These data show that higher NAAG levels are associated with better cognition, suggesting that increasing NAAG levels through FOLH1/GCPII inhibition may improve cognition. Additionally, NAAG levels measured by MRS and cortical efficiency during working memory measured by fMRI have the potential to be neuroimaging biomarkers for future clinical trials.

The distribution of 5-HT1A receptors in the subnuclei of the human caudal nucleus of solitary tract and adjacent structures in the dorsal vagal complex was studied using [3H]8-OH-DPAT, a highly selective 5-HT1A receptor agonist. The highest binding of the labeled ligand was found in the dorsal motor nucleus of the vagus, followed by the medial, intermediate, and subpostremal subnuclei of the nucleus of solitary tract. Previous animal studies suggest an important role for these structures in the regulation of visceral function, particularly for the gastrointestinal and cardiovascular systems. The results of this study suggest the possibility of an analogous role for 5-HT1A receptors in the regulation of these autonomic pathways in humans as well.

Increased dopamine (DA) signaling in the striatum has been a cornerstone hypothesis about psychosis for over 50 years. Increased dopamine release results in psychotic symptoms, while D2 dopamine receptor (DRD2) antagonists are antipsychotic. Recent schizophrenia GWAS identified risk-associated common variants near the DRD2 gene, but the risk mechanism has been unclear. To gain novel insight into risk mechanisms underlying schizophrenia, we performed a comprehensive analysis of the genetic and transcriptional landscape of schizophrenia in postmortem caudate nucleus from a cohort of 444 individuals. Integrating expression quantitative trait loci (eQTL) analysis, transcriptome wide association study (TWAS), and differential expression analysis, we found many new genes associated with schizophrenia through genetic modulation of gene expression. Using a new approach based on deep neural networks, we construct caudate nucleus gene expression networks that highlight interactions involving schizophrenia risk. Interestingly, we found that genetic risk for schizophrenia is associated with decreased expression of the short isoform of DRD2, which encodes the presynaptic autoreceptor, and not with the long isoform, which encodes the postsynaptic receptor. This association suggests that decreased control of presynaptic DA release is a potential genetic mechanism of schizophrenia risk. Altogether, these analyses provide a new resource for the study of schizophrenia that can bring insight into risk mechanisms and potential novel therapeutic targets.

Opioid abuse poses significant risk to individuals in the United States and epigenetic changes are a leading potential biomarker of abuse. Current evidence, however, is mostly limited to candidate gene analysis in whole blood. To clarify the association between opioid abuse and DNA methylation, we conducted an epigenome-wide analysis (EWAS) of DNA methylation in brains of individuals who died from opioid intoxication and controls. Tissue samples were extracted from the dorsolateral prefrontal cortex of 160 deceased individuals (Mage = 35.15, SD = 9.42 years; 62% male; 78% White). The samples included 73 individuals who died of opioid intoxication, 59 group-matched psychiatric controls, and 28 group-matched normal controls. EWAS was implemented using the Illumina Infinium MethylationEPIC BeadChip; analyses adjusted for sociodemographic characteristics, negative control and ancestry principal components, cellular composition, and surrogate variables. Epigenetic age was calculated using the Horvath and Levine clocks, and gene ontology (GO) analyses were performed. No CpG sites were epigenome-wide significant after multiple testing correction, but 13 sites reached nominal significance (p < 1.0 x 10-5). There was a significant association between opioid use and Levine phenotypic age (b = 2.24, se = 1.11, p = .045). Opioid users were approximately two years phenotypically older compared to controls. GO analyses revealed enriched pathways related to cell function and neuron differentiation, but no terms survived multiple testing correction. Results inform our understanding of the neurobiology of opioid use, and future research with larger samples across stages of opioid use will elucidate the complex genomics of opioid abuse.

Single cell/nucleus technologies are powerful tools to study cell type-specific expression in the human brain, but most large-scale efforts have focused on characterizing cortical brain regions and their constituent cell types. However, additional brain regions - particularly those embedded in basal ganglia and limbic circuits - play important roles in neuropsychiatric disorders and addiction, suggesting a critical need to better understand their molecular characteristics. We therefore created a single-nucleus RNA-sequencing (snRNA-seq) resource across five human brain regions (hippocampus, HPC; dorsolateral prefrontal cortex, DLPFC; subgenual anterior cingulate cortex, sACC; nucleus accumbens, NAc; and amygdala, AMY), with emphasis on the NAc and AMY, given their involvement in reward signaling and emotional processing. We identified distinct and potentially novel neuronal subpopulations, which we validated by smFISH for various subclasses of NAc interneurons and medium spiny neurons (MSNs). We additionally benchmarked these datasets against published datasets for corresponding regions in rodent models to define cross-species convergence and divergence across analogous cell subclasses. We characterized the transcriptomic architecture of regionally-defined neuronal subpopulations, which revealed strong patterns of similarities in specific neuronal subclasses across the five profiled regions. Finally, we measured genetic associations between risk for psychiatric disease and substance use behaviors with each of the regionally-defined cell types. This analysis further supported NAc and AMY involvement in risk for psychiatric illness by implicating specific neuronal subpopulations, and highlighted potential involvement of an MSN population associated with stress signaling in genetic risk for substance use.

DNA methylation (DNAm) regulates gene expression and may represent gene-environment interactions. Using whole genome bisulfite sequencing, we surveyed DNAm in a large sample (n=344) of human brain tissues. We identify widespread genetic influence on local methylation levels throughout the genome, with 76% of SNPs and 38% of CpGs being part of methylation quantitative trait loci (meQTLs). These associations can further be clustered into regions that are differentially methylated by a given SNP, highlighting putative functional regions that explain much of the heritability associated with risk loci. Furthermore, some CpH sites associated with genetic variation. We have established a comprehensive, single base resolution view of association between genetic variation and genomic methylation, and implicate schizophrenia GWAS-associated variants as influencing the epigenetic plasticity of the brain.One-sentence summaryMost genetic variants associated with DNA methylation levels, and implicated schizophrenia GWAS variants in the human brain.