| Nucleic Acids Research | Pages |
Expression profiling of single cells using 3 prime end amplification (TPEA) PCR
Introduction
Materials And Methods
Lymphoblastoma cell sorting
Reverse transcription and cDNA amplification
Gene-specific PCR
Extraction of neuronal contents
Results
TPEA-PCR
Reproducibility and specificity of 3[prime] amplification of cDNA from single cells
Analysis of complex cell systems: expression profiling of neurons
Discussion
Acknowledgements
References
Expression profiling of single cells using 3 prime end amplification (TPEA) PCR
ABSTRACT
INTRODUCTION
DNA sequence generated by genome and expressed sequence tag (EST) sequencing projects will provide the basis to further our understanding of the control and mode of action of individual gene products. In this respect, analysis and comparison of when, where and to what degree genes are expressed, commonly known as expression profiling, is certain to play an essential role in the functional characterization of newly identified genes (1,2).
Many tissues, such as the immune and nervous systems, are composed of highly heterogeneous cell populations. A key factor in understanding their physiology, and the role of specific gene products expressed within them, will be the ability to examine gene usage in the context of this cellular diversity. In the past, methods such as northern blotting and nuclease protection assays were employed to study gene expression. More recently, a number of powerful new technologies have been developed for assessing simultaneously the expression of large numbers of genes (3,4). All these techniques, however, require relatively large amounts of RNA and currently lack the sensitivity to analyse specimens derived from small populations of cells or indeed from an individual cell.
At present, methodologies for the analysis of gene expression within single cells are limiting. Whilst in situ hybridization provides detailed information on a gene's cellular expression pattern in intact tissue, be it whole-mounts or tissue sections, the technique is relatively laborious to perform and is incapable of analysing more than a very small number of transcripts in a single preparation. In recent years, the polymerase chain reaction (PCR) has been used successfully to investigate gene expression in cytoplasmic samples derived from single cells (5). The nested-primer approach has good sensitivity, but the analysis is restricted to a small number of closely related genes from specific gene families (6). Other techniques capable of detecting the expression of unrelated genes in a single cell include T7 RNA polymerase amplification of mRNA (7,8), and PCR after prior homopolymeric tailing of the first strand cDNA (9-11). However, neither of these approaches have been demonstrated to be capable of analysing more than a small number of genes and are also not widely used. The former is technically difficult, whilst the latter may be biased against long transcripts and often requires subsequent cloning of the amplified products. There exists therefore a great need for improved methodologies that can quickly, easily and sensitively analyse gene expression at the cellular level.
Here we report the development of an approach, called 3 prime end amplification (TPEA) PCR, for the analysis of gene expression in small samples or single cells, which overcomes many of the technical limitations of existing methodologies. In this method, the 3[prime] region of mRNA is amplified arbitrarily by PCR using a novel combination of primers. The amplified cDNA, which represents the most diverse region of gene sequence, can then be analysed by a second round of PCR using gene-specific primers. We demonstrate that it is possible to analyse the expression of up to 40 genes (20 in duplicate) in single human lymphoblastoma cells. We also show that this method is ideally suited to the analysis of genes expressed at low levels in small populations of cells, by demonstrating the expression of the adenosine A2a receptor in cholinergic neurons of the rat striatum.
MATERIALS AND METHODS
Lymphoblastoma cell sorting
An Epstein Barr virus transformed lymphoblastoid cell line (HRC575, ECACC, Porton Down, UK) was maintained in log phase growth in RPMI 1640 medium supplemented with 16% foetal calf serum, 2 mM l-glutamine and penicillin-streptomycin (100 U/ml and 100 mg/ml, respectively). Cells at ~106/ml were stained with the bisbenzimadazole dye Hoecsht 33342 (Sigma, Poole, UK) at 1 µg/ml for 30 min at 37°C. Cells were sorted by using the Autoclone attachment of a Coulter Elite ESP flow cytometer, 300 mW of all lines UV from a Coherent 306 laser and by using single drop and complete abort sorting settings. Time of flight, forward and right angle scatter and Hoechst fluorescence peak and area measurements were used to ensure the sorting of single cells. The accuracy of sorting (both spatial and numerical) was tested by sorting single fluorescent beads (DNA Check, Coulter Corp.) into 96-well plates and viewing the plates on a fluorescence microscope.
Reverse transcription and cDNA amplification
Lymphoblastoma cells were sorted into 96-well plates containing 7 µl of freshly prepared lysis buffer [50 mM Tris-HCl (pH 8.3), 75 mM KCl, 3 mM MgCl2, 5 mM NP-40 (Sigma) and 1.5 U of RNase inhibitor (Pharmacia, Milton Keynes, UK)]. This buffer leaves the nucleus intact (11). After 5 min on ice, nuclei were removed by centrifugation (8000 g, 5 min at 4°C), the supernatant aspirated and the RNA reverse transcribed in a reaction volume of 10 µl containing; 1× first-strand buffer, 200 U M-MLV reverse transcriptase (Gibco BRL, Paisley, UK), 0.5 ng reverse transcription primer for 60 min at 37°C. The RT primer was composed of an anchored oligo(dT) primer with a specific 5[prime] heel sequence: CTC TCA AGG ATC TTA CCG CTT TTT TTT TTT TTT TTT (A,G,C).
Second-strand cDNA synthesis was initiated by incubation of the first-strand cDNA with 1 ng of a primer consisting of (5[prime] to 3[prime]) a 20 base sequence selected due to its absence from the mammalian databases, a stretch of five random nucleotides and a defined pentameric sequence (CTG CAT CTA TCT AAT GCT CCNN NNN CGAGA where N represents C, G, T or A) for 15 min at 50°C under amplification conditions described below. Although this primer will undoubtedly prime second-strand DNA synthesis at many sites on the first strand cDNA, the subsequent PCR between the heel sequence of the oligo(dT) primer and the arbitrary primer closest to the 5[prime] end, ensures amplification of cDNA sequences complementary to the 3[prime] ends of the polyA tail. After allowing the second-strand primer to anneal, primer extension was performed at 72°C for 10 min using AmpliTaq DNA polymerase (0.35 U, Applied Biosystems, Warrington, UK) in PCR buffer containing 67 mM Tris-HCl (pH 8.3), 4.5 mM MgCl2 and 0.5 mM dNTPs. Subsequently, 0.4 ng of 3[prime] heel primer (CTC TCA AGG ATC TTA CCG C) was added and the reaction subjected to 10 cycles of 92°C for 2.5 min, 60°C for 1.5 min and 72°C extension for 1 min, followed by a final 10 min extension. A further 125 ng of second-strand primer and 50 ng of 3[prime] heel primer were then added in 10 µl of PCR reaction mix. After 15 cycles (as before), a further 10 µl of PCR reaction mix containing 125 ng second-strand primer and 50 ng of 3[prime] heel primer were added to the reaction and subjected to another 15 rounds of PCR. The final reaction mixture was then diluted to 200 µl with 10 mM Tris-0.1 mM EDTA (pH 8.1). Five microlitre samples were used for subsequent gene-specific PCR assays.
Gene-specific PCR
Samples (5 µl) of amplified cDNA were subjected to hot-lid PCR carried out in 1× PCR buffer (3.5 mM MgCl2, pH 8.8) containing 12.5% sucrose, 0.1 mM cresol red, 12 mM [beta]-mercaptoethanol, 0.5 mM dNTPs (Pharmacia), 0.6 U AmpliTaq DNA polymerase (Applied Biosystems), and primers at 100 ng/reaction. Amplifications were carried out on PTC-225 thermal cyclers (Tetrad, MJ Research, MA, US). Following an initial 2 min denaturing step (92°C), each PCR cycle consisted of 30 s denaturing (92°C), 90 s annealing (55°C) and 60 s elongation (72°C). After the final cycle the reaction was held for 10 min at 72°C. The PCR products were then separated on a 2.5% agarose gel, stained with ethidium bromide and photographed. All gene-specific primers are listed in Table 1.
![]() |
Table 1. Primer sequences used for gene-specific assays
Extraction of neuronal contents
Coronal slices (300 µm) from 14-28 day-old male Sprague-Dawley rats containing the striatum were viewed with a Zeiss Axioskop microscope (Carl Zeiss Ltd, Welwyn Garden City, UK) fitted with a ×64 water-immersion objective lens together with gradient contrast optics (Luigs and Neumann, Ratingen, Germany). Light in the infrared range (>740 nm) was used in conjunction with a contrast-enhancing Newvicon camera (Hamamatsu, Hamamatsu City, Japan) to resolve individual neurones within slices (27). The cytoplasm from large cells (>30 µm in one dimension) was gently aspirated under visual control into a patch-clamp recording electrode until at least 40% of the somatic cytoplasm had been collected. The electrode was then withdrawn from the cell to form an outside-out patch which prevented contamination when the electrode was subsequently withdrawn from the slice. The contents of the electrode were forced into a microtube and reverse transcribed, subjected to 3[prime] cDNA amplification, and 2.5% of the product used in each gene specific PCR reaction.
RESULTS
TPEA-PCR
Assay development was performed on lymphoblastoma cells in the G0/G1 phase of the cell cycle. Groups of 100, 10 and single cells were flow sorted into wells containing lysis buffer and the mRNA reverse transcribed. A proportion of the sorted cells then underwent 3[prime] end amplification, as described in Materials and Methods (Fig.
Reproducibility and specificity of 3[prime] amplification of cDNA from single cells
The reproducibility of this technique was assessed when the expression of 20 genes was examined in duplicate on single lymphoblastoma cells. Duplicate PCR assays revealed almost identical expression profiles both between assays on a single cell and between different cells (Fig.
Analysis of complex cell systems: expression profiling of neurons
Having shown that TPEA-PCR permits expression profiling of single cells, we wished to apply it to the analysis of a complex cell system in vivo. Striatal cholinergic interneurons are readily identifiable in rat brain slices due to their large size (>30 µm diameter) when compared to surrounding cell types, which are predominately medium spiny neurons (<15 µm diameter). After electrophysical characterisation of the cells (14), cytoplasmic samples were taken using a patch pipette, reverse transcribed and the cDNA amplified. The expression of a variety of genes of interest was then investigated (two representative expression profiles of cholinergic interneurons are shown in Fig.
Figure 1. Outline of TPEA-PCR. Details of the protocol are given in Materials and Methods, as are the sequences of the anchored oligo(dT) primer and the partially degenerate second strand primer.
DISCUSSION
Current methodologies for analysing gene expression in small samples of RNA, particularly from single cells, are severely limited in terms of the number and diversity of genes that it is possible to analyse, and the difficulty of experimental procedures. The methodology outlined here overcomes many of these limitations.
Sequence diversity between genes is at its greatest in the 3[prime] untranslated region (16) and this region provides the most unique target for gene-specific assays; this is especially important when wishing to differentiate between closely related members of a gene family. Our aim in designing this procedure has therefore been to amplify preferentially this portion of the mRNA sequences. cDNA synthesis by reverse transcriptase is initiated by an anchored oligo-dT priming so that the 3[prime] region of all genes is represented in the resulting single-stranded cDNA. A 5[prime]-specific heel is incorporated into this primer for use in the subsequent amplification procedure. For second-strand synthesis, it was seen as desirable that ~1 kb from the 3[prime] end of each gene could be selected and amplified by PCR. Assuming a completely random length of nucleotide sequence, it would be expected that a given five base sequence would appear every 1024 bases (45), even though some nucleotide sequences are more common than others (17). A pentameric sequence was selected that would enable the primer to initiate second-strand synthesis in an arbitrary manner within 1 kb of the 3[prime] end of the mRNA. A search of 30 gene sequences revealed that at least one copy of this five base sequence was present in this region of each gene. 5[prime] to this, five bases of random sequence (N5) were incorporated in order to stabilize the interaction of the arbitrary pentameric sequence, which in turn was flanked by a specific heel sequence. A similar primer design strategy has been demonstrated to be useful for DOP-PCR primers (18). After a single round of second-strand synthesis, each cDNA contains a specific priming site 5[prime] and 3[prime] to the region of interest, thus allowing amplification of the intervening sequence. At present, the exact representation of each mRNA species in the amplified cDNA relative to the starting material is unclear. As far as we have been able to ascertain, however, the majority of the mRNA species represented in the first-strand cDNA pool before amplification, as detected by conventional RT-PCR, were also detectable after amplification.
Figure 2. Demonstration of cDNA amplification. cDNA was prepared from varying numbers of sorted lymphoblastoma cells, and analysed by RT-PCR for the expression of eight `housekeeping' genes, before (left) and after (right) TPEA-PCR. After cDNA amplification, expression of all eight genes assessed could be detected after carrying out gene-specific PCR on only 5% of the amplified cDNA generated from a single cell. Without cDNA amplification, however, expression of none of the genes assayed could be detected at this level. Genes assayed: RPL21 (riboprotein L21), RP27a (riboprotein 27a), RPL28 (riboprotein L28), RPS5 (riboprotein S5), HSKPQZ7 (housekeeping protein), ACTB (Cytoplasmic beta-actin), G-3-PDH (Glyceraldehyde-3-phosphate dehydrogenase) and EF1 (Elongation factor 1). Problems that can be inherent in the detection of gene expression at the single cell level include the possibility of contamination by other cell types and the interpretation of negative results, which could be due to cell heterogeneity or an artefact of the detection system. In order to minimize these concerns, we chose to develop the system on sorted lymphoblastoma cells. These cells might be expected to be relatively homogeneous, so allowing an assessment of assay reproducibility by comparison of expression profiles between different cells and would also avoid contamination by other cell types. The use of sorted cells also provides a ready source of single cell material. Figure 3. Multiple gene expression analysis in single lymphoblastoma cells. Four cells (1-4) were lysed, reverse transcribed, the cDNA amplified and gene-specific PCR performed on the product. CD2, SI and intron primer pairs serve as negative controls to check for genomic contamination of the samples. Eleven of the other genes are expressed in all four cells in duplicate, while JUND, EF1 CDC25B and CD19 expression was not consistent between cells. Genes assayed: RPL5 (riboprotein L5), RPL21 (riboprotein L21), RP27a (riboprotein 27a), RPL28 (riboprotein L28), RPS5 (riboprotein S5), RPS9 (riboprotein S9), RPS10 (riboprotein S10), RPS29 (riboprotein S29), HSKPQZ7 (Housekeeping protein), ACTB (Cytoplasmic beta-Actin), G-3-PDH (Glyceraldehyde-3-phosphate dehydrogenase), EF1 (Elongation factor 1), JUND (JUND), CDC25B (cell cycle factor CDC25b) and the cell surface antigens CD19, CD79a, CD2, IGM (immunoglobulin IgM), SI (the intestine-specific enzyme, sucrase-isomaltase). Figure Figure 4. Demonstration of adenosine 2a receptor expression in striatal cholinergic interneurons. (a) An infrared video image of a rat striatal cholinergic interneuron during electrophysiological characterization and (b) after aspiration of cytoplasm. The expression of four housekeeping genes, the transmitter synthesising enzymes choline acetyltransferase (found only in cholinergic neurons) and glutamic acid decarboxylase (found in GABAergic, medium spiny neurons), three tachykinin (NK) receptors and the adenosine A2a receptor was assessed in 26 striatal cholinergic neurons. Two representative neurons are shown; 1 expresses the A2a receptor; 2 does not, while the expression of the other genes tested are the same in both. The power of this technique lies in its potential to facilitate expression profiling of cells derived from complex cell populations, even when they form only a small proportion of the population as a whole, and in its ability to detect low abundance transcripts. We therefore investigated the expression of the neurokinin (NK) receptors and the adenosine A2a receptor in single striatal cholinergic interneurons which constitute a small fraction of the total cellular mass of the striatum. Of the neurokinin receptors, the NK1 receptor is widely accepted as being expressed in these cells (20), so the expression of this gene was examined as an example of an mRNA species expressed at far lower levels than housekeeping genes. In contrast there is considerable controversy as to whether the A2a receptor is expressed in cholinergic interneurons (15,21-23), suggesting that the corresponding mRNA species may be present at low levels as suggested by one in situ hybridisation study (24) or not present at all, as suggested by other studies (25). Expression of the A2a receptor in these cells has important implications for the mechanism of action of adenosine and because of the potential of the A2a receptor to act as a target for novel drugs for the amelioration of Parkinson's Disease (22). Using a patch pipette it was possible to harvest upwards of an estimated 40% of the cellular contents of these neurons (Fig. An important advantage of TPEA-PCR is the relative ease of the methodology. Other protocols can be time-consuming and complex, involving DNA purification and precipitation from one step to another. The current cDNA amplification technique, however, can be carried out in a single tube with a need for only limited manual intervention. This therefore makes it possible to amplify large numbers of samples relatively easily. The ability to then analyse the expression of many genes of unrelated sequence, both at high and low abundance, in samples taken from as little as a single cell, will potentially allow it to be used in high throughput screening systems. We have demonstrated that this approach can be used to analyse gene expression in samples from as little as a single cell (Fig. In summary, this report describes a rapid, robust and reproducible procedure, called TPEA-PCR, capable of amplifying 3[prime] fragments of cDNA prior to analysis by other techniques. We believe that this will greatly facilitate the analysis of gene expression in samples obtained from small samples of tissue or single cells. In so doing, it will allow the utilisation of the wealth of new sequence data now available, to further our understanding of disease processes and the cellular physiology of complex tissues.
ACKNOWLEDGEMENTS
We would like to thank Dr David Bentley for all his support and constructive criticism. This work has been funded by the Wellcome Trust and Parke-Davis Neuroscience Research Centre. This technique is the subject of UK Patent application: 98 17055.8.
REFERENCES
This article has been cited by other articles:
This page is run by Oxford University Press, Great Clarendon Street, Oxford OX2 6DP, as part of the OUP Journals
Comments and feedback: www-admin{at}oup.co.uk
Last modification: 19 Sep 1998
Copyright©Oxford University Press, 1998.
![]()
CiteULike
Connotea
Del.icio.us What's this?
![]()
![]()

![]()
![]()
![]()
C. Brochier, M.-C. Gaillard, E. Diguet, N. Caudy, C. Dossat, B. Segurens, P. Wincker, E. Roze, J. Caboche, P. Hantraye, et al.
Quantitative gene expression profiling of mouse brain regions reveals differential transcripts conserved in human and affected in disease models
Physiol Genomics,
April 1, 2008;
33(2):
170 - 179.
[Abstract]
[Full Text]
[PDF]
![]()
![]()
![]()

![]()
![]()
![]()
M. Fry, P. M. Smith, T. D. Hoyda, M. Duncan, R. S. Ahima, K. A. Sharkey, and A. V. Ferguson
Area postrema neurons are modulated by the adipocyte hormone adiponectin.
J. Neurosci.,
September 20, 2006;
26(38):
9695 - 9702.
[Abstract]
[Full Text]
[PDF]
![]()
![]()
![]()

![]()
![]()
![]()
A. Blaschitz, H. Juch, A. Volz, H. Hutter, C. Daxboeck, G. Desoye, and G. Dohr
The soluble pool of HLA-G produced by human trophoblasts does not include detectable levels of the intron 4-containing HLA-G5 and HLA-G6 isoforms
Mol. Hum. Reprod.,
October 1, 2005;
11(10):
699 - 710.
[Abstract]
[Full Text]
[PDF]
![]()
![]()
![]()

![]()
![]()
![]()
S. P. Brandt
Microgenomics: gene expression analysis at the tissue-specific and single-cell levels
J. Exp. Bot.,
February 1, 2005;
56(412):
495 - 505.
[Abstract]
[Full Text]
[PDF]
![]()
![]()
![]()

![]()
![]()
![]()
M. Toledo-Rodriguez, B. Blumenfeld, C. Wu, J. Luo, B. Attali, P. Goodman, and H. Markram
Correlation Maps Allow Neuronal Electrical Properties to be Predicted from Single-cell Gene Expression Profiles in Rat Neocortex
Cereb Cortex,
December 1, 2004;
14(12):
1310 - 1327.
[Abstract]
[Full Text]
[PDF]
![]()
![]()
![]()

![]()
![]()
![]()
A. Peixoto, M. Monteiro, B. Rocha, and H. Veiga-Fernandes
Quantification of Multiple Gene Expression in Individual Cells
Genome Res.,
October 1, 2004;
14(10a):
1938 - 1947.
[Abstract]
[Full Text]
[PDF]
![]()
![]()
![]()

![]()
![]()
![]()
M. van den Top, M. F Nolan, K. Lee, P. J Richardson, R. M Buijs, C. H Davies, and D. Spanswick
Orexins induce increased excitability and synchronisation of rat sympathetic preganglionic neurones
J. Physiol.,
June 15, 2003;
549(3):
809 - 821.
[Abstract]
[Full Text]
[PDF]
![]()
![]()
![]()

![]()
![]()
![]()
A. Volz, A. Ehlers, R. Younger, S. Forbes, J. Trowsdale, D. Schnorr, S. Beck, and A. Ziegler
Complex Transcription and Splicing of Odorant Receptor Genes
J. Biol. Chem.,
May 23, 2003;
278(22):
19691 - 19701.
[Abstract]
[Full Text]
[PDF]
![]()
![]()
![]()

![]()
![]()
![]()
L. Smith and A. Greenfield
DNA microarrays and development
Hum. Mol. Genet.,
April 2, 2003;
12(90001):
R1 - 8.
[Abstract]
[Full Text]
[PDF]
![]()
![]()
![]()

![]()
![]()
![]()
J. M. Gale, C. P. Romero, G. B. Tafoya, and J. Conia
Application of Optical Trapping for Cells Grown on Plates: Optimization of PCR and Fidelity of DNA Sequencing of p53 Gene from a Single Cell
Clin. Chem.,
March 1, 2003;
49(3):
415 - 424.
[Abstract]
[Full Text]
[PDF]
![]()
![]()
![]()

![]()
![]()
![]()
L. Smith, P. Underhill, C. Pritchard, Z. Tymowska-Lalanne, S. Abdul-Hussein, H. Hilton, L. Winchester, D. Williams, T. Freeman, S. Webb, et al.
Single primer amplification (SPA) of cDNA for microarray expression analysis
Nucleic Acids Res.,
February 1, 2003;
31(3):
e9 - e9.
[Abstract]
[Full Text]
[PDF]
![]()
![]()
![]()

![]()
![]()
![]()
S. Ivashuta, K. Uchiyama, M. Gau, and Y. Shimamoto
Linear amplification coupled with controlled extension as a means of probe amplification in a cDNA array and gene expression analysis during cold acclimation in alfalfa (Medicago sativa L.)
J. Exp. Bot.,
February 1, 2002;
53(367):
351 - 359.
[Abstract]
[Full Text]
[PDF]
![]()
![]()
![]()

![]()
![]()
![]()
K Lee, A K Dixon, I Gonzalez, E B Stevens, S McNulty, R Oles, P J Richardson, R D Pinnock, and L Singh
Bombesin-like peptides depolarize rat hippocampal interneurones through interaction with subtype 2 bombesin receptors
J. Physiol.,
August 1, 1999;
518(3):
791 - 802.
[Abstract]
[Full Text]
[PDF]
![]()
![]()
![]()

![]()
![]()
![]()
K Lee, A K Dixon, P J Richardson, and R D Pinnock
Glucose-receptive neurones in the rat ventromedial hypothalamus express KATP channels composed of Kir6.1 and SUR1 subunits
J. Physiol.,
March 1, 1999;
515(2):
439 - 452.
[Abstract]
[Full Text]
[PDF]
![]()
![]()
![]()

![]()
![]()
![]()
K. J. Alden, P. H. Goldspink, S. W. Ruch, P. M. Buttrick, and J. Garcia
Enhancement of L-type Ca2+ current from neonatal mouse ventricular myocytes by constitutively active PKC-beta II
Am J Physiol Cell Physiol,
April 1, 2002;
282(4):
C768 - C774.
[Abstract]
[Full Text]
[PDF]
![]()
This Article ![]()
![]()
Abstract
![]()
Print PDF (315K)
![]()
Alert me when this article is cited
![]()
Alert me if a correction is posted
![]()
Services ![]()
![]()
Email this article to a friend
![]()
Similar articles in this journal
![]()
Similar articles in ISI Web of Science
![]()
Similar articles in PubMed
![]()
Alert me to new issues of the journal
![]()
Add to My Personal Archive
![]()
Download to citation manager
![]()
Search for citing articles in:
ISI Web of Science (67)
![]()
Request Permissions ![]()
Commercial Re-use Guidelines
for Open Access NAR Content
![]()
Google Scholar ![]()
![]()
Articles by Dixon, A. K.
![]()
Articles by Freeman, T. C.
![]()
Search for Related Content
![]()
PubMed ![]()
![]()
PubMed Citation
![]()
Articles by Dixon, A. K.
![]()
Articles by Freeman, T. C.
![]()
Social Bookmarking ![]()
![]()
What's this?

