Single-Cell Transcriptomics and In Vitro Lineage Tracing Reveals Differential Susceptibility of Human iPSC-Derived Midbrain Dopaminergic Neurons in a Cellular Model of Parkinson's Disease.

Advances in stem cell technologies open up new avenues for modelling development and diseases. The success of these pursuits, however, relies on the use of cells most relevant to those targeted by the disease of interest, for example, midbrain dopaminergic neurons for Parkinson's disease. In the present study, we report the generation of a human induced pluripotent stem cell (iPSC) line capable of purifying and tracing nascent midbrain dopaminergic progenitors and their differentiated progeny via the expression of a Blue Fluorescent Protein (BFP). This was achieved by CRISPR/Cas9-assisted knock-in of BFP and Cre into the safe harbour locus AAVS1 and an early midbrain dopaminergic lineage marker gene LMX1A, respectively. Immunocytochemical analysis and single-cell RNA sequencing of iPSC-derived neural cultures confirm developmental recapitulation of the human fetal midbrain and high-quality midbrain cells. By modelling Parkinson's disease-related drug toxicity using 1-Methyl-4-phenylpyridinium (MPP+), we showed a preferential reduction of BFP+ cells, a finding demonstrated independently by cell death assays and single-cell transcriptomic analysis of MPP+ treated neural cultures. Together, these results highlight the importance of disease-relevant cell types in stem cell modelling.

Single-Cell Transcriptomics Reveals Conserved Regulatory Networks in Human and Mouse Interneuron Development.

Inhibitory GABAergic interneurons originate in the embryonic medial ganglionic eminence (MGE) and control network activity in the neocortex. Dysfunction of these cells is believed to lead to runaway excitation underlying seizure-based neurological disorders such as epilepsy, autism, and schizophrenia. Despite their importance in heath and disease, our knowledge about the development of this diverse neuronal population remains incomplete. Here we conducted single-cell RNA sequencing (scRNA-seq) of human foetal MGE from 10 to 15 weeks post conception. These MGE tissues are composed of largely cycling progenitors and immature post-mitotic interneurons with characteristic regional marker expression. Analysis of integrated human and mouse MGE data revealed species-conserved transcriptomic profiles and regulatory programs. Moreover, we identified novel candidate transcription regulators for human interneuron differentiation. These findings provide a framework for in vitro modelling of interneuron development and a strategy for potentially enhancing interneuron production from human pluripotent stem cells.

Lipopolysaccharide distinctively alters human microglia transcriptomes to resemble microglia from Alzheimer's disease mouse models.

Alzheimer's disease (AD) is the most common form of dementia, and risk-influencing genetics implicates microglia and neuroimmunity in the pathogenesis of AD. Induced pluripotent stem cell (iPSC)-derived microglia (iPSC-microglia) are increasingly used as a model of AD, but the relevance of historical immune stimuli to model AD is unclear. We performed a detailed cross-comparison over time on the effects of combinatory stimulation of iPSC-microglia, and in particular their relevance to AD. We used single-cell RNA sequencing to measure the transcriptional response of iPSC-microglia after 24 h and 48 h of stimulation with prostaglandin E2 (PGE2) or lipopolysaccharide (LPS)+interferon gamma (IFN-γ), either alone or in combination with ATPγS. We observed a shared core transcriptional response of iPSC-microglia to ATPγS and to LPS+IFN-γ, suggestive of a convergent mechanism of action. Across all conditions, we observed a significant overlap, although directional inconsistency to genes that change their expression levels in human microglia from AD patients. Using a data-led approach, we identify a common axis of transcriptomic change across AD genetic mouse models of microglia and show that only LPS provokes a transcriptional response along this axis in mouse microglia and LPS+IFN-γ in human iPSC-microglia. This article has an associated First Person interview with the first author of the paper.

Deep phenotyping for precision medicine in Parkinson's disease.

