Next Article in Journal
Dynamics of the Ligand Excited States Relaxation in Novel β-Diketonates of Non-Luminescent Trivalent Metal Ions
Next Article in Special Issue
The Association of Hippocampal Long-Term Potentiation-Induced Gene Expression with Genetic Risk for Psychosis
Previous Article in Journal
Reference Values on Children’s Hair for 28 Elements (Heavy Metals and Essential Elements) Based on a Pilot Study in a Representative Non-Contaminated Local Area
Previous Article in Special Issue
Extracellular Matrix Regulation in Physiology and in Brain Disease
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

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

1
Neuroscience and Mental Health Research Institute, School of Medicine, Cardiff University, Cardiff CF24 4HQ, UK
2
UK Dementia Research Institute Cardiff, Cardiff University, Cardiff CF24 4HQ, UK
3
Division of Neuroscience, School of Bioscience, Cardiff University, Cardiff CF10 3AX, UK
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Int. J. Mol. Sci. 2023, 24(9), 8122; https://doi.org/10.3390/ijms24098122
Submission received: 28 March 2023 / Revised: 20 April 2023 / Accepted: 27 April 2023 / Published: 1 May 2023
(This article belongs to the Special Issue Gene Regulation in Brain Development and Physiology)

Abstract

:
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.

1. Introduction

The gamma-aminobutyric acid containing (GABAergic) inhibitory interneurons represent approximately 20% of the entire neuronal population in the human neocortex. Despite the small proportion, it is apparent that interneurons play a pivotal role in most, if not all, cortical functions [1,2,3]. Interneurons are required for maintaining the appropriate excitation–inhibition balance within the neocortex while the connectivity of a single interneuron onto multiple post-synaptic targets allows a synchronised, rhythmic pattern of firing [2,3,4]. The synchronised firing pattern generates network oscillations necessary for stable, long-distance transmission [5,6]. Unsurprisingly, loss of this homeostasis has been implicated in a number of neurological disorders, including autism spectrum disorders, schizophrenia, and epilepsies [2,3,4,5,7,8].
Interneurons are primarily born in two neighbouring structures of the developing ventral telencephalon, namely the medial and caudal ganglionic eminences (MGE and CGE), respectively [9,10,11]. They consist of three major subtypes based on their expression of somatostatin (SST), parvalbumin (PV), and ionotropic serotonin receptor 5HT3a (5HT3aR), respectively, and are highly heterogeneous in morphology and function [12,13]. Of these, the MGE-derived PV interneurons appear to be preferentially affected in neurodevelopmental disorders [1,14,15,16,17,18,19]. Thus, there has been growing interest in generating cortical interneurons from human pluripotent stem cells for disease modelling and drug screening [20,21,22,23,24,25]. These efforts are facilitated by recent advances in single-cell transcriptomics analysis of the developing and adult interneurons which have begun the identification of transcription factors contributing to cortical interneuron subtype fate determination in mice [26,27,28]. However, our understanding about the transcriptional program controlling interneuron differentiation remains limited.
Here we present a single-cell transcriptome analysis for human embryonic MGE of post conceptional weeks (PCWs) 10 to 15. Our finding provides support for the emergence of interneuron subtype diversification during human MGE development. Moreover, integrating mouse and human MGE transcriptional profiles identified a cross-species concordance in MGE development and identified potential novel regulators for interneuron differentiation.

2. Results

2.1. Developmental Dynamics of Human Fetal MGE

We employed the Fluidigm C1 platform and SMART-Seq2 technology to study the transcriptomic profile of individual cells of human foetal MGE (Figure 1A). Of the 358 qualifying cells corresponding to six embryos of 10 to 15 PCWs, we obtained a mean of ~2 million counts per cell and detected a median of 6208 protein-coding genes per cell (Figure S1A–C). To visualise the dynamics of these cells, we performed principal component analysis using the top 1000 most variable genes (see Methods). Uniform Manifold Approximation and Projection (UMAP) for dimensionality reduction revealed a weak segregation amongst cells of different embryonic stages, suggesting MGE cells of 10–15 PCWs are transcriptionally similar with each other (Figure 1B). However, there was a clear segregation by cell cycle phases with subtle differences in the proportion of cells at each phase between the 10 and 15 PCW samples (Figure 1C and Figure S1D).
We then analysed the gene expression using an alternative workflow of Seurat to remove the effect of cell cycle-dependent genes on clustering and identified three main clusters between the four developmental stages of foetal MGE cells (Figure 1D). Expression of neural progenitor mark genes was enriched in cluster 0 and 2, while neuronal gene markers were preferentially detected in cluster 1 (Figure 1E). This differential gene expression pattern is consistent with the cell cycle status where cluster 0 and cluster 2 cells are mostly in S or G2/M phases while cluster 1 is constituted mainly by cells in G0/1 phases (Figure S1E). Between the two progenitor populations, cluster 0 was marked by genes involved in DNA replication (including Helicase, HELLS; Proliferating Cell Nuclear Antigen, PCNA; and Ran-specific GTPase-activating protein, RANBP1), while cell division and cell cycle-related genes, such as Karyopherin Subunit Alpha 2 (KPNA2) and Ubiquitin Conjugating Enzyme E2 C (UBE2C), characterised cluster 1 (Figure 1F and Figure S2). In contrast, genes involved in terminal neuronal differentiation and synaptogenesis such as Doublecortin (DCX), Stathmin 2 (STMN2), and Synaptotagmin 1 (SYT1) were enriched in cluster 1 cells, consistent with cell differentiation and nervous system development being the most significantly enriched biological process of this cluster (Figure 1F and Figure S2).

2.2. Transcriptional Heterogeneity and Emergence of Interneuron Subtypes

