Cellular Diversity in Human Subgenual Anterior Cingulate and Dorsolateral Prefrontal Cortex by Single-Nucleus RNA-Sequencing.

Regional cellular heterogeneity is a fundamental feature of the human neocortex; however, details of this heterogeneity are still undefined. We used single-nucleus RNA-sequencing to examine cell-specific transcriptional features in the dorsolateral PFC (DLPFC) and the subgenual anterior cingulate cortex (sgACC), regions implicated in major psychiatric disorders. Droplet-based nuclei-capture and library preparation were performed on replicate samples from 8 male donors without history of psychiatric or neurologic disorder. Unsupervised clustering identified major neural cell classes. Subsequent iterative clustering of neurons further revealed 20 excitatory and 22 inhibitory subclasses. Inhibitory cells were consistently more abundant in the sgACC and excitatory neuron subclusters exhibited considerable variability across brain regions. Excitatory cell subclasses also exhibited greater within-class transcriptional differences between the two regions. We used these molecular definitions to determine which cell classes might be enriched in loci carrying a genetic signal in genome-wide association studies or for differentially expressed genes in mental illness. We found that the heritable signals of psychiatric disorders were enriched in neurons and that, while the gene expression changes detected in bulk-RNA-sequencing studies were dominated by glial cells, some alterations could be identified in specific classes of excitatory and inhibitory neurons. Intriguingly, only two excitatory cell classes exhibited concomitant region-specific enrichment for both genome-wide association study loci and transcriptional dysregulation. In sum, by detailing the molecular and cellular diversity of the DLPFC and sgACC, we were able to generate hypotheses on regional and cell-specific dysfunctions that may contribute to the development of mental illness.SIGNIFICANCE STATEMENT Dysfunction of the subgenual anterior cingulate cortex has been implicated in mood disorders, particularly major depressive disorder, and the dorsolateral PFC, a subsection of the PFC involved in executive functioning, has been implicated in schizophrenia. Understanding the cellular composition of these regions is critical to elucidating the neurobiology underlying psychiatric and neurologic disorders. We studied cell type diversity of the subgenual anterior cingulate cortex and dorsolateral PFC of humans with no neuropsychiatric illness using a clustering analysis of single-nuclei RNA-sequencing data. Defining the transcriptomic profile of cellular subpopulations in these cortical regions is a first step to demystifying the cellular and molecular pathways involved in psychiatric disorders.

3D ATAC-PALM: super-resolution imaging of the accessible genome.

To image the accessible genome at nanometer scale in situ, we developed three-dimensional assay for transposase-accessible chromatin-photoactivated localization microscopy (3D ATAC-PALM) that integrates an assay for transposase-accessible chromatin with visualization, PALM super-resolution imaging and lattice light-sheet microscopy. Multiplexed with oligopaint DNA-fluorescence in situ hybridization (FISH), RNA-FISH and protein fluorescence, 3D ATAC-PALM connected microscopy and genomic data, revealing spatially segregated accessible chromatin domains (ACDs) that enclose active chromatin and transcribed genes. Using these methods to analyze genetically perturbed cells, we demonstrated that genome architectural protein CTCF prevents excessive clustering of accessible chromatin and decompacts ACDs. These results highlight 3D ATAC-PALM as a useful tool to probe the structure and organizing mechanism of the genome.

The relationship between case-control differential gene expression from brain tissue and genetic associations in schizophrenia.

Large numbers of genetic loci have been identified that are known to contain common risk alleles for schizophrenia, but linking associated alleles to specific risk genes remains challenging. Given that most alleles that influence liability to schizophrenia are thought to do so by altered gene expression, intuitively, case-control differential gene expression studies should highlight genes with a higher probability of being associated with schizophrenia and could help identify the most likely causal genes within associated loci. Here, we test this hypothesis by comparing transcriptome analysis of the dorsolateral prefrontal cortex from 563 schizophrenia cases and 802 controls with genome-wide association study (GWAS) data from the third wave study of the Psychiatric Genomics Consortium. Genes differentially expressed in schizophrenia were not enriched for common allelic association statistics compared with other brain-expressed genes, nor were they enriched for genes within associated loci previously reported to be prioritized by genetic fine-mapping. Genes prioritized by Summary-based Mendelian Randomization were underexpressed in cases compared to other genes in the same GWAS loci. However, the overall strength and direction of expression change predicted by SMR were not related to that observed in the differential expression data. Overall, this study does not support the hypothesis that genes identified as differentially expressed from RNA sequencing of bulk brain tissue are enriched for those that show evidence for genetic associations. Such data have limited utility for prioritizing genes in currently associated loci in schizophrenia.

