Single cell sequencing

Single cell sequencing examines the sequence information from individual cells with optimized next generation sequencing (NGS) technologies, providing a higher resolution of cellular differences and a better understanding of the function of an individual cell in the context of its microenvironment.[1]

Background

A typical human cell consists of about 6 billion base pairs of DNA and 600 million bases of mRNA. With such huge amount of sequence, it is expensive and time-consuming to sequence by traditional Sanger sequencing. By using deep sequencing of DNA and RNA from single cell, cellular functions can be investigated extensively.[1] Like typical NGS experiments, the protocols of a single cell sequencing generally contain the following steps: isolation of single cell, nucleic acids extraction and amplification, sequencing library preparation, sequencing and bioinformatic data analysis. It is more challenging to perform single cell sequencing in comparison with sequencing from cells in bulk. The minimal amount of starting materials from a single cell make degradation, sample loss and contamination exert pronounced effects on quality of sequencing data. In addition, due to the picogram level of the amount of nucleic acids used,[2] heavy amplification is often needed during sample preparation of single cell sequencing, resulting in the uneven coverage, noise and inaccurate quantification of sequencing data.

Recent technical improvements make single cell sequencing a promising tool for approaching a set of seemly inaccessible problems. For example, heterogeneous samples, rare cell types, cell lineage relationships, mosaicism of somatic tissues, analyses of microbes that cannot be cultured, and disease evolution can all be elucidated through single cell sequencing.[3] Single cell sequencing was selected as the method of the year 2013 by Nature Publishing Group.[4]

Single cell genome (DNA) sequencing

Single cell DNA genome sequencing involves isolating a single cell, performing whole-genome-amplification (WGA), constructing sequencing libraries and then sequencing the DNA using a next-generation sequencer (ex. Ion Torrent, Illumina). It can be used in metagenomics studies and when sequencing the first time from novel species. In addition, it can be united with high throughput cell sorting of microorganisms and cancer. One popular method used for single cell genome sequencing is multiple displacement amplification and this enables research into various areas such as microbial genetics, ecology and infectious diseases. Furthermore, data obtained from microorganisms might establish processes for culturing in the future.[5] Some of the tools that can be used for single cell genome sequencing include: SPAdes, IDBA-UD, Cortex and HyDA.[6]

Method

This figure illustrates steps involved in workflow of single cell genome sequencing. MDA stands for multiple displacement amplification.

Multiple displacement amplification (MDA) is widely used technique, enabling amplifying femtograms of DNA from bacterium to micrograms for the use of sequencing. Reagents required for MDA reactions include: random primers and DNA polymerase from bacteriophage phi29. In 30 degree isothermal reaction, DNA is amplified with included reagents. As the polymerases manufacture new strands, a strand displacement reaction takes place, synthesizing multiple copies from each template DNA. At the same time, the strands that were extended antecedently will be displaced. MDA products result in a length of about 12 kb and ranges up to around 100 kb, enabling its use in DNA sequencing.[5] Other method includes MALBAC.[7]

Limitations

In comparison to PCR, MDA has higher amplification bias, either overrepresenting or underrepresenting different regions of the template, resulting in loss of some sequences. Because of this, amplification bias from one MDA reaction will be represented to the next reaction as well. Two ways to improve genome coverage include: pooling single-cell MDA reactions from the same cell type and pooling before the reaction is performed. Several ways to identify cells of the same strain are: fluorescent in situ hybridization (FISH) and conformation characteristics.[5]

Single-nucleotide polymorphisms (SNPs), which are a big part of genetic variation in the human genome, and copy number variation (CNV), pose problems in single cell sequencing, as well as the limited amount of DNA extracted from a single cell. Due to scant amounts of DNA, accurate analysis of DNA poses problems even after amplification since coverage is low and susceptible to errors. With MDA, average genome coverage is less than 80% and SNPs that are not covered by sequencing reads will be opted out. In addition, MDA shows a high ratio of allele dropout, not detecting alleles from heterozygous samples. Various SNP algorithms are currently in use but none are specific to single cell sequencing. MDA with CNV also poses the problem of identifying false CNVs that conceal the real CNVs. To solve this, when patterns can be generated from false CNVs, algorithms can detect and eradicate this noise to produce true variants.[7]

Applications

Microbiomes are the major targets of single cell genomics due to its difficulty for culturing. Single cell genomics is the one way to identify microbiomes’ identities and its genomes. The first microorganism used for single cell sequencing was a bacterium. When the data will be assembled in the near future, several new functions of these organisms might be discovered and might provide pros and cons regarding human health.[8]

Cancer sequencing is also an emerging application of scDNAseq. Fresh or frozen tumors may be analyzed and categorized with respect to SCNAs, SNVs, and rearrangements quite well using whole genome DNAS approaches [9] Cancer scDNAseq is particularly useful for examining the depth of complexity and compound mutations present in amplified therapeutic targets such as receptor tyrosine kinase genes (EGFR, PDGFRA etc.) where conventional population-level approaches of the bulk tumor are not able to resolve the co-occurrence patterns of these mutations within single cells of the tumor. Such overlap may provide redundancy of pathway activation and tumor cell resistance.

Single cell DNA methylome sequencing

One method for single cell DNA methylation sequencing.[10]

Single cell DNA genome sequencing quantifies DNA methylation. This is similar to single cell genome sequencing, but with the addition of a bisulfite treatment before sequencing. Forms include whole genome bisulfite sequencing,[10][11] and reduced representation bisulfite sequencing [12][13]

Comparison of single cell methylation sequencing methods in terms of coverage as at 2015 on Mus musculus

Single-cell RNA sequencing (scRNA-seq)

Current methods for quantifying molecular states of cells, from microarray to standard RNA-seq analysis, mostly depend on estimating the mean value from millions of cells by averaging the signal of individual cells. Given the heterogeneity of cell population, measurement of the mean values of signals overlooks the internal interactions and differences within a cell population that may be crucial for maintaining normal tissue functions and facilitating disease progression. Thus the cell-averaging experiments provide only partial information of the molecular state of the system.[14][15]

Single-cell RNA sequencing (scRNA-seq) provides the expression profile of individual cells. Through genes clustering analyses, rare cell types within a cell population can be identified, thereby making characterization of the subpopulation structure of a heterogeneous cell population become available. While tumor heterogeneity can be attributed to accumulated mutations, even genetically identical cells, under the same environment, display high variability of gene and protein expression levels.[16] However, RNA with low copy number, which may exert important functions in the cells, is usually undetectable or regarded as noise in traditional cell-averaging methods. Single-cell RNA sequencing on a large number of single cells can identify such uncommon RNA and also reveal the copy-number distribution of the whole mRNA population in individual cells. Knowledge about the shape of distribution can be used to understand the mechanisms of transcription regulation.[15][17]

Experimental procedures

Single-cell RNA sequencing workflow

Despite the advances in sequencing technologies, it is still unattainable to sequence RNA directly from single cell. Thus, in the current scRNA-seq protocols, RNA still needs to be converted to cDNA for sequencing. Principally, the current scRNA-seq methods contain the following steps: isolation of single cell and RNA, reverse transcription (RT), amplification, library generation and sequencing.

The ideal scRNA-seq preserves and accurately quantifies the initial relative abundance of mRNA in a cell, covers the entire transcript lengths with equal representation at each position, and retains strand information.[17] Nevertheless, a variety of noise and bias may be introduced in various steps of scRNA-seq protocol. For example, the step of reverse transcription is critical as the efficiency of the RT reaction determines the percentage of a cell’s RNA population that is eventually analyzed by the sequencer. The processivity of reverse transcriptases and the priming strategies used will affect full-length cDNA production and the generation of libraries biased toward 3’ or 5' end of genes.

In the amplification step, either PCR or in vitro transcription (IVT) is currently used to amplify cDNA. One of the advantages of PCR-based methods is able to generate full-length cDNA. However, different PCR efficiency on particular sequences (for instance, GC content and snapback structure) will also be exponentially amplified, producing libraries with uneven coverage. On the other hand, while libraries generated by IVT can avoid PCR-induced sequence bias, specific sequences may be transcribed inefficiently, thus causing sequence drop-out or generating incomplete sequences.[1][14] Several scRNA-seq protocols have been published: Tang et al.,[18] STRT,[19] SMART-seq,[20] CEL-seq[21] and Quartz-seq.[22]

Applications

The number of circulating tumor cells (CTC) in peripheral blood of cancer patients has been shown to correlate to prognosis.[23] However, it is challenging to enumerate and characterize the isolated CTCs as they are often contaminated with a large number of leukocytes and erythrocytes. Single cell RNA-seq could be applied to differentiate cancer cells from normal blood cells and obtain the expression profiles of tumor cells at the same time. Similarly, single cell RNA-seq can also be used to analyze rare cell types in early human embryo and adult stem cells, both of which exist transiently and difficult to be characterized with current technologies. Finally, single cell analysis can be applied to the study of infectious diseases.[24]

Considerations

Isolation of single cells

There is currently no standardized technique for single-cell isolation. Individual cells can be collected by micromanipulation, for example by serial dilution or by using a patch pipette or nanotube to harvest a single cell.[25][26] The advantages of micromanipulation are ease and low cost, but they are laborious and susceptible to misidentification of cell types under microscope. Laser-capture microdissection (LCM) can also be used for collecting single cells. Although LCM preserves the knowledge of the spatial location of a sampled cell within a tissue, it is hard to capture a whole single cell without also collecting the materials from neighboring cells.[14][27][28] High-throughput methods for single cell isolation include fluorescence-activated cell sorting (FACS) and microfluidics. Both of FACS and microfluidics are accurate, automatic and capable of isolating unbiased samples. However, both methods require detaching cells from their microenvironments first, thereby causing perturbation to the transcriptional profiles in RNA expression analysis.[29][30]

Number of cells to be analyzed

scRNA-Seq

Generally speaking, for a typical bulk cell RNA-sequencing (RNA-seq) experiment, ten million reads are generated and a gene with higher than the threshold of 50 reads per kb per million reads (RPKM) is considered expressed. For a gene that is 1kb long, this corresponds to 500 reads and a minimum coefficient of variation (CV) of 4% under the assumption of the Poisson distribution. For a typical mammalian cell containing 200,000 mRNA, sequencing data from at least 50 single cells need to be pooled in order to achieve this minimum CV value. However, due to the efficiency of reverse transcription and other noise introduced in the experiments, more cells are required for accurate expression analyses and cell type identification.[14]

References

  1. 1 2 3 Eberwine J, Sul JY, Bartfai T, Kim J (January 2014). "The promise of single-cell sequencing". Nat. Methods. 11 (1): 25–27. doi:10.1038/nmeth.2769. PMID 24524134.
  2. ↑ Shintaku H, Nishikii H, Marshall LA, Kotera H, Santiago JG (2014). "On-Chip Separation and Analysis of RNA and DNA from Single Cells". Anal. Chem. 86 (4): 1953–7. doi:10.1021/ac4040218. PMID 24499009.
  3. ↑ Nawy, T. (2013). "Single-cell sequencing". Nat Methods. 11 (1): 18. doi:10.1038/nmeth.2771. PMID 24524131.
  4. ↑ "Method of the Year 2013". Nat Methods. 11 (1). 2014. doi:10.1038/nmeth.2801.
  5. 1 2 3 "Lasken, RS. (2007). "Single-cell Genomic Sequencing using Multiple Displacement Amplification". Curr Opin Microbiol. 10 (5): 510–16. doi:10.1016/j.mib.2007.08.005. PMID 17923430."
  6. ↑ Taghavi, Z.; Movahedi NS; Draghici S; Chitsaz H. (2013). "Distilled single-cell genome sequencing and de novo assembly for sparse microbial communities.". Bioinformatics. 29 (19): 2395–2401. doi:10.1093/bioinformatics/btt420.
  7. 1 2 "Ning, L.; Liu G; Li G; Hou Y; Tong Y; He J. (2014). "Current Challenges in the Bioinformatics of Single Cell Genomics.". Front Oncol. 4 (7). doi:10.3389/fonc.2014.00007. PMC 3902584Freely accessible. PMID 24478987."
  8. ↑ Blainey, PC.; Quake SR. (2014). "Dissecting genomic diversity, one cell at a time". Nat Methods. 11 (1): 19–21. doi:10.1038/nmeth.2783. PMC 3947563Freely accessible. PMID 24524132.
  9. ↑ (Francis J, Zheng CZ, Maire C et al. Cancer Discovery 2014).
  10. 1 2 Farlik, M; Sheffield, NC; Nuzzo, A; Datlinger, P; Schönegger, A; Klughammer, J; Bock, C (3 March 2015). "Single-cell DNA methylome sequencing and bioinformatic inference of epigenomic cell-state dynamics.". Cell reports. 10 (8): 1386–97. doi:10.1016/j.celrep.2015.02.001. PMC 4542311Freely accessible. PMID 25732828.
  11. ↑ Smallwood, SA; Lee, HJ; Angermueller, C; Krueger, F; Saadeh, H; Peat, J; Andrews, SR; Stegle, O; Reik, W; Kelsey, G (August 2014). "Single-cell genome-wide bisulfite sequencing for assessing epigenetic heterogeneity.". Nature Methods. 11 (8): 817–20. doi:10.1038/nmeth.3035. PMC 4117646Freely accessible. PMID 25042786.
  12. ↑ Guo, H; Zhu, P; Wu, X; Li, X; Wen, L; Tang, F (December 2013). "Single-cell methylome landscapes of mouse embryonic stem cells and early embryos analyzed using reduced representation bisulfite sequencing.". Genome Research. 23 (12): 2126–35. doi:10.1101/gr.161679.113. PMC 3847781Freely accessible. PMID 24179143.
  13. ↑ Guo, H; Zhu, P; Guo, F; Li, X; Wu, X; Fan, X; Wen, L; Tang, F (May 2015). "Profiling DNA methylome landscapes of mammalian cells with single-cell reduced-representation bisulfite sequencing.". Nature protocols. 10 (5): 645–59. doi:10.1038/nprot.2015.039. PMID 25837417.
  14. 1 2 3 4 "Shapiro, E.; Biezuner T; Linnarsson S. (2013). "Single-cell sequencing-based technologies will revolutionize whole-organism science". Nat Rev Genet. 14 (9): 618–30. doi:10.1038/nrg3542. PMID 23897237."
  15. 1 2 Kolodziejczyk, Aleksandra A.; Kim, Jong Kyoung; Svensson, Valentine; Marioni, John C.; Teichmann, Sarah A. (May 2015). "The Technology and Biology of Single-Cell RNA Sequencing". Molecular Cell. 58 (4): 610–620. doi:10.1016/j.molcel.2015.04.005.
  16. ↑ Munsky B, Neuert G, van Oudenaarden A (2012). "Using gene expression noise to understand gene regulation". Science. 336 (6078): 183–7. doi:10.1126/science.1216379. PMC 3358231Freely accessible. PMID 22499939.
  17. 1 2 "Hebenstreit, D. (2012). "Methods, Challenges and Potentials of Single Cell RNA-seq.". Biology. 1 (3): 658–667. doi:10.3390/biology1030658."
  18. ↑ Tang, F.; Barbacioru C; Wang Y; et al. (2009). "mRNA-Seq whole-transcriptome analysis of a single cell". Nat Methods. 6 (5): 377–82. doi:10.1038/NMETH.1315. PMID 19349980.
  19. ↑ Islam, S.; Kjällquist U; Moliner A; et al. (2011). "Characterization of the single-cell transcriptional landscape by highly multiplex RNA-seq.". Genome Res. 21 (7): 1160–7. doi:10.1101/gr.110882.110. PMC 3129258Freely accessible. PMID 21543516.
  20. ↑ Ramsköld, D.; Luo S; Wang Y-C; et al. (2012). "Full-length mRNA-Seq from single-cell levels of RNA and individual circulating tumor cells.". Nat. Biotechnol. 30 (8): 777–82. doi:10.1038/nbt.2282. PMC 3467340Freely accessible. PMID 22820318.
  21. ↑ ,Hashimshony, T.; Wagner F; Sher N; Yanai I. (2012). "CEL-Seq: single-cell RNA-Seq by multiplexed linear amplification.". Cell Rep. 2 (3): 666–73. doi:10.1016/j.celrep.2012.08.003. PMID 22939981.
  22. ↑ Sasagawa, Y.; Nikaido I; Hayashi T; et al. (2013). "Quartz-Seq: a highly reproducible and sensitive single-cell RNA sequencing method, reveals non-genetic gene-expression heterogeneity.". Genome Biol. 14 (4): R31. doi:10.1186/gb-2013-14-4-r31. PMC 4054835Freely accessible. PMID 23594475.
  23. ↑ Olmos, D.; Arkenau H-T; Ang JE; et al. (2009). "Circulating tumour cell (CTC) counts as intermediate end points in castration-resistant prostate cancer (CRPC): a single-centre experience.". Ann. Oncol. 20 (1): 27–33. doi:10.1093/annonc/mdn544. PMID 18695026.
  24. ↑ Avraham, R; Haseley, N; Brown, D; Penaranda, C; Jijon, HB; Trombetta, JJ; Satija, R; Shalek, AK; Xavier, RJ; Regev, A; Hung, DT (10 September 2015). "Pathogen Cell-to-Cell Variability Drives Heterogeneity in Host Immune Responses.". Cell. 162 (6): 1309–21. doi:10.1016/j.cell.2015.08.027. PMID 26343579.
  25. ↑ Zong, C; Lu S; Chapman AR; Xie XS. (2012). "Genome-wide detection of single-nucleotide and copy-number variations of a single human cell". Science. 338 (6114): 1622–6. doi:10.1126/science.1229164. PMC 3600412Freely accessible. PMID 23258894.
  26. ↑ Kurimoto, K.; Yabuta Y; Ohinata Y; Saitou M. (2007). "Global single-cell cDNA amplification to provide a template for representative high-density oligonucleotide microarray analysis". Nat Protoc. 2 (3): 739–52. doi:10.1038/nprot.2007.79. PMID 17406636.
  27. ↑ Yachida, S.; Jones S; Bozic I; et al. (2010). "Distant metastasis occurs late during the genetic evolution of pancreatic cancer". Nature. 467 (7319): 1114–7. doi:10.1038/nature09515. PMC 3148940Freely accessible. PMID 20981102.
  28. ↑ Frumkin, D.; Wasserstrom A; Itzkovitz S; Harmelin A; Rechavi G; Shapiro E. (2008). "Amplification of multiple genomic loci from single cells isolated by laser micro-dissection of tissues.". BMC Biotechnol. 8 (17). doi:10.1186/1472-6750-8-17. PMC 2266725Freely accessible. PMID 18284708.
  29. ↑ Dalerba, P.; Kalisky T; Sahoo D; et al. (2011). "Single-cell dissection of transcriptional heterogeneity in human colon tumors.". Nat. Biotechnol. 29 (12): 1120–7. doi:10.1038/nbt.2038. PMC 3237928Freely accessible. PMID 22081019.
  30. ↑ White, AK.; VanInsberghe M; Petriv OI; et al. (2011). "High-throughput microfluidic single-cell RT-qPCR". Proc Natl Acad Sci U S A. 108 (34): 13999–4004. doi:10.1073/pnas.1019446108. PMC 3161570Freely accessible. PMID 21808033.
This article is issued from Wikipedia - version of the 9/26/2016. The text is available under the Creative Commons Attribution/Share Alike but additional terms may apply for the media files.