Transcription factor NKX2.1 is specifically expressed in MGE progenitors and post-mitotic MGE-derived interneurons that reside in the striatum. NKX2.1+ cells were highly enriched in samples of all ages and were more abundant in clusters 0 and 2 (Figure 2A). Consistent with the MGE identity of our cells, transcripts that are highly specific to LGE and CGE or cortex, such as ZNF503 and PROX1, respectively, were rarely detected (Figure 2A). Moreover, expression of other MGE and interneuron developmental regulators correlated well with differentiation and cell cycle status. For example, FOXG1+, DLX2+ and DLX5+ cells were distributed across 10–15 PCW MGE samples while SOX6, GAD1 and GAD2 were preferentially expressed in the neuronal cluster (Figure 2A).
In mice, Nkx2.1+ progenitors give rise to cholinergic interneurons (ChAT+) under the control of Lhx8 and GABAergic interneurons marked by Lhx6, which later differentiate into parvalbumin (PVALB) and somatostatin (SST) subtypes. maf and mef2c have been suggested to have preferential regulatory roles for SST and PV fates, respectively [27,28]. MEF2C and MAF were detected in a proportion of LHX6+ cells mostly in the neuronal cluster where SST+ and neuropeptide Y (NPY+) cells were also found (Figure 2B). On the other hand, cells expressing LHX8 and other cholinergic regulators (eg. ISL1, ZIC1) were detected albeit in fewer number of cells compared to LHX6 (Figure 2B) [29,30]. However, ChAT+ neurons were not detected. These observations either reflect a later onset of cholinergic development or the relative representation of the LHX6 and LHX8 interneuron lineages.
Repression of Nkx2.1 by Zeb2 in mouse MGE progenitors confers cortical interneuron fate, while in its absence these progenitors progress into striatal interneurons retaining Nkx2.1 expression [31]. ZEB2 exhibits a heterogeneous expression pattern across the clusters while a mosaic expression of striatal interneuron indicator genes such as ETV1 and NFIA was also observed. Together, these findings indicate the emergence of cortical interneuron subtypes and spatial specification.
The single-cell gene expression presented above is supported by immunostaining of class-defining markers in dissociated 15 PCW MGE tissue (Figure S3). After 2 days in culture, a significant proportion of the population was found to express progenitor markers FOXG1 and NKX2.1 (42 ± 15% and 25 ± 4.1%, respectively). Cells marked by post-mitotic interneuron subtype markers, such as SST (0.4% ± 0.1), ChAT (1.1 ± 0.08%) and Calretinin (encoded by CALB2, 18 ± 3%), were also readily detected along with cells stained positive for pan-neuronal marker MAP2 and NeuN (25 ± 4.5%, 17 ± 3.3%, respectively). Consistent with scRNA-seq data, no PVALB-positive cells were identified.

2.3. Conserved Biological Processes between the Human and Mouse MGE

We next integrated the human MGE scRNAseq data with that previously published from the mouse MGE based on the shared one-to-one ortholog (see Methods) [28]. The integrated dataset consisted of 1614 cells and 10,362 genes from the MGE across two species. The mouse MGE cells from embryonic day 12.5 (E12.5) and E14.5 presented the largest variance and clustered separately. The 10–15 PCW human MGE cells were, however, more evenly distributed (Figure 3A). A strong segregation of cycling and non-cycling cells based on their cell cycle phase is consistent between the two species (Figure 3B) and is in agreement with differential expression of NES and MAP2 (Figure 3C).
Seven cell clusters were identified in the integrated human and mouse MGE dataset with a variable number of cells from each species (Figure 3D). We observed a good correspondence between the clusters previously identified in the human MGE alone (Figure 1D) and those in the integrated dataset (Figure 3E). For example, human cluster 1 mostly corresponds to the integrated cluster 0 (i-cluster)/post-mitotic neurons. Using the FindMarkers function of Seurat, we identified conserved genes between the human and mouse in a given cell cluster (Table S1), all of which were detected in at least 25% of cells. Gene ontology (GO) enrichment analysis on the conserved genes of i-cluster 0 highlighted genes involved in dendrite morphogenesis, regulation of neuron projection development, axonogenesis, cell adhesion and neuron migration (Figure S4A). In turn, integrated clusters 2 and 3 are marked by the enrichment of genes involved in DNA replication (Figure S4B,C). Mitochondrial/metabolic-associated genes were heavily enriched in cluster 4 (Figure S4D). Finally, species integrated cluster 5 contained cells highly expressing mitotic cell cycle-associated genes (Figure S4E). There were no significantly enriched biological processes among the 42 and 41 conserved genes of the integrated clusters 1 and 6 when using the applied threshold.

2.4. Gene Regulatory Network Inference with SCENIC Reveals Novel Human Candidate Interneuron Regulators

