Skip to main content

Revealing cellular and molecular complexity of the central nervous system using single cell sequencing


The mammalian central nervous system (CNS) is one of the most complex systems, with thousands of cell types and subtypes with distinct and unique morphology and gene expression profiles. Based on classic histological methods and conventional cellular and molecular approaches, single cell sequencing is becoming a powerful tool to uncover the complexity of the CNS. In this review, we summarize the principle of single cell sequencing and highlight its use for studying the development of neural stem cells, neural progenitors, and distinct neurons. By revealing transcriptomes in each individual cell using single cell sequencing, we are now able to dissect the cellular heterogeneity of a hundred billion cells in the CNS and comprehensively investigate mechanisms of brain development and function at the cellular and molecular levels.


Advancing our knowledge in neuroscience relies on the fast development of modern technologies, such as next-generation sequencing (NGS), optogenetic modulation, and CRISPR-Cas9 [1,2,3]. These technologies have been used to investigate brain development and function, for example, brain morphology and electrophysiology. Recently, single cell sequencing has explored new aspects of stem cell biology and neuroscience and generated exciting discoveries based on traditional classification of cell types and subtypes in the central nervous system (CNS).

In this review, we summarize the principle of single cell sequencing and highlight its application in neuroscience. We first introduce methods of single cell sequencing, such as single cell isolation, whole-genome amplification (WGA), and whole-transcriptome amplification (WTA). We next reveal the application of single cell sequencing for classifying cell types in the CNS, for understanding molecular mechanisms of development of neural stem cells and neural progenitors in human brains, and for modeling human brain formation and disorders.

The principle of single cell sequencing

The general procedure of single cell sequencing consists of six steps: isolation of single cells; cell lysis to obtain DNA or RNA; addition of barcodes in single cells; amplification of DNA and RNA for sequencing; library preparation and sequencing; and data analysis (Fig. 1). Hierarchical clustering and principle component analysis (PCA) have been used to verify novel cell populations and distinct cell types through identification of new markers in the single cell transcriptomes.

Fig. 1
figure 1

Single cell sequencing flow chart. Brain tissues from the brain region of interest are collected, then single cells are captured by fluorescence-activated cell sorting (FACS), laser capture microdissection (LCM), and microfluidics (Fluidigm C1 and 10× Genomics Chromium). Once single cells are captured, DNA or RNA are amplified for the next sequencing step. Principle component analysis (PCA) and hierarchical clustering are ideal tools to analyze the single cell sequencing data and to classify cell types

Single cell isolation

Single cell isolation is the first key step of single cell sequencing. Several useful methods are available to capture single cells from tissues or cell cultures, such as fluorescence-activated cell sorting (FACS), laser capture microdissection (LCM), manual cell picking, microfluidic, and the use of microwells.

FACS is a classic method to isolate individual cells from brain tissue and cultured cells based on their expression of a reporter gene or their antigens against specific antibodies [4, 5] (Fig. 1). Because FACS requires a large number of cells in suspension as starting material, it is not suitable to capture low-abundance cell subpopulations.

The LCM method can capture cells under direct microscopic visualization, which permits rapid one-step sectioning of complex, heterogeneous tissue to achieve spatial information [6] (Fig. 1).

Microfluidics is a newly developed, highly integrated system that allows sequential processing of small volumes of fluids in channels with dimensions of tens to hundreds of micrometers to achieve single cell culture and sequencing [7]. Several microfluidics platform are available, such as the Fluidigm C1, Drop-seq, and 10× Genomics Chromium. The procedure of the Fluidigm C1 platform includes cell capture, verification, cell lysis, reverse transcription and pre-amplification. Because the capture sizes are tuned to 5–10, 10–17, and17–25 μm in diameter, this system is suitable only for cells that are relatively homogeneous in size. The low capture efficiency also limits its use in capturing and sequencing sticky or non-spherical cells [8]. The Drop-seq and 10× Genomics Chromium systems are droplet-based techniques that enable processing of tens of thousands of cells in a single experiment (Fig. 2). The Drop-seq device allows co-flow of two aqueous solutions across an oil channel: one flow contains microparticle beads that bind barcoded oligonucleotides and the other contains single cells for analysis (Fig. 2) [9]. In the 10× Genomics Chromium system, cells combined with reverse transcription reagents in one channel and gel beads loaded with barcoded oligonucleotides from another channel are mixed to form Gel bead in Emulsion (GEM) (Fig. 2) [10]. After droplet formation, cells are lysed and mRNAs released, which are then hybridized with oligonucleotides based on oligo dT. Droplets are then broken and mRNAs reverse-transcribed in bulk and amplified for sequencing using PCR (Fig. 2).

Fig. 2
figure 2

Flow chart of droplet-based single cell RNA sequencing. In the Drop-seq, one channel contains single cells for analysis and the other contains microparticle beads. The surface of a microparticle bead binds oligonucleotides that consist of oligo dT (green), a unique molecular identifier (UMI; red), a cell barcode (blue), and a PCR primer (brown). Immediately after droplet formation, cells are lysed and mRNAs released and then hybridized with oligonucleotides on the surface of the microparticle beads based on oligo dT binding. Droplets are then broken and mRNAs are reverse-transcribed in bulk and amplified for sequencing using PCR. Moreover, in the 10× Genomics platform, one channel contains single cells for analysis and the other contains gel beads mixed with oligonucleotides that consist of oligo dT, UMI, cell barcode, and a PCR primer. Cells and reagents are next mixed with gel beads. After cell lysis, their mRNAs are released and hybridized with oligonucleotides based on oligo dT binding, and are next reverse-transcribed in bulk and amplified for sequencing using PCR. P1 and P2 are PCR primers for establishing libraries for Illumina sequencing

In microwell sequencing, individual cells are trapped in an agarose microarray and mRNAs subsequently captured on magnetic beads for sequencing [11]. In addition, split-pool ligation-based transcriptome sequencing (SPLiT-seq) eliminates the need to separate individual cells by adding different barcodes to cells over several rounds, so each cell has a unique combination of barcodes for sequencing [12].

Adding barcodes in single cells

