Individual variation in the emergence of anterior-to-posterior neural fates from human pluripotent stem cells.

Variability between human pluripotent stem cell (hPSC) lines remains a challenge and opportunity in biomedicine. In this study, hPSC lines from multiple donors were differentiated toward neuroectoderm and mesendoderm lineages. We revealed dynamic transcriptomic patterns that delineate the emergence of these lineages, which were conserved across lines, along with individual line-specific transcriptional signatures that were invariant throughout differentiation. These transcriptomic signatures predicted an antagonism between SOX21-driven forebrain fates and retinoic acid-induced hindbrain fates. Replicate lines and paired adult tissue demonstrated the stability of these line-specific transcriptomic traits. We show that this transcriptomic variation in lineage bias had both genetic and epigenetic origins, aligned with the anterior-to-posterior structure of early mammalian development, and was present across a large collection of hPSC lines. These findings contribute to developing systematic analyses of PSCs to define the origin and consequences of variation in the early events orchestrating individual human development.

Transcriptomic analysis of the human habenula in schizophrenia.

Pathophysiology of many neuropsychiatric disorders, including schizophrenia (SCZD), is linked to habenula (Hb) function. While pharmacotherapies and deep brain stimulation targeting the Hb are emerging as promising therapeutic treatments, little is known about the cell type-specific transcriptomic organization of the human Hb or how it is altered in SCZD. Here we define the molecular neuroanatomy of the human Hb and identify transcriptomic changes in individuals with SCZD compared to neurotypical controls. Utilizing Hb-enriched postmortem human brain tissue, we performed single nucleus RNA-sequencing (snRNA-seq; n=7 neurotypical donors) and identified 17 molecularly defined Hb cell types across 16,437 nuclei, including 3 medial and 7 lateral Hb populations, several of which were conserved between rodents and humans. Single molecule fluorescent in situ hybridization (smFISH; n=3 neurotypical donors) validated snRNA-seq Hb cell types and mapped their spatial locations. Bulk RNA-sequencing and cell type deconvolution in Hb-enriched tissue from 35 individuals with SCZD and 33 neurotypical controls yielded 45 SCZD-associated differentially expressed genes (DEGs, FDR < 0.05), with 32 (71%) unique to Hb-enriched tissue. eQTL analysis identified 717 independent SNP-gene pairs (FDR < 0.05), where either the SNP is a SCZD risk variant (16 pairs) or the gene is a SCZD DEG (7 pairs). eQTL and SCZD risk colocalization analysis identified 16 colocalized genes. These results identify topographically organized cell types with distinct molecular signatures in the human Hb and demonstrate unique genetic changes associated with SCZD, thereby providing novel molecular insights into the role of Hb in neuropsychiatric disorders. One Sentence Summary: Transcriptomic analysis of the human habenula and identification of molecular changes associated with schizophrenia risk and illness state.

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.