We then employed the SCENIC (Single-Cell rEgulatory Network Inference) tool to construct gene regulatory networks (or regulons) consisting of a set of transcription factors and their target genes that are most likely active within the cell clusters defined in the integrated mouse and human MGE [32]. This revealed a network of SMAD3 containing regulons in progenitor cluster 1 (Figure 4A). Interestingly, TGFβ signalling has been shown previously to be required for regulating NKX2.1+ progenitor cell cycle exit and their differentiation into LHX6+ interneurons [25], hence providing direct experimental validation for this prediction.
Moreover, SCENIC identified a neuronal dominant cluster (i-cluster 0) containing regulons that correspond to 11 transcription factors and their target genes (Figure 4A). These include known interneuron specification genes, such as the SRY-Box Transcription Factor 6 (SOX6), LIM Homeobox 6 (LHX6), Distal-Less Homeobox 5 (DLX5) and Visual System Homeobox 1 (VSX1), as well as transcription factors known to be important in neurogenesis, such as SRY-Box Transcription Factor 4/11 (SOX4 and SOX11). In addition, we identified transcription factors B-Cell Lymphoma/Leukaemia 11A (BCL11A/CTIP1), AT-Rich Interaction Domain 5B (ARID5B), Zinc Finger Homeobox 3 (ZFHX3) and Kruppel-Like Factor 7 (KLF7) as novel candidate regulators for interneuron differentiation and/or specification. The inferred gene regulatory network containing only high confidence links for the conserved gene markers of cluster 0 and their immediate upstream regulators is shown in Figure 4B. When compared to the inferred regulatory network based on mouse data only, we also identified Sox6, Lhx6 and Klf7 as active regulators on the cells corresponding to the neuronal integrated cluster 0, indicating conserved regulatory elements in interneuron specification (Figure S5).

3. Discussion

An in-depth understanding of human cortical and striatal interneuron development is fundamental to decipher their function in the normal brain and how aberrant development and dysfunction of these cells may lead to neurodevelopmental disorders. This single-cell RNAseq study illustrates the heterogeneity and dynamic transcriptomic landscape of human foetal MGE cells and provides evidence on the conserved developmental trajectory of human and mouse MGE interneurons. Moreover, our study identified a novel set of transcription factors as potential regulatory networks governing human MGE interneuron differentiation.
With several competing theories in the literature, details concerning interneuron diversification into specific subtypes remain to be conclusively determined [33,34]. With the caveat of moderate cell numbers that limit prospective subtype clustering, cells expressing prospective MGE-derived interneuron subtype regulators or predictive genes were detected in our MGE samples of all ages. Moreover, some subtype-committed cell populations (e.g., NPY and SST) were detected. Thus, data described here support the notion that, as in mice, lineage-dictated subpopulations emerge within the MGE [27,28,35,36]. However, we did not find PVALB+ and ChAT+ neurons in human MGE of 10–15 PCWs. This may be a reflection of the developmental period studied or due to them existing in very low numbers that we were not able to detect in our current analysis [34,37]. Alternatively, these subtypes may emerge primarily within the region of destination under the influence of local extrinsic cues [34,38].
Using the SCENIC tool on the integrated human and mouse MGE single-cell profiles, we identified several sets of transcription factors characteristic to defined i-Clusters. Albeit with a moderate regulon specificity score (RSS), implicating low cell type specificity, a SMAD3-containing regulatory network was predicted for neural progenitor i-Cluster 1. This cluster of cells is enriched with a DLX2-expressing progenitor population with another SMAD (SMAD1) as one of its signature genes (Table S1). Interestingly, TGFβ signalling, which acts through SMAD proteins, has been demonstrated as a novel regulatory pathway in modulating the maintenance of human ESC-derived NKX2.1+ progenitors and their differentiation into LHX6+ interneurons [25]. This prior study thus validates the SCENIC prediction. On the other hand, the transcription factors identified in this set provide a foundation for further investigation into detailed molecular mechanisms downstream of TGFβ control of interneuron progenitor expansion and cell cycle exit.
Moreover, the regulatory network inference identified a set of species-conserved developmental regulators such as LHX6, SOX6 and DLX5 being active in the neuronal dominant i-Cluster 9 (human cluster 1). These transcription factors have previously been applied successfully to generate interneurons from human pluripotent stem cells via direct programming [39]. Significantly, the regulons contain transcription factors not previously implicated in interneuron development, such as KLF7, BCL11A and ZFHX3. KLF7 has been reported to be capable of promoting axon outgrowth and olfactory sensory neurogenesis while BCL11A is required for spinal cord sensory circuit formation [40,41,42]. ZFHX3 encodes a transcription factor that has been shown to be capable of promoting neurite outgrowth [43]. Interestingly, de novo loss-of-function mutation of this gene has been suggested to cause intellectual disability [44]. The identification of ZFHX3 as a candidate interneuron regulator thus provides further support that abnormal interneuron development contributes to the pathogenesis of neural developmental disorders. Moreover, our findings present an opportunity to generate a more authentic subtype enriched with in vitro interneuron differentiation culture through the engineered expression of one or more of these transcription factors.

4. Materials and Methods

4.1. Tissue Preparation