A major challenge in medical genomics is to understand why individuals with the same disorder have different clinical symptoms and why those who carry the same mutation may be affected by different disorders. In every complex disorder, identifying the contribution of different genetic and non-genetic risk factors is a key obstacle to understanding disease mechanisms. Genetic studies rely on precise phenotypes and are unable to uncover the genetic contributions to a disorder when phenotypes are imprecise. To address this challenge, deeply phenotyped cohorts have been developed for which detailed, fine-grained data have been collected. These cohorts help us to investigate the underlying biological pathways and risk factors to identify treatment targets, and thus to advance precision medicine. The neurodegenerative disorder Parkinson's disease has a diverse phenotypical presentation and modest heritability, and its underlying disease mechanisms are still being debated. As such, considerable efforts have been made to develop deeply phenotyped cohorts for this disorder. Here, we focus on Parkinson's disease and explore how deep phenotyping can help address the challenges raised by genetic and phenotypic heterogeneity. We also discuss recent methods for data collection and computation, as well as methodological challenges that have to be overcome.

Combining multiomics and drug perturbation profiles to identify muscle-specific treatments for spinal muscular atrophy.

Spinal muscular atrophy (SMA) is a neuromuscular disorder caused by loss of survival motor neuron (SMN) protein. While SMN restoration therapies are beneficial, they are not a cure. We aimed to identify potentially novel treatments to alleviate muscle pathology combining transcriptomics, proteomics, and perturbational data sets. This revealed potential drug candidates for repurposing in SMA. One of the candidates, harmine, was further investigated in cell and animal models, improving multiple disease phenotypes, including lifespan, weight, and key molecular networks in skeletal muscle. Our work highlights the potential of multiple and parallel data-driven approaches for the development of potentially novel treatments for use in combination with SMN restoration therapies.

Emergence of co-expression in gene regulatory networks.

Transcriptomes are known to organize themselves into gene co-expression clusters or modules where groups of genes display distinct patterns of coordinated or synchronous expression across independent biological samples. The functional significance of these co-expression clusters is suggested by the fact that highly coexpressed groups of genes tend to be enriched in genes involved in common functions and biological processes. While gene co-expression is widely assumed to reflect close regulatory proximity, the validity of this assumption remains unclear. Here we use a simple synthetic gene regulatory network (GRN) model and contrast the resulting co-expression structure produced by these networks with their known regulatory architecture and with the co-expression structure measured in available human expression data. Using randomization tests, we found that the levels of co-expression observed in simulated expression data were, just as with empirical data, significantly higher than expected by chance. When examining the source of correlated expression, we found that individual regulators, both in simulated and experimental data, fail, on average, to display correlated expression with their immediate targets. However, highly correlated gene pairs tend to share at least one common regulator, while most gene pairs sharing common regulators do not necessarily display correlated expression. Our results demonstrate that widespread co-expression naturally emerges in regulatory networks, and that it is a reliable and direct indicator of active co-regulation in a given cellular context.

TREM2 Alzheimer's variant R47H causes similar transcriptional dysregulation to knockout, yet only subtle functional phenotypes in human iPSC-derived macrophages.

BACKGROUND: TREM2 is a microglial cell surface receptor, with risk mutations linked to Alzheimer's disease (AD), including R47H. TREM2 signalling via SYK aids phagocytosis, chemotaxis, survival, and changes to microglial activation state. In AD mouse models, knockout (KO) of TREM2 impairs microglial clustering around amyloid and prevents microglial activation. The R47H mutation is proposed to reduce TREM2 ligand binding. We investigated cell phenotypes of the R47H mutant and TREM2 KO in a model of human microglia, and compared their transcriptional signatures, to determine the mechanism by which R47H TREM2 disrupts function. METHODS: We generated human microglia-like iPSC-macrophages (pMac) from isogenic induced pluripotent stem cell (iPSC) lines, with homozygous R47H mutation or TREM2 knockout (KO). We firstly validated the effect of the R47H mutant on TREM2 surface and subcellular localization in pMac. To assess microglial phenotypic function, we measured phagocytosis of dead neurons, cell morphology, directed migration, survival, and LPS-induced inflammation. We performed bulk RNA-seq, comparing significant differentially expressed genes (DEGs; p < 0.05) between the R47H and KO versus WT, and bioinformatically predicted potential upstream regulators of TREM2-mediated gene expression. RESULTS: R47H modified surface expression and shedding of TREM2, but did not impair TREM2-mediated signalling, or gross phenotypes that were dysregulated in the TREM2 KO (phagocytosis, motility, survival). However, altered gene expression in the R47H TREM2 pMac overlapped by 90% with the TREM2 KO and was characterised by dysregulation of genes involved with immunity, proliferation, activation, chemotaxis, and adhesion. Downregulated mediators of ECM adhesion included the vitronectin receptor αVß3, and consequently, R47H TREM2 pMac adhered weakly to vitronectin compared with WT pMac. To counteract these transcriptional defects, we investigated TGFß1, as a candidate upstream regulator. TGFß1 failed to rescue vitronectin adhesion of pMac, although it improved αVß3 expression. CONCLUSIONS: The R47H mutation is not sufficient to cause gross phenotypic defects of human pMac under standard culture conditions. However, overlapping transcriptional defects with TREM2 KO supports the hypothesised partial loss-of-function effects of the R47H mutation. Furthermore, transcriptomics can guide us to more subtle phenotypic defects in the R47H cells, such as reduced cell adhesion, and can be used to predict targets for therapeutic intervention.