Two strategies are most frequently used to add barcodes in single cells in order to distinguish individual cells (Fig. 3). One method is to use Tn5 transposase carrying a specific barcode to add a barcode after amplification of cDNA using oligo dT and unique molecular identifiers (UMI) (Fig. 3a). Another method is to design a primer containing an oligo dT, barcode, and PCR primer which adds a cell-unique barcode when the first cDNA strand is synthesized (Fig. 3b). Once a barcode is added, DNA and cDNA in a single cell are ready for amplification.

Fig. 3
figure 3

Two methods to add barcode in a single cell. a cDNA is reverse-transcribed and amplified using the oligo dT primer (green) in the presence of a template-switching oligonucleotide that contains a unique molecular identifier (UMI; red) and a PCR primer for Illumina sequencing (brown). After amplification, tagmentation is performed using Tn5 transposase carrying a specific cell barcode (blue) and a PCR primer. b cDNA is reverse-transcribed and amplified using an oligonucleotide that consists of oligo dT, UMI, a cell barcode, and a PCR primer, and then tagmented after amplification. P1 and P2 are PCR primers for establishing libraries for Illumina sequencing

Single cell DNA sequencing

To meet the demands of next-generation sequencing, the amount of DNA in a single cell (approximately 6 pg) needs to be amplified using whole-genome amplification (WGA) [13]. Three methods have been applied in WGA: degenerate oligonucleotide-primed PCR (DOP-PCR), multiple displacement amplification (MDA), and multiple annealing and looping-based amplification cycles (MALBAC).

DOP-PCR is widely used in WGA. This method first amplifies the DNA template using a low annealing degenerate primer extension on the DNA template and then amplifies the previous products at a high annealing temperature [14] (Fig. 4a). Because the characteristics of PCR magnify the diversity of different sequences in the genome, DOP-PCR has a low physical coverage of the genome (approximately 10%). This method can accurately retain copy number levels, which makes it an ideal method to detect single cell copy-number variants (CNVs) [15, 16].

Fig. 4
figure 4

Whole-genome amplification methods for single cell sequencing. a Degenerate oligonucleotide-primed PCR (DOP-PCR). The 3′ end of the degenerate oligonucleotide primer (the random six nucleotides) are annealed to the genomic template, allowing the primer to initiate PCR, and PCR fragments are generated to contain the full length of the oligonucleotide primer at one end and the complementary sequence at the other end. Subsequently, the temperature is increased to amplify the DNA fragments. b Multiple displacement amplification (MDA). Double-stranded DNA are melted and random primers are bound to the DNA strand. Branched structures are produced under isothermal condition using the phi29 DNA polymerase. c Multiple annealing and looping-based amplification cycles (MALBAC). Double-stranded DNA is denatured into single strands at 94 °C, then single-stranded DNA templates are homogeneously bound with random primers at 0 °C. Semi-amplified products (n) are further amplified to produce full amplicons during a subsequent five temperature cycles (m), and the complete amplification products are 5′ and 3′ complementary to each other. Cyclization of the complete amplification product is performed with the temperature dropped to 58 °C to prevent further amplification and hybridization of the sequence. Semi-amplification products and genomic DNA continue to circulate to generate the complete amplification product

MDA is a useful method for identification of SNPs, having a high physical coverage of a single cell genome (approximately 90%) due to the use of special DNA polymerase. For instance, Phi29 DNA polymerase is isolated from Bacillus subtilis and Bst DNA polymerase is isolated from Bacillus stearothermophilus [17,18,19]. Phi29 DNA polymerase is considered a better choice for MDA because it shows higher efficiency and affinity and lower error rates than Bst DNA polymerase [20]. MDA produces branched structures under the isothermal conditions through extending random primers (Fig. 4b). Because it may lead to uneven genome coverage, it is not suitable for CNV validation [19].

MALBAC uses PCR amplification of circular DNA fragments with Bst DNA polymerase [21] (Fig. 4c). MALBAC is used only to make copies of the original genomic DNA by protecting the amplification products, in which the 3′ end is complementary to the sequence on the 5′ end in full amplicons [16]. The two ends hybridize to form looped DNA, which can efficiently prevent the full amplicon from being used as a template, therefore enabling close to linear amplification [21]. Thus, MALBAC not only has accuracy for CNV detection, but also features a low false negative rate for SNV detection.

All these methods provide exponential amplification, which creates amplification bias and errors. A new WGA method, linear amplification via transposon insertion (LIANTI), can reduce exponential amplification bias and errors by combining Tn5 transposition and T7 in vitro transcription [22].

Moreover, single nucleus exome sequencing (SNES) is a newly developed method to capture single cells in G1/0 or G2/M phase in a cell cycle. After amplification of single cell genomes using MDA, SNES can reduce technical errors during the single cell sequencing of exomes and genomes [23, 24].

Single cell RNA sequencing

Due to the tiny amount of total RNA (approximately 10 pg) and mRNA (approximately 0.1 pg) in a single cell [25], to perform single cell RNA sequencing, RNA must be amplified by whole transcriptome amplification (WTA) using either the oligo dT-anchor approach (Tang2009 method) [26,27,28,29,30], template-switching approach (SMART-seq/SMART-seq2 method) [9, 31,32,33,34,35,36,37], in vitro transcription (IVT)-based linear amplification (CEL-seq/CEL-seq2 method) [38,39,40,41], or designed primer-based PCR [42, 43].

In the Tang method, mRNAs are reverse-transcribed into cDNAs using a poly(T) primer with an anchor sequence, then poly(A) tails are added to the first strand cDNAs at the 3′ end. The second-strand cDNAs are synthesized using poly(T) primers with another anchor sequence [25, 26] (Fig. 5a). The drawback of this method is that the limited length of the reverse-transcription product causes a strong 3′ mRNA bias [44].

Fig. 5
figure 5

Whole transcriptome amplification methods for single cell sequencing. a Oligo dT-anchor approach. mRNA is reverse-transcribed with an oligo-dT anchored primer, and poly(A) is added to the 3′ end of the first-strand cDNA to synthesize the second strand cDNA. b Template-switching approach. mRNA is reverse-transcribed by Moloney murine leukemia virus (MMLV) reverse transcriptase, and a few non-templated C nucleotides are added to the 3′ end of the first-strand cDNA to synthesize the second strand cDNA

The SMART-seq method can mitigate this bias by adding a new non-templated C nucleotide to the 5′ end of cDNA, which uses Moloney murine leukemia virus (MMLV) reverse transcriptase to amplify the full-length cDNA [33] (Fig. 5b). SMART-seq2 can obtain increased cDNA yields from single cells and also has higher sensitivity, fewer technical biases, and less variability by more refined reverse transcription, template switching, and pre-amplification [34]. However, exponential amplification may lead to gene expression profiles represented by shorter and less enriched G-C-bias amplicons [45].

The CEL-seq method can overcome this problem to achieve linear amplification. By using a primer designed with an anchored poly(T), a unique barcode, the 5′ Illumina sequencing adaptor, and a T7 promoter, the second-strand synthesis is performed and then the cDNA samples are pooled and the IVT reaction performed [38]. RNAs are then fragmented and ligated with an Illumina 3′ adaptor and reverse transcribed to DNA for sequencing. In the CEL-seq2 method, the sensitivity of the CEL-seq is improved by performing another round of reverse transcription with random primers [39]. Additionally, the UMI labeling technique is applied in single cell RNA sequencing, which can increase the efficiency of cDNA synthesis and quantitative accuracy [46].

Application of single cell sequencing in the nervous system

Single cell sequencing has changed our fundamental understanding of tumor heterogeneity, clonal evolution, and metastatic dissemination [23, 47,48,49,50,51,52]. It has also been used to study the immune system, to define immune cell subpopulations, and to examine gene expression variability, splice variants, and gene-regulatory networks [53,54,55]. In addition to the study of tumors and immunology, neuroscience is another field where single cell sequencing technology has a great advantage.

Classification of cell types in the CNS

The classification of cell types in the CNS has challenged the greatest minds in the history of neuroscience. Early neuroscientists used histological knowledge (morphological characteristics and localization) to classify cell types in the CNS. Ramón y Cajal, the pioneer in neuroscience, systematically described cellular structures in the brain using the Golgi stain method, which has guided neuroscience studies for more than a century [56]. Afterwards, cellular biology approaches have been used to label specific cell surface and nuclear antigens with antibodies to distinguish cell types, and then DNA and RNA sequencing has been used to profile differential gene expression from bulk tissues or cultured cells. However, few known cell surface and nuclear protein markers for specific cells and limited distinction of cell morphology have restricted comprehensive understanding of cellular complexity in the CNS [57]. Single cell sequencing has allowed characterization of various cell types in the CNS at the molecular level (Table 1).

Table 1 Summary of cell types in the CNS identified by single cell sequencing

Cerebral cortex

The cerebral cortex is a central brain region that controls complex cognitive behaviors, such as sensorimotor integration, memory, and social behaviors [58]. Recently, single cell RNA sequencing has been used to identify different neuronal types and subtypes and to discover novel cell-specific markers in the mouse cortex [31, 59, 60]. One study analyzed 3,005 cells isolated from the mouse somatosensory cortex and hippocampal CA1 [59]. These 3,005 cells have been classified into nine molecularly distinct classes by clustering cells such as pyramidal neurons, interneurons, and oligodendrocytes based on identification of the most specific markers using sequencing. Repeating bi-clustering on each of the nine major classes identified 47 molecularly distinct subclasses, for example, seven subclasses of pyramidal neurons, 16 subclasses of interneurons, and six subclasses of oligodendrocytes [59]. Another study analyzed 1,679 cells from the primary visual cortex in the adult mouse and identified 49 transcriptomic cell types, such as 23 GABAergic interneurons and 19 glutamatergic neurons [60]. While conventional research has only identified four interneuron subclasses, the single cell sequencing method has classified interneurons into 16 subclasses in the somatosensory cortex and 23 subclasses in the visual cortex, which allows further investigation of the functional complexity of these cells and brain regions [61].

Furthermore, patch-sequencing, which combines the patch-clamp technique with single cell RNA sequencing, allows us to investigate the link between molecularly defined neuronal subtypes and their corresponding electrophysiological and morphological counterparts [31]. Using this approach, 58 cells from layer 1 (L1) of the mouse neocortex have been classified into two cell types: elongated neurogliaform cells (eNGCs) and single bouquet cells (SBCs). A number of the SBCs express genes that are associated with human disease, for example Npas1 and Npas3, have been implicated in autism spectrum disorders. Moreover, a scalable approach to sequence and quantify RNA molecules from isolated neuronal nuclei from a postmortem human brain was developed, generating 3,227 sets of single-neuron data from six distinct regions of the cortex and identifying 16 neuronal subtypes (eight excitatory and eight inhibitory) [62]. In addition, 35,289 single nuclei from the visual cortex, frontal cortex, and the cerebellar hemisphere of a postmortem human brain were analyzed using an improved high-throughput method for single-nucleus droplet-based sequencing (snDrop-seq), providing a more complete brain cell atlas [63].

These reports have demonstrated the advantage of single cell sequencing to identify complex individual cell types in rodent and human brains at a more precise molecular level, which could not be done previously using conventional approaches (Table 1).


The mammalian retina consists of more than 60 distinct cell types, each playing a specific role in processing visual images [64]. Retina neurons were previously defined by their morphology and physiology, and their complexity has been underestimated. Using Drop-seq analysis, about 44,808 mouse retinal cells were analyzed and 39 transcriptionally distinct cell classes identified [9]. Even though amacrine cells account for only about 10% of retina cells, the dataset has revealed that amacrine cells display the most complex subclasses [9]. From another study, retina bipolar cells (BCs) were sorted out from a transgenic mouse line in which GFP is expressed in all BCs. Using Drop-seq, 15 types of BCs were identified and two novel types were discovered according to the molecular diversity of BCs [65, 66]. Thus, single cell RNA sequencing has revealed the unexpected molecular complexity of retina cells that comprise relatively small populations, such as amacrine cells and bipolar cells, and helped further understanding of the function of the retina (Table 1).


