Abstract
Microbiologists traditionally study population rather than individual cells, as it is generally assumed that the status of individual cells will be similar to that observed in the population. However, the recent studies have shown that the individual behavior of each single cell could be quite different from that of the whole population, suggesting the importance of extending traditional microbiology studies to single-cell level. With recent technological advances, such as flow cytometry, next-generation sequencing (NGS), and microspectroscopy, single-cell microbiology has greatly enhanced the understanding of individuality and heterogeneity of microbes in many biological systems. Notably, the application of multiple ‘omics’ in single-cell analysis has shed light on how individual cells perceive, respond, and adapt to the environment, how heterogeneity arises under external stress and finally determines the fate of the whole population, and how microbes survive under natural conditions. As single-cell analysis involves no axenic cultivation of target microorganism, it has also been demonstrated as a valuable tool for dissecting the microbial ‘dark matter.’ In this review, current state-of-the-art tools and methods for genomic and transcriptomic analysis of microbes at single-cell level were critically summarized, including single-cell isolation methods and experimental strategies of single-cell analysis with NGS. In addition, perspectives on the future trends of technology development in the field of single-cell analysis was also presented.
Introduction
Microbiologists usually study microorganisms by deciphering their physiology, internal interactions, and even genetic information. Traditionally, these studies are all carried out at the population level, typically using millions to billions of cells for analysis in bulk, and assuming the status of individual cells is similar to that observed in the population. Although these results are, no doubt, informative, they often neglect any heterogeneity that is possibly present in the population. Meanwhile, the recent studies have shown that cell-to-cell heterogeneity at both cellular and molecular levels in isogenic population could be an order of magnitude greater than previously thought (Lidstrom and Meldrum, 2003), suggesting the importance of extending traditional microbiology studies to the single-cell level. It is now increasingly accepted that conclusions based on conventional average molecular or phenotypic measurements of a population could be biased, as the patterns of distinct sub-populations cannot be revealed (Wang et al., 2015).
Heterogeneities could result from either phenotypic difference between isogenic cells or genetic diversity at population level (Davis and Isberg, 2016). Mechanisms responsible for the cell-to-cell variation could be classified into four categories: stochastic gene expression, phenotypic plasticity, genotypic plasticity, and reversible genotypic variation (Roberfroid et al., 2016). While stochastic gene expression and phenotypic plasticity only lead to phenotypic differences, genotypic plasticity and reversible genotypic variation could introduce heterogeneity to an isogenic population at the genetic level. Stochastic gene expression widely exists in both prokaryotic and eukaryotic populations and is not exclusively driven by genomic information. Noise is one of the mechanisms of stochastic variability, which could be independent of environmental signals. This variation, either triggered by intrinsic or extrinsic noise, is usually unimodal. However, a unimodal noisiness of gene expression may trigger a bimodal behavior of downstream gene expression, leading to bimodality and bistability in a population (Dubnau and Losick, 2006; Veening et al., 2008). As different subpopulations co-exist, some individual cells may express genes that allow them to survive stresses prior to environmental changes. By using this mechanism, microbes could ensure that some individuals will survive under harsh conditions (Veening et al., 2008). Phenotypic plasticity is a kind of environmental-driven viability and could make it possible for cells to adapt to the fluctuations in the environment (Viney and Reece, 2013). An example is that in heterogeneous environmental conditions such as biofilm, isogenic cells could differentiate into various phenotypes and form several sub-populations for adapting to their local environmental conditions (van Gestel et al., 2015). Genotypic plasticity usually occurs in populations subjected to the experimental evolution. Driven by clonal evolution, clonal cells could evolve and finally result in genotypic diversification (Korona et al., 1994; Rainey and Travisano, 1998). Various mechanisms of genotypic diversification, such as clonal interference (Barrick and Lenski, 2013), niche construction, and niche partitioning (Barrick and Lenski, 2009), have been reported for both well-mixed and spatially structured environments (Roberfroid et al., 2016). As the final cause of variation, reversible genotypic variations are driven by random site-specific recombination, gene conversion, or epigenetic modification, leading to phase variations that play important roles in the virulence of some pathogens, and causing increased heterogeneity in the population (Davis and Isberg, 2016; Roberfroid et al., 2016).
Another major shortage of traditional microbiology approaches is the dependence on establishing laboratory culture for studying targeted microbes. Meanwhile, it is well-known that so far only a small number of microbial species in the biosphere could be cultivated successfully in the laboratory, leaving a great deal of microbial information untouched (Cardenas and Tiedje, 2008; Rinke et al., 2013). The hidden information, also known as microbial ‘dark matter,’ has drawn great interests recently and provides potential solutions for several critical issues, such as new drugs and antibiotics discovery (Ling et al., 2015), toxic chemicals degradation (Jiang et al., 2016), understanding pathogen virulence and disease mechanisms (Omsland et al., 2009), and revealing the human microbiome (Browne et al., 2016). Although obtaining axenic culture from natural isolates remains important, it is usually labor-intensive (Connon and Giovannoni, 2002), having a low success rate, and might be biased (Wu et al., 2009). In addition, comparing with the axenic cultures in the laboratory, microorganisms usually live in a more complex and barren environment in nature, making it unable to present the original state of microorganisms in the laboratory (Stewart, 2012). In recent years, many attempts have been employed for analyzing the microbe without axenic culture. For example, metagenomics and metatranscriptomics have been widely used for studying microbial community (Venter et al., 2004; Tringe et al., 2005; Mason et al., 2012, 2014; Meng et al., 2014). However, metagenomics and metatranscriptomics are not well-suited to reveal unambiguous information about the organization of discovered genes within genomes, evolutionary histories of specific organisms, and in situ interactions among organisms (Yoon et al., 2011; Stepanauskas, 2012). Genomic information, such as genome rearrangements, gene insertions, duplications and loss, is hard to obtain from metagenomic analysis since the assembled results could be mosaics of DNA from cells sharing high-homology regions but vary in genome-wide similarity (Stepanauskas, 2012).
Single-cell analysis can be effective for addressing these issues and providing better and in-depth understanding of the status of microbial cells. As it starts from only one cell, single-cell analysis could reveal information about individual cell without laboratory cultivation. With the help of high-throughput sequencing, it is possible to obtain functional genomics information of each single cell in its natural environment, so that its original genetic and functional status in a complex community can be revealed globally, quantitatively, and absolutely. Several reports using single-cell analysis have successfully revealed information like coexisting subpopulations, organismal interactions, and new metabolic pathways from uncultivated samples, which could hardly be obtained by traditional approaches (Marcy et al., 2007b; Hess et al., 2011; Siegl et al., 2011; Yoon et al., 2011; Martinez-Garcia et al., 2012a; Kashtan et al., 2014). In recent years, significant progress has been made to apply metagenomics and metatranscriptomics approaches to reveal the genetic information and gene expression patterns of cells in a population, and uncover microbial species and gene diversity in a community (Bowler et al., 2009). However, as metagenomics and metatranscriptomics could not reveal the information such as repetitive regions or strain heterogeneity comprehensively in a complex population, single-cell-based analysis has been proposed as a valuable supplement to the efficient identification of novel microbial species and the accurate interpretation of the metagenomics and metatranscriptomics results (Massana et al., 2014; Vannier et al., 2016; Ji et al., 2017). In this review, we summarize current state-of-the-art tools and methods for genomic and transcriptomic analysis of microbes at single-cell level, including single-cell isolation methods and experimental strategies of single-cell analysis with NGS, and provide some perspectives on the future trends of technology development in single-cell analysis field.
Tools for Single-Cell Isolation
Single-cell isolation is the very first step in the single-cell analysis process (Figure 1). The major challenge of this step is: how to isolate cells of interest accurately in a high-throughput manner and without causing any genetic or physiological change to the target cells. Basically, methods being applied for microbial single-cell isolation could be classified as two principal approaches: micromanipulation and random encapsulation (Blainey, 2013). Micromanipulation methods, including micropipette and optical tweezer approaches, are carried out under high-resolution microscope. These methods offer a great confidence that every single cell can be observed, captured and delivered to the next step. Traditional micropipette method could be easily applied in any laboratory on an inverted microscope with mechanical liquid handling. Although very labor-consuming and low-throughput, approximately in the order of 50 cells/h and person (Picelli, 2016), it is the first choice if only a small number of cells are required for the next step analysis (Qi et al., 2014, 2016; Wang et al., 2015). In addition, commercial robotic manipulation system for automated single-cell selection has also been developed and applied for microbial single-cell analysis (Anis et al., 2008; Merza et al., 2009; Gao et al., 2011; Banerjee et al., 2014), making it possible for relatively high-throughput single-cell isolation. Optical tweezer approaches are implemented by tightly focusing a laser beam for trapping cells in solution (Ashkin et al., 1986, 1987; Ashkin and Dziedzic, 1987). Usually, by using near-infrared wavelengths of light, cells could be easily handled without any harm (Neuman et al., 1999; Ericsson et al., 2000). This method has been successfully applied in many microbial isolation experiments, including filamentous bacteria (Pamp et al., 2012) and even virus (Ashkin and Dziedzic, 1987). More detailed reviews of this method have been published recently (Moffitt et al., 2008; Hashemi Shabestari et al., 2017), interested readers could refer to these articles.
FIGURE 1
Flow cytometry and microfluidic device are the most widely used random encapsulation approaches in recent years. Flow cytometry and FACS have a much higher throughput, and have been demonstrated as an effective platform for single-cell analysis in microbial cells (Raghunathan et al., 2005; Stepanauskas and Sieracki, 2007; Woyke et al., 2009; Swan et al., 2011; Martinez-Garcia et al., 2012b; Field et al., 2015). Although FACS is fundamentally based on random encapsulation, flow cytometers can monitor several parameters of individuals, which means single cells can be sorted according to their size, morphology, spontaneous fluorescence, fluorescence-labeled antibodies and staining dyes simultaneously, making it possible to sort even rare cell types. In addition, it is easy to sort single cells directly into 96- and 384-well plates using commercial instruments, which means single-cell analysis workflow could be entirely performed using automated liquid handling robots. However, cells are typically subjected to physical stresses during the sorting, such as fluidic pressure, laser beam, electrostatic charges, voltage fields, and collisions with container surfaces, which could significantly affect the cell physiology and even the recovery rate during the cultivation (Marie et al., 2017). In the case when the sorted cells are used for gene expression or transcriptome analysis, proper RNA protectant needs to be added (Qi et al., 2014, 2016; Wang et al., 2015); while in the case when the sorted cells are used for clonal cultivation, extra efforts to carefully optimize the cultivation conditions are necessary in order to maximize the success rates (Marie et al., 2017). Under the conception of ‘Lab-on-a-chip,’ microfluidic devices have become the most popular method for single-cell isolation. With these devices, researchers could integrate single-cell analysis process from cell isolation to sequencing library preparation in only a coin-sized microchip, which could be either purchased from commercial manufacturers or designed and fabricated using materials such as PDMS in the laboratory. Combining with detection technologies, such as fluorescence spectroscopy (Wolff et al., 2003) or raman spectroscopy (Song et al., 2016), microfluidic devices could perform specific sorting while encapsulating cells with reagents for cell lysis and sequencing library preparation at nanoliter volume with high-throughput (Klein et al., 2015). Comparing with traditional tube-based reactions, microfluidic devices require few manual liquid handling, leading to a significant decrease in contamination and less variations among samples (de Bourcy et al., 2014). Notably, less contamination with nanoliter reaction volume means a higher concentration of substrates, resulting in better uniformity amplification (Fu et al., 2015; Leung et al., 2016). In addition, comparing with other methods, microfluidic devices cause less physical stresses to cells, leading to more accurate physiological analysis and high success rates of further cultivation analysis (Jiang et al., 2016; Kim et al., 2017; Song et al., 2017; Zhang et al., 2017). In general, both flow cytometry and microfluidic devices could provide high-throughput and accurate single-cell sorting. Flow cytometers could monitor multiple parameters and are capable for rare cell detection and sorting, but are usually expensive and require skilled operators. Before sorting, it is also necessary to prepare a sterile system for flow cytometers to prevent contamination. By contrast, microfluidic devices can be designed and made in the laboratory. They can be disposable in order to minimize contamination, and are easy to operate. Besides single-cell sorting, microfluidic devices could offer integrative single-cell analysis including cell culture and tracking (Yu et al., 2017), digital PCR (Ottesen et al., 2006), and sequencing library preparation (Hosokawa et al., 2017; Lan et al., 2017). With several advantages mentioned above, in recent years, microfluidic devices tended to be used as an analytic platform rather than just an isolation method for single-cell analysis (Marshall et al., 2012; Zhang et al., 2015, 2017; Jiang et al., 2016; Haliburton et al., 2017; Hosokawa et al., 2017; Kim et al., 2017; Lan et al., 2017; Shahi et al., 2017; Song et al., 2017). Several detailed reviews have been published recently on microfluidic devices (Wen et al., 2016; Caen et al., 2017; Prakadan et al., 2017; Xi et al., 2017), and interested readers could refer to these articles.
Tools for Genomic Analysis at Single-Cell Level
A single microbial cell usually contains picogram to femtogram level of genomic DNA (Kim et al., 2017). Sequencing technologies, up to now, are still unable to sequence such a low amount of nucleic acids directly without any amplification. Therefore, researches have been applying WGA methods since 1990 (Lichter et al., 1990). However, as amplification is conducted using DNA polymerases, the amplified products could contain genetic information of the original cell as well as some artifacts, such as genome fragment loss, amplification bias, mutations, and chimeras. Over the past 20 years, WGA methods have been optimized with substantial progress, including less contamination and better amplification performance (Blainey, 2013). In general, amplification methods could be classified into three categories: pure PCR-based amplification, isothermal amplification, and hybrid methods (Gawad et al., 2016).
Pure PCR-based WGA methods are the primary methods at early stage in the single-cell genomic analysis. Early approaches with specific primers, such as linker-adapter (also known as ligation-anchored) PCR (LA-PCR) (Troutt et al., 1992; Klein et al., 1999) and IRS-PCR (Lengauer et al., 1990; Lichter et al., 1990), require ligation reaction or prior knowledge of the target sequence. Later, methods with random primers, including primer extension pre-amplification PCR (PEP-PCR) (Hubert et al., 1992) and DOP-PCR (Telenius et al., 1992) were introduced. As the most representative method in this category, DOP-PCR typically contains two stages, with the first facilitating random primer extension on the template genome DNA and the second favoring amplicon replication with specific primer (Telenius et al., 1992).
The second category of WGA is isothermal amplification, which was first reported in Dean et al. (2001); Zhang et al. (2001) and has been demonstrated as a powerful tool in microbial single-cell genomic analysis, especially with MDA (Lasken, 2012). Unlike PCR-based methods, isothermal amplification methods use polymerases with strong strand displacement activity, such as φ29 polymerase, and 6-mer 3′-protected random primers for isothermal extension (Dean et al., 2001; Zhang et al., 2001). During extension, polymerase creates and displaces synthesized products from single-stranded DNA template, and the displaced DNA is the template for further priming and synthesis (Dean et al., 2001; Zhang et al., 2001). Compared with PCR-based methods, MDA shows higher genome coverage, lower error rates and much longer extension length over 10,000 nt (Blanco et al., 1989; de Bourcy et al., 2014). However, the loci amplified first are typically found to be overrepresented, indicating non-uniformity of MDA (de Bourcy et al., 2014). Recently, a novel primer-free method called TruePrime was reported and has been successfully used for the amplification of genomic DNA from single human HEK293 cells (Picher et al., 2016). In this method, an enzyme called TthPrimPol, which has a wide range of template specificity, serves as primase for φ29 polymerase mediated MDA. During the reaction, TthPrimPol binds to the denatured DNA and synthesizes short DNA primers. The DNA primers are recognized and extended by φ29 polymerase. Then, TthPrimPol catalyzes new rounds of priming on the elongated single-strand DNA, followed by further rounds of strand-displacement synthesis and resulting in exponential amplification (Picher et al., 2016). Another isothermal amplification method, called SPIA, could achieve linear amplification under isothermal conditions by using a specific DNA/RNA hybrid primer, together with RNase H and a strand-displacing DNA polymerase (Kurn et al., 2005). In SPIA method, strand-displacement only occurs at the DNA/RNA hybrid primer site of the amplicons, preventing the exponential amplification in MDA. Recently, a new method called LIANTI was reported (Chen et al., 2017). As an isothermal amplification method, this approach depends on RNA polymerase but not DNA polymerase for linear amplification. In this method, genomic DNA from a single cell was fragmented and tagged by Tn5 transposon with a T7 promoter, then linear amplified with T7 RNA polymerase, and finally converted to DNA by reverse transcription for further library preparation (Chen et al., 2017).
Two similar hybrid methods, displacement DOP-PCR (D-DOP-PCR, also known as PicoPLEX or GenomePlex) (Langmore, 2002) and MALBAC (Lu et al., 2012; Zong et al., 2012), were recently developed to overcome the low coverage of PCR-based methods and the non-uniformity of MDA. These two methods both use isothermal amplification followed by PCR amplification, but different primers for extension. D-DOP-PCR uses degenerated primers in the first step adding an anchor sequence with isothermal amplification and then using PCR amplification for the second step (Langmore, 2002). MALBAC, however, uses a random primer with a designed anchor which could promote looping of the isothermal amplification products to prevent further amplification before the second PCR step, suggesting a more uniform amplification (Lu et al., 2012; Zong et al., 2012).
In practice, isothermal and hybrid methods are currently the most commonly used approaches, as they show better performance comparing with pure PCR-based methods. Several groups have compared these methods using both microbial and mammalian cells (Chen et al., 2014; de Bourcy et al., 2014; Deleye et al., 2015; Hou et al., 2015; Ning et al., 2015). These reports have drawn similar conclusions that MDA has significantly higher genome coverage breadth and lower false-positive rates, while hybrid methods demonstrate better coverage uniformity (Chen et al., 2014; de Bourcy et al., 2014; Deleye et al., 2015; Hou et al., 2015; Ning et al., 2015). For example, one report showed that MDA has better coverage breadth than MALBAC (84% vs. 52%), resulting in higher detection rates of SNVs (88% vs. 52%) in human cells (Hou et al., 2015). Another report showed that hybrid methods has better coverage uniformity than MDA, suggesting that hybrid methods have better performance than MDA in detecting CNVs (Ning et al., 2015). In the report, the researchers also found that MALBAC tended to over-amplify genomic regions with a high-GC content (Ning et al., 2015). The average GC content of amplified DNA using GenomePlex (41.6%) was very close to the reference genome (41.9%), while the average GC contents of amplified DNA regions by MDA and MALBAC were 43.4 and 46.6%, respectively (Ning et al., 2015). However, after a GC-correction, the correlation of read abundance between MALBAC and bulk-cell samples (R2 = 0.53) was nearly the same as GenomePlex (R2 = 0.56), while MDA gave a very poor correlation (R2 = 0.02) (Ning et al., 2015). The TruePrime method was reported to have better coverage uniformity than the primer-based MDA, leading to an improved CNV detection accuracy, thus an advantage over the traditional primer-based MDA protocol (Picher et al., 2016). In addition, by using human genomic DNA as input, TruePrime could amplify as low as 1 fg DNA, which is about 100-fold more sensitive than the primer-based MDA (Picher et al., 2016). This superior sensitivity could be very valuable for microbial single-cell genomic analysis, as most microbes are much smaller and contain less DNA than eukaryotic cells. Notably, the most recent LIANTI method exhibited significantly improved amplification uniformity and genome coverage over the previous methods on all scales, and was capable for both high accuracy of CNV detection and low SNV false-positive rate (Chen et al., 2017). As new invented approaches, both TruePrime and LIANTI have the potential but still need more evaluation to demonstrate their performance on microbial single-cell genomic analysis. In conclusion, there is no clear winner in performance between MDA and hybrid methods yet, and researchers should choose methods depending on the metric of their interest (Gawad et al., 2016). As microbial single-cell analysis usually focuses on elucidating the genomic information of the microbial ‘dark matter,’ genome coverage is the key to be concerned. Therefore, MDA method has been far more widely used for microbial single-cell analysis rather than the others.
Besides amplification methods, previous reports also found that by using microfluidic devices, microbial single-cell analysis could obtain a better performance comparing with the traditional tube-based approach (de Bourcy et al., 2014). With higher mapping ratio and better repeatability, microfluidic devices could also reduce the contamination especially from the experiment operator (de Bourcy et al., 2014). Recently, two groups independently reported high-throughput microbial single-cell analysis protocols based on self-designed microfluidic devices (Hosokawa et al., 2017; Lan et al., 2017). These two protocols shared some similarities in single microbe encapsulation and lysis protocols. However, one protocol involved sorting the positive amplification droplets and re-amplification of the DNA for further analysis by NGS and qPCR (Hosokawa et al., 2017), while the other protocol used a strategy of labeling DNA fragments from the same cell with a barcode, and then pooling and sequencing of the barcoded fragments of all cells (Lan et al., 2017). These protocols could provide reliable pipelines for analyzing 10s of 1000s of single microbial cells within a couple of hours with a comparable performance to the conventional techniques. In conclusion, with further improvements on both amplification methods and microfluidic devices, microbial single-cell genomic analysis will be more efficient, reliable, and convenient in the near future.
Tools for Transcriptomic Analysis at Single-Cell Level
Prior to whole-genome transcriptomic analysis, relative quantification methodologies have been developed to measure expression of small number of genes at single-cell level. For example, methods using fluorescent reporter proteins coupling with high-throughput data acquisition approaches such as flow cytometry have been widely applied for detecting gene expression heterogeneities within the microbial population (Taniguchi et al., 2010; Roberfroid et al., 2016). In addition, methods using RT-qPCR for detecting gene expression in single cells have also been reported and successfully applied to several types of microbes for heterogeneity analysis (Gao et al., 2011; Shi et al., 2013; Qi et al., 2014, 2016; Wang et al., 2015; Thompson et al., 2017; Turkarslan et al., 2017). However, these methods could only reveal gene expression patterns of a very limited number of genes, while not able to uncover global information in a cell. Moreover, application of such approaches typically requires genetic engineering tools and genomic information of the target microorganisms, limiting the application to only model organisms.
Global transcriptomic analysis could circumvent the above drawbacks and even possible for unknown species without genome information using de novo NGS approach. Compared to genomic analysis, transcriptomic analysis for microbes at single-cell level is much more challenging for several reasons. First, microbial cells usually contain picogram to femtogram level RNA molecules (de Bekker et al., 2011; Kang et al., 2011; Wang et al., 2015), while mammalian cells could have up to nanogram level RNA molecules (Picelli, 2016). Besides the low-RNA content, RNA molecules of prokaryotic cells are less stable than DNA and could be degraded by ribonucleases that are widely existing and hard to be deactivated. Moreover, rRNA and tRNA molecules usually represent over 90% of total RNA, but offer limited biological information and should be excluded in the amplification process, as most researches focus on mRNA and other rare molecules. With a complicated cell wall, harsher conditions are typically required to lyze a microbial cell, which may lead to damage or loss of RNA, and accuracy and efficiency of the downstream transcriptomic analysis (Khan and Yadav, 2004; Hall et al., 2013; Heera et al., 2015; He et al., 2016). More importantly, unlike genomic analysis, in which the methods for mammalian cells could be also readily applied to prokaryotic microbes, not all methods for mammalian single-cell transcriptomic analysis could be used to microbes. This is simply because of the structure differences of mRNA molecules between eukaryotic and prokaryotic cells. Currently, most of the mammalian single-cell transcriptomic analysis approaches use oligo(dT) primers in the first cDNA synthesis step. This is based on the 3′ poly(A) structure of mRNA molecules from eukaryotic cells, which makes them easier to be enriched from rRNA and tRNA. However, mRNA molecules from prokaryotic cells usually lack the poly(A) tail, and require random primers for cDNA synthesis. By using random primers, both rRNA and tRNA will also be included in the resulting transcriptome library, thus being sequenced together, leading to a low coverage of the target mRNA. In addition, application of random primers for cDNA synthesis causes losses of 3′ sequence information, as they are usually unable to obtain the full-length transcripts. Hence, so far only a few reports on prokaryotic single-cell transcriptomic analysis have been reported (Kang et al., 2011, 2015; Wang et al., 2015). Even for eukaryotic microbes that could be analyzed with well-developed approaches for mammalian cells, only one report analyzing single-cell transcriptomics of neighboring hyphae of Aspergillus niger was reported (de Bekker et al., 2011). Concerning this circumstance, we summarized below all the state-of-the-art tools in single-cell transcriptomic analysis and discussed possibilities for their microbial applications, especially for prokaryotic microorganisms.
To our knowledge, the earliest study of single-cell transcriptomics was reported in Eberwine et al. (1992). In this work, mRNA molecules from single-living neurons were reverse transcribed to cDNA using oligo(dT)-T7 primer. Then, the synthesized double-stranded cDNA molecules with T7 promoter were used as templates for IVT with T7 RNA polymerase for producing amplified RNA. Next, the amplified RNA molecules were used as templates for the second turn of reverse transcription. After this process, over a million-fold amplification of the original RNA was achieved. Although this report only used ISH for accessing gene expression, it reveals the possibility about analyzing gene expression at a single-cell level. Based on the concept of this study, several studies have successfully analyzed the whole transcriptome of single mammalian cells (Morris et al., 2011; Hashimshony et al., 2012, 2016; Jaitin et al., 2014).
In the past decade, several new approaches were developed, leading to a tremendous progress in mammalian single-cell RNA-seq (Tang et al., 2009, 2010; Islam et al., 2011, 2012; Goetz and Trimarchi, 2012; Hashimshony et al., 2012, 2016; Picelli et al., 2013, 2014; Sasagawa et al., 2013; Jaitin et al., 2014; Soumillon et al., 2014; Fan H.C. et al., 2015; Klein et al., 2015; Macosko et al., 2015). The most widely used single-cell RNA-seq methods are characterized in Table 1. As these methods have been well reviewed in several excellent articles (Saliba et al., 2014; Chen et al., 2015; Kolodziejczyk et al., 2015; Picelli, 2016), we will focus only on some newly developed methods here. Among the methods listed in Table 1, Smart-seq/Smart-seq2 and Quartz-seq use a method called ‘template switch’ for the second strand cDNA synthesis, generating full-length double-stranded cDNA comparing with the others (Goetz and Trimarchi, 2012; Picelli et al., 2013, 2014; Sasagawa et al., 2013). Smart-seq, CEL-seq2 and STRT-seq are compatible with Fluidigm C1 Single-Cell Auto Prep system, which is an automated platform and captured using integrated fluidic circuits (Ziegenhain et al., 2017). For amplification types, Tang’s method, Smart-seq/Smart-seq2, Quartz-seq, and STRT-seq are all based on PCR amplification, while CEL-seq/CEL-seq2 and MARS-seq are based on IVT. The advantage of IVT is that the amplification efficiency is sequence independent. However, as it requires a second time of reverse transcription, there is 3′ coverage bias of the sequencing results (Kolodziejczyk et al., 2015). When choosing an appropriate single-cell RNA-seq method, transcript coverage, strand specificity, position bias, and UMI compatibility should be concerned depending on the purpose of the research (Table 1). For example, full-length RNA-seq methods such as Smart-seq/Smart-Seq2, and Quartz-seq could sequence the transcripts in their entirety, and thus are suggested for de novo sequencing and the detection of SNPs and mutations. However, these methods are not compatible with strand-specific protocol and UMI. Methods such as STRT-seq, CEL-seq/CEL-seq2, MARS-seq, and Drop-seq are all compatible with strand-specific protocol and UMI, although they tend to be 5′ or 3′ end biased. While UMI approach in single-cell sequencing could reduce amplification noise and provide more accurate expression quantification, strand-specific sequencing could provide more information for antisense transcript discovery, genome annotation, and expression profiling. In conclusion, current single-cell RNA-seq methods are still facing a trade-off between coverage and uniformity (Picelli, 2016). Notably, none of these methods has been evaluated in eukaryotic microbes, suggesting further optimization and development are needed for microbial cells. Recently, a systematically evaluation of six prominent single-cell RNA-seq methods has been reported (Ziegenhain et al., 2017), and the results indicated that Smart-seq2 had the best coverage because of its full-length synthesis ability. However, as Smart-seq2 is incompatible with UMIs, all methods using UMIs have less amplification noise. In addition, power simulations showed that Drop-seq is more cost-efficient for analyzing a large number of cells, while Smart-seq2, MARS-seq, and SCRB-seq are more efficient with the analysis of a small number of cells (Ziegenhain et al., 2017).
Table 1
| Name | Transcript coverage | Position bias | Strand specificity | UMI compatible | Key reference |
|---|---|---|---|---|---|
| Tang’s method | Nearly full-length | Strongly 3′ | No | No | Tang et al., 2009, 2010 |
| Quartz-seq | Full-length | Weakly 3′ | No | No | Sasagawa et al., 2013 |
| Smart-seq/Smart-seq2 | Full-length | Weakly 3′ | No | No | Goetz and Trimarchi, 2012; Picelli et al., 2013, 2014 |
| STRT-seq | 5′ only | 5′ only | Yes | Yes | Islam et al., 2011, 2012 |
| CEL-seq/CEL-seq2 | 3′ only | 3′ only | Yes | Yes | Hashimshony et al., 2012, 2016 |
| MARS-seq | 3′ only | 3′ only | Yes | Yes | Jaitin et al., 2014 |
| SCRB-seq | 3′ only | 3′ only | Yes | Yes | Soumillon et al., 2014 |
| Drop-seq/InDrop | 3′ only | 3′ only | Yes | Yes | Klein et al., 2015; Macosko et al., 2015 |
| Cyto-seq | Pre-defined genes only | 3′ only | Yes | Yes | Fan H.C. et al., 2015 |
Characteristics of several widely used single-cell RNA-seq methods.
Besides the above methods, several new methods have also been developed recently. Some of them have already been utilized for single-cell RNA-seq in prokaryotic cells (Kang et al., 2011, 2015; Wang et al., 2015). The first case of single-cell microbial transcriptomic analysis, to our knowledge, was reported in Kang et al. (2011). In this report, transcriptome of single bacterium Burkholderia thailandensis was analyzed using microarray through amplification of RNA molecules by rolling circle amplification. In this report, bacterial cells were first lysed with Triton X-100 and lysozyme, and then the lysate was used for direct cDNA synthesis with random primers. After genomic DNA degradation, single-stranded cDNA molecules were self-ligated and then used as the template for multiply primed rolling circle amplification using φ29 polymerase with random primers. The result showed low fold-change bias and only less than 6% drop-outs with no contamination. In addition, this method also preferred an optional rRNA/tRNA elimination step for deep sequencing. By using 5′-phosphate-dependent exonuclease, rRNA and tRNA molecules, which have the 5′-phosphate structure, will be specifically degraded, leaving the mRNA molecules which lacking the 5′-phosphate structure for the next cDNA synthesis step. This is also the only report we could find, which has successfully depleted rRNA from single microbial cells, indicating the requirements for further innovation of other effective rRNA depletion and mRNA enrichment methods for microbes. Later in the same year, another approach using Ribo-SPIA method, that is, derived from SPIA method for amplification, has successfully been employed to analyze transcriptomics of neighboring hyphae of the eukaryotic fungus A. niger using microarray (de Bekker et al., 2011). In this article, total RNA from different 5 hyphal tips were isolated using a column based kit and amplified using the WT-Ovation One-Direct RNA Amplification System (Nugen) with both oligo(dT) and random primers. Microarray analysis resulted in a present call for 4–7% of the A. niger genes, of which 12% showed heterogeneous RNA levels, indicating the feasibility of using this method for microbial transcriptomic analysis. In another study with prokaryotic cells, Wang et al. (2015) successfully conducted single-cell RNA-seq in single cyanobacterium Synechocystis sp. PCC 6803 cells with Ribo-SPIA method. To determine the heterogeneity upon environmental stress, this method was applied to Synechocystis single cells at 24 and 72 h after nitrogen starvation treatment. With up to 98% of all putative Synechocystis genes identified in single cells, a possible increasing gene-expression heterogeneity from 24 to 72 h after nitrogen starvation stress was also found, indicating the method could achieve good identification of the transcripts in single bacterial cells (Wang et al., 2015). More recently, a technology for targeted depletion of abundant transcripts was developed by Nugen (Armour et al., 2015). Unlike the exonuclease-based depletion method that Kang et al. (2011) reported, this method depleted the unwanted sequences after cDNA synthesis using probes that target unwanted sequences. However, the information of the unwanted sequences is required, making it impossible for de novo single-cell RNA-seq. SUPeR-seq (Fan X. et al., 2015) is another method to sequence both polyadenylated and non-polyadenylated RNAs, suggesting its possible application to prokaryotic microbes. This method shares some similarities to Tang’s method (Tang et al., 2009, 2010), but used a primer containing an anchor sequence (AnchorX), 15-mer dT sequence and 6-mer random sequence for simultaneous detection of both polyadenylated and non-polyadenylated RNA molecules and synthesizing the first strand cDNA. After poly(A) tailing for the first strand cDNA, a primer containing another anchor sequence (AnchorY) and 24-mer dT sequence was used for second strand cDNA synthesis, and then the double-stranded cDNA molecules were amplified by PCR using AnchorX and AnchorY primers. With this approach, the researchers discovered 2891 circRNAs in mouse preimplantation embryos. Like other methods using random primers, rRNA could not be excluded with this method. However, this method provides another possible method for single-cell RNA-seq in prokaryotic microbes, especially with the rRNA depletion methods mentioned above.
Future Perspectives
Current genomic and transcriptomic analysis of single microbial cells share several similar challenges. Cell lysis is a major challenge for single-cell analysis. As microbes typically contain complicated structure of cell walls, appropriate lysis strategies need to be chosen carefully without damaging the DNA/RNA inside. In addition, in the case if the lysate is directly used for amplification without purification, the lysis condition should also be carefully optimized to minimize the influence of lysis related reagents to the downstream reactions. Alternatively, a method called FluidFM might be a promising approach for DNA/RNA isolation from microbial cells, as it used a ‘nanosyringe’ to extract cytoplasmic and nucleoplasmic fractions from single live cells rather than lysis the cell (Meister et al., 2009; Guillaume-Gentil et al., 2016). Contamination is another key challenge in single-cell analysis. As low-input and high-fold amplification are required for sequencing, single-cell analysis is very sensitive to contamination, either from the laboratory environment or reagents and instruments used for sample preparation. Several approaches have been applied to minimize contamination, including reducing the reaction volume of lysis and amplification reaction to nanoliter scale in a sealed, disposable microfluidic device (Marcy et al., 2007a,b), using UV exposure to inactivate contaminates in reagents (Zhang et al., 2006; Woyke et al., 2011), and disposable plasticware produced from virgin materials (Blainey and Quake, 2011). Another challenge for microbial single-cell analysis is the ultra-low nuclei acids content in a single microbial cell. Current microbial single-cell sequencing methods were all modified from those developed for mammalian cells, as they contain more nuclei acids. While using these methods in microbial cells, nuclei acids template could be a 1000-fold less than using a mammalian cell. With a much lower concentration of templates, the amplification process could be more sensitive to any contamination and non-specific amplification. In addition, less input may also challenge the sensitivity of the polymerase used for the amplification process (Picher et al., 2016). Using microfluidic devices for amplification could significantly solve these problems (de Bourcy et al., 2014). Moreover, the low input also influences the uniformity of the amplification. Even for single-cell analysis of mammalian cells, the amplification uniformity is still not comparable with that at the bulk-cell level. Therefore, the sequencing depth could be a critical factor to ensure good genome coverage, especially for unculturable microbes with unknown genome sizes.
Current single-cell sequencing methods all require amplification of DNA/RNA from a single cell for NGS sequencing, which will inevitably introduce bias and loss. Recently, new sequencing platforms such as true single molecule sequencing (tSMS, Helicos, now SeqLL), SMRT sequencing (PacBio), and nanopore sequencing (Oxford Nanopore) could sequence DNA/RNA molecules at single-molecule level and prove to be possible to sequence DNA/RNA molecules directly from bulk-cells without pre-amplification (Ozsolak et al., 2009; Ozsolak and Milos, 2011; Coupland et al., 2012; Ayub et al., 2013). Although directly sequencing a single cell without pre-amplification is still challenging, further innovation of these new technologies and sequencing platforms could eventually make it possible for single-cell analysis without any amplification.
Conclusion
As a rapidly growing field, single-cell analysis plays a significant role in extending our understanding of microorganisms by revealing how individual cells perceive, respond and adapt to the environment, and determine the fate of the whole population. The key drivers of new technology for single-cell analysis will be advancement in throughput, integration of isolation and amplification, and integrated analysis with multiple ‘omics.’ Even with many challenges still ahead, we believe that this field will receive a tremendous boost with progress of several related fields, such as microfluidic devices and new sequencing platforms.
Statements
Author contributions
WZ envisioned this project. ZC, LC, and WZ wrote the manuscript. All authors have read and agreed on the manuscript.
Acknowledgments
The research was supported by grants from the Natural Science Foundation of China (NSFC) (Nos. 31470217, 91751102, and 21621004), National Basic Research Program of China (National ‘973’ program, project No. 2014CB745101), and the Tianjin Municipal Science and Technology Commission (No. 15JCZDJC32500).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Abbreviations
- CNV
copy-number variation
- D-DOP-PCR
displacement DOP-PCR
- DOP-PCR
degenerate oligonucleotide-primed PCR
- FACS
fluorescence-activated cell sorting
- FluidFM
fluidic force microscope
- IPS-PCR
interspersed repetitive sequence PCR
- ISH
in situ hybridization
- IVT
in vitro transcription
- LA-PCR
linker-adapter or ligation-anchored PCR
- LIANTI
linear amplification via transposon insertion
- MALBAC
multiple annealing and looping-based amplification cycles
- MDA
multiple displacement amplification
- NGS
next-generation sequencing
- PDMS
polydimethylsiloxane
- PEP-PCR
preamplification PCR
- poly(A)
polyadenylated
- RNA-seq
RNA sequencing
- SMRT
single molecule real-time
- SNV
single nucleotide variant
- SPIA
single primer isothermal amplification
- tSMS
true single molecule sequencing
- UMIs
unique molecular identifiers
- WGA
whole genomic amplification.
References
1
AnisY. H.HollM. R.MeldrumD. R. (2008). “Automated vision-based selection and placement of single cells in microwell array formats,” inIEEE International Conference on Automation Science and Engineering, Xian. 10.1109/COASE.2008.4626575
2
ArmourC.AmoreseD.LiB.KurnN. (2015). Compositions and methods for negative selection of non-desired nucleic acid sequences. U.S. Patent No 20150299767. Washington, DC: U.S. Patent and Trademark Office.
3
AshkinA.DziedzicJ. M. (1987). Optical trapping and manipulation of viruses and bacteria.Science2351517–1520. 10.1126/science.3547653
4
AshkinA.DziedzicJ. M.BjorkholmJ. E.ChuS. (1986). Observation of a single-beam gradient force optical trap for dielectric particles.Opt. Lett.11:288. 10.1364/OL.11.000288
5
AshkinA.DziedzicJ. M.YamaneT. (1987). Optical trapping and manipulation of single cells using infrared laser beams.Nature330769–771. 10.1038/330769a0
6
AyubM.HardwickS. W.LuisiB. F.BayleyH. (2013). Nanopore-based identification of individual nucleotides for direct RNA sequencing.Nano Lett.136144–6150. 10.1021/nl403469r
7
BanerjeeS.WangW.GutschmidtS. (2014). “Design and analysis of actuated microneedles for robotic manipulation,” inASME International Mechanical Engineering Congress and Exposition, Montreal, QC, 8–13. 10.1115/imece2014-39308
8
BarrickJ. E.LenskiR. E. (2009). Genome-wide mutational diversity in an evolving population of Escherichia coli.Cold Spring Harb. Symp. Quant. Biol.74119–129. 10.1101/sqb.2009.74.018
9
BarrickJ. E.LenskiR. E. (2013). Genome dynamics during experimental evolution.Nat. Rev. Genet.14827–839. 10.1038/nrg3564
10
BlaineyP. C. (2013). The future is now: single-cell genomics of bacteria and archaea.FEMS Microbiol. Rev.37407–427. 10.1111/1574-6976.12015
11
BlaineyP. C.QuakeS. R. (2011). Digital MDA for enumeration of total nucleic acid contamination.Nucleic Acids Res.39 e19. 10.1093/nar/gkq1074
12
BlancoL.BernadA.LazaroJ. M.MartinG.GarmendiaC.SalasM. (1989). Highly efficient DNA synthesis by the phage phi 29 DNA polymerase. Symmetrical mode of DNA replication.J. Biol. Chem.2648935–8940.
13
BowlerC.KarlD. M.ColwellR. R. (2009). Microbial oceanography in a sea of opportunity.Nature459180–184. 10.1038/nature08056
14
BrowneH. P.ForsterS. C.AnonyeB. O.KumarN.NevilleB. A.StaresM. D.et al (2016). Culturing of ’unculturable’ human microbiota reveals novel taxa and extensive sporulation.Nature533543–546. 10.1038/nature17645
15
CaenO.LuH.NizardP.TalyV. (2017). Microfluidics as a strategic player to decipher single-cell omics?Trends Biotechnol.35713–727. 10.1016/j.tibtech.2017.05.004
16
CardenasE.TiedjeJ. M. (2008). New tools for discovering and characterizing microbial diversity.Curr. Opin. Biotechnol.19544–549. 10.1016/j.copbio.2008.10.010
17
ChenC.XingD.TanL.LiH.ZhouG.HuangL.et al (2017). Single-cell whole-genome analyses by Linear Amplification via Transposon Insertion (LIANTI).Science356189–194. 10.1126/science.aak9787
18
ChenK.DongX.WuJ. (2015). “The application of single-cell sequencing in dynamic transcriptomes,” in Single Cell Sequencing and Systems Immunology, ed.WangX. (Dordrecht: Springer), 41–63. 10.1007/978-94-017-9753-5_4
19
ChenM.SongP.ZouD.HuX.ZhaoS.GaoS.et al (2014). Comparison of multiple displacement amplification (MDA) and multiple annealing and looping-based amplification cycles (MALBAC) in single-cell sequencing.PLOS ONE9:e114520. 10.1371/journal.pone.0114520
20
ConnonS. A.GiovannoniS. J. (2002). High-throughput methods for culturing microorganisms in very-low-nutrient media yield diverse new marine isolates.Appl. Environ. Microbiol.683878–3885. 10.1128/AEM.68.8.3878-3885.2002
21
CouplandP.ChandraT.QuailM.ReikW.SwerdlowH. (2012). Direct sequencing of small genomes on the Pacific Biosciences RS without library preparation.Biotechniques53365–372. 10.2144/000113962
22
DavisK. M.IsbergR. R. (2016). Defining heterogeneity within bacterial populations via single cell approaches.Bioessays38782–790. 10.1002/bies.201500121
23
de BekkerC.BruningO.JonkerM. J.BreitT. M.WostenH. A. (2011). Single cell transcriptomics of neighboring hyphae of Aspergillus niger.Genome Biol.12:R71. 10.1186/gb-2011-12-8-r71
24
de BourcyC. F.De VlaminckI.KanbarJ. N.WangJ.GawadC.QuakeS. R. (2014). A quantitative comparison of single-cell whole genome amplification methods.PLOS ONE9:e105585. 10.1371/journal.pone.0105585
25
DeanF. B.NelsonJ. R.GieslerT. L.LaskenR. S. (2001). Rapid amplification of plasmid and phage DNA using Phi 29 DNA polymerase and multiply-primed rolling circle amplification.Genome Res.111095–1099. 10.1101/gr.180501
26
DeleyeL.De ConinckD.ChristodoulouC.SanteT.DheedeneA.HeindryckxB.et al (2015). Whole genome amplification with SurePlex results in better copy number alteration detection using sequencing data compared to the MALBAC method.Sci. Rep.5:11711. 10.1038/srep11711
27
DubnauD.LosickR. (2006). Bistability in bacteria.Mol. Microbiol.61564–572. 10.1111/j.1365-2958.2006.05249.x
28
EberwineJ.YehH.MiyashiroK.CaoY.NairS.FinnellR.et al (1992). Analysis of gene expression in single live neurons.Proc. Natl. Acad. Sci. U.S.A.893010–3014. 10.1073/pnas.89.7.3010
29
EricssonM.HanstorpD.HagbergP.EngerJ.NystromT. (2000). Sorting out bacterial viability with optical tweezers.J. Bacteriol.1825551–5555. 10.1128/jb.182.19.5551-5555.2000
30
FanH. C.FuG. K.FodorS. P. (2015). Expression profiling. Combinatorial labeling of single cells for gene expression cytometry.Science347:1258367. 10.1126/science.1258367
31
FanX.ZhangX.WuX.GuoH.HuY.TangF.et al (2015). Single-cell RNA-seq transcriptome analysis of linear and circular RNAs in mouse preimplantation embryos.Genome Biol.16148. 10.1186/s13059-015-0706-1
32
FieldE. K.SczyrbaA.LymanA. E.HarrisC. C.WoykeT.StepanauskasR.et al (2015). Genomic insights into the uncultivated marine Zetaproteobacteria at Loihi Seamount.ISME J.9857–870. 10.1038/ismej.2014.183
33
FuY.LiC.LuS.ZhouW.TangF.XieX. S.et al (2015). Uniform and accurate single-cell sequencing based on emulsion whole-genome amplification.Proc. Natl. Acad. Sci. U.S.A.11211923–11928. 10.1073/pnas.1513988112
34
GaoW.ZhangW.MeldrumD. R. (2011). RT-qPCR based quantitative analysis of gene expression in single bacterial cells.J. Microbiol. Methods85221–227. 10.1016/j.mimet.2011.03.008
35
GawadC.KohW.QuakeS. R. (2016). Single-cell genome sequencing: current state of the science.Nat. Rev. Genet.17175–188. 10.1038/nrg.2015.16
36
GoetzJ. J.TrimarchiJ. M. (2012). Transcriptome sequencing of single cells with Smart-Seq.Nat. Biotechnol.30763–765. 10.1038/nbt.2325
37
Guillaume-GentilO.GrindbergR. V.KoogerR.Dorwling-CarterL.MartinezV.OssolaD.et al (2016). Tunable single-cell extraction for molecular analyses.Cell166506–516. 10.1016/j.cell.2016.06.025
38
HaliburtonJ. R.ShaoW.DeutschbauerA.ArkinA.AbateA. R. (2017). Genetic interaction mapping with microfluidic-based single cell sequencing.PLOS ONE12:e0171302. 10.1371/journal.pone.0171302
39
HallE. W.KimS.AppadooV.ZareR. N. (2013). Lysis of a single cyanobacterium for whole genome amplification.Micromachines4321–332. 10.3390/mi4030321
40
Hashemi ShabestariM.MeijeringA. E.RoosW. H.WuiteG. J.PetermanE. J. (2017). Recent advances in biological single-molecule applications of optical tweezers and fluorescence microscopy.Methods Enzymol.58285–119. 10.1016/bs.mie.2016.09.047
41
HashimshonyT.SenderovichN.AvitalG.KlochendlerA.de LeeuwY.AnavyL.et al (2016). CEL-Seq2: sensitive highly-multiplexed single-cell RNA-Seq.Genome Biol.1777. 10.1186/s13059-016-0938-8
42
HashimshonyT.WagnerF.SherN.YanaiI. (2012). CEL-Seq: single-cell RNA-Seq by multiplexed linear amplification.Cell Rep.2666–673. 10.1016/j.celrep.2012.08.003
43
HeJ.DuS.TanX.ArefinA.HanC. S. (2016). Improved lysis of single bacterial cells by a modified alkaline-thermal shock procedure.Biotechniques60129–135. 10.2144/000114389
44
HeeraR.SivachandranP.ChinniS. V.MasonJ.CroftL.RavichandranM.et al (2015). Efficient extraction of small and large RNAs in bacteria for excellent total RNA sequencing and comprehensive transcriptome analysis.BMC Res. Notes8:754. 10.1186/s13104-015-1726-3
45
HessM.SczyrbaA.EganR.KimT. W.ChokhawalaH.SchrothG.et al (2011). Metagenomic discovery of biomass-degrading genes and genomes from cow rumen.Science331463–467. 10.1126/science.1200387
46
HosokawaM.NishikawaY.KogawaM.TakeyamaH. (2017). Massively parallel whole genome amplification for single-cell sequencing using droplet microfluidics.Sci. Rep.75199. 10.1038/s41598-017-05436-4
47
HouY.WuK.ShiX.LiF.SongL.WuH.et al (2015). Comparison of variations detection between whole-genome amplification methods used in single-cell resequencing.Gigascience437. 10.1186/s13742-015-0068-3
48
HubertR.WeberJ. L.SchmittK.ZhangL.ArnheimN. (1992). A new source of polymorphic DNA markers for sperm typing: analysis of microsatellite repeats in single cells.Am. J. Hum. Genet.51985–991.
49
IslamS.KjallquistU.MolinerA.ZajacP.FanJ. B.LonnerbergP.et al (2011). Characterization of the single-cell transcriptional landscape by highly multiplex RNA-seq.Genome Res.211160–1167. 10.1101/gr.110882.110
50
IslamS.KjallquistU.MolinerA.ZajacP.FanJ. B.LonnerbergP.et al (2012). Highly multiplexed and strand-specific single-cell RNA 5’ end sequencing.Nat. Protoc.7813–828. 10.1038/nprot.2012.022
51
JaitinD. A.KenigsbergE.Keren-ShaulH.ElefantN.PaulF.ZaretskyI.et al (2014). Massively parallel single-cell RNA-seq for marker-free decomposition of tissues into cell types.Science343776–779. 10.1126/science.1247651
52
JiP.ZhangY.WangJ.ZhaoF. (2017). MetaSort untangles metagenome assembly by reducing microbial community complexity.Nat. Commun.8:14306. 10.1038/ncomms14306
53
JiangC. Y.DongL.ZhaoJ. K.HuX.ShenC.QiaoY.et al (2016). High throughput single-cell cultivation on microfluidic streak plates.Appl. Environ. Microbiol.822210–2218. 10.1128/AEM.03588-15
54
KangY.McMillanI.NorrisM. H.HoangT. T. (2015). Single prokaryotic cell isolation and total transcript amplification protocol for transcriptomic analysis.Nat. Protoc.10974–984. 10.1038/nprot.2015.058
55
KangY.NorrisM. H.Zarzycki-SiekJ.NiermanW. C.DonachieS. P.HoangT. T. (2011). Transcript amplification from single bacterium for transcriptome analysis.Genome Res.21925–935. 10.1101/gr.116103.110
56
KashtanN.RoggensackS. E.RodrigueS.ThompsonJ. W.BillerS. J.CoeA.et al (2014). Single-cell genomics reveals hundreds of coexisting subpopulations in wild Prochlorococcus.Science344416–420. 10.1126/science.1248575
57
KhanI. U. H.YadavJ. S. (2004). Development of a single-tube, cell lysis-based, genus-specific PCR method for rapid identification of Mycobacteria: optimization of cell lysis, PCR primers and conditions, and restriction rattern analysis.J. Clin. Microbiol.42453–457. 10.1128/jcm.42.1.453-457.2004
58
KimS.De JongheJ.KulesaA. B.FeldmanD.VatanenT.BhattacharyyaR. P.et al (2017). High-throughput automated microfluidic sample preparation for accurate microbial genomics.Nat. Commun.8:13919. 10.1038/ncomms13919
59
KleinA. M.MazutisL.AkartunaI.TallapragadaN.VeresA.LiV.et al (2015). Droplet barcoding for single-cell transcriptomics applied to embryonic stem cells.Cell1611187–1201. 10.1016/j.cell.2015.04.044
60
KleinC. A.Schmidt-KittlerO.SchardtJ. A.PantelK.SpeicherM. R.RiethmullerG. (1999). Comparative genomic hybridization, loss of heterozygosity, and DNA sequence analysis of single cells.Proc. Natl. Acad. Sci. U.S.A.964494–4499. 10.1073/pnas.96.8.4494
61
KolodziejczykA. A.KimJ. K.SvenssonV.MarioniJ. C.TeichmannS. A. (2015). The technology and biology of single-cell RNA sequencing.Mol. Cell58610–620. 10.1016/j.molcel.2015.04.005
62
KoronaR.NakatsuC. H.ForneyL. J.LenskiR. E. (1994). Evidence for multiple adaptive peaks from populations of bacteria evolving in a structured habitat.Proc. Natl. Acad. Sci. U.S.A.919037–9041. 10.1073/pnas.91.19.9037
63
KurnN.ChenP.HeathJ. D.Kopf-SillA.StephensK. M.WangS. (2005). Novel isothermal, linear nucleic acid amplification systems for highly multiplexed applications.Clin. Chem.511973–1981. 10.1373/clinchem.2005.053694
64
LanF.DemareeB.AhmedN.AbateA. R. (2017). Single-cell genome sequencing at ultra-high-throughput with microfluidic droplet barcoding.Nat. Biotechnol.35640–646. 10.1038/nbt.3880
65
LangmoreJ. P. (2002). Rubicon Genomics, Inc.Pharmacogenomics3557–560. 10.1517/14622416.3.4.557
66
LaskenR. S. (2012). Genomic sequencing of uncultured microorganisms from single cells.Nat. Rev. Microbiol.10631–640. 10.1038/nrmicro2857
67
LengauerC.RiethmanH.CremerT. (1990). Painting of human chromosomes with probes generated from hybrid cell lines by PCR with Alu and L1 primers.Hum. Genet.861–6. 10.1007/BF00205163
68
LeungK.KlausA.LinB. K.LaksE.BieleJ.LaiD.et al (2016). Robust high-performance nanoliter-volume single-cell multiple displacement amplification on planar substrates.Proc. Natl. Acad. Sci. U.S.A.1138484–8489. 10.1073/pnas.1520964113
69
LichterP.LedbetterS. A.LedbetterD. H.WardD. C. (1990). Fluorescence in situ hybridization with Alu and L1 polymerase chain reaction probes for rapid characterization of human chromosomes in hybrid cell lines.Proc. Natl. Acad. Sci. U.S.A.876634–6638. 10.1073/pnas.87.17.6634
70
LidstromM. E.MeldrumD. R. (2003). Life-on-a-chip.Nat. Rev. Microbiol.1158–164. 10.1038/nrmicro755
71
LingL. L.SchneiderT.PeoplesA. J.SpoeringA. L.EngelsI.ConlonB. P.et al (2015). A new antibiotic kills pathogens without detectable resistance.Nature517455–459. 10.1038/nature14098
72
LuS.ZongC.FanW.YangM.LiJ.ChapmanA. R.et al (2012). Probing meiotic recombination and aneuploidy of single sperm cells by whole-genome sequencing.Science3381627–1630. 10.1126/science.1229112
73
MacoskoE. Z.BasuA.SatijaR.NemeshJ.ShekharK.GoldmanM.et al (2015). Highly parallel genome-wide expression profiling of individual cells using nanoliter droplets.Cell1611202–1214. 10.1016/j.cell.2015.05.002
74
MarcyY.IshoeyT.LaskenR. S.StockwellT. B.WalenzB. P.HalpernA. L.et al (2007a). Nanoliter reactors improve multiple displacement amplification of genomes from single cells.PLOS Genet. 3:1702–1708. 10.1371/journal.pgen.0030155
75
MarcyY.OuverneyC.BikE. M.LosekannT.IvanovaN.MartinH. G.et al (2007b). Dissecting biological “dark matter” with single-cell genetic analysis of rare and uncultivated TM7 microbes from the human mouth.Proc. Natl. Acad. Sci. U.S.A.10411889–11894. 10.1073/pnas.0704662104
76
MarieD.Le GallF.EdernR.GourvilP.VaulotD. (2017). Improvement of phytoplankton culture isolation using single cell sorting by flow cytometry.J. Phycol.53271–282. 10.1111/jpy.12495
77
MarshallI. P.BlaineyP. C.SpormannA. M.QuakeS. R. (2012). A single-cell genome for Thiovulum sp.Appl. Environ. Microbiol.788555–8563. 10.1128/AEM.02314-12
78
Martinez-GarciaM.BrazelD. M.SwanB. K.ArnostiC.ChainP. S.ReitengaK. G.et al (2012a). Capturing single cell genomes of active polysaccharide degraders: an unexpected contribution of Verrucomicrobia.PLOS ONE7:e35314. 10.1371/journal.pone.0035314
79
Martinez-GarciaM.SwanB. K.PoultonN. J.GomezM. L.MaslandD.SierackiM. E.et al (2012b). High-throughput single-cell sequencing identifies photoheterotrophs and chemoautotrophs in freshwater bacterioplankton.ISME J.6113–123. 10.1038/ismej.2011.84
80
MasonO. U.HazenT. C.BorglinS.ChainP. S.DubinskyE. A.FortneyJ. L.et al (2012). Metagenome, metatranscriptome and single-cell sequencing reveal microbial response to Deepwater Horizon oil spill.ISME J.61715–1727. 10.1038/ismej.2012.59
81
MasonO. U.ScottN. M.GonzalezA.Robbins-PiankaA.BaelumJ.KimbrelJ.et al (2014). Metagenomics reveals sediment microbial community response to Deepwater Horizon oil spill.ISME J.81464–1475. 10.1038/ismej.2013.254
82
MassanaR.del CampoJ.SierackiM. E.AudicS.LogaresR. (2014). Exploring the uncultured microeukaryote majority in the oceans: reevaluation of ribogroups within stramenopiles.ISME J.8854–866. 10.1038/ismej.2013.204
83
MeisterA.GabiM.BehrP.StuderP.VorosJ.NiedermannP.et al (2009). FluidFM: combining atomic force microscopy and nanofluidics in a universal liquid delivery system for single cell applications and beyond.Nano Lett.92501–2507. 10.1021/nl901384x
84
MengJ.XuJ.QinD.HeY.XiaoX.WangF. (2014). Genetic and functional properties of uncultivated MCG archaea assessed by metagenome and gene expression analyses.ISME J.8650–659. 10.1038/ismej.2013.174
85
MerzaS. A.PosnerJ. D.AnisY. H.YoungA. C.JohnsonR. H.MeldrumD. R. (2009). “Microfluidic device for transport and observation of single cells,” in ASME 2009 International Mechanical Engineering Congress and Exposition, New York City, NY. 10.1115/IMECE2009-13019
86
MoffittJ. R.ChemlaY. R.SmithS. B.BustamanteC. (2008). Recent advances in optical tweezers.Annu. Rev. Biochem.77205–228. 10.1146/annurev.biochem.77.043007.090225
87
MorrisJ.SinghJ. M.EberwineJ. H. (2011). Transcriptome analysis of single cells.J. Vis. Exp.2011:e2634. 10.3791/2634
88
NeumanK. C.ChaddE. H.LiouG. F.BergmanK.BlockS. M. (1999). Characterization of photodamage to Escherichia coli in optical traps.Biophys. J.772856–2863. 10.1016/S0006-3495(99)77117-1
89
NingL.LiZ.WangG.HuW.HouQ.TongY.et al (2015). Quantitative assessment of single-cell whole genome amplification methods for detecting copy number variation using hippocampal neurons.Sci. Rep.5:11415. 10.1038/srep11415
90
OmslandA.CockrellD. C.HoweD.FischerE. R.VirtanevaK.SturdevantD. E.et al (2009). Host cell-free growth of the Q fever bacterium Coxiella burnetii.Proc. Natl. Acad. Sci. U.S.A.1064430–4434. 10.1073/pnas.0812074106
91
OttesenE. A.HongJ. W.QuakeS. R.LeadbetterJ. R. (2006). Microfluidic digital PCR enables multigene analysis of individual environmental bacteria.Science3141464–1467. 10.1126/science.1131370
92
OzsolakF.MilosP. M. (2011). Single-molecule direct RNA sequencing without cDNA synthesis.Wiley Interdiscip. Rev. RNA2565–570. 10.1002/wrna.84
93
OzsolakF.PlattA. R.JonesD. R.ReifenbergerJ. G.SassL. E.McInerneyP.et al (2009). Direct RNA sequencing.Nature461814–818. 10.1038/nature08390
94
PampS. J.HarringtonE. D.QuakeS. R.RelmanD. A.BlaineyP. C. (2012). Single-cell sequencing provides clues about the host interactions of segmented filamentous bacteria (SFB).Genome Res.221107–1119. 10.1101/gr.131482.111
95
PicelliS. (2016). Single-cell RNA-sequencing: the future of genome biology is now.RNA Biol.14637–650. 10.1080/15476286.2016.1201618
96
PicelliS.BjorklundA. K.FaridaniO. R.SagasserS.WinbergG.SandbergR. (2013). Smart-seq2 for sensitive full-length transcriptome profiling in single cells.Nat. Methods101096–1098. 10.1038/nmeth.2639
97
PicelliS.FaridaniO. R.BjorklundA. K.WinbergG.SagasserS.SandbergR. (2014). Full-length RNA-seq from single cells using Smart-seq2.Nat. Protoc.9171–181. 10.1038/nprot.2014.006
98
PicherA. J.BudeusB.WafzigO.KrugerC.Garcia-GomezS.Martinez-JimenezM. I.et al (2016). TruePrime is a novel method for whole-genome amplification from single cells based on TthPrimPol.Nat. Commun.7:13296. 10.1038/ncomms13296
99
PrakadanS. M.ShalekA. K.WeitzD. A. (2017). Scaling by shrinking: empowering single-cell ’omics’ with microfluidic devices.Nat. Rev. Genet.18345–361. 10.1038/nrg.2017.15
100
QiZ.PeiG.ChenL.ZhangW. (2014). Single-cell analysis reveals gene-expression heterogeneity in syntrophic dual-culture of Desulfovibrio vulgaris with Methanosarcina barkeri.Sci. Rep.4:7478. 10.1038/srep07478
101
QiZ. H.ChenL.ZhangW. W. (2016). Comparison of transcriptional heterogeneity of eight genes between batch Desulfovibrio vulgaris biofilm and planktonic culture at a single-cell level.Front. Microbiol.7:597. 10.3389/fmicb.2016.00597
102
RaghunathanA.FergusonH. R.Jr.BornarthC. J.SongW.DriscollM.LaskenR. S. (2005). Genomic DNA amplification from a single bacterium.Appl. Environ. Microbiol.713342–3347. 10.1128/AEM.71.6.3342-3347.2005
103
RaineyP. B.TravisanoM. (1998). Adaptive radiation in a heterogeneous environment.Nature39469–72. 10.1038/27900
104
RinkeC.SchwientekP.SczyrbaA.IvanovaN. N.AndersonI. J.ChengJ. F.et al (2013). Insights into the phylogeny and coding potential of microbial dark matter.Nature499431–437. 10.1038/nature12352
105
RoberfroidS.VanderleydenJ.SteenackersH. (2016). Gene expression variability in clonal populations: causes and consequences.Crit. Rev. Microbiol.42969–984. 10.3109/1040841X.2015.1122571
106
SalibaA. E.WestermannA. J.GorskiS. A.VogelJ. (2014). Single-cell RNA-seq: advances and future challenges.Nucleic Acids Res.428845–8860. 10.1093/nar/gku555
107
SasagawaY.NikaidoI.HayashiT.DannoH.UnoK. D.ImaiT.et al (2013). Quartz-Seq: a highly reproducible and sensitive single-cell RNA sequencing method, reveals non-genetic gene-expression heterogeneity.Genome Biol.14:R31. 10.1186/gb-2013-14-4-r31
108
ShahiP.KimS. C.HaliburtonJ. R.GartnerZ. J.AbateA. R. (2017). Abseq: ultrahigh-throughput single cell protein profiling with droplet microfluidic barcoding.Sci. Rep.7:44447. 10.1038/srep44447
109
ShiX.GaoW.ChaoS. H.ZhangW.MeldrumD. R. (2013). Monitoring the single-cell stress response of the diatom Thalassiosira pseudonana by quantitative real-time reverse transcription-PCR.Appl. Environ. Microbiol.791850–1858. 10.1128/AEM.03399-12
110
SieglA.KamkeJ.HochmuthT.PielJ.RichterM.LiangC.et al (2011). Single-cell genomics reveals the lifestyle of Poribacteria, a candidate phylum symbiotically associated with marine sponges.ISME J.561–70. 10.1038/ismej.2010.95
111
SongY.KasterA. K.VollmersJ.SongY.DavisonP. A.FrentrupM.et al (2017). Single-cell genomics based on Raman sorting reveals novel carotenoid-containing bacteria in the Red Sea.Microb. Biotechnol.10125–137. 10.1111/1751-7915.12420
112
SongY.YinH.HuangW. E. (2016). Raman activated cell sorting.Curr. Opin. Chem. Biol.331–8. 10.1016/j.cbpa.2016.04.002
113
SoumillonM.CacchiarelliD.SemrauS.van OudenaardenA.MikkelsenT. S. (2014). Characterization of directed differentiation by high-throughput single-cell RNA-Seq.bioRxiv.10.1101/003236
114
StepanauskasR. (2012). Single cell genomics: an individual look at microbes.Curr. Opin. Microbiol.15613–620. 10.1016/j.mib.2012.09.001
115
StepanauskasR.SierackiM. E. (2007). Matching phylogeny and metabolism in the uncultured marine bacteria, one cell at a time.Proc. Natl. Acad. Sci. U.S.A.1049052–9057. 10.1073/pnas.0700496104
116
StewartE. J. (2012). Growing unculturable bacteria.J. Bacteriol.1944151–4160. 10.1128/JB.00345-12
117
SwanB. K.Martinez-GarciaM.PrestonC. M.SczyrbaA.WoykeT.LamyD.et al (2011). Potential for chemolithoautotrophy among ubiquitous bacteria lineages in the dark ocean.Science3331296–1300. 10.1126/science.1203690
118
TangF.BarbacioruC.NordmanE.LiB.XuN.BashkirovV. I.et al (2010). RNA-Seq analysis to capture the transcriptome landscape of a single cell.Nat. Protoc.5516–535. 10.1038/nprot.2009.236
119
TangF.BarbacioruC.WangY.NordmanE.LeeC.XuN.et al (2009). mRNA-Seq whole-transcriptome analysis of a single cell.Nat. Methods6377–382. 10.1038/nmeth.1315
120
TaniguchiY.ChoiP. J.LiG. W.ChenH.BabuM.HearnJ.et al (2010). Quantifying E. coli proteome and transcriptome with single-molecule sensitivity in single cells.Science329533–538. 10.1126/science.1188308
121
TeleniusH.CarterN. P.BebbC. E.NordenskjoldM.PonderB. A.TunnacliffeA. (1992). Degenerate oligonucleotide-primed PCR: general amplification of target DNA by a single degenerate primer.Genomics13718–725. 10.1016/0888-7543(92)90147-K
122
ThompsonA. W.TurkarslanS.ArensC. E.Lopez Garcia de LomanaA.RamanA. V.StahlD. A.et al (2017). Robustness of a model microbial community emerges from population structure among single cells of a clonal population.Environ. Microbiol.193059–3069. 10.1111/1462-2920.13764
123
TringeS. G.von MeringC.KobayashiA.SalamovA. A.ChenK.ChangH. W.et al (2005). Comparative metagenomics of microbial communities.Science308554–557. 10.1126/science.1107851
124
TrouttA. B.McHeyzer-WilliamsM. G.PulendranB.NossalG. J. (1992). Ligation-anchored PCR: a simple amplification technique with single-sided specificity.Proc. Natl. Acad. Sci. U.S.A.899823–9825. 10.1073/pnas.89.20.9823
125
TurkarslanS.RamanA. V.ThompsonA. W.ArensC. E.GillespieM. A.von NetzerF.et al (2017). Mechanism for microbial population collapse in a fluctuating resource environment.Mol. Syst. Biol.13919. 10.15252/msb.20167058
126
van GestelJ.VlamakisH.KolterR. (2015). Division of labor in biofilms: the ecology of cell differentiation.Microbiol. Spectr.3:MB-0002-2014. 10.1128/microbiolspec.MB-0002-2014
127
VannierT.LeconteJ.SeeleuthnerY.MondyS.PelletierE.AuryJ. M.et al (2016). Survey of the green picoalga Bathycoccus genomes in the global ocean.Sci. Rep.6:37900. 10.1038/srep37900
128
VeeningJ. W.SmitsW. K.KuipersO. P. (2008). Bistability, epigenetics, and bet-hedging in bacteria.Annu. Rev. Microbiol.62193–210. 10.1146/annurev.micro.62.081307.163002
129
VenterJ. C.RemingtonK.HeidelbergJ. F.HalpernA. L.RuschD.EisenJ. A.et al (2004). Environmental genome shotgun sequencing of the Sargasso Sea.Science30466–74. 10.1126/science.1093857
130
VineyM.ReeceS. E. (2013). Adaptive noise.Proc. Biol. Sci.280:20131104. 10.1098/rspb.2013.1104
131
WangJ.ChenL.ChenZ.ZhangW. (2015). RNA-seq based transcriptomic analysis of single bacterial cells.Integr. Biol. (Camb)71466–1476. 10.1039/c5ib00191a
132
WenN.ZhaoZ.FanB.ChenD.MenD.WangJ.et al (2016). Development of droplet microfluidics enabling high-throughput single-cell analysis.Molecules21:881. 10.3390/molecules21070881
133
WolffA.Perch-NielsenI. R.LarsenU. D.FriisP.GoranovicG.PoulsenC. R.et al (2003). Integrating advanced functionality in a microfabricated high-throughput fluorescent-activated cell sorter.Lab Chip322–27. 10.1039/b209333b
134
WoykeT.SczyrbaA.LeeJ.RinkeC.TigheD.ClingenpeelS.et al (2011). Decontamination of MDA reagents for single cell whole genome amplification.PLOS ONE6:e26161. 10.1371/journal.pone.0026161
135
WoykeT.XieG.CopelandA.GonzalezJ. M.HanC.KissH.et al (2009). Assembling the marine metagenome, one cell at a time.PLOS ONE4:e5299. 10.1371/journal.pone.0005299
136
WuD.HugenholtzP.MavromatisK.PukallR.DalinE.IvanovaN. N.et al (2009). A phylogeny-driven genomic encyclopaedia of Bacteria and Archaea.Nature4621056–1060. 10.1038/nature08656
137
XiH. D.ZhengH.GuoW.Ganan-CalvoA. M.AiY.TsaoC. W.et al (2017). Active droplet sorting in microfluidics: a review.Lab Chip17751–771. 10.1039/c6lc01435f
138
YoonH. S.PriceD. C.StepanauskasR.RajahV. D.SierackiM. E.WilsonW. H.et al (2011). Single-cell genomics reveals organismal interactions in uncultivated marine protists.Science332714–717. 10.1126/science.1203163
139
YuF. B.WillisL.ChauR. M.ZambonA.HorowitzM.BhayaD.et al (2017). Long-term microfluidic tracking of coccoid cyanobacterial cells reveals robust control of division timing.BMC Biol.15:11. 10.1186/s12915-016-0344-4
140
ZhangD. Y.BrandweinM.HsuihT.LiH. B. (2001). Ramification amplification: a novel isothermal DNA amplification method.Mol. Diagn.6141–150. 10.1054/modi.2001.25323
141
ZhangK.MartinyA. C.ReppasN. B.BarryK. W.MalekJ.ChisholmS. W.et al (2006). Sequencing genomes from single cells by polymerase cloning.Nat. Biotechnol.24680–686. 10.1038/nbt1214
142
ZhangQ.WangT.ZhouQ.ZhangP.GongY.GouH.et al (2017). Development of a facile droplet-based single-cell isolation platform for cultivation and genomic analysis in microorganisms.Sci. Rep.7:41192. 10.1038/srep41192
143
ZhangQ.ZhangP.GouH.MouC.HuangW. E.YangM.et al (2015). Towards high-throughput microfluidic Raman-activated cell sorting.Analyst1406163–6174. 10.1039/c5an01074h
144
ZiegenhainC.ViethB.ParekhS.ReiniusB.Guillaumet-AdkinsA.SmetsM.et al (2017). Comparative analysis of single-cell RNA sequencing methods.Mol. Cell65631–643.e4. 10.1016/j.molcel.2017.01.023
145
ZongC.LuS.ChapmanA. R.XieX. S. (2012). Genome-wide detection of single-nucleotide and copy-number variations of a single human cell.Science3381622–1626. 10.1126/science.1229164
Summary
Keywords
single-cell analysis, microbes, heterogeneity, genomics, transcriptomics, next-generation sequencing
Citation
Chen Z, Chen L and Zhang W (2017) Tools for Genomic and Transcriptomic Analysis of Microbes at Single-Cell Level. Front. Microbiol. 8:1831. doi: 10.3389/fmicb.2017.01831
Received
25 April 2017
Accepted
06 September 2017
Published
20 September 2017
Volume
8 - 2017
Edited by
Dimitris G. Hatzinikolaou, National and Kapodistrian University of Athens, Greece
Reviewed by
Anne-Kristin Kaster, Karlsruhe Institute of Technology, Germany; Gwenael Piganeau, FR3724 Observatoire Océanologique de Banyuls sur Mer (OOB), France
Updates
Copyright
© 2017 Chen, Chen and Zhang.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Lei Chen, lchen@tju.edu.cn
This article was submitted to Systems Microbiology, a section of the journal Frontiers in Microbiology
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.