<title>Background</title> <p>With the recent technological advances, different genome-wide studies have uncovered an unprecedented number of structural variants in the human genome [##REF##15286789##1##, ####REF##15273396##2##, ##REF##15895083##3##, ##REF##15918152##4##, ##REF##16826518##5##, ##REF##17122850##6##, ##REF##17160897##7####17160897##7##], mainly in the form of copy number variations (CNVs). As much as 12% of the human genome has been reported to be variable among different individuals [##REF##17122850##6##]. The important number of genes and other regulatory elements encompassed by those variable regions, makes it very likely for them to have functional and, ultimately, phenotypical consequences [##REF##16418744##8##,##REF##17289997##9##]. In fact, several publications have already correlated the number of copies of different genes with different degrees of disease predisposition [##REF##15637236##10##, ####REF##16369530##11##, ##REF##17072318##12####17072318##12##]. Therefore, the identification of DNA copy number is important in understanding genesis and progression of human diseases.</p> <p>Several techniques and platforms have been developed for genome-wide analysisi of DNA copy number, such as array-based comparative genomic hybridization (aCGH). The goal of the analysis of DNA copy number data using this approach is to partition the whole genome into segments where copy numbers changes between contiguous segments may be present. The ability of aCGH to discern between different number of copies is very limited, thus different kinds of techniques have been developed for targeted, and more precise, analyses of genomic regions. Multiplex Ligation-dependent Probe Amplification (MLPA) [##REF##12060695##13##] is one of the most used technologies among other existing ones such us Quantitative Multiplex PCR of Short Fluorescent (QMPSF) [##REF##10850409##14##,##REF##12203994##15##] or Multiplex Amplifiable Probe Hybridization (MAPH) [##REF##10606661##16##].</p> <p>MLPA is a recently developed semi-quantitative method that aims to detect copy number alterations at the genomic level (gains or loses) in a test DNA with respect to a control. Due to its low cost, reliability and ease of implementation it has become very popular both as a research and a diagnostic tool. After hybridization, ligation and multiplex amplification of specific probes targeting different genomic regions, the probes are electrophoresed and analyzed using a DNA analyzer [##REF##12060695##13##]. Each specific probe migrates according to its molecular weight and the resulting pherograms show specific peaks that correspond to each probe (Figure ##FIG##0##1##). Relative dosage information can be obtained after the comparison of peak intensities (height or area) of the different probes between test and control samples.</p> <p>Due to the variation of PCR efficiencies across probes (due to their different size and nature) and across samples, a normalization method is needed before comparing dosage quotients. This step is crucial because the variation in experimental conditions may lead to differences of measured values between samples, thus hampering the correct interpretation of results. After normalization, which eliminates possible differences introduced during the experimental process, we aim to find biological differences in gene dosage (copy numbers). That is, the normalization tries to minimize the amount of systematic non-biological variation among samples. Different normalization methods have been used for analyzing MLPA data. The most common method divides each intensity by the sum of all intensities in each sample (see, for instance [##REF##16413013##17##] and [##REF##16542387##18##]). Other alternative approaches are based on regression methods that account for the amplification decay of larger probes. One of such regression methods uses internal control probes [##UREF##0##19##], while the another normalizes intensities based upon the statistically most probable median peak intensities using a median filter [##UREF##1##20##]. Other authors suggest to normalize peak intensities by using 4 separate peak groups according to increasing frament sizes ot the peaks (see [##REF##16240299##21##] for further details). Finally, a similar approach using the mean intensity of control probes inside a normalization group as the dividing factor was considered in [##REF##17438389##22##].</p> <p>After an intensity normalization, the next key point is the determination of genes that are significantly altered. As an example, the peak intensity of exon 13 in Figure ##FIG##0##1## seems to be lower in patient than in control indicating a deletion. Again, different approaches have been used to find probes that are altered. In this work we are considering ratio- [##REF##16542387##18##,##REF##16251890##23##] and regression-based [##UREF##0##19##,##UREF##1##20##] methods. In both approaches, the basis for the analysis is the comparison among normalized peak intensities from patient and control samples using a dosage ratio. In the ratio-based approach, deletions and duplications are given as outliers from the data set after defining a "biological" threshold (i.e. in a diploid genome, two copies exists of each gene). We assume that the simplest, most likely scenario, is an heterozygous gain or a loss of the material, under which ideal ratio values of 1.5 and 0.5 are expected, respectively. In such a simple scenario, and taking into account experimental noise, it is generally accepted in the literature that values below 0.7 and above 1.33 are indicative of loss and gain of a genetic material, respectively [##REF##12060695##13##,##REF##16542387##18##]. The regression-based approach is based on fitting several linear regression models, considering the normalized peak intensity of a given patient as the dependent variable and the normalized peak intensities of the mean control sample as independent variables. Another method based on taking into account the individual noise for each probe (i.e., standard deviation) for all control samples is considered in [##REF##16278905##24##].</p> <p>Both the ratio- and regression-based methods have general drawbacks. The ratio-based approach always considers the same gain/loss threshold without taking into account the different variability among experiments and among probes. On the other hand, although the regression-based approaches solve this problem, fitting regression models without considering that the independent variable (normalized peak intensity in controls) is also subject to error may lead to wrong conclusions. In addition, the regression method uses confidence intervals for predictions which are built under the normal theory. The normality assumption may be untenable, particularly because of the limited amount of control probes that are typically used in routine experiments (&lt; 10 probes). Finally, none of these procedures considers replicates for each individual, reducing information by averaging the replicate values.</p> <p>To overcome these difficulties, we propose alternative methods to deal with the normalization and the discovery of alterations when analyzing MLPA data. Our approach is based on both nonlinear (normalization procedure) and linear (to determine alteration) mixed models. The main advantages of our methods are: (1) the variation among individuals is modelled explicitly using a random effect; and (2) the replicates may be easily incorporated in the model. The deletions and duplications will be determined using tolerance intervals. We validate and compare our methods with the existing methods using a controlled experiment. We evaluate the presence of gains and loses in a set of DNAs from normal individuals and individuals affected by known genetic disorders. The probemix we used included probes of single copy number regions located both inside and outside of the genomic disorder regions. A total of 30 samples, including 24 controls and 6 patients, were analyzed blindly.</p>
<title>Methods</title> <p>To further illustrate the statistical methods described in this section, we are using real data provided by the United Kingdom National Genetics Reference Laboratory of Manchester [##UREF##3##26##]. In particular we are analyzing a dataset from a breast cancer study (called P002 BRCA1), in which copy number changes of different genes have previously been reported. In this example, 9 control probes and 25 analytical probes were analyzed in a total of 5 controls and 8 different test samples.</p> <title>Normalization</title> <p>Herein, we will focus on methods based on considering the total peak intensities and regression-based approaches. The first and the simplest method of normalization, determines the normalized signal of each probe in a given sample by dividing peak intensity by the sum of all peak intensity of the sample. This method was initially proposed by [##REF##12060695##13##] and is the recommended method by MRC-Holland [##UREF##4##27##] which provides commercial probes for the copy number analysis of different genome regions for diagnostic purposes. Let us begin by giving some notations to illustrate how this method works. Let <italic>H</italic><sub><italic>ipk </italic></sub>be the peak intensity for the <italic>i</italic>-th individual, <italic>i </italic>= 1, ..., <italic>C</italic>, <italic>C </italic>+ 1, ..., <italic>C </italic>+ <italic>T </italic>(<italic>C </italic>is the number of control samples and <italic>T </italic>denotes the number of test samples), the <italic>p</italic>-th probe, <italic>p </italic>= 1, ..., <italic>P</italic>, and the <italic>k</italic>-th replicate, <italic>k </italic>= 1, ..., <italic>K</italic>. The existing simple method designed for normalizing data considers the average of these replicates, , and divides each measured peak intensity by the sum of all peaks,  of that sample, i.e. .</p> <p>In addition to experimental variability, as Figure ##FIG##0##1## shows, peak intensity exhibits a negative correlation with probe size. To account for this effect, regression methods may be employed. If internal control probes are provided in the MLPA experiment (dark vertical lines in Figure ##FIG##1##2A##), a regression-based method using the control probes may also be employed to normalize the data. Different methods have been proposed. The first approach, known as "slope correction" [##UREF##0##19##], model the negative correlation between probe size and peak intensities using a linear regression model following the formula , where <italic>X</italic><sub><italic>pc </italic></sub>denotes the size of the <italic>p</italic><sub><italic>c</italic></sub>-th internal control probe for a given individual (Figure ##FIG##1##2B## solid blue line), or using a quadratic model  (Figure ##FIG##1##2B##, dotted blue line).</p> <p>Another possibility proposed by [##UREF##1##20##] is to fit an exponential decay model, . In that case, the authors propose to consider all probes <italic>p </italic>and call this method "population normalization". Notice that in that case <italic>β </italic>encodes the <italic>rate </italic>of descend. [##UREF##1##20##] argued that the trend of peak intensities varies greatly between samples regarding the control probes. This is the reason why they propose to consider all probes (both analytical and control ones) after applying a median filter approach to remove outliers. However this method has important inconveniences. The smoothed data, obtained after applying the median filter, violates some of the assumptions of nonlinear regression such as that residuals are no longer independent or that errors are not Gaussian, among others [##UREF##5##28##].</p> <p>Neither <italic>slope correction nor population normalization </italic>approaches consider replicates, e.g. they use  instead of <italic>H</italic><sub><italic>ipk</italic></sub>. In addition, these methods do not model the variability among individuals either. That is, using these methods, the authors allow only <italic>β </italic>to vary among individuals fitting different models for each one. However, the pattern of decay may vary highly between subjects, which is not captured by varying <italic>β </italic>'s only. As an example, the height of the peak intensity for the first control probe (top parameter in Figure ##FIG##1##2C##) and the asymptotic value for larger probe sizes (bottom parameter in Figure ##FIG##1##2C##) are very different among test samples, as we observe in Figure ##FIG##2##3##.</p> <p>To circumvent these problems, we propose the following nonlinear mixed model to normalize the data which includes replicates and variability among individuals:</p> <p></p> <p>where <italic>X</italic><sub><italic>p </italic></sub>denotes the size for the <italic>p</italic>-th probe (note that it only depends on the probe since it is the same for each individual and each replicate) <italic>φ</italic><sub><italic>i</italic>1 </sub>is the maximum peak intensity for control probes, <italic>φ</italic><sub>2<italic>i</italic></sub>, is the asymptotic peak intensity, and <italic>φ</italic><sub>3<italic>i </italic></sub>is the reciprocal of the rate decay constant. The term off<sub><italic>i </italic></sub>provides a more stable parametrization for the data and corresponds to the average value of the peak size for the first control probe, off<sub><italic>i </italic></sub>= , <italic>i </italic>= 1, ..., <italic>C</italic>, <italic>C </italic>+ 1, ..., <italic>C </italic>+ <italic>T</italic>. The fixed effects, <italic>β </italic>represent the population average of subject-specific parameters, <italic>φ</italic><sub><italic>i</italic></sub>, and the random effects, <italic>b</italic><sub><italic>i</italic></sub>, represents the deviations of the <italic>φ</italic><sub><italic>i</italic></sub>'s from their population average. The random effects are assumed to be independent and the within-group errors <italic>ε</italic><sub><italic>ipk </italic></sub>are assumed to be independent for different <italic>i</italic>, <italic>p </italic>and to be independent of the random effects. The model parameters are fitted by maximizing the restricted log-likelihood (RMLE) of the data using the R library nlme [##UREF##6##29##].</p> <p>After estimating this generalized exponential model, the normalization procedure is performed by dividing the peak intensities, <italic>H</italic><sub><italic>ipk</italic></sub>, by the regression estimate, , obtaining a normalized, size-adjusted peak intensity for every probe, <italic>Y</italic><sub><italic>ipk </italic></sub>= <italic>H</italic><sub><italic>ipk</italic></sub> (Figure ##FIG##1##2D##). For the control samples we consider a unique regression line corresponding to the average of all control samples. Nonetheless, for the test samples the normalization is performed individually for each subject.</p> <title>Statistical significance of dosage ratios</title> <title>Ratio approach</title> <p>The normalized peak intensities may be directly compared between control and test samples peak by peak. The first method for determining genes that are altered (gains or loses) is based on a threshold rule. The probes whose dosage ratios are outside of the threshold lines are considered altered. The thresholds are normally 0.7 for deletions (note that 0.5 would correspond to a half reduction in the dosage but this value is considered too stringent by some researchers) and 1.33 for duplications, although other cut-points have been used in other studies [##REF##16542387##18##].</p> <title>Iterative regression</title> <p>The main disadvantage of defining a threshold is that the specific variability of each experiment is not considered. Some authors have proposed alternative methods using linear regression models. [##UREF##0##19##] stated that, after a square-root transformation of the data, the plot of normalized peak intensities for test samples versus control samples should be around the diagonal if no dosage imbalances are present; otherwise, probes above or below would correspond to duplications or deletions, respectively. Their proposal was to fit a linear regression without intercept using only control probes and use the standard error of the regression to determine confidence limits for outlier detections. As in the case of normalization procedures, however, a small number of control probes is a limiting factor for adopting this approach since, as the author stated, the standard error of the regression underestimates the normal variability. Thus, [##UREF##0##19##] proposed an iterative regression algorithm for this problem. Starting with the control probes, the linear regression is fitted. Then, a prediction interval is computed for every probe. If the normalized peak intensity of a non-control probe falls into this interval, the probe is considered to belong to the non-altered population and it is added to the control probes in the next step. After a few iterations, a final model is reached including only those probes that are considered as preserved. In practice, those probes that are not included in this final model are considered to be either deleted or duplicated. The authors call this method Regression-Enhanced MLPA analysis (REX-MLPA). Figure ##FIG##3##4## illustrates the confidence limits of the final regression model. [##UREF##1##20##] proposed a similar approach but starting with all probes and retaining and rejecting points at each iteration with a given level of confidence.</p> <p>The REX-MLPA approach assumes that the X variable (e.g. squared-root of mean normalized peak intensities in controls) in the linear regression model is measured without error. This assumption is not tenable since this variable corresponds to the normalized peak intensity measured in multiple control samples. Violations of this assumption are a serious problem when trying to accurately predict Y (squared-root of normalized peak intensities in test sample) from X.</p> <title>The probe-specific mixture model</title> <p>The REX-MLPA approach has the limitations discussed above, specifically, the regression models are fitted without considering that the independent variable (normalized peak intensity in controls) is also subject to error. In addition, confidence intervals for predictions are built assuming the normal theory. This assumption may not hold when only control probes are considered since we typically have no more than 10 probes. Replicates for each individual are also not considered.</p> <p>Based on these limitations, we propose the following model to compare controls with a given test sample:</p> <p></p> <p>where <italic>μ </italic>is the mean normalized peak intensity across individuals and probes, <italic>α</italic><sub><italic>t </italic></sub>is a fixed effect representing a different mean for all controls (t = 0) and a given test sample (t = 1), <italic>β</italic><sub><italic>p </italic></sub>is a random effect for the <italic>p</italic>th probe, <italic>γ</italic><sub><italic>tp </italic></sub>is is the probe-test interaction random effect, and <italic>ε</italic><sub><italic>tpk </italic></sub>is a random error variable. The index <italic>p </italic>denotes the probe, and <italic>K </italic>is the number of replicates. We assume that <italic>β</italic><sub><italic>p</italic></sub>, <italic>γ</italic><sub><italic>tp</italic></sub>, and <italic>ε</italic><sub><italic>tpk </italic></sub>are independent, homocedastic, and normally distributed random variables with mean zero. The variances are denoted by  for the, <italic>β</italic><sub><italic>p</italic></sub>, or the "between-probe" variability,  for <italic>γ</italic><sub><italic>tp</italic></sub>, or "probe-test" variability, and  for the <italic>ε</italic><sub><italic>tpk</italic></sub>, or "within-probe" variability. Considering the assumptions of the proposed model, we need to verify that the within-groups errors are independent and identically normally distributed, with mean zero and variance , and that they are also independent of the random effect. It can be easily checked using different plots. This assumption makes sense for biologists since probes that behave abnormally in replicate studies or when hybridised onto DNAs that have a different quality are typically removed or redesigned to ensure assay robutness and reliability. The model parameters are fitted by maximizing the restricted log-likelihood (RMLE) using the R library lme (Pinheiro and Bates, 2000).</p> <title>The criterion</title> <p>For the sake of simplicity, we take control samples as the reference, i.e. <italic>α</italic><sub>0 </sub>= 0 and <italic>α</italic><sub>1 </sub>= <italic>α</italic>, and <italic>γ</italic><sub>0<italic>p </italic></sub>= 0. Therefore, <italic>α </italic>may be interpreted as an average deviation between test and control samples across all probes, and <italic>γ</italic><sub>1<italic>p </italic></sub>as the deviation control-test for the <italic>p</italic>-th probe.</p> <p>The criterion to determine those probes that are altered for each test sample is based on the observed differences among controls and its differences with controls samples. Conditioned to <italic>p</italic>-th probe, the difference between two control measurements is distributed as</p> <p></p> <p>Conversely, the difference between a test, <italic>i'</italic>, and a control sample, conditioned to a given probe, is distributed as</p> <p></p> <p>where  corresponds to the variance of the error defined in equation (2). Then, we consider that a probe is not altered if the mean of the difference between test and control samples is included within a probability interval defined by differences between control samples distribution. Hence, the criterion consists in checking whether the difference between test and control samples is greater than that expected between two control samples. That is, whether the difference is included in the interval</p> <p></p> <p>where Z<sub>1-<italic>κ</italic>/2 </sub>is the 1 - <italic>κ</italic>/2 percentile from the standard normal distribution. Hence, 1 - <italic>κ </italic>is the proportion of difference control-control that the difference test-control must exceed to declare that a probe is altered. Nonetheless, when using sample data to estimate such intervals, the uncertainty associated with the estimation estimating them has to be considered. Therefore, we propose to estimate the interval (5) through a tolerance interval over the control-control differences distribution. To build such intervals, it is assumed that , where <italic>ν </italic>are the residual degrees of freedom, and  is the REML estimator of . Therefore, the tolerance interval that contains (1 - <italic>κ</italic>)% proportion of the control-control differences estimated with a confidence of (1 - <italic>α</italic>)% is</p> <p></p> <p>Once the limits of the criterion are established, the mean difference test-control, <italic>δ</italic><sub><italic>p </italic></sub>for each probe, <italic>p</italic>, is estimated using the estimate of <italic>α </italic>+ <italic>γ</italic><sub><italic>i'p</italic></sub>, where <italic>γ</italic><sub><italic>p</italic>1 </sub>is estimated through the best linear unbiased predictor (BLUP). Then a probe for a given test sample is considered as altered when <italic>δ</italic><sub><italic>p </italic></sub>does not fall into (6). Figure ##FIG##4##5## shows those probes that are imbalanced for each individual considering the threshold obtained using tolerance limits obtained using equation (6). In this example those limits are 0.78 for deletions and 1.29 for duplications.</p>
<title>Results</title> <title>Simulation study</title> <p>In this section, we investigate the performance of the two existing methods and our proposed approach through a number of simulations. We simulated an hypothetical assay with 45 probes. We also simulated internal control probes (i.e. non-altered) to be used in the iterative regression approach by setting <italic>σ</italic><sub><italic>γ </italic></sub>= 0. Two different scenarios were simulated in order to investigate the performance of REX-MLPA approach given the number of control probes: 10% and 20% of the 45 probes. Finally, the ratios of the rest of probes were simulated as altered (gains or loses) depending on a different percentages: 5%, 10%, 20% and 50% of the probes. These values were generated from a normal distribution under different scenarios varying between probe variability (<italic>σ</italic><sub><italic>β </italic></sub>∈ {0.2, 0.4, 0.6}), probe-test variability (<italic>σ</italic><sub><italic>γ </italic></sub>∈ {0.2, 0.4, 0.6, 1.0, 1.5, 4.0}) and within-probe or random variability (<italic>σ</italic><sub><italic>ε </italic></sub>∈ {0.05, 0.08}). We simulated 3 replicates of each peak intensity for test and control samples. Note that the average of these replicates were considered in the case of using threshold and REX-MLPA approaches. We summarized our simulations by computing the mean number of altered probed simulated in each run, the empirical type-I error and the empirical power to detect gains and loses. These results are based on 1,000 simulations.</p> <p>Tables ##TAB##0##1## and ##TAB##1##2## show the results for the case of having 10% of true altered probes when 10% and 20% of probes have been generated as a internal control probes, respectively. The results for the REX-MLPA method were based on building confidence limits at 99% at each iteration. Regarding the empirical type-I error, the REX-MLPA approach usually overestimate the expected 5%, while the threshold method clearly underestimates type-I error since in all simulations the simulated type-I error was closer to 0% (results omitted in the tables). On the other hand, the probe-specific model shows a good performance. As expected, the power of the three methods increases when the probe-test variability increases. That is, the power depends on the size effect of those altered probes. In general, the probe-specific mixed-model outperforms both REX-MLPA and threshold approaches. The threshold approach is clearly underpowered for those cases in which the magnitude of the effect is not large enough. Finally, the same pattern is observed when the random variability is increased (e.g. within-probe variability) but showing lower power. The simulation results for the case of having 5%, 20% and 50% of true altered probes showed similar behaviour and they can be found in the Additional file ##SUPPL##0##1##.</p> <title>Validation study</title> <p>In order to validate our proposed method, and compare its performace with those of the existing methods, we have designed a controlled MLPA probemix in which we have included 34 different oligonucleotides, corresponding to 16 different targeted regions with single copy number in normal individuals. Eight of the targeted regions are variable in copy number in individuals suffering from different genomic disorders (see Table ##TAB##2##3##). A total of 30 DNAs from different individuals were assayed in triplicate in the study. 24 of the DNAs came from unrelated HapMap individuals (i.e. normal general population), 2 DNAs from Prader-Willi syndrome patients (15q11-q13 deletion), 2 DNAs from DiGeorge syndrome patients (22q11 deletion) and 2 DNAs from autistic patients with a duplication in 10q region. After hybridization, ligation and PCR, probes corresponding to the different individuals were electrophoresed and peak intensities were recorded and blindly analyzed using the different methods. MLPA hybridization, ligation and PCR with universal MLPA primers was performed as described elsewhere [##REF##12060695##13##]. PCR products were loaded on a capillar DNA analyzer and electrophoresed. Genescan software was used to analyze the runs and to retrieve peak intensities corresponding to each probe in the different samples.</p> <title>Comparison of methods of copy number estimation</title> <p>Before analyzing the data we illustrate that our proposed method based on a mixed-model can be applied. Supplementary Figure S1 shows the residual error for each individual and probe indicating that the residuals are centered at zero and that the variability does not change with individuals. Table ##TAB##3##4## shows the results obtained from the 9 test samples used in the validation study using the threshold, iterative regression and proposed probe-specific mixed-model method. For the two patients with Prader-Willi syndrome, we observed that all methods were able to detect the three deletions in 15q11-q13 region. The same conclusion was observed regarding patients diagnosed with DiGeorge syndrome. Considering patients with Autism, all methods detected gains in the 10q region. Only in one case threshold method had a false negative result and the REX-MLPA method indicated a deletion in the UBE3A, a false positive finding. These results clearly agree with our simulation studies in which we found that iterative regression has a increased type-I error rate and the threshold approach a low power. For instance, the ratio between the second individual diagnosed with Autism and the mean controls was 1.27 at probe named ZWINT, so the threshold method did not indicate that this was a probe altered in this individual. This result is due to the fact that the ratio between the probe-test variability and the random error variability was around 3 (3.09 = ) and the simulation study showed that in this case the power of threshold approach is very poor. Finally, regarding HapMap individuals who were considered as negative controls, we observed that any method found a false positive result.</p>
<title>Conclusion</title> <p>The MLPA method is a potentially useful semi-quantitative method to detect copy number alterations in targeted regions obtained after performing genome-wide screening using comparative methods usually based on genomic hybridization. MLPA is based on comparing peak intensities of the different probes between test and control samples through its ratio. Before determining the statistical significance of dosage ratios, a normalization procedure is needed to control variation of PCR efficiencies across probes. We have proposed a non-linear mixed-model to perform such normalization, and two other existing methods have been discussed. One of them is based on considering the sum of peak intensities in each sample, while another uses internal control probes to fit a regression line which is used as a reference.</p> <p>So far, the widely used method to detect statistical significance of dosage ratio is to consider as altered those probes that are outside a given threshold. In this paper, we have described a new approach to detect altered probes based on a mixed-model. The threshold is then established by using a tolerance interval which accommodates the specific random error variability observed in each test sample. We have also discussed another approach based on an iterative regression. Through simulation studies and a controlled MLPA assay, we have shown that our proposed method outperforms the existing ones. Another important advantage by comparing our approach and the REX-MLPA method is that our method can be used even when no internal control probes are available in the assay.</p> <p>A novelty of our algorithm is that it can handle information from replicate experiments. Replicates of MLPA assays are typically not necessary for commercial kits since in most cases more than one probe is used to interrogate the same region and concordance between those probes is considered to provide enough reliability of the existence of a copy number variation. Nevertheless, in cases where a single probe per region is placed in the assay, it is definitely useful and desirable (if not necessary) to have replicates that help minimise possible variations arising during the long MLPA protocol and to ensure reliability of the findings. Thus, having a method that is able to handle such information is very useful and provides extra statistical robustness to the findings.</p> <p>The use of control probes as references for data normalization is one of the methods recommended by the manufacturer. Commercial and home-made MLPA assays are clearly different in terms of quality of the resulting pherograms, mainly due to the process of producing the specific probes that target regions of interest. While commercial applications produce very high, sharp and neat peaks, home-made assays typically do not perform that well. Nevertheless, in both cases, the use of control probes is always desirable from a statistical point of view and it is so regardless of the number of probes that target a determined loci. In exploratory experiments with MLPA designs targeting extremely variable and polymorphic CNVs, adding control probes whose behaviour is proven unvariable and can help in assessing reliability of single experiments. In such cases, we understand that the inclusion of evenly spaced control probes is a good practice. One key point is to know the number of replicates that should be included in a MLPA assay. Our simulations indicates that 3 replicates was enough to reach a power of 90% keeping a value of alfa around 5% in those combinations where the probe-test variability was over 8 times greater than that of random error (combinations where <italic>σ</italic><sub><italic>γ </italic></sub>= 0.4 and <italic>σ</italic><sub><italic>ε </italic></sub>= 0.05). Thus, if the probe-test versus random error variability ratio was lower, a lesser power would be expected. In that case the researcher could improve the expected power by increasing the number of replicates.</p> <p>Although our proposed method outperforms other existing ones, it has to be kept in mind that in the case of having probes targeting regions that are highly variable among the population (i.e. in an extreme case where half individuals are normal homozygous and half are homozygously deleted), the tolerance intervals calculated for those probes might collect all existing variability and become very broad, thus making it impossible to distinguish the existence of copy number alteration. Nevertheless, the other existing methods would also show inconsistent copy number calls when analysing such regions. In statistical terms, this means that we should have homogeneity of the random error variance through the probes. This is the reason why a visual inspection of the pherograms as well as manual curation of the results might be the most reliable way to proceed. However, we can easily accomodate cases where the homogeneity of random error variance assumption cannot be assumed by introducing an interaction probe-random error allowing a different random error variance for each probe. In this case, different tolerance bands would be therefore estimated for each probe.</p> <p>In conclusion, we have proposed a mixed-model which is able to determine thresholds to decide whether a region is altered. These threholds are specific for each individual, incorporating experimental variability, resulting in improved sensitivity and specificity as the examples with real data have revealed. An R language package for the three approaches discussed in this paper and data analyzed will be freely available at our web page [##UREF##2##25##].</p>
<p>This is an Open Access article distributed under the terms of the Creative Commons Attribution License (<ext-link ext-link-type="uri" xlink:href="http://creativecommons.org/licenses/by/2.0"/>), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.</p> <title>Background</title> <p>MLPA method is a potentially useful semi-quantitative method to detect copy number alterations in targeted regions. In this paper, we propose a method for the normalization procedure based on a non-linear mixed-model, as well as a new approach for determining the statistical significance of altered probes based on linear mixed-model. This method establishes a threshold by using different tolerance intervals that accommodates the specific random error variability observed in each test sample.</p> <title>Results</title> <p>Through simulation studies we have shown that our proposed method outperforms two existing methods that are based on simple threshold rules or iterative regression. We have illustrated the method using a controlled MLPA assay in which targeted regions are variable in copy number in individuals suffering from different disorders such as Prader-Willi, DiGeorge or Autism showing the best performace.</p> <title>Conclusion</title> <p>Using the proposed mixed-model, we are able to determine thresholds to decide whether a region is altered. These threholds are specific for each individual, incorporating experimental variability, resulting in improved sensitivity and specificity as the examples with real data have revealed.</p>
<title>Abbreviations</title> <p>CNV: Copy Number Variant, MLPA: Multiplex Ligation-dependent Probe Amplification, DNA: Deoxyribonucleic acid, aCGH: array-based Comparative Genomic Hybridization, QMPSF: Quantitative Multiplex PCR of Short Fluorescent, PCR: Polymerase Chain Reaction, MAPH: Multiplex Amplifiable Probe Hybridization, BRCA1: Breast Cancer 1, RMLE: Restricted maximum log-likelihood, REX-MLPA: Regression-Enhanced MLPA, BLUP: Best linear unbiased predictor, PEMM: probe-specific mixture model.</p> <title>Authors' contributions</title> <p>JRG and JLC developed the new statistical methods and performed the simulation studies. JRG wrote the R functions and the main text of the manuscript. LA and SV performed the MPLA experiment, interpreted the results and tested the programs functionality. LA, LJ and YY proposed abundant suggestions for improving the implementation of the models and participated in the design and discussion of this study. YY and XE reviewed the paper and revised its framework. All authors have read, and approved the final manuscript.</p> <title>Supplementary Material</title>
<title>Acknowledgements</title> <p>We would like to thank the two reviewers for their constructive comments, which have led to improvements in the manuscript. The laboratory of X.E is supported by the Departament dEducació i Universitats and the Departament de Salut of the Catalan Autonomous Government ("Generalitat de Catalunya"); the Ministry of Health, and the Ministry of Education and Science of the Spanish Government; and the European Union Sixth Framework Programme. CIBERESP is supported by the Instituto de Salud Carlos III of the Ministry of Health of the Spanish Government.</p>
<fig position="float" id="F1"><label>Figure 1</label><caption><p><bold>Example of a deletion in an MLPA assay</bold>. Panels A and B show pherograms corresponding to the electrophoresis of an MLPA assay. In the Y-asix are depicted the intensity signals (peak heights) for each probe that are depicted in the X-axis according to their length (probe size). Peaks marked with a <bold>C </bold>correspond to control probes and peaks numbered from 1 through 24 correspond to region-specific probes. Panel A corresponds to a normal individual, while panel B corresponds to an individual with a deletion at probe #13 as visible by the reduced peak intensity in this pherogram.</p></caption></fig> <fig position="float" id="F2"><label>Figure 2</label><caption><p><bold>MLPA normalization procedure</bold>. Panel A shows the peak intensities depending on probe sizes. Panel B shows the regression estimates for different methods: linear regression model using control probes (solid blue line), quadratic model using control probes (dotted blue line), nonlinear mixed model (red lines) using control probes (solid line) or median filter approach (dotted line). Panel C illustrates the parameters involved in the nonlinear mixed model. Panel D shows the size-adjusted normalized peak intensities prepared to compute the dosage quotient. In all panels dark lines represents control probes while light lines are for analytical probes.</p></caption></fig> <fig position="float" id="F3"><label>Figure 3</label><caption><p><bold>Regression estimates for the test samples from BRCA1 data set</bold>. Regression estimates for each of the 8 test samples given in the example provided by the NGRL-Manchester called P002 BRCA1. Red lines are estimated using the nonlinear mixed model. The solid lines are estimated using control probes, while dotted lines are obtained after using the median filter approach. The dotted blue lines are showing the regression estimates using quadratic model. These regression lines are then used to normalize the peak intensities.</p></caption></fig> <fig position="float" id="F4"><label>Figure 4</label><caption><p><bold>Normal control comparison procedure using iterative regression</bold>. Plot of normalized peak intensities in controls (calculated as the mean of all controls) against normalized peak intensities in test sample number 1 using iterative regression procedure. Dotted line correspond to linear regression among those probes that are considered non altered after applying the iterative procedure. Solid lines determine upper and lower boundaries which are used to indicate whether or not a given probe is a duplication or a deletion, respectively.</p></caption></fig> <fig position="float" id="F5"><label>Figure 5</label><caption><p><bold>Normal control comparison procedure using tolerance interval criteria</bold>. Plot of ratio between normalized peak intensities in controls against normalized peak intensities in each of the 8 test samples (vertical lines). Horizontal gray lines indicate lower and upper boundaries obtained using the linear mixed model and tolerance interval criteria. Vertical green lines are indicating that those probes are duplicated for the corresponding probe.</p></caption></fig>
<table-wrap position="float" id="T1"><label>Table 1</label><caption><p>Empirical type-I error and power obtained in 1,000 simulations using the three different approaches: REX (iterative regression), PEMM (probe-specific mixture model) and threshold.</p></caption><table frame="hsides" rules="groups"><thead><tr><td/><td/><td/><td/><td align="center" colspan="2">Type-I error</td><td align="center" colspan="3">Power (gains)</td><td align="center" colspan="3">Power (loses)</td></tr><tr><td colspan="4"/><td colspan="2"><hr/></td><td colspan="3"><hr/></td><td colspan="3"><hr/></td></tr><tr><td align="center"><italic>σ</italic><sub><italic>β</italic></sub></td><td align="center"><italic>σ</italic><sub><italic>γ</italic></sub></td><td align="center"><italic>σ</italic><sub><italic>ε</italic></sub></td><td align="center"></td><td align="center">REX</td><td align="center">PEMM</td><td align="center">REX</td><td align="center">PEMM</td><td align="center">thres</td><td align="center">REX</td><td align="center">PEMM</td><td align="center">thres</td></tr></thead><tbody><tr><td align="center">0.2</td><td align="center">0.2</td><td align="center">0.05</td><td align="center">3.7</td><td align="center">0.073</td><td align="center">0.031</td><td align="center">0.815</td><td align="center">0.745</td><td align="center">0.000</td><td align="center">0.819</td><td align="center">0.751</td><td align="center">0.000</td></tr><tr><td align="center">0.2</td><td align="center">0.2</td><td align="center">0.08</td><td align="center">3.6</td><td align="center">0.071</td><td align="center">0.027</td><td align="center">0.689</td><td align="center">0.546</td><td align="center">0.000</td><td align="center">0.681</td><td align="center">0.541</td><td align="center">0.000</td></tr><tr><td align="center">0.2</td><td align="center">0.4</td><td align="center">0.05</td><td align="center">3.5</td><td align="center">0.081</td><td align="center">0.044</td><td align="center">0.923</td><td align="center">0.891</td><td align="center">0.024</td><td align="center">0.930</td><td align="center">0.908</td><td align="center">0.029</td></tr><tr><td align="center">0.2</td><td align="center">0.4</td><td align="center">0.08</td><td align="center">3.5</td><td align="center">0.079</td><td align="center">0.038</td><td align="center">0.870</td><td align="center">0.811</td><td align="center">0.028</td><td align="center">0.863</td><td align="center">0.814</td><td align="center">0.057</td></tr><tr><td align="center">0.2</td><td align="center">0.6</td><td align="center">0.05</td><td align="center">3.6</td><td align="center">0.081</td><td align="center">0.052</td><td align="center">0.965</td><td align="center">0.955</td><td align="center">0.131</td><td align="center">0.956</td><td align="center">0.956</td><td align="center">0.161</td></tr><tr><td align="center">0.2</td><td align="center">0.6</td><td align="center">0.08</td><td align="center">3.6</td><td align="center">0.086</td><td align="center">0.044</td><td align="center">0.922</td><td align="center">0.889</td><td align="center">0.140</td><td align="center">0.922</td><td align="center">0.900</td><td align="center">0.184</td></tr><tr><td align="center">0.2</td><td align="center">1.0</td><td align="center">0.05</td><td align="center">3.6</td><td align="center">0.081</td><td align="center">0.057</td><td align="center">0.980</td><td align="center">0.983</td><td align="center">0.360</td><td align="center">0.977</td><td align="center">0.972</td><td align="center">0.411</td></tr><tr><td align="center">0.2</td><td align="center">1.0</td><td align="center">0.08</td><td align="center">3.6</td><td align="center">0.085</td><td align="center">0.052</td><td align="center">0.955</td><td align="center">0.954</td><td align="center">0.396</td><td align="center">0.947</td><td align="center">0.939</td><td align="center">0.449</td></tr><tr><td align="center">0.2</td><td align="center">1.5</td><td align="center">0.05</td><td align="center">3.6</td><td align="center">0.081</td><td align="center">0.060</td><td align="center">0.982</td><td align="center">0.978</td><td align="center">0.536</td><td align="center">0.981</td><td align="center">0.983</td><td align="center">0.582</td></tr><tr><td align="center">0.2</td><td align="center">1.5</td><td align="center">0.08</td><td align="center">3.6</td><td align="center">0.084</td><td align="center">0.056</td><td align="center">0.981</td><td align="center">0.986</td><td align="center">0.623</td><td align="center">0.980</td><td align="center">0.975</td><td align="center">0.630</td></tr><tr><td align="center">0.2</td><td align="center">4.0</td><td align="center">0.05</td><td align="center">3.5</td><td align="center">0.078</td><td align="center">0.064</td><td align="center">0.995</td><td align="center">0.993</td><td align="center">0.815</td><td align="center">0.993</td><td align="center">0.994</td><td align="center">0.840</td></tr><tr><td align="center">0.2</td><td align="center">4.0</td><td align="center">0.08</td><td align="center">3.5</td><td align="center">0.079</td><td align="center">0.060</td><td align="center">0.996</td><td align="center">0.993</td><td align="center">0.848</td><td align="center">0.993</td><td align="center">0.995</td><td align="center">0.863</td></tr><tr><td colspan="12"><hr/></td></tr><tr><td align="center">0.4</td><td align="center">0.2</td><td align="center">0.05</td><td align="center">3.7</td><td align="center">0.079</td><td align="center">0.032</td><td align="center">0.822</td><td align="center">0.729</td><td align="center">0.000</td><td align="center">0.827</td><td align="center">0.737</td><td align="center">0.001</td></tr><tr><td align="center">0.4</td><td align="center">0.2</td><td align="center">0.08</td><td align="center">3.7</td><td align="center">0.070</td><td align="center">0.027</td><td align="center">0.690</td><td align="center">0.522</td><td align="center">0.000</td><td align="center">0.697</td><td align="center">0.537</td><td align="center">0.000</td></tr><tr><td align="center">0.4</td><td align="center">0.4</td><td align="center">0.05</td><td align="center">3.7</td><td align="center">0.081</td><td align="center">0.046</td><td align="center">0.918</td><td align="center">0.897</td><td align="center">0.027</td><td align="center">0.926</td><td align="center">0.905</td><td align="center">0.044</td></tr><tr><td align="center">0.4</td><td align="center">0.4</td><td align="center">0.08</td><td align="center">3.7</td><td align="center">0.078</td><td align="center">0.038</td><td align="center">0.878</td><td align="center">0.807</td><td align="center">0.035</td><td align="center">0.871</td><td align="center">0.825</td><td align="center">0.069</td></tr><tr><td align="center">0.4</td><td align="center">0.6</td><td align="center">0.05</td><td align="center">3.6</td><td align="center">0.080</td><td align="center">0.052</td><td align="center">0.958</td><td align="center">0.949</td><td align="center">0.139</td><td align="center">0.952</td><td align="center">0.942</td><td align="center">0.191</td></tr><tr><td align="center">0.4</td><td align="center">0.6</td><td align="center">0.08</td><td align="center">3.5</td><td align="center">0.083</td><td align="center">0.044</td><td align="center">0.926</td><td align="center">0.903</td><td align="center">0.150</td><td align="center">0.913</td><td align="center">0.881</td><td align="center">0.206</td></tr><tr><td align="center">0.4</td><td align="center">1.0</td><td align="center">0.05</td><td align="center">3.7</td><td align="center">0.080</td><td align="center">0.057</td><td align="center">0.981</td><td align="center">0.982</td><td align="center">0.353</td><td align="center">0.973</td><td align="center">0.977</td><td align="center">0.425</td></tr><tr><td align="center">0.4</td><td align="center">1.0</td><td align="center">0.08</td><td align="center">3.5</td><td align="center">0.081</td><td align="center">0.050</td><td align="center">0.949</td><td align="center">0.953</td><td align="center">0.392</td><td align="center">0.950</td><td align="center">0.934</td><td align="center">0.424</td></tr><tr><td align="center">0.4</td><td align="center">1.5</td><td align="center">0.05</td><td align="center">3.6</td><td align="center">0.083</td><td align="center">0.063</td><td align="center">0.989</td><td align="center">0.985</td><td align="center">0.553</td><td align="center">0.987</td><td align="center">0.989</td><td align="center">0.613</td></tr><tr><td align="center">0.4</td><td align="center">1.5</td><td align="center">0.08</td><td align="center">3.6</td><td align="center">0.080</td><td align="center">0.056</td><td align="center">0.976</td><td align="center">0.973</td><td align="center">0.568</td><td align="center">0.966</td><td align="center">0.967</td><td align="center">0.603</td></tr><tr><td align="center">0.4</td><td align="center">4.0</td><td align="center">0.05</td><td align="center">3.6</td><td align="center">0.085</td><td align="center">0.066</td><td align="center">0.997</td><td align="center">0.996</td><td align="center">0.856</td><td align="center">0.995</td><td align="center">0.995</td><td align="center">0.855</td></tr><tr><td align="center">0.4</td><td align="center">4.0</td><td align="center">0.08</td><td align="center">3.6</td><td align="center">0.086</td><td align="center">0.062</td><td align="center">0.994</td><td align="center">0.991</td><td align="center">0.841</td><td align="center">0.989</td><td align="center">0.990</td><td align="center">0.864</td></tr><tr><td colspan="12"><hr/></td></tr><tr><td align="center">0.6</td><td align="center">0.2</td><td align="center">0.05</td><td align="center">3.6</td><td align="center">0.075</td><td align="center">0.031</td><td align="center">0.828</td><td align="center">0.740</td><td align="center">0.000</td><td align="center">0.824</td><td align="center">0.741</td><td align="center">0.001</td></tr><tr><td align="center">0.6</td><td align="center">0.2</td><td align="center">0.08</td><td align="center">3.6</td><td align="center">0.071</td><td align="center">0.027</td><td align="center">0.701</td><td align="center">0.529</td><td align="center">0.000</td><td align="center">0.690</td><td align="center">0.520</td><td align="center">0.002</td></tr><tr><td align="center">0.6</td><td align="center">0.4</td><td align="center">0.05</td><td align="center">3.7</td><td align="center">0.082</td><td align="center">0.047</td><td align="center">0.932</td><td align="center">0.916</td><td align="center">0.032</td><td align="center">0.914</td><td align="center">0.902</td><td align="center">0.063</td></tr><tr><td align="center">0.6</td><td align="center">0.4</td><td align="center">0.08</td><td align="center">3.6</td><td align="center">0.077</td><td align="center">0.037</td><td align="center">0.864</td><td align="center">0.815</td><td align="center">0.041</td><td align="center">0.874</td><td align="center">0.835</td><td align="center">0.070</td></tr><tr><td align="center">0.6</td><td align="center">0.6</td><td align="center">0.05</td><td align="center">3.7</td><td align="center">0.082</td><td align="center">0.051</td><td align="center">0.964</td><td align="center">0.949</td><td align="center">0.142</td><td align="center">0.955</td><td align="center">0.944</td><td align="center">0.174</td></tr><tr><td align="center">0.6</td><td align="center">0.6</td><td align="center">0.08</td><td align="center">3.6</td><td align="center">0.078</td><td align="center">0.044</td><td align="center">0.924</td><td align="center">0.885</td><td align="center">0.177</td><td align="center">0.918</td><td align="center">0.887</td><td align="center">0.206</td></tr><tr><td align="center">0.6</td><td align="center">1.0</td><td align="center">0.05</td><td align="center">3.6</td><td align="center">0.081</td><td align="center">0.058</td><td align="center">0.977</td><td align="center">0.970</td><td align="center">0.388</td><td align="center">0.972</td><td align="center">0.974</td><td align="center">0.435</td></tr><tr><td align="center">0.6</td><td align="center">1.0</td><td align="center">0.08</td><td align="center">3.6</td><td align="center">0.080</td><td align="center">0.050</td><td align="center">0.946</td><td align="center">0.935</td><td align="center">0.380</td><td align="center">0.960</td><td align="center">0.945</td><td align="center">0.451</td></tr><tr><td align="center">0.6</td><td align="center">1.5</td><td align="center">0.05</td><td align="center">3.6</td><td align="center">0.087</td><td align="center">0.062</td><td align="center">0.982</td><td align="center">0.981</td><td align="center">0.556</td><td align="center">0.982</td><td align="center">0.985</td><td align="center">0.590</td></tr><tr><td align="center">0.6</td><td align="center">1.5</td><td align="center">0.08</td><td align="center">3.6</td><td align="center">0.080</td><td align="center">0.056</td><td align="center">0.972</td><td align="center">0.972</td><td align="center">0.597</td><td align="center">0.975</td><td align="center">0.973</td><td align="center">0.622</td></tr><tr><td align="center">0.6</td><td align="center">4.0</td><td align="center">0.05</td><td align="center">3.5</td><td align="center">0.080</td><td align="center">0.064</td><td align="center">0.995</td><td align="center">0.995</td><td align="center">0.837</td><td align="center">0.994</td><td align="center">0.993</td><td align="center">0.855</td></tr><tr><td align="center">0.6</td><td align="center">4.0</td><td align="center">0.08</td><td align="center">3.6</td><td align="center">0.082</td><td align="center">0.064</td><td align="center">0.988</td><td align="center">0.991</td><td align="center">0.850</td><td align="center">0.990</td><td align="center">0.990</td><td align="center">0.857</td></tr></tbody></table></table-wrap> <table-wrap position="float" id="T2"><label>Table 2</label><caption><p>Empirical type-I error and power obtained in 1,000 simulations using the three different approaches: REX (iterative regression), PEMM (probe-specific mixture model) and threshold.</p></caption><table frame="hsides" rules="groups"><thead><tr><td/><td/><td/><td/><td align="center" colspan="2">Type-I error</td><td align="center" colspan="3">Power (gains)</td><td align="center" colspan="3">Power (loses)</td></tr><tr><td colspan="4"/><td colspan="2"><hr/></td><td colspan="3"><hr/></td><td colspan="3"><hr/></td></tr><tr><td align="center"><italic>σ</italic><sub><italic>β</italic></sub></td><td align="center"><italic>σ</italic><sub><italic>γ</italic></sub></td><td align="center"><italic>σ</italic><sub><italic>ε</italic></sub></td><td align="center"></td><td align="center">REX</td><td align="center">PEMM</td><td align="center">REX</td><td align="center">PEMM</td><td align="center">thres</td><td align="center">REX</td><td align="center">PEMM</td><td align="center">thres</td></tr></thead><tbody><tr><td align="center">0.2</td><td align="center">0.2</td><td align="center">0.05</td><td align="center">3.9</td><td align="center">0.044</td><td align="center">0.034</td><td align="center">0.750</td><td align="center">0.765</td><td align="center">0.000</td><td align="center">0.724</td><td align="center">0.747</td><td align="center">0.000</td></tr><tr><td align="center">0.2</td><td align="center">0.2</td><td align="center">0.08</td><td align="center">3.8</td><td align="center">0.043</td><td align="center">0.029</td><td align="center">0.581</td><td align="center">0.584</td><td align="center">0.000</td><td align="center">0.544</td><td align="center">0.564</td><td align="center">0.000</td></tr><tr><td align="center">0.2</td><td align="center">0.4</td><td align="center">0.05</td><td align="center">4.0</td><td align="center">0.048</td><td align="center">0.045</td><td align="center">0.885</td><td align="center">0.919</td><td align="center">0.021</td><td align="center">0.896</td><td align="center">0.922</td><td align="center">0.041</td></tr><tr><td align="center">0.2</td><td align="center">0.4</td><td align="center">0.08</td><td align="center">3.9</td><td align="center">0.049</td><td align="center">0.038</td><td align="center">0.806</td><td align="center">0.827</td><td align="center">0.025</td><td align="center">0.791</td><td align="center">0.820</td><td align="center">0.053</td></tr><tr><td align="center">0.2</td><td align="center">0.6</td><td align="center">0.05</td><td align="center">4.0</td><td align="center">0.054</td><td align="center">0.054</td><td align="center">0.933</td><td align="center">0.954</td><td align="center">0.130</td><td align="center">0.938</td><td align="center">0.957</td><td align="center">0.178</td></tr><tr><td align="center">0.2</td><td align="center">0.6</td><td align="center">0.08</td><td align="center">4.0</td><td align="center">0.051</td><td align="center">0.045</td><td align="center">0.875</td><td align="center">0.903</td><td align="center">0.142</td><td align="center">0.870</td><td align="center">0.888</td><td align="center">0.186</td></tr><tr><td align="center">0.2</td><td align="center">1.0</td><td align="center">0.05</td><td align="center">3.9</td><td align="center">0.053</td><td align="center">0.058</td><td align="center">0.968</td><td align="center">0.976</td><td align="center">0.355</td><td align="center">0.962</td><td align="center">0.978</td><td align="center">0.429</td></tr><tr><td align="center">0.2</td><td align="center">1.0</td><td align="center">0.08</td><td align="center">4.0</td><td align="center">0.045</td><td align="center">0.055</td><td align="center">0.924</td><td align="center">0.954</td><td align="center">0.382</td><td align="center">0.926</td><td align="center">0.954</td><td align="center">0.444</td></tr><tr><td align="center">0.2</td><td align="center">1.5</td><td align="center">0.05</td><td align="center">4.0</td><td align="center">0.047</td><td align="center">0.061</td><td align="center">0.974</td><td align="center">0.983</td><td align="center">0.526</td><td align="center">0.974</td><td align="center">0.986</td><td align="center">0.597</td></tr><tr><td align="center">0.2</td><td align="center">1.5</td><td align="center">0.08</td><td align="center">4.0</td><td align="center">0.052</td><td align="center">0.059</td><td align="center">0.952</td><td align="center">0.971</td><td align="center">0.580</td><td align="center">0.960</td><td align="center">0.974</td><td align="center">0.664</td></tr><tr><td align="center">0.2</td><td align="center">4.0</td><td align="center">0.05</td><td align="center">3.9</td><td align="center">0.053</td><td align="center">0.065</td><td align="center">0.991</td><td align="center">0.993</td><td align="center">0.826</td><td align="center">0.993</td><td align="center">0.997</td><td align="center">0.860</td></tr><tr><td align="center">0.2</td><td align="center">4.0</td><td align="center">0.08</td><td align="center">4.0</td><td align="center">0.049</td><td align="center">0.063</td><td align="center">0.982</td><td align="center">0.990</td><td align="center">0.834</td><td align="center">0.989</td><td align="center">0.993</td><td align="center">0.868</td></tr><tr><td colspan="12"><hr/></td></tr><tr><td align="center">0.4</td><td align="center">0.2</td><td align="center">0.05</td><td align="center">4.0</td><td align="center">0.044</td><td align="center">0.035</td><td align="center">0.721</td><td align="center">0.754</td><td align="center">0.000</td><td align="center">0.726</td><td align="center">0.752</td><td align="center">0.000</td></tr><tr><td align="center">0.4</td><td align="center">0.2</td><td align="center">0.08</td><td align="center">3.9</td><td align="center">0.045</td><td align="center">0.030</td><td align="center">0.568</td><td align="center">0.544</td><td align="center">0.000</td><td align="center">0.528</td><td align="center">0.550</td><td align="center">0.001</td></tr><tr><td align="center">0.4</td><td align="center">0.4</td><td align="center">0.05</td><td align="center">3.9</td><td align="center">0.048</td><td align="center">0.044</td><td align="center">0.898</td><td align="center">0.912</td><td align="center">0.023</td><td align="center">0.884</td><td align="center">0.919</td><td align="center">0.040</td></tr><tr><td align="center">0.4</td><td align="center">0.4</td><td align="center">0.08</td><td align="center">4.2</td><td align="center">0.050</td><td align="center">0.040</td><td align="center">0.810</td><td align="center">0.823</td><td align="center">0.041</td><td align="center">0.799</td><td align="center">0.830</td><td align="center">0.062</td></tr><tr><td align="center">0.4</td><td align="center">0.6</td><td align="center">0.05</td><td align="center">4.0</td><td align="center">0.047</td><td align="center">0.050</td><td align="center">0.929</td><td align="center">0.945</td><td align="center">0.142</td><td align="center">0.915</td><td align="center">0.940</td><td align="center">0.176</td></tr><tr><td align="center">0.4</td><td align="center">0.6</td><td align="center">0.08</td><td align="center">4.0</td><td align="center">0.046</td><td align="center">0.046</td><td align="center">0.869</td><td align="center">0.895</td><td align="center">0.152</td><td align="center">0.859</td><td align="center">0.894</td><td align="center">0.175</td></tr><tr><td align="center">0.4</td><td align="center">1.0</td><td align="center">0.05</td><td align="center">4.0</td><td align="center">0.051</td><td align="center">0.059</td><td align="center">0.968</td><td align="center">0.983</td><td align="center">0.362</td><td align="center">0.961</td><td align="center">0.976</td><td align="center">0.428</td></tr><tr><td align="center">0.4</td><td align="center">1.0</td><td align="center">0.08</td><td align="center">3.9</td><td align="center">0.050</td><td align="center">0.052</td><td align="center">0.937</td><td align="center">0.956</td><td align="center">0.380</td><td align="center">0.937</td><td align="center">0.951</td><td align="center">0.470</td></tr><tr><td align="center">0.4</td><td align="center">1.5</td><td align="center">0.05</td><td align="center">4.0</td><td align="center">0.049</td><td align="center">0.062</td><td align="center">0.976</td><td align="center">0.981</td><td align="center">0.553</td><td align="center">0.967</td><td align="center">0.978</td><td align="center">0.599</td></tr><tr><td align="center">0.4</td><td align="center">1.5</td><td align="center">0.08</td><td align="center">4.1</td><td align="center">0.050</td><td align="center">0.056</td><td align="center">0.951</td><td align="center">0.969</td><td align="center">0.574</td><td align="center">0.955</td><td align="center">0.971</td><td align="center">0.620</td></tr><tr><td align="center">0.4</td><td align="center">4.0</td><td align="center">0.05</td><td align="center">3.9</td><td align="center">0.047</td><td align="center">0.065</td><td align="center">0.991</td><td align="center">0.995</td><td align="center">0.835</td><td align="center">0.991</td><td align="center">0.997</td><td align="center">0.848</td></tr><tr><td align="center">0.4</td><td align="center">4.0</td><td align="center">0.08</td><td align="center">4.0</td><td align="center">0.050</td><td align="center">0.064</td><td align="center">0.983</td><td align="center">0.988</td><td align="center">0.832</td><td align="center">0.990</td><td align="center">0.997</td><td align="center">0.863</td></tr><tr><td colspan="12"><hr/></td></tr><tr><td align="center">0.6</td><td align="center">0.2</td><td align="center">0.05</td><td align="center">4.0</td><td align="center">0.045</td><td align="center">0.037</td><td align="center">0.720</td><td align="center">0.748</td><td align="center">0.001</td><td align="center">0.740</td><td align="center">0.775</td><td align="center">0.004</td></tr><tr><td align="center">0.6</td><td align="center">0.2</td><td align="center">0.08</td><td align="center">4.0</td><td align="center">0.044</td><td align="center">0.029</td><td align="center">0.565</td><td align="center">0.584</td><td align="center">0.001</td><td align="center">0.533</td><td align="center">0.572</td><td align="center">0.003</td></tr><tr><td align="center">0.6</td><td align="center">0.4</td><td align="center">0.05</td><td align="center">4.0</td><td align="center">0.051</td><td align="center">0.047</td><td align="center">0.880</td><td align="center">0.896</td><td align="center">0.034</td><td align="center">0.892</td><td align="center">0.915</td><td align="center">0.055</td></tr><tr><td align="center">0.6</td><td align="center">0.4</td><td align="center">0.08</td><td align="center">4.0</td><td align="center">0.047</td><td align="center">0.041</td><td align="center">0.800</td><td align="center">0.824</td><td align="center">0.060</td><td align="center">0.792</td><td align="center">0.818</td><td align="center">0.079</td></tr><tr><td align="center">0.6</td><td align="center">0.6</td><td align="center">0.05</td><td align="center">4.1</td><td align="center">0.051</td><td align="center">0.054</td><td align="center">0.919</td><td align="center">0.938</td><td align="center">0.139</td><td align="center">0.926</td><td align="center">0.947</td><td align="center">0.186</td></tr><tr><td align="center">0.6</td><td align="center">0.6</td><td align="center">0.08</td><td align="center">3.9</td><td align="center">0.047</td><td align="center">0.046</td><td align="center">0.877</td><td align="center">0.893</td><td align="center">0.179</td><td align="center">0.889</td><td align="center">0.924</td><td align="center">0.230</td></tr><tr><td align="center">0.6</td><td align="center">1.0</td><td align="center">0.05</td><td align="center">4.0</td><td align="center">0.050</td><td align="center">0.060</td><td align="center">0.959</td><td align="center">0.976</td><td align="center">0.370</td><td align="center">0.956</td><td align="center">0.971</td><td align="center">0.435</td></tr><tr><td align="center">0.6</td><td align="center">1.0</td><td align="center">0.08</td><td align="center">4.0</td><td align="center">0.048</td><td align="center">0.054</td><td align="center">0.930</td><td align="center">0.947</td><td align="center">0.445</td><td align="center">0.922</td><td align="center">0.948</td><td align="center">0.452</td></tr><tr><td align="center">0.6</td><td align="center">1.5</td><td align="center">0.05</td><td align="center">4.0</td><td align="center">0.049</td><td align="center">0.063</td><td align="center">0.975</td><td align="center">0.983</td><td align="center">0.549</td><td align="center">0.977</td><td align="center">0.985</td><td align="center">0.601</td></tr><tr><td align="center">0.6</td><td align="center">1.5</td><td align="center">0.08</td><td align="center">3.9</td><td align="center">0.049</td><td align="center">0.059</td><td align="center">0.962</td><td align="center">0.973</td><td align="center">0.588</td><td align="center">0.960</td><td align="center">0.965</td><td align="center">0.652</td></tr><tr><td align="center">0.6</td><td align="center">4.0</td><td align="center">0.05</td><td align="center">3.9</td><td align="center">0.053</td><td align="center">0.065</td><td align="center">0.988</td><td align="center">0.993</td><td align="center">0.814</td><td align="center">0.992</td><td align="center">0.994</td><td align="center">0.845</td></tr><tr><td align="center">0.6</td><td align="center">4.0</td><td align="center">0.08</td><td align="center">3.9</td><td align="center">0.055</td><td align="center">0.066</td><td align="center">0.979</td><td align="center">0.987</td><td align="center">0.835</td><td align="center">0.984</td><td align="center">0.995</td><td align="center">0.859</td></tr></tbody></table></table-wrap> <table-wrap position="float" id="T3"><label>Table 3</label><caption><p>MLPA probemix composition</p></caption><table frame="hsides" rules="groups"><thead><tr><td align="left">Gene/Region</td><td align="left">Genomic location</td><td align="left">Probe size</td><td align="left">Band</td><td align="left">Comments</td></tr></thead><tbody><tr><td align="left">ENm323</td><td align="left">chr6:108,723,531–108,723,590</td><td align="left">126</td><td align="left">6q21</td><td align="left">single copy number region</td></tr><tr><td align="left">ENm013</td><td align="left">chr7:90,250,124–90,250,183</td><td align="left">105</td><td align="left">7q21.13</td><td align="left">single copy number region</td></tr><tr><td align="left">ENm014</td><td align="left">chr7:126,866,339–126,866,393</td><td align="left">99</td><td align="left">7q31.33</td><td align="left">single copy number region</td></tr><tr><td align="left">RNAseP (RPP30)</td><td align="left">chr10:92,621,710–92,621,757</td><td align="left">90</td><td align="left">10q23.31</td><td align="left">single copy number region</td></tr><tr><td align="left">ENr222</td><td align="left">chr10:92,621,710–92,621,757</td><td align="left">147</td><td align="left">6q23.2</td><td align="left">single copy number region</td></tr><tr><td align="left">ENr111</td><td align="left">chr13:29,519,123–29,519,182</td><td align="left">123</td><td align="left">13q12.3</td><td align="left">single copy number region</td></tr><tr><td align="left">ENr233</td><td align="left">chr15:41,662,068–41,662,127</td><td align="left">141</td><td align="left">15q15.3</td><td align="left">single copy number region</td></tr><tr><td align="left">ENm313</td><td align="left">chr16:61,141,268–61,141,327</td><td align="left">114</td><td align="left">16q21</td><td align="left">single copy number region</td></tr><tr><td align="left">ENr213</td><td align="left">chr18:24,170,726–24,170,785</td><td align="left">132</td><td align="left">18q12.1</td><td align="left">single copy number region</td></tr><tr><td align="left">RP11–71N21</td><td align="left">chr10:51,942,608–51,942,682</td><td align="left">144</td><td align="left">10q11.23</td><td align="left">10q duplication</td></tr><tr><td align="left">ZWINT</td><td align="left">chr10:57,789,521–57,789,580</td><td align="left">120</td><td align="left">10q21.1</td><td align="left">10q duplication</td></tr><tr><td align="left">PHYLIP</td><td align="left">chr10:60,674,910–60,674,995</td><td align="left">138</td><td align="left">10q21.1</td><td align="left">10q duplication</td></tr><tr><td align="left">SNRPN</td><td align="left">chr15:22,764,255–22,764,313</td><td align="left">111</td><td align="left">15q11.2</td><td align="left">15q11 deletion</td></tr><tr><td align="left">UBEA3A</td><td align="left">chr15:23,133,602–23,133,661</td><td align="left">129</td><td align="left">15q11.2</td><td align="left">15q11 deletion</td></tr><tr><td align="left">UBE3A</td><td align="left">chr15:23,171,946–23,171,999</td><td align="left">96</td><td align="left">15q11.2</td><td align="left">15q11 deletion</td></tr><tr><td align="left">HIRA</td><td align="left">chr22:17,698,971–17,699,021</td><td align="left">93</td><td align="left">22q11.21</td><td align="left">22q11 deletion</td></tr></tbody></table></table-wrap> <table-wrap position="float" id="T4"><label>Table 4</label><caption><p>MLPA results from the validation study.</p></caption><table frame="hsides" rules="groups"><tbody><tr><td/><td align="center" colspan="6">Prader-Willi</td><td align="center" colspan="6">DiGeorge</td><td align="center" colspan="6">Autism</td><td align="center" colspan="9">HapMap</td></tr><tr><td/><td colspan="6"><hr/></td><td colspan="6"><hr/></td><td colspan="6"><hr/></td><td colspan="9"><hr/></td></tr><tr><td align="left">Gene/Region</td><td align="center" colspan="3">#1</td><td align="center" colspan="3">#2</td><td align="center" colspan="3">#1</td><td align="center" colspan="3">#2</td><td align="center" colspan="3">#1</td><td align="center" colspan="3">#2</td><td align="center" colspan="3">#1</td><td align="center" colspan="3">#2</td><td align="center" colspan="3">#3</td></tr><tr><td colspan="1"><hr/></td><td colspan="3"><hr/></td><td colspan="3"><hr/></td><td colspan="3"><hr/></td><td colspan="3"><hr/></td><td colspan="3"><hr/></td><td colspan="3"><hr/></td><td colspan="3"><hr/></td><td colspan="3"><hr/></td><td colspan="3"><hr/></td></tr><tr><td align="left">ENm323</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td></tr><tr><td align="left">ENm013</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td></tr><tr><td align="left">ENm014</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td></tr><tr><td align="left">RNAseP (RPP30)</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td></tr><tr><td align="left">ENr222</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td></tr><tr><td align="left">ENr111</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td></tr><tr><td align="left">ENr233</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td></tr><tr><td align="left">ENm313</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td></tr><tr><td align="left">ENr213</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td></tr><tr><td align="left">RP11-71N21</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">1</td><td align="center">1</td><td align="center">1</td><td align="center">1</td><td align="center">1</td><td align="center">1</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td></tr><tr><td align="left">ZWINT</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">1</td><td align="center">1</td><td align="center">1</td><td align="center"><bold>1</bold></td><td align="center"><bold>0</bold></td><td align="center"><bold>1</bold></td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td></tr><tr><td align="left">PHYLIP</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">1</td><td align="center">1</td><td align="center">1</td><td align="center">1</td><td align="center">1</td><td align="center">1</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td></tr><tr><td align="left">SNRPN</td><td align="center">-1</td><td align="center">-1</td><td align="center">-1</td><td align="center">-1</td><td align="center">-1</td><td align="center">-1</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td></tr><tr><td align="left">UBEA3A</td><td align="center">-1</td><td align="center">-1</td><td align="center">-1</td><td align="center">-1</td><td align="center">-1</td><td align="center">-1</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td></tr><tr><td align="left">UBE3A</td><td align="center">-1</td><td align="center">-1</td><td align="center">-1</td><td align="center">-1</td><td align="center">-1</td><td align="center">-1</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center"><bold>0</bold></td><td align="center"><bold>0</bold></td><td align="center"><bold>-1</bold></td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td></tr><tr><td align="left">HIRA</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">-1</td><td align="center">-1</td><td align="center">-1</td><td align="center">-1</td><td align="center">-1</td><td align="center">-1</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td></tr></tbody></table></table-wrap>
<inline-formula><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M1" name="1471-2105-9-261-i1" overflow="scroll"><mml:semantics><mml:mrow><mml:msub><mml:mover accent="true"><mml:mi>x</mml:mi><mml:mo>¯</mml:mo></mml:mover><mml:mrow><mml:mi>a</mml:mi><mml:mi>l</mml:mi><mml:mi>t</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:semantics></mml:math></inline-formula> <inline-formula><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M2" name="1471-2105-9-261-i1" overflow="scroll"><mml:semantics><mml:mrow><mml:msub><mml:mover accent="true"><mml:mi>x</mml:mi><mml:mo>¯</mml:mo></mml:mover><mml:mrow><mml:mi>a</mml:mi><mml:mi>l</mml:mi><mml:mi>t</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:semantics></mml:math></inline-formula> <inline-formula><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M3" name="1471-2105-9-261-i1" overflow="scroll"><mml:semantics><mml:mrow><mml:msub><mml:mover accent="true"><mml:mi>x</mml:mi><mml:mo>¯</mml:mo></mml:mover><mml:mrow><mml:mi>a</mml:mi><mml:mi>l</mml:mi><mml:mi>t</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:semantics></mml:math></inline-formula> <inline-formula><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M4" name="1471-2105-9-261-i1" overflow="scroll"><mml:semantics><mml:mrow><mml:msub><mml:mover accent="true"><mml:mi>x</mml:mi><mml:mo>¯</mml:mo></mml:mover><mml:mrow><mml:mi>a</mml:mi><mml:mi>l</mml:mi><mml:mi>t</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:semantics></mml:math></inline-formula> <inline-formula><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M5" name="1471-2105-9-261-i2" overflow="scroll"><mml:semantics><mml:mrow><mml:mfrac><mml:mrow><mml:mn>0.1423</mml:mn></mml:mrow><mml:mrow><mml:mn>0.0461</mml:mn></mml:mrow></mml:mfrac></mml:mrow></mml:semantics></mml:math></inline-formula> <inline-formula><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M6" name="1471-2105-9-261-i3" overflow="scroll"><mml:semantics><mml:mrow><mml:msub><mml:mover accent="true"><mml:mi>H</mml:mi><mml:mo>˜</mml:mo></mml:mover><mml:mrow><mml:mi>i</mml:mi><mml:mi>p</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mstyle displaystyle="true"><mml:msubsup><mml:mo>∑</mml:mo><mml:mrow><mml:mi>k</mml:mi><mml:mo>=</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mi>K</mml:mi></mml:msubsup><mml:mrow><mml:msub><mml:mi>H</mml:mi><mml:mrow><mml:mi>i</mml:mi><mml:mi>p</mml:mi><mml:mi>k</mml:mi></mml:mrow></mml:msub><mml:mo>/</mml:mo><mml:mi>K</mml:mi></mml:mrow></mml:mstyle></mml:mrow></mml:semantics></mml:math></inline-formula> <inline-formula><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M7" name="1471-2105-9-261-i4" overflow="scroll"><mml:semantics><mml:mrow><mml:msub><mml:mi>S</mml:mi><mml:mi>i</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:mstyle displaystyle="true"><mml:msubsup><mml:mo>∑</mml:mo><mml:mrow><mml:mi>p</mml:mi><mml:mo>=</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mi>P</mml:mi></mml:msubsup><mml:mrow><mml:msub><mml:mover accent="true"><mml:mi>H</mml:mi><mml:mo>˜</mml:mo></mml:mover><mml:mrow><mml:mi>i</mml:mi><mml:mi>p</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:mstyle></mml:mrow></mml:semantics></mml:math></inline-formula> <inline-formula><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M8" name="1471-2105-9-261-i5" overflow="scroll"><mml:semantics><mml:mrow><mml:msub><mml:mover accent="true"><mml:mi>Y</mml:mi><mml:mo>˜</mml:mo></mml:mover><mml:mrow><mml:mi>i</mml:mi><mml:mi>p</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:msub><mml:mover accent="true"><mml:mi>H</mml:mi><mml:mo>˜</mml:mo></mml:mover><mml:mrow><mml:mi>i</mml:mi><mml:mi>p</mml:mi></mml:mrow></mml:msub><mml:mo>/</mml:mo><mml:msub><mml:mi>S</mml:mi><mml:mi>i</mml:mi></mml:msub></mml:mrow></mml:semantics></mml:math></inline-formula> <inline-formula><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M9" name="1471-2105-9-261-i6" overflow="scroll"><mml:semantics><mml:mrow><mml:msub><mml:mover accent="true"><mml:mi>H</mml:mi><mml:mo>˜</mml:mo></mml:mover><mml:mrow><mml:msub><mml:mi>p</mml:mi><mml:mi>c</mml:mi></mml:msub></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mi>α</mml:mi><mml:mo>+</mml:mo><mml:mi>β</mml:mi><mml:msub><mml:mi>X</mml:mi><mml:mrow><mml:msub><mml:mi>p</mml:mi><mml:mi>c</mml:mi></mml:msub></mml:mrow></mml:msub></mml:mrow></mml:semantics></mml:math></inline-formula> <inline-formula><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M10" name="1471-2105-9-261-i7" overflow="scroll"><mml:semantics><mml:mrow><mml:msub><mml:mover accent="true"><mml:mi>H</mml:mi><mml:mo>˜</mml:mo></mml:mover><mml:mrow><mml:msub><mml:mi>p</mml:mi><mml:mi>c</mml:mi></mml:msub></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mi>α</mml:mi><mml:mo>+</mml:mo><mml:msub><mml:mi>β</mml:mi><mml:mn>1</mml:mn></mml:msub><mml:msub><mml:mi>X</mml:mi><mml:mrow><mml:msub><mml:mi>p</mml:mi><mml:mi>c</mml:mi></mml:msub></mml:mrow></mml:msub><mml:mo>+</mml:mo><mml:msub><mml:mi>β</mml:mi><mml:mn>2</mml:mn></mml:msub><mml:msubsup><mml:mi>X</mml:mi><mml:mrow><mml:msub><mml:mi>p</mml:mi><mml:mi>c</mml:mi></mml:msub></mml:mrow><mml:mn>2</mml:mn></mml:msubsup></mml:mrow></mml:semantics></mml:math></inline-formula> <inline-formula><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M11" name="1471-2105-9-261-i8" overflow="scroll"><mml:semantics><mml:mrow><mml:msub><mml:mover accent="true"><mml:mi>H</mml:mi><mml:mo>˜</mml:mo></mml:mover><mml:mi>p</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:mi>α</mml:mi><mml:msup><mml:mtext>e</mml:mtext><mml:mrow><mml:mo>−</mml:mo><mml:mi>β</mml:mi><mml:msub><mml:mi>X</mml:mi><mml:mi>p</mml:mi></mml:msub></mml:mrow></mml:msup></mml:mrow></mml:semantics></mml:math></inline-formula> <inline-formula><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M12" name="1471-2105-9-261-i9" overflow="scroll"><mml:semantics><mml:mrow><mml:msub><mml:mover accent="true"><mml:mi>H</mml:mi><mml:mo>˜</mml:mo></mml:mover><mml:mrow><mml:mi>i</mml:mi><mml:mi>p</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:semantics></mml:math></inline-formula> <disp-formula id="bmcM1"><label>(1)</label><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M13" name="1471-2105-9-261-i10" overflow="scroll"><mml:semantics><mml:mrow><mml:mtable><mml:mtr><mml:mtd><mml:mrow><mml:msub><mml:mi>H</mml:mi><mml:mrow><mml:mi>i</mml:mi><mml:mi>p</mml:mi><mml:mi>k</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mo stretchy="false">(</mml:mo><mml:msub><mml:mi>φ</mml:mi><mml:mrow><mml:mn>1</mml:mn><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mo>−</mml:mo><mml:msub><mml:mi>φ</mml:mi><mml:mrow><mml:mn>2</mml:mn><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mo stretchy="false">)</mml:mo><mml:mi>exp</mml:mi><mml:mo>⁡</mml:mo><mml:mo stretchy="false">[</mml:mo><mml:mo>−</mml:mo><mml:mfrac><mml:mn>1</mml:mn><mml:mrow><mml:msub><mml:mi>φ</mml:mi><mml:mrow><mml:mn>3</mml:mn><mml:mi>i</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:mfrac><mml:mo stretchy="false">(</mml:mo><mml:msub><mml:mi>X</mml:mi><mml:mi>p</mml:mi></mml:msub><mml:mo>−</mml:mo><mml:msub><mml:mrow><mml:mtext>off</mml:mtext></mml:mrow><mml:mi>i</mml:mi></mml:msub><mml:mo stretchy="false">)</mml:mo><mml:mo stretchy="false">]</mml:mo><mml:mo>+</mml:mo><mml:msub><mml:mi>φ</mml:mi><mml:mrow><mml:mn>2</mml:mn><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mo>+</mml:mo><mml:msub><mml:mi>ε</mml:mi><mml:mrow><mml:mi>i</mml:mi><mml:mi>p</mml:mi><mml:mi>k</mml:mi></mml:mrow></mml:msub><mml:mo>,</mml:mo></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mrow><mml:msub><mml:mi>φ</mml:mi><mml:mi>i</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:mrow><mml:mo>[</mml:mo><mml:mrow><mml:mtable><mml:mtr><mml:mtd><mml:mrow><mml:msub><mml:mi>φ</mml:mi><mml:mrow><mml:mn>1</mml:mn><mml:mi>i</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mrow><mml:msub><mml:mi>φ</mml:mi><mml:mrow><mml:mn>2</mml:mn><mml:mi>i</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mrow><mml:msub><mml:mi>φ</mml:mi><mml:mrow><mml:mn>3</mml:mn><mml:mi>i</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:mrow><mml:mo>]</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:mrow><mml:mo>[</mml:mo><mml:mrow><mml:mtable><mml:mtr><mml:mtd><mml:mrow><mml:msub><mml:mi>β</mml:mi><mml:mn>1</mml:mn></mml:msub></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mrow><mml:msub><mml:mi>β</mml:mi><mml:mn>2</mml:mn></mml:msub></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mrow><mml:msub><mml:mi>β</mml:mi><mml:mn>3</mml:mn></mml:msub></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:mrow><mml:mo>]</mml:mo></mml:mrow><mml:mo>+</mml:mo><mml:mrow><mml:mo>[</mml:mo><mml:mrow><mml:mtable><mml:mtr><mml:mtd><mml:mrow><mml:msub><mml:mi>b</mml:mi><mml:mrow><mml:mn>1</mml:mn><mml:mi>i</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mrow><mml:msub><mml:mi>b</mml:mi><mml:mrow><mml:mn>2</mml:mn><mml:mi>i</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mrow><mml:msub><mml:mi>b</mml:mi><mml:mrow><mml:mn>3</mml:mn><mml:mi>i</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:mrow><mml:mo>]</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:mi>β</mml:mi><mml:mo>+</mml:mo><mml:msub><mml:mstyle mathvariant="bold" mathsize="normal"><mml:mi>b</mml:mi></mml:mstyle><mml:mi>i</mml:mi></mml:msub><mml:mo>,</mml:mo></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mrow><mml:mtable><mml:mtr><mml:mtd><mml:mrow><mml:msub><mml:mstyle mathvariant="bold" mathsize="normal"><mml:mi>b</mml:mi></mml:mstyle><mml:mi>i</mml:mi></mml:msub><mml:mo>~</mml:mo><mml:mi mathvariant="script">N</mml:mi><mml:mo stretchy="false">(</mml:mo><mml:mstyle mathvariant="bold" mathsize="normal"><mml:mn>0</mml:mn></mml:mstyle><mml:mo>,</mml:mo><mml:mi>Ψ</mml:mi><mml:mo stretchy="false">)</mml:mo><mml:mo>,</mml:mo></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:msub><mml:mi>ε</mml:mi><mml:mrow><mml:mi>i</mml:mi><mml:mi>j</mml:mi><mml:mi>k</mml:mi></mml:mrow></mml:msub><mml:mo>~</mml:mo><mml:mi mathvariant="script">N</mml:mi><mml:mo stretchy="false">(</mml:mo><mml:mn>0</mml:mn><mml:mo>,</mml:mo><mml:msup><mml:mi>σ</mml:mi><mml:mn>2</mml:mn></mml:msup><mml:mo stretchy="false">)</mml:mo><mml:mo>,</mml:mo></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:mrow></mml:semantics></mml:math></disp-formula> <inline-formula><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M14" name="1471-2105-9-261-i11" overflow="scroll"><mml:semantics><mml:mrow><mml:mstyle displaystyle="true"><mml:msubsup><mml:mo>∑</mml:mo><mml:mrow><mml:mi>k</mml:mi><mml:mo>=</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mi>K</mml:mi></mml:msubsup><mml:mrow><mml:msub><mml:mi>H</mml:mi><mml:mrow><mml:mi>i</mml:mi><mml:mn>1</mml:mn><mml:mi>k</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:mstyle></mml:mrow></mml:semantics></mml:math></inline-formula> <inline-formula><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M15" name="1471-2105-9-261-i12" overflow="scroll"><mml:semantics><mml:mrow><mml:msub><mml:mover accent="true"><mml:mi>H</mml:mi><mml:mo>^</mml:mo></mml:mover><mml:mrow><mml:mi>i</mml:mi><mml:mi>p</mml:mi><mml:mi>k</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:semantics></mml:math></inline-formula> <inline-formula><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M16" name="1471-2105-9-261-i12" overflow="scroll"><mml:semantics><mml:mrow><mml:msub><mml:mover accent="true"><mml:mi>H</mml:mi><mml:mo>^</mml:mo></mml:mover><mml:mrow><mml:mi>i</mml:mi><mml:mi>p</mml:mi><mml:mi>k</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:semantics></mml:math></inline-formula> <disp-formula id="bmcM2"><label>(2)</label><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M17" name="1471-2105-9-261-i13" overflow="scroll"><mml:semantics><mml:mrow><mml:mtable><mml:mtr><mml:mtd><mml:mrow><mml:mi>log</mml:mi><mml:mo>⁡</mml:mo><mml:mo stretchy="false">(</mml:mo><mml:msub><mml:mi>Y</mml:mi><mml:mrow><mml:mi>t</mml:mi><mml:mi>p</mml:mi><mml:mi>k</mml:mi></mml:mrow></mml:msub><mml:mo stretchy="false">)</mml:mo><mml:mo>=</mml:mo><mml:mi>μ</mml:mi><mml:mo>+</mml:mo><mml:msub><mml:mi>α</mml:mi><mml:mi>t</mml:mi></mml:msub><mml:mo>+</mml:mo><mml:msub><mml:mi>β</mml:mi><mml:mi>p</mml:mi></mml:msub><mml:mo>+</mml:mo><mml:msub><mml:mi>γ</mml:mi><mml:mrow><mml:mi>t</mml:mi><mml:mi>p</mml:mi></mml:mrow></mml:msub><mml:mo>+</mml:mo><mml:msub><mml:mi>ε</mml:mi><mml:mrow><mml:mi>t</mml:mi><mml:mi>p</mml:mi><mml:mi>k</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mrow><mml:mtable><mml:mtr><mml:mtd><mml:mrow><mml:msub><mml:mi>β</mml:mi><mml:mi>p</mml:mi></mml:msub><mml:mo>~</mml:mo><mml:mi mathvariant="script">N</mml:mi><mml:mo stretchy="false">(</mml:mo><mml:mn>0</mml:mn><mml:mo>,</mml:mo><mml:msubsup><mml:mi>σ</mml:mi><mml:mi>β</mml:mi><mml:mn>2</mml:mn></mml:msubsup><mml:mo stretchy="false">)</mml:mo><mml:mo>,</mml:mo></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:mi>p</mml:mi><mml:mo>=</mml:mo><mml:mn>1</mml:mn><mml:mo>,</mml:mo><mml:mn>...</mml:mn><mml:mo>,</mml:mo><mml:mi>P</mml:mi></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mrow><mml:mtable><mml:mtr><mml:mtd><mml:mrow><mml:msub><mml:mi>γ</mml:mi><mml:mrow><mml:mi>t</mml:mi><mml:mi>p</mml:mi></mml:mrow></mml:msub><mml:mo>~</mml:mo><mml:mi mathvariant="script">N</mml:mi><mml:mo stretchy="false">(</mml:mo><mml:mn>0</mml:mn><mml:mo>,</mml:mo><mml:msubsup><mml:mi>σ</mml:mi><mml:mi>γ</mml:mi><mml:mn>2</mml:mn></mml:msubsup><mml:mo stretchy="false">)</mml:mo><mml:mo>,</mml:mo></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:mi>t</mml:mi><mml:mo>=</mml:mo><mml:mn>0</mml:mn><mml:mo>,</mml:mo><mml:mn>1</mml:mn><mml:mo>,</mml:mo></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:mi>p</mml:mi><mml:mo>=</mml:mo><mml:mn>1</mml:mn><mml:mo>,</mml:mo><mml:mn>...</mml:mn><mml:mo>,</mml:mo><mml:mi>P</mml:mi></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mrow><mml:mtable><mml:mtr><mml:mtd><mml:mrow><mml:msub><mml:mi>ε</mml:mi><mml:mrow><mml:mi>t</mml:mi><mml:mi>p</mml:mi><mml:mi>k</mml:mi></mml:mrow></mml:msub><mml:mo>~</mml:mo><mml:mi mathvariant="script">N</mml:mi><mml:mo stretchy="false">(</mml:mo><mml:mn>0</mml:mn><mml:mo>,</mml:mo><mml:msubsup><mml:mi>σ</mml:mi><mml:mi>ε</mml:mi><mml:mn>2</mml:mn></mml:msubsup><mml:mo stretchy="false">)</mml:mo><mml:mo>,</mml:mo></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:mi>t</mml:mi><mml:mo>=</mml:mo><mml:mn>0</mml:mn><mml:mo>,</mml:mo><mml:mn>1</mml:mn><mml:mo>,</mml:mo></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:mtable><mml:mtr><mml:mtd><mml:mrow><mml:mi>p</mml:mi><mml:mo>=</mml:mo><mml:mn>1</mml:mn><mml:mo>,</mml:mo><mml:mn>...</mml:mn><mml:mo>,</mml:mo><mml:mi>P</mml:mi><mml:mo>,</mml:mo></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:mi>k</mml:mi><mml:mo>=</mml:mo><mml:mn>1</mml:mn><mml:mo>,</mml:mo><mml:mn>...</mml:mn><mml:mo>,</mml:mo><mml:mi>K</mml:mi></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:mrow></mml:semantics></mml:math></disp-formula> <inline-formula><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M18" name="1471-2105-9-261-i14" overflow="scroll"><mml:semantics><mml:mrow><mml:msubsup><mml:mi>σ</mml:mi><mml:mi>β</mml:mi><mml:mn>2</mml:mn></mml:msubsup></mml:mrow></mml:semantics></mml:math></inline-formula> <inline-formula><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M19" name="1471-2105-9-261-i15" overflow="scroll"><mml:semantics><mml:mrow><mml:msubsup><mml:mi>σ</mml:mi><mml:mi>γ</mml:mi><mml:mn>2</mml:mn></mml:msubsup></mml:mrow></mml:semantics></mml:math></inline-formula> <inline-formula><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M20" name="1471-2105-9-261-i16" overflow="scroll"><mml:semantics><mml:mrow><mml:msubsup><mml:mi>σ</mml:mi><mml:mi>ε</mml:mi><mml:mn>2</mml:mn></mml:msubsup></mml:mrow></mml:semantics></mml:math></inline-formula> <inline-formula><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M21" name="1471-2105-9-261-i16" overflow="scroll"><mml:semantics><mml:mrow><mml:msubsup><mml:mi>σ</mml:mi><mml:mi>ε</mml:mi><mml:mn>2</mml:mn></mml:msubsup></mml:mrow></mml:semantics></mml:math></inline-formula> <disp-formula id="bmcM3"><label>(3)</label><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M22" name="1471-2105-9-261-i17" overflow="scroll"><mml:semantics><mml:mrow><mml:mi>log</mml:mi><mml:mo>⁡</mml:mo><mml:mo stretchy="false">(</mml:mo><mml:msub><mml:mi>Y</mml:mi><mml:mrow><mml:mn>0</mml:mn><mml:mi>p</mml:mi><mml:mi>l</mml:mi></mml:mrow></mml:msub><mml:mo stretchy="false">)</mml:mo><mml:mo>−</mml:mo><mml:mi>log</mml:mi><mml:mo>⁡</mml:mo><mml:mo stretchy="false">(</mml:mo><mml:msub><mml:mi>Y</mml:mi><mml:mrow><mml:mn>0</mml:mn><mml:mi>p</mml:mi><mml:mi>k</mml:mi></mml:mrow></mml:msub><mml:mo stretchy="false">)</mml:mo><mml:mo>~</mml:mo><mml:mi mathvariant="script">N</mml:mi><mml:mo stretchy="false">(</mml:mo><mml:mn>0</mml:mn><mml:mo>,</mml:mo><mml:mn>2</mml:mn><mml:msubsup><mml:mi>σ</mml:mi><mml:mi>ε</mml:mi><mml:mn>2</mml:mn></mml:msubsup><mml:mo stretchy="false">)</mml:mo><mml:mo>.</mml:mo></mml:mrow></mml:semantics></mml:math></disp-formula> <disp-formula id="bmcM4"><label>(4)</label><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M23" name="1471-2105-9-261-i18" overflow="scroll"><mml:semantics><mml:mrow><mml:mi>log</mml:mi><mml:mo>⁡</mml:mo><mml:mo stretchy="false">(</mml:mo><mml:msub><mml:mi>Y</mml:mi><mml:mrow><mml:mn>0</mml:mn><mml:mi>p</mml:mi><mml:mi>k</mml:mi></mml:mrow></mml:msub><mml:mo stretchy="false">)</mml:mo><mml:mo>−</mml:mo><mml:mi>log</mml:mi><mml:mo>⁡</mml:mo><mml:mo stretchy="false">(</mml:mo><mml:msub><mml:mi>Y</mml:mi><mml:mrow><mml:msup><mml:mi>i</mml:mi><mml:mo>′</mml:mo></mml:msup><mml:mi>p</mml:mi><mml:mi>k</mml:mi></mml:mrow></mml:msub><mml:mo stretchy="false">)</mml:mo><mml:mo>~</mml:mo><mml:mi mathvariant="script">N</mml:mi><mml:mo stretchy="false">(</mml:mo><mml:mi>α</mml:mi><mml:mo>+</mml:mo><mml:msub><mml:mi>γ</mml:mi><mml:mrow><mml:msup><mml:mi>i</mml:mi><mml:mo>′</mml:mo></mml:msup><mml:mi>p</mml:mi></mml:mrow></mml:msub><mml:mo>,</mml:mo><mml:mn>2</mml:mn><mml:msubsup><mml:mi>σ</mml:mi><mml:mi>ε</mml:mi><mml:mn>2</mml:mn></mml:msubsup><mml:mo stretchy="false">)</mml:mo><mml:mo>,</mml:mo></mml:mrow></mml:semantics></mml:math></disp-formula> <inline-formula><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M24" name="1471-2105-9-261-i16" overflow="scroll"><mml:semantics><mml:mrow><mml:msubsup><mml:mi>σ</mml:mi><mml:mi>ε</mml:mi><mml:mn>2</mml:mn></mml:msubsup></mml:mrow></mml:semantics></mml:math></inline-formula> <disp-formula id="bmcM5"><label>(5)</label><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M25" name="1471-2105-9-261-i19" overflow="scroll"><mml:semantics><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mo>−</mml:mo><mml:msub><mml:mi>Z</mml:mi><mml:mrow><mml:mn>1</mml:mn><mml:mo>−</mml:mo><mml:mi>κ</mml:mi><mml:mo>/</mml:mo><mml:mn>2</mml:mn></mml:mrow></mml:msub><mml:msqrt><mml:mrow><mml:mn>2</mml:mn><mml:msubsup><mml:mi>σ</mml:mi><mml:mi>ε</mml:mi><mml:mn>2</mml:mn></mml:msubsup></mml:mrow></mml:msqrt><mml:mo>,</mml:mo><mml:msub><mml:mi>Z</mml:mi><mml:mrow><mml:mn>1</mml:mn><mml:mo>−</mml:mo><mml:mi>κ</mml:mi><mml:mo>/</mml:mo><mml:mn>2</mml:mn></mml:mrow></mml:msub><mml:msqrt><mml:mrow><mml:mn>2</mml:mn><mml:msubsup><mml:mi>σ</mml:mi><mml:mi>ε</mml:mi><mml:mn>2</mml:mn></mml:msubsup></mml:mrow></mml:msqrt><mml:mo stretchy="false">)</mml:mo><mml:mo>,</mml:mo></mml:mrow></mml:semantics></mml:math></disp-formula> <inline-formula><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M26" name="1471-2105-9-261-i20" overflow="scroll"><mml:semantics><mml:mrow><mml:mfrac><mml:mrow><mml:mi>ν</mml:mi><mml:msubsup><mml:mover accent="true"><mml:mi>σ</mml:mi><mml:mo>^</mml:mo></mml:mover><mml:mi>e</mml:mi><mml:mn>2</mml:mn></mml:msubsup></mml:mrow><mml:mrow><mml:msubsup><mml:mi>σ</mml:mi><mml:mi>e</mml:mi><mml:mn>2</mml:mn></mml:msubsup></mml:mrow></mml:mfrac><mml:mo>~</mml:mo><mml:msubsup><mml:mi>χ</mml:mi><mml:mi>ν</mml:mi><mml:mn>2</mml:mn></mml:msubsup></mml:mrow></mml:semantics></mml:math></inline-formula> <inline-formula><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M27" name="1471-2105-9-261-i21" overflow="scroll"><mml:semantics><mml:mrow><mml:msubsup><mml:mover accent="true"><mml:mi>σ</mml:mi><mml:mo>^</mml:mo></mml:mover><mml:mi>e</mml:mi><mml:mn>2</mml:mn></mml:msubsup></mml:mrow></mml:semantics></mml:math></inline-formula> <inline-formula><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M28" name="1471-2105-9-261-i22" overflow="scroll"><mml:semantics><mml:mrow><mml:msubsup><mml:mi>σ</mml:mi><mml:mi>e</mml:mi><mml:mn>2</mml:mn></mml:msubsup></mml:mrow></mml:semantics></mml:math></inline-formula> <disp-formula id="bmcM6"><label>(6)</label><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M29" name="1471-2105-9-261-i23" overflow="scroll"><mml:semantics><mml:mrow><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mo>−</mml:mo><mml:msub><mml:mi>Z</mml:mi><mml:mrow><mml:mn>1</mml:mn><mml:mo>−</mml:mo><mml:mi>κ</mml:mi><mml:mo>/</mml:mo><mml:mn>2</mml:mn></mml:mrow></mml:msub><mml:msqrt><mml:mrow><mml:mfrac><mml:mrow><mml:mi>ν</mml:mi><mml:msubsup><mml:mover accent="true"><mml:mi>σ</mml:mi><mml:mo>^</mml:mo></mml:mover><mml:mi>ε</mml:mi><mml:mn>2</mml:mn></mml:msubsup></mml:mrow><mml:mrow><mml:msubsup><mml:mi>χ</mml:mi><mml:mrow><mml:mi>α</mml:mi><mml:mo>,</mml:mo><mml:mi>ν</mml:mi></mml:mrow><mml:mn>2</mml:mn></mml:msubsup></mml:mrow></mml:mfrac></mml:mrow></mml:msqrt><mml:mo>,</mml:mo><mml:msub><mml:mi>Z</mml:mi><mml:mrow><mml:mn>1</mml:mn><mml:mo>−</mml:mo><mml:mi>κ</mml:mi><mml:mo>/</mml:mo><mml:mn>2</mml:mn></mml:mrow></mml:msub><mml:msqrt><mml:mrow><mml:mfrac><mml:mrow><mml:mi>ν</mml:mi><mml:msubsup><mml:mover accent="true"><mml:mi>σ</mml:mi><mml:mo>^</mml:mo></mml:mover><mml:mi>ε</mml:mi><mml:mn>2</mml:mn></mml:msubsup></mml:mrow><mml:mrow><mml:msubsup><mml:mi>χ</mml:mi><mml:mrow><mml:mi>α</mml:mi><mml:mo>,</mml:mo><mml:mi>ν</mml:mi></mml:mrow><mml:mn>2</mml:mn></mml:msubsup></mml:mrow></mml:mfrac></mml:mrow></mml:msqrt></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:mo>.</mml:mo></mml:mrow></mml:semantics></mml:math></disp-formula>
<supplementary-material content-type="local-data" id="S1"><caption><title>Additional file 1</title><p>mix_model_MLPA.pdf, 173.4K.</p></caption></supplementary-material>
<table-wrap-foot><p>These results are for the case of having 20% of probes as a internal control probes (needed for the REX approach) and 10% of probes as a true altered probes. The results are given for different scenarios between probe variability (<italic>σ</italic><sub><italic>β</italic></sub>), probe-test variability (<italic>σ</italic><sub><italic>γ</italic></sub>) and within-probe variability (<italic>σ</italic><sub><italic>ε</italic></sub>). The column  indicates the mean number of simulated altered probes.</p></table-wrap-foot> <table-wrap-foot><p>These results are for the case of having 10% of probes as a internal control probes (needed for the REX approach) and 10% of probes as a true altered probes. The results are given for different scenarios between probe variability (<italic>σ</italic><sub><italic>β</italic></sub>), probe-test variability (<italic>σ</italic><sub><italic>γ</italic></sub>) and within-probe variability (<italic>σ</italic><sub><italic>ε</italic></sub>). The column  indicates the mean number of simulated altered probes.</p></table-wrap-foot> <table-wrap-foot><p>The three colums for each individual indicates the result obtained using mixed-model, threshold and REX-MLPA approaches, respectively. The code for the results are the following: -1: relative loss, 0:normal, 1:relative gain. Those result where a disagreement between the three methods is observed are in bold face.</p></table-wrap-foot>
<graphic xlink:href="1471-2105-9-261-1"/> <graphic xlink:href="1471-2105-9-261-2"/> <graphic xlink:href="1471-2105-9-261-3"/> <graphic xlink:href="1471-2105-9-261-4"/> <graphic xlink:href="1471-2105-9-261-5"/>
<media xlink:href="1471-2105-9-261-S1.pdf" mimetype="application" mime-subtype="pdf"><caption><p>Click here for file</p></caption></media>
[{"surname": ["Mavrogiannis", "Cockburn"], "given-names": ["LA", "DJ"], "article-title": ["Regression-Enhanced Analysis of Multiplex Ligation-Dependent Probe Amplification (REX-MLPA)"], "source": ["Tech rep, Yorkshire Regional DNA Laboratory"], "year": ["2004"]}, {"surname": ["Kellander", "Riley", "Liu"], "given-names": ["M", "M", "C"], "article-title": ["GeneMarker Software for Multiplex Ligation-dependent Probe Amplification (MLPA)"], "source": ["Tech rep, SoftGenetics LLC"], "year": ["2006"]}, {"article-title": ["CREAL's web-page"]}, {"article-title": ["NGRL-Manchester"]}, {"article-title": ["MRC-Holland"]}, {"surname": ["Motulsky", "Christopoulos"], "given-names": ["H", "A"], "source": ["Fitting Models to Biological Data Using Linear and Nonlinear Regression: a practical guide to curve fitting"], "year": ["2004"], "publisher-name": ["Oxford: Oxford Univesity Press"]}, {"surname": ["Pinheiro", "Bates"], "given-names": ["JC", "DM"], "source": ["Mixed-Effects Models in S and S-plus"], "year": ["2000"], "publisher-name": ["New York: Springer Verlag"]}]
BMC Bioinformatics. 2008 Jun 4; 9:261