An integrated Bayesian analysis of LOH and copy number data
 Paola MV Rancoita^{1, 2, 3}Email author,
 Marcus Hutter^{4},
 Francesco Bertoni^{2} and
 Ivo Kwee^{1, 2}
DOI: 10.1186/1471210511321
© Rancoita et al; licensee BioMed Central Ltd. 2010
Received: 4 November 2009
Accepted: 15 June 2010
Published: 15 June 2010
Abstract
Background
Cancer and other disorders are due to genomic lesions. SNPmicroarrays are able to measure simultaneously both genotype and copy number (CN) at several Single Nucleotide Polymorphisms (SNPs) along the genome. CN is defined as the number of DNA copies, and the normal is two, since we have two copies of each chromosome. The genotype of a SNP is the status given by the nucleotides (alleles) which are present on the two copies of DNA. It is defined homozygous or heterozygous if the two alleles are the same or if they differ, respectively. Loss of heterozygosity (LOH) is the loss of the heterozygous status due to genomic events.
Combining CN and LOH data, it is possible to better identify different types of genomic aberrations. For example, a long sequence of homozygous SNPs might be caused by either the physical loss of one copy or a uniparental disomy event (UPD), i.e. each SNP has two identical nucleotides both derived from only one parent. In this situation, the knowledge of the CN can help in distinguishing between these two events.
Results
To better identify genomic aberrations, we propose a method (called gBPCR) which infers the type of aberration occurred, taking into account all the possible influence in the microarray detection of the homozygosity status of the SNPs, resulting from an altered CN level. Namely, we model the distributions of the detected genotype, given a specific genomic alteration and we estimate the parameters involved on public reference datasets. The estimation is performed similarly to the modified Bayesian Piecewise Constant Regression, but with improved estimators for the detection of the breakpoints.
Using artificial and real data, we evaluate the quality of the estimation of gBPCR and we also show that it outperforms other wellknown methods for LOH estimation.
Conclusions
We propose a method (gBPCR) for the estimation of both LOH and CN aberrations, improving their estimation by integrating both types of data and accounting for their relationships. Moreover, gBPCR performed very well in comparison with other methods for LOH estimation and the estimated CN lesions on real data have been validated with another technique.
Background
Although most of the human genome is identical among individuals, there are about 10 million single nucleotide polymorphisms (SNPs) which distinguish us [1]. SNPs are single basepair loci where the nucleotides can assume two possible values (called alleles) among the four bases (thymine, adenine, cytosine, guanine). In general, since we have two copies of each chromosome, the genotype at any SNP can be: AA, BB or AB, where A and B represent the two alleles. Moreover, a SNP can be classified as homozygous (i.e., AA or BB) or heterozygous (i.e., AB), whether or not its genotype consists of two equal alleles. Cancer and several human diseases are caused by genomic aberrations, which can affect the homozygous status and/or the DNA copy number (the normal copy number, CN, is two since we have two copies of each chromosome, except for the chromosomes X and Y). The former type of aberrations is often displayed by unusual long stretches of homozygous SNPs, called loss of heterozygosity (LOH) region. The latter type of aberrations consists in genomic regions with DNA copy number different from two.
In general, LOH can arise by several mechanisms, such as deletions and germline or somatic recombinations. When the LOH occurs without a change in copy number, it is referred to as copyneutral LOH (or sometimes run of homozygosity, ROH). In the past, copyneutral LOH regions were usually explained as a consequence of an uniparental disomy event (UPD), see [2]. Recently, long homozygous segments have also been detected in genomes of normal individuals, supporting the hypothesis that some copyneutral LOH segments might represent autozygosity (see, for example, [3–5]). In the literature, it has been shown a relationship between some tumors and both types of aberrant events (see, for example, [6–8]).
Uniparental disomy (UPD) occurs when both homologues of a part of a chromosome are inherited from only one parent. It can be divided in: uniparental isodisomy (when the two copies are two replicates of one homologue of one parent) and the uniparental heterodisomy (when both homologues are inherited from the same parent). Because of meiotic recombination, a mixture of both events is also possible, and similar events can also happen during the mitosis. Moreover, in cancer cells, the uniparental isodisomy can also occur when a homologue of a part of a chromosome is lost and the remaining homologue is duplicated. The autozygosity describes a situation where the homologues are identical by descendent (IBD), because they are inherited from a common ancestor. Inbreeding is usually uncommon because of laws and social conventions, but it does occur in small isolated populations.
SNP microarrays are able to measure simultaneously both the DNA copy number and the genotype at each SNP position considered [9]. We call LOH data the homozygous status of the SNPs deduced from the genotyping data. By integrating the information given by both LOH and copy number data, we can better identify several types of lesions of the genome (regarding combinations of both DNA copy number and LOH aberrations). For example, when one copy of a chromosomal segment is deleted, we usually detect a long stretch of homozygous SNPs (since the genotype calling algorithm is unable to distinguish between the presence of only one copy and the presence of two equal copies), but the same homozygous status can also occur for other reasons, such as uniparental disomy. In this situation, the knowledge of both types of data can lead to the correct interpretation of the phenomenon, while with only one type of data it would not be possible. Another example is when an amplified genomic segment is present: if one of the two copies of the segment is highly amplified, then even the heterozygous SNPs will be likely detected as homozygous, because the DNA quantity of one allele is much higher than the other one. In this case again, the integration of both types of data is able to better identify the dosage of the DNA aberration.
Many methods have been developed for the estimation of the copy number profile (see, for example, [10–14]) and others for the discovery of LOH regions from the genotyping data, without distinguishing if they are caused by either the loss of one copy or other genomic events like uniparental disomy or autozygosity (see, for example, [15, 16]). To the best of our knowledge, only one method integrates these two types of data for the estimation of both copy number aberrations and copyneutral LOH regions and it uses a hidden Markov model (HMM) [17]. Other statistical procedures use the information regarding both the total and the allelic copy number to infer these kind of lesions (for example, [18–24]) and some of these algorithms are available only for the analysis of data coming from Illumina Beadarrays. Finally, in [25] the authors describe an HMM with the same purpose, which employs the allelic copy number data from a tumor sample and the genotyping data from the matched normal sample.
Here, we propose a method which estimates the copy number changes and the copyneutral LOH regions at the same time, using both LOH and DNA copy number data. The estimation procedure consists of a Bayesian piecewise constant regression, thus we call our algorithm genomic Bayesian Piecewise Constant Regression (gBPCR). Our model is more general than [17], since the latter cannot be applied to data, whose DNA sample come from a mixture of cell populations (which is usually the case for samples of patients affected by cancer). Moreover, the algorithm in [17] needs the specification of some parameters by the user and is sensitive to their values.
Our method was implemented in R and is freely available at http://www.idsia.ch/~paola/gBPCR/ or in Additional file 1. Furthermore, an R package will be soon available.
Methods
Because of the complexity of the biological model, we first describe a preliminary simplified model (called Model 1), which only estimates the copy number events exploiting the relationship between copy number and LOH data. Therefore, it does not identify copyneutral LOH regions (called IBD/UPD regions), which are due to events like uniparental disomy, and it does not distinguish the normal regions from the gained one (because we suppose that the capability of detection of the homozygous status is the same in these two types of regions). In the subsequent subsections, we add to the model the detection of copyneutral LOH regions (Model 2) and of gained ones (Model 3). Therefore, the explanation is structured in the following way:

Model 1: relationship between LOH and copy number data to detect copy number changes (apart from the gained regions);

Model 2: addition to Model 1 of the IBD/UPD region detection (i.e. determination of copyneutral LOH regions);

Model 3: addition to Model 2 of the gained region detection.
Each of the three models is contained in the subsequent. The final model (Model 3) represents our algorithm for the estimation of both copy number changes and copyneutral LOH regions and we call it genomic Bayesian Piecewise Constant Regression (gBPCR).
Model 1: relationship between LOH and copy number data
Although in nature the copy number is an integer, the raw copy number values detected by the microarray are usually continuous values, due to technical procedures. Moreover, the samples often contain also a percentage of normal cells.
It is common practice to treat copy number data in a log_{2}ratio scale (where the ratio is defined with respect to a normal reference dataset) which makes the errors approximately normally distributed. Then, the copy number profile is estimated as a piecewise constant function (i.e. the genome is divided in regions of constant copy number), where the levels assume real values. For the purpose of our model, we estimate this profile by mBPCR, which is a Bayesian piecewise constant regression procedure [14]. It has been shown in [14] that this method outperformed wellknown other methods on several datasets.
Commonly, in biomedical/cancer research, after estimating the log_{2}ratio profile, the copy number aberrations are defined as those regions with values outside an interval around zero (notice that, in the log_{2}ratio scale, zero represents CN = 2, i.e. a normal copy number). Often, the interval is a statistical confidence interval computed on the basis of the samples of the whole dataset.
In Model 1, our aim is to classify better the copy number changes, trying to reduce the number of false positives, by exploiting the relationship between copy number and LOH data.
Mathematical model of the biology mechanism
The aim of Model 1 is to obtain a better estimation of the true underlying copy number events, using both the information given by copy number and LOH data. In a genomic region, a copy number event is defined as a particular class of copy number values. The definition of the categories into which the copy number values are divided will follow from the description of the LOH data.
For the purpose of better identifying the copy number events, we can consider two classes of SNP values: Heterozygosity (Het) and Homozygosity (Hom). Thus, the LOH data are deduced from the genotyping data, by mapping the genotypes AA and BB into Hom and the genotype AB into Het, for all SNPs. The genotype calling algorithm (e.g. BRLMM [26]) is unable to distinguish between a homozygosity due to the presence of two equal nucleotides or the one due to the loss or high amplification of one of them. Hence, the presence of heterozygosity can ensure that the copy number is normal or gained with a high probability, while the homozygosity can be due to different events. It follows that there are only four relevant classes of copy number events that can be distinguished by looking at the LOH data. Therefore, if we call the random variable which represents a copy number event at SNP i, it can assume only the following values:

= 2, when CN > 4 (amplification);

= 0, when 1 < CN ≤ 4 (normal or gain);

= 1, when CN = 1 (loss);

= 2, when CN = 0 (homozygous deletion).
The homozygous deletion corresponds to the loss of both copies of a genomic region. Ideally, the genotype calling algorithm should detect a NoCall genotype at the corresponding SNP position (i.e. it should not be able to identify the genotype of the SNP). Although not common, since cancer DNA samples usually contain a mixture of normal and tumor cells (with also different cancer cell subpopulations), the information given by the NoCall genotype can be used to better distinguish between a monoallelic deletion and a biallelic (homozygous) deletion.
Therefore, three different LOH variables are present in the model: the true homozygous status in normal cells (X^{ N }), the homozygous status in abnormal cells (X), which is the consequence of copy number changes (in Model 1 we do not consider other biological events), and the homozygous status detected by the genotype calling algorithm (Y). The components of the first two random vectors can assume only values in = {Het, Hom} and ^{*} = {∅, Het, Hom}, respectively, and we suppose that they are independently distributed as Bernoulli random variable. The components of Y can assume values in = {NoCall, Het, NHet} (NHet stands for "not heterozygous", since the genotype calling algorithm cannot distinguish between two equal nucleotides, i.e. homozygosity, and the loss of one copy).
Notations.
Het  heterozygous 

Hom  homozygous 
NHet  not heterozygous (is used when we cannot distinguish between two equal nucleotides, i.e. homozygosity, and the loss of one copy) 
 {Het, Hom} 
^{*}  {ø, Het, Hom} 
 {Het, NHet, NoCall} 
X ^{ N }  true genotypes in normal cells ( ∈ ) 
X  true genotypes in abnormal cells (X_{ i }∈ ^{*}) 
Y  genotypes detected by the genotype calling algorithm (Y_{ i }∈ ) 
Y ^{ cn }  raw copy number data 
 copy number events/aberrations 
 occurrence of copyneutral LOH (i.e. IBD/UPD event) 
 IBD/UPD & copy number aberrations ({ = w} = { , = u} for some w, z, u) 
