Skip to main content

Single-cell sequencing in stem cell biology


Cell-to-cell variation and heterogeneity are fundamental and intrinsic characteristics of stem cell populations, but these differences are masked when bulk cells are used for omic analysis. Single-cell sequencing technologies serve as powerful tools to dissect cellular heterogeneity comprehensively and to identify distinct phenotypic cell types, even within a ‘homogeneous’ stem cell population. These technologies, including single-cell genome, epigenome, and transcriptome sequencing technologies, have been developing rapidly in recent years. The application of these methods to different types of stem cells, including pluripotent stem cells and tissue-specific stem cells, has led to exciting new findings in the stem cell field. In this review, we discuss the recent progress as well as future perspectives in the methodologies and applications of single-cell omic sequencing technologies.


An individual cell is the smallest functional and universal unit of organisms. Gene expression is regulated within or between individual cells, and so, ideally, analyses of gene expression would be performed using single cells; but owing to technical limitations, such as the tiny size of an individual cell, nearly all of the gene-expression studies described in the literature (especially those at a whole-genome scale) have been performed using bulk samples of thousands or even millions of cells. The data based on these ensemble analyses are valid; but the gene expression heterogeneity between individual cells, especially at the whole-genome scale, is still largely unexplored.

Cellular heterogeneity is a general feature of biological tissues that is influenced by both physiological and pathological conditions. Even a ‘pure’ cell type will have heterogeneous gene expression because individual cells may occur in a range of extrinsic microenvironments and niches that influence gene expression, because gene expression may differ throughout the cell cycle, and because of the intrinsic stochastic nature of gene-expression systems [14]. By definition, a stem cell is characterized as both being capable of unlimited self-renewal and having the potential to differentiate into specialized types of cells. Stem cells are generally classified into pluripotent stem cells, which can give rise to cells of all three germ layers (the ectoderm, mesoderm and endoderm), and tissue-specific stem cells, which play essential roles in the development of embryonic tissues and the homeostasis of adult tissues. Pluripotent stem cells in a mammalian early embryo are few in number; tissue-specific stem cells always form a minor proportion of the cell population of a particular tissue or organ. These minor cell populations are thus intermingled with a variety of differentiated and intermediate cell types in the embryonic or adult tissues, forming heterogeneous populations. Single-cell sequencing provides powerful tools for characterizing the omic-scale features of heterogeneous cell populations, including those of stem cells. The beauty of single-cell sequencing technologies is that they permit the dissection of cellular heterogeneity in a comprehensive and unbiased manner, with no need of any prior knowledge of the cell population.

In this review, we discuss the methodologies of recently developed single-cell omic sequencing methods, which include single-cell transcriptome, epigenome, and genome sequencing technologies, and focus on their applications in stem cells, both pluripotent and tissue-specific stem cells. Finally, we briefly discuss the future of methodologies and applications for single-cell sequencing technologies in the stem cell field.

Single-cell RNA-sequencing (RNA-seq) technologies

Introduction of single-cell RNA-seq technologies

RNA-seq technology provides an unbiased view of the transcriptome at single-base resolution. It has been shown that the transcriptome of a mammalian cell can accurately reflect its pluripotent or differentiated status, and it will be of great interest to explore the transcriptome diversity and dynamics of self-renewing and differentiating stem cells at single-cell resolution. The first method for single-cell RNA-seq was reported in 2009, only 2 years after standard RNA-seq technology using millions of cells was developed [5]. Subsequently, many other single-cell RNA-seq methods based on different cell capture, RNA capture, cDNA amplification, and library establishment strategies were reported, including Smart-seq/Smart-seq2 [6, 7], CEL-seq [8], STRT-seq [9, 10], Quartz-seq [11], multiple annealing and looping-based amplification cycles (MALBAC)-RNA [12], Phi29-mRNA amplification (PMA), Semirandom primed polymerase chain reaction (PCR)-based mRNA amplification (SMA) [13], transcriptome in vivo analysis (TIVA) [14], fixed and recovered intact single-cell RNA (FRISCR) [15], Patch-seq [16, 17], microfluidic single-cell RNA-seq [18, 19], massively parallel single-cell RNA-sequencing (MARS-seq) [20], CytoSeq [21], Drop-seq [22] and inDrop [23].

Methods allowing in situ single-cell RNA sequencing or highly multiplexed profiling have also been developed recently [24, 25]. Furthermore, methods for three-dimensional reconstructed RNA-seq at single-cell resolution have also been developed [2628]. A summary of these methods can be found in Table 1, and detailed descriptions of them can also be seen in other recent reviews [2931]. All of these methods detect only poly(A)-plus RNAs from an individual cell and thus miss the important poly(A)-minus RNAs. Recently, we developed the SUPeR-seq technique, which detects both poly(A)-plus and poly(A)-minus RNAs from an individual cell, and we used it to discover several thousands of circular RNAs with no poly(A) tail as well as hundreds of poly(A)-minus linear RNAs in mouse pre-implantation embryos [32].

Table 1 Summary of single-cell RNA-seq technologies

To obtain a comprehensive view of the heterogeneity of a complex population of cells, a large number of individual cells must be sequenced. During the past several years, the throughput of the single-cell RNA-seq technologies has been greatly improved. The microfluidic and robotic systems provide high-throughput strategies that can handle hundreds of individual cells [1821]. Notably, two recently reported methods, Drop-seq and inDrop, dramatically improve the throughput to thousands or even tens of thousands of individual cells for each experimental run by using a combination of the one-bead–one-cell droplet and an unique barcoding strategy [22, 23]. Very different cell types can be distinguished by sequencing as few as 50,000 reads for each cell [33, 34], though deeper sequencing may be necessary to discriminate between types of cells that have relatively subtle differences, such as mouse embryonic stem cells and epiblast stem cells.

Many bioinformatics tools that were designed for bulk RNA-seq analyses are also applicable to single-cell RNA-seq data; further tools have been designed specifically for analyses of single-cell RNA-seq data. An in-depth review of these approaches can be seen elsewhere [35]. Bioinformatics tools have been used in the stem cell field to identify different cell types and sub-populations, as well as their marker genes, from the relatively noisy dataset. Determining sub-populations of stem cells within a dataset is achieved by methods for unbiased clustering and differential gene expression analysis. Zeisel et al. [36] recently described a biclustering-based algorithm called BackSPIN that increases the accuracy of identifying cell types from single-cell RNA-seq data. Grun et al. [37] developed another algorithm called RaceID, which is based on a feature of the single-cell RNA-seq technique that creates extremely low false-positive errors if cross contamination is carefully controlled, especially when unique molecular identifiers (UMIs) are applied. It does, however, generate a high number of false negative errors, where a gene is expressed in a cell, but missed by this technique. These and other methods have greatly improved the analyses of single-cell RNA-seq data in stem cells or embryos. In addition, bioinformatic analysis algorithms such as Monocle and Waterfall have been developed to provide a time-serial reconstruction of a developmental or differentiation process, also using single cell RNA-seq datasets [38, 39]. These algorithms produce a ‘pseudotime’ trajectory through a reduced dimension data space by calculating a minimum spanning tree.