GRK3 deficiency elicits brain immune activation and psychosis.

The G protein-coupled receptor kinase (GRK) family member protein GRK3 has been linked to the pathophysiology of schizophrenia and bipolar disorder. Expression, as well as protein levels, of GRK3 are reduced in post-mortem prefrontal cortex of schizophrenia subjects. Here, we investigate functional behavior and neurotransmission related to immune activation and psychosis using mice lacking functional Grk3 and utilizing a variety of methods, including behavioral, biochemical, electrophysiological, molecular, and imaging methods. Compared to wildtype controls, the Grk3-/- mice show a number of aberrations linked to psychosis, including elevated brain levels of IL-1ß, increased turnover of kynurenic acid (KYNA), hyper-responsiveness to D-amphetamine, elevated spontaneous firing of midbrain dopamine neurons, and disruption in prepulse inhibition. Analyzing human genetic data, we observe a link between psychotic features in bipolar disorder, decreased GRK expression, and increased concentration of CSF KYNA. Taken together, our data suggest that Grk3-/- mice show face and construct validity relating to the psychosis phenotype with glial activation and would be suitable for translational studies of novel immunomodulatory agents in psychotic disorders.

Melanin Concentrating Hormone Signaling Deficits in Schizophrenia: Association With Memory and Social Impairments and Abnormal Sensorimotor Gating.

BACKGROUND: Evidence from anatomical, pharmacological, and genetic studies supports a role for the neuropeptide melanin concentrating hormone system in modulating emotional and cognitive functions. Genome-wide association studies revealed a potential association between the melanin concentrating hormone receptor (MCHR1) gene locus and schizophrenia, and the largest genome-wide association study conducted to date shows a credible genome-wide association. METHODS: We analyzed MCHR1 and pro-melanin concentrating hormone RNA-Seq expression in the prefrontal cortex in schizophrenia patients and healthy controls. Disruptions in the melanin concentrating hormone system were modeled in the mouse brain by germline deletion of MCHR1 and by conditional ablation of melanin concentrating hormone expressing neurons using a Cre-inducible diphtheria toxin system. RESULTS: MCHR1 expression is decreased in the prefrontal cortex of schizophrenia samples (false discovery rate (FDR) P < .05, CommonMind and PsychEncode combined datasets, n = 901) while pro-melanin concentrating hormone is below the detection threshold. MCHR1 expression decreased with aging (P = 6.6E-57) in human dorsolateral prefrontal cortex. The deletion of MCHR1 was found to lead to behavioral abnormalities mimicking schizophrenia-like phenotypes: hyperactivity, increased stereotypic and repetitive behavior, social impairment, impaired sensorimotor gating, and disrupted cognitive functions. Conditional ablation of pro-melanin concentrating hormone neurons increased repetitive behavior and produced a deficit in sensorimotor gating. CONCLUSIONS: Our study indicates that early disruption of the melanin concentrating hormone system interferes with neurodevelopmental processes, which may contribute to the pathogenesis of schizophrenia. Further neurobiological research on the developmental timing and circuits that are affected by melanin concentrating hormone may lead to a therapeutic target for early prevention of schizophrenia.

Antipsychotic drug use complicates assessment of gene expression changes associated with schizophrenia.

Recent postmortem transcriptomic studies of schizophrenia (SCZ) have shown hundreds of differentially expressed genes. However, the extent to which these gene expression changes reflect antipsychotic drug (APD) exposure remains uncertain. We compared differential gene expression in the prefrontal cortex of SCZ patients who tested positive for APDs at the time of death with SCZ patients who did not. APD exposure was associated with numerous changes in the brain transcriptome, especially among SCZ patients on atypical APDs. Brain transcriptome data from macaques chronically treated with APDs showed that APDs affect the expression of many functionally relevant genes, some of which show expression changes in the same directions as those observed in SCZ. Co-expression modules enriched for synaptic function showed convergent patterns between SCZ and some of the APD effects, while those associated with inflammation and glucose metabolism exhibited predominantly divergent patterns between SCZ and APD effects. In contrast, major cell-type shifts inferred in SCZ were primarily unaffected by APD use. These results show that APDs may confound SCZ-associated gene expression changes in postmortem brain tissue. Disentangling these effects will help identify causal genes and improve our neurobiological understanding of SCZ.

