ABSTRACT
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.
ABSTRACT
Human induced pluripotent stem cells (hiPSCs) are a powerful model of neural differentiation and maturation. We present a hiPSC transcriptomics resource on corticogenesis from 5 iPSC donor and 13 subclonal lines across 9 time points over 5 broad conditions: self-renewal, early neuronal differentiation, neural precursor cells (NPCs), assembled rosettes, and differentiated neuronal cells. We identify widespread changes in the expression of both individual features and global patterns of transcription. We next demonstrate that co-culturing human NPCs with rodent astrocytes results in mutually synergistic maturation, and that cell type-specific expression data can be extracted using only sequencing read alignments without cell sorting. We lastly adapt a previously generated RNA deconvolution approach to single-cell expression data to estimate the relative neuronal maturity of iPSC-derived neuronal cultures and human brain tissue. Using many public datasets, we demonstrate neuronal cultures are maturationally heterogeneous but contain subsets of neurons more mature than previously observed.
Subject(s)
Cell Differentiation/genetics , Induced Pluripotent Stem Cells/cytology , Induced Pluripotent Stem Cells/physiology , Neural Stem Cells/physiology , Transcriptome , Algorithms , Animals , Astrocytes/cytology , Cells, Cultured , Cerebral Cortex/cytology , Coculture Techniques , Databases, Genetic , Gene Expression Regulation , Humans , Models, Neurological , Neural Stem Cells/cytology , Neurons/cytology , Neurons/physiology , RatsABSTRACT
Transcriptome compartmentalization by the nuclear membrane provides both stochastic and functional buffering of transcript activity in the cytoplasm, and has recently been implicated in neurodegenerative disease processes. Although many mechanisms regulating transcript compartmentalization are also prevalent in brain development, the extent to which subcellular localization differs as the brain matures has yet to be addressed. To characterize the nuclear and cytoplasmic transcriptomes during brain development, we sequenced both RNA fractions from homogenate prenatal and adult human postmortem cortex using poly(A)+ and Ribo-Zero library preparation methods. We find that while many genes are differentially expressed by fraction and developmental expression changes are similarly detectable in nuclear and cytoplasmic RNA, the compartmented transcriptomes become more distinct as the brain matures, perhaps reflecting increased utilization of nuclear retention as a regulatory strategy in adult brain. We examined potential mechanisms of this developmental divergence including alternative splicing, RNA editing, nuclear pore composition, RNA-binding protein motif enrichment, and RNA secondary structure. Intron retention is associated with greater nuclear abundance in a subset of transcripts, as is enrichment for several splicing factor binding motifs. Finally, we examined disease association with fraction-regulated gene sets and found nuclear-enriched genes were also preferentially enriched in gene sets associated with neurodevelopmental psychiatric disorders. These results suggest that although gene-level expression is globally comparable between fractions, nuclear retention of transcripts may play an underappreciated role in developmental regulation of gene expression in brain, particularly in genes whose dysregulation is related to neuropsychiatric disorders.
Subject(s)
Cell Nucleus/metabolism , Cerebral Cortex/metabolism , Cytoplasm/metabolism , Genetic Predisposition to Disease , Mental Disorders/genetics , Mental Disorders/psychology , Transcriptome , Age Factors , Alternative Splicing , Computational Biology/methods , Gene Expression Profiling , Gene Expression Regulation , Genetic Association Studies , High-Throughput Nucleotide Sequencing , Humans , Molecular Sequence Annotation , RNA EditingABSTRACT
BACKGROUND: RNA sequencing offers advantages over other quantification methods for microRNA (miRNA), yet numerous biases make reliable quantification challenging. Previous evaluations of these biases have focused on adapter ligation bias with limited evaluation of reverse transcription bias or amplification bias. Furthermore, evaluations of the quantification of isomiRs (miRNA isoforms) or the influence of starting amount on performance have been very limited. No study had yet evaluated the quantification of isomiRs of altered length or compared the consistency of results derived from multiple moderate starting inputs. We therefore evaluated quantifications of miRNA and isomiRs using four library preparation kits, with various starting amounts, as well as quantifications following removal of duplicate reads using unique molecular identifiers (UMIs) to mitigate reverse transcription and amplification biases. RESULTS: All methods resulted in false isomiR detection; however, the adapter-free method tested was especially prone to false isomiR detection. We demonstrate that using UMIs improves accuracy and we provide a guide for input amounts to improve consistency. CONCLUSIONS: Our data show differences and limitations of current methods, thus raising concerns about the validity of quantification of miRNA and isomiRs across studies. We advocate for the use of UMIs to improve accuracy and reliability of miRNA quantifications.
Subject(s)
Sequence Analysis, RNA/standards , Animals , Bias , Humans , Mice , RNA Isoforms , RNA, Viral , Rats , Reproducibility of Results , Sequence Analysis, RNA/methodsABSTRACT
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.
Subject(s)
Frontal Lobe , Gene Expression Regulation, Developmental/physiology , Hippocampus , Schizophrenia/genetics , Schizophrenia/pathology , Adolescent , Adult , Aged , Aged, 80 and over , Female , Frontal Lobe/embryology , Frontal Lobe/growth & development , Frontal Lobe/metabolism , Gene Ontology , Genetic Predisposition to Disease , Genome-Wide Association Study , Hippocampus/embryology , Hippocampus/growth & development , Hippocampus/metabolism , Humans , Male , Middle Aged , Young AdultABSTRACT
Brain-derived neurotrophic factor (Bdnf) transcription is controlled by several promoters, which drive expression of multiple transcripts encoding an identical protein. We previously reported that BDNF derived from promoters I and II is highly expressed in hypothalamus and is critical for regulating aggression in male mice. Here we report that BDNF loss from these promoters causes reduced sexual receptivity and impaired maternal care in female mice, which is concomitant with decreased oxytocin (Oxt) expression during development. We identify a novel link between BDNF signaling, oxytocin, and maternal behavior by demonstrating that ablation of TrkB selectively in OXT neurons partially recapitulates maternal care impairments observed in BDNF-deficient females. Using translating ribosome affinity purification and RNA-sequencing we define a molecular profile for OXT neurons and delineate how BDNF signaling impacts gene pathways critical for structural and functional plasticity. Our findings highlight BDNF as a modulator of sexually-dimorphic hypothalamic circuits that govern female-typical behaviors.
Subject(s)
Brain-Derived Neurotrophic Factor/metabolism , Neurons/metabolism , Oxytocin/metabolism , Receptor, trkB/metabolism , Signal Transduction , Animals , Estrous Cycle , Female , Gene Expression Regulation , Maternal Behavior , Mice, Inbred C57BL , Neuronal Plasticity , Oxytocin/genetics , Postpartum Period/metabolism , Promoter Regions, Genetic/genetics , Protein Biosynthesis , RNA, Messenger/genetics , RNA, Messenger/metabolism , Reproducibility of Results , Ribosomes/metabolism , Sexual Behavior, AnimalABSTRACT
Accurate, RNA-seq based, microRNA (miRNA) expression estimates from primary cells have recently been described. However, this in vitro data is mainly obtained from cell culture, which is known to alter cell maturity/differentiation status, significantly changing miRNA levels. What is needed is a robust method to obtain in vivo miRNA expression values directly from cells. We introduce expression microdissection miRNA small RNA sequencing (xMD-miRNA-seq), a method to isolate cells directly from formalin fixed paraffin-embedded (FFPE) tissues. xMD-miRNA-seq is a low-cost, high-throughput, immunohistochemistry-based method to capture any cell type of interest. As a proof-of-concept, we isolated colon epithelial cells from two specimens and performed low-input small RNA-seq. We generated up to 600,000 miRNA reads from the samples. Isolated epithelial cells, had abundant epithelial-enriched miRNA expression (miR-192; miR-194; miR-200b; miR-200c; miR-215; miR-375) and overall similar miRNA expression patterns to other epithelial cell populations (colonic enteroids and flow-isolated colon epithelium). xMD-derived epithelial cells were generally not contaminated by other adjacent cells of the colon as noted by t-SNE analysis. xMD-miRNA-seq allows for simple, economical, and efficient identification of cell-specific miRNA expression estimates. Further development will enhance rapid identification of cell-specific miRNA expression estimates in health and disease for nearly any cell type using archival FFPE material.
Subject(s)
Colon/cytology , Epithelial Cells/metabolism , Gene Expression Regulation , MicroRNAs/metabolism , Microdissection , Sequence Analysis, RNA , Cell Aggregation , Cell Separation , Cells, Cultured , Humans , MicroRNAs/geneticsABSTRACT
RNA editing is increasingly recognized as a molecular mechanism regulating RNA activity and recoding proteins. Here we surveyed the global landscape of RNA editing in human brain tissues and identified three unique patterns of A-to-I RNA editing rates during cortical development: stable high, stable low and increasing. RNA secondary structure and the temporal expression of adenosine deaminase acting on RNA (ADAR) contribute to cis- and trans-regulatory mechanisms of these RNA editing patterns, respectively. Interestingly, the increasing pattern was associated with neuronal maturation, correlated with mRNA abundance and potentially influenced miRNA binding energy. Gene ontology analyses implicated the increasing pattern in vesicle or organelle membrane-related genes and glutamate signaling pathways. We also found that the increasing pattern was selectively perturbed in spinal cord injury and glioblastoma. Our findings reveal global and dynamic aspects of RNA editing in brain, providing new insight into epitranscriptional regulation of sequence diversity.
Subject(s)
Adenosine Deaminase/metabolism , Brain/growth & development , Gene Ontology , RNA Editing/genetics , RNA-Binding Proteins/metabolism , Adenosine/metabolism , Diagnosis , Humans , Inosine/metabolism , RNA-Binding Proteins/geneticsABSTRACT
Chromosomal amplifications are among the most common genetic alterations found in human cancers. However, experimental systems to study the processes that lead to specific, recurrent amplification events in human cancers are lacking. Moreover, some common amplifications, such as that at 8p11-12 in breast cancer, harbor multiple driver oncogenes, which are poorly modeled by conventional overexpression approaches. We sought to develop an experimental system to model recurrent chromosomal amplification events in human cell lines. Our strategy is to use homologous-recombination-mediated gene targeting to deliver a dominantly selectable, amplifiable marker to a specified chromosomal location. We used adeno-associated virus vectors to target human MCF-7 breast cancer cells at the ZNF703 locus, in the recurrent 8p11-12 amplicon, using the E. coli inosine monophosphate dehydrogenase (IMPDH) enzyme as a marker. We applied selective pressure using IMPDH inhibitors. Surviving clones were found to have increased copy number of ZNF703 (average 2.5-fold increase) by droplet digital PCR and FISH. Genome-wide array comparative genomic hybridization confirmed that amplifications had occurred on the short arm of chromosome 8, without changes on 8q or other chromosomes. Patterns of amplification were variable and similar to those seen in primary human breast cancers, including "sawtooth" patterns, distal copy number loss, and large continuous regions of copy number gain. This system will allow study of the cis- and trans-acting factors that are permissive for chromosomal amplification and provide a model to analyze oncogene cooperativity in amplifications harboring multiple candidate driver genes.
Subject(s)
Breast Neoplasms/genetics , Breast Neoplasms/pathology , Gene Amplification , Gene Targeting , Genetic Loci , Carrier Proteins/genetics , Cell Line, Tumor , Chromosomes, Human, Pair 8 , Comparative Genomic Hybridization , DNA Copy Number Variations , Female , Gene Expression , Gene Order , Humans , In Situ Hybridization, FluorescenceABSTRACT
PURPOSE: Activating mutations in the phosphoinositide-3-kinase (PI3K)/AKT/mTOR pathway are present in the majority of breast cancers and therefore are a major focus of drug development and clinical trials. Pathway mutations have been proposed as predictive biomarkers for efficacy of PI3K-targeted therapies. However, the precise contribution of distinct PI3K pathway mutations to drug sensitivity is unknown. EXPERIMENTAL DESIGN: We describe the creation of a physiologic human luminal breast cancer cell line model to study the phenotype of these mutations using the MCF-7 cell line. We used somatic cell gene targeting to "correct" PIK3CA E545K-mutant alleles in MCF-7 cells to wild-type sequence. The AKT1 E17K hotspot mutation was knocked in on this wild-type background. RESULTS: Loss of mutant PIK3CA dramatically reduced phosphorylation of AKT proteins and several known AKT targets, but other AKT target proteins and downstream effectors of mTOR were not affected. PIK3CA wild-type cells exhibited reduced proliferation in vitro and in vivo. Knockin of the AKT1 E17K hotspot mutation on this PIK3CA wild-type background restored pathway signaling, proliferation, and tumor growth in vivo. PIK3CA, but not AKT1 mutation, increased sensitivity to the PI3K inhibitor GDC-0941 and the allosteric AKT inhibitor MK-2206. CONCLUSIONS: AKT1 E17K is a bona fide oncogene in a human luminal breast cancer context. Distinct PI3K pathway mutations confer differential sensitivity to drugs targeting the pathway at different points and by distinct mechanisms. These findings have implications for the use of tumor genome sequencing to assign patients to targeted therapies.
Subject(s)
Breast Neoplasms/genetics , Breast Neoplasms/metabolism , Mutation , Phosphatidylinositol 3-Kinases/genetics , Protein Kinase Inhibitors/pharmacology , Proto-Oncogene Proteins c-akt/genetics , Signal Transduction/drug effects , Animals , Breast Neoplasms/pathology , Cell Line, Tumor , Cell Proliferation/drug effects , Class I Phosphatidylinositol 3-Kinases , Disease Models, Animal , Drug Resistance, Neoplasm , Female , Heterografts , Humans , MCF-7 Cells , Mice , Phosphoinositide-3 Kinase Inhibitors , Phosphorylation/drug effects , Proto-Oncogene Proteins c-akt/antagonists & inhibitors , Tumor Burden/drug effects , Tumor Burden/geneticsABSTRACT
The phosphatidylinositol 3-kinase (PI3-kinase)-Akt-mTOR pathway is mutated at high frequency in human breast cancer, and this pathway is the focus of active drug discovery and clinical investigation. Trials of personalized cancer therapy seek to leverage knowledge of cancer gene mutations by using mutations to guide the choice of targeted therapies. At the same time, cancer genome sequencing studies are identifying low frequency variants of unknown significance in known cancer genes, as well as genes of unknown function. We have performed functional analysis of six non-hotspot AKT1 pleckstrin homology domain mutants identified in recent large-scale breast cancer sequencing studies. Three of these mutants cause constitutive activation of Akt1 in the absence of growth factors, leading to phosphorylation of downstream target proteins. Like the hotspot E17K mutation, these mutants confer constitutive membrane localization of Akt1. Finally, the same three mutants showed oncogenic activity in a cellular transformation assay. The other three mutants were inactive in all assays. These findings validate novel driver mutations in AKT1, and extend the number and type of mutations that activate the PI3-kinase pathway in human breast cancers.
Subject(s)
Amino Acid Substitution , Mutation , Proto-Oncogene Proteins c-akt/genetics , Signal Transduction/genetics , Animals , Binding Sites/genetics , Breast Neoplasms/genetics , Breast Neoplasms/metabolism , Breast Neoplasms/pathology , Cell Line , Cell Membrane/metabolism , Cell Transformation, Neoplastic/genetics , Female , Gene Expression Regulation, Neoplastic , Green Fluorescent Proteins/genetics , Green Fluorescent Proteins/metabolism , Humans , Immunoblotting , MCF-7 Cells , Microscopy, Confocal , Phosphatidylinositol 3-Kinases/genetics , Phosphatidylinositol 3-Kinases/metabolism , Proto-Oncogene Proteins c-akt/metabolism , Rats , Reverse Transcriptase Polymerase Chain Reaction , TOR Serine-Threonine Kinases/genetics , TOR Serine-Threonine Kinases/metabolismABSTRACT
Deletion of the lissencephaly disease gene LIS-1 in humans causes an extreme disorganization of the brain associated with significant reduction in cortical neurons. Here we show that deletion or RNA interference (RNAi) of Caenorhabditis elegans lis-1 results in a reduction in germline nuclei and causes a variety of cellular, developmental, and neurological defects throughout development. Our analysis of the germline defects suggests that the reduction in nuclei number stems from dysfunctional mitotic spindles resulting in cell cycle arrest and eventually programmed cell death (apoptosis). Deletion of the spindle checkpoint gene mdf-1 blocks lis-1(lf)-induced cell cycle arrest and germline apoptosis, placing the spindle checkpoint pathway upstream of the programmed cell death pathway. These results suggest that apoptosis may contribute to the cell-sparse pathology of lissencephaly.