Identifying single-cell molecular programs by stochastic profiling

Cells in tissues can be morphologically indistinguishable yet show molecular expression patterns that are remarkably heterogeneous. Here we describe an approach to comprehensively identify co-regulated, heterogeneously expressed genes among cells that otherwise appear identical. The technique, called stochastic profiling, involves repeated, random selection of very small cell populations via laser-capture microdissection followed by a customized single-cell amplification procedure and transcriptional profiling. Fluctuations in the resulting gene-expression measurements are then analyzed statistically to identify transcripts that are heterogeneously coexpressed. We stochastically profiled matrix-attached human epithelial cells in a three-dimensional culture model of mammary-acinar morphogenesis. Of 4,557 transcripts, we identified 547 genes with strong cell-to-cell expression differences. Clustering of this heterogeneous subset revealed several molecular 'programs' implicated in protein biosynthesis, oxidative-stress responses and NF-κB signaling, which we independently confirmed by RNA fluorescence in situ hybridization. Thus, stochastic profiling can reveal single-cell heterogeneities without the need to measure expression in individual cells.

[1]  T. Hunter,et al.  PI-3 kinase activity is necessary for ERK1/2-induced disruption of mammary epithelial architecture , 2009, Breast Cancer Research.

[2]  E. Petricoin,et al.  Laser Capture Microdissection , 1996, Science.

[3]  G Brady,et al.  Construction of cDNA libraries from single cells. , 1993, Methods in enzymology.

[4]  David J. Sheskin,et al.  Handbook of Parametric and Nonparametric Statistical Procedures , 1997 .

[5]  A. E. Hirsh,et al.  Noise Minimization in Eukaryotic Gene Expression , 2004, PLoS biology.

[6]  J. Sedat,et al.  Fluorescence microscopy: reduced photobleaching of rhodamine and fluorescein protein conjugates by n-propyl gallate. , 1982, Science.

[7]  Tomoharu Kajiyama,et al.  Quantitative analysis of gene expression in a single cell by qPCR , 2009, Nature Methods.

[8]  Marc S Halfon,et al.  A combinatorial code for pattern formation in Drosophila oogenesis. , 2008, Developmental cell.

[9]  R. Sakai,et al.  A Novel RNA-Binding Protein, Ossa/C9orf10, Regulates Activity of Src Kinases To Protect Cells from Oxidative Stress-Induced Apoptosis , 2008, Molecular and Cellular Biology.

[10]  Esteban O. Mazzoni,et al.  Stochastic spineless expression creates the retinal mosaic for colour vision , 2006, Nature.

[11]  M. Karin,et al.  Phosphorylation meets ubiquitination: the control of NF-[kappa]B activity. , 2000, Annual review of immunology.

[12]  Liliána Z Fehér,et al.  Real-time polymerase chain reaction-based exponential sample amplification for microarray gene expression profiling. , 2005, Analytical biochemistry.

[13]  Mario Cortina-Borja,et al.  Handbook of Parametric and Nonparametric Statistical Procedures, 5th edn , 2012 .

[14]  David M. Stern,et al.  Stress-Associated Endoplasmic Reticulum Protein 1 (Serp1)/Ribosome-Associated Membrane Protein 4 (Ramp4) Stabilizes Membrane Proteins during Stress and Facilitates Subsequent Glycosylation , 1999, The Journal of cell biology.

[15]  W. Staines,et al.  Reduction of Lipofuscin-like Autofluorescence in Fluorescently Labeled Tissue , 1999, The journal of histochemistry and cytochemistry : official journal of the Histochemistry Society.

[16]  Ian O Ellis,et al.  Basal-like breast cancer: a critical review. , 2008, Journal of clinical oncology : official journal of the American Society of Clinical Oncology.

[17]  P. Lásló,et al.  Multilineage Transcriptional Priming and Determination of Alternate Hematopoietic Cell Fates , 2006, Cell.

[18]  J. Derisi,et al.  Single-cell proteomic analysis of S. cerevisiae reveals the architecture of biological noise , 2006, Nature.

[19]  J. Debnath,et al.  Akt activation disrupts mammary acinar architecture and enhances proliferation in an mTOR-dependent manner , 2003, The Journal of cell biology.

[20]  A. T. McKay,et al.  Distribution of the Coefficient of Variation and the Extended “T” Distribution , 1932 .

[21]  Jan Vijg,et al.  Increased cell-to-cell variation in gene expression in ageing mouse heart , 2006, Nature.

[22]  Ronald N Germain,et al.  Modeling T Cell Antigen Discrimination Based on Feedback Control of Digital ERK Responses , 2005, PLoS biology.

[23]  Steven A. Julious,et al.  Using confidence intervals around individual means to assess statistical significance between two means , 2004 .

[24]  C. Klein,et al.  Gene expression profiling of single cells on large-scale oligonucleotide arrays , 2006, Nucleic acids research.

[25]  Daniel Birnbaum,et al.  ALDH1 is a marker of normal and malignant human mammary stem cells and a predictor of poor clinical outcome. , 2007, Cell stem cell.

[26]  Kazuki Kurimoto,et al.  An improved single-cell cDNA amplification method for efficient high-density oligonucleotide microarray analysis , 2006, Nucleic acids research.

[27]  Hannah H. Chang,et al.  Transcriptome-wide noise controls lineage choice in mammalian progenitor cells , 2008, Nature.

[28]  E. Cox,et al.  Real-Time Kinetics of Gene Activity in Individual Bacteria , 2005, Cell.

[29]  David Bryder,et al.  Transcription factor profiling in individual hematopoietic progenitors by digital RT-PCR , 2006, Proceedings of the National Academy of Sciences.

[30]  M. Birnbaum,et al.  The role of FoxO in the regulation of metabolism , 2008, Oncogene.

[31]  Tony Hunter,et al.  Real-time imaging reveals that noninvasive mammary epithelial acini can contain motile cells , 2007, The Journal of cell biology.

[32]  J. Avery Critical review. , 2006, The Journal of the Arkansas Medical Society.

[33]  Wen-Lin Kuo,et al.  A collection of breast cancer cell lines for the study of functionally distinct cancer subtypes. , 2006, Cancer cell.

[34]  P. Rorsman,et al.  Gene expression profiling in single cells from the pancreatic islets of Langerhans reveals lognormal distribution of mRNA levels. , 2005, Genome research.

[35]  Douglas A. Lauffenburger,et al.  Common effector processing mediates cell-specific responses to stimuli , 2007, Nature.

[36]  J. Rihel,et al.  Single-Cell Transcriptional Analysis of Neuronal Progenitors , 2003, Neuron.

[37]  G. Nolan,et al.  Mapping normal and cancer cell signalling networks: towards single-cell proteomics , 2006, Nature Reviews Cancer.

[38]  Tobias Schmelzle,et al.  Functional role and oncogene-regulated expression of the BH3-only factor Bmf in mammary epithelial anoikis and morphogenesis , 2007, Proceedings of the National Academy of Sciences.

[39]  J E Ferrell,et al.  The biochemical basis of an all-or-none cell fate switch in Xenopus oocytes. , 1998, Science.

[40]  Jayanta Debnath,et al.  Morphogenesis and oncogenesis of MCF-10A mammary epithelial acini grown in three-dimensional basement membrane cultures. , 2003, Methods.

[41]  Jayanta Debnath,et al.  Modelling glandular epithelial cancers in three-dimensional cultures , 2005, Nature Reviews Cancer.