Skip to main content
  • Methodology article
  • Open access
  • Published:

MASQOT: a method for cDNA microarray spot quality control

Abstract

Background

cDNA microarray technology has emerged as a major player in the parallel detection of biomolecules, but still suffers from fundamental technical problems. Identifying and removing unreliable data is crucial to prevent the risk of receiving illusive analysis results. Visual assessment of spot quality is still a common procedure, despite the time-consuming work of manually inspecting spots in the range of hundreds of thousands or more.

Results

A novel methodology for cDNA microarray spot quality control is outlined. Multivariate discriminant analysis was used to assess spot quality based on existing and novel descriptors. The presented methodology displays high reproducibility and was found superior in identifying unreliable data compared to other evaluated methodologies.

Conclusion

The proposed methodology for cDNA microarray spot quality control generates non-discrete values of spot quality which can be utilized as weights in subsequent analysis procedures as well as to discard spots of undesired quality using the suggested threshold values. The MASQOT approach provides a consistent assessment of spot quality and can be considered an alternative to the labor-intensive manual quality assessment process.

Background

At present, the DNA microarray technology allows simultaneous monitoring of the expression levels of thousands of genes. The technique produces large and complex datasets that are relatively easy to generate but non-trivial to analyze and extract information from. So far, much of the data mining efforts have been focused on the statistical analysis (see, for instance [14]) and less on acquiring high quality data from the image analysis. Image analysis is the process of extracting information from the scanned microarray images, which is an important step due to the sequential nature of the microarray analysis [5]. Consequently, problems in the initial steps have large impact on the interpretation of the final results from the experiment.

A number of technical issues during microarray preparation potentially affect the spot quality.

  • Low signal intensity is perhaps the most generally acknowledged property that affects spot quality due to the subsequent problems in distinguishing signal from noise for spots with weak signals. Weak signal intensities should result from physiologically low expression levels but might also be related to surface properties of the slide, signal bleaching, scanner problems or incomplete or irregular hybridization.

  • Intensity distribution issues appear as regions of pixels containing signals that clearly deviate from the average signal, typically as distinct sub-areas within the foreground region. As the signal of any given spot on a microarray slide is expected to be uniform over the entire spot foreground area, intensity distribution issues are usually a consequence of non-specific binding or irregular distribution of the printed DNA on the slide.

  • Morphological issues refer to unexpected shape-related variations of the spot foreground region. This includes very small or very large spot sizes, low spot circularities or spot mixing. Size aberrations might be a consequence of precipitates or impurities in the printing solution or needle clogging during printing. Spots are expected to be roughly circular in shape, but manufacturing issues might result in deviation from the circularity norm. Furthermore, imperfections on the slide or washing problems might cause the dye from several spots to mix, referred to as bleeding, making the separation of these signals difficult or even impossible.

  • Background issues appear as intensity fluctuations in the local background region immediately surrounding the foreground region. An increase in local background intensity or variance compared to the global slide background typically result from dye contaminants due to non-specific binding or incomplete washing.

Microarray spot quality control is essentially the identification and removal of spots with properties that cause the subsequent interpretation of the signal from these spots to be unreliable or misleading. Analogously, it is the recognition of characteristics that enable dependable interpretations and conclusions from the generated data. The field of microarray spot quality control has been largely neglected in the past but has recently become an area of interest. Existing documented semi-automatic methodologies include Bayesian networks [6] as well as linear combinations of quality parameters allegedly related to the quality of the spot [710]. Interestingly, manual evaluation is still a common quality assessment procedure, despite the labor-intensive nature of visually inspecting spots in the order of hundreds of thousands. Increasing availability and usage of microarrays in the field of transcriptomics has caused experiments to involve an escalating number of slides, which in turn has highlighted the primary bottleneck of manual quality assessment.

Partial least squares (PLS) is a generalized regression method which aims to maximize the covariance between the X (descriptor) and Y (response) matrices. PLS can handle large data sets of multi-collinear and noisy data with moderate amounts of missing data in both X and Y. PLS-DA can be seen as a special case of PLS where the response matrix Y is categorical (numerically represented as 0 or 1) and determines class belonging of observations. PLS-DA has been widely applied in microarray analysis (see, for instance [11, 12]) as well as other areas of life science (see, for instance [13, 14]). For a more detailed description of the properties of PLS, please consult [1517] and references therein.

Here, we propose the m icroa rray s pot q uality co nt rol (MASQOT) methodology for assessment of cDNA microarray spot quality, outlined in figure 1. A set of existing and novel spot descriptors were identified that aimed to characterize spot quality in terms of physical attributes of the spot. Prior to the extraction of descriptors, manual assessment of the spot quality was performed independently by three experienced microarray users on roughly eighty thousand spots in order to provide a sufficiently large data set of known quality.

Figure 1
figure 1

Flowchart of the classification procedure. The classification process involves an 8-bit image, optimized for segmentation, as well as a 32-bit image, used for information extraction. During the training phase, visual classification results are required while this is not necessary for external data.

Spot descriptors were subsequently subjected to multivariate discriminant analysis by means of PLS-DA with the aim to categorize spots of low quality and spots of high quality by treating these spots as separate classes. The utilized descriptors aimed to describe foreground and background irregularity measures, spot morphology and foreground density attributes that were potentially useful for discriminating between the reliable (not bad) and unreliable (bad) spots. For instance, the circularity measure is a descriptor ranging from 0% to 100%. If the circularity descriptor approaches its minimum value, the predicted class belonging should typically be higher (closer to 1) for the bad class compared to the not bad class. The coefficient of variation for the foreground and background regions is an example of an employed descriptor that illustrates reverse characteristics; higher values should typically provide greater class conformity with the bad class compared to the not bad class. However, all employed descriptors together contribute to the regression model and consequently to the final class determination at varying degrees depending on the properties of the spot.

The MASQOT approach aims to provide a consistent assessment of spot quality, applicable to various types of microarray data, thus avoiding the labor-intensive manual quality assessment process. The methodology generates continuous values of spot quality which can be utilized to discard spots of undesired quality or used as weights in subsequent analysis procedures.

Results

Five cDNA microarray slides using the Populus second generation microarray slide layout (POP2) where the samples originate from a previous investigation of leaf development (Sjödin et al, in preparation) were used for classification training. Five additional POP2 slides, not included in the training set, were employed for external validation. Segmentation of raw images was performed using an implementation of the Seeded Region Growing (SRG) algorithm [18]. The properties of each spot were subsequently characterized using a large set of descriptors allegedly linked to spot quality. These properties include foreground and background variability measures, spot morphology and foreground intensity distribution measures. Please consult Additional file 1 for a complete list of all the utilized descriptors.

Following segmentation, spots were inspected by three experienced microarray users and independently assigned to the two quality categories {bad, not bad}. Spots in the bad category consisted of all the spots that were classified as bad by at least one of the experienced users while the remaining spots were categorized as not bad. For classification and evaluation purposes, the spots in the bad category were subsequently partitioned into different sub-classes based on visual properties as described in table 1. This can be seen as characterizing each spot as exhibiting

Table 1 The different sub-classes of bad spots.
  1. 1.

    No issues (not bad); or

  2. 2.

    Foreground issues (FI); or

  3. 3.

    Background issues (BI); or

  4. 4.

    Any combination of 2 and 3

To avoid confounding of properties, only spots displaying pure issues (entries 1–3) were used in the classification training, although all spots were used in the model evaluation. The three-class problem was subjected to multivariate analysis by means of PLS [16] regression coupled with discriminant analysis (PLS-DA) with the aim to discriminate not bad spots from FI spots and BI spots. The result from the PLS-DA regression model is a predicted class conformity (CC) value for each of the classes: not bad (CCnb), foreground issues (CCFI) and background issues (CCBI) with the added restriction that CCnb = 1 - (CCFI + CCBI). Due to this restriction, only the conformity value of the not bad spots (CCnb) will be interpreted in the upcoming sections. A CCnb value approaching 1 denotes high compliance with the not bad class, which can be interpreted as a quality measure of the spot. Spots visually categorized as bad should thus exhibit a value of CCnb close to 0 whereas spots categorized as not bad should exhibit a value of CCnb close to 1.

Receiver Operating Characteristics (ROC) plot of the classifications of the POP2 training set and POP2 test set, respectively, is available in figure 2. A density plot of the CCnb value for the bad and not bad spots in the POP2 training set is shown in figure 3, showing a partial overlap between the discrimination of the two classes. Due to this overlap, discrimination accuracy was dependent on a threshold value t denoting the separation point between the two classes.

Figure 2
figure 2

Receiver Operating Characteristics (ROC) plot. The relation between true positives (bad spots classified as bad) and false positives (not bad spots classified as bad) for the training and test data. The solid line denotes training data whereas the dashed line denotes test data.

Figure 3
figure 3

Density plot of the predicted class conformity of the not bad class. A class conformity value of 1 signifies perfect class conformity while a value of 0 signifies no class conformity. The dashed line illustrates the density for the prediction of the bad spots in the POP2 training set whereas the solid line illustrates the density of the prediction of the not bad spots in the POP2 training set.

Spots with a predicted CCnb value below t were classified as bad whereas the remaining spots were classified as not bad. The threshold value can be set more or less stringently depending on the quality filtering requirements. This is illustrated in figures 4a–b, depicting different views of interpreting the classification accuracy of the POP2 training set and the POP2 test set. The threshold values were set either to maximize the overall classification accuracy or to maximize the class-wise classification accuracy. The overall classification accuracies for the POP2 training set (38 627 spots) and POP2 test set (39 421 spots) were calculated using equation 1 for all t in the interval (0,1) using CCnb values for the {bad, not bad} spots. The classification accuracy peak at a level of approximately 98% where t = 0.4 (see figure 4a). Predicted class-wise accuracies were calculated using equation 1 employing CCnb values for the bad spots and CCnb values for the not bad spots separately. The predicted class-wise accuracies intersect at a level of 95% (see figure 4b) for t = 0.5. Exact classification accuracies per sub-class, based on the intersection threshold value t = 0.5 as illustrated in figure 4b, are available in table 2.

Figure 4
figure 4

Relationship between classification accuracy and threshold value for the POP2 data. The threshold value t defines the boundary between bad and not bad spots for the POP2 training set (38 627 spots) and the POP2 test set (39 421 spots). Spots with a predicted class conformity value for the not bad class (CCnb) below the threshold value t are classified as bad while the remaining spots are classified as not bad. a) Overall classification accuracy vs. threshold value calculated as the fraction of correctly classified spots in the data set for a given threshold value. The solid line represents the POP2 training set whereas the dashed line represents the POP2 test set. The dotted vertical line at threshold value t = 0.4 illustrates an approximate maximum. b) Classification accuracy of the bad and not bad spots vs. threshold value. For the POP2 training set, the solid line represents the classification accuracy of the not bad spots and the dashed line represents the classification accuracy of the bad spots. For the POP2 test set, the dot-dashed line represents the classification accuracy of the not bad spots and the long-dashed line represents the classification accuracy of the bad spots. The dotted vertical line at threshold value t = 0.5 denotes the intersection point.

Table 2 Classification accuracy of the POP2 training data. The classification accuracy for each sub-class as calculated using threshold value t = 0.5.

The presented MASQOT approach was compared to three existing quality control methods: the composite quality score qcom proposed by Wang et al [8], the mean-median correlation factor mmcorr evaluated by Tran et al [9] and the coefficient of variation (CV) parameter CVspot evaluated by Sauer et al [10]. Threshold values for all quality control parameters were set to achieve maximum overall classification accuracy. The result, shown in table 3, demonstrates that the MASQOT approach provides a greater level of class discrimination for the POP2 test set compared to the remaining evaluated quality control methods.

Table 3 Comparison to other quality control methods. The presented quality control parameter CCnb was compared to the composite quality score qcom, the mean-median correlation factor mmcorr and the CVspot value. Threshold values for all quality control parameters were set to maximize overall classification accuracy. The classification accuracy was determined from classification of the POP2 test set.

Discussion

cDNA microarray spot quality control is, in many aspects, a complex problem. Naturally, the automatic assessment of quality of each spot is highly reliant on the characterization of the spot. Incorrect approximations of the spatial location will affect the properties of the segmented foreground region, which in turn will influence the values of the quality control descriptors. In such a sequential process, where each step is dependent on the preceding steps, errors will propagate down-stream at a high rate. However, the most striking intricacy is perhaps the visual assessment, which is the foundation of this computer-based classification, where even experienced microarray users tend to disagree. As shown from the results presented here, it is the spots with unanimous visual quality assessment that are the most complicated to reproduce accurately. This disagreement stems from the more fundamental issue of defining 'quality' and in understanding the basal aspects that affect this quality.

The approach described here aims to assess the technical precision of each spot, which is believed to be linked to the biological accuracy (see [19] for a discussion regarding precision and accuracy in the microarray field). It should thus be stated, in this context, that lack of precision in a microarray spot measurement does not necessarily infer lack of accuracy. However, it is arguably reasonable to handle spots of questionable precision with specific care during the analysis procedure to aid the concluding biological interpretations.

The spot quality control assessment is commonly treated as a discrete problem (essentially, separating 'bad' spots from 'good' spots) but the spot quality varies on a continuous scale, ranging from very bad to very good. Instead of discarding spots, one might weight the spots according to the quality assessment. The concept of relative spot weights has previously been acknowledged in microarray normalization techniques (see, for instance [2, 20]) but might also prove to be valuable in additional analysis steps. However, evaluation of the usage and validity of spot weights based on the quality assessments provided here remains the scope of a future paper.

The rate of accuracy in prediction of the true positives (bad spots predicted as bad) and the true negatives (not bad spots predicted as not bad) are illustrated separately since these accuracies are not consistently of equal importance. For instance, depending on the user and the question at hand, it might be more important to avoid the risk of removing spots of decent quality than to eliminate all of the bad spots from the data set. Simply using the overall classification accuracy could be rather illusive, merely since the number of bad spots is much lower than the number of not bad spots in a typical data set.

The methodology presented here is based on the scaled sum of the intensities from both channels but can, with minor adjustments, also be based on single channel intensities. By using the scaled channel intensity levels, one avoids the risk of drowning information, in particular when there is a great difference in intensity level between the channels. In addition, the presented approach greatly resembles the visual illustrations of the spots which, by design, will provide an advantage in finding correlation between the visual quality assessment and the spot descriptors. Furthermore, it is more feasible for the average user to achieve a per-spot quality measure than a per-channel quality measure since this avoids raising questions with regard to what to do when only one channel is of low or moderate quality.

The recent advances in spot quality control have clearly shown that a good explanation of training data is possible using several different methodologies of varying complexity. However, very few efforts have been made to evaluate further aspects of the quality (for instance, more refined descriptors) and, most importantly, the reproducibility of the classification on external data. External reproducibility has been the major aim here, partly overshadowing the aim of internal reproducibility on the training data, which is shown by the clear agreement in accuracy between the independent POP2 training set and the POP2 test set.

Conclusion

The presented MASQOT technique provides a robust methodology for semi-automated cDNA microarray spot quality control with high accuracy of training data as well as external data compared to other evaluated methods. The MASQOT methodology generates non-discrete values of spot quality which can be utilized as weights in subsequent analysis procedures as well as to discard spots of undesired quality using the proposed threshold values.

Methods

Microarray preparation

Samples for the microarray slides used in this paper originate from an experiment of Populus tremula leaves, investigating regulation of leaf development (Sjödin et al, in preparation). The utilized microarray layout, referred to as POP2, consist of 25 278 single spotted cDNA clones from a recent assembly of more than 100 000 expressed sequence tags (ESTs) from the Populus genus [21]. All sequence information is available in the online sequence resource PopulusDB [22] and a full array layout is available for download from the online microarray resource UPSC-BASE [23].

Ten out of a total of 28 POP2 slides were randomly chosen for classification and were subsequently grouped into two equally large sets of five slides each; the POP2 training set and POP2 external test set. See additional file 2 for a complete list of the POP2 microarray slides used here. All POP2 slides were printed using a QArray arrayer (Genetix, Hampshire, U.K.). The preparation, labeling and hybridization of cDNA clones and mRNA samples were carried out according to the protocol described by Smith et al [24]. The arrays were scanned on a ScanArray 4000 (Perkin-Elmer Wellesley, MA) at 5 μm resolution to obtain raw image files for the red-fluorescent dye Cy5 and the green-fluorescent dye Cy3. All POP2 raw image files are available online for download at the UPSC-BASE microarray database [23] from experiment number 0013.

Image analysis

The workflow from scanned cDNA images to computer-based classification was separated into seven sub-procedures, outlined below and illustrated in figure 1.

1. Image merging generated a combined image from the intensity measurements of both the red-fluorescent dye Cy5 and the green-fluorescent dye Cy3, which was used in the subsequent gridding and segmentation steps.

2. Gridding attempted to identify the precise spatial center of the spots on the scanned microarray images.

3. Segmentation classified the pixels as either representing the cDNA expression level (foreground pixels) or an estimation of the local noise level (background pixels). In addition, a thin strip of pixels in the boundary region between the two segments (border pixels) was identified.

4. Information extraction refers to the characterization of the foreground and background regions from the segmentation process. In general terms, information extraction should provide a description of each region that is relevant in some sense (for instance, the spatial location of the foreground region or the foreground intensity level.) The focus here was on features that captured the overall quality of the spot.

5. Manual classification provided a measure of the spot quality by means of visual inspection carried out by three experienced microarray users.

6. Computer-based classification of spot quality (the training phase) generated a model for the differences between the spot quality classes using discriminant analysis based on the PLS regression method (PLS-DA).

7. Verification of the computer-based classification (the test phase) validated the predictive ability of the model using processed data not included in the training phase.

Image merging

Both segmentation and gridding were based on a combined eight-bit image constructed from the intensity measurements of the red-fluorescent dye Cy5 and the green-fluorescent dye Cy3. The merged eight-bit image lacks some details of the original images but is computationally efficient, in particular concerning memory requirements. Details of the utilized damping and scaling procedures are described by Yang et al [5] and briefly outlined below.

  • The intensity levels in both images were square-root transformed. The square-root transform utilizes damping, which ensures that the relative impact of high-intensity pixels is decreased during gridding and segmentation.

  • Median intensity values were computed from the transformed images.

  • A joint intensity value was calculated using the sum of the square-root transformed intensities from both channels scaled according to the median values, respectively.

  • Intensity values greater than 255 were truncated.

Gridding

Approximate spatial centers of each spot, referred to as the grid points, were manually located using an in-house developed Java application. This procedure is the only step in the classification process that requires user intervention. A more precise midpoint of the foreground region was found using a square pixel mask with the expected spot diameter (100 μm for the POP2 data) surrounding the initial grid point. The pixel mask was spatially reallocated in all directions, deviating at most 30 μm from the initial grid point, and the center position of the square pixel mask containing the highest total sum of intensities was selected as seed point.

Segmentation

The employed segmentation method was an implementation of the seeded region growing (SRG) algorithm, initially proposed by Adams and Bischof [18]. The SRG method has earlier been utilized in microarray spot segmentation by Yang et al [5]. Implementation details are available in additional file 3. The result from the segmentation process was a pixel mask categorizing each pixel into one of the four groups {foreground, border, background, un-assigned}. Each spot thus consisted of a distinct foreground region with the following characteristics:

  • All pixels within the foreground region were spatially connected.

  • No pixels overlapped with the foreground region of another spot.

  • Minor fluctuations in intensity level within the region were accepted.

  • The maximum Euclidean distance between any two pixels in the foreground region was restricted.

  • Spot circularity was not assumed.

Information extraction

The data utilized in the information extraction originates from the sum of the raw intensities of both the Cy5 and the Cy3 channels scaled according to the respective median intensity value. The scaling was applied to decrease the impact of the channel demonstrating the highest median intensity value. A large set of different features were extracted which were believed to be linked to spot quality and these were subsequently used in the upcoming computer-based classification. For purposes of repeatability and applicability to various types of microarray slides, all descriptors were corrected by an approximation of the slide background mean based on the mean intensity level of the local background regions from all spots on the slide.

Furthermore, spots where the saturation contents in at least one of the channels exceeded 10% of the total number of pixels, as suggested by Wang et al [8], were not included in the classification.

A complete table of all extracted descriptors including a description or definition is available in additional file 1. The descriptors aimed to capture foreground and background variability properties, spot morphology and foreground intensity and density properties.

Manual classification

The spots from ten POP2 slides were independently inspected by three experienced microarray users and assigned to the two quality categories {bad, not bad}. The spots in the bad category consisted of all the spots that were classified as bad by at least one of the experienced users while the remaining spots were categorized as not bad.

During the visual classification, the experienced users worked according to four basic rules of thumb related to the technical precision of each spot.

  • The signal within the foreground region should have low variability.

  • The foreground region should be circular.

  • The foreground region should be spatially located at the expected position.

  • The background region should have low variability and low intensity level compared to the global slide background.

The relation between these, that is, how much each factor was allowed to deviate, alone and in combination with other factors, was the task for the multivariate classification model. The utilized data sets, subsequent to segmentation, are available at additional file 4 (training set) and additional file 5 (test set). A summary of the manual classifications as performed by the experienced users is available in additional file 6.

The POP2 slides were randomly partitioned into two equally large sets of five slides each; the POP2 training set and the POP2 test set. For classification and evaluation purposes, the bad spots of the training set were subsequently divided into different sub-classes based on visual properties as described in table 1. The HIFI and LIFI sub-classes were used as representatives of the pure foreground issues (FI) during classification training. Analogously, the HIBI and LIBI sub-classes were used as representatives of the pure background issues (BI) during classification training. Typical examples of the described sub-classes can be found in additional file 7.

Computer-based classification

The computer-based classification was performed using PLS-DA as implemented in SIMCA-P+ 10.0 (Umetrics AB, Umeå, Sweden). Cross-validation [25] with seven groups was used to determine the number of latent variables. Prior to analysis, all descriptors were column-wise mean-centered and scaled to unit variance (UV) by dividing each descriptor with the standard deviation of the descriptor. The UV scaling procedure in combination with mean-centering translates the distribution of each descriptor to unit variance. Results and model statistics from the PLS-DA training phase are described in additional file 8.

Classification training was based on discriminant analysis of subsets of the foreground issues (FI) class, the background issues (BI) class and the not bad class. See table 1 for a more detailed description of the available sub-classes. Prior to discriminant analysis, a representative subset of each class consisting of 355 spots each was selected using D-optimal design [26] in order to eliminate the large differences in data set size between the three classes. See additional file 9 for the designed data set and additional file 10 for details regarding the D-optimal design.

Classification accuracies were calculated using equation 1, where n is the number of observations; t is a threshold value and x the predicted class conformity values for a given set of spots. Equation 1 utilizes the corr pred (i, y, t) function that returns 1 if i {bad} and y < t or if i {not bad} and y ≥ t; or 0 otherwise.

p r e d a c c ( x , t ) = 100 n i = 1 n c o r r p r e d ( i , x i , t ) ( 1 ) MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciGacaGaaeqabaqabeGadaaakeaacqWGWbaCcqWGYbGCcqWGLbqzcqWGKbazdaWgaaWcbaGaemyyaeMaem4yamMaem4yamgabeaakiabcIcaOiabdIha4jabcYcaSiabdsha0jabcMcaPiabg2da9maalaaabaGaeGymaeJaeGimaaJaeGimaadabaGaemOBa4gaamaaqahabaGaem4yamMaem4Ba8MaemOCaiNaemOCai3aaSbaaSqaaiabdchaWjabdkhaYjabdwgaLjabdsgaKbqabaGccqGGOaakcqWGPbqAcqGGSaalcqWG4baEdaWgaaWcbaGaemyAaKgabeaakiabcYcaSiabdsha0jabcMcaPaWcbaGaemyAaKMaeyypa0JaeGymaedabaGaemOBa4ganiabggHiLdGccaWLjaGaaCzcamaabmaabaGaeGymaedacaGLOaGaayzkaaaaaa@6051@

Abbreviations

PLS :

Partial Least Squares

PLS-DA :

Partial Least Squares Discriminant Analysis

SRG :

Seeded Region Growing

ROC :

Receiver Operating Characteristics

EST :

Expressed Sequence Tag

References

  1. Kerr MK, Martin M, Churchill GA: Analysis of variance for gene expression microarray data. J Comput Biol 2000, 7(6):819–837. 10.1089/10665270050514954

    Article  CAS  PubMed  Google Scholar 

  2. Smyth GK, Speed T: Normalization of cDNA microarray data. Methods 2003, 31(4):265–273. 10.1016/S1046-2023(03)00155-5

    Article  CAS  PubMed  Google Scholar 

  3. Wolfinger RD, Gibson G, Wolfinger ED, Bennett L, Hamadeh H, Bushel P, Afshari C, Paules RS: Assessing gene significance from cDNA microarray expression data via mixed models. J Comput Biol 2001, 8(6):625–637. 10.1089/106652701753307520

    Article  CAS  PubMed  Google Scholar 

  4. Efron B, Tibshirani R, Storey JD, Tusher V: Empirical Bayes analysis of a microarray experiment. J Am Stat Assoc 2001, 96(456):1151–1160. 10.1198/016214501753382129

    Article  Google Scholar 

  5. Yang YH, Buckley MJ, Dudoit S, Speed TP: Comparison of methods for image analysis on cDNA microarray data. J Comput Graph Stat 2002, 11(1):108–136. 10.1198/106186002317375640

    Article  Google Scholar 

  6. Hautaniemi S, Edgren H, Vesanen P, Wolf M, Jarvinen AK, Yli-Harja O, Astola J, Kallioniemi O, Monni O: A novel strategy for microarray quality control using Bayesian networks. Bioinformatics 2003, 19(16):2031–2038. 10.1093/bioinformatics/btg275

    Article  CAS  PubMed  Google Scholar 

  7. Chen YD, Kamat V, Dougherty ER, Bittner ML, Meltzer PS, Trent JM: Ratio statistics of gene expression levels and applications to microarray data analysis. Bioinformatics 2002, 18(9):1207–1215. 10.1093/bioinformatics/18.9.1207

    Article  CAS  PubMed  Google Scholar 

  8. Wang X, Ghosh S, Guo SW: Quantitative quality control in microarray image processing and data acquisition. Nucleic Acids Res 2001, 29(15):e75. 10.1093/nar/29.15.e75

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  9. Tran PH, Peiffer DA, Shin Y, Meek LM, Brody JP, Cho KW: Microarray optimizations: increasing spot accuracy and automated identification of true microarray signals. Nucleic Acids Res 2002, 30(12):e54. 10.1093/nar/gnf053

    Article  PubMed Central  PubMed  Google Scholar 

  10. Sauer U, Preininger C, Hany-Schmatzberger R: Quick and simple: quality control of microarray data. Bioinformatics 2005, 21(8):1572–1578. 10.1093/bioinformatics/bti238

    Article  CAS  PubMed  Google Scholar 

  11. Man MZ, Dyson G, Johnson K, Liao B: Evaluating methods for classifying expression data. J Biopharm Stat 2004, 14(4):1065–1084. 10.1081/BIP-200035491

    Article  PubMed  Google Scholar 

  12. Perez-Enciso M, Tenenhaus M: Prediction of clinical outcome with microarray data: a partial least squares discriminant analysis (PLS-DA) approach. Hum Genet 2003, 112(5–6):581–592.

    PubMed  Google Scholar 

  13. Eriksson L, Antti H, Gottfries J, Holmes E, Johansson E, Lindgren F, Long I, Lundstedt T, Trygg J, Wold S: Using chemometrics for navigating in the large data sets of genomics, proteomics, and metabonomics (gpm). Anal Bioanal Chem 2004, 380(3):419–429. 10.1007/s00216-004-2783-y

    Article  CAS  PubMed  Google Scholar 

  14. Jonsson P, Bruce SJ, Moritz T, Trygg J, Sjöström M, Plumb R, Granger J, Maibaum E, Nicholson JK, Holmes E, Antti H: Extraction, interpretation and validation of information for comparing samples in metabolic LC/MS data sets. Analyst 2005, 130(5):701–707. 10.1039/b501890k

    Article  CAS  PubMed  Google Scholar 

  15. Martens H, Naes T: Multivariate Calibration. Chichester , John Wiley & Sons; 1992.

    Google Scholar 

  16. Burnham AJ, MacGregor JF, Viveros R: Latent variable multivariate regression modeling. Chemom Intell Lab Syst 1999, 48(2):167–180. 10.1016/S0169-7439(99)00018-0

    Article  CAS  Google Scholar 

  17. Eriksson L, Johansson E, Kettaneh-Wold N, Wold S: Multi- and Megavariate Data Analysis: Principles and Applications. Umeå , Umetrics Academy; 2001.

    Google Scholar 

  18. Adams R, Bischof L: Seeded Region Growing. IEEE T Pattern Anal 1994, 16(6):641–647. 10.1109/34.295913

    Article  Google Scholar 

  19. van Bakel H, Holstege FC: In control: systematic assessment of microarray performance. EMBO Rep 2004, 5(10):964–969. 10.1038/sj.embor.7400253

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  20. 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.

    Google Scholar 

  21. Sterky F, Bhalerao RR, Unneberg P, Segerman B, Nilsson P, Brunner AM, Charbonnel-Campaa L, Lindvall JJ, Tandre K, Strauss SH, Sundberg B, Gustafsson P, Uhlen M, Bhalerao RP, Nilsson O, Sandberg G, Karlsson J, Lundeberg J, Jansson S: A Populus EST resource for plant functional genomics. Proc Natl Acad Sci USA 2004, 101(38):13951–13956. 10.1073/pnas.0401641101

    Article  PubMed Central  PubMed  Google Scholar 

  22. Sterky F, Bhalerao RR, Unneberg P, Segerman B, Nilsson P, Brunner AM, Charbonnel-Campaa L, Lindvall JJ, Tandre K, Strauss SH, Sundberg B, Gustafsson P, Uhlen M, Bhalerao RP, Nilsson O, Sandberg G, Karlsson J, Lundeberg J, Jansson S: PopulusDB: A Populus EST resource for plant functional genomics.[http://www.populus.db.umu.se/]

  23. Sjödin A, Bylesjö M: UPSC-BASE: Populus transcriptomics online.[http://www.upscbase.db.umu.se/]

  24. Smith CM, Rodriguez-Buey M, Karlsson J, Campbell MM: The response of the poplar transcriptome to wounding and subsequent infection by a viral pathogen. New Phytol 2004, 164(1):123–136. 10.1111/j.1469-8137.2004.01151.x

    Article  CAS  Google Scholar 

  25. Wold S: Cross Validatory Estimation of the Number of Components in Factor and Principal Components Models. Technometrics 1978, 20: 397–406.

    Article  Google Scholar 

  26. deAguiar PF, Bourguignon B, Khots MS, Massart DL, PhanThanLuu R: D-optimal designs. Chemom Intell Lab Syst 1995, 30(2):199–210. 10.1016/0169-7439(94)00076-X

    Article  Google Scholar 

Download references

Acknowledgements

This work was supported by grants from

• The Swedish Foundation for Strategic Research (MB, HA)

• The Knut and Alice Wallenberg Foundation (JT)

• The European Commission through the Directorate General Research within the Fifth Framework for Research – Quality of Life and Management of the Living Resources Programme, contract No.QLK5-CT-2002-00953 coordinated by the University of Southampton (AS, SJ).

• The Swedish Research Council (AS, SJ, DE, MS)

• EU-strategic funding (DE)

• The Functional Genomics Initiative at Swedish University of Agricultural Sciences (DE)

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Max Bylesjö.

Additional information

Authors' contributions

MB implemented the segmentation and gridding tools, performed visual classification, implemented the D-optimal design, conceived and generated the classification model and drafted the manuscript. DE generated the microarray slides, performed visual classification and helped to draft the manuscript. AS conceived the study, collected leaf samples, performed visual classification and helped to draft the manuscript. MS, SJ, HA and JT supervised the project. All authors read and approved the final manuscript.

Electronic supplementary material

12859_2005_575_MOESM1_ESM.pdf

Additional File 1: Definition of the employed spot descriptors. Provides a definition of the employed spot descriptors used to assess the quality of each spot. (PDF 121 KB)

Additional File 2: A list of the employed POP2 slides. Provides a list of the employed POP2 slides. (TXT 499 bytes)

12859_2005_575_MOESM3_ESM.pdf

Additional File 3: Implementation details of the segmentation process. Provides in-depth information regarding the implementation of the seeded region growing (SRG) algorithm. (PDF 161 KB)

12859_2005_575_MOESM4_ESM.zip

Additional File 4: The processed POP2 training data. Provides the processed POP2 training data set, which contains per-spot values of all descriptors employed here. (ZIP 7 MB)

12859_2005_575_MOESM5_ESM.zip

Additional File 5: The processed POP2 test data. Provides the processed POP2 test data set, which contains per-spot values of all descriptors employed here. (ZIP 7 MB)

12859_2005_575_MOESM6_ESM.zip

Additional File 6: Manual quality assessments. Provides a summary of the quality assessments as classified by the three experienced microarray users. (ZIP 193 KB)

12859_2005_575_MOESM7_ESM.pdf

Additional File 7: Visual representations of the sub-classes of bad spots. Provides images of typical examples of the 4 main sub-classes of bad spots. (PDF 50 KB)

12859_2005_575_MOESM8_ESM.pdf

Additional File 8: Details of the PLS-DA model. Provides details and statistics from the utilized PLS-DA model. (PDF 132 KB)

12859_2005_575_MOESM9_ESM.zip

Additional File 9: The designed subset of the POP2 training data. Provides the processed and filtered POP2 training data set, containing only the spots from the not bad, FI and BI classes which were selected according to the D-optimal design. (ZIP 206 KB)

12859_2005_575_MOESM10_ESM.pdf

Additional File 10: Description of the utilized D-optimal design. Provides information regarding generation of the D-optimal design used in to select subsets of the three classes. (PDF 75 KB)

Authors’ original submitted files for images

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://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

Bylesjö, M., Eriksson, D., Sjödin, A. et al. MASQOT: a method for cDNA microarray spot quality control. BMC Bioinformatics 6, 250 (2005). https://doi.org/10.1186/1471-2105-6-250

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1471-2105-6-250

Keywords