Polygenic scores for psychiatric disorders in a diverse postmortem brain tissue cohort.

A new era of human postmortem tissue research has emerged thanks to the development of 'omics technologies that measure genes, proteins, and spatial parameters in unprecedented detail. Also newly possible is the ability to construct polygenic scores, individual-level metrics of genetic risk (also known as polygenic risk scores/PRS), based on genome-wide association studies, GWAS. Here, we report on clinical, educational, and brain gene expression correlates of polygenic scores in ancestrally diverse samples from the Human Brain Collection Core (HBCC). Genotypes from 1418 donors were subjected to quality control filters, imputed, and used to construct polygenic scores. Polygenic scores for schizophrenia predicted schizophrenia status in donors of European ancestry (p = 4.7 × 10-8, 17.2%) and in donors with African ancestry (p = 1.6 × 10-5, 10.4% of phenotypic variance explained). This pattern of higher variance explained among European ancestry samples was also observed for other psychiatric disorders (depression, bipolar disorder, substance use disorders, anxiety disorders) and for height, body mass index, and years of education. For a subset of 223 samples, gene expression from dorsolateral prefrontal cortex (DLPFC) was available through the CommonMind Consortium. In this subgroup, schizophrenia polygenic scores also predicted an aggregate gene expression score for schizophrenia (European ancestry: p = 0.0032, African ancestry: p = 0.15). Overall, polygenic scores performed as expected in ancestrally diverse samples, given historical biases toward use of European ancestry samples and variable predictive power of polygenic scores across phenotypes. The transcriptomic results reported here suggest that inherited schizophrenia genetic risk influences gene expression, even in adulthood. For future research, these and additional polygenic scores are being made available for analyses, and for selecting samples, using postmortem tissue from the Human Brain Collection Core.

Altered transcriptomes, cell type proportions, and dendritic spine morphology in hippocampus of suicide deaths.

Suicide is a condition resulting from complex environmental and genetic risks that affect millions of people globally. Both structural and functional studies identified the hippocampus as one of the vulnerable brain regions contributing to suicide risk. Here, we have identified the hippocampal transcriptomes, gene ontology, cell type proportions, dendritic spine morphology, and transcriptomic signature in iPSC-derived neuronal precursor cells (NPCs) and neurons in postmortem brain tissue from suicide deaths. The hippocampal tissue transcriptomic data revealed that NPAS4 gene expression was downregulated while ALDH1A2, NAAA, and MLXIPL gene expressions were upregulated in tissue from suicide deaths. The gene ontology identified 29 significant pathways including NPAS4-associated gene ontology terms "excitatory post-synaptic potential", "regulation of postsynaptic membrane potential" and "long-term memory" indicating alteration of glutamatergic synapses in the hippocampus of suicide deaths. The cell type deconvolution identified decreased excitatory neuron proportion and an increased inhibitory neuron proportion providing evidence of excitation/inhibition imbalance in the hippocampus of suicide deaths. In addition, suicide deaths had increased dendric spine density, due to an increase of thin (relatively unstable) dendritic spines, compared to controls. The transcriptomes of iPSC-derived hippocampal-like NPCs and neurons revealed 31 and 33 differentially expressed genes in NPC and neurons, respectively, of suicide deaths. The suicide-associated differentially expressed genes in NPCs were RELN, CRH, EMX2, OXTR, PARM1 and IFITM2 which overlapped with previously published results. The previously-known suicide-associated differentially expressed genes in differentiated neurons were COL1A1, THBS1, IFITM2, AQP1, and NLRP2. Together, these findings would help better understand the hippocampal neurobiology of suicide for identifying therapeutic targets to prevent suicide.

Functional impairment of cortical AMPA receptors in schizophrenia.

