An experimental loop design for the detection of constitutional chromosomal aberrations by array CGH
BMC Bioinformatics volume 10, Article number: 380 (2009)
Comparative genomic hybridization microarrays for the detection of constitutional chromosomal aberrations is the application of microarray technology coming fastest into routine clinical application. Through genotype-phenotype association, it is also an important technique towards the discovery of disease causing genes and genomewide functional annotation in human. When using a two-channel microarray of genomic DNA probes for array CGH, the basic setup consists in hybridizing a patient against a normal reference sample. Two major disadvantages of this setup are (1) the use of half of the resources to measure a (little informative) reference sample and (2) the possibility that deviating signals are caused by benign copy number variation in the "normal" reference instead of a patient aberration. Instead, we apply an experimental loop design that compares three patients in three hybridizations.
We develop and compare two statistical methods (linear models of log ratios and mixed models of absolute measurements). In an analysis of 27 patients seen at our genetics center, we observed that the linear models of the log ratios are advantageous over the mixed models of the absolute intensities.
The loop design and the performance of the statistical analysis contribute to the quick adoption of array CGH as a routine diagnostic tool. They lower the detection limit of mosaicisms and improve the assignment of copy number variation for genetic association studies.
Array Comparative Genomic Hybridization (array CGH) [1, 2], also called molecular karyotyping , detects copy number aberrations and variations at high resolution on a genomewide scale [4, 5]. Genomewide array CGH has been applied to detect chromosomal imbalances in patients with congenital anomalies and mental retardation [6–12]. An illustrative example is presented in Figure 1. Array CGH is a highly effective technique that is entering routine clinical use much faster than other microarray technologies. Indeed, compared to, for example, expression microarrays, array CGH enjoys several technical advantages: (1) genomic DNA samples are less prone to degradation than mRNA samples, (2) genomic DNA samples show much less variation between biological replicates than mRNA samples, and (3) interpretation of chromosomal imbalances is much easier than that of expression fingerprints. These advantages explain why array CGH for the diagnosis of constitutional anomalies is progressing faster towards the clinic than expression microarrays for the prediction of clinical outcome (e.g., in cancer), for which a few applications are now entering clinical practice [13–15].
Array CGH mostly competes with and is complementary to conventional karyotyping and Fluorescent In Situ Hybridization (FISH). Compared to conventional karyotyping, it offers a resolution between 10 kb and 1 Mb, instead of about 5 Mb, and detects at least twice as many aberrations . Furthermore, it does not require the use of metaphase chromosomes, which makes it faster and less labor intensive. However, current array CGH techniques cannot detect balanced translocations, while this is straightforward with conventional karyotyping. Compared to FISH, array CGH provides genomewide coverage, instead of covering only a limited set of probes--so, it does not require prior knowledge of which aberration might be present (based on the phenotype of the patient).
The most frequent experimental setup for array CGH consists in comparing genomic DNA of a patient (test) with that of a normal individual (reference) using a two-channel microarray consisting of DNA segments spread across the whole genome. In the case of our clinical platform, the DNA segments consist in PCR-amplified BAC clones. However, the discussion applies equally to spotted long oligo platforms. So, we will refer to our probes as reporters. DNA from the test and reference samples is extracted, labeled with different fluorescent dyes (usually Cy3 and Cy5), hybridized to the microarray, and then scanned by two-channel laser scanner. Aneuploid chromosomal regions are detected as probes with a deviant log ratio of the intensities of the test against reference signal (approximately log2(1/2) for a deletion and log2(3/2) for a duplication). Usually the experiment is repeated in a dye-swap with the uorescent labeling of test and reference exchanged. The signals are then averaged over the dye-swap replicates to reduce the signal-to-noise ratio.
An alternative design  is a loop design in which three hybridizations are carried out with three test patients that are compared with each other: Patient 1 versus Patient 2, Patient 2 versus Patient 3, and Patient 3 versus Patient 1, as shown in Figure 2. This design measures the intensities of three test samples in a statistically balanced way and requires no normal reference sample. Hence, only three arrays are used to analyze three patients and to obtain two measurements for each of them. For the classical dye-swap design, half of the resources are consumed to measure the reference sample of a normal individual and, therefore, six arrays would be necessary to obtain as many measurements from the test samples.
Extensive genomic variation (called copy number variation (CNV)) is also present in normal individuals [16–18]. The extent of this variation is surprising (covering at least 10% of the genome) and likely to have major implications for human variation and disease. In the classical dye-swap design, a deviant log ratio for one reporter in the test sample could just as well be associated with a variation in the reference sample. The difficulty in disambiguating deviations between the test and reference sample prevents us also from replacing the reference sample with a second test sample in the dye-swap design. The loop design, on the contrary, unambiguously associates a deviation to the correct sample by looking for a unique pattern of log ratios. For example, a duplication in Patient 1 will be associated approximately to a positive log ratio in the Patient 1 vs. Patient 2 hybridization, a negative log ratio in the Patient 3 vs. Patient 1 hybridization, and a null log ratio in the Patient 2 vs. Patient 3 hybridization. No deletion or duplication in another patient will display the same pattern, so the association is unambiguous. Another way to elevate this issue in dye swap experiments would be to use a DNA sample, pooled from several individuals. However, for frequently occurring CNVs, the intensity ratios will be reduced and, therefore, pooling will rather complicate data interpretation instead of simplifying .
For the statistical analysis, we consider two approaches for array CGH: linear modeling of the log ratios and mixed modeling of the absolute signal intensities. We compare both methods on a test data set consisting of 27 patients (9 loop designs) and we implement the method with the best signal-to-noise ratio as a user-friendly web application. Both methods analyze the data in a clone by clone way. On high resolution arrays, the resulting estimates can however be used as input for segment wise analysis techniques. Experimental designs that make the best use of available resources are essential for the widespread adoption of array CGH as a routine clinical tool for the diagnosis of constitutional chromosomal aberrations. Reduction of false positives and negatives guarantees the best service to the patient and the best use of economic resources, which are key factors in a clinical environment. Furthermore, correct assignment of "benign" copy number variations to the right sample will be important in upcoming studies of association between copy number variations and disease. Finally, this design will also lower the detection limit for mosaicism (i.e., chromosomal aberrations present in only a percentage of the cells).
Array CGH for the detection of congenital chromosomal aberrations is also a key method for the genomewide discovery of gene function. Patients with chromosomal aberrations provide a natural form of forward genetics screen. Through genotype-phenotype correlations, phenotypes can be associated to chromosomal regions containing only a few to a few tens of candidate genes . Further prioritization of candidate genes using bioinformatics approaches and validation in small animal models (for example, zebrafish or fruit fly) allows rapid identification of disease causing genes .
Results and Discussion
New statistical models for the analysis of loop designs for array CGH
Two main philosophies vie for dominance for the statistical modeling of microarray signal intensities: (1) linear models of log ratios of intensities and (2) mixed models of absolute intensities. Two-channel microarrays were originally developed so that taking the ratio between the Cy3 and Cy5 intensity of a spot would eliminate multiple sources of variations (in particular, the amount of DNA material per spot). Linear models of log ratios keep in track with this philosophy and extend it by observing that log ratio measurements from different hybridizations containing equivalent samples are interdependent. They formulate these dependencies as a set of linear relations that are then inverted to obtain tighter estimates for a smaller set of independent statistical effects (which remain in essence similar to log ratios). By contrast, mixed models of absolute intensities aim at disentangling the signals from the Cy3 and Cy5 channels by expressing the intensities as the sum of an extensive set of fixed and random effects dependent on many factors that systematically affect the microarray measurements (dye effect, array effect, spot effect, etc.). Although the mixed models we consider are linear mixed models and thus also linear models, we will for simplicity refer in this paper to the two classes of methods as being "linear models" vs. "mixed models", respectively.
Mixed model of absolute signal intensities
We will apply a mixed model as proposed by . These models were originally applied to cDNA microarrays, but they can be tailored to the analysis of array CGH. Before applying a mixed model, we want to correct hybridization signals for possible spatial effects (which cannot be easily corrected with the mixed model). We therefore apply a 2D spatial loess correction to the hybridization intensities and obtain loess corrected log ratios, from which we can derive corrected intensities. The mixed models proposed by  consist of two successive models: the normalization and the reporter-specific model. The normalization model corrects for array, dye, and patient effects. The fitted model can be written as y cij = μ + P i + A j + PA ij + r cij , where y cij are the Cy3 and Cy5 intensities for Clone c, Patient i, and Array j, μ is the overall average, P i is the fixed patient effect with three levels (i = Patient 1, 2, or 3), A j estimates the random array effect (also with three levels j = Array 1, 2, or 3), and PA ij fits the interaction effect between the patient and array effect, and in this way it also corrects for the dye effect. For each reporter, we extract the residuals r cij from the normalization model and fit a reporter-specific model r cij = C c + CP ci + CA cj + ε cij , where r cij are the residuals obtained from the normalization model, C c is the overall average for Reporter or Clone c, CP ci is the fixed patient effect for Reporter c with three levels (i = Patient 1, 2, or, 3), CA cj estimates the random array effect for Reporter c (also with three levels j = Array 1, 2, or 3), and ε cij fits the random error effect. Our main interest is in the estimates of the CP ci effects, which reflect the difference between the patients for Reporter c. Specifically, we assess whether the contrasts Patient 2 vs. Patient 1 (= CPc 2- CPc 1) and Patient 1 vs. Patient 3 (= CPc 1- CPc 3) are equal to zero with a Wald's F-test. In the case where the contrast is significantly larger than zero for a chosen significance level α, we call this contrast positive. In the case where it is smaller than zero, it is called negative. Else we assign 0. Based on both hypothesis tests, the reporters are classified as duplicated or deleted according to the classification shown in Table 1. For example, if the contrast Patient 1 vs. Patient 3 (= CPc 1- CPc 3) is positive and the contrast Patient 2 vs. Patient 1 (= CPc 2- CPc 1) is negative for a reporter, then this reporter is likely to be duplicated for Patient 1. In some rare cases, we obtain as result a reporter that has, for example, a negative value for both contrasts CPc 1CPc 3and CPc 2CPc 1, which is none of the combinations in Table 1. In this case, we call the reporter strange.
Linear model of log ratios
An alternative statistical tool is a linear model of the log ratios, as proposed by . In contrast to the mixed model, this technique fits the 2D spatial loess corrected log2 ratios directly. In this particular case, we can choose the following contrasts Cc 1= log2(Pc 1/Pc 3) and Cc 2= log2(Pc 2/Pc 1), where P ci corresponds to the true underlying signal for Reporter c for Patient i. These contrasts correspond to the samples that were directly compared on the first two slides in Figure 2 and the observed log ratios on these slides should on average be equal to the contrast. The data of the third slide should then correspond on average to Cc 3= log2(Pc 3/Pc 2) = Cc 1-Cc 2. The linear model that fits the data can be written as
where E denotes the expectation of a random variable, X is the matrix of linear dependencies, C c is the vector of contrasts for Reporter c, and y ci denotes the log2 ratio for Reporter c measured on the ithslide. For each reporter, the least squares estimates of the three contrasts are obtained. To classify the contrasts as significantly duplicated, deleted, or not changed in copy-number, we apply the moderated t-statistic as implemented in the Bioconductor package limma, which implements linear models for microarray data analysis . The p-values from the moderated t-test were corrected to control the false discovery rate with Benjamini-Hochberg . Similarly to the mixed model, we can detect reporters that are duplicated or deleted for a patient, based on the p-values of the contrasts. For a chosen cut-off value α, we decide whether a reporter is not differentially expressed (0), upregulated (positive), or downregulated (negative) for a contrast. Based on the two contrasts, we can again classify a reporter as duplicated or deleted for a patient according to Table 1. Again, we can on rare occasions obtain strange reporters.
The mixed model and linear model provide two distinct ways to analyze the loop design experiments. To decide which method is preferable, we will first check which estimation method best separates the aberrant from the non-aberrant reporters. This will already give an indication to which method is preferred. Secondly, we will compare the false positive and false negative rates for a number of cut-off values α. Based on this information, we will decide which method to use and choose a cut-off value to call a reporter significantly duplicated or deleted. For the comparison of the analysis approaches, we consider a data set consisting out of nine loop designs from patients seen at our genetics center (Center for Human Genetics, U.Z.Leuven; see Table 2 and Methods section). Fifteen of the 27 patients involved in the nine loop designs showed one or more confirmed deletions or duplications of a chromosomal segment. Two experiments (Experiments 1 and 9) in our test data set include a sex mismatch. As for these experiments, the Y chromosome is absent for at least one of the patients, the measurements on the Y chromosome were excluded for both experiments from all computations. Because the X chromosome has regions with chromosome Y homology, the intensity ratios of chromosome X reporters are also more variable, and hence the X chromosome was also excluded from the computations for Experiments 1 and 9. For each of the nine loop designs, we classify the reporter as being deleted/duplicated if the reporter was classified as deleted or duplicated for one of the patients in the loop design. If the reporter was not deleted/duplicated for any of the patients in the loop design, the reporter was classified as non-aberrant. In total, this data set comprises 328 aberrant reporters: 116 deleted and 212 duplicated reporters. Over all nine experiments, we have a set of 30,668 measurements for non-aberrant reporters for any of the three patients in the loop.
Assessing which method is best capable of distinguishing between the intensities of aberrant and non-aberrant reporters can be done by computing signal-to-noise-ratios (SN) (for both deletions and duplications separately) as
where alteration type is deletion or duplication. As we have collected a data set with 212 duplicated reporters, 116 deleted, and 30,668 non-aberrant reporters, we can compute the SN values based on the absolute values of the contrasts Patient 1 vs. Patient 3 and Patient 2 vs. Patient 1, for both the linear model and the mixed model. The results are shown in Table 3. As for deleted clones the difference between their log ratios (± ) and log ratios of non-aberrant clones is larger than for the log ratios derived from duplicated clones (± ), the signal-to-noise-ratio is of course larger for the deleted reporters than for the duplicated reporters. The linear model leads to a significant reduction in the noise, especially for the non-aberrant reporters, and this results in a larger signal-to-noise-ratio. Therefore, these statistics are favorable to the linear model.
True positive and false positive rate
First, we compute for a number of significance levels α, the percentage of the duplicated and deleted reporters that are correctly classified as duplicated and deleted, respectively, in our test data set, according to both methods. The true positive (TP) rates are shown for the mixed model and the linear model in Figure 3 in function of the significance level α. For the linear model, the TP rate reaches a maximum of 0.954 for significance level α = 0.009. For higher α, surprisingly, the TP rate drops a little. This effect results from the fact that some reporters become strange reporters for larger significance levels α. For the classification with the mixed model, the TP rates grow slowly as the significance level increases. Within this range of significance levels α, it never reaches the maximum TP rate value obtained with the linear model. Perhaps it comes closer to the result obtained with the linear model if we allow for even larger significance levels α, but this will increase the FP rate and the number of strange reporters to an unacceptable level. Outside the duplicated and deleted regions, other reporters were also classified as duplicated or deleted. These positives can be false positives, due to technical artifacts, or they can indicate true biological variations. At this point, we will not make the distinction between both kinds of aberrant reporters, as it does not affect the method comparison, and we will refer to this set of positives as non-confirmed positives. At a later stage, this set of non-confirmed positives will be examined in depth, for one method and one significance level α. The number of these non-confirmed positives is shown in Figure 4. The figure shows that there is no clear difference in the non-confirmed positive rates between both methods. For low significance levels α the linear model has a slightly smaller number of non-confirmed positives. The combined results on the signal-to-noise ratio, the TP, and non-confirmed positive rate lead to the conclusion that the linear model is the preferable method.
Optimization of the linear model
In the previous section, we focused on how well the different methods fit the measurements by assessing their capabilities to divide the non-aberrant reporters from the deviating reporters and by comparing the FP and TP rates. This indicated that the linear model was best suited to distinguish these groups of reporters, although also the linear model has a fairly high FP rate. However, we did not yet benefit from all available information.
Completely and partially deleted reporters
To further optimize the reporter classification, we will use the fact that if, for example, a reporter of Patient 1 is deleted or duplicated, its contrast Patient 2 vs. Patient 1 should theoretically be equal to log2 () = 1 or log2 () = -0.58, respectively. However, nonlinear saturation effects in the signals cause a deviation from these values. Instead of taking the theoretically expected values (i.e., ± 1 and ± 0.58), we will estimate the expected values based on the linear model estimates of the contrasts for the group of confirmed deletions and duplications, after exclusion of the deletions and duplications on the X and Y chromosome. This results in an average for the absolute log ratios of 0.86 for the deleted reporters and 0.54 for the duplicated reporters. Therefore, if we detect with the linear model a reporter that is likely to be duplicated or deleted, we extract its contrasts Cc 1and Cc 2. If their absolute value is not larger than 0.54 or 0.86, respectively, we use an adapted version of the moderated t-test, as implemented in the limma package. We use the same standard deviation of the contrasts, as computed within the previous limma procedure, and test one-sidedly the hypothesis H0 : |C| = 0.54 versus H a : |C| < 0.54 for the duplicated reporters and H0 : |C| = 0.86 versus H a : |C| < 0.86 for the deleted reporters. If we cannot reject the hypothesis at a significance level αpartial for both tests, we call the reporter completely deleted or duplicated, else we call the reporter partially deviating (i.e., only a part of the clone has been deleted/duplicated instead of the complete clone). As a significance level, we choose αpartial = 0.01. For this significance level, the non-confirmed positives restricted to the reporters that are completely deleted are plotted as a blue line in Figure 4. This non-confirmed positives rate ranges between 0.001 and 0.002.
The non-confirmed positives
The non-confirmed positives rate obtained in the previous section is not a direct indication of the false positives rate, as they can comprise not only false positives, but also both true positives or polymorphic reporters. In our data set of nine loops, we saw that after extraction of clones that were proven to be copy number variable regions, 15 single clones with complete deletions or duplications are picked up by our method, at a significance level of α = 0.001. Seven of these clones have been previously described as polymorphic . One clone (RP1-93N13) overlaps with a CNV present in the normal population (Genomic Database of Variants; http://projects.tcag.ca/variation/). To determine whether the remaining reporters are false positives or true CNVs, qPCR or FISH  was performed. For six clones, the deletion or duplication was confirmed. For one clone we could not determine the copy number variation due to lack of DNA. The remaining four clones were all measured in Experiment 1 and caused by low spot quality, as they had a large discrepancy between the mean and the median pixel intensities and a large standard deviation of the pixel intensity. Hence, our method reliably detects aberrant clones. An overview of these results is shown in Table 4.
At α = 0.001, also 13 clones were classified as being strange. Ten of these clones overlapped with CNVs that are present in the normal population, according to the Genomic Database of Variants. (Note also that not all normal CNVs have already been identified.) Therefore, we suspect that a strange clone can often be explained as a polymorphic clone that is shared by two patients. Taken together these results indicate that at α = 0.001 our procedure has both low false negatives (TP rate around 95%) and essentially no false positives for completely deleted or duplicated clones. Most clones detected as positive or strange outside the regions known in the benchmark have been confirmed by qPCR or FISH; several of them being known normal CNVs.
The method is implemented as a web application and is available at http://www.esat.kuleuven.be/loop. A demo account with test data is available. Currently, the application and the statistical analysis have been tested and refined on an in-house series of over 400 patients. Details on the implementation and the use of the web application can be found in Additional file 1.
A first point of discussion is to which platforms the proposed statistical models are applicable. Some array CGH platforms are available using single-channel microarrays. In this situation, the loop design is simply not applicable. The loop design is directly dependent on using a two-channel microarray.
A second point of discussion is a main assumption that underlies our analysis, that two patients never share the same aberration. If two patients have an aberration for the same reporter, the statistical model cannot correctly interpret the result|resulting mostly into a strange reporter. For the detection of congenital anomalies, this effect is prevented by putting into a loop design patients who have clearly distinct phenotypic patterns. The rare cases where a common aberration is still present can be rescued at the validation stage. Incorrect assignment to the third patient would not be validated by FISH or qPCR and, in this case, the possibility of a common aberration for the two other patients should be kept in mind. Reporters flagged as strange can be validated in all three patients to clear out the situation.
At our genetics center, the pick-up rate (i.e., the percentage of patients that show an chromosomal aberration) is around 20%. Hence, a loop design consisting out of three patients with a CNVs occurs only in one out of 125 cases. The likelihood of having patients that share an abnormality is even smaller. However, this assumption is unrealistic when studying tumors or when doing preimplantation diagnosis. In the case of tumors, many chromosomal regions can be affected and overlap between patients is essentially unavoidable. In the case of preimplantation diagnosis, we have developed a procedure where one of the eight blastomere cells of an 8-cell embryo from in vitro fertilization is assessed by single-cell array CGH . In this procedure it is essentially impossible to guarantee that the aberrations from different blastomeres do not overlap (for example, some aberrant embryos have chaotic genomes where many chromosomes are affected).
In some rare cases, deviating reporters cannot be assigned as duplicated or deleted, but are instead labeled as strange. We observed that the majority of these strange reporters corresponded to polymorphic reporters. This can correspond to the situation where two patients share a complex polymorphic reporter  and will be investigated more deeply in a subsequent study.
Our analysis was characterized by the fact that for our optimized threshold we have a low number of false negatives and close to no false positives. Given that a secondary validation is available through FISH, which can catch false positives but not false negatives, we would want to increase our threshold further to decrease the number of false negatives at the cost of a few false positives that can then be caught by the FISH validation. However, increasing the threshold further results rather in more reporters being labeled as strange, which we prefer to avoid.
At a significance level of 0.001, for 0.251% of the clones the intensity ratios were suggestive of partial deletions of duplications. Most likely, the intermediate intensity ratios reflect the presence of subclonal CNVs (i.e., only a part of the clone has been deleted/duplicated instead of the complete clone). Recently, several articles show high levels of small CNVs in the human genome [16–18, 26]. Moreover, we performed qPCRs on some of these clones using several PCR primers and could confirm such partial deletions (data not shown). However, since no systematic analysis of all these clones has been performed, we do not know the number of false positives.
In our specific setting and in our hands, linear models clearly outperformed mixed models. We do not draw any conclusion with regard to the suitability of one class of method versus the other in general. In our setting, we hypothesize that the mixed model was less robust to deviation from the underlying normality assumptions or that the more compact (fewer parameters) estimation procedure of the linear model increased its robustness. The improvement of the mixed model in this setting appears to be an interesting research direction.
Microarray CGH for the diagnosis of congenital chromosomal aberrations is progressing rapidly from the research lab to the clinic. What value do such improved statistical procedures add to the diagnosis? First of all, increasing the signal-to-noise ratio between aberrant and normal clones through our statistical procedure is likely to improve the detection of low-grade mosaicism (i.e., chromosomal aberrations present in only a percentage of all cells), which is currently difficult to detect by array CGH. In a previous study , we presented a power analysis model for the detection of low-grade mosaicism where the signal-to-noise ratio was a critical factor in determining the limit of detection for mosaicism.
Secondly, CNVs are hard to assign with certainty to the patient instead of the "normal" reference. While these CNVs are currently largely handled as noise, some of them are likely to act as genetic modifiers that are risk factors for disease or modulate its penetrance and phenotypic spectrum . Association studies are currently being started to evaluate the importance of such CNVs in both congenital and acquired disorders. Unambiguous assignment of the "normal" CNVs to the correct sample is therefore paramount. And as a last point, there is the issue of cost per experiment. To keep cost down, some laboratories do not perform a dye swap when using the basic array CGH setup (patient vs. normal reference). In this case, a single array is used per patient. However, only one measurement per patient will be available in this case instead of two with our design, which means a higher level of false positives and negatives. False positives will be caught by the secondary FISH or qPCR validation, but false negatives will not. This results in lower quality diagnostics for the patient.
In conclusion, our results indicate that the experimental loop design, together with a statistical analysis by a linear model, provides an efficient procedure for the detection of chromosomal aberrations in congenital anomalies by array CGH. It is significantly superior to the classical setup by doubling the use of resources and unambiguously assigning variation to the correct patient. These improvements have a direct impact on the diagnosis offered to the patient using the microarray technology that is closest to routine clinical use.
A 1 Mb resolution BAC array was performed as described in . In short, 3500 BAC en PAC clones from the Welcome Trust Sanger Institute were amplified by two rounds of DOP-PCR . The purified aminolinked PCR products were spotted in duplicate on 3-D CodeLink Bioarray System slides (Amersham Biosciences, Piscataway, NJ). 150 ng of patient DNA was labeled by random prime labeling system (BioPrime Array CGH Genomic Labeling System, Invitrogen) using Cy3- and Cy5-dCTPs (Amersham Biosciences). Probe concentration and labeling efficiencies were measured using the nanodrop ND-1000 spectrophotometer (Nanodrop Technologies, Rockland, DE). The probe was placed on the slide under a glass cover slip (24 × 24 mm). The slides were incubated for 42-72 hours at 37°C under humified conditions.
Post hybridization washing was performed by soaking in 1 × PBS for 10 minutes and followed by a stringency wash on 42°C for 30 minutes in 50% formamide 2× SSC solution and 10 minutes in 1× PBS at room temperature. To finalise, the slides were dried by centrifugation for at least 1 minute.
Image and data analysis
After the washing, a two-channel scan was performed with an axon laser scanner GenePix 4000B (Molecular Devices, Union City, CA) at 532 nm and 635 nm using GenePixPro 6.0 program. The results are converted into GPR format and can be directly uploaded to the web-based application. The median spot intensities were corrected with the local median background, and only those spots with a signal above background (i.e., foreground intensity larger than local background intensity plus twice the local standard deviation of the background) were retained for the analysis. In this way, only few spots are lost, as almost all spots are above background (on average 96.6%). The ratios of the Cy5 to the Cy3 intensities were computed for each reporter and base 2 log transformed. The log ratios are normalized using a 2D spatial loess normalization, in which one applies a loess regression to fit the log2-ratios (M-values) on the coordinates on the slide as predictor variables.
Benchmark data set
For the comparison of the analysis approaches, we consider a data set consisting out of nine loop designs or 27 patients with mental retardation (MR) and multiple congenital anomalies (MCA). The patients were seen at our genetics center (Center for Human Genetics, U.Z.Leuven). Conventional karyotyping showed chromosomal imbalances in 11 patients. Analysis of the patients was carried out with a 1 Mb BAC array. In first instance, the data analysis strategy that we had previously developed was applied to this data. In this procedure, a region is called aberrant, if one clone passes the threshold of 4 × SD and if two or more flanking clones were passing the threshold of log2 () - 2 × SD as described in . If a deletion or a duplication larger then 3 Mb was detected, FISH was performed to confirm the results of the array. In case of a duplication smaller than 3 Mb, we performed quantitative PCR (qPCR) . 16 out of 27 patients show one or multiple clone anomalies, whereas 10 patients are apparently normal, at least according to the results of the array. One of the patients was a carrier of an inversion and another patient had a balanced translocation. Both aberrations cannot be detected by array CGH, but with conventional karyotyping; as such these patients did not contribute data to the benchmark. A short summary of the data set is shown in Table 2. In total, this data set comprises 635 aberrant clones: 274 deleted and 361 duplicated clones.
The study was approved by the institutional review board and appropriate informed consent was obtained from human subjects. The data set has been uploaded to the Gene Expression Omnibus (Accession number GSE6538) and is publicly available.
Written informed consent was obtained from the patient for publication of accompanying images. A copy of the written consent is available for review by the Editor-in-Chief of this journal.
Pinkel D, Segraves R, Sudar D, Clark S, Poole I, Kowbel D, Collins C, Kuo WL, Chen C, Zhai Y, Dairkee SH, Ljung BM, Gray JW, Albertson DG: High resolution analysis of DNA copy number variation using comparative genomic hybridization to microarrays. Nat Genet 1998, 20(2):207–211. 10.1038/2524
Solinas-Toldo S, Lampel S, Stilgenbauer S, Nickolenko J, Benner A, Dohner H, Cremer T, Lichter P: Matrix-based comparative genomic hybridization: biochips to screen for genomic imbalances. Genes Chromosomes Cancer 1997, 20(4):399–407. 10.1002/(SICI)1098-2264(199712)20:4<399::AID-GCC12>3.0.CO;2-I
Vermeesch JR, Rauch A: Reply to Hochstenbach et al. 'Molecular karyotyping'. Eur J Hum Genet 2006, 14(10):1063–1064. Comment Comment 10.1038/sj.ejhg.5201663
Snijders AM, Nowak N, Segraves R, Blackwood S, Brown N, Conroy J, Hamilton G, Hindle AK, Huey B, Kimura K, Law S, Myambo K, Palmer J, Ylstra B, Yue JP, Gray JW, Jain AN, Pinkel D, Albertson DG: Assembly of microarrays for genome-wide measurement of DNA copy number. Nat Genet 2001, 29(3):263–264. 10.1038/ng754
Pollack JR, Perou CM, Alizadeh AA, Eisen MB, Pergamenschikov A, Williams CF, Jeffrey SS, Botstein D, Brown PO: Genome-wide analysis of DNA copy-number changes using cDNA microarrays. Nat Genet 1999, 23: 41–46. 10.1038/14385
Shaw-Smith C, Redon R, Rickman L, Rio M, Willatt L, Fiegler H, Firth H, Sanlaville D, Winter R, Colleaux L, Bobrow M, Carter NP: Microarray based comparative genomic hybridisation (array-CGH) detects submicroscopic chromosomal deletions and duplications in patients with learning disability/mental retardation and dysmorphic features. J Med Genet 2004, 41(4):241–248. [Evaluation Studies] [Evaluation Studies] 10.1136/jmg.2003.017731
Rosenberg C, Knijnenburg J, Bakker E, Vianna-Morgante AM, Sloos W, Otto PA, Kriek M, Hansson K, Krepischi-Santos ACV, Fiegler H, Carter NP, Bijlsma EK, van Haeringen A, Szuhai K, Tanke HJ: Array-CGH detection of micro rearrangements in mentally retarded individuals: clinical significance of imbalances present both in affected children and normal parents. J Med Genet 2006, 43(2):180–186. [Letter] [Letter] 10.1136/jmg.2005.032268
Schoumans J, Ruivenkamp C, Holmberg E, Kyllerman M, Anderlid BM, Nordenskjold M: Detection of chromosomal imbalances in children with idiopathic mental retardation by array based comparative genomic hybridisation (array-CGH). J Med Genet 2005, 42(9):699–705. [Case Reports] [Case Reports] 10.1136/jmg.2004.029637
Vissers LELM, de Vries BBA, Osoegawa K, Janssen IM, Feuth T, Choy CO, Straatman H, Vliet W, Huys EHLPG, van Rijk A, Smeets D, van Ravenswaaij-Arts CMA, Knoers NV, Burgt I, de Jong PJ, Brunner HG, van Kessel AG, Schoenmakers EFPM, Veltman JA: Array-based comparative genomic hybridization for the genomewide detection of submicroscopic chromosomal abnormalities. Am J Hum Genet 2003, 73(6):1261–1270. 10.1086/379977
Menten B, Maas N, Thienpont B, Buysse K, Vandesompele J, Melotte C, de Ravel T, Van Vooren S, Balikova I, Backx L, Janssens S, De Paepe A, De Moor B, Moreau Y, Marynen P, Fryns JP, Mortier G, Devriendt K, Speleman F, Vermeesch JR: Emerging patterns of cryptic chromosomal imbalance in patients with idiopathic mental retardation and multiple congenital anomalies: a new series of 140 patients and review of published reports. J Med Genet 2006, 43(8):625–633. 10.1136/jmg.2005.039453
Friedman JM, Baross A, Delaney AD, Ally A, Arbour L, Asano J, Bailey DK, Barber S, Birch P, Brown-John M, Cao M, Chan S, Charest DL, Farnoud N, Fernandes N, Flibotte S, Go A, Gibson WT, Holt RA, Jones SJM, Kennedy GC, Krzywinski M, Langlois S, Li HI, McGillivray BC, Nayar T, Pugh TJ, Rajcan-Separovic E, Schein JE, Schnerch A, Siddiqui A, Van Allen MI, Wilson G, Yong SL, Zahir F, Eydoux P, Marra MA: Oligonucleotide microarray analysis of genomic imbalance in children with mental retardation. Am J Hum Genet 2006, 79(3):500–513. 10.1086/507471
de Vries BBA, Pfundt R, Leisink M, Koolen DA, Vissers LELM, Janssen IM, Reijmersdal Sv, Nillesen WM, Huys EHLPG, Leeuw Nd, Smeets D, Sistermans EA, Feuth T, van Ravenswaaij-Arts CMA, van Kessel AG, Schoenmakers EFPM, Brunner HG, Veltman JA: Diagnostic genome profiling in mental retardation. Am J Hum Genet 2005, 77(4):606–616. 10.1086/491719
Bogaerts J, Cardoso F, Buyse M, Braga S, Loi S, Harrison JA, Bines J, Mook S, Decker N, Ravdin P, Therasse P, Rutgers E, van 't Veer LJ, Piccart M: Gene signature evaluation as a prognostic tool: challenges in the design of the MINDACT trial. Nat Clin Pract Oncol 2006, 3(10):540–551. 10.1038/ncponc0591
Buyse M, Loi S, van't Veer L, Viale G, Delorenzi M, Glas AM, d'Assignies MS, Bergh J, Lidereau R, Ellis P, Harris A, Bogaerts J, Therasse P, Floore A, Amakrane M, Piette F, Rutgers E, Sotiriou C, Cardoso F, Piccart MJ: Validation and clinical utility of a 70-gene prognostic signature for women with node-negative breast cancer. J Natl Cancer Inst 2006, 98(17):1183–1192.
Shi L, Reid LH, Jones WD, Shippy R, Warrington JA, Baker SC, Collins PJ, de Longueville F, Kawasaki ES, Lee KY, Luo Y, Sun YA, Willey JM, Setterquist RA, Fischer GM, Tong W, Dragan YP, Dix DJ, Frueh FW, Goodsaid FM, Herman D, Jensen RV, Johnson CD, Lobenhofer EK, Puri RK, Schrf U, Thierry-Mieg J, Wang C, Wilson M, Wolber PK, Zhang L, Slikker WJ, Shi L, Reid LH: The MicroArray Quality Control (MAQC) project shows inter- and intraplatform reproducibility of gene expression measurements. Nat Biotechnol 2006, 24(9):1151–1161. [Comparative Study] [Comparative Study] 10.1038/nbt1239
Redon R, Ishikawa S, Fitch KR, Feuk L, Perry GH, Andrews TD, Fiegler H, Shapero MH, Carson AR, Chen W, Cho EK, Dallaire S, Freeman JL, Gonzalez JR, Gratacos M, Huang J, Kalaitzopoulos D, Komura D, MacDonald JR, Marshall CR, Mei R, Montgomery L, Nishimura K, Okamura K, Shen F, Somerville MJ, Tchinda J, Valsesia A, Woodwark C, Yang F, Zhang J, Zerjal T, Zhang J, Armengol L, Conrad DF, Estivill X, Tyler-Smith C, Carter NP, Aburatani H, Lee C, Jones KW, Scherer SW, Hurles ME: Global variation in copy number in the human genome. Nature 2006, 444(7118):444–454. 10.1038/nature05329
Sebat J, Lakshmi B, Troge J, Alexander J, Young J, Lundin P, Maner S, Massa H, Walker M, Chi M, Navin N, Lucito R, Healy J, Hicks J, Ye K, Reiner A, Gilliam TC, Trask B, Patterson N, Zetterberg A, Wigler M: Large-scale copy number polymorphism in the human genome. Science 2004, 305(5683):525–528. 10.1126/science.1098918
Iafrate AJ, Feuk L, Rivera MN, Listewnik ML, Donahoe PK, Qi Y, Scherer SW, Lee C: Detection of large-scale variation in the human genome. Nat Genet 2004, 36(9):949–951. 10.1038/ng1416
Van Buggenhout G, Melotte C, Dutta B, Froyen G, Van Hummelen P, Marynen P, Matthijs G, de Ravel T, Devriendt K, Fryns JP, Vermeesch JR: Mild Wolf-Hirschhorn syndrome: micro-array CGH analysis of atypical 4p16.3 deletions enables refinement of the genotype-phenotype map. J Med Genet 2004, 41(9):691–698. [Case Reports] [Case Reports] 10.1136/jmg.2003.016865
Aerts S, Lambrechts D, Maity S, Van Loo P, Coessens B, De Smet F, Tranchevent L, De Moor B, Marynen P, Hassan B, Carmeliet P, Moreau Y: Gene prioritization through genomic data fusion. Nature Biotechnology 2006, 24(5):537–544. 10.1038/nbt1203
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
Smyth G: Linear models and empirical bayes methods for assessing differential expression in microarray experiments. Statistical Applications in Genetics and Molecular Biology 2004., 3:
Benjamini Y, Hochberg Y: Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc Ser B 1995, 57: 289–300.
Vermeesch JR, Melotte C, Froyen G, Van Vooren S, Dutta B, Maas N, Vermeulen S, Menten B, Speleman F, De Moor B, Van Hummelen P, Marynen P, Fryns JP, Devriendt K: Molecular karyotyping: array CGH quality criteria for constitutional genetic diagnosis. J Histochem Cytochem 2005, 53(3):413–422. 10.1369/jhc.4A6436.2005
Le Caignec C, Spits C, Sermon K, De Rycke M, Thienpont B, Debrock S, Staessen C, Moreau Y, Fryns JP, Van Steirteghem A, Liebaers I, Vermeesch JR: Single-cell chromosomal imbalances detection by array CGH. Nucleic Acids Res 2006, 34(9):e68. 10.1093/nar/gkl336
Tuzun E, Sharp AJ, Bailey JA, Kaul R, Morrison VA, Pertz LM, Haugen E, Hayden H, Albertson D, Pinkel D, Olson MV, Eichler EE: Fine-scale structural variation of the human genome. Nat Genet 2005, 37(7):727–732. 10.1038/ng1562
Fiegler H, Carr P, Douglas EJ, Burford DC, Hunt S, Scott CE, Smith J, Vetrie D, Gorman P, Tomlinson IPM, Carter NP: DNA microarrays for comparative genomic hybridization based on DOP-PCR amplification of BAC and PAC clones. Genes Chromosomes Cancer 2003, 36(4):361–374. 10.1002/gcc.10155
The authors wish to thank the MicroArray Facility, Flanders Institute for Biotechnology (VIB) for their help in the spotting of the arrays and the Mapping Core and Map Finishing groups of the Wellcome Trust Sanger Institute for the initial clone supply and verification. This work was made possible by grants GOA/2006/12 and Centre of Excellence SymBioSys (Research Council K.U.Leuven EF/05/007), Catholic University of Leuven.
YM, JV, and BDM coordinated the study and formulated the methodology and research schema. JA performed the statistical analysis. JA and SVV contributed to the implementation. FH carried out the microarray experiments and did the validation experiments. JA, FH, YM, and JV interpreted the results. JA and YM wrote the manuscripts. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 1: Documentation on the web application. This additional file provides a description of the web application, implementing the algorithm presented in this paper. (PDF 32 KB)
About this article
Cite this article
Allemeersch, J., Van Vooren, S., Hannes, F. et al. An experimental loop design for the detection of constitutional chromosomal aberrations by array CGH. BMC Bioinformatics 10, 380 (2009). https://doi.org/10.1186/1471-2105-10-380
- Copy Number Variation
- Chromosomal Aberration
- True Positive Rate
- Loop Design
- Conventional Karyotyping