Systematic analysis of T7 RNA polymerase basedin vitrolinear RNA amplification for use in microarray experiments

  • Jörg Schneider1,

    Affiliated with

    • Andreas Buneß1,

      Affiliated with

      • Wolfgang Huber1,

        Affiliated with

        • Joachim Volz2,

          Affiliated with

          • Petra Kioschis3,

            Affiliated with

            • Mathias Hafner3,

              Affiliated with

              • Annemarie Poustka1 and

                Affiliated with

                • Holger Sültmann1Email author

                  Affiliated with

                  BMC Genomics20045:29

                  DOI: 10.1186/1471-2164-5-29

                  Received: 30 October 2003

                  Accepted: 30 April 2004

                  Published: 30 April 2004

                  Abstract

                  Background

                  The requirement of a large amount of high-quality RNA is a major limiting factor for microarray experiments using biopsies. An average microarray experiment requires 10–100 μg of RNA. However, due to their small size, most biopsies do not yield this amount. Several different approaches for RNA amplificationin vitrohave been described and applied for microarray studies. In most of these, systematic analyses of the potential bias introduced by the enzymatic modifications are lacking.

                  Results

                  We examined the sources of error introduced by the T7 RNA polymerase based RNA amplification method through hybridisation studies on microarrays and performed statistical analysis of the parameters that need to be evaluated prior to routine laboratory use. The results demonstrate that amplification of the RNA has no systematic influence on the outcome of the microarray experiment. Although variations in differential expression between amplified and total RNA hybridisations can be observed, RNA amplification is reproducible, and there is no evidence that it introduces a large systematic bias.

                  Conclusions

                  Our results underline the utility of the T7 based RNA amplification for use in microarray experiments provided that all samples under study are equally treated.

                  Background

                  The utility of microarrays for disease classification, prognosis and progression, or the identification of target genes for novel therapeutic approaches is well documented [15] and is likely to change our views of disease development [6,7]. The major bottlenecks for these experiments are the limited availability and low quality of tissue or RNA. The issue of quality can be easily solved using appropriate tissue handling techniques. However, solid tumours are usually too small to yield enough RNA for direct use in microarray experiments. Therefore, amplification techniques have to be applied to increase the amount of available RNA and minimise the required starting material. Two basically different approaches for RNA amplification have been used by various laboratories: In the first, the polymerase chain reaction (PCR) is performed to increase the amount of sample either exponentially [8,9] or linearly [10]. The second approach appliesin vitrotranscription with T7 RNA polymerase [1114] for linear amplification of the sample. Both approaches differ in the lengths of the double-stranded cDNA molecules generated in the reverse transcription process prior to amplification. Since enzymatic modifications on a highly complex mixture of RNA/cDNA molecules are performed during the reaction, it can be presumed that noise and systematic biases are introduced. Thus, it is essential to quantify the effect of amplification and examine the reproducibility of the method.

                  To study genetic changes during breast cancer, we have developed a cDNA microarray consisting of 7,347 genes and expressed sequence tags (ESTs). We have established a robust method to amplify total RNA via a T7 RNA polymerase based process [11] and used this method for the amplification of RNA from different tissues and hybridised the corresponding labeled cDNA products on cDNA microarrays. Our goals were to assess the reproducibility of the amplifications, compare the results between unamplified and amplified RNA (aRNA), and quantify the bias introduced to the data by the RNA amplification process.

                  Results

                  Amplification factors

                  We used the linear amplification technique of van Gelder et al. [11] which is based on a double-stranded cDNA synthesis with an oligo-dT primer coupled to the T7 RNA polymerase promoter and subsequent in vitro transcription into aRNA using T7 RNA polymerase. In 30 separate reactions, we achieved amplification factors of 150–560 when 100–3,200 ng of total RNA was used as starting material (table1). The amplification of no more than 100–200 ng of total RNA (14 reactions) resulted in a 368 fold average increase of poly (A)+RNA equivalents whereas the amplification of 1,000–3,200 ng (16 reactions) resulted in an average amplification factor of 253. All amplification factors were calculated according to the assumption that 5% of total RNA correspond to poly (A)+RNA. The quality of all total and amplified RNA was checked using the Agilent 2100 Bioanalyzer. Only high-quality RNA was used for amplification because low-quality RNA can influence the outcome of amplification. The measurement of the amplified RNA showed a length reduction of the first-round amplified RNA compared to the total RNA (not shown). The sizes of the amplified RNA molecules were distributed between 100 and 4,500 bases. In total, 69 microarray hybridisations were used for the analysis of the different parameters and were performed as summarised in tables2and3. All microarrays used in these experiments were taken from two spotting runs in which a total of 200 arrays had been produced.
                  Table 1

                  Amplification factors

                  tot RNA starting amount [ng]

                  Repetitions

                  Range of amplified RNA [μg]

                  Average amount of amplified RNA [μg]

                  Range of amplification factors

                  Average amplification factor*

                  100–200

                  14

                  1–4.4

                  2.49

                  200–560

                  368

                  1,000–3,200

                  16

                  12–36

                  22

                  150–500

                  253

                  * calculated under the assumption that 5% of total RNA corresponds to poly (A)+RNA

                  Table 2

                  Microarray hybridisations used for the data analysis

                  Parameter

                  Figure

                  Microarrays used

                  Total Nr. of arrays

                  Normalisation group

                  Comparability

                  1a

                  Ta: 1,2

                  10

                  α

                    

                  Tb: 3,4

                    
                    

                  Tc: 5–7

                    
                    

                  Td: 8–10

                    
                    

                  Aa: 11,12

                  10

                   
                    

                  Ab: 13,14

                    
                    

                  Ac: 15–17

                    
                    

                  Ad: 18–20

                    
                   

                  1b

                  A: 21–35

                  29

                  β

                    

                  B: 36–43

                    
                    

                  C: 44–49

                    

                  Linearity

                  2a

                  5–7,50–53 (Cy5)

                  14

                  γ

                    

                  15–17,54–57 (Cy5)

                    
                   

                  -

                  5–10 (Cy3)

                  12

                  δ

                    

                  15–20 (Cy3)

                    

                  Correlation

                  -

                  44–49

                  6

                  ε

                  Reproducibility

                  3

                  58–69

                  12

                  ζ

                  Correspondence

                  4

                  50–57

                  8

                  η

                  Table 3

                  Experimental design of microarray hybridizations

                  Array Nr.

                  Cy3 RNA source*

                  Cy5 RNA source*

                  Array Nr.

                  Cy3 RNA source*

                  Cy5 RNA source*

                  Array Nr.

                  Cy3 RNA source*

                  Cy5 RNA source*

                  1

                  b

                  a

                  24

                  a

                  A2

                  47

                  B1

                  B2

                  2

                  b

                  a

                  25

                  b

                  B1

                  48

                  B1

                  B3

                  3

                  c

                  b

                  26

                  b

                  B1

                  49

                  B1

                  B4

                  4

                  c

                  b

                  27

                  b

                  B2

                  50

                  d

                  pla

                  5

                  ref2

                  pla

                  28

                  b

                  B2

                  51

                  d

                  pla

                  6

                  ref2

                  pla

                  29

                  a

                  A3

                  52

                  d

                  pla

                  7

                  ref2

                  pla

                  30

                  a

                  A3

                  53

                  d

                  pla

                  8

                  ref2

                  d

                  31

                  a

                  A4

                  54

                  D

                  PLA

                  9

                  ref2

                  d

                  32

                  b

                  B3

                  55

                  D

                  PLA

                  10

                  ref2

                  d

                  33

                  b

                  B3

                  56

                  D

                  PLA

                  11

                  B1

                  A1

                  34

                  b

                  B4

                  57

                  D

                  PLA

                  12

                  B1

                  A1

                  35

                  b

                  B4

                  58

                  REF1

                  A1

                  13

                  C4

                  B4

                  36

                  A1

                  A1

                  59

                  REF1

                  A2

                  14

                  C4

                  B4

                  37

                  A2

                  A2

                  60

                  REF1

                  A3

                  15

                  REF2

                  PLA

                  38

                  A3

                  A3

                  61

                  REF1

                  A4

                  16

                  REF2

                  PLA

                  39

                  A4

                  A4

                  62

                  REF1

                  A1

                  17

                  REF2

                  PLA

                  40

                  B1

                  B1

                  63

                  REF1

                  A1

                  18

                  REF2

                  D

                  41

                  B2

                  B2

                  64

                  REF1

                  A2

                  19

                  REF2

                  D

                  42

                  B3

                  B3

                  65

                  REF1

                  A2

                  20

                  REF2

                  D

                  43

                  B4

                  B4

                  66

                  REF1

                  A3

                  21

                  a

                  A1

                  44

                  A1

                  A2

                  67

                  REF1

                  A3

                  22

                  a

                  A1

                  45

                  A1

                  A3

                  68

                  REF1

                  A4

                  23

                  a

                  A2

                  46

                  A1

                  A4

                  69

                  REF1

                  A4

                  * lowercase letters: total RNA; capital letters: amplified RNA. A, a, B, b, C, c: RNA from human (untreated) mcf7 cell tumours induced in different SCID mice; D, d: adriamycin-treated mcf7 cell tumours grown in a SCID mouse. PLA, pla: human placenta; REF1, ref1, REF2, ref2: different pools of a, b and c.

                  Comparability of data sets obtained from amplified and non-amplified RNA

                  In 49 hybridisation experiments (tables2and3), we compared the generalised log ratios of 7,347 genes using amplified or total RNA for labelling (figure1). When labelled samples derived from total (Cy-3) vs. total (Cy-5) RNA were hybridised on 10 different microarrays, the boxplots, indicating the variation of the generalised log ratios, were generally wider in comparison to those of a second set of 10 microarrays with samples derived from amplified (Cy-3) vs. amplified (Cy-5) RNA (figure1A). This result holds for the cases where differential expression is not expected (figure1A, Ta and Tb compared to Aa and Ab) due to the use of the same tissue source for RNA isolation or whether it is anticipated (figure1A, Tc and Td compared to Ac and Ad) since different tissues were used. These findings indicate a considerably lower variability of generalised log ratios of aRNA vs. aRNA hybridisations compared to total vs. total RNA hybridisations.
                  http://static-content.springer.com/image/art%3A10.1186%2F1471-2164-5-29/MediaObjects/12864_2003_Article_131_Fig1_HTML.jpg
                  Figure 1

                  Comparison of generalised log ratios in hybridisations on 7,347-gene microarraysEach boxplot (x-axis) characterises the distribution of the generalised log ratios (y-axis) ofA)total RNA vs. total RNA and aRNA vs. aRNA from the same origin (duplicate hybridisations of Ta, Tb and Aa, Ab, respectively) and from different tissue sources (triplicate hybridisations of Tc, Td and Ac, Ad, respectively) andB)15 microarray hybridisations from the same tissue sources using total RNA vs. aRNA (A) and 14 hybridisations using aRNA vs. aRNA (B and C).

                  In a second experiment (figure1B), we hybridised total RNA against aRNA from the same tissue samples on 15 microarrays. The result (figure1B, group A) shows that the degree of variation of generalised log ratios is even larger than in the case of the total vs. total RNA hybridisation (figure1A). Although no differentially expressed gene is expected in these experiments, the high variation of the generalised log ratios indicates decreased sensitivity of the measurement if total RNA is compared to aRNA on the same array. Therefore, different types of RNA (amplified and total RNA) should not be compared within one experiment. In contrast, only small variations were observed for the generalised log ratios of the corresponding 14 microarray hybridisations with aRNA vs. aRNA (figure1B, groups B and C). In summary, the variability of the generalised log ratios of the 49 microarray hybridisations decreased from (1) total vs. aRNA, (2) total vs. total RNA to (3) aRNA vs. aRNA.

                  Linearity of amplification reactions

                  Prior to performing large scale microarray gene expression studies using RNA amplification, it is important to address the question as to whether the enzymatic modification of the RNA introduces a systematic bias to the list of differentially expressed genes. In other words: Is the amplification "linear"? To approach this problem, we hybridised 20 microarrays with samples derived from either total vs. total RNA or aRNA vs. aRNA from the same sources (table2). The arithmetic mean of the resulting values from all hybridisations was determined for each gene and each channel separately (figure2, Cy-3 channel not shown). In the ideal situation (100% correlation, identity line), the overall correlation coefficient should be equal to 1. The scattering of values around the identity line reflects the influences of gene specific RNA amplification factors and noise. For both channels, the scattering of values is homogeneously distributed on both sides of the identity line, which indicates an overall linear amplification. The correlation coefficients between total and aRNA calculated from the experiments were 0.87 and 0.89 for the individual Cy-5 and Cy-3 channels, respectively. Thus, there is a clear correlation between signal intensities resulting from non-amplified RNA as compared to amplified RNA, suggesting that aRNA hybridisations correspond to those obtained with unmodified samples.
                  http://static-content.springer.com/image/art%3A10.1186%2F1471-2164-5-29/MediaObjects/12864_2003_Article_131_Fig2_HTML.jpg
                  Figure 2

                  Correlation of normalised signal intensitiesAveragedvsn-normalised signal intensities [quanta] of unamplified samples (x-axis) were compared to amplified ones (y-axis) for the Cy-5 channel. Fourteen microarray hybridisations were used (for details, see tables2and3). A symmetric scatter around the identity line indicates an overall linear amplification. Similar results were obtained for the Cy-3 channel.

                  Reproducibility of RNA amplification

                  Several other groups [8,10,14,15] have applied Pearson correlation coefficients between log ratios in order to show the reproducibility of the RNA amplifications. High correlation coefficients are 0.8 or greater. We performed pairwise calculations of correlation coefficients using six hybridisations with aRNA samples (table2, correlation). In result, the coefficients varied from 0.95 to 0.99 indicating a very high correlation between the individual amplification reactions.

                  A statistically solid way to estimate the reproducibility of hybridisations is theFstatistic in which the ratio of variancesbetweenthe different amplification groups is compared to thatwithinthe amplification groups. When thebetweenvariation is almost the same as thewithinvariation, the logarithmicF-value comes to lie around zero, indicating high reproducibility; when the variationbetweenis larger or smaller thanwithinhybridisations, theF-value becomes larger or smaller than zero, respectively. To estimate the ratios, we performed four amplification reactions from the same total RNA source. The labelled products were hybridised in triplicate against an amplified reference sample on 12 microarrays (table2, reproducibility). These experiments were carried out on two different days (table2; day 1, arrays 58–61; day 2, arrays 62–69). For each gene from these hybridisations, anF-value was calculated and plotted in a histogram (figure3). The negative median logF-value (-0.66) indicated a higher variancewithinthanbetweenthe amplification groups. When only hybridisations performed on day 2 were considered, a median logFof 0.20 was found. These results suggest that other sources of variation (e.g. day-to-day variation, labelling, chip-batches) are of greater influence on the total variation of microarray hybridisations rather than the RNA amplification itself.
                  http://static-content.springer.com/image/art%3A10.1186%2F1471-2164-5-29/MediaObjects/12864_2003_Article_131_Fig3_HTML.jpg
                  Figure 3

                  Reproducibility of RNA amplificationsThe distribution (frequency, y-axis) is plotted against the gene specific logFvalues (x-axis). aRNA samples from four distinct amplification reactions were used in 12 separate microarray hybridisations.

                  Correspondence of gene expression data obtained with amplified and total RNA

                  For the biological interpretation of the data, the crucial question concerning the use of RNA amplification methods is: Which proportion of the differentially expressed genes identified by hybridisation of total vs. total RNA corresponds to those found by hybridisation of aRNA vs. aRNA? To examine this question, we used eight microarray hybridisations (table2), four of which were performed using total RNA vs. total RNA from different tissues, the other four with aRNA vs. aRNA derived from these total RNA samples. The generalised log ratios of four hybridisations were averaged and plotted against each other. For the highest 2,000 out of 7,347 genes, the correlation coefficient was 0.85 (figure4). When low intensity genes were included, the correlation dropped to 0.72 (not shown), owing to the lower signal to noise ratios for these genes. The distribution of generalised log ratios in the aRNA hybridisations appears to be compressed (slope of 0.54 for the regression line), indicating a decreased sensitivity of the aRNA results compared to total RNA hybridisation data. Within the group of 2,000 highest expressed genes we also looked in each case, total RNA and aRNA experiments, at the 500 highest differentially expressed genes. The minimum fold change was 2.6 for total RNA and 1.8 for amplified RNA experiments. Out of 344 common genes (69%) all showed the same direction of regulation.
                  http://static-content.springer.com/image/art%3A10.1186%2F1471-2164-5-29/MediaObjects/12864_2003_Article_131_Fig4_HTML.jpg
                  Figure 4

                  Correspondence of gene expression data between amplified and non-amplified samplesThe averaged generalised log ratios of the total RNA hybridisations (x-axis) are plotted against those of the aRNA hybridisations (y-axis). Only the 2,000 highest intensity genes were considered.

                  Discussion

                  In vitrotranscription of cDNA for the amplification of RNA is a commonly used method that requires low amounts of input total RNA, typically less than 1 μg, often as little as 100 ng. This amount can be reliably reduced further to the level of RNA from a few hundred cells when microdissected tissue is used for RNA isolation. Thus, the method is an absolute prerequisite for gene expression microarray experiments in which the amount of sample is the major limiting factor. We have obtained average amplification factors of more than 350 poly (A)+equivalents in a single amplification round. With these values, the amount of aRNA derived from 1 μg total RNA is sufficient to perform at least eight microarray hybridisations with 2 μg aRNA each. The benefit is twofold: First, it is possible to perform multiple experiments with minute amounts of sample and perform a statistical analysis of the data. Second, additional data validation experiments (e.g. quantitative RT-PCR) can be performed using the same physical sample. These benefits will further increase in the future when array based diagnostics become available.

                  Despite the obvious advantages of the amplification method, systematic analyses of the impact of this enzymatic modification on the microarray gene expression results are scarce. From the analysis of the variability of the results, we conclude that it is vital to use equally treated samples for any particular experiment: If one sample requires amplification, all other samples should be amplified as well. Comparison of data from amplified and non-amplified RNA results in the widest log ratio distribution of all combinations of aRNA and total RNA samples. This may lead to the identification of too many false positive differentially expressed genes. The opposite effect was seen when amplified vs. amplified RNA was used: The log ratio plots are slimmer when compared to total vs. total RNA hybridisations. Although this leads to a loss of sensitivity of the experiment, the criteria for the identification of differentially expressed genes are more stringent when aRNA rather than total RNA samples are used. RNA amplification does not introduce a systematic bias to the gene expression data: Microarray hybridisation results obtained with aRNA are comparable with those from total RNA. The relationship is linear over the whole dynamic range. Despite the linearity for most genes, however, one should consider the possibility that, due to sequence or structure-specific properties, a minority of genes have gene-specific amplification factors different to the rest. We cannot exclude the existence of such genes; however, we are not aware of a systematic study that has addressed this question down to the level of individual genes.

                  RNA amplification reactions are very reproducible. We found high correlation coefficients between repeated microarray hybridisations indicating that there is no major bias due to the amplification reaction that affects the results obtained from separate amplification reactions. Thus, in the light of all other well known technologically caused sources of error in microarray experiments which have to be considered carefully when the experiments are designed, the error introduced by RNA amplification is small. The differentially expressed genes identified by hybridisation of aRNA vs. aRNA correspond well with those found by hybridisation of total RNA. For the differentially expressed genes with the 27% highest signal intensities, the correlation coefficient was 0.85.

                  Conclusions

                  We showed that the use of the T7 RNA polymerase based amplification prior to microarray experiments is well-grounded because no systematic bias is introduced. Although variations in differential expression between amplified and total RNA hybridisations can be observed, the amplifications are reproducible, and the results of aRNA correspond to those obtained with total RNA. The slightly lower sensitivity of aRNA to detect differential gene expression can be compensated by the use of additional technical improvements or biological repetitions coupled to statistical data analysis.

                  Methods

                  Tissue samples and RNA isolation

                  Tissue samples: Intraperitoneal injection of the human breast cancer cell line mcf7 into severe combined immunodeficiency (SCID) mice resulted in the growth of a solid tumour mass [1618] varying in diameter from 1 to 7 mm. Six to nine weeks post injection, mice were sacrificed by application of CO2, and the tumour mass was immediately frozen in liquid nitrogen. Human placental tissue was frozen as soon as possible. Total cellular RNA was isolated with the Trizol-method (TriFast, peqlab, Erlangen, Germany), following the manufacturers instructions, after homogenisation with a Mikro-Dismembrator S (Braun Biotech, Melsungen, Germany). The quality of the RNA samples was checked by the Agilent 2100 bioanalyzer (Agilent Technologies, Waldbronn, Germany). Only high-quality RNA samples (ratio 28S/18S rRNA > 1.8) were selected for the experiments.

                  RNA amplification

                  The cDNA synthesis was performed using a cDNA Synthesis System Kit (Roche Diagnostics, Mannheim, Germany) according to the manufacturers instructions. The ds cDNA was purified by phenol/chloroform extraction [19] precipitated in 5 M NH4OAc [19] and resuspended in 8 μl RNase-free water.

                  In vitro transcription [11] was performed using the AmpliScribeTM T7 High Yield Transcription Kit (Epicentre Technologies, Madison, WI, USA) according to the manufacturers instructions with the exception of increasing the incubation time to 4 hours. The newly synthesized aRNA was extracted in phenol/chloroform/isoamylalcohol as described above. Unincorporated dNTPs were removed by chromatography with MobiSpin S-300 columns (MoBiTec GmbH, Göttingen, Germany). Ethanol precipitation was performed as described in [19] and the resulting pellet was resuspended in 20–50 μl RNase-free water.

                  Microarray spotting

                  Glass slides used for this study carried 7,347 breast specific cDNA clones selected from the Human UniGene 1 clone set (German Resource Centre for Genome Research, Berlin, Germany). The clones were amplified by PCR and spotted in a "replicate array" design on silane-prep™ slides (Sigma Diagnostics, St. Louis, USA) to produce 14,694 individual spots. The slides were rehydrated, and the DNA was denatured with boiling water treatment prior to washing with 0.2% SDS, MilliQ- H2O, 95% and 100% Ethanol. After the washing procedure, the microarrays were dried with compressed air.

                  RNA labelling and hybridisation

                  2 μg amplified antisense RNA (aRNA) were mixed with 500 ng random hexamer primers, incubated at 70°C for 10 min and cooled on ice. Alternatively, 10 μg total RNA were mixed with 500 ng (dT)17primer, incubated at 70°C for 10 min and cooled on ice. The labelling reaction was performed in a 12.5 μl reaction volume using 2.5 μl 5x RT buffer (Invitrogen, Karlsruhe, Germany), 1.25 μl 0.1 M DTT (Invitrogen), 1 μl dNTP-mix (dGAT, 5 mM each), 0.5 μl dCTP (3 mM), 0.5 μl RNasin (40 U/μl), 0.5 μl Cy-3 or Cy-5 labelled dCTP (1 mM, Amersham Biosciences, Freiburg, Germany) and 1 μl Superscript II reverse transcriptase (100 U/μl, Invitrogen). The mixture was incubated for 1 h at 42°C, and the reaction was stopped by addition of 1.25 μl 50 mM EDTA (pH 8). The input RNA was removed by hydrolysis with 5 μl 1 M NaOH at 65°C for 10 min, followed by neutralisation with 1 μl 5 M acetic acid. Cy-3 and Cy-5 labelled probes were combined, precipitated with isopropanol, washed with 70% ethanol, air dried and resuspended in 30 μl hybridisation buffer. The hybridisation buffer consisted of 350 μl 2x DIG-EASY (Roche Diagnostics), 210 μl MilliQ-H2O, 70 μl 50x Denhardt's solution [19] and 70 μl Cot1-DNA (10 ng/μl). The 30 μl probe was heat denatured at 65°C for 2 min and hybridised to cDNA glass microarrays in a hybridisation chamber (Corning Inc., Acton, USA) over night at 37°C. After hybridisation the microarrays were washed at 25°C in a water bath with 1x SSC containing 0.1% SDS once for 15 min, and twice for 10 min. This was followed by washing twice with 0.1x SSC containing 0.1% SDS for 10 min. Washing steps were completed in 70% and 95% ethanol before the microarrays were dried.

                  Image quantification

                  The hybridised arrays were scanned with the GenePix 4000 B microarray scanner (Axon Instruments Inc., Union City, CA, USA), and the scanned images were analysed using GenePix Pro 4.0 software (Axon Instruments). Spot intensities were obtained by subtracting the median brightness of a region around the spot from the median of the region within.

                  Data analysis

                  For each of the datasets α ... η (table2), the spot intensities were calibrated and transformed by the VSN method [20]. Differences between the resulting values are referred to as "generalised log ratios". All 14.694 spots from each slide were analysed, without filtering or thresholding procedures.

                  The reproducibility of different amplification reactions was measured using theFstatistic to compare the variance for each genewithinandbetweengroups of hybridisations using amplified samples. TheFvalue is defined as the ratio of thebetweenandwithingroups mean squaresmSS B andmSS W , whereas the mean square is defined as the sum of squares divided by the degrees of freedom:

                  http://static-content.springer.com/image/art%3A10.1186%2F1471-2164-5-29/MediaObjects/12864_2003_Article_131_Equa_HTML.gif

                  given http://static-content.springer.com/image/art%3A10.1186%2F1471-2164-5-29/MediaObjects/12864_2003_Article_131_IEq1_HTML.gif and http://static-content.springer.com/image/art%3A10.1186%2F1471-2164-5-29/MediaObjects/12864_2003_Article_131_IEq2_HTML.gif ,

                  wherey ikl denotes the measured value of thelth repetition of thekth amplification reaction for RNA transcripti.Kdenotes the number of all amplification reactions andL k the number of replicas in groupk.

                  Declarations

                  Acknowledgements

                  We thank Patricia McCabe and Ruprecht Kuner for critical reading of the manuscript and helpful comments. This work was supported by the German Federal Ministry for Science and Education (BMBF; FKZ: 1706300).

                  Authors’ Affiliations

                  (1)
                  Department of Molecular Genome Analysis, German Cancer Research Centre
                  (2)
                  Department of Gynecology and Obstetrics, Städtische Kliniken Bielefeld
                  (3)
                  Institute of Molecular Biology and Cell Culture Technology, Mannheim University of Applied Sciences

                  References

                  1. Beer DG, Kardia SL, Huang CC, Giordano TJ, Levin AM, Misek DE, Lin L, Chen G, Gharib TG, Thomas DG, Lizyness ML, Kuick R, Hayasaka S, Taylor JM, Iannettoni MD, Orringer MB, Hanash S:Gene-expression profiles predict survival of patients with lung adenocarcinoma. Nat Med2002,8:816–824.PubMed
                  2. Dyrskjot L, Thykjaer T, Kruhoffer M, Jensen JL, Marcussen N, Hamilton-Dutoit S, Wolf H, Orntoft TF:Identifying distinct classes of bladder carcinoma using microarrays. Nat Genet2003,33:90–96.View ArticlePubMed
                  3. Hedenfalk I, Ringner M, Ben-Dor A, Yakhini Z, Chen Y, Chebil G, Ach R, Loman N, Olsson H, Meltzer P, Borg A, Trent J:Molecular classification of familial non-BRCA1/BRCA2 breast cancer. Proc Natl Acad Sci U S A2003,100:2532–2537.View ArticlePubMed
                  4. Nutt CL, Mani DR, Betensky RA, Tamayo P, Cairncross JG, Ladd C, Pohl U, Hartmann C, McLaughlin ME, Batchelor TT, Black PM, Von Deimling A, Pomeroy SL, Golub TR, Louis DN:Gene Expression-based Classification of Malignant Gliomas Correlates Better with Survival than Histological Classification. Cancer Res2003,63:1602–1607.PubMed
                  5. van 't Veer LJ, Dai H, van de Vijver MJ, He YD, Hart AA, Mao M, Peterse HL, van der Kooy K, Marton MJ, Witteveen AT, Schreiber GJ, Kerkhoven RM, Roberts C, Linsley PS, Bernards R, Friend SH:Gene expression profiling predicts clinical outcome of breast cancer. Nature2002,415:530–536.View ArticlePubMed
                  6. Bernards R, Weinberg RA:A progression puzzle. Nature2002,418:823.View ArticlePubMed
                  7. van de Vijver MJ, He YD, van't Veer LJ, Dai H, Hart AA, Voskuil DW, Schreiber GJ, Peterse JL, Roberts C, Marton MJ, Parrish M, Atsma D, Witteveen A, Glas A, Delahaye L, van der Velde T, Bartelink H, Rodenhuis S, Rutgers ET, Friend SH, Bernards R:A gene-expression signature as a predictor of survival in breast cancer. N Engl J Med2002,347:1999–2009.View ArticlePubMed
                  8. 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 Biotechnol2002,20:940–943.View ArticlePubMed
                  9. Puskas LG, Zvara A, Hackler L., Jr., Van Hummelen P:RNA amplification results in reproducible microarray data with slight ratio bias. Biotechniques2002,32:1330–4, 1336, 1338, 1340.PubMed
                  10. Smith L, Underhill P, Pritchard C, Tymowska-Lalanne Z, Abdul-Hussein S, Hilton H, Winchester L, Williams D, Freeman T, Webb S, Greenfield A:Single primer amplification (SPA) of cDNA for microarray expression analysis. Nucleic Acids Res2003,31:e9.View ArticlePubMed
                  11. Van Gelder RN, von Zastrow ME, Yool A, Dement WC, Barchas JD, Eberwine JH:Amplified RNA synthesized from limited quantities of heterogeneous cDNA. Proc Natl Acad Sci U S A1990,87:1663–1667.View ArticlePubMed
                  12. Hu L, Wang J, Baggerly K, Wang H, Fuller GN, Hamilton SR, Coombes KR, Zhang W:Obtaining reliable information from minute amounts of RNA using cDNA microarrays. BMC Genomics2002,3:16.View ArticlePubMed
                  13. Pabon C, Modrusan Z, Ruvolo MV, Coleman IM, Daniel S, Yue H, Arnold L. J., Jr.:Optimized T7 amplification system for microarray analysis. Biotechniques2001,31:874–879.PubMed
                  14. Wang E, Miller LD, Ohnmacht GA, Liu ET, Marincola FM:High-fidelity mRNA amplification for gene profiling. Nat Biotechnol2000,18:457–459.View ArticlePubMed
                  15. Zhao H, Hastie T, Whitfield ML, Borresen-Dale AL, Jeffrey SS:Optimization and evaluation of T7 based RNA linear amplification protocols for cDNA microarray analysis. BMC Genomics2002,3:31.View ArticlePubMed
                  16. Ripamonti M, Canevari S, Menard S, Mezzanzanica D, Miotti S, Orlandi R, Rilke F, Tagliabue E, Colnaghi MI:Human carcinoma cell lines xenografted in athymic mice: biological and antigenic characteristics of an intraabdominal model. Cancer Immunol Immunother1987,24:13–18.View ArticlePubMed
                  17. Volz J, Koster S, Spacek Z, Paweletz N:The influence of pneumoperitoneum used in laparoscopic surgery on an intraabdominal tumor growth. Cancer1999,86:770–774.View ArticlePubMed
                  18. Volz J, Volz-Koster S, Kanis S, Klee D, Ahlert C, Melchert F:Modulation of tumor-induced lethality after pneumoperitoneum in a mouse model. Cancer2000,89:262–266.View ArticlePubMed
                  19. Sambrook J, Russel DW:Molecular Cloning: A Laboratory Manuel. 3rd Edition EditionCold Spring Harbor, NY, Cold Spring Harbor Laboratory Press 2001.
                  20. Huber W, Von Heydebreck A, Sultmann H, Poustka A, Vingron M:Variance stabilization applied to microarray data calibration and to the quantification of differential expression. Bioinformatics2002,18 Suppl 1:S96-S104.PubMed

                  Copyright

                  © Schneider et al 2004

                  This article is published under license to BioMed Central Ltd. This is an Open Access article: verbatim copying and redistribution of this article are permitted in all media for any purpose, provided this notice is preserved along with the article's original URL.

                  Advertisement