A highly sensitive and specific system for large-scale gene expression profiling

Background Rapid progress in the field of gene expression-based molecular network integration has generated strong demand on enhancing the sensitivity and data accuracy of experimental systems. To meet the need, a high-throughput gene profiling system of high specificity and sensitivity has been developed. Results By using specially designed primers, the new system amplifies sequences in neighboring exons separated by big introns so that mRNA sequences may be effectively discriminated from other highly related sequences including their genes, unprocessed transcripts, pseudogenes and pseudogene transcripts. Probes used for microarray detection consist of sequences in the two neighboring exons amplified by the primers. In conjunction with a newly developed high-throughput multiplex amplification system and highly simplified experimental procedures, the system can be used to analyze >1,000 mRNA species in a single assay. It may also be used for gene expression profiling of very few (n = 100) or single cells. Highly reproducible results were obtained from duplicate samples with the same number of cells, and from those with a small number (100) and a large number (10,000) of cells. The specificity of the system was demonstrated by comparing results from a breast cancer cell line, MCF-7, and an ovarian cancer cell line, NCI/ADR-RES, and by using genomic DNA as starting material. Conclusion Our approach may greatly facilitate the analysis of combinatorial expression of known genes in many important applications, especially when the amount of RNA is limited.


Background
Biological processes are underlain by interactions between various genes and their products through defined pathways in the molecular network, in which molecules cross communicate in hitherto unknown ways under both healthy and disease conditions. Learning gene expression patterns on a genomic scale would substantially help deconvolute these complex processes. Exhaustive identification of human genes during the Human Genome Project has made such studies possible. By global gene expression profiling in cells and tissues under either physiological or in vitro conditions, our understanding of the correlation between gene functions and their phenotypic effects could be significantly enhanced.
The advent of the microarray-based high-throughput RNA detection system [1,2] has made it possible to profile gene expression patterns for the entire transcriptome. However, to detect gene transcripts very specifically, one needs to discriminate them from closely related sequences including: (1) the corresponding gene sequences. Although contamination of gene sequences may not be a concern for applications using purified mRNA, gene sequences must be taken into consideration for applications directly using cell lysate without RNA extraction. This becomes especially important when the studied transcripts are present at low abundance; (2) pseudogenes and their possible transcripts. The number of pseudogenes in the human genome was estimated to be 20,000 to 33,000, which are widely expressed [3,4]. These sequences usually share a high degree of sequence identity with the closely related genes; (3) unprocessed RNA containing the same exons as those of the corresponding mRNA. So far, no system has addressed the above issue very effectively.
Among the microarray-based platforms, GeneChip is a commonly used system and has been improved significantly since it was invented, and has contributed to understanding the complex gene expression network in a great deal. However, since this technology is limited by its high degree of nonspecificity and insensitivity, its application has been limited in molecular network integration. Results from a recent analysis [5] indicated that on the Affymetrix GeneChip U95A/Av2 array, 20,696 (10.5%) probes were nonspecific, which could cross-hybridize to multiple genes, and 18,363 (9.3%) probes missed the target transcript sequences. The numbers of nonspecific and mis-targeted probes on the U133A array were comparable, which were 29,405 (12.1%) and 19,717 (8.0%), respectively [5]. These ~20% of problematic probes certainly and substantially compromise the data accuracy, decrease the value of microarray data, and are not acceptable for the studies of molecular network integration. It was also found that some probe sets representing the same genes on Affymetrix microarrays could show significant discrepancy because of the non-specific hybridization [6,7] In most applications, gene expression profiling with microarrays including GeneChip requires amplification of sample RNA, regardless of how much material is available. Normally, 1 to 3 μg of RNA is required for each assay [8]. However, high-throughput gene expression profiling with superior sensitivity is becoming more and more demanded, and has its wide applications. For example, in breast cancer research, analysis of specimens from microdissection may provide important information about genes involved in different cancer development stages and for understanding the molecular mechanisms underlying cancer development [9]. Specimens from fine needle biopsy are also important in diagnostic procedures and in evaluating therapeutic effects. The ability to analyze a large number of genes in single cells may help understand the origin and clonality of cancer development and learn the molecular details involved in different stages of the cell cycle.
Current methodologies for gene expression profiling in small RNA samples, especially those from single cells, are very limited. Many of these protocols [2,10,11] require multiple enzymatic reactions that may seriously reduce the sensitivity and compromise the specificity. RNA preparation in most of applications also involves a number of steps, which is rather lengthy, tedious, and requires highly skilled personnel.
To solve the above problems, we have developed a highly specific and sensitive gene expression profiling system. With this system, primers are specially designed to amplify mRNA sequences very specifically. Probes used for microarray detection are designed only to hybridize to sequences amplified from mRNA. In conjunction with the high-throughput multiplex amplification protocol developed in our laboratory lately [12], a large number of mRNA species directly released from very few cells or even single cells can be amplified to a detectable amount without RNA isolation. Amplified products can then be detected by the single-base extension [13] assay on an oligonucleotide microarray [14].

