- Methodology article
- Open Access

# Improved analysis of bacterial CGH data beyond the log-ratio paradigm

- Lars Snipen
^{1}Email author, - Otto L Nyquist
^{2}, - Margrete Solheim
^{2}, - Ågot Aakra
^{2}and - Ingolf F Nes
^{2}

**10**:91

https://doi.org/10.1186/1471-2105-10-91

© Snipen et al; licensee BioMed Central Ltd. 2009

**Received:**12 November 2008**Accepted:**19 March 2009**Published:**19 March 2009

## Abstract

### Background

Existing methods for analyzing bacterial CGH data from two-color arrays are based on log-ratios only, a paradigm inherited from expression studies. We propose an alternative approach, where microarray signals are used in a different way and sequence identity is predicted using a supervised learning approach.

### Results

A data set containing 32 hybridizations of sequenced versus sequenced genomes have been used to test and compare methods. A ROC-analysis has been performed to illustrate the ability to rank probes with respect to Present/Absent calls. Classification into Present and Absent is compared with that of a gaussian mixture model.

### Conclusion

The results indicate our proposed method is an improvement of existing methods with respect to ranking and classification of probes, especially for multi-genome arrays.

## Keywords

- Positive Predict Value
- Negative Predict Value
- Comparative Genomic Hybridization
- Area Under Curve
- Probe Sequence

## Background

Microarray based comparative genomic hybridizations (CGH) is a tool for rapid investigation of the genetic content of bacteria. The technique is used for comparative genomic studies as well as screening for virulence factors or other genomic features of interest in a population [1–3]. The basic idea behind the technology is to construct microarrays from sequenced and annotated genomes, and then hybridize genomic DNA from other sources to these arrays to detect similarities and differences in genomic content. For two-color arrays DNA from some sampled genome is labeled and hybridized against labeled DNA from a reference. This reference is typically genomic DNA from one or several fully sequenced genomes, usually those from which the array was constructed.

The results obtained from such experiments can be seen as projections of the genomes in question onto the sequence space spanned by the microarray probe sequences. This probe space may vary in size, representing only a set of selected genomic features all the way up to pan-genomes. Probes may be short or long oligonucleotides, or PCR products, and we will in this paper only consider cases where the probe sequences are known exactly.

The data from these experiments are qualitatively different from those obtained in gene expression studies, where signal intensities must be seen as a continuum due to the dynamic abundance of mRNA. In bacterial CGH (bCGH) differences in signal intensities are predominately due to differences in sequence composition, copy number abberations are few and give smaller signal fluctuations. For this reason bCGH signals tend to behave more like a categorical variable with two possible outcomes, usually denoted Present and Absent. A strong signal, corresponding to Present, means the corresponding probe sequence is found, with sufficient similarity to yield hybridization, in the investigated genome. A weak signal means a too small part of the probe sequence is found in the genome to give hybridization, and the probe is called Absent.

Some methods to analyze bCGH data of this type have been proposed, and some of them are reviewed and tested in a recent publication by [4]. Most of these methods base their results on the log-ratio of signals, which is a standard adopted from the analysis of expression data. We will in this paper propose a new strategy for analyzing bCGH data, that does not rely on log-ratios, which we believe is a misleading paradigm for this type of data. Also, some previously proposed methods utilizing more than just log-ratios, like [5–7], are all unsupervised methods, not taking into account the sequence information from the reference genomes. In our approach this information is also included to aid the analysis. In the analysis of two-color microarray data demonstrated in this paper, we treat the array signals separately, almost as two single-color arrays, hence the method could easily be used for data from this technology as well. We test our method on data from *S. aureus* and *E. faecalis*, and compare our results to those achieved by other effective methods.

## Methods

### Sequence identity

Any bCGH experiment starts by performing alignments of every array probe sequence against the fully sequenced reference genomes to establish which probes are present and absent in these genomes. We use the term R-genome for a reference genome. We define the *identity* between probe and an R-genome as the number of identical bases in the best local alignment between them divided by the probe length to obtain a value between 0 and 1. We call this quantity *Rb*_{
ij
}for probe *i* against the R-genome in hybridization *j*. If *Rb*_{
ij
}= 1 it means an exact copy of the probe sequence is found in the R-genome, while if probe *i* has no significant hits in R-genome *j* we set *Rb*_{
ij
}= 0 even if all probes will of course have some very short subsequences in common with any genome.