Clinical and preclinical studies suggest that some of the behavioral alterations observed in schizophrenia (SZ) may be mechanistically linked to synaptic dysfunction of glutamatergic signaling. Recent genetic and proteomic studies suggest alterations of cortical glutamate receptors of the AMPA-type (AMPARs), which are the predominant ligand-gated ionic channels of fast transmission at excitatory synapses. The impact of gene and protein alterations on the electrophysiological activity of AMPARs is not known in SZ. In this proof of principle work, using human postmortem brain synaptic membranes isolated from the dorsolateral prefrontal cortex (DLPFC), we combined electrophysiological analysis from microtransplanted synaptic membranes (MSM) with transcriptomic (RNA-Seq) and label-free proteomics data in 10 control and 10 subjects diagnosed with SZ. We observed in SZ a reduction in the amplitude of AMPARs currents elicited by kainate, an agonist of AMPARs that blocks the desensitization of the receptor. This reduction was not associated with protein abundance but with a reduction in kainate's potency to activate AMPARs. Electrophysiologically-anchored dataset analysis (EDA) was used to identify synaptosomal proteins that linearly correlate with the amplitude of the AMPARs responses, gene ontology functional annotations were then used to determine protein-protein interactions. Protein modules associated with positive AMPARs current increases were downregulated in SZ, while protein modules that were upregulated in SZ were associated with decreased AMPARs currents. Our results indicate that transcriptomic and proteomic alterations, frequently observed in the DLPFC in SZ, converge at the synaptic level producing a functional electrophysiological impairment of AMPARs.

Increased excitatory to inhibitory synaptic ratio in parietal cortex samples from individuals with Alzheimer's disease.

Synaptic disturbances in excitatory to inhibitory (E/I) balance in forebrain circuits are thought to contribute to the progression of Alzheimer's disease (AD) and dementia, although direct evidence for such imbalance in humans is lacking. We assessed anatomical and electrophysiological synaptic E/I ratios in post-mortem parietal cortex samples from middle-aged individuals with AD (early-onset) or Down syndrome (DS) by fluorescence deconvolution tomography and microtransplantation of synaptic membranes. Both approaches revealed significantly elevated E/I ratios for AD, but not DS, versus controls. Gene expression studies in an independent AD cohort also demonstrated elevated E/I ratios in individuals with AD as compared to controls. These findings provide evidence of a marked pro-excitatory perturbation of synaptic E/I balance in AD parietal cortex, a region within the default mode network that is overly active in the disorder, and support the hypothesis that E/I imbalances disrupt cognition-related shifts in cortical activity which contribute to the intellectual decline in AD.

Patterns of cilia gene dysregulations in major psychiatric disorders.

Primary cilia function as cells' antennas to detect and transduce external stimuli and play crucial roles in cell signaling and communication. The vast majority of cilia genes that are causally linked with ciliopathies are also associated with neurological deficits, such as cognitive impairments. Yet, the roles of cilia dysfunctions in the pathogenesis of psychiatric disorders have not been studied. Our aim is to identify patterns of cilia gene dysregulation in the four major psychiatric disorders: schizophrenia (SCZ), autism spectrum disorder (ASD), bipolar disorder (BP), and major depressive disorder (MDD). For this purpose, we acquired differentially expressed genes (DEGs) from the largest and most recent publicly available databases. We found that 42%, 24%, 17%, and 15% of brain-expressed cilia genes were significantly differentially expressed in SCZ, ASD, BP, and MDD, respectively. Several genes exhibited cross-disorder overlap, suggesting that typical cilia signaling pathways' dysfunctions determine susceptibility to more than one psychiatric disorder or may partially underlie their pathophysiology. Our study revealed that genes encoding proteins of almost all sub-cilia structural and functional compartments were dysregulated in the four psychiatric disorders. Strikingly, the genes of 75% of cilia GPCRs and 50% of the transition zone proteins were differentially expressed in SCZ. The present study is the first to draw associations between cilia and major psychiatric disorders, and is the first step toward understanding the role that cilia components play in their pathophysiological processes, which may lead to novel therapeutic targets for these disorders.

Deep transcriptome sequencing of subgenual anterior cingulate cortex reveals cross-diagnostic and diagnosis-specific RNA expression changes in major psychiatric disorders.