A single-cell atlas of the human substantia nigra reveals cell-specific pathways associated with neurological disorders.

We describe a human single-nuclei transcriptomic atlas for the substantia nigra (SN), generated by sequencing approximately 17,000 nuclei from matched cortical and SN samples. We show that the common genetic risk for Parkinson's disease (PD) is associated with dopaminergic neuron (DaN)-specific gene expression, including mitochondrial functioning, protein folding and ubiquitination pathways. We identify a distinct cell type association between PD risk and oligodendrocyte-specific gene expression. Unlike Alzheimer's disease (AD), we find no association between PD risk and microglia or astrocytes, suggesting that neuroinflammation plays a less causal role in PD than AD. Beyond PD, we find associations between SN DaNs and GABAergic neuron gene expression and multiple neuropsychiatric disorders. Conditional analysis reveals that distinct neuropsychiatric disorders associate with distinct sets of neuron-specific genes but converge onto shared loci within oligodendrocytes and oligodendrocyte precursors. This atlas guides our aetiological understanding by associating SN cell type expression profiles with specific disease risk.

Pro-maturational Effects of Human iPSC-Derived Cortical Astrocytes upon iPSC-Derived Cortical Neurons.

Astrocytes influence neuronal maturation and function by providing trophic support, regulating the extracellular environment, and modulating signaling at synapses. The emergence of induced pluripotent stem cell (iPSC) technology offers a human system with which to validate and re-evaluate insights from animal studies. Here, we set out to examine interactions between human astrocytes and neurons derived from a common cortical progenitor pool, thereby recapitulating aspects of in vivo cortical development. We show that the cortical iPSC-derived astrocytes exhibit many of the molecular and functional hallmarks of astrocytes. Furthermore, optogenetic and electrophysiological co-culture experiments reveal that the iPSC-astrocytes can actively modulate ongoing synaptic transmission and exert pro-maturational effects upon developing networks of iPSC-derived cortical neurons. Finally, transcriptomic analyses implicate synapse-associated extracellular signaling in the astrocytes' pro-maturational effects upon the iPSC-derived neurons. This work helps lay the foundation for future investigations into astrocyte-to-neuron interactions in human health and disease.

Human-Specific Transcriptome of Ventral and Dorsal Midbrain Dopamine Neurons.

OBJECTIVE: Neuronal loss in the substantia nigra pars compacta (SNpc) in Parkinson disease (PD) is not uniform, as dopamine neurons from the ventral tier are lost more rapidly than those of the dorsal tier. Identifying the intrinsic differences that account for this differential vulnerability may provide a key for developing new treatments for PD. METHODS: Here, we compared the RNA-sequenced transcriptomes of ~100 laser captured microdissected SNpc neurons from each tier from 7 healthy controls. RESULTS: Expression levels of dopaminergic markers were similar across the tiers, whereas markers specific to the neighboring ventral tegmental area were virtually undetected. After accounting for unwanted sources of variation, we identified 106 differentially expressed genes (DEGs) between the SNpc tiers. The genes higher in the dorsal/resistant SNpc tier neurons displayed coordinated patterns of expression across the human brain, their protein products had more interactions than expected by chance, and they demonstrated evidence of functional convergence. No significant shared functionality was found for genes higher in the ventral/vulnerable SNpc tier. Surprisingly but importantly, none of the identified DEGs was among the familial PD genes or genome-wide associated loci. Finally, we found some DEGs in opposite tier orientation between human and analogous mouse populations. INTERPRETATION: Our results highlight functional enrichments of vesicular trafficking, ion transport/homeostasis and oxidative stress genes showing higher expression in the resistant neurons of the SNpc dorsal tier. Furthermore, the comparison of gene expression variation in human and mouse SNpc populations strongly argues for the need of human-focused omics studies. ANN NEUROL 2020;87:853-868.

