|
|
|
|
Published online before print
June 29, 2006, 10.1101/gr.4999906 Genome Res. 16:1041-1045, 2006 ©2006 by Cold Spring Harbor Laboratory Press; ISSN 1088-9051/06 $5.00 OPEN ACCESS ARTICLE
Methods RNA expression profiling at the single molecule level1 Biophysics Institute, Johannes Kepler University Linz, A-4040 Linz, Austria; 2 Division of Genomics, Department of Molecular Biology, University of Salzburg, A-5020 Salzburg, Austria; 3 Center for Biomedical Nanotechnology, Upper Austrian Research GmbH, A-4020 Linz, Austria; 4 Department of Knowledge-based Mathematical Systems, Johannes Kepler University Linz, A-4040 Linz, Austria
We developed a microarray platform for PCR amplification-independent expression profiling of minute samples. A novel scanning system combined with specialized biochips enables detection down to individual fluorescent oligonucleotide molecules specifically hybridized to their complementary sequence over the entire biochip surface of cm2 size. A detection limit of 1.3 fM target oligonucleotide concentrationcorresponding to only 39,000 molecules in the sample solutionand a dynamic range of 4.7 orders of magnitude have been achieved. The applicability of the system to PCR amplification-independent gene-expression profiling of minute samples was demonstrated by complex hybridization of cDNA derived from the equivalent of only 104 cells, which matches results obtained in ensemble studies on large samples. By counting each hybridized molecule on the microarray, the method is insusceptible to gene-specific variations of the labeling, thereby representing a principle advance to conventional ensemble-based microarray analysis.
Microarray analysis of the global transcriptional state of a small cell number represents a major target for both medical diagnostics and basic research. In particular, precise characterization of minute patients samples such as malignant and multidrug-resistant cancer stem cells is decisive for a detailed molecular understanding of the oncogenic mechanisms driving tumor growth as well as for diagnosis and individualized therapy; furthermore, global molecular characterization of small subregions within heterogeneous tissues or of tiny biological units such as stem-cell niches is of fundamental importance for a more defined understanding of key biological processes such as tissue homeostasis in normal and diseased states.
Recent technological advances allow for routine isolation of small, yet defined subpopulations of cells such as cancer stem cells from heterogeneous fluid or tissue samples using Fluorescence Activated Cell Sorting (FACS) and Laser Capture Microdissection (LCM). Subsequent purification and biochemical processing can be performed with minimum loss of material via microfluidic devices (Paegel et al. 2003
In conventional DNA microarray readout, the sensitivity is limited by standard formats of biochip substrates. Their thickness of
Setting up a robust scanning technology for cm2 areas with single-molecule sensitivity is a challenging task. To reliably identify single dye molecules, imaging at the resolution limit is inevitable; the concomitant short depth of focus of We report here the development of a system for PCR amplification-independent DNA microarray analysis down to the single-molecule level, which includes an ultrasensitive detection device and accessory biochips. Hybridization of fluorescent oligonucleotides was used to characterize the platform in terms of sensitivity and dynamic range. Complex hybridization of biological material corresponding to only 104 cells demonstrates the applicability of the system for biomedical diagnostics.
Biochip surfaces of arbitrary size were imaged at the sensitivity of single-dye molecules using a home-built scanning system (Hesse et al. 2004
We first tested the binding specificity by hybridizing a Cy5-labeled 60mer target oligonucleotide (c = 100 fM) to its complementary sequence (probe oligonucleotide) bound to the biochip surface (Fig. 1). Homogenously distributed diffraction-limited peaks were observed over the entire image; the corresponding surface density was 1000-fold higher when compared with a surface functionalized with control oligonucleotide. To ascertain whether each fluorescence peak corresponds to a single Cy5-oligonucleotide molecule, we determined the brightness of individual peaks as a measure of the number of colocalized dye molecules; fitting with a two-dimensional Gaussian profile (Schmidt et al. 1996
To determine the detection limit and the dynamic range of the platform, we varied the target concentration c from 100 aM to 1 nM. Hybridization experiments were performed on microarrays with a regular checkerboard pattern of probe and control oligonucleotides using a spot size of 100 x 100 µm2. Figure 2A shows a microarray hybridized with 100 fM Cy5-labeled target oligonucleotide. Binding efficiency was homogenous over the chip surface, characterized by a spot-to-spot signal variation of only 13%. Single oligonucleotide molecules bound to the chip surface become visible when the data are plotted at higher magnification (Fig. 2B). At low surface densities of N 1000 molecules per spot, more than 95% of all molecules are well separated in the image and were counted via automated algorithms. At higher surface densities, peaks begin to overlap; in such cases, the number of hybridized molecules per spot was inferred from the corresponding integrated fluorescence signal and the single molecule brightness B. Residual signals due to unspecific adsorption (Nunspec) on control spots were subtracted for calculating the specific signal Nspec. A constant concentration of Cy3-labeled target oligonucleotide (c = 10 fM) was included in all experiments; the fluorescence signal SCy3 is a measure of the binding capacity of each spot and was used to correct for chip-to-chip variations with (SCy3(c)/ SCy3 all conc.) as linear normalization. Figure 2C shows the binding curve for specific hybridization (red) and unspecific adsorption (black). A linear relation between specifically hybridized oligonucleotides and sample concentration was found over six orders of magnitude. At all applied concentrations, the average distance between two neighboring molecules is >30 nm; therefore, fluorescence quenching has no significant effect on the molecular brightness.
A constant unspecific signal of Nunspec = 13.2 molecules per spot was observed for c < 1 pM, due to residual surface and buffer impurities. At sample concentrations c 1 pM unspecific oligonucleotide binding resulted in a linear increase of the unspecific signal, approximately four orders of magnitude below the specific signal.
For the definition of the detection limit, two major noise contributions have to be considered; first, the signal noise, which includes variations in the brightness of the fluorophore and of the background; and second, variations due to the stochastic nature of the binding process. While for conventional microarray analysis the first term is dominant and limits the accuracy of signal quantification, its influence on single-molecule counting is negligible due to the high-detection reliability. Here, the ultimate statistical limit of the binding process defines the sensitivity of the system. Assuming Poissonian-binding statistics, the signal-to-noise ratio (SNR) is described by SNR = Nspec/ The high instrumental dynamic range of the developed system enables signal quantification in a range of from 1.3 fM to 1 nM. The operational dynamic range, however, is currently limited by unspecific adsorption of fluorescent oligonucleotides to the target control spots. At c = 67 pM, the unspecific signal equals the specific signal recorded at the detection limit of c = 1.3 fM, yielding an operational dynamic range of 4.7 orders of magnitude.
To demonstrate the applicability of the system to gene-expression profiling, we used arrays of 125 distinct 6770mer oligonucleotides spotted in eight replicates (1000 probe features per array). Complex cDNA was synthesized from 200 ng of total RNA from the equivalent of only 104 cells of a human keratinocyte cell line (HaCaT) and used for hybridization (Fig. 3A). Samples were scanned at a speed of 2.4 sec/spot; in total, the analysis took
For validation, we compared the expression profiles obtained with the platform described above to those of a conventional microarray platform. A 100-fold higher amount of cDNA from the same RNA preparation was hybridized to commercially available aldehyde substrates (Schott Nexterion) under identical conditions and scanned with a conventional confocal microarray scanner (ScanArray 4000, Perkin Elmer). The signals derived from both measurements are in good agreement, yielding a Pearson correlation coefficient of 0.85. The high consistency of the data obtained here using different microarray platforms and different sample amounts appears remarkable, as typically a considerable divergence of gene-expression results has been observed when comparing commercial products (Tan et al. 2003 0.910.98). The intermediate result revealed in our study reflects the fact that different platforms containing the same oligoset have been investigated.
The system presented enables expression profiling of minute amounts of sample material. Material from the equivalent of only 104 cells can be reliably analyzed without prior PCR amplification steps. Connection to sample isolation and microfluidic purification techniques offers the exciting perspective to directly analyze minute cell populations in biomedical applications. In addition to sensitivity, 4.7 orders of magnitude dynamic range significantly exceeds specifications of conventional microarray platforms and approaches the complete biological range of mRNA expression levels of up to six orders of magnitude (Holland 2002
Single molecule biochip reader The detection system was set up on an epifluorescence microscope (Axiovert 200, Zeiss) (Hesse et al. 2004 -Fluar, Zeiss). After appropriate filtering using standard Cy3 and Cy5 filter sets (Chroma Technology Corp.), fluorescence was imaged onto a back-illuminated CCD camera (SPEC10:100B, Princeton Instruments; quantum efficiency = 90%, gain = 0.77counts/e-). Biochip readout was performed in time-delay and integration (TDI) mode. For this, samples were shifted using a motorized xy-stage (Scan IM 120 x 100, Märzhäuser) synchronized to the line-shift of the camera. During the TDI-scan, data were streamed into the RAM of the PC and stored on the hard disk in the course of repositioning for the next scan. The reader was equipped with an automated focus hold system operating during the scanning process; the back-reflected laser beam is imaged on a two-segment photodiode, and the differential signal is used to control a z-piezo (PIFOC, Physik Instrumente) for fast refocusing.
Microarray preparation Target oligonucleotides labeled with Cy3 or Cy5 at the 5' end (Cy3/Cy5-GGTGGCACAAGCGCTACAGCTCGGTTTC AATCTGCTAAGAGATCCACGCTGCATTTCTAT) (VBC Genomics) were used for evaluation of sensitivity and dynamic range of the developed platform. For RNA profiling, 96 oligonucleotides of the Human Genome Oligo Set V3 (Operon) and 29 oligonucleotides designed using OligoWiz 1.0 (http://www.cbs.dtu.dk/services/OligoWiz/) were printed. Oligonucleotides were dissolved in 3xSSC/1.5 M betaine to a final concentration of 20 µM. After printing, the slides were stored under dry conditions for 1 d at room temperature. Stable covalent coupling to the chip surface and removal of unbound probe oligonucleotides was achieved by extensive washing (1 h) with blocking buffer (0.1 M Na-carbonate, 0.1 M 4-aminobutanoic acid, 10 µM NaBH3CN at pH 8.0), followed by a washing step with 1xSSC buffer. Directly after washing, slides were used for hybridization experiments.
RNA preparation
Hybridization Conventional microarrays were hybridized with labeled cDNA synthesized from 20 µg of total RNA. Following precipitation with EtOH, cDNA was dissolved in 3xSSC/0.1% SDS, heated to 96°C for 3 min, and hybridized overnight at 50°C in a semiautomated hybridization station (SlideBooster, Advalytix AG) at a mixing:pause ratio of 3:7, mixing power: 27. Slides were washed at 30°C in a semiautomated wash station (Advalytix AG) two times for 10 min in 2xSSC/0.1% SDS, two times for 10 min in 0.2xSSC/0.1%SDS, and two times for 5 min in 0.2xSSC. Slides were dried by spinning for 3 min at 250g and subsequently scanned on a ScanArray 4000 confocal microarray scanner (Perkin Elmer).
Data analysis
Diffraction-limited peaks are small, approximately symmetric features over the background noise. For counting the total number of peaks per spot, N, we used the à trous wavelet filter bank method as described previously (Starck et al. 1998
b0, where b0 = [1/16,1/4,3/8,1/4,1/16]. Low-pass filtering is iterated with kernels bj, where bj+1 is obtained by inserting lines and columns of zeros between the lines and columns of bj, yielding
j. We estimate the variance of the noise via robust statistics methods j2 = (x,ymedian|Wj(x,y) m|)2, where m is the median of Wj(x, y). j allows separation of the significant wavelet coefficients from those due to noise:
This study was supported by the GEN-AU program of the Austrian Federal Ministry of Education, Science, and Culture, by the Austrian Research Fund, by the state of Upper Austria, and by the Priority Program "Life Sciences and Health" of the University of Salzburg.
5 Corresponding authors. E-mail gerhard.schuetz{at}jku.at; fax 43-732-3468-29284.
E-mail annemarie.frischauf{at}sbg.ac.at; fax 43-662-8044-183. [Supplemental material is available online at www.genome.org.] Article published online before print. Article and publication date are at http://www.genome.org/cgi/doi/10.1101/gr.4999906
Baugh L.R., Hill A.A., Brown E.L., Hunter C.P. 2001. Quantitative analysis of mRNA amplification by in vitro transcription. Nucleic Acids Res. 29: E29. Hesse J., Sonnleitner M., Sonnleitner A., Freudenthaler G., Jacak J., Höglinger O., Schindler H., Schutz G.J. 2004. Single-molecule reader for high-throughput bioanalysis. Anal. Chem. 76: 59605964.[Medline] Holland M.J. 2002. Transcript abundance in yeast varies over six orders of magnitude. J. Biol. Chem. 277: 1436314366. Hong J.W., Studer V., Hang G., Anderson W.F., Quake S.R. 2004. A nanoliter-scale nucleic acid processor with parallel architecture. Nat. Biotechnol. 22: 435439.[CrossRef][Medline] Hu L., Wang J., Baggerly K., Wang H., Fuller G.N., Hamilton S.R., Coombes K.R., Zhang W. 2002. Obtaining reliable information from minute amounts of RNA using cDNA microarrays. BMC Genomics 3: 16.[CrossRef][Medline] Irizarry R.A., Warren D., Spencer F., Kim I.F., Biswal S., Frank B.C., Gabrielson E., Garcia J.G.N., Geoghegan J., Germino G. et al. 2005. Multiple-laboratory comparison of microarray platforms. Nat. Methods 2: 345349.[CrossRef][Medline] Kacharmina J.E., Crino P.B., Eberwine J. 1999. Preparation of cDNA from single cells and subcellular regions. Methods Enzymol. 303: 318.[Medline] Mahadevappa M. and Warrington J.A. 1999. A high-density probe array sample preparation method using 10- to 100-fold fewer cells. Nat. Biotechnol. 17: 11341136.[CrossRef][Medline] Olivo-Marin J.C. 2002. Extraction of spots in biological images using multiscale products. Pattern Recognit. 35: 19891996.[CrossRef] Paegel B.M., Blazej R.G., Mathies R.A. 2003. Microfluidic devices for DNA sequencing: Sample preparation and electrophoretic analysis. Curr. Opin. Biotechnol. 14: 4250.[CrossRef][Medline] Schlapak R., Pammer P., Armitage D., Zhu R., Hinterdorfer P., Vaupel M., Frühwirt T., Howorka S. 2005. Glass surfaces grafted with high-density poly(ethylene glycol) as substrates for DNA oligonucleotide microarrays. Langmuir 22: 277285.[CrossRef] Schmidt T., Schutz G.J., Baumgartner W., Gruber H.J., Schindler H. 1996. Imaging of single molecule diffusion. Proc. Natl. Acad. Sci. 93: 29262929. Sonnleitner M., Freudenthaler G., Hesse J., Schütz G.J. 2005. High-throughput scanning with single-molecule sensitivity. Proc. SPIE 5699: 202210.[CrossRef] Starck J.-L., Murtagh F., Bijaoui A. 1998. Image processing and data analysisThe multiscale approach. Cambridge University Press, Cambridge, UK. Tan P.K., Downey T.J., Spitznagel E.L., Xu P., Fu D., Dimitrov D.S., Lempicki R.A., Raaka B.M., Cam M.C. 2003. Evaluation of gene expression measurements from commercial microarray platforms. Nucleic Acids Res. 31: 56765684. Taylor B.T., Nambiar P.R., Raja R., Cheung E., Rosenberg D.W., Anderegg B. 2004. Microgenomics: Identification of new expression profiles via small and single-cell sample analyses. Cytometry A 59A: 254261.[CrossRef][Medline] Van Gelder R.N., von Zastrow M.E., Yool A., Dement W.C., Barchas J.D., Eberwine J.H. 1990. Amplified RNA synthesized from limited quantities of heterogeneous cDNA. Proc. Natl. Acad. Sci. 87: 16631667. Wang E., Miller L.D., Ohnmacht G.A., Liu E.T., Marincola F.M. 2000. High-fidelity mRNA amplification for gene profiling. Nat. Biotechnol. 18: 457459.[CrossRef][Medline]
Received December 2, 2005; accepted in revised format May 4, 2006. This article has been cited by other articles:
|
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||