The categorical response Present or Absent is coded as 1 or 0, respectively. This require, however, that intermediate *Rb*-values must be rounded to either 1 or 0, i.e. we need some *a priori* threshold that specify the sequence identity needed to be Present. Our analysis approach does not require categorical responses, and intermediate *Rb*-values can be used as is. However, if the ultimate goal is to classify between Present and Absent, the analysis is usually favored by having only 1 and 0 as responses from the start.

A sampled, un-sequenced, genome we call a sample-genome or S-genome. Corresponding to *Rb*_{
ij
}for the R-genome, we also have a similar *Sb*_{
ij
}for the S-genome. The motivation behind the entire bCGH experiment is to say something about this *Sb*_{
ij
}, i.e. the sequence identity between probe *i* and the S-genome in hybridization *j*.

### Preprocessing

For each array in the experiment, we assume background correction and within-array normalization has been done. We have employed standard methods in the LIMMA package [8] in R [9], available from the Bioconductor [10]. Normalization of CGH arrays has recently been discussed by [11] and [12], and nothing in our downstream analysis prevent the use of these or other approaches.

The flagging of low quality spots should be done very careful for bCGH analyses. In standard procedures for expression data, spots with low signals are removed. For bCGH data these spots turn out informative, because they span the range of array signals. Especially negative control probes, e.g. spots with alien or no DNA, are important since they carry information about which signals to expect when no hybridization takes place.

On microbial arrays probes are usually spotted multiple times (replicates). We will only consider the median value of these replicates on each array, but the number of replicates for each probe is kept as a weight in the final prediction, i.e. probes with more replicates have larger impact. Let
and
be the median preprocessed log-transformed signals from the R- and S-genome channel for probe *i* in hybridization *j*.

*I*

_{j 0}= {

*i*|

*Rb*

_{ ij }< 0.1}, i.e. the set of probes with R-genome sequence identity less than 0.1. Also, let

*I*

_{j 1}= {

*i*|

*Rb*

_{ ij }> 0.9}. Let

*Ra*

_{j 0}and

*Ra*

_{j 1}be the median of the

*Ra*

_{ ij }values for the probes in

*I*

_{j 0}and

*I*

_{j 1}, respectively. Then the between-array normalized R-signal is

*Sa*

_{ ij }. Notice that this procedure requires a significant number of probes to have low (less than 0.1) sequence identity with the R-genome, i.e. negative control probes are essential here. The effect of this normalization can be seen in Figure 1.

### Probe bias

Given a sequence identity *Rb*_{
ij
}, the corresponding array signals *Ra*_{
ij
}will in general correlate in a positive way, i.e. stronger sequence identity yields stronger array signal, and a similar relation we assume also holds between *Sa*_{
ij
}and the the unknown *Sb*_{
ij
}. However, probes with similar *Rb*-value may show consistently different *Ra*-values. This reflects a variable signal potential for the different probes due to sequence composition and/or bias during construction of the arrays. We refer to this as the probe bias. The same probe bias we assume is also present in the relation between *Sb*_{
ij
}and *Sa*_{
ij
}.

The *Rb*-values take on *L* discrete values between 0 and 1, and consider subsets of probes with similar *Rb*-value, i.e. ℐ_{
l
}= {*i*|*Rb*_{
ij
}= *l*} for *l* = 0,...,1. We assume for all *i* ∈ ℐ_{
l
}and hybridization *j* the linear model

*Ra*_{
ij
}= *μ*_{
lj
}+ *B*_{
ij
}

where *μ*_{
lj
}is the unconditional expected array signal at *Rb*-value *l* and *B*_{
ij
}is the probe bias for probe *i* in hybridization *j*. From this we get estimates of the probe bias for each hybridization
.

*i*we can get a pooled estimate of probe bias by averaging over the

*J*hybridizations, i.e.

*ω*∈ [0,1] and use as the final estimate of probe bias

Choosing *ω* close to 1 means information is 'borrowed' across hybridizations.

### Predicting sequence identity

The basic idea is, for each array, to fit a function that describes how *Rb*-values depend on bias-corrected *Ra*-values, and then use the same function to predict *Sb*-values from bias-corrected *Sa*-values.

We will make few assumptions about the shape of the function *f*_{
j
}, but we will require it to be monotonously increasing, since an increased array signal should always indicate stronger sequence identity.

*f*

_{ j }by a weighted running mean, where probes are weighted by their number of within-array replicates. For notational simplicity, let

*x*

_{ ij }=

*Ra*

_{ ij }- . The range of the function is divided into

*N*equally spaced

*knots*,

*x*

_{1},...,

*x*

_{ N }, and let

*D*be the width between two knots. For knot

*n*, let

*C*

_{ n }be the data subset {

*x*

_{ ij },

*Rb*

_{ ij }} whose value of

*x*

_{ ij }falls within

*x*

_{ n }± 3

*D*/2. Finding

*f*

_{ j }(

*x*

_{1}),...,

*f*

_{ j }(

*x*

_{ N }) leads to the constrained optimization problem

This problem can be solved by first computing the unconstrained optimum (weighted running mean), and then resolving the violated constraints in a recursive way. If the initial estimate of *f*_{
j
}(*x*_{n+1}) is smaller than that of *f*_{
j
}(*x*_{
n
}), both are replaced by the weighted average of them, weighted by the number of data points behind each initial estimate. This may again violate the constraints on the estimates of *f*_{
j
}(*x*_{n+2}) and/or *f*_{
j
}(*x*_{n-1}), and hence the recursion.

*f*

_{ j }(

*x*

_{1}),...,

*f*

_{ j }(

*x*

_{ N }) the estimated function value at any point within the range is found by linear interpolation between the knots. Let denote this estimated function for array

*j*. Figure 2 illustrates how fits a typical data set.

*j*, the prediction of the sequence identity for probe

*i*is now given as

It is not uncommon to repeat experiments, i.e. hybridize the same S-genome to several arrays. In this case it is natural to first analyze each array separately, obtain predictions from each array, and in the end average these
for each S-genome. A description of uncertainty in the prediction is best achieved by constructing a confidence interval for *Sb*_{
ij
}. Since this variable is trapped between 0 and 1 it seems reasonable to avoid inference based on specific distributions, and instead rely on some non-parametric approach. In case of a categorical response (Present/Absent), majority vote should be used instead of average, and statements concerning uncertainty should be put forward as some estimate of posterior probability of Present. The proportion of Present-votes for each probe is the maximum likelihood estimate of this probability, assuming the repeated experiments are independent.

### Data

In order to test methods we performed bCGH experiments using only sequenced genomes, i.e. the *Sb*-values are, contrary to a real situation, all known. Two different arrays were used, one representing 6 genomes of *Staphylococcus aureus* available from J. Craig Venter Institute [13] (JCVI), and one representing the genome of *Enterococcus faecalis* strain V583. In both cases probes are 70-mer oligonucleotides. The *S. aureus* array contains 5057 different probes spotted six times, where 4515 are ordinary probes representing genomes, and the remaining 542 negative control probes include various alien DNA and the 'empty probe' (no DNA). The *E. faecalis* array contains 3218 probes representing genes in the genome of V583, 10 probes representing the enterococcal pathogenicity island of strain MMH594 and 15 negative controls, giving a total of 3243 probes, spotted three times each.

*S. auerus*strains COL, N315, Mu50, NCTC8325 and RF122 and

*E. faecalis*strains V583 and OG1RF, whose genome sequences are available at NCBI [14]. For the

*S. aureus*experiments seven different pairs of genomes were selected for hybridization, and for each pair a dye-swap was performed. For each of these 14 hybridizations both genomes involved can play the role as R-genome and S-genome, hence there are altogether 28 different

*S. aureus*data sets where we can compare predicted and true sequence identity. Two hybridizations of V583 versus OG1RF were conducted (dye swap), and again both genomes can play the role as R-genome and S-genome, giving 4 additional

*E. faecalis*data sets. In order to compare our method against other methods we use a categorical response, i.e. each probe is classified as Present (1) or Absent (0). This means we have assigned a threshold to the

*Rb*- and

*Sb*-values in order to round each value to 1 or 0. We have used the threshold 0.7 (70% identity), i.e. an

*Rb*- or

*Sb*-value above 0.7 corresponds to Present and is rounded to 1 and values below 0.7 is rounded to 0. The threshold is chosen on the basis of the histogram in Figure 3. The

*S. aureus*arrays contain probes representing genes in 6 different strains. By BLASTing the probe sequences against the genome sequences of these strains, the identities distribute as indicated in Figure 3. Thus, it seems that probes matching with approximately 70% identity or more are considered Present in the genome by JCVI who designed the arrays. This also corresponds well with our experience regarding the degree of match giving hybridizations. This threshold will in general depend on array design and hybridization conditions, and a proper value must be decided upon for each experiment separately. Our method is independent of this choice as long as it is a reasonable value for the experiments analyzed. Table 1 show the percent of truly Present/Absent probes in each of the genomes using our probe set and threshold.

Genomes and microarrays

Genome | Size (Mb) | Present | Absent |
---|---|---|---|

| 2.81 | 74% | 26% |

| 2.84 | 74% | 26% |

| 2.90 | 76% | 24% |

| 2.82 | 74% | 26% |

| 2.74 | 70% | 30% |

| 3.36 | 99% | 1% |

| 2.73 | 73% | 27% |

As previously mentioned, we advocate a weak flagging of array spots during the preprocessing of the data. This means only truly damaged spots should be flagged, and spots with weak signals or negative controls, should be part of the data set through the entire analysis. When comparing our proposed method against other approaches, we used both 'hard' and 'weak' flagging of spots to illustrate the differences between these strategies. By 'hard' flagging we mean removing all negative controls as well as all spots flagged by the image analysis software, i.e. in our case all spots with negative flag value from GenePix. In the 'weak' flagging only manually discarded spots were removed, i.e. only spots with flag value -100 from GenePix.

## Results

Our proposed method predicts probe sequence similarity to a sampled genome based on a biased-corrected array signal. Based on observed array signal and probe sequence similarities to the reference genome, we estimate a probe bias for each probe. Then, correcting for this probe bias, we fit a non-parametric function describing the relation between array signal and probe sequence similarity for the reference genome. Finally, we use this function to predict probe sequence similarity from observed array signals for the sampled genome. If a categorical response (Present or Absent) is desired this is coded as Present = 1 and Absent = 0. Comparison to other approaches are here made on data sets where true sequence similarities (Present/Absent status) are known.

### ROC-analysis

*E. faecalis*data sets completely lacked absent probes, and hence no AUC-values could be computed for these data sets. Thus, only 2 of the 4

*E. faecalis*data sets were included in the ROC-analysis. Figure 4 shows the AUC-values for both ranking criteria. An AUC-value of 1.0 means perfect separation of classes, while a value close to 0.5 means ranking is completely random, i.e. both classes are mixed in the ranked list. In this analysis we used the weight

*ω*= 0.75 to compute the probe bias effect. Other choices of these weights produced very similar AUC-values, and did not alter the big picture.

### Effect of bias weight *ω*

*ω*from (4). A weight close to 1 means information is borrowed between arrays when it comes to estimating the probe bias. To get an impression of the effect of this constant, we varied it systematically over the interval [0, 1], and for each weight classified all probes in all data sets. For each data set we computed the classification error as the geometric average [16]. This is the square root of the product of sensitivity (probability of classifying as Present when truly Present) and specificity (probability of classifying as Absent when truly Absent). Figure 5 illustrate how the geometric average varies for different choices of

*ω*over the

*S. aureus*and

*E. faecalis*data sets.

### Comparing classification results

Classification results

Array | Flagging | # probes | Method | Sens. | Spec. | PPV | NPV |
---|---|---|---|---|---|---|---|

| Weak | 4515 | BCSP | 0.968 | 0.938 | 0.987 | 0.861 |

Mix.mod. | 0.945 | 0.901 | 0.979 | 0.776 | |||

p-value | 4.1·10 | 7.4·10 | 1.7·10 | 2.5·10 | |||

Hard | 3539 | BCSP | 0.977 | 0.864 | 0.989 | 0.754 | |

Mix.mod. | 0.960 | 0.825 | 0.985 | 0.607 | |||

p-value | 3.2·10 | 0.02 | 0.003 | 3.7·10 | |||

| Weak | 3228 | BCSP | 0.989 | 0.953 | 0.983 | 0.560 |

Mix.mod. | 0.946 | 0.890 | 0.989 | 0.458 | |||

Hard | 3145 | BCSP | 0.992 | 0.0 | 0.869 | 0.0 | |

Mix.mod. | 0.955 | 0.474 | 0.990 | 0.409 |

### Prediction error

*Rb*-values should not be rounded to 0 or 1, but used as is, as illustrated in the left panel of Figure 2. However, since the large majority of probes are either completely present or absent, predicting an intermediate sequence identity is usually a sign of uncertainty of the probes actual status. This is reflected in Figure 6, where we have indicated the average absolute error | -

*Sb*

_{ ij }| for the different predicted values .

## Discussion

There is at present no standard approach for analyzing bacterial CGH data, and the methods reviewed by [4] are only a selection of approaches employed in recent bCGH-publications, e.g. see [18] and [19]. Common to the large majority of these methods is the use of the log-ratio for ranking and classifying probes. In our notation it means sequence identity *Sb* is predicted from array signal *Sa* - *Ra*. This is a paradigm inherited from the analysis of expression data. However, for bCGH data it is actually possible to test how informative this quantity is, since we can perform experiments with one sequenced genome against another. This was done by [6] and [7], and from both publications we may conclude that combining *Sa* and *Ra* in other ways than just subtracting one from the other, is superior. In this paper we have a much larger data set, and the results from Figure 3 clearly show the same picture. For both weak and hard flagging ranking by the bias-corrected S-signal produce larger AUC values than ranking by log-ratio. Hence, we can extract more information from array-signals than just the log-ratios.

In our present approach we have also utilized the sequence information *Rb* directly in the prediction of *Sb*. This seems like a new idea, even if [20] has utilized sequence information in the analysis of single-channel CGH data. When predicting the sequence identity of the S-genome, *Sb*, we first consider how sequence identity *Rb* and array signal *Ra* relates to each other, and then use this to predict *Sb* from *Sa*. The reason a rather obvious approach like this has not be tried out long ago must be due to the tunnel-vision imposed by the log-ratio paradigm. In our approach we treat signals from dual-dye arrays almost as if they were from two single channel arrays, and then use the signal-genotype relation on one array to predict the signal-genotype relation on the other. For this reason the implementation of our method for single channel arrays is straightforward. The only requirement is that for each sample-genome investigated there is also a set of reference signals, i.e. at least one array must be used to hybridize an already sequenced genome to obtain these reference signals.

An argument for using log-ratios is that probe signal biases are canceled. Since we do not use log-ratios, we compensate for this effect by estimating a probe bias from Eq. 4 and then subtract it in Eq. 6. Figure 4 indicates that the weight *ω* should be large, somewhere between 0.7 and 1.0. However, the differences in geometric average are small for various choices of *ω*, and even at *ω* = 0 it is well above 0.9. The values at *ω* = 0 also indicates the precision we get for analyzing a single array, because here we do not borrow any information across arrays. Hence, these results indicates only a small gain in performing a batch of hybridizations, and analyze all arrays together compared to doing it array-by-array.

In Table 2 the results for classification in all 32 data sets are displayed. For the 28 *S. aureus* data sets the picture is clear: Our proposed method, denoted BCSP, performs better than the log-ratio-based mixture model, which is the 'winner' in [4]. For all four criteria sensitivity, specificity, positive predicted value and negative predicted value, the BCSP method gives significant improvement to the mixture model method (small p-values). Noticeable is also the difference between weak and hard flagging. By hard flagging around 1000 ordinary probes are removed from the data set (in addition to all negative control probes), while with weak flagging none are removed. Sensitivity is always improved by hard flagging, but specificity is poorer. The latter means Absent probes become more difficult to detect after hard flagging. This is natural for the BCSP method, since the informative negative controls are no longer available. In general, hard flagging means there are fewer data with small *Ra* and *Rb*-values, and the shape of the functions displayed in Figure 2 become more uncertain and difficult to estimate. Given the excellent results for weak flagging, we can think of no good reason to throw away a large proportion of the probes in a hard flagging procedure. For the *E. faecalis* data the results are more unclear. For weak flagging the BCSP method gives better sensitivity, specificity and NPV, but slightly poorer PPV. No differences are significant, basically because there are only 4 data sets. For hard flagging BCSP produce absolute no specificity, i.e. no truly Absent probes are classified as absent! This illustrates the dramatic effect of losing all information about negative controls and other probes with *Rb*-value equal to 0. Also the mixture model behaves poorly for hard flagging, and again this support a weak flagging strategy.

A difference between the *S. aureus* and *E. faecalis* data is that the *S. aureus* array contain probes representing features in several genomes, a multi-genome array, while the *E. faecalis* array contain little more than what is found in the strain V583. Hence, in the *S. aureus* case there is always a large number of probes that should not hybridize against a specific *S. aureus* genome used for reference. This situation is ideal for our proposed method because there will always be a good balance between probes with small and large *Rb*-values. In a recent publication [21] argues that for multi-genome arrays a mixture of all genomes represented on the array should be used as the reference DNA pool. Their conclusion is based on an analysis of log-ratios. For our supervised learning approach, this strategy should clearly be avoided. If you want to discriminate between Present and Absent in the S-genome channel, you must make certain you have data that show the difference between Present and Absent in the R-genome channel as well. Hence, there should always be a substantial amount of probes against which a reference does not hybridize. Figure 6 illustrate that reliable predictions of sequence identity can only be given for very low or very high identities, i.e. for probes who are either more or less completely Absent or Present. Thus, even if our proposed method opens up the possibility to use and predict any sequence identity, intermediate identities always introduce difficulties. Thus, predicting an identity around 0.5 can be seen as an indication of a large uncertainty.

## Conclusion

We have proposed a method for analyzing bacterial CGH data that seems to be a significant improvement compared to any log-ratio based approach, as indicated by the ROC-analysis. For actual classification we also tend to get improved results compared to the log-ratio based mixture model approach, which was the 'winner' in the survey of [4]. Instead of forming log-ratios, we employ a supervised learning approach where sequence identities are predicted from bias-corrected array signals in each channel separately. The proposed method require a substantial number of probes with little or no sequence identity to the reference genome used in the hybridization. Thus, the method is particulary well suited for data from multi-genome arrays.

## Availability

R code for handling bCGH data using this method, as well as other approaches, is freely available from the corresponding author.

## Declarations

### Acknowledgements

OLN was financially supported by a research grant from the Norwegian University of Life Sciences. MS was financially supported by the European Union 6th Framework Programme "Approaches to Control multi-resistant Enterococci: Studies on molecular ecology, horizontal gene transfer, fitness and prevention". ÅA were supported by grants from the Research Council of Norway. *S. aureus* micorarrays were kindly provided by the Pathogen Functional Genomics Resource Center (PFGRC) at the J. Craig Venter Institute (JCVI), Rockville, MD, USA. We acknowledge Aksel Flack, The Norwegian Microarray Consortium, Oslo, for printing of the *E. faecalis* microarray slides.

## Authors’ Affiliations

## References

