- Methodology article
- Open Access
Increased DNA microarray hybridization specificity using sscDNA targets
BMC Genomicsvolume 6, Article number: 57 (2005)
The most widely used amplification method for microarray analysis of gene expression uses T7 RNA polymerase-driven in vitro transcription (IVT) to produce complementary RNA (cRNA) that can be hybridized to arrays. However, multiple rounds of amplification are required when assaying very small amounts of starting RNA. Moreover, certain cRNA-DNA mismatches are more stable than the analogous cDNA-DNA mismatches and this might increase non-specific hybridization. We sought to determine whether a recently developed linear isothermal amplification method (ribo-SPIA) that produces single stranded cDNA would offer advantages over traditional IVT-based methods for microarray-based analyses of transcript expression.
A single round of ribo-SPIA amplification produced sufficient sscDNA for hybridizations when as little as 5 ng of starting total RNA was used. Comparisons of probe set signal intensities obtained from replicate amplifications showed consistently high correlations (r = 0.99). We compared gene expression in two different human RNA samples using ribo-SPIA. Compared with one round IVT, ribo-SPIA had a larger dynamic range and correlated better with quantitative PCR results even though we used 1000-fold less starting RNA. The improved dynamic range was associated with decreases in hybridization to mismatch control probes.
The use of amplified sscDNA may offer substantial advantages over IVT-based amplification methods, especially when very limited amounts of starting RNA are available. The use of sscDNA targets instead of cRNA targets appears to improve hybridization specificity.
DNA microarrays are a powerful tool for global analysis of gene transcript expression. The initial studies using arrays required large amounts of starting material in order to reliably detect sample signals. Since that time, improvements in sample preparation, amplification and labeling methods [1–5] have reduced the starting material requirement to ~1–5 μg of total RNA . Efforts to use smaller amounts of starting material have focused on PCR [7, 8] and multiple rounds of T7 RNA polymerase in vitro transcription [IVT] [9–12]. PCR based methods have been successfully used to amplify subnanogram quantities of RNA from as little as a single cell [13, 14], but these approaches have not been widely adopted. Most attempts to perform arrays using submicrogram amounts of RNA have relied on 2 or 3 rounds of linear amplification using IVT, but this approach has proven to be time consuming and technically demanding. In our hands, two round IVT is necessary to prepare samples from 5–50 ng total RNA and the amplification typically takes 4–5 days to complete. Others have reported a 10% decrease in sensitivity in detection of differentially expressed genes with the addition of a second IVT round .
A new single primer, isothermal linear amplification method (ribo-SPIA) has been specifically developed for amplification of very small samples for use on DNA microarrays [16, 17]. With this method (Figure 1), small amounts of total RNA are reverse transcribed into cDNA using a chimeric RNA/DNA primer containing oligo(dT) and a unique RNA sequence tag at the 5' end. Linear amplification requires the addition of RNase H, DNA polymerase and excess chimeric primer. The RNase H digests RNA from RNA/DNA hybrids thus exposing a single stranded binding site where a new copy of the primer anneals and the DNA polymerase initiates synthesis of a fresh copy of cDNA, displacing the original antisense strand of the cDNA. A single isothermal linear amplification reaction rapidly generates sufficient single-stranded cDNA (sscDNA) for multiple hybridization reactions. sscDNA samples are fragmented to provide sscDNA fragments of ~50–200 bp, end labeled with biotin and used for microarray hybridization. Approximately 100,000-fold amplification is typical for a single amplification step. The ribo-SPIA method is potentially attractive because the amplification can be completed in a single day and there are no purification steps until after completion of amplification, thus reducing the risk of losing sample during handling. In this study, we have investigated the utility of ribo-SPIA generated sscDNA for DNA microarray analysis of small starting samples. We assessed yield and reproducibility of sscDNA, and compared sscDNA-based microarray results with microarray results obtained using IVT amplification and with results of quantitative PCR in order to assess if this method conferred new advantages for use with DNA microarrays.
Results and discussion
Yield of sscDNA
We used the ribo-SPIA method to amplify several different total RNA samples. RNA was obtained from a pool of human tissues (Clontech Universal Human Reference RNA, or cUHR, Experiment 1), mouse liver (Experiment 2), a second pool of human cells (Stratagene Universal Human Reference RNA, or sUHR, Experiment 3), and K562 human erythroleukemia cells (also Experiment 3). The amount of starting total RNA ranged from 5–100 ng and yields were in the range of ~6–12 μg of sscDNA (Table I). The somewhat lower yields seen in Experiment 1 are likely attributable to our unfamiliarity with the protocol, and yields improved in subsequent experiments. There was no clear relationship between the amount of starting RNA and the sscDNA yield. To determine how much sscDNA product was produced in the absence of any template, we performed two additional amplifications with no input RNA (Experiment 4). Some sscDNA was produced, although the amount was substantially less than that seen when input RNA was present (Table 1). To determine the possible impact of this template-independent product on microarray results, we hybridized the entire sscDNA product from a template-independent reaction to a U95Av2 microarray. This resulted in low overall signal intensity with only 0.6% of probe sets yielding "present" calls.
Size of sscDNA products
sscDNA preparations were analyzed by electrophoresis using an Agilent 2100 BioAnalyzer. sscDNAs ranged widely in size and the median size was typically slightly greater than 1 kb (data not shown). sscDNAs were fragmented in preparation for hybridization resulting in fragments of ~50–200 bp. These results are similar to those previously obtained using this method [16, 17].
Reproducibility of microarray hybridization results
Each of the experiments included replicate amplifications (independent amplifications of aliquots of the same starting material). We hybridized each replicate to a separate microarray and calculated intensities for each probe set. Fig. 2A shows an example of one pair of replicate hybridizations from Experiment 1, each performed using 20 ng of cUHR RNA as starting material. Pairwise comparisons of probe set intensities for replicate hybridizations (Table 2) produced very high correlations (r = 0.983–0.996) across the entire range of starting RNA amounts used in the three experiments (5–100 ng). One previous study also showed high correlations (r ~ 0.97–99) between replicate array data produced using the ribo-SPIA method . To assess how the amount of input RNA affects microarray results, we compared intensities found using 5 versus 100 ng of murine liver RNA (Experiment 2). When a 5 ng sample was compared to a 100 ng sample, there was a similar strong correlation (r = 0.987), indicating that large differences (20-fold) in starting material between reactions had small effects on measurements of gene expression (Fig. 2B). Previous reports have used correlation values between replicates to assess the reproducibility of other amplification methods. Those studies involve different laboratories and a wide range of microarray platforms, which makes direct comparison challenging. However, the correlations that we obtained using the ribo-SPIA method compare favorably with those reported for two rounds of IVT amplification (r = 0.92–0.98) [18–21], SMART amplification (r = 0.85–0.97) [22–24], and a PCR-based form of amplification (r = 0.97) .
Differential gene expression measurement
To test the ability to detect differential gene expression using amplified sscDNAs, we compared gene expression in two different RNA samples. We chose to compare K562 cell and sUHR RNAs since we have previously used these two RNAs to compare the performance of single round IVT-based amplification with other methods . In Experiment 3, we did three separate sscDNA amplifications of K562 and sUHR RNAs. We used the same RNA preparations as for the previously reported IVT-based amplification experiments, but started with 1000-fold less material (10 ng instead of 10 μg). For each probe set, we used intensity values from replicate arrays to calculate relative gene expression (M, defined as log2 [mean K562 intensity/mean sUHR intensity]) and average signal intensity (A, defined as 1/2 log2 K562 mean intensity + 1/2 log2 mean sUHR intensity). M and A values obtained using sscDNA are shown in Fig. 3B and those obtained using IVT-generated cRNA are shown in Fig. 3A. The range of signal intensities was similar for the two methods, although the mean intensity was lower for sscDNA hybridizations (A = 5.22 for sscDNA, A = 6.12 for cRNA). In contrast, the range of M values was somewhat larger with sscDNA. In particular, the sscDNA method identified several transcripts that were more than 27-fold higher in sUHR than K562 cell RNA (M < -7), but no differences of this magnitude were identified using cRNA. The number of probe sets associated with greater than 2-fold differences in expression (|M| > 1) was 1518 for sscDNA and 1043 for cRNA. 51% of the genes with >2-fold differences in expression on sscDNA arrays were not detected as >2-fold on cRNA arrays, but only 25% of the genes that were >2-fold different on cRNAs were not detected as >2-fold different on sscDNA arrays (Fig. 3).
The new observation that sscDNA gave a wider range of relative expression (M) values despite lower average intensity (A) values could be explained by improved hybridization specificity under the conditions used in this study. This is plausible because the binding energy for DNA-DNA interactions is more sensitive to base pair mismatching than the binding energy for DNA-RNA interactions [26, 27]. To look for further evidence about specificity of hybridization, we took advantage of the mismatched (MM) probes included on the arrays. For each perfect match (PM) GeneChip 25 mer probe, there is a corresponding MM probe with a single base mismatch at base 13. The MM probes were included in the probe set design to allow adjustments for nonspecific hybridization. Under ideal conditions, MM probes would never give signals higher than PM probes, although in practice this does sometimes occur. MM probes would be more likely to give stronger signals than PM probes if there was more non-specific hybridization of off-target sequences to the probes. We found that MM intensities exceeded PM intensities less frequently when we used sscDNA as compared to cRNA. When sUHR RNA was used as starting material, the average number of probe sets where MM intensity exceeded PM intensity was 2247 for cRNA versus 1671 for sscDNA (34% higher, p = 0.008). MM intensity also exceeded PM intensity more frequently with cRNA probes for K562 RNA arrays (2903 vs. 2482 probe sets, 17% higher, p = 0.017). When we looked at raw signal intensity for all MM probes, we found that the cRNA MM intensity distributions were skewed compared to the sscDNA MM distribution (Fig. 4). A closer examination of these distributions revealed that the use of sscDNA instead of cRNA resulted in a substantial reduction in the number of MM probes that gave relatively high intensity signals (Table 3). These findings strongly suggest that hybridization specificity is better for sscDNA than for cRNA. In a related study, Gingeras and coworkers  observed that increased nonspecific hybridization was observed when using directly labeled E. coli RNA as compared to cDNA. The increased nonspecificity was attributed to the presence of large amounts of rRNA in the samples. In our study however, both target preparations were prepared using oligo(dT) primers for the synthesis of first strand cDNA, so this explanation is less likely.
Comparison of expression measurements made with sscDNA, cRNA, and qPCR
We wished to compare how measurements made using amplified sscDNA and microarrays compared with measurements made using other approaches. We began by comparing results obtained using sscDNA and cRNA microarray hybridizations for all 12,625 probe sets on the arrays. Since the sscDNA and cRNA methods would be expected to introduce different systematic biases, we were not surprised that direct correlations between signal intensities obtained with the two different methods showed show relatively poor agreement (r = 0.72–0.75 for K562 and r = 0.68–0.70 for sUHR, as opposed to r = 0.98–0.99 between replicates performed using the same sample preparation method). The finding indicates that it will not be useful to directly compare one array hybridized with sscDNA to another one hybridized with cRNA.
We next compared differential gene expression measures (M values) determined using sscDNA with those determined using cRNA. There was a clear correlation (r = 0.83, Fig. 5A). We expected that probe sets associated with low intensity signals would give less reliable measures of gene expression and when we removed these probe sets from the analysis the correlation improved (r = 0.90, Fig. 5B). On average, the estimated M values were slightly larger (~1.2 times higher) when sscDNA was used instead of cRNA.
Next we generated another set of expression measurements that could be used as a basis for comparison for the sscDNA and cRNA array results. qPCR is typically used as "gold standard" to confirm putative differentially expressed genes detected with microarrays. Since we saw a subset of genes for which expression differed between sscDNA and cRNA targets, we next assessed if either method tracked more closely to qPCR. We chose qPCR primers and probe sets from a large group of >1000 sets that have been developed for various studies. From these, primers and probes for four subsets of genes were selected for qPCR. The first set included all genes with >4 fold difference in expression between K562 and sUHR samples as determined using the sscDNA method, the cRNA method, or both methods (53 primer/probe sets). The second set included all other genes in which the two methods disagreed by more than 2-fold (29 primer probe sets). The third set consisted of a group of 33 empirically-derived 'housekeeping genes.' These were all genes that were nearly equally expressed (|M| < 0.1) in K562 and sUHR samples according to both the sscDNA and cRNA methods and gave strong signals (A > 5 for both methods). The fourth set included 8 housekeeping genes that had been previously validated as controls for qPCR in other experiments. We determined the gene copy number for each qPCR primer and probe set and then calculated a measure of relative expression, M = log2 (K562 copy number)/(sUHR copy number), that could be directly compared to M values from arrays. 37 putative duplicate probe sets from 17 genes probe sets were hand-curated to confirm that they would correspond to the predicted qPCR product. In two cases probe sets were found to be misidentified in the GeneChip annotation and were removed from the analysis. In the remaining cases of duplication, the qPCR and microarray values were averaged across the duplicates. The final set of 106 curated genes and the associated data can be found at http://asthmagenomics.ucsf.edu.
There were clear correlations between qPCR M values and array M values obtained using sscDNA (Fig. 6A) or cRNA (Fig. 6B). When all 106 genes were included, qPCR results correlated slightly better with sscDNA (r = 0.72) than with cRNA (r = 0.66). When we included only the 29 genes for which sscDNA and cRNA methods disagreed by more than 2-fold, the difference between the two sample preparation methods became more pronounced (r = 0.75 for sscDNA vs. r = 0.57 for cRNA). Not surprisingly, both array-based methods tended to give smaller estimates of M than qPCR; this relative underestimation was somewhat less marked for the sscDNA than the cRNA method. In summary, results from arrays hybridized with sscDNA samples amplified using the ribo-SPIA method tracked with qRT-PCR more closely than did results from arrays hybridized with cRNA samples prepared using the traditional IVT method.
We examined the suitability of a new isothermal linear amplification method for application to Affymetrix GeneChip microarrays. We performed a series of tests using starting amounts of RNA ranging from 5 to 100 ng for amplification yield and reproducibility. The amplification reactions consistently produced sufficient sscDNA for multiple array hybridizations. Pairwise comparison of technical replicates hybridized to microarrays by regression analysis showed excellent consistency. When we used sscDNA to analyze differential gene expression between two samples, we found a larger dynamic range than that obtained with cRNA hybridizations. The improved performance appears to be related to increased sscDNA hybridization specificity. The data obtained using this new method also more closely matched the results from qRT-PCR than data obtained using standard IVT reactions, even though the amount of starting RNA used was 1000-fold less. This new amplification method is a useful alternative approach for preparing targets that is especially well-suited for experiments involving small amounts of starting material.
Clontech Human Universal Reference Pool total RNA (cUHR), derived by pooling RNA from a variety of human tissues, was purchased from BD Biosciences and used in Experiment 1. Mouse total liver RNA was isolated by standard methods from C57/BL6 mice according to procedures approved by the UCSF Committee on Animal Research and used in Experiment 2. For Experiment 3, we used Stratagene Human Universal Reference Pool and K562 erythroleukemia total RNAs from the same batches used in a previous study . All samples were assessed for size and integrity using the Agilent 2100 BioAnalyzer RNA 6000 Nano LabChip assay. RNA and DNA samples were quantified using a NanoDrop ND-1000 spectrophotometer.
sscDNA samples were prepared using the NuGEN Technologies Ovation RNA amplification and Biotin Labeling system (Version 1.0) according to the manufacturer's directions from the indicated amount of starting RNA (5–100 ng). All reactions were performed in 0.2 ml strip PCR tubes in an MJ GeneWorks PTC-100 thermocycler using recommended programs. Since the seal for PCR tubes and caps tends to deteriorate with repeated use, we replaced the caps for each tube before each resealing step in the protocol. Following amplification, sscDNA product was purified using QIAquick PCR purification kits (Qiagen). Samples were fragmented and end labeled with biotin. After stopping, each reaction was concentrated in a Microcon YM-3 column to a final volume of ~20 μl. The concentrated material was purified using a Centri-Sep 100 spin column (Princeton Separations). Negative control reactions were prepared by replacing input RNA with the appropriate volume of RNase free water.
All samples were placed in standard Affymetrix hybridization buffer. The sample denaturation time for the sscDNA samples was reduced from 5 to 2 minutes and hybridization time increased from 16 to 20 hours as recommended by NuGEN Technologies. cUHR gene expression was assayed using Affymetrix Human Genome U133A GeneChip arrays (Experiment 1). Mouse liver RNA was assayed using Murine Genome Mu6500A arrays (Experiment 2). K562 and sUHR RNAs were assayed using Human Genome U95Av2 arrays (Experiment 3). One template independent sample was also analyzed using a Human Genome U95Av2 array (Experiment 4). Arrays were stained with phycoerythrin-streptavidin according to the manufacturer's instructions. Metrics for all sample hybridizations including scaling factors, mean background intensities, and percent present calls have been provided (see Additional File 1). Each set of data was normalized independently using RMAExpress software http://stat-www.berkeley.edu/~bolstad/RMAExpress/RMAExpress.html. K562 and sUHR microarray data were also analyzed using Microarray Suite 4.0 in order to calculate PM-MM values for each transcript probe set. Probe level analyses were performed using the BioConductor  affy analysis package . All microarray data have been deposited in the Gene Expression Omnibus (GEO) database under the accession numbers GSM41384 – GSM41393, GSM41433 – GSM41438 and GSM4843 – GSM4847.
Real-time (RT) PCR was used to measure the expression of selected genes in sUHR and K562 cells. Gene-specific primers for multiplex real time RT-PCR were designed for each gene of interest using "Primer Express" software (Perkin-Elmer) and purchased from Biosearch Technologies. Sequence data for all oligonucleotides primers has been provided (see Additional File 2). First strand cDNA synthesis was performed using total RNA, Powerscript reverse transcriptase (BD Biosciences), and random hexamer primers. Real time amplification was performed using an ABI Prizm7900 and Invitrogen Universal Master Mix. Relative gene copy numbers (GCN) were calculated as described previously . GeNorm  was used to select the two most stable housekeeping genes across all specimens for normalization.
Luo L, Salunga RC, Guo H, Bittner A, Joy KC, Galindo JE, Xiao H, Rogers KE, Wan JS, Jackson MR, Erlander MG: Gene expression profiles of laser-captured adjacent neuronal subtypes. Nat Med. 1999, 5: 117–122-10.1038/4806.
Ohyama H, Zhang X, Kohno Y, Alevizos I, Posner M, Wong DT, Todd R: Laser capture microdissection-generated target sample for high-density oligonucleotide array hybridization. Biotechniques. 2000, 29: 530–536-
Luzzi V, Holtschlag V, Watson MA: Expression profiling of ductal carcinoma in situ by laser capture microdissection and high-density oligonucleotide arrays. Am J Pathol. 2001, 158: 2005–2010-
Dunmire V, Wu C, Symmans WF, Zhang W: Increased yield of total RNA from fine-needle aspirates for use in expression microarray analysis. Biotechniques. 2002, 33: 890–896-
Luzzi V, Mahadevappa M, Raja R, Warrington JA, Watson MA: Accurate and reproducible gene expression profiles from laser capture microdissection, transcript amplification, and high density oligonucleotide microarray analysis. J Mol Diagn. 2003, 5: 9–14-
Eberwine J, Kacharmina JE, Andrews C, Miyashiro K, McIntosh T, Becker K, Barrett T, Hinkle D, Dent G, Marciano P: mRNA expression analysis of tissue sections and single cells. J Neurosci. 2001, 21: 8310–8314-
Livesey FJ, Furukawa T, Steffen MA, Church GM, Cepko CL: Microarray analysis of the transcriptional network controlled by the photoreceptor homeobox gene Crx. Curr Biol. 2000, 10: 301–310-10.1016/S0960-9822(00)00379-1.
Iscove NN, Barbara M, Gu M, Gibson M, Modi C, Winegarden N: Representation is faithfully preserved in global cDNA amplified exponentially from sub-picogram quantities of mRNA. Nat Biotechnol. 2002, 20: 940–943-10.1038/nbt729.
Eberwine J, Yeh H, Miyashiro K, Cao Y, Nair S, Finnell R, Zettel M, Coleman P: Analysis of gene expression in single live neurons. Proc Natl Acad Sci U S A. 1992, 89: 3010–3014-
Baugh LR, Hill AA, Brown EL, Hunter CP: Quantitative analysis of mRNA amplification by in vitro transcription. Nucleic Acids Res. 2001, 29: e29-10.1093/nar/29.5.e29.
Kamme F, Salunga R, Yu J, Tran DT, Zhu J, Luo L, Bittner A, Guo HQ, Miller N, Wan J, Erlander M: Single-cell microarray analysis in hippocampus CA1: Demonstration and validation of cellular heterogeneity. J Neurosci. 2003, 23: 3607–3615-
Scherer A, Krause A, Walker JR, Korn A, Niese D, Raulf F: Early prognosis of the development of renal chronic allograft rejection by gene expression profiling of human protocol biopsies. Transplantation. 2003, 75: 1323–1330-
Chiang MK, Melton DA: Single-cell transcript analysis of pancreas development. Dev Cell. 2003, 4: 383–393-10.1016/S1534-5807(03)00035-2.
Tietjen I, Rihel JM, Cao Y, Koentges G, Zakhary L, Dulac C: Single-cell transcriptional analysis of neuronal progenitors. Neuron. 2003, 38: 161–175-10.1016/S0896-6273(03)00229-0.
Gold D, Coombes K, Medhane D, Ramaswamy A, Ju Z, Strong L, Koo JS, Kapoor M: A comparative analysis of data generated using two different target preparation methods for hybridization to high-density oligonucleotide microarrays. BMC Genomics. 2004, 5: 2 (http://www.biomedcentral.com/1471-2164/5/2). 10.1186/1471-2164-5-2.
Dafforn A, Chen P, Deng G, Herrler M, Iglehart D, Koritala S, Lato S, Pillarisetty S, Purohit R, Wang M, Wang S, Kurn N: Linear mRNA amplification from as little as 5 ng total RNA for global gene expression analysis. Biotechniques. 2004, 37: 854–857-
Singh R, Maganti RJ, Jabba SV, Wang M, Deng G, Heath JD, Kurn N, Wangemann P: Microarray based comparison of three amplification methods for nanogram amounts of total RNA. Am J Physiol Cell Physiol. (in press)-
Wilson CL, Pepper SD, Hey Y, Miller CJ: Amplification protocols introduce systematic but reproducible errors into gene expression studies. Biotechniques. 2004, 36: 498–506-
Klur S, Toy K, Williams MP, Certa U: Evaluation of procedures for amplification of small-size samples for hybridization on microarrays. Genomics. 2004, 83: 508–517-10.1016/j.ygeno.2003.09.005.
Kenzelmann M, Klären R, Hergenhahn M, Bonrouhi M, Gröne HJ, Schmid W, Schütz G: High-accuracy amplification of nanogram total RNA amounts for gene profiling. Genomics. 2003, 83: 550–558-
Feldman AL, Costouros NG, Wang E, Qian M, Marincola FM, Alexander HR, Libutti SK: Advantages of mRNA amplification for microarray analysis. Biotechniques. 2002, 33: 906–914-
Puskás LG, Zvara , Hackler LJ, Van Hummelen P: RNA amplification results in reproducible microarray data with slight ratio bias. Biotechniques. 2002, 32: 1330–1340-
Wang J, Hu L, Hamilton SR, Coombes KR, Zhang W: RNA amplification strategies for cDNA microarray experiments. Biotechniques. 2003, 34: 394-400.
Saghizadeh M, Brown DJ, Tajbakhsh J, Chen Z, Kenney MC, Farber DB, Nelson SF: Evaluation of techniques using amplified nucleic acid probes for gene expression profiling. Biomol Eng. 2003, 20: 97–106-
Barczak A, Rodriguez MW, Hanspers K, Koth LL, Tai YC, Bolstad BM, Speed TP, Erle DJ: Spotted long oligonucleotide arrays for human gene expression analysis. Genome Res. 2003, 13: 1775–1785-10.1101/gr.1048803.
Casey J, Davidson N: Rates of formation and thermal stabilities of RNA:DNA and DNA:DNA duplexes at high concentrations of formamide. Nucleic Acids Res. 1977, 4: 1539–1552-
Sugimoto N, Nakano M, Nakano SI: Thermodynamics—structure relationship of single mismatches in RNA/DNA duplexes. Biochemistry. 2000, 39: 11270-11281. 10.1021/bi000819p.
Rosenow C, Saxena RM, Durst M, Gingeras TR: Prokaryotic RNA preparation methods useful for high density array analysis: Comparison of two approaches. Nucleic Acids Res. 2001, 29: e112-10.1093/nar/29.22.e112.
Gentleman RC, Carey VJ, Bates DM, Bolstad B, Dettling M, Dudoit S, Ellis B, Gautier L, Ge Y, Gentry J, Hornik K, Hothorn T, Huber W, Iacus S, Irizarry R, Leisch F, Li C, Maechler M, Rossini AJ, Sawitzki G, Smith C, Smyth G, Tierney L, Yang JYH, Zhang J: Bioconductor: Open software development for computational biology and bioinformatics. Genome Biol. 2004, 5: R80-10.1186/gb-2004-5-10-r80.
Gautier L, Cope L, Bolstad BM, Irizarry RA: Affy—analysis of Affymetrix GeneChip data at the probe level. Bioinformatics. 2004, 20: 306–315-10.1093/bioinformatics/btg405.
Dolganov GM, Woodruff PG, Novikov AA, Zhang Y, Ferrando RE, Szubin R, Fahy JV: A novel method of gene transcript profiling in airway biopsy homogenates reveals increased expression of a Na+-K+-Cl– cotransporter (NKCC1) in asthmatic subjects. Genome Res. 2001, 11: 1473–1483-10.1101/gr.191301.
Vandesompele J, De Preter K, Pattyn F, Poppe B, Van Roy N, De Paepe A, Speleman F: Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes. Genome Biol. 2002, 3: research0034.1–0034.11. Available at: http://genomebiology.com/2002/3/7/research/0034.-10.1186/gb-2002-3-7-research0034.
This study was supported in part in the General Clinical Research Center at San Francisco General Hospital and supported by Grant 5-MO1-RR00083 from the Division of Research Resources, National Institutes of Health, the UCSF NHLBI Shared Microarray Facility (NIH grant HL072301) and the UCSF Sandler Center for Basic Research in Asthma. We thank Yanxia Hao, Jennifer Gregg, Agnes Paquet and Michael Salazar for excellent technical assistance, Ron Davis for helpful discussions and Karena Essex for assistance in preparation of the manuscript.
CSB conceived of the study, participated in the design and coordination, and wrote the manuscript. CG participated in the design of the study and both CG and KH performed the microarray experiments. GMD carried out the real time PCR studies. JYHY performed the statistical analyses. DJE participated in the study design and analysis and revised the manuscript. All authors read and approved the final manuscript.
Electronic supplementary material
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.