Transcriptional, Behavioral and Biochemical Profiling in the 3xTg-AD Mouse Model Reveals a Specific Signature of Amyloid Deposition and Functional Decline in Alzheimer's Disease.

Alzheimer's disease (AD)-related degenerative decline is associated to the presence of amyloid beta (Aß) plaque lesions and neuro fibrillary tangles (NFT). However, the precise molecular mechanisms linking Aß deposition and neurological decline are still unclear. Here we combine genome-wide transcriptional profiling of the insular cortex of 3xTg-AD mice and control littermates from early through to late adulthood (2-14 months of age), with behavioral and biochemical profiling in the same animals to identify transcriptional determinants of functional decline specifically associated to build-up of Aß deposits. Differential expression analysis revealed differentially expressed genes (DEGs) in the cortex long before observed onset of behavioral symptoms in this model. Using behavioral and biochemical data derived from the same mice and samples, we found that down but not up-regulated DEGs show a stronger average association with learning performance than random background genes in control not seen in AD mice. Conversely, these same genes were found to have a stronger association with Aß deposition than background genes in AD but not in control mice, thereby identifying these genes as potential intermediaries between abnormal Aß/NFT deposition and functional decline. Using a complementary approach, gene ontology analysis revealed a highly significant enrichment of learning and memory, associative, memory, and cognitive functions only among down-regulated, but not up-regulated, DEGs. Our results demonstrate wider transcriptional changes triggered by the abnormal deposition of Aß/NFT occurring well before behavioral decline and identify a distinct set of genes specifically associated to abnormal Aß protein deposition and cognitive decline.

CSF1R inhibitor JNJ-40346527 attenuates microglial proliferation and neurodegeneration in P301S mice.

Neuroinflammation and microglial activation are significant processes in Alzheimer's disease pathology. Recent genome-wide association studies have highlighted multiple immune-related genes in association with Alzheimer's disease, and experimental data have demonstrated microglial proliferation as a significant component of the neuropathology. In this study, we tested the efficacy of the selective CSF1R inhibitor JNJ-40346527 (JNJ-527) in the P301S mouse tauopathy model. We first demonstrated the anti-proliferative effects of JNJ-527 on microglia in the ME7 prion model, and its impact on the inflammatory profile, and provided potential CNS biomarkers for clinical investigation with the compound, including pharmacokinetic/pharmacodynamics and efficacy assessment by TSPO autoradiography and CSF proteomics. Then, we showed for the first time that blockade of microglial proliferation and modification of microglial phenotype leads to an attenuation of tau-induced neurodegeneration and results in functional improvement in P301S mice. Overall, this work strongly supports the potential for inhibition of CSF1R as a target for the treatment of Alzheimer's disease and other tau-mediated neurodegenerative diseases.

Postmitotic cell longevity-associated genes: a transcriptional signature of postmitotic maintenance in neural tissues.