Human foetal tissue was collected from medical terminations of pregnancy with full donor consent, through the SWIFT human foetal tissue bank (http://www.biobankswales.org.uk/swift-research-tissue-bank), under UK Human Tissue Authority research license (No. 12457) held by Cardiff University and with ethical approval of the project from the Bro Taf local research ethics committee. Gestational age was estimated through ultrasound scan prior to the procedure in combination with measurement of foetal regions [45]. MGE tissue was harvested from seven foetuses corresponding to approximately 10 and 12–15 post conception weeks (PCWs). MGEs were individually incubated at 37 °C for 20 min in TrypLE Express containing Donase-α. Tissue was dissociated by trituration in 200 μL of DMEM/D to obtain a quasi-single-cell suspension. The cell concentration was calculated by manual counting and viability was measured by Trypan Blue assay.

4.2. Single-Cell Capture, Library Preparation and Sequencing

Single cells were captured according to the manufacturer’s protocol on a C1 Single-Cell IFC microfluidic chip designed for 10 to 17 μm cells (Fluidigm, Inc., South San Francisco, CA, USA). This captured one cell in each of up to 96 capture chambers. Capture efficiency and cell viability were determined by microscopy before cDNA amplification and library preparation using the SMART-Seq 46 v4 Ultra Low Input RNA kit chemistry (Clontech, Mountain View, CA, USA) and Nextera XT DNA Library Prep Index Kit (Illumina, San Diego, CA, USA). Single-cell sequencing was carried out on the HiSeq 4000 Illumina platform with the following parameters: paired-end reads, read length 2 × 75 bp and read depth of approximately 2 million reads per cell.

4.3. Data Processing

The quality of the sequences was assessed using FastQC (version 0.11.8) and results were summarised using MultiQC. After removing adapters by TrimGalore (version 0.5.0), transcript abundances were quantified using Kallisto (version 0.46.0) [46]. The reference index used include all cDNA and ncRNA sequences from Ensembl (release 99) as well as the Array Control RNA spike sequences from Thermo Fisher. Sequences in scaffold chromosomes were removed and we kept only chromosomes 1–22, X and Y. We removed cells with less than 25% of reads mapping to protein-coding genes and retained cells with at least half a million reads. We further removed an outlier sample from 12 PCWs. We filtered out lowly expressed genes, keeping only those with more than 5 counts in at least 1% of the cells. The filtered dataset consisted of 14,673 protein-coding genes across 358 cells.

4.4. Dimensionality Reduction

Gene expression data were logNormalized using Seurat function and a scale factor equal to 10,000. Principal component analysis was based on the top 1000 most variable features (vst method). We identify cell clusters based on the first 10 principal components using the FindCluster function with a resolution equal to one in first instance and equal to 0.5 after regressing out the difference in cell cycle scores. Similarly, Uniform Manifold Approximation and Projection (UMAP) was based on the first 10 principal components. We followed the alternative workflow to keep differences between cycling and non-cycling cells by regressing out the differences between the S and the G2M scores, based on the updated (2019) cell cycle gene annotations from Seurat.

4.5. Differential Expression Analysis

We used the FindMarkers function from Seurat to obtain gene markers for each of the four human MGE clusters, including only positive markers detected in at least 25% of cells and with a minimum logFC of 0.25. For further analysis, we only considered gene markers with an adjusted p value < 0.05.

4.6. Gene Ontology Enrichment Analysis

GO annotations were retrieved from Ensembl through biomaRt, and we segregated annotations according to three main branching categories: biological processes, cellular component and molecular function. We only focused on GO terms with at least 20 genes in our gene background population (n = 14,673), and to assess GO enrichment we used a hypergeometric test. We adjusted for multiple testing and only GO terms with a false discovery rate < 0.05 were considered enriched.

4.7. Integration with Mouse MGE Dataset

Gene expression data from the mouse MGE were obtained from the GEO record (GSE109796), and cell annotations were kindly provided by Zhen Li and further analysis was based only on the MGE cells from their reported filtered dataset [28]. In order to integrate human-on-mouse gene expression data from the MGE, we focused only on one-to-one gene orthologs that were present in both datasets. This led us to a combined dataset of 1256 mouse and 358 human MGE cells across 10,362 genes. We used SCTransform from Seurat to integrate both datasets, based on the top 3000 most variable features and the first 30 principal components [47]. For conserved gene markers in the integrated dataset we used the FindConservedMarkers function form Seurat with species as the grouping variable, focusing only in positive markers, a min percentage of 25% cells and a logFC threshold of 0.25. We considered only conserved markers, those with a maximum combined p-value < 0.05.

4.8. Gene Regulatory Network Inference

To infer the gene regulatory network, we ran SCENIC [32] in R separately for the human and mouse MGEs [32]. Briefly, we used the GENIE3 algorithm to infer the gene regulatory network from the human MGE expression data, and regulons for a particular transcription factor were based on DNA motif analysis using RcisTarget. We calculated the regulon specificity scores (RSS) [48] to identify specific regulators within each of the species-conserved cell clusters based on the individual AUCell scores.

Supplementary Materials

The following supporting information can be downloaded at https://www.mdpi.com/article/10.3390/ijms24098122/s1.

Author Contributions

F.K. and M.L. conceived the study. F.K. carried out all experimental work. J.M.-S. analysed the sequencing data with guidance from C.W. A.E.R. supplied the foetal tissue. F.K., J.M.-S. and M.L. interpreted the data and wrote the manuscript with comments from all authors. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by a Wellcome trust Ph.D. studentship to F.K., a grant to M.L. from the UK Medical Research Council (MR/L020807/1) and funding to C.W. from UK Dementia Research Institute which is funded by the UK Medical Research Council, Alzheimer’s Society and Alzheimer’s Research UK (MC_PC_17112). The SWIFT tissue bank is funded by the MRC (MR/M02475X/1) and the Welsh government via Cardiff Brain Unit to A.E.R.

Informed Consent Statement

The informed consent procedures are approved by the national research ethics service.

Data Availability Statement

The data discussed in this publication are deposited in NCBI’s Gene Expression Omnibus and are accessible through GEO Series accession number GSE230467 (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE230467).

Acknowledgments

We would like to thank Sophie Rowlands for collection and dissection of foetal tissues. Thanks also to Nenad Sestan and Zhen Li for providing cell annotations of the mouse MGE dataset. We acknowledge the support of the Supercomputing Wales project, which is part-funded by the European Regional Development Fund (ERDF) via the Welsh Government.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Lewis, D.A.; Hashimoto, T.; Volk, D.W. Cortical inhibitory neurons and schizophrenia. Nat. Rev. Neurosci. 2005, 6, 312–324. [Google Scholar] [CrossRef] [PubMed]
  2. Yizhar, O.; Fenno, L.E.; Prigge, M.; Schneider, F.; Davidson, T.J.; O’Shea, D.J.; Sohal, V.S.; Goshen, I.; Finkelstein, J.; Paz, J.T.; et al. Neocortical excitation/inhibition balance in information processing and social dysfunction. Nature 2011, 477, 171–178. [Google Scholar] [CrossRef] [PubMed]
  3. Marin, O. Interneuron dysfunction in psychiatric disorders. Nat. Rev. Neurosci. 2012, 13, 107–120. [Google Scholar] [CrossRef]
  4. Rubenstein, J.L.; Merzenich, M.M. Model of autism: Increased ratio of excitation/inhibition in key neural systems. Genes Brain Behav. 2003, 2, 255–267. [Google Scholar] [CrossRef]
  5. Hirano, S.; Hirano, Y.; Maekawa, T.; Obayashi, C.; Oribe, N.; Kuroki, T.; Kanba, S.; Onitsuka, T. Abnormal neural oscillatory activity to speech sounds in schizophrenia: A magnetoencephalography study. J. Neurosci. 2008, 28, 4897–4903. [Google Scholar] [CrossRef]
  6. Sohal, V.S.; Zhang, F.; Yizhar, O.; Deisseroth, K. Parvalbumin neurons and gamma rhythms enhance cortical circuit performance. Nature 2009, 459, 698–702. [Google Scholar] [CrossRef]
  7. Powell, E.M.; Campbell, D.B.; Stanwood, G.D.; Davis, C.; Noebels, J.L.; Levitt, P. Genetic disruption of cortical interneuron development causes region- and GABA cell type-specific deficits, epilepsy, and behavioral dysfunction. J. Neurosci. 2003, 23, 622–631. [Google Scholar] [CrossRef] [PubMed]
  8. Paterno, R.; Casalia, M.; Baraban, S.C. Interneuron deficits in neurodevelopmental disorders: Implications for disease pathology and interneuron-based therapies. Eur. J. Paediatr. Neurol. 2020, 24, 81–88. [Google Scholar] [CrossRef] [PubMed]
  9. Zecevic, N.; Hu, F.; Jakovcevski, I. Interneurons in the developing human neocortex. Dev. Neurobiol. 2011, 71, 18–33. [Google Scholar] [CrossRef]
  10. Hansen, D.V.; Lui, J.H.; Flandin, P.; Yoshikawa, K.; Rubenstein, J.L.; Alvarez-Buylla, A.; Kriegstein, A.R. Non-epithelial stem cells and cortical interneuron production in the human ganglionic eminences. Nat. Neurosci. 2013, 16, 1576–1587. [Google Scholar] [CrossRef]
  11. Al-Jaberi, N.; Lindsay, S.; Sarma, S.; Bayatti, N.; Clowry, G.J. The early fetal development of human neocortical GABAergic interneurons. Cereb. Cortex 2015, 25, 631–645. [Google Scholar] [CrossRef] [PubMed]
  12. Whittington, M.A.; Traub, R.D. Interneuron diversity series: Inhibitory interneurons and network oscillations in vitro. Trends Neurosci. 2003, 26, 676–682. [Google Scholar] [CrossRef]
  13. Tremblay, R.; Lee, S.; Rudy, B. GABAergic Interneurons in the Neocortex: From Cellular Properties to Circuits. Neuron 2016, 91, 260–292. [Google Scholar] [CrossRef] [PubMed]
  14. Reynolds, G.P.; Beasley, C.L.; Zhang, Z.J. Understanding the neurotransmitter pathology of schizophrenia: Selective deficits of subtypes of cortical GABAergic neurons. J. Neural. Transm. 2002, 109, 881–889. [Google Scholar] [CrossRef]
  15. Lin, L.C.; Sibille, E. Reduced brain somatostatin in mood disorders: A common pathophysiological substrate and drug target? Front. Pharmacol. 2013, 4, 110. [Google Scholar] [CrossRef] [PubMed]
  16. Ito-Ishida, A.; Ure, K.; Chen, H.; Swann, J.W.; Zoghbi, H.Y. Loss of MeCP2 in Parvalbumin-and Somatostatin-Expressing Neurons in Mice Leads to Distinct Rett Syndrome-like Phenotypes. Neuron 2015, 88, 651–658. [Google Scholar] [CrossRef]
  17. Darmanis, S.; Sloan, S.A.; Zhang, Y.; Enge, M.; Caneda, C.; Shuer, L.M.; Hayden Gephart, M.G.; Barres, B.A.; Quake, S.R. A survey of human brain transcriptome diversity at the single cell level. Proc. Natl. Acad. Sci. USA 2015, 112, 7285–7290. [Google Scholar] [CrossRef]
  18. Lake, B.B.; Ai, R.; Kaeser, G.E.; Salathia, N.S.; Yung, Y.C.; Liu, R.; Wildberg, A.; Gao, D.; Fung, H.L.; Chen, S.; et al. Neuronal subtypes and diversity revealed by single-nucleus RNA sequencing of the human brain. Science 2016, 352, 1586–1590. [Google Scholar] [CrossRef]
  19. Boldog, E.; Bakken, T.E.; Hodge, R.D.; Novotny, M.; Aevermann, B.D.; Baka, J.; Borde, S.; Close, J.L.; Diez-Fuertes, F.; Ding, S.L.; et al. Transcriptomic and morphophysiological evidence for a specialized human cortical GABAergic cell type. Nat. Neurosci. 2018, 21, 1185–1195. [Google Scholar] [CrossRef]
  20. Cambray, S.; Arber, C.; Little, G.; Dougalis, A.G.; de Paola, V.; Ungless, M.A.; Li, M.; Rodriguez, T.A. Activin induces cortical interneuron identity and differentiation in embryonic stem cell-derived telencephalic neural precursors. Nat. Commun. 2012, 3, 841. [Google Scholar] [CrossRef]
  21. Nicholas, C.R.; Chen, J.; Tang, Y.; Southwell, D.G.; Chalmers, N.; Vogt, D.; Arnold, C.M.; Chen, Y.J.; Stanley, E.G.; Elefanty, A.G.; et al. Functional maturation of hPSC-derived forebrain interneurons requires an extended timeline and mimics human neural development. Cell Stem Cell 2013, 12, 573–586. [Google Scholar] [CrossRef] [PubMed]
  22. Kim, T.G.; Yao, R.; Monnell, T.; Cho, J.H.; Vasudevan, A.; Koh, A.; Peeyush, K.T.; Moon, M.; Datta, D.; Bolshakov, V.Y.; et al. Efficient specification of interneurons from human pluripotent stem cells by dorsoventral and rostrocaudal modulation. Stem Cells 2014, 32, 1789–1804. [Google Scholar] [CrossRef] [PubMed]
  23. Maroof, A.M.; Keros, S.; Tyson, J.A.; Ying, S.W.; Ganat, Y.M.; Merkle, F.T.; Liu, B.; Goulburn, A.; Stanley, E.G.; Elefanty, A.G.; et al. Directed differentiation and functional maturation of cortical interneurons from human embryonic stem cells. Cell Stem Cell 2013, 12, 559–572. [Google Scholar] [CrossRef] [PubMed]
  24. Noakes, Z.; Keefe, F.; Tamburini, C.; Kelly, C.M.; Cruz Santos, M.; Dunnett, S.B.; Errington, A.C.; Li, M. Human Pluripotent Stem Cell-Derived Striatal Interneurons: Differentiation and Maturation In Vitro and in the Rat Brain. Stem Cell Rep. 2019, 12, 191–200. [Google Scholar] [CrossRef]
  25. Cruz-Santos, M.; Li, M. Identification of TGFβ signaling as a regulator of interneuron neurogenesis in a human pluripotent stem cell model. Neuronal. Signal. 2021, 5, NS20210020. [Google Scholar] [CrossRef]
  26. Chen, Y.J.; Friedman, B.A.; Ha, C.; Durinck, S.; Liu, J.; Rubenstein, J.L.; Seshagiri, S.; Modrusan, Z. Single-cell RNA sequencing identifies distinct mouse medial ganglionic eminence cell types. Sci. Rep. 2017, 7, 45656. [Google Scholar] [CrossRef]
  27. Mayer, C.; Hafemeister, C.; Bandler, R.C.; Machold, R.; Batista Brito, R.; Jaglin, X.; Allaway, K.; Butler, A.; Fishell, G.; Satija, R. Developmental diversification of cortical inhibitory interneurons. Nature 2018, 555, 457–462. [Google Scholar] [CrossRef]
  28. Mi, D.; Li, Z.; Lim, L.; Li, M.; Moissidis, M.; Yang, Y.; Gao, T.; Hu, T.X.; Pratt, T.; Price, D.J.; et al. Early emergence of cortical interneuron diversity in the mouse embryo. Science 2018, 360, 81–85. [Google Scholar] [CrossRef]
  29. Chen, L.; Chatterjee, M.; Li, J.Y. The mouse homeobox gene Gbx2 is required for the development of cholinergic interneurons in the striatum. J. Neurosci. 2010, 30, 14824–14834. [Google Scholar] [CrossRef]
  30. Fragkouli, A.; van Wijk, N.V.; Lopes, R.; Kessaris, N.; Pachnis, V. LIM homeodomain transcription factor-dependent specification of bipotential MGE progenitors into cholinergic and GABAergic striatal interneurons. Development 2009, 136, 3841–3851. [Google Scholar] [CrossRef]
  31. McKinsey, G.L.; Lindtner, S.; Trzcinski, B.; Visel, A.; Pennacchio, L.A.; Huylebroeck, D.; Higashi, Y.; Rubenstein, J.L. Dlx1&2-dependent expression of Zfhx1b (Sip1, Zeb2) regulates the fate switch between cortical and striatal interneurons. Neuron 2013, 77, 83–98. [Google Scholar] [CrossRef] [PubMed]
  32. Aibar, S.; Gonzalez-Blas, C.B.; Moerman, T.; Huynh-Thu, V.A.; Imrichova, H.; Hulselmans, G.; Rambow, F.; Marine, J.C.; Geurts, P.; Aerts, J.; et al. SCENIC: Single-cell regulatory network inference and clustering. Nat. Methods 2017, 14, 1083–1086. [Google Scholar] [CrossRef] [PubMed]
  33. Hu, J.S.; Vogt, D.; Sandberg, M.; Rubenstein, J.L. Cortical interneuron development: A tale of time and space. Development 2017, 144, 3867–3878. [Google Scholar] [CrossRef] [PubMed]
  34. Wamsley, B.; Fishell, G. Genetic and activity-dependent mechanisms underlying interneuron diversity. Nat. Rev. Neurosci. 2017, 18, 299–309. [Google Scholar] [CrossRef]
  35. Shi, Y.; Wang, M.; Mi, D.; Lu, T.; Wang, B.; Dong, H.; Zhong, S.; Chen, Y.; Sun, L.; Zhou, X.; et al. Mouse and human share conserved transcriptional programs for interneuron development. Science 2021, 374, eabj6641. [Google Scholar] [CrossRef] [PubMed]
  36. Zhao, Z.; Zhang, D.; Yang, F.; Xu, M.; Zhao, S.; Pan, T.; Liu, C.; Liu, Y.; Wu, Q.; Tu, Q.; et al. Evolutionarily conservative and non-conservative regulatory networks during primate interneuron development revealed by single-cell RNA and ATAC sequencing. Cell Res. 2022, 32, 425–436. [Google Scholar] [CrossRef]
  37. Bandler, R.C.; Mayer, C.; Fishell, G. Cortical interneuron specification: The juncture of genes, time and geometry. Curr. Opin. Neurobiol. 2017, 42, 17–24. [Google Scholar] [CrossRef]
  38. Kessaris, N.; Magno, L.; Rubin, A.N.; Oliveira, M.G. Genetic programs controlling cortical interneuron fate. Curr. Opin. Neurobiol. 2014, 26, 79–87. [Google Scholar] [CrossRef]
  39. Colasante, G.; Lignani, G.; Rubio, A.; Medrihan, L.; Yekhlef, L.; Sessa, A.; Massimino, L.; Giannelli, S.G.; Sacchetti, S.; Caiazzo, M.; et al. Rapid Conversion of Fibroblasts into Functional Forebrain GABAergic Interneurons by Direct Genetic Reprogramming. Cell Stem Cell 2015, 17, 719–734. [Google Scholar] [CrossRef]
  40. Laub, F.; Aldabe, R.; Friedrich, V., Jr.; Ohnishi, S.; Yoshida, T.; Ramirez, F. Developmental expression of mouse Kruppel-like transcription factor KLF7 suggests a potential role in neurogenesis. Dev. Biol. 2001, 233, 305–318. [Google Scholar] [CrossRef]
  41. Blackmore, M.G.; Wang, Z.; Lerch, J.K.; Motti, D.; Zhang, Y.P.; Shields, C.B.; Lee, J.K.; Goldberg, J.L.; Lemmon, V.P.; Bixby, J.L. Kruppel-like Factor 7 engineered for transcriptional activation promotes axon regeneration in the adult corticospinal tract. Proc. Natl. Acad. Sci. USA 2012, 109, 7517–7522. [Google Scholar] [CrossRef] [PubMed]
  42. John, A.; Brylka, H.; Wiegreffe, C.; Simon, R.; Liu, P.; Juttner, R.; Crenshaw, E.B., 3rd; Luyten, F.P.; Jenkins, N.A.; Copeland, N.G.; et al. Bcl11a is required for neuronal morphogenesis and sensory circuit formation in dorsal spinal cord development. Development 2012, 139, 1831–1841. [Google Scholar] [CrossRef]
  43. Hashimoto, R.; Nakazawa, T.; Tsurusaki, Y.; Yasuda, Y.; Nagayasu, K.; Matsumura, K.; Kawashima, H.; Yamamori, H.; Fujimoto, M.; Ohi, K.; et al. Whole-exome sequencing and neurite outgrowth analysis in autism spectrum disorder. J. Hum. Genet. 2016, 61, 199–206. [Google Scholar] [CrossRef]
  44. Wang, Y.; Koh, K.; Ichinose, Y.; Yasumura, M.; Ohtsuka, T.; Takiyama, Y. A de novo mutation in the NALCN gene in an adult patient with cerebellar ataxia associated with intellectual disability and arthrogryposis. Clin. Genet. 2016, 90, 556–557. [Google Scholar] [CrossRef] [PubMed]
  45. Evtouchenko, L.; Studer, L.; Spenger, C.; Dreher, E.; Seiler, R.W. A mathematical model for the estimation of human embryonic and fetal age. Cell Transpl. 1996, 5, 453–464. [Google Scholar] [CrossRef] [PubMed]
  46. Bray, N.L.; Pimentel, H.; Melsted, P.; Pachter, L. Near-optimal probabilistic RNA-seq quantification. Nat. Biotechnol. 2016, 34, 525–527. [Google Scholar] [CrossRef]
  47. Hao, Y.; Hao, S.; Andersen-Nissen, E.; Mauck, W.M., 3rd; Zheng, S.; Butler, A.; Lee, M.J.; Wilk, A.J.; Darby, C.; Zager, M.; et al. Integrated analysis of multimodal single-cell data. Cell 2021, 184, 3573–3587.e29. [Google Scholar] [CrossRef]
  48. Suo, S.; Zhu, Q.; Saadatpour, A.; Fei, L.; Guo, G.; Yuan, G.C. Revealing the Critical Regulators of Cell Identity in the Mouse Cell Atlas. Cell Rep. 2018, 25, 1436–1445.e3. [Google Scholar] [CrossRef]
Figure 1. Transcriptional landscape of human MGE. (A), Single-cell RNA sequencing of human MGE, up to two samples were taken from 10 to 15 PCWs as illustrated in the timeline. A total of 358 cells passed the quality control and were analysed. (B), UMAP visualisation of cells coloured by post conceptional week. (C), UMAP visualisation of cells coloured by cell cycle phase. (D), UMAP after the difference between cell cycle scores (S-G2M) was regressed out from the data, colours indicate each of the three unbiased clusters. (E), Expression profile of neural progenitor and neuronal genes visualised by UMAP. (F), Top four biological processes enriched among the 64, 34 and 198 gene markers that characterise cluster 0, 1 and 2, respectively (from left to right). Gene Ontology (GO) terms were ordered based on their lowest adjusted p value. Dot size indicates the number of genes annotated to each GO term while dashed line indicates an adjusted p value threshold of 0.05.
Figure 1. Transcriptional landscape of human MGE. (A), Single-cell RNA sequencing of human MGE, up to two samples were taken from 10 to 15 PCWs as illustrated in the timeline. A total of 358 cells passed the quality control and were analysed. (B), UMAP visualisation of cells coloured by post conceptional week. (C), UMAP visualisation of cells coloured by cell cycle phase. (D), UMAP after the difference between cell cycle scores (S-G2M) was regressed out from the data, colours indicate each of the three unbiased clusters. (E), Expression profile of neural progenitor and neuronal genes visualised by UMAP. (F), Top four biological processes enriched among the 64, 34 and 198 gene markers that characterise cluster 0, 1 and 2, respectively (from left to right). Gene Ontology (GO) terms were ordered based on their lowest adjusted p value. Dot size indicates the number of genes annotated to each GO term while dashed line indicates an adjusted p value threshold of 0.05.
Ijms 24 08122 g001
Figure 2. Expression of regional and interneuron subtype genes. (A), UMAP coloured by the gene expression levels of MGE, LGE and CGE markers. (B), Heatmap shows gene expression patterns of selected genes across PCW time points. (C), UMAP coloured by the gene expression of markers that characterise interneuron LHX6 and LXH8 populations. (D), Gene expression patterns of post-mitotic interneuron markers NPY, SST, CALB2. Bar plot indicates the percentage of neurons that express each marker (relative to the MAP2-expressing cells).
Figure 2. Expression of regional and interneuron subtype genes. (A), UMAP coloured by the gene expression levels of MGE, LGE and CGE markers. (B), Heatmap shows gene expression patterns of selected genes across PCW time points. (C), UMAP coloured by the gene expression of markers that characterise interneuron LHX6 and LXH8 populations. (D), Gene expression patterns of post-mitotic interneuron markers NPY, SST, CALB2. Bar plot indicates the percentage of neurons that express each marker (relative to the MAP2-expressing cells).
Ijms 24 08122 g002
Figure 3. Integration of human and mouse MGE gene expression data and identification of conserved cell clusters. UMAP was calculated based on the first 30 principal components of the top 3000 most variable genes. (A), UMAP of the integrated human and mouse data has been coloured by time points on embryonic days E12.5 and E14.5 (mouse) and 10–15 PCWs (human), different shapes denote each species. (B), UMAP coloured by cell cycle phase with shapes identifying each species. (C), UMAP coloured by the gene expression levels of NES and MAP2 on the SCT normalised dataset, plot has been divided into a single panel per species. (D), UMAP coloured by the integrated cell clusters from the MGE (mouse and human, denoted as i-clusters), again each species has been plotted in a separate panel. (E), Alluvial plot showing the correspondence between the cell clusters identified only in the human dataset (h-clusters) and their correspondence to the integrated clusters based on the combined analysis of human and mouse data.
Figure 3. Integration of human and mouse MGE gene expression data and identification of conserved cell clusters. UMAP was calculated based on the first 30 principal components of the top 3000 most variable genes. (A), UMAP of the integrated human and mouse data has been coloured by time points on embryonic days E12.5 and E14.5 (mouse) and 10–15 PCWs (human), different shapes denote each species. (B), UMAP coloured by cell cycle phase with shapes identifying each species. (C), UMAP coloured by the gene expression levels of NES and MAP2 on the SCT normalised dataset, plot has been divided into a single panel per species. (D), UMAP coloured by the integrated cell clusters from the MGE (mouse and human, denoted as i-clusters), again each species has been plotted in a separate panel. (E), Alluvial plot showing the correspondence between the cell clusters identified only in the human dataset (h-clusters) and their correspondence to the integrated clusters based on the combined analysis of human and mouse data.
Ijms 24 08122 g003
Figure 4. Human active regulons within the integrated MGE clusters. We inferred the gene regulatory network based on co-expression and TF binding site enrichment separately for human MGE using SCENIC. (A), The regulon specific scores (RSS) allow us to identify active regulons in the integrated MGE cell clusters. (B), Inferred regulatory network showing the high confidence links among the conserved markers in the integrated neuronal cluster 0 and their immediate upstream regulators. Genes are coloured in red if they were identified as conserved gene markers in the integrated mouse/human MGE cluster 0. Transcription factors are shaped in diamonds.
Figure 4. Human active regulons within the integrated MGE clusters. We inferred the gene regulatory network based on co-expression and TF binding site enrichment separately for human MGE using SCENIC. (A), The regulon specific scores (RSS) allow us to identify active regulons in the integrated MGE cell clusters. (B), Inferred regulatory network showing the high confidence links among the conserved markers in the integrated neuronal cluster 0 and their immediate upstream regulators. Genes are coloured in red if they were identified as conserved gene markers in the integrated mouse/human MGE cluster 0. Transcription factors are shaped in diamonds.
Ijms 24 08122 g004
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Keefe, F.; Monzón-Sandoval, J.; Rosser, A.E.; Webber, C.; Li, M. Single-Cell Transcriptomics Reveals Conserved Regulatory Networks in Human and Mouse Interneuron Development. Int. J. Mol. Sci. 2023, 24, 8122. https://doi.org/10.3390/ijms24098122

AMA Style

Keefe F, Monzón-Sandoval J, Rosser AE, Webber C, Li M. Single-Cell Transcriptomics Reveals Conserved Regulatory Networks in Human and Mouse Interneuron Development. International Journal of Molecular Sciences. 2023; 24(9):8122. https://doi.org/10.3390/ijms24098122

Chicago/Turabian Style

Keefe, Francesca, Jimena Monzón-Sandoval, Anne E. Rosser, Caleb Webber, and Meng Li. 2023. "Single-Cell Transcriptomics Reveals Conserved Regulatory Networks in Human and Mouse Interneuron Development" International Journal of Molecular Sciences 24, no. 9: 8122. https://doi.org/10.3390/ijms24098122

APA Style

Keefe, F., Monzón-Sandoval, J., Rosser, A. E., Webber, C., & Li, M. (2023). Single-Cell Transcriptomics Reveals Conserved Regulatory Networks in Human and Mouse Interneuron Development. International Journal of Molecular Sciences, 24(9), 8122. https://doi.org/10.3390/ijms24098122

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop