|
|
|
|
Genome Res. 15:1388-1392, 2005 ©2005 by Cold Spring Harbor Laboratory Press; ISSN 1088-9051/05 $5.00 Letter Gene expression profiling in single cells from the pancreatic islets of Langerhans reveals lognormal distribution of mRNA levels1 Department of Experimental Medical Science, Lund University, 221 84 Lund, Sweden 2 Department of Chemistry and Biosciences, Molecular Biotechnology, Chalmers University of Technology and TATAA Biocenter, Lundberg Laboratory, 405 30 Göteborg, Sweden 3 The Oxford Centre for Diabetes, Endocrinology and Metabolism, The Churchill Hospital, Oxford, OX3 7LJ, United Kingdom
The transcriptional machinery in individual cells is controlled by a relatively small number of molecules, which may result in stochastic behavior in gene activity. Because of technical limitations in current collection and recording methods, most gene expression measurements are carried out on populations of cells and therefore reflect average mRNA levels. The variability of the transcript levels between different cells remains undefined, although it may have profound effects on cellular activities. Here we have measured gene expression levels of the five genes ActB, Ins1, Ins2, Abcc8, and Kcnj11 in individual cells from mouse pancreatic islets. Whereas Ins1 and Ins2 expression show a strong cellcell correlation, this is not the case for the other genes. We further found that the transcript levels of the different genes are lognormally distributed. Hence, the geometric mean of expression levels provides a better estimate of gene activity of the typical cell than does the arithmetic mean measured on a cell population.
A typical eukaryotic cell contains
Here we have studied the expression of five genes in individual mouse pancreatic islet cells and MIN6 insulinoma cells (Miyazaki et al. 1990
Out of a total of 169 mouse islet cells, 84 were incubated in the presence of 5 mM glucose (low) and 85 in 20 mM glucose (high). In Table 1, the arithmetic means of the expression level of the five genes indicate that whereas the number of insulin transcripts per cell is in the order of several thousand copies, ActB and Abcc8 transcripts are present in a few hundred copies. The number of transcripts of the KATP-channel subunit Kcnj11 is <30 copies per cell. Based on the presence of Ins1 or Ins2 transcripts, it was concluded that at least 123 cells (73%) were -cells. The fraction of -cells in islet preparations is known to show large variations, and the average is 70%80% (Barg et al. 2000
To visualize the gene expression profile in a population of cells, distribution plots were used. Figure 1 shows histograms of ActB expression levels in both logarithmic and linear scale. As confirmed by the Shapiro-Wilk normality test, the transcript distribution is lognormal at 95% significance level (P = 0.05) (Table 1). The transcript distribution of the other genes is also lognormal at the same significance level. Corresponding data for the MIN6 cells are in Supplemental Table 1.
Lognormal distributions are common in nature. They are occasionally mistaken for normal (Gaussian) distributions, although the difference is fundamental (Limpert et al. 2001
We ascertained that the observed lognormal distribution of expression levels reflects true biological variability and is not an artifact of the technology or the approach used (see Supplemental information; Methods; Supplemental Fig. 1). The finding that cellular transcript levels are lognormally distributed has implications on the interpretation of gene expression data in general. If mRNA expression levels among cells are lognormally rather than normally distributed, then the average expression measured on a cell population does not reflect the expression in the typical cell in the population. The average value is strongly biased by a small population of cells with very active transcription of the particular gene. Accordingly, it may not be valid to extrapolate results of gene expression measurements on cell populations to the single-cell level. We analyzed this aspect by measuring the distribution of the ratios between the expression levels at high and low glucose concentration for Ins1 and Ins2 (see Table 3). Glucose stimulation has been reported to increase Ins1 and Ins2 expression two- to fivefold (Nielsen et al. 1985
The expression levels of Ins1 and Ins2 are affected similarly by glucose, suggesting that the two genes are regulated by similar or perhaps even a common mechanism (Wicksteed et al. 2001
Figure 2 shows the logarithm of Ins1 and Ins2 expression at low and high glucose levels in histograms. The large shift toward higher expression levels at elevated glucose concentration reflects 100-fold increase in insulin expression triggered by the sugar. Asymmetry, or skewness, of the distribution constitutes evidence of deviation from lognormal behavior. The skewness values for the genes analyzed, quantified as (X µ)3/(N 1) 3, are presented in Table 1. A positive value suggests that the distribution is skewed toward higher values, and vice versa. Ins1 distributions exhibited high skewness of opposite signs at low and high glucose concentrations. This may suggest the existence of two subpopulations of cells: one active to secrete insulin and one dormant, exhibiting a bimodal distribution. The bimodal gene induction model (Ko 1992 -cells has enhanced transcription levels characterized by a high mean value, while another subpopulation has low transcriptional activity characterized by a much lower mean value. Elevated glucose level thus increases the probability that individual cells are activated. The currently prevailing model of the action of enhancers is consistent with the binary model for gene activation (Fiering et al. 2000
An alternative interpretation of the data shown in Figure 2 is that high glucose "locks" the insulin gene in a high-expression state, whereas its expression at low glucose concentrations is more random. It should be noted that a fairly large portion of the cells at low glucose concentration exhibit as high insulin gene expression as that observed for most cells during stimulation with 20 mM glucose. When the expression levels of two genes are correlated, a common regulatory mechanism is often assumed. It may be through a mechanism that actually affects the two genes the same way, such as a common transcription factor, resulting in correlation at the single cell level, but it can also be a general increase in transcriptional activity due to, for example, environmental factors. The latter would also give rise to correlation in expression between genes, but not necessarily at the level of the individual cell. In our system, Ins1 and Ins2 are correlated on the cell level, while Ins1/Ins2 and ActB are correlated only on the population level. Our technology offers means to distinguish between these two cases and is expected to become especially useful for studies of molecular mechanisms underlying complex biological processes as well as disease.
Preparation and culture of cells Animals used in this study were healthy female National Maritime Research Institute (NMRI) mice aged 34 mo that were obtained from a commercial breeder (Bomholtgaard, Ry, Denmark) and fed a normal diet ad libitum. Care and use of animals were approved by the ethical committee of Lund University. The mice were sacrificed by cervical dislocation, and pancreatic islets were isolated by collagenase P digestion (Roche) (Olofsson et al. 2002
Mouse insulinoma MIN6-cells (passage 30 and above) were cultured in DMEM medium (10 mM glucose, Invitrogen) supplemented with 10% fetal calf serum, 100 U/mL penicillin, and 10 µg/mL streptomycin (all from Invitrogen) to
Single cell isolation and cDNA synthesis
Quantitative real-time PCR
Some cells appeared to lack at least one transcript; 96% of all collected samples contained detectable transcripts from at least one gene, and 83% had detectable levels of ActB transcript. For the analysis of primary cells in this article, only the insulin-producing
We thank Kristina Borglid, Britt-Marie Nilsson, and Lena Thiman for technical assistance with cell preparations. Research was supported by the Crafoord Foundation, the Royal and Hvitfeldtska Foundation, Swedish Research Council (8647), the Swedish Diabetes Association, the NovoNordisk Foundation, the Swedish Strategic Foundation, the Goran Gustafsson Foundation for Research in the Natural Sciences and Medicine, and the Medical Faculty Lund University. We also thank Prof. Miyazaki for providing us with the MIN6 cell line. P.R. is a Royal Society-Wolfson Research Fellow.
Article and publication are at http://www.genome.org/cgi/doi/10.1101/gr.3820805.
4 Corresponding author. [Supplemental material is available online at www.genome.org.]
Barg, S., Galvanovskis, J., Gopel, S.O., Rorsman, P., and Eliasson, L. 2000. Tight coupling between electrical activity and exocytosis in mouse glucagon-secreting Bishop, A.E. and Polak, J.M. 1997. The anatomy, organisation and ultrastructure of the islets of Langerhans. In Textbook of diabetes (eds. J. Pickup and G. Williams), pp. 6.1-6.16. Blackwell Science, Oxford, UK. Blake, W.J., Kaern, M., Cantor, C.R., and Collins, J.J. 2003. Noise in eukaryotic gene expression. Nature 422: 633-637.[CrossRef][Medline] Bustin, S.A. 2000. Absolute quantification of mRNA using real-time reverse transcription polymerase chain reaction assays. J. Mol. Endocrinol. 25: 169-193.[Abstract] Castano, J.P., Kineman, R.D., and Frawley, L.S. 1996. Dynamic monitoring and quantification of gene expression in single, living cells: A molecular basis for secretory cell heterogeneity. Mol. Endocrinol. 10: 599-605.[Abstract] Chiang, M.K. and Melton, D.A. 2003. Single-cell transcript analysis of pancreas development. Dev. Cell 4: 383-393.[CrossRef][Medline] Eliasson, L., Renstrom, E., Ding, W.G., Proks, P., and Rorsman, P. 1997. Rapid ATP-dependent priming of secretory granules precedes Ca(2+)-induced exocytosis in mouse pancreatic B-cells. J. Physiol. 503: 399-412.[CrossRef][Medline]
Elowitz, M.B., Levine, A.J., Siggia, E.D., and Swain, P.S. 2002. Stochastic gene expression in a single cell. Science 297: 1183-1186. Fiering, S., Whitelaw, E., and Martin, D.I. 2000. To be or not to be active: The stochastic nature of enhancer action. Bioessays 22: 381-387.[CrossRef][Medline] Gehwolf, R., Griessner, M., Pertl, H., and Obermeyer, G. 2002. First patch, then catch: Measuring the activity and the mRNA transcripts of a proton pump in individual Lilium pollen protoplasts. FEBS Lett. 512: 152-156.[CrossRef][Medline]
Kamme, F., Salunga, R., Yu, J., Tran, D.T., Zhu, J., Luo, L., Bittner, A., Guo, H.Q., Miller, N., Wan, J., et al. 2003. Single-cell microarray analysis in hippocampus CA1: Demonstration and validation of cellular heterogeneity. J. Neurosci. 23: 3607-3615. Ko, M.S. 1992. Induction mechanism of a single gene molecule: Stochastic or deterministic? Bioessays 14: 341-346.[CrossRef][Medline] Koch, A.L. 1966. The logarithm in biology, 1: Mechanisms generating the log-normal distribution exactly. J. Theor. Biol. 12: 276-290.[CrossRef][Medline]
Kuznetsov, V.A., Knott, G.D., and Bonner, R.F. 2002. General statistics of stochastic process of gene expression in eukaryotic cells. Genetics 161: 1321-1332. Levsky, J.M. and Singer, R.H. 2003. Gene expression and the myth of the average cell. Trends Cell Biol. 13: 4-6.[CrossRef][Medline]
Levsky, J.M., Shenoy, S.M., Pezo, R.C., and Singer, R.H. 2002. Single-cell gene expression profiling. Science 297: 836-840. Limpert, E., Stahel, W.A., and Abbt, M. 2001. Log-normal distributions across the sciences: Keys and clues. Bioscience 51: 341-352.[CrossRef] Liss, B., Franz, O., Sewing, S., Bruns, R., Neuhoff, H., and Roeper, J. 2001. Tuning pacemaker frequency of individual dopaminergic neurons by Kv4.3L and KChip3.1 transcription. EMBO J. 20: 5715-5724.[CrossRef][Medline] Luo, L., Salunga, R.C., Guo, H., Bittner, A., Joy, K.C., Galindo, J.E., Xiao, H., Rogers, K.E., Wan, J.S., Jackson, M.R., et al. 1999. Gene expression profiles of laser-captured adjacent neuronal subtypes. Nat. Med. 5: 117-122.[CrossRef][Medline]
McAdams, H.H. and Arkin, A. 1997. Stochastic mechanisms in gene expression. Proc. Natl. Acad. Sci. 94: 814-819.
Miyazaki, J., Araki, K., Yamato, E., Ikegami, H., Asano, T., Shibasaki, Y., Oka, Y., and Yamamura, K. 1990. Establishment of a pancreatic
Nielsen, D.A., Welsh, M., Casadaban, M.J., and Steiner, D.F. 1985. Control of insulin gene expression in pancreatic Olofsson, C.S., Gopel, S.O., Barg, S., Galvanovskis, J., Ma, X., Salehi, A., Rorsman, P., and Eliasson, L. 2002. Fast insulin secretion reflects exocytosis of docked granules in mouse pancreatic B-cells. Pflugers Arch. 444: 43-51.[CrossRef][Medline] Ozbudak, E.M., Thattai, M., Kurtser, I., Grossman, A.D., and van Oudenaarden, A. 2002. Regulation of noise in the expression of a single gene. Nat. Genet. 31: 69-73.[CrossRef][Medline] Paldi, A. 2003. Stochastic gene expression during cell differentiation: Order from disorder? Cell. Mol. Life Sci. 60: 1775-1778.[CrossRef][Medline]
Peixoto, A., Monteiro, M., Rocha, B., and Veiga-Fernandes, H. 2004. Quantification of multiple gene expression in individual cells. Genome Res. 14: 1938-1947.
Raser, J.M. and O'Shea, E.K. 2004. Control of stochasticity in eukaryotic gene expression. Science 304: 1811-1814. Stahlberg, A., Zoric, N., Aman, P., and Kubista, M. 2005. Quantitative real-time PCR for cancer detection: The lymphoma case. Expert Rev. Mol. Diagn. 5: 221-230.[CrossRef][Medline] Sucher, N.J., Deitcher, D.L., Baro, D.J., Warrick, R.M., and Guenther, E. 2000. Genes and channels: Patch/voltage-clamp analysis and single-cell RT-PCR. Cell Tissue Res. 302: 295-307.[CrossRef][Medline] Tietjen, I., Rihel, J.M., Cao, Y., Koentges, G., Zakhary, L., and Dulac, C. 2003. Single-cell transcriptional analysis of neuronal progenitors. Neuron 24: 161-175. Vandesompele, J., De Preter, K., Pattyn, F., Poppe, B., Van Roy, N., De Paepe, A., and Speleman, F. 2002. Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes. Genome Biol. 3: research0034.
Walters, M.C., Fiering, S., Eidemiller, J., Magis, W., Groudine, M., and Martin, D.I. 1995. Enhancers increase the probability but not the level of gene expression. Proc. Natl. Acad. Sci. 92: 7125-7129.
Wicksteed, B., Herbert, T.P., Alarcon, C., Lingohr, M.K., Moss, L.G., and Rhodes, C.J. 2001. Cooperativity between the preproinsulin mRNA untranslated regions is necessary for glucose-stimulated translation. J. Biol. Chem. 276: 22553-22558.
Received February 11, 2005; accepted in revised format June 28, 2005. This article has been cited by other articles:
|
|||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||