- Dorrell N, Champion OL, Wren BW: Application of DNA Microarrays for Comparative and Evolutionary Genomics. Methods in Microbiology. 2002, 33: 121-136.View ArticleGoogle Scholar
- Lindsay JA, Moore CE, Day NP, Peacock SJ, Witney AA, Stabler RA, Husain PDSE, Butcher JH: Microarrays Reveal that Each of the Ten Dominant Lineages of Staphylococcus aureus Has a Unique Combination of Surface-Associated and Regulatory Genes. Journal of Bacteriology. 2006, 188 (2): 669-676.PubMed CentralView ArticlePubMedGoogle Scholar
- Willenbrock H, Petersen A, Sekse C, Kiil K, Wasteson Y, Ussery DW: Design of a Seven-Genome Escherichia coli Microarray for Comparative Genomic Profiling. Journal of Bacteriology. 2006, 188 (22):Google Scholar
- Carter B, Wu G, Woodward MJ, Anjum MF: A process for analysis of microarray comparative genomics hybridisation studies for bacterial genomes. BMC Genomics. 2008, 9 (53):Google Scholar
- Repsilber D, Mira A, Lindroos H, Andersson S, Ziegler A: Data rotation improves genomotyping efficiency. Biometrical Journal. 2005, 47 (4): 585-598.View ArticlePubMedGoogle Scholar
- Snipen L, Repsilber D, Nyquist L, Ziegler A, Aakra Å, Aastveit A: Detection of divergent genes in microbial aCGH experiments. BMC Bioinformatics. 2006, 7 (181):Google Scholar
- Feten G, Almøy T, Snipen L, Aakra Å, Nyquist OL, Aastveit AH: Mixture Models as a Method to Find Present and Divergent Genes in Comparative Genomic Hybridization Studies on Bacteria. Biometrical journal. 2007, 49 (2): 242-258.View ArticlePubMedGoogle Scholar
- Smyth GK, Speed TP: Normalization of cDNA microarray data. Methods. 2003, 31: 265-273.View ArticlePubMedGoogle Scholar
- The R project. [http://www.r-project.org/]
- The Bioconductor. [http://www.bioconductor.org/]
- van Hijum SAFT, Baerends RJS, Zomer AL, Karsens HA, Martin-Requena V, Trelles O, Kok J, Kuipers OP: Supervised Lowess normalization of comparative genome hybridization data – application to lactococcal strain comparisons. BMC Bioinformatics. 2008, 9: 93-PubMed CentralView ArticlePubMedGoogle Scholar
- Staaf J, Jonsson G, Ringner M, Vallon-Christersson J: Normalization of array-CGH data: influence of copy number imbalances. BMC Genomics. 2007, 8: 382-PubMed CentralView ArticlePubMedGoogle Scholar
- The J. Craig Venter Institute. [http://www.jcvi.org/]
- GenBank. [http://www.ncbi.nlm.nih.gov/Genomes/]
- Hanley JA, McNeil BJ: The meaning and use of the area under a receiver operating characteristic (ROC) curve. Radiology. 1982, 143: 29-36.View ArticlePubMedGoogle Scholar
- Kubat M, Holte R, Matwin S: Machine learning for the detection of oil spills in satellite radar images. Machine Learning. 1998, 30: 195-215.View ArticleGoogle Scholar
- McLachlan GJ, Peel D: Finite Mixture Models. 2000, New York: John Wiley & SonsView ArticleGoogle Scholar
- da Silva VS, Shida CS, Rodrigues FB, Ribeiro DCD, de Souza AA, Coletta-Fiho HD, Machada MA, Nunes LR, de Oliveira RC: Comparative genomic characterization of citrus-associated Xylella fastidiosa strains. BMC Genomics. 2007, 8 (474):Google Scholar
- Jayapal KP, Lian W, Glod F, Sherman DH, Hu WS: Comparative genomic hybridizations reveal absence of large Streptomyces coelicolor genomic islands in Streptomyces lividans. BMC Genomics. 2007, 8 (229):Google Scholar
- Schuster EF, Blanc E, Partridge L, Thornton J: Correcting for sequence biases in present/absent calls. Genome Biology. 2007, 8: R125-PubMed CentralView ArticlePubMedGoogle Scholar
- Pinto FR, Aguiar SI, Melo-Cristino J, Ramirez M: Optimal control and analysis of two-color genomotyping experiments using bacterial multistrain arrays. BMC Genomics. 2008, 9 (230):Google Scholar

## Copyright

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.