Quantitative assessment of the current single-cell RNA-seq methods shows that these methods have a capture efficiency ranging between 5 % and 60 % [10, 18, 19, 40, 41]. Owing to the biases of molecular capture and amplification, current methods for the sequencing of single cells still have relatively high technical noise, which is acceptable when studying highly expressed genes but which masks the biological variations of genes that are expressed at low levels. Several studies have made great efforts to improve signal-to-noise performance by optimizing the efficiency of reverse transcription and PCR amplification [7], by performing the reactions in nanoliter volumes in a microfluidic system instead of in microliter volumes in tubes [18, 19], through the use of UMIs [10, 33], or by using spike-in of reference mRNAs to discriminate the technical noise and real biological variation signals [42]; nevertheless, there is still much room for improvement.

In the past several years, single-cell RNA-seq methods have been applied to a wide variety of systems, including early mammalian embryos [4348], developing tissues [33, 4951], adult tissues [22, 36, 37, 52, 53], immune cells [20, 21, 5456], cancer cells [6, 5759], and stem cells that are either isolated in vivo [39, 6063] or cultured in vitro [23, 38, 6467]. A flowchart of a typical single-cell RNA-seq project is shown in Fig. 1. The work of Zeisel et al. is an excellent and representative example of these studies, showing that single-cell RNA-seq can identify numerous sub-populations of cells that would be missed if bulk RNA-seq were performed instead [36]. These authors unbiasedly sequenced the transcriptomes of 3005 single cells isolated from the mouse primary somatosensory cortex (S1) and the hippocampal CA1 region. A total of 47 molecularly distinct subclasses of cells were identified, comprising nine major cell types including S1 and CA1 pyramidal neurons, interneurons, oligodendrocytes, astrocytes, microglia, vascular endothelial cells, mural cells, and ependymal cells. This and other studies demonstrate that the current single-cell RNA-seq technology, even with much room for improvement, has become an established and powerful tool that has practical applications in a wide variety of biological fields.

Fig. 1
figure 1

Flowchart of a typical single cell omic sequencing project. A typical single-cell sequencing project comprises four major steps: experimental design, performing the experiments, bioinformatic analyses and verification of the discoveries. Here we use a single-cell RNA-seq project as an example. Note that if the project fails at any step, researchers should go back to previous steps to identify the cause of the failure and re-design accordingly. In a real project, this process may need to be repeated several times. FACS fluorescence-activated cell sorting, MACS magnetic-activated cell sorting, qPCR quantitative polymerase chain reaction

Pluripotent stem cells

Pre-implantation development

Mammalian pre-implantation development represents the start of a new life and involves global gene expression changes during this process. Because the cell numbers during this developmental process are very limited, single-cell RNA-seq provides an unprecedented opportunity to decipher gene expression dynamics during this process. Comprehensive sets of transcriptome profiles from both human and mouse cells undergoing pre-implantation development have been generated [4345]. The gene expression features of the maternal-zygotic transition have been accurately captured. Although cells of the same stage are relatively similar, there is evidence that inter-blastomere differences occur as early as the four-cell stage of mouse embryos [46, 68]. These differences may be functionally relevant to the first cell-fate decision event of the pre-implantation embryo, which is the segregation between the trophectoderm (TE) and the inner cell mass (ICM). Later, the ICM further segregates into primitive endoderm (PE) and pluripotent epiblast (EPI) that give rise to all the cell lineages of the embryo proper. Single-cell RNA-seq analysis offers a comprehensive view of the transcriptome of these divergent cell lineages. It has been shown that Id2 and Sox2 are two early markers that are strongly activated in TE and ICM cells, respectively, during the 16- to 32-cell stage of the mouse embryo [69]. In addition, before the segregation of PE and EPI, a precursor cell expresses both the PE and EPI markers, leading to a model of stochastic cell-to-cell expression heterogeneity that is followed by signal reinforcement and commitment of cell-fate determination [70].

Conservation of gene expression and differences between human and mouse pre-implantation development have been identified on the basis of single-cell RNA-seq data. One documented difference between human and mouse development is the timing of zygotic genome activation, which has been clearly verified using single-cell transcriptome data. In a mixed background mouse (CAST/EiJ × C57BL/6 J), embryos showed rapid maternal transcript clearance and zygotic genome activation at the two-cell stage, as well as significant allele-specific gene expression [45]. In humans, principal component analysis (PCA) and differential gene expression analysis confirmed that zygotic genome activation occurs between the four- and eight-cell stages [44]. A careful comparative analysis revealed many other important differences between human and mouse development [47]. For example, the transcription factor KLF17 is exclusively expressed in the human EPI, and key components of the transforming growth factor (TGF)-β signaling pathway are highly enriched in human, but not mouse, embryos. In addition, the key factors Id2, Elf5, and Eomes are exclusively expressed in TE cells in the mouse, but not in humans.

Embryonic stem cells

Both mouse and human embryonic stem cells (ESCs) serve as excellent in vitro models for studying the self-renewal ability and differentiation potential of pluripotent stem cells. The ICM of blastocysts can form ESCs when cultured in proper pluripotency maintenance conditions, and the derivation of both human and mouse ESCs have been traced using single-cell RNA-seq methods [44, 64]. These studies demonstrate that the outgrowth process is associated with prominent expression changes for transcriptional regulators and for genes that are associated with pluripotency. A comparison between human embryonic stem cells (hESCs) and EPI showed that genes that are involved in pluripotency are conserved, but enriched for different pathways [44, 47]. Human EPI is enriched for oxidative phosphorylation signaling, reflecting a difference in growth environment from that of hESCs, which are cultured in vitro under oxygen-rich conditions and preferentially switch to glycolytic metabolism. The hESCs are enriched for the regulation of cell proliferation and genes involved in the fibroblast growth factor (FGF), MAPK and Wnt signaling pathways, suggesting that the EPI and hESCs have distinct mechanisms for maintaining the pluripotency state.

Although ESCs are relatively homogeneous, they still contain different sub-populations. Single-cell RNA-seq analysis has revealed that many genes have variable expression among individual mouse embryonic stem cells (mESCs) [18, 19] and, importantly, has identified sub-populations that have distinct transcriptomes [23, 65, 66]. By sequencing nearly 1000 individual mESCs using the droplet-barcoding approach, Klein et al. [23] characterized several minor sub-populations, including an epiblast-like sub-population, a Prdm1-high sub-population and an Hsp90-high sub-population. The same study also sequenced thousands of cells to examine the differentiation of mESCs after withdrawal of leukemia inhibitory factor (LIF), and characterized the dynamic changes during differentiation in several sub-populations that do not map to any known cell type.

Primordial germ cells

Primordial germ cells (PGCs) are precursors of mature germ cells—the oocyte and sperm. Single-cell RNA-seq datasets of human PGCs from the migrating stage to the gonadal stage have been created and reveal the dynamic and balanced expression of both pluripotency genes and germline-specific genes during PGC development [62]. Cell populations of early PGCs in mitosis are relatively homogenous, whereas the later female PGCs are highly heterogeneous during meiotic arrest, even within the same embryo. This finding suggests that entry into meiotic arrest is unsynchronized for human female PGCs in vivo. Unique features that distinguish human PGCs from those of mice were also systematically explored. For example, human early PGCs highly expressed SOX15 and SOX17, whereas those from mice express Sox2.