The hypothalamus is capable of sensing light, glucose, and blood osmolarity, and is also one of the brain regions with the highest diversity of neurons [67]. Understanding the cell composition and features of cell type-specific transcription in the hypothalamus is essential to explore the function of this important brain region. One study analyzed 3,131 hypothalamus cells and identified 62 subtypes, such as glutamatergic, dopaminergic, and GABAergic neurons, using single cell RNA sequencing [68]. Another study defined 11 non-neuronal and 34 neuronal cell clusters in the adult mouse hypothalamus [69]. Most of the non-neuronal cell types identified in this study are similar to those found in the cerebral cortex [59]. In contrast, the neuronal cell types are predominantly hypothalamus-specific, indicating that different neuronal components are the basis of different functionalities in different brain regions.

Olfactory system

The mammalian olfactory system is capable of detecting and discriminating between a tremendous number of odorants by following the “one neuron–one receptor” rule. Olfactory sensory neurons (OSNs) detect odorants by expressing one olfactory receptor gene (Olfr) in each neuron [70]. A study analyzed 178 cells from the main olfactory epithelium of adult and newborn mice using single cell RNA sequencing [71]. Surprisingly, a subset of cells were found to express multiple Olfrs, and most of these cells are developmentally immature OSNs [71]. Another study also used single cell RNA sequencing to analyze transcriptomes of single epithelial neurons during development by incorporating cell stage markers in individual transcriptomes which reflect the developmental progression from progenitors to mature OSNs. This study showed that immature OSNs express low levels of multiple Olfrs. As development proceeds, this ability declines and mature OSNs typically express high levels of a single Olfr [72]. Moreover, immature OSNs are not restricted to expressing Olfrs from a single chromosomal region [72]. These studies have revealed a new olfactory development model in which the maturation of OSNs is associated with transition of their expression of low levels of multiple Olfrs to high levels of a single Olfr (Table 1).

Exploring human brain development at the single cell level

The human brain is a complex and highly evolving structure. Single cell sequencing has allowed us for the first time to dissect the cell type diversity and pedigree trajectory of the human brain at the molecular level.

Neural stem cells and progenitors in the human brain

Radial glia cells (RGCs) are neural stem cells that generate cortical neurons and serve as scaffolds for migrating neurons in the mammalian cortex [58]. RGCs have been classified into apical RGCs (aRGCs) that reside in the ventricular zone (VZ) and basal RGCs (bRGCs) that are located in the outer subventricular zone (OSVZ) in the human brain cortex. Although aRGCs and bRGCs have different morphologies, the diversity of their molecular features is unknown. It has been proposed that bRGCs drive the dramatic cortical expansion observed in gyrified brains, such as the human brain [73, 74]. Single cells were captured from microdissected specimens from the VZ and SVZ of the human cortex at gestational weeks 16–18 (GW16–18) for single cell RNA sequencing [75]. bRGCs were found to preferentially express genes related to extracellular matrix formation, migration, and stemness of neural stem cells, including TNC, PTPRZ1, FAM107A, HOPX, and LIFR [75]. Gene enrichment analysis identified the selective expression of genes in the LIFR/STAT3 self-renewal pathway in bRGCs, which suggests that expansion of bRGCs in the OSVZ relies on local production of growth factors and activation of the LIFR/STAT3-mediated signaling pathway.

In addition, a method called fixed and recovered intact single cell RNA (FRISCR) has been developed. It enables targeted mRNA purification and single cell transcriptomic profiling of fixed cells without compromising data quality compared with live cells [76]. Using FRISCR, aRGCs have been found to express ANXA1 and CRYAB, while bRGCs express HOPX, which further indicates complex lineage relationships in the human brain [76]. Moreover, a recent study analyzed 2,309 single cells from the human fetal prefrontal cortex (PFCs) between gestational week 8 (GW8) and GW26 and classified PFCs into six major cell types, such as neural progenitor cells (NPCs), excitatory neurons, and interneurons [77]. This study has found that bRGCs actively proliferate from GW9 to GW16 and undergo neuronal differentiation from GW16 to GW26 [77]. Interestingly, this report has also revealed that intermediate progenitor cells (IPCs) comprise two types produced at two key time points; one is mainly produced by aRGCs at GW10 and the other is produced by bRGCs at GW16 [77].

Modeling human brain development

To study the complexity of human brains, systems for modeling human brain development using human embryonic stem cells (hESCs) and induced pluripotent stem cells (iPSCs) in 2D or 3D culture systems have been developed [78,79,80]. However, whether these in vitro models can faithfully recapitulate human brain development at the cellular and molecular levels remains unclear.

A recent study analyzed 4,368 progenitor-enriched and neuron-enriched cells, which differentiate from hESCs, after in vitro culturing at multiple points such as day 0 (D0), D12, D26, and D54 using single cell sequencing [81]. Clustering single cell transcriptomic data have identified 41 distinct populations, for example, 14 types of progenitors and 19 types of neurons [81]. Moreover, comparisons with primary mouse and human brain cell expression datasets and creation of lineage trees using the Bayesian algorithm have revealed the importance of Wnt/β-catenin signaling in autonomously controlling cell fate commitment of hESC-derived neural progenitors [81, 82]. In another study, human cerebral organoids were generated from ESCs and iPSCs. After single cell sequencing, single cell transcriptomic data from organoids were analyzed in parallel with molecular markers of primary human fetal cortex [83]. Interestingly, organoid cortical cells have been found to display remarkably similar gene expression profiles to those in the human fetal cortex [83]. Because samples of primary human fetal brains are extremely rare and difficult to collect, the organoid culture system, in combination with single cell sequencing, has provided a valuable vehicle to understand early patterning and regional formation of the human fetal brain and to reveal critical genes that are responsible for early human brain malformation.


Single cell sequencing is a powerful tool in the fields of stem cell biology and neuroscience. Single cell sequencing allows classification of cell types and subtypes in the CNS at the transcriptome level, which can help to solve the mystery of the heterogeneity and complex function of the nervous system and uncover mechanisms of brain formation at different developmental stages. Single cell sequencing also helps to understand the molecular relationship between disease status and cellular response and reveals specific neuronal cell types in association with different neurological diseases.

Technically, though single cell sequencing provides useful information on gene expression in individual cells, it has relatively low detection efficiencies and requires cells to be removed from their native environment, resulting in the loss of spatial information [84]. Thus, in situ single cell sequencing is a method to overcome these drawbacks [84,85,86].

