Skip to main content

Evaluating methods for ranking differentially expressed genes applied to microArray quality control data

Abstract

Background

Statistical methods for ranking differentially expressed genes (DEGs) from gene expression data should be evaluated with regard to high sensitivity, specificity, and reproducibility. In our previous studies, we evaluated eight gene ranking methods applied to only Affymetrix GeneChip data. A more general evaluation that also includes other microarray platforms, such as the Agilent or Illumina systems, is desirable for determining which methods are suitable for each platform and which method has better inter-platform reproducibility.

Results

We compared the eight gene ranking methods using the MicroArray Quality Control (MAQC) datasets produced by five manufacturers: Affymetrix, Applied Biosystems, Agilent, GE Healthcare, and Illumina. The area under the curve (AUC) was used as a measure for both sensitivity and specificity. Although the highest AUC values can vary with the definition of "true" DEGs, the best methods were, in most cases, either the weighted average difference (WAD), rank products (RP), or intensity-based moderated t statistic (ibmT). The percentages of overlapping genes (POGs) across different test sites were mainly evaluated as a measure for both intra- and inter-platform reproducibility. The POG values for WAD were the highest overall, irrespective of the choice of microarray platform. The high intra- and inter-platform reproducibility of WAD was also observed at a higher biological function level.

Conclusion

These results for the five microarray platforms were consistent with our previous ones based on 36 real experimental datasets measured using the Affymetrix platform. Thus, recommendations made using the MAQC benchmark data might be universally applicable.

Background

Identification of differentially expressed genes (DEGs) under different conditions is an important goal in microarray-based gene expression analysis. For this identification, new gene ranking methods have been developed and comparative studies have been performed [1–8]. Several evaluation metrics for comparing these methods have also been used [8–17]. A major metric is the area under the curve (AUC), which evaluates both sensitivity and specificity of the methods simultaneously [8, 14–17]. A high AUC value indicates that a set of true DEGs are highly concentrated at the top in ranked gene lists.

Currently, by virtue of the MicroArray Quality Control (MAQC) project [18, 19], researchers can comprehensively evaluate many competing methods with regard to reproducibility, another important metric [9–14]. The MAQC study provides a large number of benchmark datasets measured using different microarray platforms and at different test sites for a set of common samples (so-called "Samples A-D"; for details, see Materials and Methods). This enables us to evaluate inter-site (or intra-platform) and inter-platform reproducibility. The study concluded that a fold-change-based method, an average difference (AD) statistic for log-scale data, has the most reproducible results of ranked gene lists when inter-site reproducibility for the DEG lists was evaluated [18]. However, many other existing methods, such as rank products (RP) [4], were not compared in the original study. Moreover, new methods such as weighted average difference (WAD) [8] have been developed since the MAQC study. Evaluations for those methods are therefore important for creating up-to-date guidelines.

We recently reported that WAD outperformed AD, which was recommended by the MAQC study with regard to inter-site reproducibility. We also reported that the use of WAD or RP, in conjunction with suitable preprocessing algorithms dedicated to the Affymetrix (AFX) GeneChip data, can increase both sensitivity and specificity of the results [14]. However, as the datasets we used were derived from only one platform (i.e., the AFX data), the applicability of the conclusions to other platforms, such as Agilent One-Color (AG1) and Applied Biosystems (ABI) data, should be determined.

In this study, we analyzed the MAQC data [18, 19] measured using five platforms: ABI, AG1, GE Healthcare (GEH), and Illumina (ILM), as well as AFX. We report the suitable gene ranking methods having high sensitivity, specificity, and inter-site and inter-platform reproducibility for these data. We also report the capability of the recommended methods for generating consistent results at higher biological function levels as well as lower DEG list levels. Our recommendations for the five platforms are essentially the same as the previous ones, which were based on 36 real experimental datasets measured using only the AFX platform [14].

Methods

MAQC expression data

Expression data derived from five platforms, i.e., AFX, ABI, AG1, GEH, and ILM, were obtained from the MAQC website [20] and analyzed. The MAQC project produced four sample types: Sample A, a universal human reference RNA; Sample B, a human brain reference RNA; Sample C, consisting of 75 and 25 % Samples A and B respectively; and Sample D, consisting of 25 and 75% Samples A and B respectively [18]. Five replicate experiments for each of the four sample types at three (for ABI, AG1, GEH, and ILM) and six (for AFX) test sites were conducted. The downloaded data were pre-processed using the individual manufacturer protocols. More detail about pre-processing is given in Ref. [18]. The gene ranking methods were independently applied for "Sample A versus B" and "Sample C versus D" comparisons.

For the GEH data, signals from multiple probes with the same ID (e.g., 12 probes with the ID "GE1071362" were present) were averaged. This produced 53,517 unique probes from the original 54,359 probes. Then, signals under a minimum value for signals flagged as 'Good' quality were set to the minimum value (i.e., 0.169). This substitution of a floor value was made for 21.75% (698,424/(53,517 probes × 60 samples)) of the signals. For reproducing the research, the R-code for the preprocessing of GEH data is available in Additional file 1. After the preprocessing of the GEH data, log2 transformation was performed. The GEH data, i.e., 53,517 probes, were used for analyzing sensitivity and specificity. For the other platforms (ABI, AFX, AG1, and ILM), the expression data on a log2 scale were analyzed.

MAQC RT-PCR data

Quantitative real-time polymerase chain reaction (RT-PCR) data were obtained from the MAQC website [20] and used for defining "true" DEGs. The original MAQC study [19] produced TaqMan (TAQ) assay data on 997 genes, StaRT-PCR (GEX) assay data on 205 genes, and QuantiGene (QGN) assay data on 244 genes. Of these, we only used the TAQ assay data on 905 genes and the GEX assay data on 161 genes because those genes were common among the five different microarray platforms.

To determine the DEGs, we applied two metrics for the log2 transformed data: a t-test and an AD statistic (i.e., log-fold-change). The false discovery rate (FDR) for the p-values from the t-test was calculated using the Benjamini and Hochberg approach [21]. DEGs were selected based on FDR < 0.05 or an absolute AD value > 1. For comparing "Sample A versus B" using the TAQ data, the FDR threshold defined 817 genes as DEGs (TAQ_AB_FDR = 817, denoted as "Platform_Comparison_Metric ") and the AD threshold produced 553 DEGs (TAQ_AB_AD = 553). Similarly, 767, 351, 144, 106, 114, and 60 DEGs were determined by TAQ_CD_FDR, TAQ_CD_AD, GEX_AB_FDR, GEX_AB_AD, GEX_CD_FDR, and GEX_CD_AD, respectively. The R-codes for analyzing the TAQ and GEX data are available in Additional files 2 and 3, respectively.

Gene ranking

The genes were ranked according to the differential expression using R (ver. 2.11.0) [22] and Bioconductor [23]. Following on from our previous study [8, 14], eight methods for analyzing two-class data were compared: the WAD [8], AD [8], fold change (FC), RP [4], moderated t statistic (modT) [3], significance analysis of microarrays (samT) [1], shrinkage t statistic (shrT) [7], and intensity-based moderated t statistic (ibmT) [6]. A higher absolute value for all methods except for RP indicates a higher degree of differential expression. The genes were ranked in descending order of the absolute value.

RP is only a rank-based non-parametric method. Different from the other methods, RP independently handles up-regulated genes and down-regulated genes under one class and therefore produces two separate ranked gene lists. In other words, RP outputs two values per gene. To obtain one RP value per gene for comparison with the results of the other methods, we defined a lower value as a net value for the gene. A small net value for RP is therefore evidence of differential expression. This procedure is the same as that described in Ref. [14]. The genes were ranked in ascending order of the net value.

For evaluation of the methods based both on sensitivity and specificity, an AUC value was calculated for each ranked gene list, where true DEGs were mainly defined by two metrics (i.e., FDR and AD) as mentioned above. For evaluation based on inter-site and inter-platform reproducibility, the percentages of overlapping genes (POGs) and Spearman's rank correlation coefficients between different ranked gene lists were calculated.

Functional analysis

Analyzing gene set enrichment is a common approach in functional analysis where functionally or structurally related sets of genes regarding differential expression are assessed. Methods for this purpose include Gene Set Enrichment Analysis (GSEA) [24] and Parametric Analysis of Gene Set Enrichment (PAGE) [25]. These methods calculate the differences between two classes of samples (e.g., Sample A versus B) by using associative statistics, such as a signal-to-noise metric [24] or AD [25], and output one score per gene set. Different associative statistics with a particular method (e.g., GSEA with AD and GSEA with samT) could result in different top-rank lists of gene sets.

To compare the performances of the eight gene ranking methods as simply as possible, we calculate the average of the ranks of the genes belonging to a gene set as an "enrichment score". Consider a gene ranking method that produces m associative statistics S= {s 1 , s 2 , ..., s m } composed of m genes. Given R= {r 1 , r 2 , ..., r m } as the corresponding ranks of the values of the statistics and G= {g 1 , g2, ..., g n } as a gene set composed of n genes, the enrichment score for the gene set, E G , is simply calculated as. For example, consider three genes in a gene set ranked as 2nd, 8th, and 14th respectively. The enrichment score is simply calculated as E G = (2+8+14)/3 = 8. Accordingly, if a gene set has a small E G value, it is expected to be enriched (or de-regulated) in the two classes analyzed.

The gene sets were obtained from the Molecular Signatures Databases (MSigDB) ver. 3.0 [24]. We used 186 gene sets associated with the Kyoto Encyclopedia of Genes and Genomes (KEGG) Pathway. Since the analysis of gene set enrichment provides an enrichment score for each gene set, we obtained 186 enrichment scores for each gene ranking method. The pathways were analyzed for each of the eight gene ranking methods, each of the two comparisons ("Sample A versus B" and "C versus D"), each of the three test sites, and each of the five microarray platforms. Accordingly, the analysis produced a total of (8 × 2 × 3 × 5 =) 240 lists, each of which had 186 scores. Since the comparisons were performed for each platform, each list was denoted as "Method_Comparison_Site ".

Gene set information is available as gene symbols, while expression data are measured on the probe basis. Inevitably, many multiple probes are mapped to the same gene on existing microarray platforms. For example, three probes ('201127_s_at', '201128_s_at', and '210337_s_at') are mapped to the ACLY gene on the AFX platform, while only one probe is mapped to that gene on the other platforms ('102218' for ABI, '17711' for AG1, '199102' for GEH, and 'GI_38569420-S' for ILM). In this case, expression signals for the three probes on the AFX platform should be collapsed to avoid inflating the enrichment scores and to facilitate biological interpretation of the analysis results [24]. We averaged the signals from the probes mapped to the same gene.

The mapping information for the five platforms was obtained from the Gene Expression Omnibus [26] as 'GPL570' for AFX, 'GPL2986' for ABI, 'GPL1708' for AG1, 'GPL2895' for GEH, and 'GPL2507' for ILM. For the GEH data, the original data composed of 54,359 probes were collapsed. After this, expression data composed of 20,723 unique genes for AFX, 16,275 for ABI, 19,053 for AG1, 19,352 for GEH, and 16,992 for ILM on a log2 scale were analyzed. The POGs and Spearman's rank correlation coefficients for the 186 enrichment scores between test sites (or different platforms) were calculated for estimating inter-site (or inter-platform) reproducibility of the individual gene ranking methods at the function level.

Results

Investigation of suitable gene ranking methods having high sensitivity, specificity, and (inter-site and inter-platform) reproducibility is the main purpose of this work. Similar to our previous study on the AFX platform [14], eight gene ranking methods (WAD, AD, FC, RP, modT, samT, shrT, and ibmT) were evaluated. Compared to the abundant real experimental data measured using the AFX platform in our previous study, the data for other platforms are considerably limited. We therefore analyzed the MAQC data measured using five platforms (AFX, ABI, AG1, GEH, and ILM). The results for the MAQC AFX data were compared with our previous results, where 36 real experimental datasets measured using the AFX platform were analyzed [14]. We will follow the terminology used in our previous study [14] for consistency.

Sensitivity and specificity

We and other researchers believe the human samples compared in the MAQC study, cell lines (sample A) versus brain tissues (sample B), are clearly distinct, i.e., such comparisons are unrealistic in practical biological research [27–29]. This may weaken the confidence of any results obtained using the MAQC benchmark data. To mitigate this concern, we first analyzed the MAQC data measured using the AFX platform and compared the current results to our previous ones (where the best method was either WAD or RP) [14]. The evaluation frameworks were set to be as similar as possible to those in the previous study [14]. Two assay technologies for RT-PCR data (TAQ and GEX) [19] were used to determine true DEGs, and an AUC value was used to compare sensitivity and specificity simultaneously [30–32]. Currently, there is no convincing rationale for defining DEGs from the RT-PCR data. We here define DEGs by two metrics: a FDR threshold (FDR < 0.05) and an absolute AD threshold (|AD| > 1; AD corresponds to "log-fold-change").

The average AUC values for the AFX data across six test sites are shown in Table 1. WAD and RP, recommended by Ref. [14], show the best performance in six of the eight cases. Note that the ibmT method also shows high AUC values as a whole; it is always the second best when WAD is the best. These observations are similar to those in our previous study [14], suggesting that findings derived from the MAQC data measured using the other four platforms (ABI, AG1, GEH, and ILM) could be applicable for practical microarray research.

Table 1 Average AUC values for AFX data

The average AUC values for ABI, AG1, GEH, and ILM data across three test sites are shown in Table 2. Similar to the results for the AFX data, the best was either WAD, RP, or ibmT in all but one (i.e., TAQ_AB_AD ) case. The ibmT method may be suitable for ABI data. RP was overall good when DEGs were determined by the AD metric for TAQ RT-PCR data. WAD was overall good when DEGs were determined by the FDR metric from a t-test, regardless of the choice of RT-PCR technology (i.e., TAQ or GEX).

Table 2 Average AUC values for ABI, AG1, GEH, and ILM data

Note that these results were based on only two sets of DEGs arbitrarily defined by commonly used procedures (i.e., FDR < 0.05 and |AD| > 1). As the AUC values vary with the thresholds for individual metrics, the effects of changing the parameters should also be analyzed. Table 3 shows the gene ranking methods having the highest AUC values when the top X genes are defined as DEGs (X = 100, 200, 300, ..., and 800 for TAQ data and 50, 100, and 150 for GEX data). We confirmed the superiority of the WAD, RP, and ibmT methods against various thresholds.

Table 3 List of best methods with different threshold settings for defining DEGs

For reproducing the research, the R-codes for obtaining the AUC values shown in Tables 1 and 2 are given in Additional files 4, 5, 6, 7 and 8. The raw AUC values for all methods shown in Table 3 and the R-code for obtaining a fraction of the values are given in Additional files 9 and 10, respectively.

Inter-site reproducibility at lower DEG list level

Next, we examined the inter-site reproducibility of ranked gene lists among three test sites for each platform. Our previous study using AFX data reports that WAD is the most reproducible method when POG is evaluated (Figure one in Ref. [14]). The POG values for the 100 top-ranked genes among three test sites for the five platforms AFX, ABI, AG1, GEH, and ILM are shown in Table 4. The complete POG values are given in Additional files 11 and 12.

Table 4 Inter-site reproducibility at lower DEG list level

The POG values for WAD were clearly higher than those for the other seven methods, especially in the "Sample C versus D" comparison. The results of that comparison should take precedence over those of "Sample A versus B". This is because sample C is a 3:1 and sample D is a 1:3 mixture of samples A and B, respectively, and real biological samples consist of various kinds of tissues and cells. These results suggest that the use of WAD can be promising for users concerned about inter-site reproducibility at the DEG list level.

Note that the POG values shown in Table 4 only describe inter-site (or intra-platform) reproducibility among three test sites for each gene ranking method. Accordingly, there may be concern about the potential difference of the ranked gene lists between WAD and the other methods. Figure 1 shows the results of hierarchical clustering applied to a total of 48 ranked gene lists (8 methods × 2 comparisons × 3 test sites = 48 lists; each list is therefore denoted as "Method_Comparison_Site "). Each list consists of (a) 54,676 gene ranks for AFX data, (b) 32,879 for ABI, (c) 41,059 for AG1, (d) 53,517 for GEH, and (e) 47,294 for ILM. When we focus on the dendrograms for the ILM data (Figure 1e), six distinct clusters can be seen, each consisting of eight lists with the same sample comparison and the same test site but different gene ranking methods. We also observe this trend for the ABI data (Figure 1b). Different from those for the ILM and ABI data, the dendrogram for the AG1 data shows a distinct WAD cluster (Figure 1c). This cluster is clearly located in a larger cluster mainly consisting of gene lists from the "Sample A versus B" comparison. These results indicate that the ranked gene lists produced from WAD are undoubtedly not outliers.

Figure 1
figure 1

Dendrograms of average-linkage hierarchical clustering. Dendrograms for (a) AFX, (b) ABI, (c) AG1, (d) GEH, and (e) ILM data. Spearman's rank correlation coefficient is used as a similarity metric; left-hand scale represents (1- correlation coefficient). Each gene list is denoted as "Method_Comparison_Site ". For example, a ranked list calculated from samT when comparing C versus D for site 2 data is denoted as "samT_CD_2".

Similar to the POG values in Table 4 the inter-site reproducibility can also be evaluated by calculating the Spearman's rank correlation coefficients of two ranked gene lists. The high inter-site reproducibility of WAD was also confirmed using this metric. The correlation coefficients between any two test sites (i.e., "Site 1 versus 2", "1 versus 3", and "2 versus 3") are given in Additional file 13.

Inter-site reproducibility at higher biological function level

In general, the higher the reproducibility of a gene ranking method at a lower DEG list level, the higher the reproducibility the method should have at a higher biological function level. To clarify this idea, we performed a functional analysis by estimating gene set enrichment. There are many sophisticated methods for analyzing gene set enrichment, such as GSEA [24] and PAGE [25]. We calculated the average rank for genes belonging to each gene set as the "enrichment score." This was to manipulate the eight gene ranking methods equally and to perform the comparison as simply as possible (for details, see Materials and Methods).

We performed the functional analysis for 186 gene sets associated with KEGG Pathway. The POG values for different settings of X top-ranked gene sets among three test sites for the five platforms (AFX, ABI, AG1, GEH, and ILM) are shown in Figure 2. While the best methods depend on the number of top-ranked gene sets, WAD seems to have relatively good performance compared to the other methods. We also observed this trend when Spearman's rank correlation coefficients were evaluated as another metric (see Additional file 14). The raw data for Figure 2 are available in Additional file 15. The R-code for making the left side of Figure 2(a) is also available in Additional file 16.

Figure 2
figure 2

Inter-site reproducibility at higher biological function level. POG values for different settings of X top-ranked gene sets among three test sites. For AFX data, POG values were calculated based on first three test sites.

Inter-platform reproducibility at higher biological function level

Note that the numbers of genes belonging to each gene set differ across platforms. Only 5 of the 186 KEGG Pathway gene sets (STEROID_BIOSYNTHESIS, VALINE_LEUCINE_AND_ISOLEUCINE_BIOSYNTHESIS, FOLATE_BIOSYNTHESIS, TERPENOID_BACKBONE_BIOSYNTHESIS, and RENIN_ANGIOTENSIN_SYSTEM ) were identical across the 5 platforms. To estimate the overall similarity among platforms, we calculated an intersection-union ratio for each gene set. For example, there are 61, 61, 60, 54, and 61 genes in the respective platforms AFX, ABI, AG1, GEH, and ILM in the GLYCOLYSIS_GLUCONEOGENESIS pathway. As the intersection and union across the 5 platforms are 53 and 61, respectively, the ratio is 53/61 = 0.87. Detailed information on the 186 gene sets is available in Additional file 17. The average of the 186 intersection-union ratios was 0.83. Arguably, we assert this high value enables a comparison between platforms.

The POG values for different settings of X top-ranked gene sets among the five platforms are shown in Figure 3. Similar to the results for inter-site reproducibility shown in Figure 2, the results for inter-platform reproducibility do not indicate the best gene ranking method because one method does not always outperform the others, especially when comparing ~25% of the top-ranked gene sets. This result is somewhat disappointing because only ~20 top-ranked gene sets have been evaluated at most [13, 33]. Nevertheless, as a whole, WAD outperforms the others. We also observed this trend for the 825 Gene Ontology gene sets (file name "c5.bp.v3.0.symbols.gmt") and for all 3,272 curated gene sets ("c2.all.v3.0.symbols.gmt"), obtained from MSigDB ver. 3.0 [24] (data not shown). The raw data for Figure 3 are available in Additional file 18.

Figure 3
figure 3

Inter-platform reproducibility at higher biological function level. POG values for different settings of X top-ranked gene sets among five platforms. Values correspond to intersections among five platforms shown in Figure 2 when comparing (a) A versus B and (b) C versus D.

Discussion

We investigated suitable gene ranking methods across various microarray platforms using a large set of MAQC microarray data. In accordance with previous evaluation frameworks [8, 14], sensitivity and specificity of the methods were evaluated using an AUC metric where genes were determined to be DEGs or non-DEGs based on RT-PCR data. Reproducibility was evaluated using a POG metric and Spearman's correlation coefficient. The current conclusions are essentially the same as those of our previous studies [8, 14]: the best of the eight methods investigated was (i) either WAD, RP, or ibmT, when both sensitivity and specificity were evaluated and (ii) WAD when reproducibility was evaluated. As we proposed one of the methods (WAD) that is compared here, this study could be regarded as the third in a series of studies on WAD. In the first study [8], the high sensitivity and specificity of the above-mentioned methods (especially WAD) were demonstrated using two spike-in datasets [34, 35] and 36 real experimental datasets. In the second study [14], the sensitivity and specificity of the eight gene ranking methods coupled with nine pre-processing algorithms dedicated to the Affymetrix (AFX) GeneChip data were further evaluated using the 36 real experimental datasets, and high reproducibility of WAD was demonstrated using the MAQC data. However, our previous studies were based only on datasets measured using the AFX platform. We have demonstrated in this study that those recommendations are also applicable for the other four microarray platforms.

Sensitivity and specificity

The current results indicate that care should be taken when choosing both a gene ranking method for analyzing microarray data and a metric for analyzing RT-PCR data. In addition, there are more effective combinations of the two choices: the WAD method in conjunction with the FDR metric (WAD/FDR), RP/AD, and ibmT/FDR. These combinations are logical because RP/AD are both fold-change-based statistics and ibmT/FDR are both t-statistic-based ones. Although WAD is a fold-change-based method, the overall performance of WAD/FDR is superior to that of ibmT/FDR (Table 3). To our knowledge, no study has reported suitable combinations of the two choices to facilitate agreement between the results. However, note that excessive value may have been placed on the RT-PCR data when determining these combinations. Although the RT-PCR technology is generally considered the "gold-standard" assay for measuring gene expression by biologists [36], it has not been proven as the best. In this regard, our proposed combinations may be useful for discussion on the agreement or disagreement between the results obtained by the two technologies.

Reproducibility

As mentioned, our previous study reports that WAD outperforms the other methods when inter-site reproducibility is evaluated for AFX microarray data [14]. In the current study, we demonstrate that WAD has clearly higher inter-site reproducibility when DEG lists produced from different test sites are compared, irrespective of the choice of microarray platform (Table 4 and Additional files 11, 12).

Li et al. recently reported high inter-platform reproducibility at the biological function level obtained by using a commercially available tool [13]. Only the samT method [1] was used in their functional analysis, without any justification, despite citation of the original MAQC paper where AD was recommended [18]. In addition, there should be "levels" of reproducibility; the use of better gene ranking methods should increase the level of inter-platform reproducibility at the biological function level. However, comparison of many gene ranking methods at the biological function level using such specialized tools is difficult practically for the following reasons: (1) gene ranking methods of interest may not be used, (2) only a set of significant gene sets can be manipulated, and so on. In particular, related to reason (1), to our knowledge there is no method that can also analyze gene set enrichment in conjunction with non-parametric gene ranking methods such as RP [4]. In this sense, our simple non-parametric approach to estimating gene set enrichment can be deemed a new method. Our functional analysis using this simple approach demonstrates that WAD has both high intra- and inter-platform reproducibility at the biological function level (Figures 2 and 3). We observed similarity between the top-rank lists produced by this simple procedure and those by two widely used methods (i.e., GSEA [24] and PAGE [25]) when some real experimental datasets used in Refs. [24, 25, 37] were applied (data not shown), but such a comparison is out of the scope of this study. Improving the approach of using WAD for functional analysis is future work.

In this study, we could perform a more detailed analysis with regard to reproducibility by using MAQC datasets. However, the original MAQC study [18] has often been criticized for its evaluation framework, where the correlation coefficient and the POG were used as measures for reproducible ranked gene lists [27–29]. Chen et al. [29] stated that the POG metric does not reflect the accuracy (i.e., sensitivity and specificity) of a ranked gene list. This is true in a sense because, for example, POG values higher than those of WAD can easily be obtained if we rank genes according to relative average signal intensity (i.e., a weight term (w ) in a WAD statistic; see Figure one in Ref. [14]). This w statistic is clearly inferior to that obtained from the other specialized gene ranking methods such as WAD when the AUC value (both sensitivity and specificity) is evaluated [8, 14]. Evaluation based only on the reproducibility of ranked gene lists is therefore insufficient. In the current study, however, we compare the eight gene ranking methods using the AUC value (both sensitivity and specificity) as well as the POG value (reproducibility). In addition, the most reproducible method (i.e., WAD) is one of the three methods that found to have high sensitivity and specificity.

Other pre-processing methods

In this study, we analyzed publicly available datasets that were already pre-processed according to the MAQC philosophy [18]. Our additional processing for the datasets was basically only the log2 transformation (see Materials and Methods). As demonstrated by our previous study [14] and another study of the MAQC project [38], there are many other choices for pre-processing data (for example, see Refs. [39–42]) and different choices may result in different conclusions. Thus, the other pre-processing choices should be evaluated next.

Conclusion

In conclusion, the superiority of WAD was confirmed with regard to sensitivity, specificity, and reproducibility by using the MAQC datasets. We thus recommend the use of WAD for more effective transcriptome analysis.

Abbreviations

ABI:

Applied Biosystems

AD:

average difference

AFX:

Affymetrix

AG1:

Agilent one-color microarray

AUC:

area under the curve

FDR:

false discovery rate

GEH:

GE Healthcare

GEX:

StaRT-PCR

ibmT:

intensity-based moderated t-statistic

ILM:

Illumina

MAQC:

MicroArray Quality Control (project)

POG:

percentage of overlapping genes

RP:

rank products

samT:

significance analysis of microarrays

shrT:

shrinkage t-statistic

TAQ:

TaqMan

WAD:

weighted average difference (method)

References

  1. Tusher VG, Tibshirani R, Chu G: Significance analysis of microarrays applied to the ionizing radiation response. Proc Natl Acad Sci USA 2001, 98(9):5116–5121. 10.1073/pnas.091062498

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  2. Pepe MS, Longton G, Anderson GL, Schummer M: Selecting differentially expressed genes from microarray experiments. Biometrics 2003, 59(1):133–142. 10.1111/1541-0420.00016

    Article  PubMed  Google Scholar 

  3. Smyth GK: Linear models and empirical bayes methods for assessing differential expression in microarray experiments. Stat Appl Genet Mol Biol 2004., 3(1): Article 3 Article 3

  4. Breitling R, Armengaud P, Amtmann A, Herzyk P: Rank products: a simple, yet powerful, new method to detect differentially regulated genes in replicated microarray experiments. FEBS Lett 2004, 573(1–3):83–92. 10.1016/j.febslet.2004.07.055

    Article  CAS  PubMed  Google Scholar 

  5. Jeffery IB, Higgins DG, Culhane AC: Comparison and evaluation of methods for generating differentially expressed gene lists from microarray data. BMC Bioinformatics 2006, 7: 359. 10.1186/1471-2105-7-359

    Article  PubMed Central  PubMed  Google Scholar 

  6. Sartor MA, Tomlinson CR, Wesselkamper SC, Sivaganesan S, Leikauf GD, Medvedovic M: Intensity-based hierarchical Bayes method improves testing for differentially expressed genes in microarray experiments. BMC Bioinformatics 2006, 7: 538. 10.1186/1471-2105-7-538

    Article  PubMed Central  PubMed  Google Scholar 

  7. Opgen-Rhein R, Strimmer K: Accurate ranking of differentially expressed genes by a distribution-free shrinkage approach. Stat Appl Genet Mol Biol 2007, 6: Article 9.

    Google Scholar 

  8. Kadota K, Nakai Y, Shimizu K: A weighted average difference method for detecting differentially expressed genes from microarray data. Algorithm Mol Biol 2008, 3: 8. 10.1186/1748-7188-3-8

    Article  Google Scholar 

  9. Shi L, Tong W, Fang H, Scherf U, Han J, Puri RK, Frueh RK, Goodsaid FM, Guo L, Su Z, Han T, Fuscoe JC, Xu ZA, Patterson TA, Hong H, Xie Q, Perkins RG, Chen JJ, Casciano DA: Cross-platform comparability of microarray technology: intra-platform consistency and appropriate data analysis procedures are essential. BMC Bioinformatics 2005, 6(Suppl 2):S12. 10.1186/1471-2105-6-S2-S12

    Article  PubMed Central  PubMed  Google Scholar 

  10. Kadota K, Ye J, Nakai Y, Terada T, Shimizu K: ROKU: a novel method for identification of tissue-specific genes. BMC Bioinformatics 2006, 7: 294. 10.1186/1471-2105-7-294

    Article  PubMed Central  PubMed  Google Scholar 

  11. Kadota K, Konishi T, Shimizu K: Evaluation of two outlier-detection-based methods for detecting tissue-selective genes from microarray data. Gene Regul Syst Biol 2007, 1: 9–15.

    Google Scholar 

  12. Stafford P, Brun M: Three methods for optimization of cross-laboratory and cross-platform microarray expression data. Nucleic Acids Res 2007, 35: e72. 10.1093/nar/gkl1133

    Article  PubMed Central  PubMed  Google Scholar 

  13. Li Z, Su Z, Wen Z, Shi L, Chen T: Microarray platform consistency is revealed by biologically functional analysis of gene expression profiles. BMC Bioinformatics 2009, 10(Suppl 11):S12. 10.1186/1471-2105-10-S11-S12

    Article  PubMed Central  PubMed  Google Scholar 

  14. Kadota K, Nakai Y, Shimizu K: Ranking differentially expressed genes from Affymetrix gene expression data: methods with reproducibility, sensitivity, and specificity. Algorithm Mol Biol 2009, 4: 7. 10.1186/1748-7188-4-7

    Article  Google Scholar 

  15. Pearson RD: A comprehensive re-analysis of the Golden Spike data: towards a benchmark for differential expression methods. BMC Bioinformatics 2008, 9: 164. 10.1186/1471-2105-9-164

    Article  PubMed Central  PubMed  Google Scholar 

  16. Murie C, Woody O, Lee AY, Nadon R: Comparison of small n statistical tests of differential expression applied to microarrays. BMC Bioinformatics 2009, 10: 45. 10.1186/1471-2105-10-45

    Article  PubMed Central  PubMed  Google Scholar 

  17. Bullard JH, Purdom E, Hansen KD, Dudoit S: Evaluation of statistical methods for normalization and differential expression in mRNA-Seq experiments. BMC Bioinformatics 2010, 11: 94. 10.1186/1471-2105-11-94

    Article  PubMed Central  PubMed  Google Scholar 

  18. MAQC Consortium: The MicroArray Quality Control (MAQC) project shows inter- and intraplatform reproducibility of gene expression measurements. Nat Biotechnol 2006, 24(9):1151–1161. 10.1038/nbt1239

    Article  PubMed Central  Google Scholar 

  19. Canales RD, Luo Y, Willey JC, Austermiller B, Barbacioru CC, Boysen C, Hunkapiller K, Jensen RV, Knight CR, Lee KY, Ma Y, Maqsodi B, Papallo A, Peters EH, Poulter K, Ruppel PL, Samaha RR, Shi L, Yang W, Zhang L, Goodsaid FM: Evaluation of DNA microarray results with quantitative gene expression platforms. Nat Biotechnol 2006, 24(9):1115–1122. 10.1038/nbt1236

    Article  CAS  PubMed  Google Scholar 

  20. The MAQC website[http://www.fda.gov/ScienceResearch/BioinformaticsTools/MicroarrayQualityControlProject/]

  21. Benjamini Y, Hochberg Y: Controlling the False Discovery Rate: A practical and powerful approach to multiple testing. J Roy Statist Soc Ser B 1995, 57(1):289–300.

    Google Scholar 

  22. R Development Core Team: R: A language and environment for statistical computing. In R Foundation for Statistical Computing. Vienna, Australia; 2010.

    Google Scholar 

  23. 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 JY, Zhang J: Bioconductor: open software development for computational biology and bioinformatics. Genome Biol 2004, 5(10):R80. 10.1186/gb-2004-5-10-r80

    Article  PubMed Central  PubMed  Google Scholar 

  24. Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, Paulovich A, Pomeroy SL, Golub TR, Lander ES, Mesirov JP: Gene set enrichment analysis: A knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci USA 2005, 102: 15545–15550. 10.1073/pnas.0506580102

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  25. Kim SY, Volsky DJ: PAGE: Parametric Analysis of Gene Set Enrichment. BMC Bioinformatics 2005, 6: 144. 10.1186/1471-2105-6-144

    Article  PubMed Central  PubMed  Google Scholar 

  26. Barrett T, Troup DB, Wilhite SE, Ledoux P, Rudnev D, Evangelista C, Kim IF, Soboleva A, Tomashevsky M, Marshall KA, Phillippy KH, Sherman PM, Muertter RN, Edgar R: NCBI GEO: archive for high-throughput functional genomic data. Nucleic Acids Res 2009, 37: D885-D890. 10.1093/nar/gkn764

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  27. Strauss E: Arrays of Hope. Cell 2006, 127: 657–659. 10.1016/j.cell.2006.11.005

    Article  CAS  PubMed  Google Scholar 

  28. Liang P: MAQC papers over the cracks. Nature Biotechnol 2007, 25: 27–29. 10.1038/nbt0107-27

    Article  CAS  Google Scholar 

  29. Chen JJ, Hsueh HM, Delongchamp RR, Lin CJ, Tsai CA: Reproducibility of microarray data: a further analysis of microarray quality control (MAQC) data. BMC Bioinformatics 2007, 8: 412. 10.1186/1471-2105-8-412

    Article  PubMed Central  PubMed  Google Scholar 

  30. Chen Z, McGee M, Liu Q, Scheuermann RH: A distribution free summarization method for Affymetrix GeneChip arrays. Bioinformatics 2007, 23(3):321–327. 10.1093/bioinformatics/btl609

    Article  CAS  PubMed  Google Scholar 

  31. Hu J, Xu J: Density based pruning for identification of differentially expressed genes from microarray data. BMC Genomics 2010, 11: S3.

    Article  PubMed Central  PubMed  Google Scholar 

  32. Hashimoto-Torii K, Kawasawa YI, Kuhn A, Rakic P: Combined transcriptome analysis of fetal human and mouse cerebral cortex exposed to alcohol. Proc Natl Acad Sci USA 2011, 108(10):4212–4217. 10.1073/pnas.1100903108

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  33. Abatangelo L, Maglietta R, Distaso A, D'Addabbo A, Creanza TM, Mukherjee S, Ancona N: Comparative study of gene set enrichment methods. BMC Bioinformatics 2009, 10: 275. 10.1186/1471-2105-10-275

    Article  PubMed Central  PubMed  Google Scholar 

  34. Irizarry RA, Bolstad BM, Collin F, Cope LM, Hobbs B, Speed TP: Summaries of Affymetrix GeneChip probe level data. Nucleic Acids Res 2003, 31(4):e15. 10.1093/nar/gng015

    Article  PubMed Central  PubMed  Google Scholar 

  35. Cope LM, Irizarry RA, Jaffee HA, Wu Z, Speed TP: A benchmark for Affymetrix GeneChip expression measures. Bioinformatics 2004, 20(3):323–331. 10.1093/bioinformatics/btg410

    Article  CAS  PubMed  Google Scholar 

  36. Qin LX, Beyer RP, Hudson FN, Linford NJ, Morris DE, Kerr KF: Evaluation of methods for oligonucleotide array data via quantitative real-time PCR. BMC Bioinformatics 2006, 7: 23. 10.1186/1471-2105-7-23

    Article  PubMed Central  PubMed  Google Scholar 

  37. Nakai Y, Hashida H, Kadota K, Minami M, Shimizu K, Matsumoto I, Kato H, Abe K: Up-regulation of genes related to the ubiquitin-proteasome system in the brown adipose tissue of 24-h-fasted rats. Biosci Biotechnol Biochem 2008, 72(1):139–148. 10.1271/bbb.70508

    Article  CAS  PubMed  Google Scholar 

  38. Shippy R, Fulmer-Smentek S, Jensen RV, Jones WD, Wolber PK, Johnson CD, Pine PS, Boysen C, Guo X, Chudin E, Sun YA, Willey JC, Thierry-Mieg J, Thierry-Mieg D, Setterquist RA, Wilson M, Lucas AB, Novoradovskaya N, Papallo A, Turpaz Y, Baker SC, Warrington JA, Shi L, Herman D: Using RNA sample titrations to assess microarray platform performance and normalization techniques. Nat Biotechnol 2006, 24(9):1123–1131. 10.1038/nbt1241

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  39. Barbacioru CC, Wang Y, Canales RD, Sun YA, Keys DN, Chan F, Poulter KA, Samaha RR: Effect of various normalization methods on Applied Biosystems expression array system data. BMC Bioinformatics 2006, 7: 533. 10.1186/1471-2105-7-533

    Article  PubMed Central  PubMed  Google Scholar 

  40. Zahurak M, Parmigiani G, Yu W, Scharpf RB, Berman D, Schaeffer E, Shabbeer S, Cope L: Pre-processing Agilent microarray data. BMC Bioinformatics 2007, 8: 142. 10.1186/1471-2105-8-142

    Article  PubMed Central  PubMed  Google Scholar 

  41. Diez D, Alvarez R, Dopazo A: Codelink: an R package for analysis of GE healthcare gene expression bioarrays. Bioinformatics 2007, 23(9):1168–1169. 10.1093/bioinformatics/btm072

    Article  CAS  PubMed  Google Scholar 

  42. Shi W, Banerjee A, Ritchie ME, Gerondakis S, Smyth GK: Illumina WG-6 BeadChip strips should be normalized separately. BMC Bioinformatics 2009, 10: 372. 10.1186/1471-2105-10-372

    Article  PubMed Central  PubMed  Google Scholar 

Download references

Acknowledgements and Funding

This study was supported by Grants-in-Aid for Scientific Research (21710208 and 22128008) to KK from the Japanese Ministry of Education, Culture, Sports, Science and Technology (MEXT).

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Koji Kadota.

Additional information

Authors' contributions

KK performed the analyses and wrote the paper. KS supervised the critical discussion by KK. All authors read and approved the final manuscript.

Electronic supplementary material

12859_2010_4605_MOESM1_ESM.R

Additional file 1:R-code for pre-processing GEH data. Input file must be obtained from http://edkb.fda.gov/MAQC/MainStudy/upload/norm_MAQC_GEH_123_Median1.zip. After execution of this R-code, expression data consisting of 53,517 unique probes can be obtained as an output file arbitrarily named "norm_MAQC_GEH_123_Median1_unique.txt". Note that the website with the normalized MAQC data was unavailable under review. The raw data may be availabl from. http://www.fda.gov/ScienceResearch/BioinformaticsTools/Arraytrack/ or http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE5350. (R 1 KB)

12859_2010_4605_MOESM2_ESM.R

Additional file 2:R-code for analyzing TAQ data. Two input files must first be obtained from the following websites: "http://www.nature.com/nbt/journal/v24/n9/extref/nbt1239-S5.txt" and "http://edkb.fda.gov/MAQC/MainStudy/upload/norm_TAQ_1_POLR2A.zip". After execution of this R-code, FDR and AD statistics for each comparison (i.e., "Sample A versus B" and "C versus D") can be obtained in an output file arbitrarily named "data_TAQ.txt". Note that the website with the normalized MAQC data was unavailable under review. The raw data may be available from http://www.fda.gov/ScienceResearch/BioinformaticsTools/Arraytrack/ or http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE5350. (R 2 KB)

12859_2010_4605_MOESM3_ESM.R

Additional file 3:R-code for analyzing GEX data. Two input files must first be obtained from the following websites: "http://www.nature.com/nbt/journal/v24/n9/extref/nbt1239-S5.txt" and "http://edkb.fda.gov/MAQC/MainStudy/upload/norm_MAQC_GEX_1_ACTB.zip". After execution of this R-code, FDR and AD statistics for each comparison (i.e., "Sample A versus B" and "C versus D") can be obtained in an output file arbitrarily named "data_GEX.txt". Note that the website with the normalized MAQC data was unavailable under review. The raw data may be available from http://www.fda.gov/ScienceResearch/BioinformaticsTools/Arraytrack/ or http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE5350. (R 2 KB)

12859_2010_4605_MOESM4_ESM.R

Additional file 4:R-code for calculating average AUC values for AFX data across six test sites shown in Table 1. Two input files are required. One is the AFX gene expression data file obtained from http://edkb.fda.gov/MAQC/MainStudy/upload/norm_MAQC_AFX_123456_qPLIER16.zip. The other includes the DEG information derived from individual RT-PCR data. Specifically, the filename should be "data_TAQ.txt" (default) or "data_GEX.txt", generated by executing Additional file 2 or 3, respectively. After execution of this R-code, the average AUC values shown in Table 1 can be obtained in an output file arbitrarily named "result_AFX.txt". Note that the website with the normalized MAQC data was unavailable under review. The raw data may be available from http://www.fda.gov/ScienceResearch/BioinformaticsTools/Arraytrack/ or http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE5350. (R 57 KB)

12859_2010_4605_MOESM5_ESM.R

Additional file 5:R-code for calculating average AUC values for ABI data across three test sites shown in Table 2. Two input files are required. One is the ABI gene expression data file obtained from http://edkb.fda.gov/MAQC/MainStudy/upload/norm_MAQC_ABI_123_qNorm.zip. The other includes the DEG information derived from individual RT-PCR data. Specifically, the filename should be "data_TAQ.txt" (default) or "data_GEX.txt", generated by executing Additional file 2 or 3, respectively. After execution of this R-code, the average AUC values for ABI data shown in Table 2 can be obtained in an output file arbitrarily named "result_ABI.txt". Note that the website with the normalized MAQC data was unavailable under review. The raw data may be available from http://www.fda.gov/ScienceResearch/BioinformaticsTools/Arraytrack/ or http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE5350. (R 30 KB)

12859_2010_4605_MOESM6_ESM.R

Additional file 6:R-code for calculating average AUC values for AG1 data across three test sites shown in Table 2. Two input files are required. One is the AG1 gene expression data file obtained from http://edkb.fda.gov/MAQC/MainStudy/upload/norm_MAQC_AG1_123_Median1GeneSpring.zip. The other includes the DEG information derived from individual RT-PCR data. Specifically, the filename should be "data_TAQ.txt" (default) or "data_GEX.txt", generated by executing Additional file 2 or 3, respectively. After execution of this R-code, the average AUC values for AG1 data shown in Table 2 can be obtained in an output file arbitrarily named "result_AG1.txt". Note that the website with the normalized MAQC data was unavailable under review. The raw data may be available from http://www.fda.gov/ScienceResearch/BioinformaticsTools/Arraytrack/ or http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE5350. (R 30 KB)

12859_2010_4605_MOESM7_ESM.R

Additional file 7:R-code for calculating average AUC values for GEH data across three test sites shown in Table 2. Two input files are required. One is the GEH gene expression data file obtained by executing Additional file 1, and the output file should be named "norm_MAQC_GEH_123_Median1_unique.txt". The other includes the DEGs information derived from individual RT-PCR data. Specifically, the filename should be "data_TAQ.txt" (default) or "data_GEX.txt", generated by executing Additional file 2 or 3, respectively. After execution of this R-code, the average AUC values for GEH data shown in Table 2 can be obtained in an output file arbitrarily named "result_GEH.txt". (R 30 KB)

12859_2010_4605_MOESM8_ESM.R

Additional file 8:R-code for calculating average AUC values for ILM data across three test sites shown in Table 2. Two input files are required. One is the ILM gene expression data file obtained from http://edkb.fda.gov/MAQC/MainStudy/upload/norm_MAQC_ILM_123_qNorm16.zip. The other includes the information about DEGs derived from individual RT-PCR data. Specifically, the filename should be "data_TAQ.txt" (default) or "data_GEX.txt", generated by executing Additional file 2 or 3, respectively. After execution of this R-code, the average AUC values for ILM data shown in Table 2 can be obtained in an output file arbitrarily named "result_ILM.txt". Note that the website with the normalized MAQC data was unavailable under review. The raw data may be available from http://www.fda.gov/ScienceResearch/BioinformaticsTools/Arraytrack/ or http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE5350. (R 30 KB)

12859_2010_4605_MOESM9_ESM.XLS

Additional file 9:Average AUC values for all methods regarding Table 3. Table 3 was obtained based on these values. (XLS 66 KB)

12859_2010_4605_MOESM10_ESM.R

Additional file 10:R-code for obtaining average AUC values for AFX data across six test sites when top X genes are defined as DEGs. Two input files are required. One is the AFX gene expression data file obtained from http://edkb.fda.gov/MAQC/MainStudy/upload/norm_MAQC_AFX_123456_qPLIER16.zip. The other includes the DEG information derived from individual RT-PCR data. Specifically, the filename should be "data_TAQ.txt" (default) or "data_GEX.txt", generated by executing Additional file 2 or 3, respectively. After execution of this R-code, a fraction of the average AUC values for the AFX data when top X genes are defined as DEGs (X= 100 by default), shown in Additional file 9, can be obtained in an output file arbitrarily named "result_AFX_100.txt". Note that the website with the normalized MAQC data was unavailable under review. The raw data may be available from http://www.fda.gov/ScienceResearch/BioinformaticsTools/Arraytrack/ or http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE5350. (R 57 KB)

12859_2010_4605_MOESM11_ESM.PPT

Additional file 11:POG values for given numbers of top-ranked genes among three test sites. (a) AFX data, (b) ABI data, (c) AG1 data, (d) GEH data, and (e) ILM data. Number of DEGs (X) is shown on x-axis (log-scale). Percentage of genes (POG) for X top-ranked genes among three test sites is shown on y-axis. (PPT 144 KB)

12859_2010_4605_MOESM12_ESM.XLS

Additional file 12:Raw data for Additional file11. There are a total of ten sheets (five platforms × two comparisons). (XLS 4 MB)

Additional file 13:Correlation coefficients of ranked gene lists between any two test sites. (XLS 32 KB)

Additional file 14:Correlation coefficients of 186 ranked lists of gene sets between any two test sites. (XLS 33 KB)

12859_2010_4605_MOESM15_ESM.XLS

Additional file 15:Raw data for Figure 2. There are a total of ten sheets (five platforms × two comparisons). (XLS 336 KB)

12859_2010_4605_MOESM16_ESM.R

Additional file 16:R-code for making left side of Figure 2(a). The R-code assumes the input filename is "POG_AFX_KEGG_AB.txt," which contains data for the first sheet (i.e., the sheet named "AFX_AvsB") in Additional file 15. (R 2 KB)

12859_2010_4605_MOESM17_ESM.XLS

Additional file 17:Numbers of genes belonging to each of 186 KEGG Pathway gene sets in individual platforms. (XLS 54 KB)

Additional file 18:Raw data for Figure 3. There are two sheets. (XLS 86 KB)

Authors’ original submitted files for images

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Kadota, K., Shimizu, K. Evaluating methods for ranking differentially expressed genes applied to microArray quality control data. BMC Bioinformatics 12, 227 (2011). https://doi.org/10.1186/1471-2105-12-227

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1471-2105-12-227

Keywords