Experimental system used in the study
To establish a cancer gene expression array, a panel of cancer-related genes were selected based on their known functions and/or cancer-associated expression patterns from published literature [15][16][17][18][19][20][21][22][23][24][25][26][27][28]. All amplicon sequences were subjected to computational screening to ensure their uniqueness. Primers and probes were selected according to a series of criteria as specified in Materials and Methods. Most primer pairs amplify sequences in two neighboring exons separated by large introns. The intron lengths ranged from 79 bp to 90 kb with an average of 2.0 kb and 97% of the introns are longer than 200 bp. Initially 1,445 genes were used as the input for the primer and probe design program. Primers and probes were selected for 1,120 (77.5%) of these genes. The remaining 22.5% had either no introns or no suitable sequences for primers and/or probes. Fifteen of these remaining genes with important functions in cancer development were included in the panel. Primers and probes were designed based on the unique sequences in these genes, and were not required to have introns internally located within the amplified sequences. Therefore, a total of 1,135 genes were included in our multiplex assay. (Details about these genes, and their corresponding primers and probes used for the study are listed in Additional files 1 and 2.) Microarray-based single-base extension (SBE) assay has been used to genotype single nucleotide polymorphisms (SNPs) [12,29,30] in our laboratory. In the present study, SBE was adapted for gene expression profiling. To simplify the analysis, all probes were designed to terminate immediately before a 'G' base in the templates. In this way, the probes were extended by a single base, dideoxynucleoside triphosphate (ddCTP) that was fluorescently labeled. By using one color, the bias associated with different dyes was also eliminated. The detection procedure is schematically illustrated in Fig. 1. Resulting data have been deposited to the NCBI's Gene Expression Omnibus (GEO) [31] and are accessible through GEO Series accession number GSE5920.