Tissue-specific stem cells

Tissue-specific stem cells reside in developing or differentiated tissues. They also undergo self-renewal and have the potential to differentiate into a variety of specified cell types. In the past 2 years, single-cell RNA-seq methods have been applied to tissue-specific stem cells. These studies have identified novel stem cell types and have dissected cell heterogeneity within a ‘homogenous’ stem cell population.

Identification of novel stem cell types

The study by Treutlein et al. [49] on developing mouse lung epithelium provides an elegant example of how a novel stem cell type could be identified using the single-cell RNA-seq approach. The alveolar type 1 (AT1) and AT2 cells in the lung are two epithelial cell types that play crucial roles in air exchange, but the identity of the alveolar progenitors remains elusive. Treutlein et al. [49] identified five distinct cell populations through evaluation of 80 individual epithelial cells from distal lung regions of E18.5 mouse embryos, which included four known cell types: two bronchiolar lineages (Clara and ciliated cells), and the alveolar type AT1 and AT2 cells. An undefined and interesting fifth cell group co-expresses the marker genes of AT1 and AT2 and is positioned between the populations of AT1 and AT2 cells on the PCA plot, pointing to a bi-potential progenitor population for AT1 and AT2 cells. These alveolar progenitors have been verified by an independent set of experiments, including immunostaining, lineage tracing and clonal analysis [71]. Furthermore, single-cell RNA-seq data allowed Treutlein et al. [49] to characterize the developmental intermediates from the bi-potential progenitor to AT1 and AT2 cells completely, and even to reconstruct a continual route of the differentiation process to reveal the transcriptome dynamics.

Dissecting cell heterogeneity among a stem cell population

Single-cell RNA-seq has been used to dissect cellular heterogeneity within a tissue-specific stem cell population. These studies have revealed both similarities and differences in the structures of the stem cell populations of different tissue types. Hematopoietic stem cells (HSCs) generate all blood lineages. Long-term reconstituting HSCs (LT-HSCs) are at the top of the hematopoietic hierarchy and can undergo self-renewal and division to replenish committed cells, which are called short-term reconstituting HSCs (ST-HSCs). Kowalczyk et al. [60] and Tsang et al. [61] have shown that cell-cycle differences dominate the cell heterogeneity of each HSC type. The cell-cycle progression of HSCs can be re-established using single-cell transcriptome data, which provide a promising new approach for studying the characteristics of quiescent and proliferative stem cells. Analysis of non-cycling cells revealed a clear difference between LT-HSCs and ST-HSCs. Within the LT-HSCs, subgroups of cells that are associated with markers of specific lineages also exist, as revealed by analyzing only the hematopoietic genes, even though these genes may still be related to the cell cycle.

The neural stem cells (NSCs) in the subventricular zone and the subgranular zone of the dentate gyrus continually give rise to new neurons and glia in the adult mammalian brain. The neurogenesis process starts from quiescent NSCs (qNSCs), which become activated NSCs (aNSCs) and, subsequently, early intermediate progenitor cells (eIPCs). Using the ‘pseudotime’ concept to analyze the single-cell transcriptome data, Shin et al. [39] and Llorens-Bobadilla et al. [63] charted a continual developing trajectory for this early neurogenesis process. Like the work of Treutlein et al. [49], these two studies again demonstrate that the single-cell RNA-seq approach can provide a snapshot of the transcriptome dynamics of a developmental process if reasonable numbers of individual cells of the population are sequenced at a given time point.

The continual self-renewal of the intestinal epithelium is another well-established model for studying adult stem cells. Lgr5-positive cells positioned at crypt bottoms serve as the stem cells that fuel the self-renewal process. Grun et al. [37] sequenced nearly 200 green fluorescent protein (GFP)-marked Lgr5-positive cells, and found that these cells formed a single large homogenous population with a few outliers, which indicates a distinct population structure different from that of HSCs and NSCs.

Perturbation of the stem cell populations under non-physiological conditions has also been studied. Llorens-Bobadilla et al. [63] analyzed NSCs in ischemic brain injury. In NSCs under physiological conditions, these authors identified a transition from dormant NSCs to primed-quiescent NSCs and then activated NSCs. In injured NSCs, the authors found that the proportion of dormant NSCs prominently decreases, whereas the primed-quiescent and activated NSCs greatly increase. In another study, Kowalczyk et al. [60] compared young and old mice and found that ageing is associated with a decrease in the length of the G1 phase of the LT-HSCs, which should be linked to LT-HSC accumulation in older mice. In addition, they found that the transcriptome states of the ageing HSCs are inversely correlated with their differentiation states, such that the old ST-HSCs are similar to the young LT-HSCs [60]. Tsang et al. [61] investigated the knockout phenotype of the transcription factor Bcl11a and found abnormal proliferation and selective elimination of lymphoid-competent HSCs in Bcl11a-knockout HSCs [61]. Together, these studies demonstrate that single-cell RNA-seq can provide rich information on the structure of a stem cell population and its behavior under different conditions, and offer great insight into the function of tissue-specific stem cells.

Cancer stem cells

Cancer tissue usually contains sub-populations of cells that have strong phenotypic and functional heterogeneity. The cancer stem cell (CSCs) concept holds that there is a sub-population of highly malignant stem cells at the top of the tumor cell hierarchy. The existence of these CSCs, however, is still controversial in many cancer types. Single-cell RNA-seq has the potential to help identify these cells and, more generally, to provide new insight into complex intra-tumoral heterogeneity. Patel et al. [57] sequenced 672 single cells from five glioblastoma samples. Each tumor showed high intra-tumoral cell heterogeneity in many aspects, including copy number variations as well as cell cycle, hypoxia, and immune response. By examining a set of ‘stemness’ genes, Patel et al. identified continuous, rather than discrete, stemness-related expression states among the individual cells of all five tumors, reflecting the complex stem cell states within a primary tumor. Even though there have been only a few studies addressing the question of tumor transcriptome heterogeneity down to the single-cell resolution [57, 72, 73], a more complete and accurate view of heterogeneity in various cancer types, including the characteristics of the CSC, is expected to be obtained in the near future.

Single-cell epigenome sequencing technologies

The development, maintenance and differentiation of a stem cell are orchestrated by epigenetic modifications of its genome, including covalent modifications of genomic DNA and histones. Cell-to-cell epigenetic variation is an important layer of cell heterogeneity necessary for the transcriptional regulation of gene expression. Of particular interest will be the epigenome heterogeneity that underlies the transcriptome heterogeneity of cell populations such as the pluripotent and adult stem cells described above. In addition, how these heterogeneities are associated with changes in chromosome conformation in individual cells is not yet known. Conventional genome-wide epigenetic methods require millions of cells and cannot identify epigenetic heterogeneity among different individual cells, but recent studies have made great efforts in developing technologies to perform single-cell epigenome analysis (Table 2).

Table 2 Summary of single-cell epigenome sequencing technologies

DNA modifications

DNA methylation is the major DNA modification in the mammalian genome and plays important roles in many developmental processes. Recently, single-cell DNA methylome sequencing methods have been reported by our group and others [7476]. Our method (scRRBS) is based on the reduced representation bisulfite sequencing (RRBS) strategy [77], whereas the methods of Smallwood et al. (scBS-seq) [75] and Farlik et al. (scWGBS) [76] are based on a post-bisulfite adaptor tagging (PBAT) approach [78]. Using these methods, we have charted the DNA methylation landscapes of human and mouse pre-implantation development, as well as human PGC development [62, 74, 79]. These and other studies have comprehensively characterized the two global DNA demethylation waves that occur during mammalian pre-implantation and PGC development at the genome-scale and the single-base resolution [8082]. These studies have shown that human PGCs at about 10 to 11 weeks after gestation have lower methylation levels (6–8 %) than other types of cells including blastocysts (~40 %). This serial hypomethylated DNA methylome dataset of human PGCs in vivo can be used as a standard reference for assessing the quality of PGC-like cells differentiated from hESCs or human induced pluripotent stem cells (hiPSCs) in vitro. Smallwood et al. [75] demonstrated that integration of just 12 single oocyte scBS-seq datasets can largely recover the major pattern of their entire DNA methylome. Although successful, the present single-cell DNA methylome sequencing methods have much sparser coverage than bulk methods, and thus have much room for improvement. In addition to DNA methylation, recent studies have uncovered hydroxymethylation (5hmC) as well as 5-formylcytosine (5fC) and 5-carboxylcytosine (5caC) modifications on genomic DNAs [83]. Although whole genome scale methods for detecting these DNA modifications on bulk cells have been established, methods at the single-cell level still await development in the near future.

Chromatin accessibility and structure

Genomic methods for assessing the chromatin accessibility of bulk cell populations have been effective for identifying active regulatory elements. Several recent studies have adapted these methods to single-cell resolution. The methods of Buenrostro et al. [84] and Cusanovich et al. [85] (scATAC-seq) are based on ATAC-seq (assay for transposase-accessible chromatin) and rely on the ability of the prokaryotic Tn5-transposase to insert preferentially into accessible chromatin regions in the genome. The method used by Jin et al. [86] is based on the more conventional DNase sequencing approach (scDNase-seq). scDNase-seq appears to detect more open chromatin regions per individual cell than scATAC-seq. In addition, chromosome structure capture technologies have recently been adapted to single cell analysis [87, 88]. These methods, which have been shown to distinguish correctly between ESCs and other cell types at different chromatin state layers [86], should be applied to dissect the heterogeneity of chromatin states of stem cell populations in the near future.

Histone modifications

Histone modifications play essential roles in the regulation of gene expression in stem cells. Chromatin immunoprecipitation followed by sequencing (ChIP-seq) is a widely used method for mapping histone modifications at the whole-genome scale. Rotem et al. [89] recently adapted ChIP-seq to a single-cell analysis by combining droplet and barcoding strategies (Drop-ChIP). A stringent negative control using a non-specific IgG antibody was not performed side-by-side for murine embryonic fibroblasts (MEFs) or ESCs, however, leaving the potential non-specific noise in their single-cell ChIP-seq dataset unresolved. Drop-ChIP is able to detect only approximately 1000 H3K4me3 peaks per cell, corresponding to a peak detection sensitivity of approximately 5 %. Nevertheless, the method is capable of separating mouse ESCs into three sub-populations that have distinct H3K4me2 signals over loci bound by pluripotency-associated transcription factors such as Oct4, Sox2, and Nanog, and differentiation-associated transcription factors such as FoxA2, as well as epigenetic repressors including Polycomb and CoREST. The first group of cells has the highest signal for these pluripotency signature genes, the second group has intermediate signals, and the third group has the lowest signals, while H3K4me2 signals for differentiation and epigenetic repressor signature genes are reversed. Thus, these sub-populations may have distinct chromatin states that are related to pluripotency and differentiation priming. This finding implicates a new layer of cell heterogeneity in the epigenome of ESCs. Further improvement of single-cell epigenome sequencing technologies will provide a deeper understanding of the cell heterogeneity of chromatin states in ESCs and other types of stem cells.

Single-cell genome sequencing technologies

The genomes of individual cells carry another layer of information that is useful in revealing the development and heterogeneity of a stem cell population: the cell lineage. During development, one stem cell gives rise to many specialized cells through continuous cell division and differentiation. During each cell division, replication errors may occur. Although such errors (replication mutations) occur at an extremely low frequency in normal mammalian cells (0–1 mutations per cell division), any replication mutations that are detected in individual progeny cells can be used to trace the developmental lineage of those cells. A cell lineage tree, such as the detailed lineage tree that has been illustrated for Caenorhabditis elegans, can greatly help to illustrate a developmental process.

To detect replication mutations in individual cells, a single-cell whole-genome amplification is necessary in order to get enough material for sequencing analysis. This can be accomplished using methods that include degenerate oligonucleotide-primed polymerase chain reaction (DOP-PCR) [90], multiple displacement amplification (MDA) [91], MALBAC [92], microfluidics-based MDA [9395] and MDA for G2/M nuclei (Nuc-seq) [96, 97] (Table 3). Detailed and elegant reviews of these methods can also be seen elsewhere [98, 99].

Table 3 Summary of single-cell genome sequencing technologies

Single-cell genome sequencing has been applied to human germ cells for sperm and oocytes to study meiotic recombination, aneuploidy, and the mutation rate of these cells [93, 100, 101]. These studies have generated the first personal recombination maps of individual men and women and have detected aneuploidy during human gametogenesis [93, 100, 101]. Behjati et al. [102] have also applied genome sequencing for lineage tracing of the development of normal cells. In this study, the early cell lineage and the contribution of these early cells to adult tissues were elucidated by whole-genome sequencing of 25 single-cell-derived organoid lines from the mouse gut and prostate. Single-cell whole-genome sequencing has also been used to study tumor cells. Clonal evolution of a tumor can be elucidated on the basis of the copy number variation (CNV) and single-nucleotide variation (SNV) of single tumor cells [96]. How to authenticate a SNV accurately within a single cell with essentially no false-positive calls remains a challenge. Future improvement of single-cell whole-genome-sequencing technologies will help resolve this issue and will promote the application of this technology for the lineage tracing of stem cells by comprehensively identifying genomic variations within each single stem or differentiated cell in normal or cancerous tissue.


Despite the fact that single-cell sequencing methods have been widely applied to dissecting the heterogeneity of stem cells, all of the currently available single-cell omic sequencing technologies are clearly not ideal. There exist significant technical noise and amplification errors, and they provide relatively low coverage when compared to bulk sequencing methods. This is expected, as the whole field is still in its infancy beginning only 7 years ago. From this perspective, it is amazing that the single-cell omic sequencing field has already had such great influence and has contributed so tremendously to numerous biological fields. There is huge room for additional development and improvement of the technologies.

Amplification error is a crucial parameter and an issue that limits the accuracy of current single-cell omic sequencing technologies, all of which are based on the pre-amplification of the nucleic acids in individual cells before deep sequencing. After amplification, the single cell being analyzed is already ‘destroyed’; thus, the results cannot be verified in the same individual cell. Some reports use Sanger sequencing to re-sequence the amplified product from the same individual cell for selected loci at which point mutations have been called. Nevertheless, this strategy can detect only the next generation sequencing errors, leaving the single-cell amplification errors concealed and untestable. The other strategy is to use several cells to verify each other and to count only the SNVs that are called in three or more individual cells [92, 96]. Unfortunately, this approach is possible only for cells that can be cultured and amplified at a clonal level in vitro, which is very difficult, if not impossible, for the majority of types of primary cells. If such cell preparation is not possible, this approach will remove the real SNVs that are unique to an individual cell and will severely limit the applications of single-cell omic sequencing technologies. An ideal single-cell genome-sequencing technology would accurately identify both common and ‘private’ SNVs within an individual cell without any false positives resulting from amplification errors. We propose that, in the near future, better single-cell omic sequencing technologies should permit several repeated measurements of the original copy of the nucleic acids within an individual cell. In this way, the amplification errors of sequencing a single cell could be accurately and directly evaluated and determined within the same cell. This would also permit the authentic mutation in an individual cell to be firmly called and verified with essentially no false positives.

It is also important to develop a full set of new bioinformatics tools that are specifically designed for analyses of single-cell omic datasets. These bioinformatics tools should carefully consider both the cons of the single-cell omic datasets, such as high technical noise and high false-negative rates, and the pros of these datasets, such as high sampling numbers and UMIs or spike-in based absolute counting. Despite being valuable, current ‘pseudotime’ analyses have problems in resolving some of the intermediate states during differentiation, especially when these states are dramatically different from both the earlier stem cell state and the later committed state. This is due to the fact that single-cell transcriptome analysis by its nature provides only a snapshot of the gene expression profile for each individual cell, which is an intrinsic disadvantage of this technique when compared with time-lapse imaging methods. As the technique can offer a whole-genome-scale gene expression profile, and because the gene-expression changes in an individual cell at the whole-transcriptome scale can usually be assumed to be ‘continual’ and traceable within a short time interval, one possible resolution for the ‘snapshot’ problem is to sample the cell population much more intensely, ideally every hour or so. Including the following assumption into the pseudotime algorithms is also likely to be helpful: the later differentiation time point will very probably contain differentiation-delayed residual stem cells of the earlier time point. By contrast, the earlier time point is very unlikely to contain fully differentiated cells. For example, during ESC differentiation into liver cells, functional liver cells are very unlikely to be found in the population after just 1 or 2 days; but in the several-week differentiated cell population that contains functional liver cells, it will still be possible to find some residual stem-like cells. Adding this constraint will probably help to resolve the true differentiation pathway of stem cells.

Single-cell multiple omics sequencing technologies have also been developed recently. These methods are capable of simultaneously obtaining information from a single cell on the transcriptome and genome (G&T-seq) [103], or on the transcriptome and DNA methylome (scM&T-seq) [104], or even on all three of these omics (genome, DNA methylome, and transcriptome; scTrio-seq) [73]. New methods covering more layers of different omics are expected to emerge in the near future. These methods are invaluable for elucidating the relationship between different layers of omics in an individual cell. When they become routinely available, permitting the precise recovery of genome, epigenome and transcriptome information from the same individual cell, an ideal approach would be to use single-cell genome sequencing data to perform lineage tracing to reconstruct the pedigree of the cells during stem cell differentiation in vivo. Then, transcriptome data from these cells could be analyzed and used to identify different cell types or sub-populations in the complex tissue. The epigenome information from the same set of single cells could be used subsequently to investigate how different epigenetic layers regulate transcription. Finally, to build a causal relationship between genotype and phenotype, it will be ideal to knockout key component genes for stem cells in vivo using gene-editing technologies. Single-cell multiple omics sequencing at serial time points during the development and differentiation process of stem cells could then be used to reconstruct the core gene regulation network within each individual cell during the differentiation process. The phenotype–genotype relationship for each gene within each individual cell, or between different individual cells, will finally permit us to understand thoroughly the complexity and beauty of the gene regulation network under both physiological and pathological conditions, and will provide us with new insights into the biological basis of human development and diseases.



alveolar type


assay for transposase-accessible chromatin


bisulfite sequencing


chromatin immunoprecipitation


cancer stem cell


pluripotent epiblast


embryonic stem cell


human embryonic stem cell


hematopoietic stem cell


inner cell mass


long-term reconstituting hematopoietic stem cell


multiple annealing and looping-based amplification cycles


multiple displacement amplification


mouse embryonic stem cell


neural stem cell


principal component analysis


polymerase chain reaction


primitive endoderm


primordial germ cell


reduced representation bisulfite sequencing






single-nucleotide variation


short-term reconstituting hematopoietic stem cell




unique molecular identifier


whole-genome bisulfite sequencing


  1. Junker JP, van Oudenaarden A. Every cell is special: genome-wide studies add a new dimension to single-cell biology. Cell. 2014;157:8–11.

    Article  CAS  PubMed  Google Scholar 

  2. Eldar A, Elowitz MB. Functional roles for noise in genetic circuits. Nature. 2010;467:167–73.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  3. Li GW, Xie XS. Central dogma at the single-molecule level in living cells. Nature. 2011;475:308–15.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  4. Huang S. Non-genetic heterogeneity of cells in development: more than just noise. Development. 2009;136:3853–62.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  5. Tang F, Barbacioru C, Wang Y, Nordman E, Lee C, Xu N, et al. mRNA-seq whole-transcriptome analysis of a single cell. Nat Methods. 2009;6:377–82.

    Article  CAS  PubMed  Google Scholar 

  6. Ramskold D, Luo S, Wang YC, Li R, Deng Q, Faridani OR, 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  PubMed Central  CAS  Google Scholar 

  7. 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  CAS  PubMed  Google Scholar 

  8. 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  CAS  PubMed  Google Scholar 

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

    Article  CAS  PubMed  PubMed Central  Google Scholar 

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

    Article  CAS  PubMed  Google Scholar 

  11. Sasagawa Y, Nikaido I, Hayashi T, Danno H, Uno KD, Imai T, et al. 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  PubMed Central  CAS  Google Scholar 

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

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  13. Pan X, Durrett RE, Zhu H, Tanaka Y, Li Y, Zi X, et al. Two methods for full-length RNA sequencing for low quantities of cells and single cells. Proc Natl Acad Sci U S A. 2013;110:594–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Lovatt D, Ruble BK, Lee J, Dueck H, Kim TK, Fisher S, et al. Transcriptome in vivo analysis (TIVA) of spatially defined single cells in live tissue. Nat Methods. 2014;11:190–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Thomsen ER, Mich JK, Yao Z, Hodge RD, Doyle AM, Jang S, et al. Fixed single-cell transcriptomic characterization of human radial glial diversity. Nat Methods. 2016;13:87–93.

    CAS  PubMed  Google Scholar 

  16. Fuzik J, Zeisel A, Máté Z, Calvigioni D, Yanagawa Y, Szabó G, et al. Integration of electrophysiological recordings with single-cell RNA-seq data identifies neuronal subtypes. Nat Biotechnol. 2016;34:175–83.

    Article  CAS  PubMed  Google Scholar 

  17. Cadwell CR, Palasantza A, Jiang X, Berens P, Deng Q, Yilmaz M, et al. Electrophysiological, transcriptomic and morphologic profiling of single neurons using Patch-seq. Nat Biotechnol. 2016;34:199–203.

    Article  CAS  PubMed  Google Scholar 

  18. Wu AR, Neff NF, Kalisky T, Dalerba P, Treutlein B, Rothenberg ME, et al. Quantitative assessment of single-cell RNA-sequencing methods. Nat Methods. 2014;11:41–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Streets AM, Zhang X, Cao C, Pang Y, Wu X, Xiong L, et al. Microfluidic single-cell whole-transcriptome sequencing. Proc Natl Acad Sci U S A. 2014;111:7048–53.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Jaitin DA, Kenigsberg E, Keren-Shaul H, Elefant N, Paul F, Zaretsky I, et al. Massively parallel single-cell RNA-seq for marker-free decomposition of tissues into cell types. Science. 2014;343:776–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Fan HC, Fu GK, Fodor SP. Expression profiling. Combinatorial labeling of single cells for gene expression cytometry. Science. 2015;347:1258367.

    Article  PubMed  CAS  Google Scholar 

  22. Macosko EZ, Basu A, Satija R, Nemesh J, Shekhar K, Goldman M, et al. Highly parallel genome-wide expression profiling of individual cells using nanoliter droplets. Cell. 2015;161:1202–14.

    Article  CAS  PubMed  Google Scholar 

  23. Klein AM, Mazutis L, Akartuna I, Tallapragada N, Veres A, Li V, et al. Droplet barcoding for single-cell transcriptomics applied to embryonic stem cells. Cell. 2015;161:1187–201.

    Article  CAS  PubMed  Google Scholar 

  24. Lee JH, Daugharthy ER, Scheiman J, Kalhor R, Yang JL, Ferrante TC, et al. Highly multiplexed subcellular RNA sequencing in situ. Science. 2014;343:1360–3.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Chen KH, Boettiger AN, Moffitt JR, Wang S, Zhuang X. RNA imaging. Spatially resolved, highly multiplexed RNA profiling in single cells. Science. 2015;348:aaa6090.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  26. Junker JP, Noel ES, Guryev V, Peterson KA, Shah G, Huisken J, et al. Genome-wide RNA tomography in the zebrafish embryo. Cell. 2014;159:662–75.

    Article  CAS  PubMed  Google Scholar 

  27. Satija R, Farrell JA, Gennert D, Schier AF, Regev A. Spatial reconstruction of single-cell gene expression data. Nat Biotechnol. 2015;33:495–502.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Peng G, Suo S, Chen J, Chen W, Liu C, Yu F, et al. Spatial transcriptome for the molecular annotation of lineage fates and cell identity in mid-gastrula mouse embryo. Dev Cell. 2016;36:681–97.

    Article  CAS  PubMed  Google Scholar 

  29. Grun D, van Oudenaarden A. Design and analysis of single-cell sequencing experiments. Cell. 2015;163:799–810.

    Article  CAS  PubMed  Google Scholar 

  30. Shapiro E, Biezuner T, Linnarsson S. Single-cell sequencing-based technologies will revolutionize whole-organism science. Nat Rev Genet. 2013;14:618–30.

    Article  CAS  PubMed  Google Scholar 

  31. Eberwine J, Sul JY, Bartfai T, Kim J. The promise of single-cell sequencing. Nat Methods. 2014;11:25–7.

    Article  CAS  PubMed  Google Scholar 

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

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  33. Pollen AA, Nowakowski TJ, Shuga J, Wang X, Leyrat AA, Lui JH, et al. Low-coverage single-cell mRNA sequencing reveals cellular heterogeneity and activated signaling pathways in developing cerebral cortex. Nat Biotechnol. 2014;32:1053–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Streets AM, Huang Y. How deep is enough in single-cell RNA-seq? Nat Biotechnol. 2014;32:1005–6.

    Article  CAS  PubMed  Google Scholar 

  35. Stegle O, Teichmann SA, Marioni JC. Computational and analytical challenges in single-cell transcriptomics. Nat Rev Genet. 2015;16:133–45.

    Article  CAS  PubMed  Google Scholar 

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

    Article  CAS  PubMed  Google Scholar 

  37. Grun D, Lyubimova A, Kester L, Wiebrands K, Basak O, Sasaki N, et al. Single-cell messenger RNA sequencing reveals rare intestinal cell types. Nature. 2015;525:251–5.

    Article  PubMed  CAS  Google Scholar 

  38. Trapnell C, Cacchiarelli D, Grimsby J, Pokharel P, Li S, Morse M, et al. The dynamics and regulators of cell fate decisions are revealed by pseudotemporal ordering of single cells. Nat Biotechnol. 2014;32:381–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Shin J, Berg DA, Zhu Y, Shin JY, Song J, Bonaguidi MA, et al. Single-cell RNA-seq with waterfall reveals molecular cascades underlying adult neurogenesis. Cell Stem Cell. 2015;17:360–72.

    Article  CAS  PubMed  Google Scholar 

  40. Grun D, Kester L, van Oudenaarden A. Validation of noise models for single-cell transcriptomics. Nat Methods. 2014;11:637–40.

    Article  PubMed  CAS  Google Scholar 

  41. Marinov GK, Williams BA, McCue K, Schroth GP, Gertz J, Myers RM, et al. From single-cell to cell-pool transcriptomes: stochasticity in gene expression and RNA splicing. Genome Res. 2014;24:496–510.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. Brennecke P, Anders S, Kim JK, Kolodziejczyk AA, Zhang X, Proserpio V, et al. Accounting for technical noise in single-cell RNA-seq experiments. Nat Methods. 2013;10:1093–5.

    Article  CAS  PubMed  Google Scholar 

  43. Tang F, Barbacioru C, Nordman E, Bao S, Lee C, Wang X, et al. Deterministic and stochastic allele specific gene expression in single mouse blastomeres. PLoS One. 2011;6:e21208.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  44. Yan L, Yang M, Guo H, Yang L, Wu J, Li R, et al. Single-cell RNA-seq profiling of human preimplantation embryos and embryonic stem cells. Nat Struct Mol Biol. 2013;20:1131–9.

    Article  CAS  PubMed  Google Scholar 

  45. Deng Q, Ramskold D, Reinius B, Sandberg R. Single-cell RNA-seq reveals dynamic, random monoallelic gene expression in mammalian cells. Science. 2014;343:193–6.

    Article  CAS  PubMed  Google Scholar 

  46. Biase F, Cao X, Zhong S. Cell fate inclination within 2-cell and 4-cell mouse embryos revealed by single-cell RNA sequencing. Genome Res. 2014;24:1787–96.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  47. Blakeley P, Fogarty NM, Del Valle I, Wamaitha SE, Hu TX, Elder K, et al. Defining the three cell lineages of the human blastocyst by single-cell RNA-seq. Development. 2015;142:3613.

    Article  PubMed  PubMed Central  Google Scholar 

  48. Xue Z, Huang K, Cai C, Cai L, Jiang CY, Feng Y, et al. Genetic programs in human and mouse early embryos revealed by single-cell RNA sequencing. Nature. 2013;500:593–7.

    Article  CAS  PubMed  Google Scholar 

  49. Treutlein B, Brownfield DG, Wu AR, Neff NF, Mantalas GL, Espinoza FH, et al. Reconstructing lineage hierarchies of the distal lung epithelium using single-cell RNA-seq. Nature. 2014;509:371–5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  50. Hanchate NK, Kondoh K, Lu Z, Kuang D, Ye X, Qiu X, et al. Single-cell transcriptomics reveals receptor transformations during olfactory neurogenesis. Science. 2015;350:1251–5.

    Article  CAS  PubMed  Google Scholar 

  51. Brunskill EW, Park JS, Chung E, Chen F, Magella B, Potter SS. Single cell dissection of early kidney development: multilineage priming. Development. 2014;141:3093–101.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  52. Darmanis S, Sloan SA, Zhang Y, Enge M, Caneda C, Shuer LM, et al. A survey of human brain transcriptome diversity at the single cell level. Proc Natl Acad Sci U S A. 2015;112:7285–90.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  53. Li CL, Li KC, Wu D, Chen Y, Luo H, Zhao JR, et al. Somatosensory neuron types identified by high-coverage single-cell RNA-sequencing and functional heterogeneity. Cell Res. 2016;26:83–102.

    Article  CAS  PubMed  Google Scholar 

  54. Shalek AK, Satija R, Adiconis X, Gertner RS, Gaublomme JT, Raychowdhury R, et al. Single-cell transcriptomics reveals bimodality in expression and splicing in immune cells. Nature. 2013;498:236–40.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  55. Shalek AK, Satija R, Shuga J, Trombetta JJ, Gennert D, Lu D, et al. Single-cell RNA-seq reveals dynamic paracrine control of cellular variation. Nature. 2014;510:363–9.

    CAS  PubMed  PubMed Central  Google Scholar 

  56. Avraham R, Haseley N, Brown D, Penaranda C, Jijon HB, Trombetta JJ, et al. Pathogen cell-to-cell variability drives heterogeneity in host immune responses. Cell. 2015;162:1309–21.

    Article  CAS  PubMed  Google Scholar 

  57. Patel AP, Tirosh I, Trombetta JJ, Shalek AK, Gillespie SM, Wakimoto H, et al. Single-cell RNA-seq highlights intratumoral heterogeneity in primary glioblastoma. Science. 2014;344:1396–401.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  58. Cann GM, Gulzar ZG, Cooper S, Li R, Luo S, Tat M, et al. mRNA-seq of single prostate cancer circulating tumor cells reveals recapitulation of gene expression and pathways found in prostate cancer. PLoS One. 2012;7:e49144.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  59. Miyamoto DT, Zheng Y, Wittner BS, Lee RJ, Zhu H, Broderick KT, et al. RNA-seq of single prostate CTCs implicates noncanonical Wnt signaling in antiandrogen resistance. Science. 2015;349:1351–6.

    Article  CAS  PubMed  Google Scholar 

  60. Kowalczyk MS, Tirosh I, Heckl D, Rao TN, Dixit A, Haas BJ, et al. Single-cell RNA-seq reveals changes in cell cycle and differentiation programs upon aging of hematopoietic stem cells. Genome Res. 2015;25:1860–72.

    Article  PubMed  Google Scholar 

  61. Tsang JC, Yu Y, Burke S, Buettner F, Wang C, Kolodziejczyk AA, et al. Single-cell transcriptomic reconstruction reveals cell cycle and multi-lineage differentiation defects in Bcl11a-deficient hematopoietic stem cells. Genome Biol. 2015;16:178.

    Article  PubMed  PubMed Central  Google Scholar 

  62. Guo F, Yan L, Guo H, Li L, Hu B, Zhao Y, et al. The transcriptome and DNA methylome landscapes of human primordial germ cells. Cell. 2015;161:1437–52.

    Article  CAS  PubMed  Google Scholar 

  63. Llorens-Bobadilla E, Zhao S, Baser A, Saiz-Castro G, Zwadlo K, Martin-Villalba A. Single-cell transcriptomics reveals a population of dormant neural stem cells that become activated upon brain injury. Cell Stem Cell. 2015;17:329–40.

    Article  CAS  PubMed  Google Scholar 

  64. Tang F, Barbacioru C, Bao S, Lee C, Nordman E, Wang X, et al. Tracing the derivation of embryonic stem cells from the inner cell mass by single-cell RNA-seq analysis. Cell Stem Cell. 2010;6:468–78.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  65. Kumar RM, Cahan P, Shalek AK, Satija R, DaleyKeyser AJ, Li H, et al. Deconstructing transcriptional heterogeneity in pluripotent stem cells. Nature. 2014;516:56–61.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  66. Kolodziejczyk AA, Kim JK, Tsang JC, Ilicic T, Henriksson J, Natarajan KN, et al. Single cell RNA-sequencing of pluripotent states unlocks modular transcriptional variation. Cell Stem Cell. 2015;17:471–85.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  67. Kim DH, Marinov GK, Pepke S, Singer ZS, He P, Williams B, et al. Single-cell transcriptome analysis reveals dynamic changes in lncRNA expression during reprogramming. Cell Stem Cell. 2015;16:88–101.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  68. Shi J, Chen Q, Li X, Zheng X, Zhang Y, Qiao J, et al. Dynamic transcriptional symmetry-breaking in pre-implantation mammalian embryo development revealed by single-cell RNA-seq. Development. 2015;142:3468–77.

    Article  PubMed  Google Scholar 

  69. Guo G, Huss M, Tong GQ, Wang C, Li Sun L, Clarke ND, et al. Resolution of cell fate decisions revealed by single-cell gene expression analysis from zygote to blastocyst. Dev Cell. 2010;18:675–85.

    Article  CAS  PubMed  Google Scholar 

  70. Ohnishi Y, Huber W, Tsumura A, Kang M, Xenopoulos P, Kurimoto K, et al. Cell-to-cell expression variability followed by signal reinforcement progressively segregates early mouse lineages. Nat Cell Biol. 2014;16:27–37.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  71. Desai TJ, Brownfield DG, Krasnow MA. Alveolar progenitor and stem cells in lung development, renewal and cancer. Nature. 2014;507:190–4.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  72. Dalerba P, Kalisky T, Sahoo D, Rajendran PS, Rothenberg ME, Leyrat AA, et al. Single-cell dissection of transcriptional heterogeneity in human colon tumors. Nat Biotechnol. 2011;29:1120–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  73. Hou Y, Guo H, Cao C, Li X, Hu B, Zhu P, et al. Single-cell triple omics sequencing reveals genetic, epigenetic, and transcriptomic heterogeneity in hepatocellular carcinomas. Cell Res. 2016;26:304–19.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  74. Guo H, Zhu P, Wu X, Li X, Wen L, Tang F. Single-cell methylome landscapes of mouse embryonic stem cells and early embryos analyzed using reduced representation bisulfite sequencing. Genome Res. 2013;23:2126–35.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  75. Smallwood SA, Lee HJ, Angermueller C, Krueger F, Saadeh H, Peat J, et al. Single-cell genome-wide bisulfite sequencing for assessing epigenetic heterogeneity. Nat Methods. 2014;11:817–20.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  76. Farlik M, Sheffield NC, Nuzzo A, Datlinger P, Schonegger A, Klughammer J, et al. Single-cell DNA methylome sequencing and bioinformatic inference of epigenomic cell-state dynamics. Cell Rep. 2015;10:1386–97.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  77. Gu H, Smith ZD, Bock C, Boyle P, Gnirke A, Meissner A. Preparation of reduced representation bisulfite sequencing libraries for genome-scale DNA methylation profiling. Nat Protoc. 2011;6:468–81.

    Article  CAS  PubMed  Google Scholar 

  78. Miura F, Enomoto Y, Dairiki R, Ito T. Amplification-free whole-genome bisulfite sequencing by post-bisulfite adaptor tagging. Nucleic Acids Res. 2012;40:e136.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  79. Guo H, Zhu P, Yan L, Li R, Hu B, Lian Y, et al. The DNA methylation landscape of human early embryos. Nature. 2014;511:606–10.

    Article  CAS  PubMed  Google Scholar 

  80. Smith ZD, Chan MM, Humm KC, Karnik R, Mekhoubad S, Regev A, et al. DNA methylation dynamics of the human preimplantation embryo. Nature. 2014;511:611–5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  81. Tang WW, Dietmann S, Irie N, Leitch HG, Floros VI, Bradshaw CR, et al. A unique gene regulatory network resets the human germline epigenome for development. Cell. 2015;161:1453–67.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  82. Gkountela S, Zhang KX, Shafiq TA, Liao WW, Hargan-Calvopina J, Chen PY, et al. DNA demethylation dynamics in the human prenatal germline. Cell. 2015;161:1425–36.

    Article  CAS  PubMed  Google Scholar 

  83. Xu GL, Walsh CP. Enzymatic DNA oxidation: mechanisms and biological significance. BMB Rep. 2014;47:609–18.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  84. Buenrostro JD, Wu B, Litzenburger UM, Ruff D, Gonzales ML, Snyder MP, et al. Single-cell chromatin accessibility reveals principles of regulatory variation. Nature. 2015;523:486–90.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  85. Cusanovich DA, Daza R, Adey A, Pliner HA, Christiansen L, Gunderson KL, et al. Epigenetics. Multiplex single-cell profiling of chromatin accessibility by combinatorial cellular indexing. Science. 2015;348:910–4.

    Article  CAS  PubMed  Google Scholar 

  86. Jin W, Tang Q, Wan M, Cui K, Zhang Y, Ren G, et al. Genome-wide detection of DNase I hypersensitive sites in single cells and FFPE tissue samples. Nature. 2015;528:142–6.

    CAS  PubMed  Google Scholar 

  87. Nagano T, Lubling Y, Stevens TJ, Schoenfelder S, Yaffe E, Dean W, et al. Single-cell Hi-C reveals cell-to-cell variability in chromosome structure. Nature. 2013;502:59–64.

    Article  CAS  PubMed  Google Scholar 

  88. Kind J, Pagie L, de Vries SS, Nahidiazar L, Dey SS, Bienko M, et al. Genome-wide maps of nuclear lamina interactions in single human cells. Cell. 2015;163:134–47.

    Article  CAS  PubMed  Google Scholar 

  89. Rotem A, Ram O, Shoresh N, Sperling RA, Goren A, Weitz DA, et al. Single-cell ChIP-seq reveals cell subpopulations defined by chromatin state. Nat Biotechnol. 2015;33:1165–72.

    Article  CAS  PubMed  Google Scholar 

  90. 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  CAS  PubMed  Google Scholar 

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

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  92. 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  CAS  PubMed  PubMed Central  Google Scholar 

  93. Wang J, Fan HC, Behr B, Quake SR. Genome-wide single-cell analysis of recombination activity and de novo mutation rates in human sperm. Cell. 2012;150:402–12.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  94. Gole J, Gore A, Richards A, Chiu YJ, Fung HL, Bushman D, et al. Massively parallel polymerase cloning and genome sequencing of single cells using nanoliter microwells. Nat Biotechnol. 2013;31:1126–32.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  95. Fu Y, Li C, Lu S, Zhou W, Tang F, Xie XS, et al. Uniform and accurate single-cell sequencing based on emulsion whole-genome amplification. Proc Natl Acad Sci U S A. 2015;112:11923–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  96. Wang Y, Waters J, Leung ML, Unruh A, Roh W, Shi X, et al. Clonal evolution in breast cancer revealed by single nucleus genome sequencing. Nature. 2014;512:155–60.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

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

    Article  PubMed  PubMed Central  CAS  Google Scholar 

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

    Article  CAS  PubMed  Google Scholar 

  99. Gawad C, Koh W, Quake SR. Single-cell genome sequencing: current state of the science. Nat Rev Genet. 2016;17:175–88.

    Article  CAS  PubMed  Google Scholar 

  100. Lu S, Zong C, Fan W, Yang M, Li J, Chapman AR, et al. Probing meiotic recombination and aneuploidy of single sperm cells by whole-genome sequencing. Science. 2012;338:1627–30.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  101. Hou Y, Fan W, Yan L, Li R, Lian Y, Huang J, et al. Genome analyses of single human oocytes. Cell. 2013;155:1492–506.

    Article  CAS  PubMed  Google Scholar 

  102. Behjati S, Huch M, van Boxtel R, Karthaus W, Wedge DC, Tamuri AU, et al. Genome sequencing of normal cells reveals developmental lineages and mutational processes. Nature. 2014;513:422–5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  103. Macaulay IC, Haerty W, Kumar P, Li YI, Hu TX, Teng MJ, et al. G&T-seq: parallel sequencing of single-cell genomes and transcriptomes. Nat Methods. 2015;12:519–22.

    Article  CAS  PubMed  Google Scholar 

  104. Angermueller C, Clark SJ, Lee HJ, Macaulay IC, Teng MJ, Hu TX, et al. Parallel single-cell sequencing links transcriptional and epigenetic heterogeneity. Nat Methods. 2016;13:229–32.

    Article  CAS  PubMed  Google Scholar 

Download references


We apologize for the relevant papers not cited here due to strict space limitation.



Author information

Authors and Affiliations


Corresponding authors

Correspondence to Lu Wen or Fuchou Tang.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

Both authors wrote, read and approved the final manuscript.

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

Wen, L., Tang, F. Single-cell sequencing in stem cell biology. Genome Biol 17, 71 (2016).

Download citation

  • Published:

  • DOI: