Transcriptomic analysis of the human habenula in schizophrenia.

Importance: Habenula (Hb) pathophysiology is involved in many neuropsychiatric disorders, including schizophrenia. Deep brain stimulation and pharmacological targeting of the Hb are emerging as promising therapeutic treatments. However, little is known about the cell type-specific transcriptomic organization of the human Hb or how it is altered in schizophrenia. Objective: To define the molecular neuroanatomy of the human habenula and identify transcriptomic changes in individuals with schizophrenia compared to neurotypical controls. Design Setting and Participants: This study utilized Hb-enriched postmortem human brain tissue. Single nucleus RNA-sequencing (snRNA-seq) and single molecule fluorescent in situ hybridization (smFISH) experiments were conducted to identify molecularly defined Hb cell types and map their spatial location (n=3-7 donors). Bulk RNA-sequencing and cell type deconvolution were used to investigate transcriptomic changes in Hb-enriched tissue from 35 individuals with schizophrenia and 33 neurotypical controls. Gene expression changes associated with schizophrenia in the Hb were compared to those previously identified in the dorsolateral prefrontal cortex (DLPFC), hippocampus, and caudate. Main Outcomes and Measures: Semi-supervised snRNA-seq cell type clustering. Transcript visualization and quantification of smFISH probes. Bulk RNA-seq cell type deconvolution using reference snRNA-seq data. Schizophrenia-associated gene differential expression analysis adjusting for Hb and thalamus fractions, RNA degradation-associated quality surrogate variables, and other covariates. Cross-brain region schizophrenia-associated gene expression comparison. Results: snRNA-seq identified 17 cell type clusters across 16,437 nuclei, including 3 medial and 7 lateral Hb populations. Cell types were conserved with those identified in a rodent model. smFISH for cell type marker genes validated snRNA-seq Hb cell types and depicted the spatial organization of subpopulations. Bulk RNA-seq analyses yielded 45 schizophrenia-associated differentially expressed genes (FDR < 0.05), with 32 (71%) unique to Hb-enriched tissue. Conclusions: These results identify topographically organized cell types with distinct molecular signatures in the human Hb. They further demonstrate unique transcriptomic changes in the epithalamus associated with schizophrenia, thereby providing molecular insights into the role of Hb in neuropsychiatric disorders.

Identification and evaluation of midbrain specific longevity-related genes in exceptionally long-lived but healthy mice.

Brain aging is a complex biological process that is affected by both genetic background and environment. The transcriptomic analysis of aged human and rodent brains has been applied to identify age-associated molecular and cellular processes for which intervention could possibly restore declining brain functions induced by aging. However, whether these age-associated genetic alterations are indeed involved in the healthy aging of the brain remains unclear. We herein characterized a naturally occurring, extremely long-lived (34 months of age) but healthy mouse group retaining well-preserved motor functions. Strikingly, these long-lived mice maintained tyrosine hydroxylase expression and dopaminergic fiber densities, even in the presence of persistent neuroinflammation and expression of aging markers. Combined with Endeavor gene prioritization, we identified the following midbrain-specific longevity-associated genes in the midbrain of these mice: aimp2, hexb, cacybp, akt2, nrf1, axin1, wwp2, sp2, dnajb9, notch, traf7, and lrp1. A detailed biochemical analysis of the midbrain of these long-lived mice confirmed the increased expression of Nrf1 and the activation of Akt1 and 2. Interestingly, dopaminergic neuroprotective and age-associated E3 ubiquitin ligase parkin expression was retained at high levels in the aforementioned midbrains, possibly supporting the suppression of its toxic substrates AIMP2 and PARIS. In contrast, the 24-month-old mice with dopaminergic neurite deficits failed to maintain parkin expression in the midbrain. AIMP2-induced cytotoxicity, mitochondrial stress, and neurite toxicity can be prevented by overexpression of parkin, Akt1, and Nrf1 in SH-SY5Y and PC12 cells, and basal expression of parkin, Akt1, and Nrf1 is required for maintenance of mitochondrial function and neurite integrity in PC12 cells. Taken together, this longevity-associated pathway could be a potential target of intervention to maintain nigrostriatal dopaminergic fibers and motor ability to ensure healthy longevity.

Variation of Human Neural Stem Cells Generating Organizer States In Vitro before Committing to Cortical Excitatory or Inhibitory Neuronal Fates.

Better understanding of the progression of neural stem cells (NSCs) in the developing cerebral cortex is important for modeling neurogenesis and defining the pathogenesis of neuropsychiatric disorders. Here, we use RNA sequencing, cell imaging, and lineage tracing of mouse and human in vitro NSCs and monkey brain sections to model the generation of cortical neuronal fates. We show that conserved signaling mechanisms regulate the acute transition from proliferative NSCs to committed glutamatergic excitatory neurons. As human telencephalic NSCs develop from pluripotency in vitro, they transition through organizer states that spatially pattern the cortex before generating glutamatergic precursor fates. NSCs derived from multiple human pluripotent lines vary in these early patterning states, leading differentially to dorsal or ventral telencephalic fates. This work furthers systematic analyses of the earliest patterning events that generate the major neuronal trajectories of the human telencephalon.

