Isoform diversity and regulation in peripheral and central neurons revealed through RNA-Seq

Selective neuronal vulnerability to deficits in RNA processing

Emerging evidence indicates that errors in RNA processing can causally drive neurodegeneration. Given that RNA produced from expressed genes of all cell types undergoes processing (splicing, polyadenylation, 5' capping, etc.), the particular vulnerability of neurons to deficits in RNA processing calls for careful consideration. The activity-dependent transcriptome remodeling associated with synaptic plasticity in neurons requires rapid, multilevel post-transcriptional RNA processing events that provide additional opportunities for dysregulation and consequent introduction or persistence of errors in RNA transcripts. Here we review the accumulating evidence that neurons have an enhanced propensity for errors in RNA processing alongside grossly insufficient defenses to clear misprocessed RNA compared to other cell types. Additionally, we explore how tau, a microtubule-associated protein implicated in Alzheimer's disease and related tauopathies, contributes to deficits in RNA processing and clearance.

Palmitoylation-dependent control of JAK1 kinase signaling governs responses to neuropoietic cytokines and survival in DRG neurons

Janus Kinase-1 (JAK1) plays key roles during neurodevelopment and following neuronal injury, while activatory JAK1 mutations are linked to leukemia. In mice, Jak1 genetic deletion results in perinatal lethality, suggesting non-redundant roles and/or regulation of JAK1 for which other JAKs cannot compensate. Proteomic studies reveal that JAK1 is more likely palmitoylated compared to other JAKs, implicating palmitoylation as a possible JAK1-specific regulatory mechanism. However, the importance of palmitoylation for JAK1 signaling has not been addressed. Here, we report that JAK1 is palmitoylated in transfected HEK293T cells and endogenously in cultured Dorsal Root Ganglion (DRG) neurons. We further use comprehensive screening in transfected non-neuronal cells and shRNA-mediated knockdown in DRG neurons to identify the related enzymes ZDHHC3 and ZDHHC7 as dominant protein acyltransferases (PATs) for JAK1. Surprisingly, we found palmitoylation minimally affects JAK1 localization in neurons, but is critical for JAK1's kinase activity in cells and even in vitro. We propose this requirement is likely because palmitoylation facilitates transphosphorylation of key sites in JAK1's activation loop, a possibility consistent with structural models of JAK1. Importantly, we demonstrate a leukemia-associated JAK1 mutation overrides the palmitoylation-dependence of JAK1 activity, potentially explaining why this mutation is oncogenic. Finally, we show that JAK1 palmitoylation is important for neuropoietic cytokine-dependent signaling and neuronal survival and that combined Zdhhc3/7 loss phenocopies loss of palmitoyl-JAK1. These findings provide new insights into the control of JAK signaling in both physiological and pathological contexts.

mRNA isoform balance in neuronal development and disease

Balanced mRNA isoform diversity and abundance are spatially and temporally regulated throughout cellular differentiation. The proportion of expressed isoforms contributes to cell type specification and determines key properties of the differentiated cells. Neurons are unique cell types with intricate developmental programs, characteristic cellular morphologies, and electrophysiological potential. Neuron-specific gene expression programs establish these distinctive cellular characteristics and drive diversity among neuronal subtypes. Genes with neuron-specific alternative processing are enriched in key neuronal functions, including synaptic proteins, adhesion molecules, and scaffold proteins. Despite the similarity of neuronal gene expression programs, each neuronal subclass can be distinguished by unique alternative mRNA processing events. Alternative processing of developmentally important transcripts alters coding and regulatory information, including interaction domains, transcript stability, subcellular localization, and targeting by RNA binding proteins. Fine-tuning of mRNA processing is essential for neuronal activity and maintenance. Thus, the focus of neuronal RNA biology research is to dissect the transcriptomic mechanisms that underlie neuronal homeostasis, and consequently, predispose neuronal subtypes to disease. This article is categorized under: RNA in Disease and Development > RNA in Disease RNA in Disease and Development > RNA in Development.

Endogenous Glycoprotein GPM6a Is Involved in Neurite Outgrowth in Rat Dorsal Root Ganglion Neurons

The peripheral nervous system (PNS) has a unique ability for self-repair. Dorsal root ganglion (DRG) neurons regulate the expression of different molecules, such as neurotrophins and their receptors, to promote axon regeneration after injury. However, the molecular players driving axonal regrowth need to be better defined. The membrane glycoprotein GPM6a has been described to contribute to neuronal development and structural plasticity in central-nervous-system neurons. Recent evidence indicates that GPM6a interacts with molecules from the PNS, although its role in DRG neurons remains unknown. Here, we characterized the expression of GPM6a in embryonic and adult DRGs by combining analysis of public RNA-seq datasets with immunochemical approaches utilizing cultures of rat DRG explants and dissociated neuronal cells. M6a was detected on the cell surfaces of DRG neurons throughout development. Moreover, GPM6a was required for DRG neurite elongation in vitro. In summary, we provide evidence on GPM6a being present in DRG neurons for the first time. Data from our functional experiments support the idea that GPM6a could contribute to axon regeneration in the PNS.

Neuronal Redevelopment and the Regeneration of Neuromodulatory Axons in the Adult Mammalian Central Nervous System

One reason that many central nervous system injuries, including those arising from traumatic brain injury, spinal cord injury, and stroke, have limited recovery of function is that neurons within the adult mammalian CNS lack the ability to regenerate their axons following trauma. This stands in contrast to neurons of the adult mammalian peripheral nervous system (PNS). New evidence, provided by single-cell expression profiling, suggests that, following injury, both mammalian central and peripheral neurons can revert to an embryonic-like growth state which is permissive for axon regeneration. This “redevelopment” strategy could both facilitate a damage response necessary to isolate and repair the acute damage from injury and provide the intracellular machinery necessary for axon regrowth. Interestingly, serotonin neurons of the rostral group of raphe nuclei, which project their axons into the forebrain, display a robust ability to regenerate their axons unaided, counter to the widely held view that CNS axons cannot regenerate without experimental intervention after injury. Furthermore, initial evidence suggests that norepinephrine neurons within the locus coeruleus possess similar regenerative abilities. Several morphological characteristics of serotonin axon regeneration in adult mammals, observable using longitudinal in vivo imaging, are distinct from the known characteristics of unaided peripheral nerve regeneration, or of the regeneration seen in the spinal cord and optic nerve that occurs with experimental intervention. These results suggest that there is an alternative CNS program for axon regeneration that likely differs from that displayed by the PNS.

Suppression of premature transcription termination leads to reduced mRNA isoform diversity and neurodegeneration

Tight regulation of mRNA isoform expression is essential for neuronal development, maintenance, and function; however, the repertoire of proteins that govern isoform composition and abundance remains incomplete. Here, we show that the RNA kinase CLP1 regulates mRNA isoform expression through suppression of proximal cleavage and polyadenylation. We found that human stem-cell-derived motor neurons without CLP1 or with the disease-associated CLP1 p.R140H variant had distinct patterns of RNA-polymerase-II-associated cleavage and polyadenylation complex proteins that correlated with polyadenylation site usage. These changes resulted in imbalanced mRNA isoform expression of long genes important for neuronal function that were recapitulated in vivo. Strikingly, we observed the same pattern of reduced mRNA isoform diversity in 3' end sequencing data from brain tissues of patients with neurodegenerative disease. Together, our results identify a previously uncharacterized role for CLP1 in mRNA 3' end formation and reveal an mRNA misprocessing signature in neurodegeneration that may suggest a common mechanism of disease.

Phenotypic Screening Following Transcriptomic Deconvolution to Identify Transcription Factors Mediating Axon Growth Induced by a Kinase Inhibitor

After injury to the central nervous system (CNS), both neuron-intrinsic limitations on regenerative responses and inhibitory factors in the injured CNS environment restrict regenerative axon growth. Instances of successful axon regrowth offer opportunities to identify features that differentiate these situations from that of the normal adult CNS. One such opportunity is provided by the kinase inhibitor RO48, which dramatically enhances neurite outgrowth of neurons in vitro and substantially increased contralateral sprouting of corticospinal tract neurons when infused intraventricularly following unilateral pyramidotomy. The authors present here a transcriptomic deconvolution of RO48-associated axon growth, with the goal of identifying transcriptional regulators associated with axon growth in the CNS. Through the use of RNA sequencing (RNA-seq) and transcription factor binding site enrichment analysis, the authors identified a list of transcription factors putatively driving differential gene expression during RO48-induced neurite outgrowth of rat hippocampal neurons in vitro. The 82 transcription factor motifs identified in this way included some with known association to axon growth regulation, such as Jun, Klf4, Myc, Atf4, Stat3, and Nfatc2, and many with no known association to axon growth. A phenotypic loss-of-function screen was carried out to evaluate these transcription factors for their roles in neurite outgrowth; this screen identified several potential outgrowth regulators. Subsequent validation suggests that the Forkhead box (Fox) family transcription factor Foxp2 restricts neurite outgrowth, while FoxO subfamily members Foxo1 and Foxo3a promote neurite outgrowth. The authors' combined transcriptomic-phenotypic screening strategy therefore allowed identification of novel transcriptional regulators of neurite outgrowth downstream of a multitarget kinase inhibitor.

ExTraMapper: Exon- and Transcript-level mappings for orthologous gene pairs

MOTIVATION

Access to large-scale genomics and transcriptomics data from various tissues and cell lines allowed the discovery of wide-spread alternative splicing events and alternative promoter usage in mammalians. Between human and mouse, gene-level orthology is currently present for nearly 16k protein-coding genes spanning a diverse repertoire of over 200k total transcript isoforms.

RESULTS

Here, we describe a novel method, ExTraMapper, which leverages sequence conservation between exons of a pair of organisms and identifies a fine-scale orthology mapping at the exon and then transcript level. ExTraMapper identifies more than 350k exon mappings, as well as 30k transcript mappings between human and mouse using only sequence and gene annotation information. We demonstrate that ExTraMapper identifies a larger number of exon and transcript mappings compared to previous methods. Further, it identifies exon fusions, splits, and losses due to splice site mutations, and finds mappings between microexons that are previously missed. By reanalysis of RNA-seq data from 13 matched human and mouse tissues, we show that ExTraMapper improves the correlation of transcript-specific expression levels suggesting a more accurate mapping of human and mouse transcripts. We also applied the method to detect conserved exon and transcript pairs between human and rhesus macaque genomes to highlight the point that ExTraMapper is applicable to any pair of organisms that have orthologous gene pairs.

AVAILABILITY

The source code and the results are available at https://github.com/ay-lab/ExTraMapper and http://ay-lab-tools.lji.org/extramapper.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Birth of a Regulatory Long Non-coding RNA/Gene, linc-UR-UB

The origin of genes has been a major topic of research for many years, albeit in some cases, it has been a difficult process to elucidate. Insightful is a recent publication that experimentally shows how one gene, linc-UR-UB was born. This gene is regulated in a complex manner in male germ cells during spermatogenesis and is believed to participate in the regulation of levels of the ubiquitin specific peptidase 18 (USP18) mRNA. The process of formation of linc-UR-UB appears relatively simple. It involves a transcription read through from an upstream gene to a downstream functional element, the USP18 3' UTR sequence. This small element also shares the same sequence as the 3' ends of the lincRNA FAM247 family genes. In addition to linc-UR-UB, it is possible that other genes formed in a similar fashion that involves a genomic sequence read through to a functional element.

The palmitoyl acyltransferases ZDHHC5 and ZDHHC8 are uniquely present in DRG axons and control retrograde signaling via the Gp130/JAK/STAT3 pathway

Palmitoylation, the modification of proteins with the lipid palmitate, is a key regulator of protein targeting and trafficking. However, knowledge of the roles of specific palmitoyl acyltransferases (PATs), which catalyze palmitoylation, is incomplete. For example, little is known about which PATs are present in neuronal axons, although long-distance trafficking of palmitoyl-proteins is important for axon integrity and for axon-to-soma retrograde signaling, a process critical for axon development and for responses to injury. Identifying axonally targeted PATs might thus provide insights into multiple aspects of axonal biology. We therefore comprehensively determined the subcellular distribution of mammalian PATs in dorsal root ganglion (DRG) neurons and, strikingly, found that only two PATs, ZDHHC5 and ZDHHC8, were enriched in DRG axons. Signals via the Gp130/JAK/STAT3 and DLK/JNK pathways are important for axonal injury responses, and we found that ZDHHC5 and ZDHHC8 were required for Gp130/JAK/STAT3, but not DLK/JNK, axon-to-soma signaling. ZDHHC5 and ZDHHC8 robustly palmitoylated Gp130 in cotransfected nonneuronal cells, supporting the possibility that Gp130 is a direct ZDHHC5/8 substrate. In DRG neurons, Zdhhc5/8 shRNA knockdown reduced Gp130 palmitoylation and even more markedly reduced Gp130 surface expression, potentially explaining the importance of these PATs for Gp130-dependent signaling. Together, these findings provide new insights into the subcellular distribution and roles of specific PATs and reveal a novel mechanism by which palmitoylation controls axonal retrograde signaling.

isoTarget: A Genetic Method for Analyzing the Functional Diversity of Splicing Isoforms In Vivo

Protein isoforms generated by alternative splicing contribute to proteome diversity. Because of the lack of effective techniques, the isoform-specific function, expression, localization, and signaling of endogenous proteins are unknown for most genes. Here, we report a genetic method, isoTarget, for multi-purpose studies of targeted isoforms in select cells. Applying isoTarget to two isoforms of Drosophila Dscam, Dscam[TM1] and [TM2], we found that, in neurons, endogenous Dscam[TM1] is in dendrites, whereas Dscam[TM2] is in both dendrites and axons. We demonstrate that the difference in subcellular localization, rather than biochemical properties, leads to the two isoforms’ functional differences. Moreover, we show that the subcellular enrichment of functional partners results in a DLK/Wallenda-Dscam[TM2]-Dock signaling cascade in axons. We further apply isoTarget to study two isoforms of a GABA receptor to demonstrate its general applicability. isoTarget is an effective technique for studying how alternative splicing enhances proteome complexity.

Liu et al. develop a genetic method that enables the investigation of isoform-specific function, expression, localization, and signaling of endogenous proteins in select cells. Using this method, they demonstrate that the difference in subcellular localization of two isoforms of Down syndrome cell adhesion molecule leads to functional differences between them.

Spliced genes in muscle from Nelore Cattle and their association with carcass and meat quality

Transcript data obtained by RNA-Seq were used to identify differentially expressed alternatively spliced genes in ribeye muscle tissue between Nelore cattle that differed in their ribeye area (REA) or intramuscular fat content (IF). A total of 166 alternatively spliced transcripts from 125 genes were significantly differentially expressed in ribeye muscle between the highest and lowest REA groups (p ≤ 0.05). For animals selected on their IF content, 269 alternatively spliced transcripts from 219 genes were differentially expressed in ribeye muscle between the highest and lowest IF animals. Cassette exons and alternative 3′ splice sites were the most frequently found alternatively spliced transcripts for REA and IF content. For both traits, some differentially expressed alternatively spliced transcripts belonged to myosin and myotilin gene families. The hub transcripts were identified for REA (LRRFIP1, RCAN1 and RHOBTB1) and IF (TRIP12, HSPE1 and MAP2K6) have an important role to play in muscle cell degradation, development and motility. In general, transcripts were found for both traits with biological process GO terms that were involved in pathways related to protein ubiquitination, muscle differentiation, lipids and hormonal systems. Our results reinforce the biological importance of these known processes but also reveal new insights into the complexity of the whole cell muscle mRNA of Nelore cattle.

Identification of a Sacral, Visceral Sensory Transcriptome in Embryonic and Adult Mice

Visceral sensory neurons encode distinct sensations from healthy organs and initiate pain states that are resistant to common analgesics. Transcriptome analysis is transforming our understanding of sensory neuron subtypes but has generally focused on somatic sensory neurons or the total population of neurons in which visceral neurons form the minority. Our aim was to define transcripts specifically expressed by sacral visceral sensory neurons, as a step towards understanding the unique biology of these neurons and potentially leading to identification of new analgesic targets for pelvic visceral pain. Our strategy was to identify genes differentially expressed between sacral dorsal root ganglia (DRG) that include somatic neurons and sacral visceral neurons, and adjacent lumbar DRG that comprise exclusively of somatic sensory neurons. This was performed in adult and E18.5 male and female mice. By developing a method to restrict analyses to nociceptive Trpv1 neurons, a larger group of genes were detected as differentially expressed between spinal levels. We identified many novel genes that had not previously been associated with pelvic visceral sensation or nociception. Limited sex differences were detected across the transcriptome of sensory ganglia, but more were revealed in sacral levels and especially in Trpv1 nociceptive neurons. These data will facilitate development of new tools to modify mature and developing sensory neurons and nociceptive pathways.

Comparative analysis of the transcriptome of injured nerve segments reveals spatiotemporal responses to neural damage in mice

Nerves are particularly vulnerable to damage due to their unique structure with meter-long axons. In the peripheral nervous system neurons and Schwann cells can activate the injury-response program that directs axons to either regenerate or degenerate after traumatic nerve injury. However, the differences between the genetic programs driving nerve regeneration and degeneration have not yet been described extensively. To understand these differences, in this study we have compared the injury-induced transcriptomic changes between the regenerating proximal segment and the degenerating distal segment of a transected nerve, at different post-injury time points. We analyzed the spatiotemporal dynamics of the mouse transcriptome using a sciatic nerve-injury model by means of RNA sequencing. The results of the differentially regulated genes (DEGs) analysis show that some DEG groups are similarly regulated in both proximal and distal segments, and primarily display a positive correlation. However, some DEG groups are exclusively regulated in either the proximal or the distal segment, suggesting that these DEG groups constitute a genetic network for distinguishing the regenerative and degenerative responses. In addition, our gene ontology analysis revealed an enrichment of particular biological processes in different phases and locations. Thus, our data provide a spatiotemporal profile of the transcriptomes that are differentially regulated in either regenerating or degenerating nerves, in vivo. The specific biological processes enriched in the DEG groups might delineate the injury-responsive program that induces contrasting regenerative and degenerative responses in different nerve segments.

Highly conserved molecular pathways, including Wnt signaling, promote functional recovery from spinal cord injury in lampreys

In mammals, spinal cord injury (SCI) leads to dramatic losses in neurons and synaptic connections, and consequently function. Unlike mammals, lampreys are vertebrates that undergo spontaneous regeneration and achieve functional recovery after SCI. Therefore our goal was to determine the complete transcriptional responses that occur after SCI in lampreys and to identify deeply conserved pathways that promote regeneration. We performed RNA-Seq on lamprey spinal cord and brain throughout the course of functional recovery. We describe complex transcriptional responses in the injured spinal cord, and somewhat surprisingly, also in the brain. Transcriptional responses to SCI in lampreys included transcription factor networks that promote peripheral nerve regeneration in mammals such as Atf3 and Jun. Furthermore, a number of highly conserved axon guidance, extracellular matrix, and proliferation genes were also differentially expressed after SCI in lampreys. Strikingly, ~3% of differentially expressed transcripts belonged to the Wnt pathways. These included members of the Wnt and Frizzled gene families, and genes involved in downstream signaling. Pharmacological inhibition of Wnt signaling inhibited functional recovery, confirming a critical role for this pathway. These data indicate that molecular signals present in mammals are also involved in regeneration in lampreys, supporting translational relevance of the model.

Large-scale identification of human cerebrovascular proteins: Inter-tissue and intracerebral vascular protein diversity

The human cerebrovascular system is responsible for regulating demand-dependent perfusion and maintaining the blood-brain barrier (BBB). In addition, defects in the human cerebrovasculature lead to stroke, intracerebral hemorrhage, vascular malformations, and vascular cognitive impairment. The objective of this study was to discover new proteins of the human cerebrovascular system using expression data from the Human Protein Atlas, a large-scale project which allows public access to immunohistochemical analysis of human tissues. We screened 20,158 proteins in the HPA and identified 346 expression patterns correlating to blood vessels in human brain. Independent experiments showed that 51/52 of these distributions could be experimentally replicated across different brain samples. Some proteins (40%) demonstrated endothelial cell (EC)-enriched expression, while others were expressed primarily in vascular smooth muscle cells (VSMC; 18%); 39% of these proteins were expressed in both cell types. Most brain EC markers were tissue oligospecific; that is, they were expressed in endothelia in an average of 4.8 out of 9 organs examined. Although most markers expressed in endothelial cells of the brain were present in all cerebral capillaries, a significant number (21%) were expressed only in a fraction of brain capillaries within each brain sample. Among proteins found in cerebral VSMC, virtually all were also expressed in peripheral VSMC and in non-vascular smooth muscle cells (SMC). Only one was potentially brain specific: VHL (Von Hippel-Lindau tumor suppressor). HRC (histidine rich calcium binding protein) and VHL were restricted to VSMC and not found in non-vascular tissues such as uterus or gut. In conclusion, we define a set of brain vascular proteins that could be relevant to understanding the unique physiology and pathophysiology of the human cerebrovasculature. This set of proteins defines inter-organ molecular differences in the vasculature and confirms the broad heterogeneity of vascular cells within the brain.

High Content Screening of Mammalian Primary Cortical Neurons

High Content Screening (HCS) can be used to analyze the morphology of neuronal primary cultures on a large scale. When used in the field of neuronal regeneration this approach allows the screening of hundreds or thousands of perturbagens, such as miRNAs, cDNAs, or compounds, for their ability to induce neuronal growth. One of the most important steps while designing these kinds of experiments is the choice of the correct neuronal model. Testing the correct neuronal type is critical to obtain results that are biologically significant and that can later be translated to a clinical setting. For example, if the goal is identifying possible therapies for Spinal Cord Injury (SCI), a challenging target is the neuronal projection from the motor cortex to the spinal cord, the corticospinal tract. Here, we describe the experimental protocols that can be used to produce primary cortical culture from young rat cortices, electroporate the neurons to study the effect of altered gene expression on neurite growth, and immunostain to measure neurite growth parameters.

Identification of miRNAs involved in DRG neurite outgrowth and their putative targets

Peripheral neurons regenerate their axons after injury. Transcriptional regulation by microRNAs (miRNAs) is one possible mechanism controlling regeneration. We profiled miRNA expression in mouse dorsal root ganglion neurons after a sciatic nerve crush, and identified 49 differentially expressed miRNAs. We evaluated the functional role of each miRNA using a phenotypic analysis approach. To predict the targets of the miRNAs we employed RNA-Sequencing and examined transcription at the isoform level. We identify thousands of differentially expressed isoforms and bioinformatically associate the miRNAs that modulate neurite growth with their putative target isoforms to outline a network of regulatory events underlying peripheral nerve regeneration. MiR-298, let-7a, and let-7f enhance neurite growth and target the majority of isoforms in the differentially expressed network.

Phenotypic screening with primary neurons to identify drug targets for regeneration and degeneration

High-throughput, target-based screening techniques have been utilized extensively for drug discovery in the past several decades. However, the need for more predictive in vitro models of in vivo disease states has generated a shift in strategy towards phenotype-based screens. Phenotype based screens are particularly valuable in studying complex conditions such as CNS injury and degenerative disease, as many factors can contribute to a specific cellular response. In this review, we will discuss different screening frameworks and their relative utility in examining mechanisms of neurodegeneration and axon regrowth, particularly in cell-based in vitro disease models.

Transcriptome profile of the human placenta

The human placenta is a particular organ that inseparably binds the mother and the fetus. The proper development and survival of the conceptus relies on the essential interplay between maternal and fetal factors involved in cooperation within the placenta. In our study, high-throughput sequencing (RNA-seq) was applied to analyze the global transcriptome of the human placenta during uncomplicated pregnancies. The RNA-seq was utilized to identify the global pattern of the gene expression in placentas (N = 4) from women in single and twin pregnancies. During analyses, we obtained 228,044 transcripts. More than 91% of them were multi-exon, and among them 134 were potentially unknown protein coding genes. Expression levels (FPKM) were estimated for 38,948 transcriptional active regions, and more than 3000 of genes were expressed with FPKM >20 in each sample. Additionally, all unannotated transcripts with estimated FPKM values were localized on the human genome. Highly covered splice junctions unannotated in the human genome (6497) were identified, and among them 30 were novel. To gain a better understanding of the biological implications, the assembled transcripts were annotated with gene ontology (GO) terms. Single nucleotide variants were predicted for the transcripts assigned to each analyzed GO category. Our results may be useful for establishing a general pattern of the gene expression in the human placenta. Characterizing placental transcriptome, which is crucial for a pregnancy’s outcome, can serve as a basis for identifying the mechanisms underlying physiological pregnancy, as well as may be useful for an early detection of the genomic defects.

Macrophage Transcriptional Profile Identifies Lipid Catabolic Pathways That Can Be Therapeutically Targeted after Spinal Cord Injury

Although infiltrating macrophages influence many pathological processes after spinal cord injury (SCI), the intrinsic molecular mechanisms that regulate their function are poorly understood. A major hurdle has been dissecting macrophage-specific functions from those in other cell types as well as understanding how their functions change over time. Therefore, we used the RiboTag method to obtain macrophage-specific mRNA directly from the injured spinal cord in mice and performed RNA sequencing to investigate their transcriptional profile. Our data show that at 7 d after SCI, macrophages are best described as foam cells, with lipid catabolism representing the main biological process, and canonical nuclear receptor pathways as their potential mediators. Genetic deletion of a lipoprotein receptor, CD36, reduces macrophage lipid content and improves lesion size and locomotor recovery. Therefore, we report the first macrophage-specific transcriptional profile after SCI and highlight the lipid catabolic pathway as an important macrophage function that can be therapeutically targeted after SCI.SIGNIFICANCE STATEMENT The intrinsic molecular mechanisms that regulate macrophage function after spinal cord injury (SCI) are poorly understood. We obtained macrophage-specific mRNA directly from the injured spinal cord and performed RNA sequencing to investigate their transcriptional profile. Our data show that at 7 d after SCI, macrophages are best described as foam cells, with lipid catabolism representing the main biological process and canonical nuclear receptor pathways as their potential mediators. Genetic deletion of a lipoprotein receptor, CD36, reduces macrophage lipid content and improves lesion size and locomotor recovery. Therefore, we report the first macrophage-specific transcriptional profile after SCI and highlight the lipid catabolic pathway as an important macrophage function that can be therapeutically targeted after SCI.

Cell type-specific transcriptome profiling in mammalian brains

A mammalian brain contains numerous types of cells. Advances in neuroscience in the past decade allow us to identify and isolate neural cells of interest from mammalian brains. Recent developments in high-throughput technologies, such as microarrays and next-generation sequencing (NGS), provide detailed information on gene expression in pooled cells on a genomic scale. As a result, many novel genes have been found critical in cell type-specific transcriptional regulation. These differentially expressed genes can be used as molecular signatures, unique to a particular class of neural cells. Use of this gene expression-based approach can further differentiate neural cell types into subtypes, potentially linking some of them with neurological diseases. In this article, experimental techniques used to purify neural cells are described, followed by a review on recent microarray- or NGS-based transcriptomic studies of common neural cell types. The future prospects of cell type-specific research are also discussed.

Slippery signaling: Palmitoylation-dependent control of neuronal kinase localization and activity

Modification of proteins with the lipid palmitate, a process called palmitoylation, is important for the normal function of neuronal cells. However, most attention has focused on how palmitoylation regulates the targeting and trafficking of neurotransmitter receptors and non-enzymatic scaffold proteins. In this review we discuss recent studies that suggest that palmitoylation also plays additional roles in neurons by controlling the localization, interactions and perhaps even the activity of protein kinases that play key roles in physiological neuronal regulation and in neuropathological processes.

Hyperactivated Stat3 boosts axon regeneration in the CNS

Axonal regeneration after spinal cord injury (SCI) is intrinsically and extrinsically inhibited by multiple factors. One major factor contributing to intrinsic regeneration failure is the inability of mature neurons in the central nervous system (CNS) to activate regeneration-associated transcription factors (TFs) post-injury. A prior study identified TFs overexpressed in neurons of the peripheral nervous system (PNS) compared to the CNS; some of these could be involved in the ability of PNS neurons to regenerate. Of these, signal transducer and activator of transcription 3 (STAT3), as well its downstream regeneration-associated targets, showed a significant upregulation in PNS neurons relative to CNS neurons, and a constitutively active variant of Stat3 (Stat3CA) promoted neurite growth when expressed in cerebellar neurons (Lerch et al., 2012; Smith et al., 2011). To further enhance STAT3's neurite outgrowth enhancing activity, Stat3CA was fused with a viral activation domain (VP16). VP16 hyperactivates TFs by recruiting transcriptional co-factors to the DNA binding domain (Hirai et al., 2010). Overexpression of this VP16-Stat3CA chimera in primary cortical neurons led to a significant increase of neurite outgrowth as well as Stat3 transcriptional activity in vitro. Furthermore, in vivo transduction of retinal ganglion cells (RGCs) with AAV constructs expressing VP16-Stat3CA resulted in regeneration of optic nerve axons after injury, to a greater degree than for those expressing Stat3CA alone. These findings confirm and extend the concept that overexpression of hyperactivated transcription factors identified as functioning in PNS regeneration can promote axon regeneration in the CNS.

Neuronal networks in mental diseases and neuropathic pain: Beyond brain derived neurotrophic factor and collapsin response mediator proteins

The brain is a complex network system that has the capacity to support emotion, thought, action, learning and memory, and is characterized by constant activity, constant structural remodeling, and constant attempt to compensate for this remodeling. The basic insight that emerges from complex network organization is that substantively different networks can share common key organizational principles. Moreover, the interdependence of network organization and behavior has been successfully demonstrated for several specific tasks. From this viewpoint, increasing experimental/clinical observations suggest that mental disorders are neural network disorders. On one hand, single psychiatric disorders arise from multiple, multifactorial molecular and cellular structural/functional alterations spreading throughout local/global circuits leading to multifaceted and heterogeneous clinical symptoms. On the other hand, various mental diseases may share functional deficits across the same neural circuit as reflected in the overlap of symptoms throughout clinical diagnoses. An integrated framework including experimental measures and clinical observations will be necessary to formulate a coherent and comprehensive understanding of how neural connectivity mediates and constraints the phenotypic expression of psychiatric disorders.

Treatment with analgesics after mouse sciatic nerve injury does not alter expression of wound healing-associated genes

Animal models of sciatic nerve injury are commonly used to study neuropathic pain as well as axon regeneration. Administration of post-surgical analgesics is an important consideration for animal welfare, but the actions of the analgesic must not interfere with the scientific goals of the experiment. In this study, we show that treatment with either buprenorphine or acetaminophen following a bilateral sciatic nerve crush surgery does not alter the expression in dorsal root ganglion (DRG) sensory neurons of a panel of genes associated with wound healing. These findings indicate that the post-operative use of buprenorphine or acetaminophen at doses commonly suggested by Institutional Animal Care and Use Committees does not change the intrinsic gene expression response of DRG neurons to a sciatic nerve crush injury, for many wound healing-associated genes. Therefore, administration of post-operative analgesics may not confound the results of transcriptomic studies employing this injury model.

Transcriptome Analysis of the Human Striatum in Tourette Syndrome

BACKGROUND

Genome-wide association studies have not revealed any risk-conferring common genetic variants in Tourette syndrome (TS), requiring the adoption of alternative approaches to investigate the pathophysiology of this disorder.

METHODS

We obtained the basal ganglia transcriptome by RNA sequencing in the caudate and putamen of nine TS and nine matched normal control subjects.

RESULTS

We found 309 downregulated and 822 upregulated genes in the caudate and putamen (striatum) of TS individuals. Using data-driven gene network analysis, we identified 17 gene coexpression modules associated with TS. The top-scoring downregulated module in TS was enriched in striatal interneuron transcripts, which was confirmed by decreased numbers of cholinergic and gamma-aminobutyric acidergic interneurons by immunohistochemistry in the same regions. The top-scoring upregulated module was enriched in immune-related genes, consistent with activation of microglia in patients' striatum. Genes implicated by copy number variants in TS were enriched in the interneuron module, as well as in a protocadherin module. Module clustering revealed that the interneuron module was correlated with a neuronal metabolism module.

CONCLUSIONS

Convergence of differential expression, network analyses, and module clustering, together with copy number variants implicated in TS, strongly implicates disrupted interneuron signaling in the pathophysiology of severe TS and suggests that metabolic alterations may be linked to their death or dysfunction.

Transcriptomic Approaches to Neural Repair

Understanding why adult CNS neurons fail to regenerate their axons following injury remains a central challenge of neuroscience research. A more complete appreciation of the biological mechanisms shaping the injured nervous system is a crucial prerequisite for the development of robust therapies to promote neural repair. Historically, the identification of regeneration associated signaling pathways has been impeded by the limitations of available genetic and molecular tools. As we progress into an era in which the high-throughput interrogation of gene expression is commonplace and our knowledge base of interactome data is rapidly expanding, we can now begin to assemble a more comprehensive view of the complex biology governing axon regeneration. Here, we highlight current and ongoing work featuring transcriptomic approaches toward the discovery of novel molecular mechanisms that can be manipulated to promote neural repair.

SIGNIFICANCE STATEMENT

Transcriptional profiling is a powerful technique with broad applications in the field of neuroscience. Recent advances such as single-cell transcriptomics, CNS cell type-specific and developmental stage-specific expression libraries are rapidly enhancing the power of transcriptomics for neuroscience applications. However, extracting biologically meaningful information from large transcriptomic datasets remains a formidable challenge. This mini-symposium will highlight current work using transcriptomic approaches to identify regulatory networks in the injured nervous system. We will discuss analytical strategies for transcriptomics data, the significance of noncoding RNA networks, and the utility of multiomic data integration. Though the studies featured here specifically focus on neural repair, the approaches highlighted in this mini-symposium will be of broad interest and utility to neuroscientists working in diverse areas of the field.

Phenotypic Screening of Small-Molecule Inhibitors: Implications for Therapeutic Discovery and Drug Target Development in Traumatic Brain Injury

The inability of central nervous system (CNS) neurons to regenerate damaged axons and dendrites following traumatic brain injury (TBI) creates a substantial obstacle for functional recovery. Apoptotic cell death, deposition of scar tissue, and growth-repressive molecules produced by glia further complicate the problem and make it challenging for re-growing axons to extend across injury sites. To date, there are no approved drugs for the treatment of TBI, accentuating the need for relevant leads. Cell-based and organotypic bioassays can better mimic outcomes within the native CNS microenvironment than target-based screening methods and thus should speed the discovery of therapeutic agents that induce axon or dendrite regeneration. Additionally, when used to screen focused chemical libraries such as small-molecule protein kinase inhibitors, these assays can help elucidate molecular mechanisms involved in neurite outgrowth and regeneration as well as identify novel drug targets. Here, we describe a phenotypic cellular (high content) screening assay that utilizes brain-derived primary neurons for screening small-molecule chemical libraries.

Palmitoylation controls DLK localization, interactions and activity to ensure effective axonal injury signaling

Dual leucine-zipper kinase (DLK) is critical for axon-to-soma retrograde signaling following nerve injury. However, it is unknown how DLK, a predicted soluble kinase, conveys long-distance signals and why homologous kinases cannot compensate for loss of DLK. Here, we report that DLK, but not homologous kinases, is palmitoylated at a conserved site adjacent to its kinase domain. Using short-hairpin RNA knockdown/rescue, we find that palmitoylation is critical for DLK-dependent retrograde signaling in sensory axons. This functional importance is because of three novel cellular and molecular roles of palmitoylation, which targets DLK to trafficking vesicles, is required to assemble DLK signaling complexes and, unexpectedly, is essential for DLK's kinase activity. By simultaneously controlling DLK localization, interactions, and activity, palmitoylation ensures that only vesicle-bound DLK is active in neurons. These findings explain how DLK specifically mediates nerve injury responses and reveal a novel cellular mechanism that ensures the specificity of neuronal kinase signaling.

Defining the nociceptor transcriptome

Unbiased “omics” techniques, such as next generation RNA-sequencing, can provide entirely novel insights into biological systems. However, cellular heterogeneity presents a significant barrier to analysis and interpretation of these datasets. The neurons of the dorsal root ganglia (DRG) are an important model for studies of neuronal injury, regeneration and pain. The majority of investigators utilize a dissociated preparation of whole ganglia when studying cellular and molecular function. We demonstrate that the standard methods for producing these preparations gives a 10%-neuronal mixture of cells, with the remainder of cells constituting satellite glia and other non-neuronal cell types. Using a novel application of magnetic purification, we consistently obtain over 95% pure, viable neurons from adult tissue, significantly enriched for small diameter nociceptors expressing the voltage gated ion channel Na_v1.8. Using genome-wide RNA-sequencing we compare the currently used (10% neuronal) and pure (95% nociceptor) preparations and find 920 genes enriched. This gives an unprecedented insight into the molecular composition of small nociceptive neurons in the DRG, potentially altering the interpretation of previous studies performed at the tissue level, and indicating a number of novel markers of this widely-studied population of cells. We anticipate that the ease of use, affordability and speed of this technique will see it become widely adopted, delivering a greatly improved capacity to study the roles of nociceptors in health and disease.

An RNA-sequencing transcriptome and splicing database of glia, neurons, and vascular cells of the cerebral cortex

The major cell classes of the brain differ in their developmental processes, metabolism, signaling, and function. To better understand the functions and interactions of the cell types that comprise these classes, we acutely purified representative populations of neurons, astrocytes, oligodendrocyte precursor cells, newly formed oligodendrocytes, myelinating oligodendrocytes, microglia, endothelial cells, and pericytes from mouse cerebral cortex. We generated a transcriptome database for these eight cell types by RNA sequencing and used a sensitive algorithm to detect alternative splicing events in each cell type. Bioinformatic analyses identified thousands of new cell type-enriched genes and splicing isoforms that will provide novel markers for cell identification, tools for genetic manipulation, and insights into the biology of the brain. For example, our data provide clues as to how neurons and astrocytes differ in their ability to dynamically regulate glycolytic flux and lactate generation attributable to unique splicing of PKM2, the gene encoding the glycolytic enzyme pyruvate kinase. This dataset will provide a powerful new resource for understanding the development and function of the brain. To ensure the widespread distribution of these datasets, we have created a user-friendly website (http://web.stanford.edu/group/barres_lab/brain_rnaseq.html) that provides a platform for analyzing and comparing transciption and alternative splicing profiles for various cell classes in the brain.

Regulatory gene networks that shape the development of adaptive phenotypic plasticity in a cichlid fish

Phenotypic plasticity is the ability of organisms with a given genotype to develop different phenotypes according to environmental stimuli, resulting in individuals that are better adapted to local conditions. In spite of their ecological importance, the developmental regulatory networks underlying plastic phenotypes often remain uncharacterized. We examined the regulatory basis of diet-induced plasticity in the lower pharyngeal jaw (LPJ) of the cichlid fish Astatoreochromis alluaudi, a model species in the study of adaptive plasticity. Through raising juvenile A. alluaudi on either a hard or soft diet (hard-shelled or pulverized snails) for between 1 and 8 months, we gained insight into the temporal regulation of 19 previously identified candidate genes during the early stages of plasticity development. Plasticity in LPJ morphology was first detected between 3 and 5 months of diet treatment. The candidate genes, belonging to various functional categories, displayed dynamic expression patterns that consistently preceded the onset of morphological divergence and putatively contribute to the initiation of the plastic phenotypes. Within functional categories, we observed striking co-expression, and transcription factor binding site analysis was used to examine the prospective basis of their coregulation. We propose a regulatory network of LPJ plasticity in cichlids, presenting evidence for regulatory crosstalk between bone and muscle tissues, which putatively facilitates the development of this highly integrated trait. Through incorporating a developmental time-course into a phenotypic plasticity study, we have identified an interconnected, environmentally responsive regulatory network that shapes the development of plasticity in a key innovation of East African cichlids.

Identification of potential therapeutic targets in a model of neuropathic pain

Neuropathic pain (NP) is caused by damage to the nervous system, resulting in dysfunction and aberrant pain. The cellular functions (e.g., peripheral neuron spinal cord innervation, neuronal excitability) associated with NP often develop over time and are likely associated with gene expression changes. Gene expression studies on the cells involved in NP (e.g., sensory dorsal root ganglion neurons) are publically available; the mining of these studies may enable the identification of novel targets and the subsequent development of therapies that are essential for improving quality of life for the millions of individuals suffering with NP. Here we analyzed a publically available microarray dataset (GSE30165) in order to identify new RNAs (e.g., messenger RNA (mRNA) isoforms and non-coding RNAs) underlying NP. GSE30165 profiled gene expression in dorsal root ganglion neurons (DRG) and in sciatic nerve (SN) after resection, a NP model. Gene ontological analysis shows enrichment for sensory and neuronal processes. Protein network analysis demonstrates DRG upregulated genes typical to an injury and NP response. Of the top changing genes, 34 and 36% are associated with more than one protein coding isoform in the DRG and SN, respectively. The majority of genes are receptor and enzymes. We identified 15 long non-coding RNAs (lncRNAs) targeting these genes in LNCipedia.org, an online comprehensive lncRNA database. These RNAs represent new therapeutic targets for preventing NP development and this approach demonstrates the feasibility of data reanalysis for their identification.

cJun promotes CNS axon growth

A number of genes regulate regeneration of peripheral axons, but their ability to drive axon growth and regeneration in the central nervous system (CNS) remains largely untested. To address this question we overexpressed eight transcription factors and one small GTPase alone and in pairwise combinations to test whether combinatorial overexpression would have a synergistic impact on CNS neuron neurite growth. The Jun oncogene/signal transducer and activator of transcription 6 (JUN/STAT6) combination increased neurite growth in dissociated cortical neurons and in injured cortical slices. In injured cortical slices, JUN overexpression increased axon growth to a similar extent as JUN and STAT6 together. Interestingly, JUN overexpression was not associated with increased growth associated protein 43 (GAP43) or integrin alpha 7 (ITGA7) expression, though these are predicted transcriptional targets. This study demonstrates that JUN overexpression in cortical neurons stimulates axon growth, but does so independently of changes in expression of genes thought to be critical for JUNs effects on axon growth. We conclude that JUN activity underlies this CNS axonal growth response, and that it is mechanistically distinct from peripheral regeneration responses, in which increases in JUN expression coincide with increases in GAP43 expression.

Lentiviral vector-mediated RNA silencing in the central nervous system

RNA silencing is an established method for investigating gene function and has attracted particular interest because of the potential for generating RNA-based therapeutics. Using lentiviral vectors as an efficient delivery system that offers stable, long-term expression in postmitotic cells further enhances the applicability of an RNA-based gene therapy for the CNS. In this review we provide an overview of both lentiviral vectors and RNA silencing along with design considerations for generating lentiviral vectors capable of RNA silencing. We go on to describe the current preclinical data regarding lentiviral vector-mediated RNA silencing for CNS disorders and discuss the concerns of side effects associated with lentiviral vectors and small interfering RNAs and how these might be mitigated.

Chemical interrogation of the neuronal kinome using a primary cell-based screening assay

A fundamental impediment to functional recovery from spinal cord injury (SCI) and traumatic brain injury is the lack of sufficient axonal regeneration in the adult central nervous system. There is thus a need to develop agents that can stimulate axon growth to re-establish severed connections. Given the critical role played by protein kinases in regulating axon growth and the potential for pharmacological intervention, small molecule protein kinase inhibitors present a promising therapeutic strategy. Here, we report a robust cell-based phenotypic assay, utilizing primary rat hippocampal neurons, for identifying small molecule kinase inhibitors that promote neurite growth. The assay is highly reliable and suitable for medium-throughput screening, as indicated by its Z'-factor of 0.73. A focused structurally diverse library of protein kinase inhibitors was screened, revealing several compound groups with the ability to strongly and consistently promote neurite growth. The best performing bioassay hit robustly and consistently promoted axon growth in a postnatal cortical slice culture assay. This study can serve as a jumping-off point for structure activity relationship (SAR) and other drug discovery approaches toward the development of drugs for treating SCI and related neurological pathologies.

Pathway analysis of the human brain transcriptome in disease

Pathway analysis is a powerful method for discerning differentially regulated genes and elucidating their biological importance. It allows for the identification of perturbed or aberrantly expressed genes within a biological context from extensive data sets and offers a simplistic approach for interrogating such data sets. With the growing use of microarrays and RNA-Seq, data for genome-wide studies are growing at an alarming rate, and the use of deep sequencing is revealing elements of the genome previously uncharacterised. Through the employment of pathway analysis, mechanisms in complex diseases may be explored and novel causatives found primarily through differentially regulated genes. Further, with the implementation of next generation sequencing, a deeper resolution may be attained, particularly in identification of isoform diversity and SNPs. Here, we look at a broad overview of pathway analysis in the human brain transcriptome and its relevance in teasing out underlying causes of complex diseases. We will outline processes in data gathering and analysis of particular diseases in which these approaches have been successful.

Dyrk kinases regulate phosphorylation of doublecortin, cytoskeletal organization, and neuronal morphology

In a neuronal overexpression screen focused on kinases and phosphatases, one "hit" was the dual specificity tyrosine phosphorylation-regulated kinase (Dyrk4), which increased the number of dendritic branches in hippocampal neurons. Overexpression of various Dyrk family members in primary neurons significantly changed neuronal morphology. Dyrk1A decreased axon growth, Dyrk3 and Dyrk4 increased dendritic branching, and Dyrk2 decreased both axon and dendrite growth and branching. Kinase-deficient mutants revealed that most of these effects depend on kinase activity. Because doublecortin (DCX), a microtubule-binding protein, regulates cytoskeletal dynamics and neuronal morphogenesis, we investigated the possibility that DCX is a target of Dyrks. We found that overexpression of Dyrk2 and Dyrk3, but not Dyrk1A or Dyrk4, can change DCX phosphorylation status. Mutation of a consensus phosphorylation site for Dyrk kinases at Serine 306 (Ser306) in DCX indicated that this is one target site for Dyrk2 and Dyrk3. Overexpression of Dyrk2 restored altered DCX distribution in the growth cones of dendrites and axons, and partially reversed the morphological effects of DCX overexpression; some of these effects were abrogated by mutation of Ser306 to alanine. These studies implicate Dyrks in the regulation of cytoskeletal organization and process outgrowth in neurons, and suggest that DCX is one relevant Dyrk target.

Molecular control of axon growth: insights from comparative gene profiling and high-throughput screening

Axon regeneration in the mammalian adult central nervous system (CNS) is limited by an intrinsically low capacity for axon growth in many CNS neurons. In contrast, embryonic, peripheral, and many nonmammalian neurons are capable of successful regeneration. Numerous studies have compared mammalian CNS neurons to their counterparts in regenerating systems in an effort to identify candidate genes that control regenerative ability. This review summarizes work using this comparative strategy and examines our current understanding of gene function in axon growth, highlighting the emergence of genome-wide expression profiling and high-throughput screening strategies to identify novel regulators of axon growth.