Despite strong evidence of heritability and growing discovery of genetic markers for major mental illness, little is known about how gene expression in the brain differs across psychiatric diagnoses, or how known genetic risk factors shape these differences. Here we investigate expressed genes and gene transcripts in postmortem subgenual anterior cingulate cortex (sgACC), a key component of limbic circuits linked to mental illness. RNA obtained postmortem from 200 donors diagnosed with bipolar disorder, schizophrenia, major depression, or no psychiatric disorder was deeply sequenced to quantify expression of over 85,000 gene transcripts, many of which were rare. Case-control comparisons detected modest expression differences that were correlated across disorders. Case-case comparisons revealed greater expression differences, with some transcripts showing opposing patterns of expression between diagnostic groups, relative to controls. The ~250 rare transcripts that were differentially-expressed in one or more disorder groups were enriched for genes involved in synapse formation, cell junctions, and heterotrimeric G-protein complexes. Common genetic variants were associated with transcript expression (eQTL) or relative abundance of alternatively spliced transcripts (sQTL). Common genetic variants previously associated with disease risk were especially enriched for sQTLs, which together accounted for disproportionate fractions of diagnosis-specific heritability. Genetic risk factors that shape the brain transcriptome may contribute to diagnostic differences between broad classes of mental illness.

Unexpected Transcriptional Programs Contribute to Hippocampal Memory Deficits and Neuronal Stunting after Early-Life Adversity.

Early-life adversity (ELA) is associated with lifelong memory deficits, yet the responsible mechanisms remain unclear. We impose ELA by rearing rat pups in simulated poverty, assess hippocampal memory, and probe changes in gene expression, their transcriptional regulation, and the consequent changes in hippocampal neuronal structure. ELA rats have poor hippocampal memory and stunted hippocampal pyramidal neurons associated with ~140 differentially expressed genes. Upstream regulators of the altered genes include glucocorticoid receptor and, unexpectedly, the transcription factor neuron-restrictive silencer factor (NRSF/REST). NRSF contributes critically to the memory deficits because blocking its function transiently following ELA rescues spatial memory and restores the dendritic arborization of hippocampal pyramidal neurons in ELA rats. Blocking NRSF function in vitro augments dendritic complexity of developing hippocampal neurons, suggesting that NRSF represses genes involved in neuronal maturation. These findings establish important, surprising contributions of NRSF to ELA-induced transcriptional programming that disrupts hippocampal maturation and memory function.

Cortical RORß is required for layer 4 transcriptional identity and barrel integrity.

Retinoic acid-related orphan receptor beta (RORß) is a transcription factor (TF) and marker of layer 4 (L4) neurons, which are distinctive both in transcriptional identity and the ability to form aggregates such as barrels in rodent somatosensory cortex. However, the relationship between transcriptional identity and L4 cytoarchitecture is largely unknown. We find RORß is required in the cortex for L4 aggregation into barrels and thalamocortical afferent (TCA) segregation. Interestingly, barrel organization also degrades with age in wildtype mice. Loss of RORß delays excitatory input and disrupts gene expression and chromatin accessibility, with down-regulation of L4 and up-regulation of L5 genes, suggesting a disruption in cellular specification. Expression and binding site accessibility change for many other TFs, including closure of neurodevelopmental TF binding sites and increased expression and binding capacity of activity-regulated TFs. Lastly, a putative target of RORß, Thsd7a, is down-regulated without RORß, and Thsd7a knock-out alone disrupts TCA organization in adult barrels.

Novel Complex Interactions between Mitochondrial and Nuclear DNA in Schizophrenia and Bipolar Disorder.

Mitochondrial dysfunction has been associated with schizophrenia (SZ) and bipolar disorder (BD). This review examines recent publications and novel associations between mitochondrial genes and SZ and BD. Associations of nuclear-encoded mitochondrial variants with SZ were found using gene- and pathway-based approaches. Two control region mitochondrial DNA (mtDNA) SNPs, T16519C and T195C, both showed an association with SZ and BD. A review of 4 studies of A15218G located in the cytochrome B oxidase gene (CYTB, SZ = 11,311, control = 35,735) shows a moderate association with SZ (p = 2.15E-03). Another mtDNA allele A12308G was nominally associated with psychosis in BD type I subjects and SZ. The first published study testing the epistatic interaction between nuclear-encoded and mitochondria-encoded genes demonstrated evidence for potential interactions between mtDNA and the nuclear genome for BD. A similar analysis for the risk of SZ revealed significant joint effects (34 nuclear-mitochondria SNP pairs with joint effect p ≤ 5E-07) and significant enrichment of projection neurons. The mitochondria-encoded gene CYTB was found in both the epistatic interactions for SZ and BD and the single SNP association of SZ. Future efforts considering population stratification and polygenic risk scores will test the role of mitochondrial variants in psychiatric disorders.

A repeated molecular architecture across thalamic pathways.

The thalamus is the central communication hub of the forebrain and provides the cerebral cortex with inputs from sensory organs, subcortical systems and the cortex itself. Multiple thalamic regions send convergent information to each cortical region, but the organizational logic of thalamic projections has remained elusive. Through comprehensive transcriptional analyses of retrogradely labeled thalamic neurons in adult mice, we identify three major profiles of thalamic pathways. These profiles exist along a continuum that is repeated across all major projection systems, such as those for vision, motor control and cognition. The largest component of gene expression variation in the mouse thalamus is topographically organized, with features conserved in humans. Transcriptional differences between these thalamic neuronal identities are tied to cellular features that are critical for function, such as axonal morphology and membrane properties. Molecular profiling therefore reveals covariation in the properties of thalamic pathways serving all major input modalities and output targets, thus establishing a molecular framework for understanding the thalamus.

Publisher Correction: A repeated molecular architecture across thalamic pathways.

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

Mapping the transcriptional diversity of genetically and anatomically defined cell populations in the mouse brain.

Understanding the principles governing neuronal diversity is a fundamental goal for neuroscience. Here, we provide an anatomical and transcriptomic database of nearly 200 genetically identified cell populations. By separately analyzing the robustness and pattern of expression differences across these cell populations, we identify two gene classes contributing distinctly to neuronal diversity. Short homeobox transcription factors distinguish neuronal populations combinatorially, and exhibit extremely low transcriptional noise, enabling highly robust expression differences. Long neuronal effector genes, such as channels and cell adhesion molecules, contribute disproportionately to neuronal diversity, based on their patterns rather than robustness of expression differences. By linking transcriptional identity to genetic strains and anatomical atlases, we provide an extensive resource for further investigation of mouse neuronal cell types.

Inference of cell type content from human brain transcriptomic datasets illuminates the effects of age, manner of death, dissection, and psychiatric diagnosis.

Psychiatric illness is unlikely to arise from pathology occurring uniformly across all cell types in affected brain regions. Despite this, transcriptomic analyses of the human brain have typically been conducted using macro-dissected tissue due to the difficulty of performing single-cell type analyses with donated post-mortem brains. To address this issue statistically, we compiled a database of several thousand transcripts that were specifically-enriched in one of 10 primary cortical cell types in previous publications. Using this database, we predicted the relative cell type content for 833 human cortical samples using microarray or RNA-Seq data from the Pritzker Consortium (GSE92538) or publicly-available databases (GSE53987, GSE21935, GSE21138, CommonMind Consortium). These predictions were generated by averaging normalized expression levels across transcripts specific to each cell type using our R-package BrainInABlender (validated and publicly-released on github). Using this method, we found that the principal components of variation in the datasets strongly correlated with the predicted neuronal/glial content of the samples. This variability was not simply due to dissection-the relative balance of brain cell types appeared to be influenced by a variety of demographic, pre- and post-mortem variables. Prolonged hypoxia around the time of death predicted increased astrocytic and endothelial gene expression, illustrating vascular upregulation. Aging was associated with decreased neuronal gene expression. Red blood cell gene expression was reduced in individuals who died following systemic blood loss. Subjects with Major Depressive Disorder had decreased astrocytic gene expression, mirroring previous morphometric observations. Subjects with Schizophrenia had reduced red blood cell gene expression, resembling the hypofrontality detected in fMRI experiments. Finally, in datasets containing samples with especially variable cell content, we found that controlling for predicted sample cell content while evaluating differential expression improved the detection of previously-identified psychiatric effects. We conclude that accounting for cell type can greatly improve the interpretability of transcriptomic data.