Moreover, future development of single cell sequencing should aim to be more refined, high-throughput, economical, and easier to handle in standard research and clinical laboratories. Nevertheless, single cell sequencing has opened a new frontier to classify distinct cell types, identify new biomarkers, and reveal novel mechanisms of CNS development, function, and disorders. Single cell sequencing is becoming a valuable and unique tool in addition to classic morphological methods and conventional cellular and molecular approaches.



Apical radial glia cell


Bipolar cell


Basal radial glia cell


Central nervous system


Copy-number variant


Degenerate oligonucleotide-primed PCR


Elongated neurogliaform cell


Fluorescence-activated cell sorting


Fixed and recovered intact single cell RNA


Gel bead in emulsion


Gestational week


Human embryonic stem cell


Intermediate progenitor cell


Induced pluripotent stem cell


In vitro transcription


Laser capture microdissection


Linear amplification via transposon insertion


Multiple annealing and looping base amplification cycles


Multiple displacement amplification


Moloney murine leukemia virus


Next-generation sequencing


Neural progenitor cells


Olfactory receptor gene


Olfactory sensory neurons


Outer subventricular zone


Principle component analysis


Prefrontal cortex cell


Radial glia cell


Single bouquet cell


Single nucleus exome sequencing


Split-pool ligation-based transcriptome sequencing


Unique molecular identifier


Ventricular zone


Whole-genome amplification