Reproducibility of the high-throughput gene expression profiling system
To test the reproducibility of our system, gene expression was profiled for three duplicated 100-cell samples from an ovarian cancer cell line, NCI/ADR-RES [32] and two 100-cell samples from a breast cancer cell line, MCF-7.
Resulting microarray data are supplied in Additional file 3. Table 1 summarizes the numbers of gene transcripts detected from different samples. As shown, 660 (58.2%), 663 (58.4%), and 662 (58.3%) gene transcripts were detected from the three 100-cell samples of NCI/ADR-RES, respectively. Of these transcripts, 650 (>98%) were detected from all three duplicates. Signal intensities for the 1,135 genes were strongly correlated between the duplicates (Pearson's r = 0.977, 0.974, and 0.949, respectively). Fig. 2A shows a scatter plot of two duplicates. Of the 650 transcripts detected in all three NCI/ADR-RES 100-cell samples, only 6 (0.9%), 17 (2.6%), and 1 (0.2%) transcripts had their signal intensities differing by >2 fold between each two of these three duplicates. Twenty-six transcripts were detected from only one or two of the three samples. The signal intensities for these transcripts were low. Only one transcript in one sample had its signal intensity >1,000, indicating that the inconsistence among the duplicates was due to low signals of these transcripts.
For the two 100-cell samples from MCF-7, 615 (54.2%) and 614 (54.1%) gene transcripts were detected, respectively, with 597 (>97%) detected in both. Of these 597 transcripts, 562 (94.1%) had signal intensities differing less than two fold. Similar to the situation with NCI/ADR-Schematic illustration of the high-throughput gene expression profiling procedure Figure 1 Schematic illustration of the high-throughput gene expression profiling procedure. Fluorescent labeling is indicated by an asterisk. RES samples, all 34 transcripts that were detected in only one sample but not the other had low signal intensities with only nine genes whose signal intensities were >1,000 in one of the two samples.
Because samples prepared from a large number of cells are usually associated with high reliability, we further compared the microarray results of the NCI/ADR-RES 100-cell samples with those from a 10,000-cell sample of the same cell line. Resulting data also show a high degree of correlation (r = 0.961, Figure 2B). As shown in Table 1, 630 (96.7%) of the 650 gene transcripts detected from all the 100-cell samples were also detected from the 10,000-cell sample. Sixty-three gene transcripts were detected in at least one of the three 100-cell samples but not in the 10,000-cell sample, or vice versa. Of these 63 gene transcripts, 61 had signal intensities below 1,000 in all the three 100-cells. However, the change from 100 to 10,000 cells did enhance the detection of 21 gene transcripts whose signal intensities were >2 fold greater in the 10,000-cell sample than those in the 100-cell samples. Among these 21 transcripts, six had signal intensities in

Undetectable in all 442
Detectable in all 531

Detectable in some 212
Undetectable in all 392

Detectable in all non-single, not in single 61
Detectable in all single, not in non-single 27 *ADR: NCI/ADR-RES the 10,000-cell sample more than 15 fold greater than the average intensities of the corresponding genes in the three 100-cell samples, indicating that using 10,000 cells may have significantly increased the copy numbers of these transcripts or changed their absence status to presence. These data indicate that our system not only can produce very reliable results even with as few as 100 cells but also is very sensitive to the copy number change for the lowcopy-number gene transcripts.

Sensitivity of the high-throughput gene expression profiling system
To further test the sensitivity of our high-throughput gene expression profiling system, single NCI/ADR-RES cell samples were prepared and used for multiplex gene expression assay of the 1,135 mRNA species. Microarray results from three of these samples are listed in Additional file 3. The numbers of gene transcripts detected from the three single-cell samples were 590, 576, and 614, respectively. Of these transcripts, 503 were detected from all single cells. Of the 503, 463 (92.0%) were also detected from all non-single-cell (100-cell and 10,000-cell) samples, indicating a prevalent expression of these genes in most, if not all, cells at relatively high levels.
On the other hand, the detection range of gene transcripts from the three single-cell samples was wider compared to the non-single-cell samples. As shown in Table 1, 449 transcripts were undetectable in all three single cell samples, a number which is not greater than that (459) for the three 100-cell samples and is comparable to that (442) for all non-single-cell samples. The number of undetectable gene transcripts from all single and non-single cell samples is 357. This number means that from single cells, we not only detected a comparable number of genes, but also detected a new set of 449-357 = 92 genes that could not be detected with non-single-cell samples of the same cell line.
The robustness of gene expression profiling with singlecell samples was also demonstrated by the signal intensities. As described above, most transcripts that were detected from some but not all non-single-cell samples had low signal intensities and very few were >1,000. The scenario with single cells is very different. Of the 503 gene transcripts detected from all single cells, 40 were detected in one to three non-single-cell samples but not all four. All 40 but one have signal intensity >1,000 in at least one of the three single-cell samples. Of the 183 transcripts that were only detected from one or two single-cell samples, 108 (59.0%) had signal intensity >1,000. The strong and robust signal intensities detected from single-cell samples indicate that our system is very sensitive.
Unlike the gene transcripts detected from all non-singlecell samples which account for more than 95% of gene transcripts detected from each of these samples, tively, compared with = 0.949 for the non-single-cell samples. From all single and non-single-cell samples, 778 gene transcripts were detected, of which 315 (40.5%) were detected from some but not all samples. This is in contrast with the scenario of non-single-cell samples from which gene transcripts that were only detected from some but not all samples were a very small portion (Table 1). Furthermore, of these 315 transcripts, 177 (56.2%) were either detected from only single cells or from non-single cell samples.
The high degree of concordance among the results from the non-single-cell samples, and the significant differences among those from single cells, and between single cells and non-single-cell samples indicate that most, if not all, of these differences are real. As mentioned above, this is further supported by the robustness of the signal intensities detected from single-cell samples for the gene transcripts that were detected differently between the single cells and non-single-cell samples. It is conceivable that heterogeneity in clonality and/or genetic alterations in the cells of a cell line could be major factors contributing to the differences. In addition, a considerable portion of the cells may be at different cycle stages during which groups of genes are expressed differently. Therefore, while gene expression in single cells could differ in various aspects, 100 cells may well represent the entire cell population because, after all, the cell line cells are from the same tissue and the same donor. Therefore, genes that are detectable in a cell population may not be expressed or expressed at very low levels in certain single cells. Conversely, genes that are detectable in particular single cell samples may not be expressed or expressed at very low levels in the majority of the cell population.

Differential gene expression in the two cell lines, NCI/ ADR-RES and MCF-7
When the gene expression profiles of NCI/ADR-RES were compared with those of MCF-7, a considerable number of genes were shown to be expressed differentially in these two cell lines. Of the 1,135 gene products, 531 (46.8%) were detected from samples of both cell lines (not including single cell samples). Seventy-five gene transcripts were detected in all NCI/ADR-RES non-single-cell samples, but not in the MCF-7 samples, and 43 were detected in the opposite way.
Of the 118 differentially expressed genes, 69 were shown to be expressed with more than 10-fold difference ( Table  2). Of the 69 genes, 37 (53.6%) were detected as strongly or relatively strongly expressed in MCF-7, but weakly or not expressed in NCI/ADR-RES, and 32 were detected in the opposite way. To validate the gene expression data, 22 of these 69 genes, and another 46 gene transcripts detected with various microarray signal intensities differ-ent between the samples of the two cell lines were randomly selected (Table 3) and subjected to RT-PCR amplification individually. The amplified products were resolved by gel electrophoresis. The signal intensities of the respective bands were quantified with a gel documentation system. Part of results from microarrays and gel assays are shown in Fig. 3. Table 3 summarizes the results from both microarray and gel assays. Based on the results from microarray, genes in Table 3 are subdivided into four groups. Transcripts of Group I genes were detected from all samples, while no transcripts were detectable from all samples for Group IV. Transcripts of Group III genes were detected only from the NCI/ADR-RES samples but not from the MCF-7 samples, and those of Group II genes were detected in an opposite way. As shown, the signal ratios between NCI/ADR-RES and MCF-7 from microarray for the Group I genes are well in concordance with the ratios from gel assay, with a correlation coefficient of 0.679, indicating results from microarray and gel assays match very well. Because signals of genes in Groups II to IV were below the background for at least one of the two cell lines, ratio comparison may not be meaningful. However, it is clear that the microarray signals detected from MCF-7 are all greater than those from NCI/ADR-RES for the Group II genes, and vice versa for the Group III genes. Results from the gel assay are very well in concordance with this correlation. The only exception is HDAC5 whose microarray signal from MCF-7 is approximately 5 times that from NCI/ADR-RES, while its gel signal is more than 2 times that of the latter. Since the microarray signal intensity from MCF-7 for this gene is the lowest among the Group II genes, this discrepancy could be caused by wider variation of the low signal intensity. For all genes in Groups II to IV, if the microarray signals are lower than background, the corresponding gel signals are also low (<5,000) except for three genes, HDAC5 in Group II, and DAB2 and CTSZ in Group III. The fact that low or relatively low signals were detected by the gel assay for the genes whose array signals were weaker than background may be a reflection of the difference between the two assays. For microarray, all 1,135 transcripts were amplified in the same tube, while all transcripts analyzed by gel assay were amplified individually. It is known that during PCR, after the reaction reaches a saturation point, very little additional products may be generated. When the gene transcripts are amplified in a multiplex way, certain low-copy-number sequences may not be amplified to the detectable amounts when the reaction reaches a saturation point.

Specificity of the high-throughput gene expression system
The specificity of our high-throughput gene expression system was demonstrated by the results from different cell line samples and by those from different single cells as  described above. To further demonstrate the specificity of our system, human genomic DNA samples were amplified with the same multiplex RT-PCR procedure and analyzed by microarray. Very few probes (<0.2%) were shown to have signals above background (data not shown), indicating that our system is very specific and can discriminate between the target mRNA sequences from their genomic counterparts, and therefore, the unprocessed transcripts.
Our previous experience also showed that in the absence of specific templates, a few primer sets may amplify nonspecific sequences. However, such non-specific amplification may become undetectable in the presence of specific templates because the specific sequences are much stronger in competition. In addition, using specially designed probes also enhanced the specificity.

Discussion
Compared with other existing gene expression profiling methods, our approach has the following advantages: (1) Highly specific To date, no other high-throughput system has been reported to be highly discriminative of mRNA from other related DNA and RNA sequences. Using primers amplifying sequences across intron(s) and probes consisting of sequences in adjacent exons is a critical enhancement to achieve such high specificity. Furthermore, all primer, probe and amplicon sequences were subjected to exhaustive searches against the databases of the entire human genome and transcriptome to ensure these sequences are unique. Such a step was proven very effective for enhancing the specificity [12]. Experimentally, when genomic DNA was used as samples, signals were only detected for only 2 or 3 genes (0.2%) out of the 1,135 genes. Based on our previous studies, these signals may become undetectable in the presence of specific sequences which may compete out the nonspecific amplification.

(2) Highly sensitive
We showed previously that our multiplex amplification system could detect >1,000 single-copy sequences simultaneously from single haploid sperm cells [12]. The fact that >90% of these sequences are detectable indicates that with our specially designed primers, most, if not all, sequences may be well amplified parallelly with very limited, if any, interaction among the primers. Since the primers used for gene profiling are designed in the same way, it is reasonable to believe that most gene transcripts are also amplified parallelly. However, since the copy number of different gene transcripts in the cells varies in a wide range, the outcome of amplification would be differ-  Table 3. M, MCF-7, and A, NCI/ADR-RES.  ent from that using single-copy sequences. When only single-copy sequences are used in multiplex amplification, most, if not all, sequences may reach the detectable amount before the system is saturated. However, when gene transcripts are amplified, whether a transcript reaches a detectable amount before the system is saturated depends on its copy number in the sample, and not all sequences may reach a detectable amount at the end of amplification. This is probably why some sequences were undetectable by microarray but detectable by gel assay.

MCF-7 NCI/ADR-RES M A M A M A M A M A M A
With our system, a total of 686 gene transcripts were detected from three single cells, which is comparable to 676 for the three 100-cell samples and 693 for all non-single-cell samples from the same cell line. The sensitivity of our system is further proved by the facts that results from 100-cell samples are very similar to each other and to those from 10,000 cells, and that specific gene expression profiles were obtained from different cell lines using as few as 100 cells.
The sensitivity of our system is further illustrated by the results that a significant portion of transcripts that could not be detected from the NCI/ADR-RES samples but were detected from the MCF-7 samples or single cell samples, and vice versa. This also indicates that low microarray intensities for these transcripts were not false negatives, and they were either not present or present in very low abundance in the respective samples.

(3) Very simple
Unlike other methods that involve multiple steps and use multiple enzymes, our method allows a large number of gene products amplified by a single RT-PCR step directly from cell lysates without RNA extraction. In this way, a large number of samples may be analyzed easily and costeffectively. Our simple experimental procedure is also the basis of the high degree of sensitivity since it avoids complicated mRNA extraction and processing procedures before and during amplification, which may cause mRNA degradation or loss.

(4) Very safe for RNA samples
When working with RNA, one has to take extra precaution to prevent mRNA from degradation. Our method does not need RNA extraction. Once cells are lysed, RNA is directly released to the RT-PCR buffer and used as template immediately. There is almost no chance for RNase to degrade the mRNA templates.

(5) Highly flexible
Many studies may not need to analyze all genes in the human genome and may often need to focus on different gene groups. Therefore, flexibility of the experimental system would be highly desirable. With our computer pro-gram, a large number of gene products can be designed into a single multiplex group. Genes can be easily organized into different subgroups upon need, and can also be re-grouped at any time without altering the reaction conditions. New gene products can be added to an existing set easily.
The capacity of multiplex RT-PCR is another concern for high-throughput gene expression profiling because it not only makes the amplification of a large number of gene products affordable and cost-effective, but also eliminates challenges involved in quality control of RT-PCR for a large number of genes individually [33,34]. However, the capacity of multiplex amplification was limited by interaction between primers. A previous study reported a screening of 29 expressed genes using multiplex RT-PCR, but was unable to reduce the number of the reaction tubes less than eight [35]. Other studies multiplexed up to nine genes with nonspecific RT primers [36,37]. Studies using multiple sets of gene-specific primers in single reactions were also reported, but none of these generated enough products for the analysis of all expressed genes in the samples [34,38]. In the present study, we report our success with multiplex RT-PCR for 1,135 mRNA species. Such a success was based on a combination of several technological developments, including computerized primer design with predicted minimal interaction, a narrow primer T m range, small amplicon sizes, and optimization of amplification conditions based on our previous experience [12,29,30]. With our current protocol, it is possible to include two thousand or more gene transcripts in a single multiplex amplification group, and to analyze all human gene transcripts using several multiplexing amplification groups. After pooling amplified products from the multiplexing groups, all genes may be analyzed with a single microarray. With our system, large-scale gene expression profiling becomes highly affordable and cost-effective. If the primers and probes used in the high-throughput analysis are made accessible to the research community through a distribution system, large-and genome-scale gene expression profiling may be even more affordable and cost-effective.
A major limitation of our system is the requirement of presence of large introns in genes under study. When the introns are small, discrimination between mRNA and closely related DNA and RNA sequences is still possible by using probes consisting of sequences in the neighboring exons. For genes with no introns, primers and probes can be designed only to discriminate mRNA sequences from related pseudogenes and their transcripts but not the corresponding gene sequences. In this case, discrimination between mRNAs and their gene sequences is only possible when the mRNAs are present abundantly. High-throughput gene expression profiling with single cells is also interesting for most laboratories studying molecular neurophysiology, but has been hampered by the capacity of multiplex PCR. Our approach can be used to examine the expression of many genes within individual neurons or other cells. The gene expression profiles can also be correlated to the phenotypes of these cells such as morphological, electrophysiological and pharmacological features to understand the underlying molecular mechanisms.

Conclusion
This report describes a high-throughput gene expression profiling technology, which is simple, highly reproducible, specific and sensitive, and may greatly facilitate gene expression profiling of a small number of or even single cells. It may also be applicable to many applications where the amount of material is limited, and to diagnostic assays that identify the onset of cancer and monitor its progression, remediation or response to treatments.
Data discussed in this publication have been deposited in the NCBI's Gene Expression Omnibus [31] and are accessible through GEO Series accession number GSE5920.

Cell lines and single cell preparation
Human breast cancer cell line MCF-7 and ovarian cancer cell line NCI/ADR-RES were kindly provided by Drs. Jinming Yang, Hao Wu and William Hait [42]. The cell lines were maintained in RPMI 1640 medium containing 10% fetal bovine serum, 100 units/ml penicillin, and 100 μg/ ml streptomycin at 37°C in a humidified atmosphere containing 5% CO 2 . After counting with a hemacytometer, cells were suspended in PBS (phosphate buffer solution) to 1000 cells/μl or other desirable densities. Two μl was dispensed into an Eppendorf tube containing cell lysis buffer (1.5 μl RNase inhibitor, 4 μl of 5× QIAGEN OneStep RT-PCR buffer, 12.5 μl H 2 O). Single cells were prepared from a diluted cell suspension of 2 cells/μl in 1 × PBS. About 0.5 μl of the suspension was pipetted onto a small piece of glass coverslip, and was checked under a microscope. If the droplet contained only one cell, the piece of the coverslip was then transferred into an Eppendorf tube containing the cell lysis buffer. The tube was immediately frozen in an ethanol/dry ice bath and stored at -80°C until use.

Selection of genes for mRNA profiling
Genes used in the present study were selected based on previous publications [15][16][17][18][19][20][21][22][23][24][25][26][27][28], and are those involved in fundamental cell functions such as cell cycle, apoptosis, cell matrix, DNA repair, DNA replication, somatic recombination, RNA transcription and regulation, and protein translation and regulation. The borders between exons and introns for the selected genes were determined by aligning of the mRNA to genomic sequences using the BLAT program [43] maintained by the University of California, Santa Cruz.

Primer and probe design
A computer program was written for primer and probe selection. Each pair of PCR primers was designed to amplify sequences in two adjacent exons flanking a large intron and to ensure specific amplification of the desirable mRNA sequences rather than the respective gene or unprocessed RNA sequences. To enhance the amplification specificity, the program always searches for candidate amplicon sequences separated by large introns in each gene. The melting temperatures (T m 's) for all selected primers ranged from 50.1°C to 61.6°C, and the GC-contents ranged from 32% to 70%. The lengths of the amplicons ranged from 72 to 150 bases.
Each oligonucleotide probe for microarray analysis was designed to consist of sequences of two adjacent exons to specifically interrogate the cDNA from corresponding mRNA sequence, but not the corresponding gene sequences or cDNA from unprocessed RNA. To facilitate microarray analysis, the 3'-ends of all probes terminated before a "G" base in the template sequence so that they can be labeled with the same fluorescent color by incorporating fluorescently labeled Cy5-ddCTP. The lengths of the probes ranged from 22 to 31 bases, and the GC-content of the probes ranged from 30% to 70% with their T m 's from 54.4°C to 65.2°C. dom sequences, and was subtracted from the sample signals. For the comparative expression analysis of the cell lines MCF-7 and NCI/ADR-RES in Table 1, the array data were normalized by the Lowess smoothing method [45,46]. After background subtraction, genes with negative values of signal intensities in both duplicated samples were excluded for further analysis. The log ratios of the intensities of the remaining genes in two cells lines were used to make calls and to identify the differentially expressed genes in the samples.