cn  all copy number information (both raw data and estimated profile by mBPCR) 
p  vector of posterior probabilities to be a breakpoint (for all SNP positions) 
Nevertheless, the homozygous status of abnormal cells estimated by the genotype calling algorithm (Y_{ i }) is affected by several sources of errors.
Hypothesis of the model
The genome of cancer cells can be divided in subregions where the copy number is constant. Since we divided the copy number values in four classes (i.e. the copy number events), we can also consider regions with the same copy number event.
Let us consider a genomic region where the microarray measures the DNA copy number and the genotype at n SNP loci. Then, from the previous discussion, the vector of the copy number events at all positions can be seen as a piecewise constant function. This function consists of k_{0} intervals with the same copy number event and with boundaries so that , for all p = 1, ..., k_{0}. To estimate this function we use a Bayesian piecewise constant regression method, which determines the number of segments k_{0}, the boundaries ( ) and the copy number events .
where and are the domains of k and t, respectively (they will be defined later).
To specify the model (see Figure 1), we need to define the likelihood, i.e. the conditional distribution of Y, given and X^{ N }. To model it, we take into account all the variability that can affect the genotype detection, such as: the polymerase chain reaction (PCR) amplification, the presence of different cancer cell subpopulations or normal cells, and the amplification of only one copy. For example, the probabilities and are not zero, because of the error in the genotype detection even in case of a normal DNA sample. The probabilities and are related to the detection errors due to the presence of normal cells and/or different types of cancer cell subpopulations, or to PCR amplification errors, while is related to the errors that can be due to the amplification of only one allele. Also and account for the errors that can be due to the presence of cell subpopulations.
The set of conditional probabilities are considered as parameters of the model. To quantify them, we needed paired normalcancer samples, since they are related to the probability of detecting a certain homozygous status in a cancer cell, given the corresponding one in a normal cell of the same patient and under some copy number event. Therefore, to compute maximum likelihood estimates of these parameters, we used 13 samples from an available cancer dataset consisting of breast cancer cell lines [27, 28] (see Section S.1 in Additional file 2, for further explanations).
where = {1, ..., k_{max}} and is a subspace of such that t_{0} = 0, t_{ k }= n and t_{ q }∈ {1, ..., n  1} for all q = 1, ..., k  1, in an ordered way and without repetitions.
The prior probabilities of heterozygosity of the SNPs are the frequencies of heterozygosity computed on the samples of the matched race in the HapMap project [3, 29]. They are usually provided by the manufacturer in the documentation related to the microarray used. In Section "Results and Discussion", the microarray mostly employed is the GeneChip Human Mapping 250K NspI (Affymetrix, Santa Clara, CA, USA).
Z prior definition
The only prior that we have not yet defined is the one of Z. While the estimated levels of the log_{2}ratio profile are continuous variables, Z classifies the copy number as discrete events. Then, the major problem consists in mapping the continuous values into the discrete values of Z, i.e. in defining a partition of the log_{2}ratio values such that each interval corresponds to a particular copy number event.
where are, respectively, the estimated mean and variance of the normal distribution corresponding to CN = cn.
where cn represents all copy number information (both raw data and estimated profile by mBPCR) and M_{ p }is the random variable representing the log_{2}ratio value in the p^{ th }segment. From the mBPCR model, given cn, M_{ p }is normally distributed with mean and variance where ( , ) are the posterior mean and variance of M_{ p }estimated by mBPCR, respectively.
The estimation
Namely, corresponds to the positions which have the highest posterior probability to be a breakpoint. The main differences with respect to mBPCR are in the prior over K and in the estimation of K. Instead of using a uniform prior and an estimator which minimizes the posterior expected squared error, here we consider a prior similar to 1/k^{2} and an estimator which minimizes the 01 error, in order to reduce the false discovery rate (FDR) in case of few segments.
for p = 1, ..., , where and are any estimate of t^{0} and k_{0}, respectively. For the computation of all the posterior probabilities involved, we used dynamic programming as described in Section S.3 in Additional file 2.
Therefore, if the boundary estimator misses a clear boundary between and , then the probability at the denominator of Equation (5) could be zero and thus the level would not be estimated. The best way to prevent this event consists in using a good estimator for the boundaries.
Since commonly the vector of the posterior probabilities shows clearly the position of the breakpoints in correspondence to the maxima, we estimate the number of the segments and the breakpoints with the number of peaks and the locations of their maxima, respectively (see Section S.4 in Additional file 2). Essentially, the algorithm for the determination of the peaks, after applying a kernel method to reduce the noise of the function, uses two thresholds: one for the determination of the peaks (thr_{1}) and one for the definition of the values close to zero (thr_{2}). Therefore, we will denote the corresponding estimators by and .
In Section "Results and Discussion", we will consider several pairs of thresholds and we will apply the corresponding estimators to simulated data, in order to determine the best paired thresholds and to compare their performance with . We will also compare with , another boundary estimator described in [14].
Model 2: addition of the IBD/UPD region detection
LOH data are used in biology not only to better identify regions of loss and amplifications, but, especially, to detect regions of copyneutral LOH, which can be identified by unusual long stretches of homozygous SNPs, with normal copy number. In Section "Background", we explained that this type of aberrations can be a consequence of UPD (either uniparental isodisomy or heterodisomy) or autozygosity (IBD regions). From the description of these genomic events, it follows that the uniparental isodisomy and the IBD regions can be detected because they appear as a long sequence of homozygous SNPs with a low probability to occur, while the uniparental heterodisomy consists in a sequence of both homozygous and heterozygous SNPs as in a normal condition. Therefore, without the genotypes of the parents, from SNP data we can only detect the uniparental isodisomy (iUPD) and the IBD segments. In the following, we will consider only these two events, referring to them as IBD/UPD events.
Since an IBD/UPD event, by definition, only exists in regions of normal copy number (CN = 2), the only probabilities which are affected by the presence of this event are those involving {Z = 0}. Therefore, we define the following sets of conditional probabilities and , where the variable indicates if an IBD/UPD event occurred at SNP i (if it happened = 1, otherwise = 0). We can notice that, given { = 0, = 0}, the distribution of Y_{ i }is equal to the conditional distribution with respect to { = 0} in Model 1, since the latter was modeled with no possibility of an IBD/UPD event. Instead, in case of an IBD/UPD event, we do not need to condition with respect to , since, in case of a somatic iUPD event, the genotype of an iUPD region is independent of the homozygous status of the same region in a normal cell. Otherwise, in case of autozygosity or germ line iUPD, the genotypes of normal and abnormal cells are the same and it makes no sense to condition one to the other.
and we assume that P(U_{ p }= 1) =: p_{ upd }, for all p = 1, ..., .
Both and p_{ upd }are parameters of the model. For the maximum likelihood estimation of , we used 11 IBD/UPD regions previously found by us on 5 samples of patients with hairy cell leukemia [34] and on the Bcell lymphoma cell line KARPAS422. All regions were detected by dChip [16] and their width was between 3 Mb and 100 Mb (covering from 300 to 9800 SNPs), so that they were large enough to be really considered IBD/UPD regions (for further explanations, see Section S.1 in Additional file 2).
Values for the parameter p_{upd}
We expect the prior probability of an IBD/UPD event to be low. In order to estimate the order of magnitude of this parameter, we considered two studies on IBD regions: [6] and [3]. In the former, they considered as IBD regions only stretches of at least 50 homozygous SNPs (with at maximum 2% of heterozygous) longer than 4 Mb and the platform used was the Affymetrix GeneChip Human Mapping 50K Array. In the latter, a denser microarray was used and the stretches considered were longer than 1 Mb (with at least 50 probes) or longer than 3 Mb. Using the data of the former paper (only the normal samples), we estimated p_{ upd }≈ 1.7·10^{3}. Instead, with the data of the latter, we estimated p_{ upd }≈ 1.5·10^{3} considering all regions greater than 1 Mb, while p_{ upd }≈ 1.46·10^{4}, considering only the regions greater than 3 Mb. The differences in the estimated values are due to the different resolution of the technologies used (in fact, in the former the number of SNPs used was 58,960, while in the latter it was 3,107,620). Moreover, the probability depends on the minimum length allowed for these regions. The wider the regions are, the higher is the probability that the regions represent "abnormalities" and the lower becomes the probability of their occurrence (so that p_{ upd }is lower). Therefore, in the following applications (see Section "Results and Discussion"), we will use two values: p_{ upd }= 10^{3} and p_{ upd }= 10^{4}.
Another possible way to solve the problem could be to assign a prior distribution to p_{ upd }(for example, a uniform distribution over its range) and integrate it out in the equations of the model.
Model 3: addition of the gained region detection
In the following, Model 3 (which is the complete model) will be called genomic Bayesian Piecewise Constant Regression (gBPCR).
Adjustment of the parameters related to NoCall
The probabilities {P(Y_{ i }= NoCall = x, = w), x ∈ , w = 3, 2, 1, 0, 2} are related to the detection of NoCall s under some conditions. Generally, the presence of NoCall s is not only due to diffculties of the genotype calling algorithm in the detection of the genotype (technical noise) but also to the noise of the sample because of differences in the quality of extracted DNA. Therefore, we need to adjust the estimated values of these parameters on the basis of the sample noise.
In the following, we will denote the sample to estimate with s = 1 and the reference sample with s = 2.
for z = 2, 1, 0, 2. r_{1}(z) is the frequency of NoCall in regions with copy number event z of Sample 1, for z = 2, 1, 0, 2.
In Section "Results and Discussion", we will compare the estimations resulting from gBPCR with and without the adjustment of these parameters.
Results and Discussion
In this section, we apply gBPCR to both artificial and real data. First, we compare the boundary estimators (described in the previous section) on data simulated by using Model 1. Then, we evaluate the detection of IBD/UPD regions on the artificial dataset of [35], in comparison with three wellknown methods for LOH estimation. Using the same data, we also show the difference in the estimation when adjusting the parameters. Finally, we show the performance of gBPCR, when applied to real data.
With the current implementation, on a computer with dual CPU (AMD Opteron 250, 2.4 GHz) and 4 GB RAM, the algorithm needed almost two days to estimate the profile of an Affymetrix GeneChip Mapping 250K NspI sample (using k_{ max }= 50). Nevertheless, the computations can be performed by chromosome (and by arm for the longest chromosomes), reducing the computational cost. In any case, an optimized version of the code will be soon available.
Comparison of the breakpoint estimators on simulated data
In Section "Methods", we have described several possible boundary estimators: , and . The last one actually defines a class of estimators which depend on the values of the thresholds thr_{1} and thr_{2}. We tried several pairs of the following types of thresholds:

"005" := max(0.005, quantile of p at 0.95)

"01" := max(0.01, quantile of p at 0.95)

"01 90" := max(0.01, quantile of p at 0.90)

"mad" := median(p) + 3 * mad(p)
where mad is the median absolute deviation and p is the vector of posterior probabilities to be a breakpoint. All these thresholds derive from different definitions of which probability values are to be considered significant.
We assessed the quality of all the estimators of k_{0} and t^{0} considered, by applying them on two artificial datasets (called datasets A and B), each of 100 samples. We used as prior probabilities of heterozygosity the frequencies of heterozygosity (in the CEU population of the HapMap project [29]) given by the annotation file of the Affymetrix GeneChip Mapping 250K NspI microarray. Just for illustrative purpose and because of limited computational time, we considered only the SNPs of a single chromosome (chromosome 22), hence the number of data points in each sample is n = 2520. Since our complete model (Model 3) does not provide a realistic way to simulate IBD/UPD regions and the identification of gained regions depends mainly on copy number data, the samples were generated using Model 1.
Simulation description
Prior distribution of Z in the simulated dataset B.
segment  

prior  I  II  III  IV  V  VI  VII  VIII  IX  X  XI  XII  XIII  XIV  XV 
P(^{ Z }p ^{= 2)}  0  0.1  0  0.1  0.5  0.1  0  0  0.1  0.5  0  0.1  0.5  0.1  0 
P(^{ Z }p ^{= 0)}  0.1  0.6  0.1  0.6  0.4  0.6  0.1  0.1  0.6  0.4  0.1  0.6  0.4  0.6  0.1 
P(^{ Z }p ^{= 1)}  0.6  0.3  0.6  0.3  0.1  0.3  0.6  0.4  0.3  0.1  0.6  0.3  0.1  0.3  0.6 
P(^{ Z }p ^{= 2)}  0.3  0  0.3  0  0  0  0.3  0.5  0  0  0.3  0  0  0  0.3 
 1.
we generated a true profile of the homozygous status in normal cells X ^{ N }, by using the prior probabilities of heterozygosity, described previously;
 2.
we generated a true copy number event profile , by using the prior distribution of Z (notice that in some cases the final profile can have less than 15 segments, since, if consecutive segments have the same copy number value, then they are joined together);
 3.
given the true copy number event profile and the profile of the homozygous status in normal cells, we generated Y (the profile of the homozygous status in cancer cells detected by the genotype calling algorithm), by using the conditional probability distributions of Model 1.
Results of the comparisons
To evaluate the performance of the estimators, we computed several error measures regarding the estimation of the number of segments (01, absolute and squared errors), the boundary estimation (binary error, sensitivity and false discovery rate, FDR) and the profile estimation (sum of squared distance, SSQ, sum 01 error, sensitivity and FDR for all copy number events). The explanation of these error measures can be found in Section S.5 in Additional file 2.
Comparison among the breakpoint estimators with respect to error measures regarding copy number aberration detection.
dataset  method  sum 01 err  SSQ 


( , )  51.53  86.08  0.19  
A  ( , )  146.91  596.78  0.49 
( )  91.99  345.64  0.37  
( , )  421.79  1226.59  0.70  
( )  110.39  287.21  0.34  
( )  109.39  286.15  0.34  
B  ( )  141.65  370.78  0.38 
( )  154.56  424.2  0.41  
( )  109.39  286.15  0.34  
( )  111.75  283.77  0.34 
Based upon these results, we decided to not apply the estimators ( , ) on dataset B and to try other paired thresholds for , in order to reduce the FDR of the boundary estimation. By looking globally at the results of all error measures (see Table 3, Supplementary Tables S.2S.5 and Supplementary Figures S.2 and S.3 in Additional file 2), we can suggest the use of the following pairs of estimators: ( ), ( ) or ( ). Moreover, from the study of the behavior of ( ) and ( ), we can understand the role of the two thresholds in our algorithm for the determination of the maxima in a multimodal function (see Section S.4 in Additional file 2). The threshold thr_{1} is used to decide which points belong to the same peak: all the points, between two regions of points below thr_{1}, are considered in the same peak. Hence, with a low threshold, more points are considered belonging to the same peak and thus we can eliminate lot of false breakpoints (like in ( )). But, at the same time, if two true peaks are close, then it is possible that they are considered as only one peak, losing a true breakpoint (low sensitivity). Instead, the threshold thr_{2} is used to choose which estimated breakpoints are significant for the regression, i.e. if their posterior probabilities are to be considered different from zero. Therefore, using a lower value of thr_{2}, we select a higher number of breakpoints obtaining a higher percentage of both false ones (high FDR) and true ones (high sensitivity, as in ( )).
A detailed description of the results obtained in the comparison is in Section S.5 in Additional file 2.
Comparisons on simulated data with LOH regions
In order to evaluate the IBD/UPD detection of gBPCR, we applied it to simulated data of [35]. These data are based on three real samples of the HapMap dataset (see [1]), obtained with the Affymetrix GeneChip Mapping 250K NspI. For each sample and signal to noise ratio (SNR) value, the authors simulated two profiles: one with regions of copyneutral LOH and one with regions of loss. In both cases the number of regions was 50 and their width ranged from 20 SNPs to a whole chromosome. The values of SNR considered were 5, 2 and 1.25. Therefore, the total number of samples was 18, because for each normal sample we had two LOH profiles and three SNR values. The authors simulated the noise and the aberrations at probe level intensity saving the data in .CEL file format, thus we used BRLMM [26] to extract the genotyping data and CNAT 4.01 [36] for the raw copy number data.
Similar to [35], the estimation of gBPCR was compared with the ones given by three wellknown methods in the field: dChip [16], CNAT 4.01 [36] and PennCNV [24]. The evaluation has been done by computing the true positive rate (TPR) and the false positive rate (FPR), i.e. the proportion of SNPs inside the LOH regions that are correctly identified (as belonging to a LOH region) and the proportion of SNPs outside these segments that are wrongly identified (as belonging to them), respectively. We used ( ), ( ) or ( ) as paired estimators of the number of segments and the boundaries, and either p_{ upd }= 10^{3} or p_{ upd }= 10^{4} mas the prior probability of IBD/UPD.
Since CNAT does not consider the NoCall SNPs (called noninformative SNPs) for the estimation of the LOH profile, we compared the TPR and FPR computed using only either the informative or the noninformative SNPs (see Supplementary Figures S.4, S.5, S.6 and S.7 in Additional file 2).
Finally, we also evaluated the effect of the adjustment of the model parameters related to the NoCall detection (see Section "Methods"), using the same data. At low or medium noise, no significant differences in the goodness of the estimation could be observed (see, for example, Supplementary Figure S.10 in Additional file 2). Instead, in presence of high noise, the FPR regarding the IBD/UPD detection without the adjustment of the model parameters was close to one. In fact, in this situation a segment with normal copy number is more often classified as IBD/UPD, since the NoCall rate is higher and, without the correction, the IBD/UPD segments are allowed to contain a higher percentage of NoCall s with respect to the normal ones. Instead, with the adjustment, all types of regions are allowed to have a higher number of NoCall s in proportion to the noise, obtaining a less biased estimation.
In conclusion, we suggest to use ( ) or ( ) with p_{ upd }= 10^{4}, due to their good results obtained on the noninformative SNPs. A detailed description of all the results is in Section S.6 in Additional file 2.
Application to real data
In this subsection, we show the behavior of gBPCR on three real datasets. The first dataset consisted of eight paired cancer samples of patients affected by chronic lymphocytic leukemia (CLL), which then developed in diffuse large Bcell lymphoma (DLBCL), see [37, 38]. For two patients we had also a third sample, thus the total number of samples was 18. The second dataset consisted of 18 patients affected by CLL, see [39]. For both of these datasets, genomewide DNA profiles were obtained using the GeneChip Human Mapping 250K NspI (Affymetrix, Santa Clara, CA, USA). The genotype calls were calculated with BRLMM [26] using 46 Caucasian normal female samples of the HapMap Project as reference samples and the raw copy number data were retrieved using CNAT 4.01 [36]. In [37–39], the copy number of some genomic regions was also measured with fluorescent in situ hybridization (FISH). Therefore, on these regions we compared the copy number event estimated by gBPCR with the copy number measured by FISH. Moreover, since samples coming from the same patient should present the same copyneutral LOH regions (the germ line ones) for the majority of the genome, we used the two patients with three samples to evaluate the IBD/UPD detection of gBPCR.
The third dataset was a dilution series of the CRL2324 breast cancer cell line from [23]. The series comprised 12 samples, corresponding to the following percentages of tumor content: 0, 10, 14, 21, 23, 30, 34, 45, 47, 50, 79, 100. The genomewide DNA profiles were obtained using Illumina 370K BeadChips. The authors preprocessed the data with BeadStudio software (Illumina Inc.) and we used both the genotyping and logRratio data available at [40]. In [23], the authors chose eight genomic aberrations and compared the estimation given by their method (called BAFsegmentation) with the ones of the following algorithms: dChip [16], PennCNV [24], QuantiSNP [20] and SOMATICs [18]. Thus, we compared the estimations of these genomic regions given by gBPCR with the ones given by the previous methods. We also used these data to evaluate the performance of gBPCR in the detection of gains, for different values of the allelic copy numbers.
Based on the previous results on simulated data, for the analysis of these real data, we used: ( ), as paired estimators of the number of segments and the boundaries, and p_{ upd }= 10^{4} as prior probability of IBD/UPD.
Results regarding the identification of the copy number changes in CLL samples
In the samples considered for the comparison, we had a total of 169 regions measured by FISH (which provides also an estimate of the percentage of cells bearing the aberration): 38 regions were gains or losses in at least 60% of the cells (called detectable aberrations), 33 were gains or losses in less than 60% (called less detectable aberrations) and 98 were identified as normal segments. Regarding the detectable aberrations, only two copy number events were not identified by our algorithm. One loss was not found, because the estimated log_{2}ratio was very close to zero, and the other, because of a different percentage of Het SNPs from what was expected by our algorithm. We discovered 13 of the 33 less detectable copy number changes and we also detected two of the 98 normal segments as aberrations, but one of these copy number changes was equal to the one discovered in the same region of the paired sample, thus it was likely to be true.
Instead, by simply using the thresholds of the prior of Z on the profiles estimated by mBPCR for the classification of the copy number events (similarly to what is usually done), we detected one alteration less than what found by gBPCR and other 5 normal regions were seen as aberrations.
For the analysis of the results, we have to consider that the samples used for FISH came from peripheral blood, for the CLL samples, and from paraffn embedded tissues or lymph node, for the DLBCLs. Because of the consequently different cell content, in the former case, the results are better estimated. Moreover, the samples used for microarray and FISH might not be exactly the same, hence the percentage of cells which carry the aberrations can be different and a discordance between the two techniques is possible. Thus, gBPCR performed well in estimating the copy number changes on these samples.
Results regarding IBD/UPD region detection in CLL samples
For the evaluation of the IBD/UPD region detection, we considered the two patients with three samples. For the first patient (called Patient 1), we had: one matched normal DNA sample extracted from peripheral blood granulocytes (called Sample 1.1), one sample from neoplastic cells at CLL phase (called Sample 1.2) and the last one from neoplastic cells at DLBCL phase (called Sample 1.3). For the second patient (Patient 2), we had: one sample from neoplastic cells at CLL phase (called Sample 2.1), one at DLBCL phase (called Sample 2.2) and the last one from neoplastic cells at a further progression of the DLBCL (called Sample 2.3).
Applying gBPCR to the three samples of Patient 1, we found that the number of aberrations in each sample increased with the progression of the disease. The lower number of segments discovered in Sample 1.1 could also be due to a higher NoCall rate in comparison to the other samples. The same happened for Sample 2.3 of Patient 2.
We compared the IBD/UPD segments found in the three samples of each patient and we divided them into three classes (see Supplementary Table S.6 in Additional file 2):

equal regions: segments that are exactly the same in two or three samples;

overlapping regions: segments that are common in at least two samples but do not have the same boundaries;

single sample regions: the remaining segments.
Then, we defined the number of distinct regions as the sum of all these regions and the number of validated ones as the sum of all types of regions except the single sample regions. The proportions of equal and overlapping regions were similar in the two patients and the validated regions were 73% of the distinct regions detected in Patient 1 and 79% of the distinct regions in Patient 2. The single sample regions were about the 21% of the distinct regions in Patient 2, but the majority of them had length less than 50 SNPs. Instead, since the samples of Patient 1 belonged to different stages of the disease, in this patient we found a higher number of single sample regions and most of them were wider than 50 SNPs. In fact, the majority of these regions was detected in Sample 1.3, thus they were likely to be somatic.
Results regarding the identification of genomic aberrations in the dilution series
In [23], the authors observed that the BeadStudio normalization produced copy number profiles which were centered differently as the tumor content decreased and, as a consequence, many algorithms wrongly assigned the type of genetic aberration. Therefore, they evaluated the methods by considering only if they found any aberration in the eight regions considered, without looking at the type of aberration.
Due to this variation in centering the normal copy number, we estimated the histogram of the estimated log_{2}ratio values (which is used for the definition of the prior of Z), separately by using only samples with similar tumor content. Nevertheless, this shrewdness was not suffcient to well distinguish the peaks of the histogram in some cases.
For all samples, we computed the sensitivity in detecting the eight aberrations considered in [23]. For each of them, the calculations were done in two ways: by looking if gBPCR found any aberration (like in [23]) and by looking if it found the correct lesion (see Supplementary Figures S.11 and S.12 in Additional file 2). By comparing the results obtained by gBPCR using the first type of sensitivity with the ones given by the other methods in Figure 7 of [23], we can observe that gBPCR outperformed dChip and PennCNV and often also QuantiSNP. Sometimes it also performed better than BAFsegmentation and SOMATICs in the detection of the regions of gain. Moreover, occasionally gBPCR had a nonzero sensitivity in the normal sample, because it detected small IBD/UPD regions. By looking at the results obtained by gBPCR with the second type of sensitivity, we can notice that the correct aberration was usually detected in samples with at least 60% of tumor content and the sensitivity was still often higher than the one of dChip and PennCNV.
Finally, we computed the sensitivity of gBPCR for eight regions of gain, to evaluate its performance depending on the value of the copy number of the alleles. For all eight lesions, the total copy number was four. Instead, the minor allele copy number (maCN , i.e. the copy number of the allele less frequent in a normal population) changed from two to zero. The selection of these regions of gain was based on the estimated genomic profile of CRL2324, provided by The Cancer Genome Project at the Wellcome Trust Sanger Institute and available at [41]. For each aberration, the sensitivity was computed in two ways: by looking if the region was identified as a gain and by looking if it was detected as either a gain or an IBD/UPD segment (see Supplementary Figure S.13 in Additional file 2).
The differences between the two types of sensitivity were observed for some percentages of tumor content, in gains with maCN = 2 or maCN = 0. Regarding the lesions with maCN = 2, a small part of the gain was identified as IBD/UPD region in few cases with a small percentage of tumor content. This phenomenon was due to the presence of a high percentage of homozygous SNPs with copy number close to the normal copy number. For the same reason, the whole gain 6q22.31 (maCN = 0) was identified as an IBD/UPD region at 100% and 79% percentages of tumor content and the same happened also for a part of 6q15 (maCN = 0) at 79%. As we explained in Section "Method", the detection of the gains highly depends on the copy number value. Thus, if the copy number of a region of gain is close to the normal value, it is identified as either normal or IBD/UPD, depending on the homozygous status of the SNPs inside it. Therefore, the performance of gBPCR depends mainly on the quality of the copy number data and not on the value of the copy number of the alleles.
Conclusions
We have derived a new algorithm (called gBPCR) for the simultaneous estimation of copy number changes and IBD/UPD regions, by using both copy number and genotyping data. To the best of our knowledge, only one other algorithm exists which uses the same input data for the same purpose [17], but it does not appear appropriate for data coming from a DNA sample of a mixture of cell populations (like cancer DNA samples).
Our model takes into account the errors due to both the microarray procedure and the biological processes that lead to aberrations affecting the DNA copy number and the homozygous status. Because of the complexity of the algorithm and the high noise of the real data, we introduced new estimators to improve the detection of the breakpoints. On the basis of the results on simulated data, we selected the best performing one: ( ).
On the artificial dataset of [35] (and especially in samples with high noise), gBPCR outperformed three wellknown methods which estimate regions of LOH: dChip [16], CNAT [36] and PennCNV [24]. We also tested gBPCR on real data. On 36 CLL samples [37–39], we found a high agreement between the copy number changes estimated by gBPCR and the ones obtained by FISH (used as reference). Moreover, on two patients with three samples we could validate at least 73% of the identified IBD/UPD segments. On the samples of the CRL2324 dilution series of [23], we showed that in samples with at least 60% of tumor content, gBPCR was able to detect the genomic aberrations, while with less tumor content only some aberrations could be seen. Moreover, on these data gBPCR outperformed dChip [16] and PennCNV [24] and sometimes QuantiSNP [20]. Since other methods (SOMATICs [18] and BAFsegmentation [23]), which use the allelic copy number information, seemed to perform well, as future work we intend to add also this useful information in our model.
Availability and requirements
Project name: gBPCR.
Project home page: http://www.idsia.ch/~paola/gBPCR/.
Operating systems: the software should run in Linux, MacOS or Windows. Tests were performed on Windows and Linux systems.
Programming language: R.
Other requirements: none.
Licence: GNU GPL.
Any restrictions to use by nonacademics: none.
Declarations
Acknowledgements
This work was supported by Swiss National Science Foundation (grants 205321112430, 205320121886/1); Oncosuisse grants OCS193982006 and OCS  02296082008; Cantone Ticino ("Computational life science/Ticino in rete" program); Fondazione per la Ricerca e la Cura sui Linfomi (Lugano, Switzerland).
Authors’ Affiliations
References
 The International HapMap Consortium: The International HapMap Project. Nature 2003, 426: 789–796. 10.1038/nature02168View ArticleGoogle Scholar
 Kotzot D: Complex and segmental uniparental disomy (UPD): review and lessons from rare chromosomal complements. Journal of Medical Genetics 2001, 38: 497–507. 10.1136/jmg.38.8.497View ArticlePubMedPubMed CentralGoogle Scholar
 The international HapMap Consortium: A second generation human haplotype map of over 3.1 million SNPs. Nature 2007, 449: 851–862. 10.1038/nature06258View ArticlePubMed CentralGoogle Scholar
 Li LH, Ho SF, Chen CH, Wei WC, Wong CY, Li LY, Hung SI, Chung WH, Pan WH, Lee MTM, Tsai FJ, Chang CF, Wu JY, Chen YT: Long Contiguous Stretches of Homozygosity in the Human Genome. Human Mutation 2006, 27(11):1115–1121. 10.1002/humu.20399View ArticlePubMedGoogle Scholar
 Broman KW, Weber JL: Long homozygous chromosomal segments in reference families from the Centre d'Etude du Polymorphisme Humain. American journal of human genetics 1999, 65: 1493–1500. 10.1086/302661View ArticlePubMedPubMed CentralGoogle Scholar
 Bacolod MD, Schemmann GS, Wang S, Shattock R, Giardina SF, Zeng Z, Shia J, Stengel RF, Gerry N, Hoh J, Kirchhoff T, Gold B, Christman MF, Offt K, Gerald WL, Notterman DA, Ott J, Paty PB, Barany F: The Signatures of Autozygosity among Patients with Colorectal Cancer. Cancer Research 2008, 68(8):2610–2621. 10.1158/00085472.CAN075250View ArticlePubMedPubMed CentralGoogle Scholar
 Beà S, Salaverria I, Armengol L, Pinyol M, Fernández V, Hartmann EM, Jares P, Amador V, Hernández L, Navarro A, Ott G, Rosenwald A, Estivill X, Campo E: Uniparental disomies, homozygous deletions, amplifications and target genes in mantle cell lymphoma revealed by integrative highresolution whole genome profiling. Blood 2009, 113(13):3059–3069. 10.1182/blood200807170183View ArticlePubMedPubMed CentralGoogle Scholar
 Gondek LP, Tiu R, O'Keefe CL, Sekeres MA, Theil KS, Maciejewski JP: Chromosomal lesions and uniparental disomy detected by SNP arrays in MDS, MDS/MPD, and MDSderived AML. Blood 2008, 111(3):1534–1542. 10.1182/blood200705092304View ArticlePubMedPubMed CentralGoogle Scholar
 Huang J, Wei W, Zhang J, Liu G, Bignell GR, Stratton MR, Futreal PA, Wooster R, Jones KW, Shapero MH: Whole Genome DNA Copy Number Changes Identified by High Density Oligonucleotide Arrays. Human Genomics 2004, 1(4):287–299.View ArticlePubMedPubMed CentralGoogle Scholar
 Fridlyand J, Snijders AM, Pinkel D, Albertson DG, Jain AN: Hidden Markov Models approach to the analysis of array CGH data. Journal of Multivariate Analysis 2004, 90: 132–153. 10.1016/j.jmva.2004.02.008View ArticleGoogle Scholar
 Hupé P, Stransky N, Thiery JP, Radvanyi F, Barillot E: Analysis of array CGH data: from signal ratio to gain and loss of DNA regions. Bioinformatics 2004, 20(18):3413–3422. 10.1093/bioinformatics/bth418View ArticlePubMedGoogle Scholar
 Olshen AB, Venkatraman ES, Lucito R, Wigler M: Circular Binary Segmentation for the Analysis of Arraybased DNA Copy Number Data. Biostatistics 2004, 5(4):557–572. 10.1093/biostatistics/kxh008View ArticlePubMedGoogle Scholar
 Picard F, Robin S, Lavielle M, Vaisse C, Daudin JJ: A Statistical Approach for Array CGH Data Analysis. BMC Bioinformatics 2005., 6(27):Google Scholar
 Rancoita PMV, Hutter M, Bertoni F, Kwee I: Bayesian DNA copy number analysis. BMC Bioinformatics 2009., 10(10):Google Scholar
 Newton MA, Lee Y: Inferring the Location and Effect of Tumor Suppressor Genes by InstabilitySelection Modelling of AllelicLoss Data. Biometrics 2000, 56: 1088–1097. 10.1111/j.0006341X.2000.01088.xView ArticlePubMedGoogle Scholar
 Beroukhim R, Lin M, Park Y, Hao K, Zhao X, Garraway LA, Fox EA, Hochberg EP, Mellingho IK, Hofer MD, Descazeaud A, Rubin MA, Meyerson M, Wong WH, Sellers WR, Li C: Inferring LossofHeterozygosity from Unpaired Tumors Using HighDensity Oligonucleotide SNP Arrays. PLOS Computational Biology 2006, 2(5):323–332. 10.1371/journal.pcbi.0020041View ArticleGoogle Scholar
 Scharpf RB, Parmigiani G, Pevsner J, Ruczinski I: Hidden Markov Models for the assessment of chromosomal alterations using highthroughput SNP arrays. Annals of Applied Statistics 2008, 2(2):687–713. 10.1214/07AOAS155View ArticlePubMedPubMed CentralGoogle Scholar
 Assié G, LaFramboise T, Platzer P, Bertherat J, Stratakis CA, Eng C: SNP Arrays in Heterogeneous Tissue: Highly Accurate Collection of Both Germline and Somatic Genetic Information from Unpaired Single Tumor Samples. The American Journal of Human Genetics 2008, 82: 903–915. 10.1016/j.ajhg.2008.01.012View ArticlePubMedGoogle Scholar
 Attiyeh EF, Diskin SH, Attiyeh MA, Mossé YP, Hou C, Jackson EM, Kim C, Glessner J, Hakonarson H, Biegel JA, Maris JM: Genomic copy number determination in cancer cells from single nucleotide polymorphism microarrays based on quantitative genotyping corrected for aneuploidy. Genome Research 2009, 19: 276–283. 10.1101/gr.075671.107View ArticlePubMedPubMed CentralGoogle Scholar
 Colella S, Yau C, Taylor JM, Mirza G, Butler H, Clouston P, Bassett AS, Seller A, Holmes CC, Ragoussis J: QuantiSNP: an Objective Bayes HiddenMarkov Model to detect and accurately map copy number variation using SNP genotyping data. Nucleic Acids Research 2007, 35(6):2013–2025. 10.1093/nar/gkm076View ArticlePubMedPubMed CentralGoogle Scholar
 Huang J, Wei W, Chen J, Zhang J, Liu G, Di X, Mei R, Ishikawa S, Aburatani H, Jones KW, Shapero MH: CARAT: A novel method for allelic detection of DNA copy number changes using high density oligonucleotide rrays. BMC Bioinformatics 2006, 7: 83.View ArticlePubMedPubMed CentralGoogle Scholar
 LaFramboise T, Weir BA, Zhao X, Beroukhim R, Li C, Harrington D, Sellers WR, Meyerson M: AlleleSpecific Amplification in Cancer Revealed by SNP Array Analysis. PLOS Computational Biology 2005, 1(6):e65. 10.1371/journal.pcbi.0010065View ArticlePubMedPubMed CentralGoogle Scholar
 Staaf J, Lindgren D, VallonChristersson J, Isaksson A, G¨oransson H, Juliusson G, Rosenquist R, H¨oglund M, Borg A, Ringnér M: Segmentationbased detection of allelic imbalance and lossofheterozygosity in cancer cells using whole genome SNP arrays. Genome Biology 2008., 9(R136):Google Scholar
 Wang K, Li M, Hadley D, Liu R, Glessner J, Grant SFA, Hakonarson H, Bucan M: PennCNV: an integrated hidden Markov model designed for highresolution copy number variation detection in wholegenome SNP genotyping data. Genome Research 2007, 17(11):1665–1674. 10.1101/gr.6861907View ArticlePubMedPubMed CentralGoogle Scholar
 Lamy P, Andersen CL, Dyrskjot L, Torring N, Wiuf C: A Hidden Markov Model to estimate population mixture and allelic copynumbers in cancers using Affymetrix SNP arrays. BMC Bioinformatics 2007., 8(434):Google Scholar
 Affymetrix: BRLMM: an Improved Genotyping Calling Method for the GeneChip Human Mapping 500 K Array Set. 2006.Google Scholar
 Zhao X, Li C, Guillermo Paez J, Chin K, Jänne PA, Chen TH, Girard L, Minna J, Christiani D, Leo C, Gray JW, Sellers WR, Meyerson M: An Integrated View of Copy Number and Allelic Alterations in the Cancer Genome Using Single Nucleotide Polymorphism Arrays. Cancer Research 2004, 64: 3060–3071. 10.1158/00085472.CAN033308View ArticlePubMedGoogle Scholar
 Lai Y, Zhao H: A statistical method to detect chromosomal regions with DNA copy number alterations using SNParraybased CGH data. Computational Biology and Chemistry 2005, 29: 47–54. 10.1016/j.compbiolchem.2004.12.004View ArticlePubMedGoogle Scholar
 The international HapMap Consortium: A haplotype map of the human genome. Nature 2005, 437: 1299–1320. 10.1038/nature04226View ArticlePubMed CentralGoogle Scholar
 Veltman JA, Fridlyand J, Pejavar S, Olshen AB, Korkola JE, DeVries S, Carroll P, Kuo WL, Pinkel D, Albertson D, CordonCardo C, Jain AN, Waldman FM: Arraybased Comparative Genomic Hybridization for GenomeWide Screening of DNA Copy Number in Bladder Tumors. Cancer Research 2003, 63: 2872–2880.PubMedGoogle Scholar
 Nakao K, Mehta KR, Fridlyand J, Moore DH, Jain AN, Lafuente A, Wiencke JW, Terdiman JP, Waldman FM: Highresolution analysis of DNA copy number alterations in colorectal cancer by arraybased comparative genomic hybridization. Carginogenesis 2004, 25(8):1345–1357. 10.1093/carcin/bgh134View ArticleGoogle Scholar
 Capello D, Scandurra M, Poretti G, Rancoita PMV, Mian M, Gloghini A, Deambrogi C, Martini M, Rossi D, Greiner TC, Chan WC, Ponzoni M, Montes Moreno S, Piris MA, Canzonieri V, Spina M, Tirelli U, Inghirami G, Rinaldi A, Zucca E, Dalla Favera R, Cavalli F, Larocca LM, Kwee I, Carbone A, Gaidano G, Bertoni F: Genome wide DNAprofiling of HIVrelated Bcell lymphomas. British Journal of Haematology 2010, 148(2):245–255. 10.1111/j.13652141.2009.07943.xView ArticlePubMedGoogle Scholar
 Hodgson G, Hager JH, Volik S, Hariono S, Wernick M, Moore D, Nowak N, Albertson DG, Pinkel D, Collins C, Hanahan D, Gray JW: Genome scanning with array CGH delineates regional alterations in mouse islet carcinomas. Nature Genetics 2001, 29: 459–464. 10.1038/ng771View ArticlePubMedGoogle Scholar
 Forconi F, Poretti G, Kwee I, Sozzi E, Rossi D, Rancoita PMV, Capello D, Rinaldi A, Zucca E, Donatella Raspadori D, Spina V, Lauria F, Gaidano G, Bertoni F: High density genomewide DNA profiling reveals a remarkably stable profile in hairy cell leukaemia. British Journal of Haematology 2008, 141(5):622–630. 10.1111/j.13652141.2008.07106.xView ArticlePubMedGoogle Scholar
 Wu LY, Zhou X, Li F, Yang X, Chang CC: Conditional random pattern algorithm for LOH inference and segmentation. Bioinformatics 2009, 25: 61–67. 10.1093/bioinformatics/btn561View ArticlePubMedPubMed CentralGoogle Scholar
 Affymetrix: CNAT 4.0: Copy Number and Loss of Heterozygosity Estimation Algorithms for the GeneChip Human Mapping 10/50/100/250/500 K Array Set. 2007.Google Scholar
 Scandurra M, Rossi D, Deambrogi C, Rancoita PMV, Chigrinova E, Mian M, Cerri M, Rasi S, Sozzi F E Forconi, Ponzoni M, MontesMoreno S, Piris MA, Inghirami G, Zucca E, Gattei V, Rinaldi A, Kwee I, G G, Bertoni F: Genomic profiling of Richters syndrome: recurrent lesions and differences with de novo diffuse large Bcell lymphomas. Hematological Oncology 2010, 28(2):62–67.View ArticlePubMedGoogle Scholar
 Rossi D, Cerri M, Capello D, Deambrogi C, Rossi FM, Zucchetto A, De Paoli L, Cresta S, Rasi S, Spina V, Franceschetti S, Lunghi M, Vendramin C, Bomben R, Ramponi A, Monga G, Conconi A, Magnani C, Gattei V, Gaidano G: Biological and clinical risk factors of chronic lymphocytic leukaemia transformation to Richter syndrome. British Journal of Haematology 2008, 142(2):202–215. 10.1111/j.13652141.2008.07166.xView ArticlePubMedGoogle Scholar
 Forconi F, Rinaldi A, Kwee I, Sozzi E, Raspadori D, Rancoita PMV, Scandurra M, Rossi D, Deambrogi C, Capello D, Zucca E, Marconi D, Bomben R, Gattei V, Lauria F, Gaidano G, Bertoni F: Genomewide DNA analysis identifies recurrent imbalances predicting outcome in chronic lymphocytic leukaemia with 17 p deletion. British Journal of Haematology 2008, 143(4):532–536.PubMedGoogle Scholar
 BAFsegmentation[http://baseplugins.thep.lu.se/wiki/se.lu.onk.BAFsegmentation]
 The Cancer Genome Project[http://www.sanger.ac.uk/cgibin/genetics/CGP/cghviewer/CghHome.cgi]
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.