Different cell types have different postmitotic maintenance requirements. Nerve cells, however, are unique in this respect as they need to survive and preserve their functional complexity for the entire lifetime of the organism, and failure at any level of their supporting mechanisms leads to a wide range of neurodegenerative conditions. Whether these differences across tissues arise from the activation of distinct cell type-specific maintenance mechanisms or the differential activation of a common molecular repertoire is not known. To identify the transcriptional signature of postmitotic cellular longevity (PMCL), we compared whole-genome transcriptome data from human tissues ranging in longevity from 120 days to over 70 years and found a set of 81 genes whose expression levels are closely associated with increased cell longevity. Using expression data from 10 independent sources, we found that these genes are more highly coexpressed in longer-living tissues and are enriched in specific biological processes and transcription factor targets compared with randomly selected gene samples. Crucially, we found that PMCL-associated genes are downregulated in the cerebral cortex and substantia nigra of patients with Alzheimer's and Parkinson's disease, respectively, as well as Hutchinson-Gilford progeria-derived fibroblasts, and that this downregulation is specifically linked to their underlying association with cellular longevity. Moreover, we found that sexually dimorphic brain expression of PMCL-associated genes reflects sexual differences in lifespan in humans and macaques. Taken together, our results suggest that PMCL-associated genes are part of a generalized machinery of postmitotic maintenance and functional stability in both neural and non-neural cells and support the notion of a common molecular repertoire differentially engaged in different cell types with different survival requirements.

Neocortex expansion is linked to size variations in gene families with chemotaxis, cell-cell signalling and immune response functions in mammals.

Increased brain size is thought to have played an important role in the evolution of mammals and is a highly variable trait across lineages. Variations in brain size are closely linked to corresponding variations in the size of the neocortex, a distinct mammalian evolutionary innovation. The genomic features that explain and/or accompany variations in the relative size of the neocortex remain unknown. By comparing the genomes of 28 mammalian species, we show that neocortical expansion relative to the rest of the brain is associated with variations in gene family size (GFS) of gene families that are significantly enriched in biological functions associated with chemotaxis, cell-cell signalling and immune response. Importantly, we find that previously reported GFS variations associated with increased brain size are largely accounted for by the stronger link between neocortex expansion and variations in the size of gene families. Moreover, genes within these families are more prominently expressed in the human neocortex during early compared with adult development. These results suggest that changes in GFS underlie morphological adaptations during brain evolution in mammalian lineages.

Modular reorganization of the global network of gene regulatory interactions during perinatal human brain development.

BACKGROUND: During early development of the nervous system, gene expression patterns are known to vary widely depending on the specific developmental trajectories of different structures. Observable changes in gene expression profiles throughout development are determined by an underlying network of precise regulatory interactions between individual genes. Elucidating the organizing principles that shape this gene regulatory network is one of the central goals of developmental biology. Whether the developmental programme is the result of a dynamic driven by a fixed architecture of regulatory interactions, or alternatively, the result of waves of regulatory reorganization is not known. RESULTS: Here we contrast these two alternative models by examining existing expression data derived from the developing human brain in prenatal and postnatal stages. We reveal a sharp change in gene expression profiles at birth across brain areas. This sharp division between foetal and postnatal profiles is not the result of pronounced changes in level of expression of existing gene networks. Instead we demonstrate that the perinatal transition is marked by the widespread regulatory rearrangement within and across existing gene clusters, leading to the emergence of new functional groups. This rearrangement is itself organized into discrete blocks of genes, each targeted by a distinct set of transcriptional regulators and associated to specific biological functions. CONCLUSIONS: Our results provide evidence of an acute modular reorganization of the regulatory architecture of the brain transcriptome occurring at birth, reflecting the reassembly of new functional associations required for the normal transition from prenatal to postnatal brain development.

Optimization of next-generation sequencing transcriptome annotation for species lacking sequenced genomes.

Next-generation sequencing methods, such as RNA-seq, have permitted the exploration of gene expression in a range of organisms which have been studied in ecological contexts but lack a sequenced genome. However, the efficacy and accuracy of RNA-seq annotation methods using reference genomes from related species have yet to be robustly characterized. Here we conduct a comprehensive power analysis employing RNA-seq data from Drosophila melanogaster in conjunction with 11 additional genomes from related Drosophila species to compare annotation methods and quantify the impact of evolutionary divergence between transcriptome and the reference genome. Our analyses demonstrate that, regardless of the level of sequence divergence, direct genome mapping (DGM), where transcript short reads are aligned directly to the reference genome, significantly outperforms the widely used de novo and guided assembly-based methods in both the quantity and accuracy of gene detection. Our analysis also reveals that DGM recovers a more representative profile of Gene Ontology functional categories, which are often used to interpret emergent patterns in genomewide expression analyses. Lastly, analysis of available primate RNA-seq data demonstrates the applicability of our observations across diverse taxa. Our quantification of annotation accuracy and reduced gene detection associated with sequence divergence thus provides empirically derived guidelines for the design of future gene expression studies in species without sequenced genomes.