Whole-transcription amplification


  1. Stuber GD, Britt JP, Bonci A. Optogenetic modulation of neural circuits that underlie reward seeking. Biol Psychiatry. 2012;71:1061–7.

    Article  PubMed  Google Scholar 

  2. Heidenreich M, Zhang F. Applications of CRISPR-Cas systems in neuroscience. Nat Rev Neurosci. 2015;17:36–44.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  3. Mills JD, Nalpathamkalam T, Jacobs HIL, Janitz C, Merico D, Hu P, Janitz M. RNA-Seq analysis of the parietal cortex in Alzheimer's disease reveals alternatively spliced isoforms related to lipid metabolism. Neurosci Lett. 2013;536:90–5.

    Article  PubMed  CAS  Google Scholar 

  4. Schwarz JM, Smith SH, Bilbo SD. FACS analysis of neuronal-glial interactions in the nucleus accumbens following morphine administration. Psychopharmacology. 2013;230:525–35.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  5. Basu S, Campbell HM, Dittel BN, Ray A. Purification of specific cell population by fluorescence activated cell sorting (FACS). J Vis Exp. 2010. (

  6. Emmertbuck MR, Bonner RF, Smith PD, Chuaqui RF, Zhuang Z, Goldstein SR, Weiss RA, Liotta LA. Laser capture microdissection. Science. 1996;274:998–1001.

    Article  CAS  Google Scholar 

  7. Whitesides GM. The origins and the future of microfluidics. Nature. 2006;442:368–73.

    Article  PubMed  CAS  Google Scholar 

  8. Kolodziejczyk AA, Kim JK, Svensson V, Marioni JC, Teichmann SA. The technology and biology of single-cell RNA sequencing. Mol Cell. 2015;58:610–20.

    Article  PubMed  CAS  Google Scholar 

  9. Macosko EZ, Basu A, Satija R, Nemesh J, Shekhar K, Goldman M, Tirosh I, Bialas AR, Kamitaki N, Martersteck EM. Highly parallel genome-wide expression profiling of individual cells using Nanoliter droplets. Cell. 2015;161:1202–14.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  10. Zheng GX, Terry JM, Belgrader P, Ryvkin P, Bent ZW, Wilson R, Ziraldo SB, Wheeler TD, McDermott GP, Zhu J, et al. Massively parallel digital transcriptional profiling of single cells. Nat Commun. 2017;8:14049.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  11. Han X, Wang R, Zhou Y, Fei L, Sun H, Lai S, Saadatpour A, Zhou Z, Chen H, Ye F, et al. Mapping the mouse cell atlas by microwell-Seq. Cell. 2018;172:1091–1107.e17.

  12. Rosenberg AB, Roco CM, Muscat RA, Kuchina A, Sample P, Yao Z, Graybuck LT, Peeler DJ, Mukherjee S, Chen W, et al. Single-cell profiling of the developing mouse brain and spinal cord with split-pool barcoding. Science. 2018;360:176–82.

  13. Zhang X, Marjani SL, Hu Z, Weissman SM, Pan X, Wu S. Single-cell sequencing for precise cancer research: progress and prospects. Cancer Res. 2016;76:1305–12.

    Article  PubMed  CAS  Google Scholar 

  14. Telenius H, Carter NP, Bebb CE, Nordenskjold M, Ponder BA, Tunnacliffe A. Degenerate oligonucleotide-primed PCR: general amplification of target DNA by a single degenerate primer. Genomics. 1992;13:718–25.

    Article  PubMed  CAS  Google Scholar 

  15. Hou Y, Wu K, Shi X, Li F, Song L, Wu H, Dean M, Li G, Tsang S, Jiang R, et al. Comparison of variations detection between whole-genome amplification methods used in single-cell resequencing. Gigascience. 2015;4:37.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  16. Huang L, Ma F, Chapman AR, Lu S, Xie XS. Single-cell whole-genome amplification and sequencing: methodology and applications. Annu Rev Genomics Hum Genet. 2015;16:79–102.

    Article  PubMed  CAS  Google Scholar 

  17. Aliotta JM, Pelletier JJ, Ware JL, Moran LS, Benner JS, Kong H. Thermostable Bst DNA polymerase I lacks a 3′→ 5′ proofreading exonuclease activity. Genet Anal: Biomol Eng. 1996;12:185–95.

    Article  CAS  Google Scholar 

  18. Dean FB, Nelson J, Giesler TL, Lasken RS. Rapid amplification of plasmid and phage DNA using Phi29 DNA polymerase and multiply-primed rolling circle amplification. Genome Res. 2001;11:1095–9.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  19. Dean FB, Hosono S, Fang L, Wu X, Faruqi AF, Brayward P, Sun Z, Zong Q, Du Y, Du J. Comprehensive human genome amplification using multiple displacement amplification. Proc Natl Acad Sci U S A. 2002;99:5261–6.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  20. Spits C, Le Caignec C, De Rycke M, Van Haute L, Van Steirteghem A, Liebaers I, Sermon K. Optimization and evaluation of single-cell whole-genome multiple displacement amplification. Hum Mutat. 2006;27:496–503.

    Article  PubMed  CAS  Google Scholar 

  21. Zong C, Lu S, Chapman AR, Xie XS. Genome-wide detection of single-nucleotide and copy-number variations of a single human cell. Science. 2012;338:1622–6.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  22. Chen C, Xing D. Single-cell whole-genome analyses by Linear Amplification via Transposon Insertion (LIANTI). Science. 2017;356:189–94.

  23. Wang Y, Waters J, Leung ML, Unruh A, Roh W, Shi X, Chen K, Scheet P, Vattathil S, Liang H. Clonal evolution in breast cancer revealed by single nucleus genome sequencing. Nature. 2014;512:155–60.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  24. Leung ML, Wang Y, Waters J, Navin NE. SNES: single nucleus exome sequencing. Genome Biol. 2015;16:55.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  25. Wang Y, Navin NE. Advances and applications of single-cell sequencing technologies. Mol Cell. 2015;58:598–609.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  26. Tang F, Barbacioru C, Wang Y, Nordman E, Lee C, Xu N, Wang X, Bodeau J, Tuch BB, Siddiqui A, et al. mRNA-Seq whole-transcriptome analysis of a single cell. Nat Methods. 2009;6:377–82.

    Article  PubMed  CAS  Google Scholar 

  27. Sasagawa Y, Nikaido I, Hayashi T, Danno H, Uno KD, Imai T, Ueda HR. Quartz-Seq: a highly reproducible and sensitive single-cell RNA sequencing method, reveals non-genetic gene-expression heterogeneity. Genome Biol. 2013;14:R31.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  28. Nakamura T, Yabuta Y, Okamoto I, Aramaki S, Yokobayashi S, Kurimoto K, Sekiguchi K, Nakagawa M, Yamamoto T, Saitou M. SC3-seq: a method for highly parallel and quantitative measurement of single-cell gene expression. Nucleic Acids Res. 2015;43:e60.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  29. Fan X, Zhang X, Wu X, Guo H, Hu Y, Tang F, Huang Y. Single-cell RNA-seq transcriptome analysis of linear and circular RNAs in mouse preimplantation embryos. Genome Biol. 2015;16:148.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  30. Sheng K, Cao W, Niu Y, Deng Q, Zong C. Effective detection of variation in single-cell transcriptomes using MATQ-seq. Nat Methods. 2017;14:267–70.

    Article  PubMed  CAS  Google Scholar 

  31. Cadwell CR, Palasantza A, Jiang X, Berens P, Deng Q, Yilmaz M, Reimer J, Shen S, Bethge M, Tolias KF, Sandberg R. Electrophysiological, transcriptomic and morphologic profiling of single neurons using Patch-seq. Nat Biotechnol. 2016;34:199–203.

  32. Fuzik J, Zeisel A, Mate Z, Calvigioni D, Yanagawa Y, Szabo G, Linnarsson S, Harkany T. Integration of electrophysiological recordings with single-cell RNA-seq data identifies neuronal subtypes. Nat Biotechnol. 2016;34:175–83.

    Article  PubMed  CAS  Google Scholar 

  33. Ramskold D, Luo S, Wang YC, Li R, Deng Q, Faridani OR, Daniels GA, Khrebtukova I, Loring JF, Laurent LC, et al. Full-length mRNA-Seq from single-cell levels of RNA and individual circulating tumor cells. Nat Biotechnol. 2012;30:777–82.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  34. Picelli S, Bjorklund AK, Faridani OR, Sagasser S, Winberg G, Sandberg R. Smart-seq2 for sensitive full-length transcriptome profiling in single cells. Nat Methods. 2013;10:1096–8.

    Article  PubMed  CAS  Google Scholar 

  35. Picelli S, Faridani OR, Bjorklund AK, Winberg G, Sagasser S, Sandberg R. Full-length RNA-seq from single cells using Smart-seq2. Nat Protoc. 2014;9:171–81.

    Article  PubMed  CAS  Google Scholar 

  36. Islam S, Kjallquist U, Moliner A, Zajac P, Fan J, Lonnerberg P, Linnarsson S. Characterization of the single-cell transcriptional landscape by highly multiplex RNA-seq. Genome Res. 2011;21:1160–7.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  37. Kivioja T, Vaharautio A, Karlsson K, Bonke M, Enge M, Linnarsson S, Taipale J. Counting absolute numbers of molecules using unique molecular identifiers. Nat Methods. 2011;9:72–4.

    Article  PubMed  CAS  Google Scholar 

  38. Hashimshony T, Wagner F, Sher N, Yanai I. CEL-Seq: single-cell RNA-Seq by multiplexed linear amplification. Cell Rep. 2012;2:666–73.

    Article  PubMed  CAS  Google Scholar 

  39. Hashimshony T, Senderovich N, Avital G, Klochendler A, de Leeuw Y, Anavy L, Gennert D, Li S, Livak KJ, Rozenblatt-Rosen O, et al. CEL-Seq2: sensitive highly-multiplexed single-cell RNA-Seq. Genome Biol. 2016;17:77.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  40. Jaitin DA, Kenigsberg E, Keren-Shaul H, Elefant N, Paul F, Zaretsky I, Mildner A, Cohen N, Jung S, Tanay A, Amit I. Massively parallel single-cell RNA-seq for marker-free decomposition of tissues into cell types. Science. 2014;343:776–9.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  41. Klein AM, Mazutis L, Akartuna I, Tallapragada N, Veres A, Li V, Peshkin L, Weitz DA, Kirschner MW. Droplet barcoding for single-cell transcriptomics applied to embryonic stem cells. Cell. 2015;161:1187–201.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  42. Brennecke P, Anders S, Kim JK, Kolodziejczyk AA, Zhang X, Proserpio V, Baying B, Benes V, Teichmann SA, Marioni JC, Heisler MG. Accounting for technical noise in single-cell RNA-seq experiments. Nat Methods. 2013;10:1093–5.

    Article  PubMed  CAS  Google Scholar 

  43. Chapman AR, He Z, Lu S, Yong J, Tan L, Tang F, Xie XS. Single cell transcriptome amplification with MALBAC. PLoS One. 2015;10:e0120889.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  44. Liang J, Cai W, Sun Z. Single-cell sequencing technologies: current and future. J Genet Genomics. 2014;41:513–28.

    Article  PubMed  Google Scholar 

  45. Hedlund E, Deng Q. Single-cell RNA sequencing: technical advancements and biological applications. Mol Asp Med. 2018;59:36–46.

    Article  CAS  Google Scholar 

  46. Islam S, Zeisel A, Joost S, La Manno G, Zajac P, Kasper M, Lonnerberg P, Linnarsson S. Quantitative single-cell RNA-seq with unique molecular identifiers. Nat Methods. 2014;11:163–6.

    Article  PubMed  CAS  Google Scholar 

  47. Navin N, Kendall J, Troge J, Andrews PW, Rodgers L, Mcindoo J, Cook K, Stepansky A, Levy D, Esposito D. Tumour evolution inferred by single-cell sequencing. Nature. 2011;472:90–4.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  48. Kim K, Lee HW, Lee H, Kim SC, Seo YJ, Chung W, Eum HH, Nam D, Kim J, Joo KM. Single-cell mRNA sequencing identifies subclonal heterogeneity in anti-cancer drug responses of lung adenocarcinoma cells. Genome Biol. 2015;16:127.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  49. Tirosh I, Venteicher AS, Hebert C, Escalante L, Patel AP, Yizhak K, Fisher JM, Rodman C, Mount C, Filbin MG. Single-cell RNA-seq supports a developmental hierarchy in human oligodendroglioma. Nature. 2016;539:309–13.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  50. Gerber T, Willscher E, Loefflerwirth H, Hopp L, Schadendorf D, Schartl M, Anderegg U, Camp G, Treutlein B, Binder H. Mapping heterogeneity in patient-derived melanoma cultures by single-cell RNA-seq. Oncotarget. 2016;8:846–62.

    PubMed Central  Google Scholar 

  51. Ni X, Zhuo M, Su Z, Duan J, Gao Y, Wang Z, Zong C, Bai H, Chapman AR, Zhao J. Reproducible copy number variation patterns among single circulating tumor cells of lung cancer patients. Proc Natl Acad Sci U S A. 2013;110:21083–8.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  52. Lohr J, Adalsteinsson VA, Cibulskis K, Choudhury AD, Rosenberg M, Cruzgordillo P, Francis JM, Zhang C, Shalek AK, Satija R. Whole exome sequencing of circulating tumor cells provides a window into metastatic prostate cancer. Nat Biotechnol. 2014;32:479–84.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  53. Mahata B, Zhang X, Kolodziejczyk AA, Proserpio V, Haimvilmovsky L, Taylor AE, Hebenstreit D, Dingler FA, Moignard V, Gottgens B. Single-cell RNA sequencing reveals T helper cells synthesizing steroids De novo to contribute to immune homeostasis. Cell Rep. 2014;7:1130–42.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  54. Han A, Glanville J, Hansmann L, Davis MM. Linking T-cell receptor sequence to functional phenotype at the single-cell level. Nat Biotechnol. 2014;32:684–92.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  55. Shalek AK, Satija R, Adiconis X, Gertner RS, Gaublomme JT, Raychowdhury R, Schwartz S, Yosef N, Malboeuf CM, Lu D. Single-cell transcriptomics reveals bimodality in expression and splicing in immune cells. Nature. 2013;498:236–40.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  56. Compston A. Texture of the nervous system of man and the vertebrates by Santiago Ramon y Cajal. Volume 1. J Neurol Neurosurg Psychiatry. 2001;70:421.

    Article  PubMed Central  Google Scholar 

  57. Ofengeim D, Giagtzoglou N, Huh D, Zou C, Yuan J. Single-cell RNA sequencing: unraveling the brain one cell at a time. Trends Mol Med. 2017;23:563–76.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  58. Sun T, Hevner RF. Growth and folding of the mammalian cerebral cortex: from molecules to malformations. Nat Rev Neurosci. 2014;15:217–32.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  59. Zeisel A, Munoz-Manchado AB, Codeluppi S, Lonnerberg P, La Manno G, Jureus A, Marques S, Munguba H, He L, Betsholtz C, et al. Brain structure. Cell types in the mouse cortex and hippocampus revealed by single-cell RNA-seq. Science. 2015;347:1138–42.

    Article  PubMed  CAS  Google Scholar 

  60. Tasic B, Menon V, Nguyen TN. Adult mouse cortical cell taxonomy revealed by single cell transcriptomics. Nat Neurosci. 2016;19:335–46.

  61. Wamsley B, Fishell G. Genetic and activity-dependent mechanisms underlying interneuron diversity. Nat Rev Neurosci. 2017;18:299.

    Article  PubMed  CAS  Google Scholar 

  62. Lake BB, Ai R, Kaeser GE, Salathia NS, Yung YC, Liu R, Wildberg A, Gao D, Fung HL, Chen S, et al. Neuronal subtypes and diversity revealed by single-nucleus RNA sequencing of the human brain. Science. 2016;352:1586–90.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  63. Lake BB, Chen S, Sos BC, Fan J, Kaeser GE, Yung YC, Duong TE, Gao D, Chun J, Kharchenko PV, Zhang K. Integrative single-cell analysis of transcriptional and epigenetic states in the human adult brain. Nat Biotechnol. 2017;36:70.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  64. Masland RH. The neuronal organization of the retina. Neuron. 2012;76:266–80.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  65. Shekhar K, Lapan SW, Whitney IE, Tran NM, Macosko EZ, Kowalczyk MS, Adiconis X, Levin JZ, Nemesh J, Goldman M. Comprehensive classification of retinal bipolar neurons by single-cell transcriptomics. Cell. 2016;166:1308.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  66. Euler T, Haverkamp S, Schubert T, Baden T. Retinal bipolar cells: elementary building blocks of vision. Nat Rev Neurosci. 2014;15:507–19.

    Article  PubMed  CAS  Google Scholar 

  67. Alvarez-Bolado G, Grinevich V, Puelles L. Editorial: Development of the hypothalamus. Front Neuroanat. 2015;9:83.

  68. Romanov RA, Zeisel A, Bakker J, Girach F, Hellysaz A, Tomer R, Alpar A, Mulder J, Clotman F, Keimpema E, et al. Molecular interrogation of hypothalamic organization reveals distinct dopamine neuronal subtypes. Nat Neurosci. 2017;20:176–88.

  69. Chen R, Wu X, Jiang L, Zhang Y. Single-cell RNA-Seq reveals hypothalamic cell diversity. Cell Rep. 2017;18:3227–41.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  70. Mombaerts P. Odorant receptor gene choice in olfactory sensory neurons: the one receptor-one neuron hypothesis revisited. Curr Opin Neurobiol. 2004;14:31–6.

    Article  PubMed  CAS  Google Scholar 

  71. Tan L, Li Q, Xie XS. Olfactory sensory neurons transiently express multiple olfactory receptors during development. Mol Syst Biol. 2015;11:844.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  72. Hanchate NK, Kondoh K, Lu Z, Kuang D, Ye X, Qiu X, Pachter L, Trapnell C, Buck LB. Single-cell transcriptomics reveals receptor transformations during olfactory neurogenesis. Science. 2015;350:1251–5.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  73. Dehay C, Kennedy H, Kosik KS. The outer subventricular zone and primate-specific cortical Complexification. Neuron. 2015;85:683–94.

    Article  PubMed  CAS  Google Scholar 

  74. Lui JH, Hansen DV, Kriegstein AR. Development and evolution of the human Neocortex. Cell. 2011;146:18–36.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  75. Pollen AA, Nowakowski TJ, Chen J, Retallack H, Sandoval-Espinosa C, Nicholas CR, Shuga J, Liu SJ, Oldham MC, Diaz A, et al. Molecular identity of human outer radial glia during cortical development. Cell. 2015;163:55–67.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  76. Thomsen ER, Mich JK, Yao Z, Hodge RD, Doyle AM, Jang S, Shehata SI, Nelson AM, Shapovalova NV, Levi BP, Ramanathan S. Fixed single-cell transcriptomic characterization of human radial glial diversity. Nat Methods. 2016;13:87–93.

    Article  PubMed  CAS  Google Scholar 

  77. Zhong S, Zhang S, Fan X, Wu Q, Yan L, Dong J, Zhang H, Li L, Sun L, Pan N, et al. A single-cell RNA-seq survey of the developmental landscape of the human prefrontal cortex. Nature. 2018;555:524–8.

    Article  PubMed  CAS  Google Scholar 

  78. Chambers SM, Fasano CA, Papapetrou EP, Tomishima M, Sadelain M, Studer L. Highly efficient neural conversion of human ES and iPS cells by dual inhibition of SMAD signaling. Nat Biotechnol. 2009;27:275–80.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  79. Espunycamacho I, Michelsen KA, Gall D, Linaro D, Hasche A, Bonnefont J, Bali C, Orduz D, Bilheu A, Herpoel A. Pyramidal neurons derived from human pluripotent stem cells integrate efficiently into mouse brain circuits in vivo. Neuron. 2013;77:440–56.

    Article  CAS  Google Scholar 

  80. Shi Y, Kirwan P, Smith JC, Robinson HPC, Livesey FJ. Human cerebral cortex development from pluripotent stem cells to functional excitatory synapses. Nat Neurosci. 2012;15:477–86.

    Article  PubMed  CAS  Google Scholar 

  81. Yao Z, Mich JK, Ku S, Menon V, Krostag A-R, Martinez RA, Furchtgott L, Mulholland H, Bort S, Fuqua MA. A single-cell roadmap of lineage bifurcation in human ESC models of embryonic brain development. Cell Stem Cell. 2017;20:120–34.

    Article  PubMed  CAS  Google Scholar 

  82. Clevers H, Nusse R. Wnt/β-catenin signaling and disease. Cell. 2012;149:1192–205.

    Article  PubMed  CAS  Google Scholar 

  83. Camp JG, Badsha F, Florio M, Kanton S, Gerber T, Wilschbrauninger M, Lewitus E, Sykes AM, Hevers W, Lancaster MA. Human cerebral organoids recapitulate gene expression programs of fetal neocortex development. Proc Natl Acad Sci U S A. 2015;112:15672–7.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  84. Shah S, Lubeck E, Zhou W, Cai L. In situ transcription profiling of single cells reveals spatial organization of cells in the mouse hippocampus. Neuron. 2016;92:342–57.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  85. Ke R, Mignardi M, Pacureanu A, Svedlund J, Botling J, Wahlby C, Nilsson M. In situ sequencing for RNA analysis in preserved tissue and cells. Nat Methods. 2013;10:857–60.

    Article  PubMed  CAS  Google Scholar 

  86. Lee J, Daugharthy ER, Scheiman J, Kalhor R, Yang JL, Ferrante TC, Terry RC, Jeanty SSF, Li C, Amamoto R. Highly multiplexed subcellular RNA sequencing in situ. Science. 2014;343:1360–3.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

Download references


We thank members of the Sun laboratory for their valuable discussions and advice.


This work was supported by the subsidized Project for Postgraduates’ Innovative Fund in Scientific Research of Huaqiao University (Z.Z.), an R01-MH083680 grant from the NIH/NIMH (T. S.), the China Postdoctoral Science Foundation (2017 M622053), and National Natural Science Foundation of China (81471152, 31771141, and 81701132).

Availability of data and materials

Data sharing is not applicable to this article as no datasets were generated or analyzed during the current study.

Author information

Authors and Affiliations



ZZ and TS wrote the manuscript, NM revised the paper, TS conceived of the study. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Tao Sun.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Zeng, Z., Miao, N. & Sun, T. Revealing cellular and molecular complexity of the central nervous system using single cell sequencing. Stem Cell Res Ther 9, 234 (2018).

Download citation

  • Published:

  • DOI: