Abstract
Single-cell transcriptomics has recently emerged as a powerful technology to explore gene expression heterogeneity among single cells. Here we identify two major sources of technical variability: sampling noise and global cell-to-cell variation in sequencing efficiency. We propose noise models to correct for this, which we validate using single-molecule FISH. We demonstrate that gene expression variability in mouse embryonic stem cells depends on the culture condition.
This is a preview of subscription content, access via your institution
Access options
Subscribe to this journal
Receive 12 print issues and online access
$259.00 per year
only $21.58 per issue
Buy this article
- Purchase on SpringerLink
- Instant access to full article PDF
Prices may be subject to local taxes which are calculated during checkout
Similar content being viewed by others
Accession codes
References
Munsky, B., Neuert, G. & van Oudenaarden, A. Science 336, 183–187 (2012).
Eldar, A. & Elowitz, M.B. Nature 467, 167–173 (2010).
Hashimshony, T., Wagner, F., Sher, N. & Yanai, I. Cell Rep. 2, 666–673 (2012).
Sasagawa, Y. et al. Genome Biol. 14, R31 (2013).
Tang, F. et al. Nat. Methods 6, 377–382 (2009).
Ramsköld, D. et al. Nat. Biotechnol. 30, 777–782 (2012).
Islam, S. et al. Genome Res. 21, 1160–1167 (2011).
Picelli, S. et al. Nat. Methods 10, 1096–1098 (2013).
Shapiro, E., Biezuner, T. & Linnarsson, S. Nat. Rev. Genet. 14, 618–630 (2013).
Kivioja, T. et al. Nat. Methods 9, 72–74 (2012).
Shiroguchi, K., Jia, T.Z., Sims, P.A. & Xie, X.S. Proc. Natl. Acad. Sci. USA 109, 1347–1352 (2012).
Hug, H. & Schuler, R. J. Theor. Biol. 221, 615–624 (2003).
Shalek, A.K. et al. Nature 498, 236–240 (2013).
Islam, S. et al. Nat. Methods 11, 163–166 (2014).
Jaitin, D.A. et al. Science 343, 776–779 (2014).
Brennecke, P. et al. Nat. Methods 10, 1093–1095 (2013).
Ying, Q.-L. et al. Nature 453, 519–523 (2008).
The External RNA Controls Consortium. Nat. Methods 2, 731–734 (2005).
Raj, A., Peskin, C.S., Tranchina, D., Vargas, D.Y. & Tyagi, S. PLoS Biol. 4, e309 (2006).
Raj, A., van den Bogaard, P., Rifkin, S.A., van Oudenaarden, A. & Tyagi, S. Nat. Methods 5, 877–879 (2008).
Li, H. & Durbin, R. Bioinformatics 26, 589–595 (2010).
Meyer, L.R. et al. Nucleic Acids Res. 41, D64–D69 (2013).
Anders, S. & Huber, W. Genome Biol. 11, R106 (2010).
Robinson, M.D., McCarthy, D.J. & Smyth, G.K. Bioinformatics 26, 139–140 (2010).
Byrd, R.H., Lu, P., Nocedal, J. & Zhu, C. SIAM J. Sci. Comput. 16, 1190–1208 (1995).
Acknowledgements
This work was supported by a European Research Council Advanced grant (ERC-AdG 294325-GeneNoiseControl) and a Nederlandse Organisatie voor Wetenschappelijk Onderzoek (NWO) Vici award.
Author information
Authors and Affiliations
Contributions
D.G., L.K. and A.v.O. conceived the methods. D.G. developed the noise models, performed all computations and wrote the manuscript. L.K. performed all experiments and corrected the manuscript. A.v.O. guided experiments, data analysis and writing of the manuscript, and corrected the manuscript.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing financial interests.
Supplementary information
Supplementary Text and Figures
Supplementary Figures 1–15, Supplementary Table 1 and Supplementary Notes 1–4. (PDF 11686 kb)
Supplementary Table 2
GO terms enriched among genes with increased expression variability in serum versus 2i culture condition. Enriched biological processes and enriched molecular functions are given as separate lists. Only significantly enriched GO-terms (P < 0.05) were included. The lists indicate the GO-term ID, the hypergeometric P-value, the odds ratio, the expected number of genes associated with each GO-term, the observed number of genes for each GO-term, the size of the GO-term (total number of genes associated) and a short description. For the inference of over-represented GO terms, the set of differentially variable genes was compared to the universe of all genes expressed in the two conditions. The GOstats package was used to compute GO enrichment in R. (XLSX 82 kb)
Supplementary Table 3
Probe set composition of smFISH probes used. Each column represents a probe set for the gene specified in the column header. All probes were labeled on the 3' end with TMR, Alexa594 or Cy5. (XLSX 56 kb)
Rights and permissions
About this article
Cite this article
Grün, D., Kester, L. & van Oudenaarden, A. Validation of noise models for single-cell transcriptomics. Nat Methods 11, 637–640 (2014). https://doi.org/10.1038/nmeth.2930
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1038/nmeth.2930