Nociceptin/Orphanin FQ Inhibits the Survival and Axon Growth of Midbrain Dopaminergic Neurons Through a p38-MAPK Dependent Mechanism.

Nociceptin/orphanin FQ (N/OFQ) is an opioid-like neuropeptide that binds and signals through a G-protein-coupled receptor called the N/OFQ peptide (NOP) receptor. N/OFQ and the NOP receptor are expressed in the midbrain and have been implicated in the pathogenesis of Parkinson's disease (PD). Genetic removal of the N/OFQ precursor partially protects midbrain dopaminergic neurons from 1-methyl-4-phenylpyridine-induced toxicity, suggesting that endogenous N/OFQ may be detrimental to dopaminergic neurons. However, whether N/OFQ directly affects the survival and growth of dopaminergic neurons is unknown. Here, we show that N/OFQ has a detrimental effect on the survival of dopaminergic neurons and the growth of their axons in primary cultures of the E14 rat ventral mesencephalon. N/OFQ potentiates the effects of the neurotoxins 6-hydroxydopamine and 1-methyl-4-phenylpyridinium through p38-MAPK signalling. We also show that like α-synuclein, there is a significant reduction in N/OFQ messenger RNA (mRNA) expression in the midbrain of patients with Parkinson's disease. These results demonstrate for the first time that N/OFQ is detrimental to the survival and growth of dopaminergic neurons and that its expression is altered in the midbrain of patients with Parkinson's disease.

Modular and coordinated expression of immune system regulatory and signaling components in the developing and adult nervous system.

During development, the nervous system (NS) is assembled and sculpted through a concerted series of neurodevelopmental events orchestrated by a complex genetic programme. While neural-specific gene expression plays a critical part in this process, in recent years, a number of immune-related signaling and regulatory components have also been shown to play key physiological roles in the developing and adult NS. While the involvement of individual immune-related signaling components in neural functions may reflect their ubiquitous character, it may also reflect a much wider, as yet undescribed, genetic network of immune-related molecules acting as an intrinsic component of the neural-specific regulatory machinery that ultimately shapes the NS. In order to gain insights into the scale and wider functional organization of immune-related genetic networks in the NS, we examined the large scale pattern of expression of these genes in the brain. Our results show a highly significant correlated expression and transcriptional clustering among immune-related genes in the developing and adult brain, and this correlation was the highest in the brain when compared to muscle, liver, kidney and endothelial cells. We experimentally tested the regulatory clustering of immune system (IS) genes by using microarray expression profiling in cultures of dissociated neurons stimulated with the pro-inflammatory cytokine TNF-alpha, and found a highly significant enrichment of immune system-related genes among the resulting differentially expressed genes. Our findings strongly suggest a coherent recruitment of entire immune-related genetic regulatory modules by the neural-specific genetic programme that shapes the NS.

Increased brain size in mammals is associated with size variations in gene families with cell signalling, chemotaxis and immune-related functions.

Genomic determinants underlying increased encephalization across mammalian lineages are unknown. Whole genome comparisons have revealed large and frequent changes in the size of gene families, and it has been proposed that these variations could play a major role in shaping morphological and physiological differences among species. Using a genome-wide comparative approach, we examined changes in gene family size (GFS) and degree of encephalization in 39 fully sequenced mammalian species and found a significant over-representation of GFS variations in line with increased encephalization in mammals. We found that this relationship is not accounted for by known correlates of brain size such as maximum lifespan or body size and is not explained by phylogenetic relatedness. Genes involved in chemotaxis, immune regulation and cell signalling-related functions are significantly over-represented among those gene families most highly correlated with encephalization. Genes within these families are prominently expressed in the human brain, particularly the cortex, and organized in co-expression modules that display distinct temporal patterns of expression in the developing cortex. Our results suggest that changes in GFS associated with encephalization represent an evolutionary response to the specific functional requirements underlying increased brain size in mammals.