Comparison of quantitative trait loci methods: Total expression and allelic imbalance method in brain RNA-seq.

BACKGROUND: Of the 108 Schizophrenia (SZ) risk-loci discovered through genome-wide association studies (GWAS), 96 are not altering the sequence of any protein. Evidence linking non-coding risk-SNPs and genes may be established using expression quantitative trait loci (eQTL). However, other approaches such allelic expression quantitative trait loci (aeQTL) also may be of use. METHODS: We applied both the eQTL and aeQTL analysis to a biobank of deeply sequenced RNA from 680 dorso-lateral pre-frontal cortex (DLPFC) samples. For each of 340 genes proximal to the SZ risk-SNPs, we asked how much SNP-genotype affected total expression (eQTL), as well as how much the expression ratio between the two alleles differed from 1:1 as a consequence of the risk-SNP genotype (aeQTL). RESULTS: We analyzed overlap with comparable eQTL-findings: 16 of the 30 risk-SNPs known to have gene-level eQTL also had gene-level aeQTL effects. 6 of 21 risk-SNPs with known splice-eQTL had exon-aeQTL effects. 12 novel potential risk genes were identified with the aeQTL approach, while 55 tested SNP-pairs were found as eQTL but not aeQTL. Of the tested 108 loci we could find at least one gene to be associated with 21 of the risk-SNPs using gene-level aeQTL, and with an additional 18 risk-SNPs using exon-level aeQTL. CONCLUSION: Our results suggest that the aeQTL strategy complements the eQTL approach to susceptibility gene identification.

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

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

Genetic risk mechanisms of posttraumatic stress disorder in the human brain.

Posttraumatic stress disorder (PTSD) follows exposure to a traumatic event in susceptible individuals. Recently, genome-wide association studies have identified a number of genetic sequence variants that are associated with the risk of developing PTSD. To follow up on identifying the molecular mechanisms of these risk variants, we performed genotype to RNA sequencing-derived quantitative expression (whole gene, exon, and exon junction levels) analysis in the dorsolateral prefrontal cortex (DLPFC) of normal postmortem human brains. We further investigated genotype-gene expression associations within the amygdala in a smaller independent RNA sequencing (Genotype-Tissue Expression [GTEx]) dataset. Our DLPFC analyses identified significant expression quantitative trait loci (eQTL) associations for a "candidate" PTSD risk SNP rs363276 and the expression of two genes: SLC18A2 and PDZD8, where the PTSD risk/minor allele T was associated with significantly lower levels of gene expression for both genes, in the DLPFC. These eQTL associations were independently confirmed in the amygdala from the GTEx database. Rs363276 "T" carriers also showed significantly increased activity in the amygdala during an emotional face-matching task in healthy volunteers. Taken together, our preliminary findings in normal human brains represent a tractable approach to identify mechanisms by which genetic variants potentially increase an individual's risk for developing PTSD. © 2016 Wiley Periodicals, Inc.

Developmental regulation of human cortex transcription and its clinical relevance at single base resolution.

Transcriptome analysis of human brain provides fundamental insight into development and disease, but it largely relies on existing annotation. We sequenced transcriptomes of 72 prefrontal cortex samples across six life stages and identified 50,650 differentially expression regions (DERs) associated with developmental and aging, agnostic of annotation. While many DERs annotated to non-exonic sequence (41.1%), most were similarly regulated in cytosolic mRNA extracted from independent samples. The DERs were developmentally conserved across 16 brain regions and in the developing mouse cortex, and were expressed in diverse cell and tissue types. The DERs were further enriched for active chromatin marks and clinical risk for neurodevelopmental disorders such as schizophrenia. Lastly, we demonstrate quantitatively that these DERs associate with a changing neuronal phenotype related to differentiation and maturation. These data show conserved molecular signatures of transcriptional dynamics across brain development, have potential clinical relevance and highlight the incomplete annotation of the human brain transcriptome.

Recognition of partially occluded and rotated images with a network of spiking neurons.

In this paper, we introduce a novel system for recognition of partially occluded and rotated images. The system is based on a hierarchical network of integrate-and-fire spiking neurons with random synaptic connections and a novel organization process. The network generates integrated output sequences that are used for image classification. The proposed network is shown to provide satisfactory predictive performance given that the number of the recognition neurons and synaptic connections are adjusted to the size of the input image. Comparison of synaptic plasticity activity rule (SAPR) and spike timing dependant plasticity rules, which are used to learn connections between the spiking neurons, indicates that the former gives better results and thus the SAPR rule is used. Test results show that the proposed network performs better than a recognition system based on support vector machines.