<title>Background</title> <p>DNA methylation frequently occurs in mammalian DNA at the 5 position of cytosine in CpG dinucleotides. It has been estimated that over 70% of cytosines of CpG dinucleotides are methylated in the human genome. CpG dinucleotides are under-represented in the genome and methylated CpG dinucleotides predominantly occur within repetitive elements [##REF##2777259##2##]. However, there are CpG rich regions of the genome which generally remain unmethylated [##REF##2423876##3##]. These CpG rich regions are known as CpG islands and are frequently located in the promoter or the first exon regions of approximately 60% of all genes [##REF##7505451##4##]. The unmethylated status of CpG islands is thought to be a prerequisite state to maintain the linked gene in an active transcribed and transcriptional permissive state.</p> <p>Differential Methylation Hybridisation (DMH) is one of several techniques for examining CpG island methylation at a genome-wide scale that has been applied to the identification of aberrantly methylated gene promoters in various cancers [##REF##9949205##5##, ####REF##11731411##6##, ##REF##12114427##7##, ##REF##12163706##8##, ##REF##12095276##9##, ##REF##11951657##10##, ##REF##16900213##11##, ##REF##17504987##12####17504987##12##]. Nouzova <italic>et al</italic>[##REF##15302897##13##] modified the original method by using digestion with a methylation-dependent enzyme, M<italic>cr</italic>BC. This enzyme only cleaves methylated CpG DNA sequences. Within-sample comparison is applied after competitive hybridisation with M<italic>cr</italic>BC digested DNA and undigested (mock digested) DNA labelled with Cy3 and Cy5. If a locus is unmethylated the signal intensities of Cy3 and Cy5 are equivalent, while if methylated the Cy5/Cy3 (undigested/digested) ratio is greater than one. However, no common reference is generally used in the modified DMH method, and the unequal representation of methylated and unmethylated sequences due to competitive hybridisation may reduce sensitivity and specificity to detect differential methylation.</p> <p>Currently, Significance Analysis of Microarrays (SAM) [##REF##11309499##14##] and Prediction Analysis for Microarrays (PAM) [##REF##12011421##15##] are commonly applied in DNA methylation analysis. Based on the change in hybridisation relative to the standard deviation of repeated measurements, SAM assigns each gene a score that is an extension of the t-statistic. For significant genes with a score over a certain threshold, SAM uses permutations to estimate the false discovery rate (FDR). It has been implemented in many studies of gene expression data [##REF##17478416##16##, ####REF##17720884##17##, ##REF##18065728##18##, ##REF##17875689##19##, ##REF##17599952##20##, ##REF##17591626##21####17591626##21##] as well as DMH data, e.g. Wei <italic>et al</italic>. [##REF##16675572##22##] applied SAM to find the differential methylation of CpG island loci between ovarian caner patient groups with short and long progression-free survival (PFS). However, SAM assumes that the microarray data conform to approximate normality and symmetry, leading to the loss of power in the analysis of DMH data that are inherently skewed due to the biological features of DNA methylation in cancer and competitive hybridisation on DMH arrays (Figure ##FIG##0##1##).</p> <p>In the modified DMH method, the ratios of raw signal intensities (undigested/digested) greater than 1 reflect the various methylation levels [##REF##15302897##13##]. A ratio cut-off is generally used to identify the hypermethylated loci [##REF##12114427##7##]. However, this is an arbitrary value and does not necessarily accurately reflect the various sources of variation in the experiment. It is therefore desirable to develop an algorithm to more objectively assess the methylation status of loci from DMH data.</p> <p>PAM is a nearest centroid shrinkage method that identifies those genes that discriminate best between classes. This technique shrinks the class gene centroid towards the overall centroid by a "threshold" amount after standardizing each gene by its within class standard deviation. The "threshold" is identified by cross-validation. This approach was applied in the study by Wei <italic>et al</italic>. [##REF##16675572##22##] and showed certain power in the identification of differentially methylated loci, but PAM is designed for class prediction rather than class comparison. Although the class predictor used in PAM can reflect the difference between classes, a large number of loci actually differentially methylated between the classes are excluded to improve the accuracy of prediction.</p> <p>Although normalisation has become a standard procedure for the study of microarray data and is necessary for SAM and PAM analysis, unbalanced shifts in methylation status between class samples in DMH limit the use of between-class normalisation which assumes the changes are roughly symmetric. Thus, the differential methylation can be masked by the over-correction of normalisation and it would be preferable to use a method of analysis that does not require normalisation of the data.</p> <p>Since PAM and SAM may have limitations for analysing DMH data, we have developed an alternative approach based on the specific features and known biological properties of the arrays used for DMH analysis. The algorithm is named as Methylation Linear Discriminant Analysis (MLDA) and has been applied to identify a set of loci differentially methylated between ovarian cisplatin sensitive and resistant cancer cell lines.</p>
<title>Methods</title> <p>First, all intensity values were log transformed. A multiplicative background correction was applied to correct signal intensities for the background noise in each array. After background correction, the log-transformed digested and undigested intensities show three approximately parallel linear patterns (Figure ##FIG##1##2a##). The first pattern (digested/undigested is close to 1) represents the unmethylated sequences. The second pattern represents either hemi-methylated sequences or the unmethylated sequences cross-hybridised with the methylated ones on the panel. The third pattern represents the methylated sequences in target DNA. The methylated and unmethylated loci in target DNA can be characterised by a linear regression model for each pattern. The distance of each spot to the methylated and unmethylated lines respectively can then be estimated by standardised residuals. The log likelihood ratio of a locus being methylated is then proportional to the difference between the squared standardised residual from the methylated line and that from the unmethylated one. The algorithm based around this regression approach is named Methylation Linear Discriminant Analysis (MLDA) and was programmed in R version 2.7.0.</p> <title>Log-likelihood ratio transformation</title> <p>a An univariate linear regression model was constructed for the unmethylated probes (e.g. mitochondrial derived features) using formula (1) where <italic>α </italic>is the intercept, <italic>β </italic>is the slope of the model, and <italic>ξ </italic>is the error representing the unpredicted or unexplained variation in the model (Figure ##FIG##1##2b##). The parameters of regression line were estimated by the method of least squares (formula 2 and 3).</p> <p></p> <p></p> <p></p> <p><italic>k </italic>is the number of unmethylated controls on DMH array. <italic>G</italic><sub><italic>i </italic></sub>and <italic>R</italic><sub><italic>i </italic></sub>are the logarithmic-transformed digested and undigested intensities of microarray probes for mitochondrial sequences, respectively.  and  are the averaged logarithmic-transformed undigested and digested intensities of the k unmethylated controls.</p> <p>b. The scale estimate <italic>σ</italic><sub>mito </sub>associated with the error term in the linear regression model was estimated from the residuals from the observed <italic>k </italic>points to the fitted line. The most extreme 10% of residuals was omitted from either end of the distribution to minimise the impact of extreme residuals on this estimate.</p> <p>c. The standardised residuals of all the microarray probes to the unmethylation regression line were calculated as formula (4).</p> <p></p> <p>d. The point corresponding to the 97.5-quantiles residual below the unmethylation line is represented as X (R.975, G.975). The intermediate linear model (Figure ##FIG##1##2c##) was constructed through point X with a slope assumed to be 1 and the intercept estimated as formula (5).</p> <p></p> <p>e. The standardised residuals of all the microarray probes to the line with slope 1 and intercept estimated from (5) were calculated as formula (6). The variance of the residuals to the intermediate model was assumed to be similar as that in the mitochondrial model.</p> <p></p> <p>f. The microarray probes with standardised residuals less than 2 were included for later robust regression analysis. The line estimated from this regression analysis represents the methylation regression line (Figure ##FIG##1##2d##).</p> <p>g. The scale estimate <italic>σ</italic><sub>meth </sub>of the methylation regression line was estimated using only those microarray probes below the line, with the most extreme 5% removed.</p> <p>h. The standardised residuals of all the microarray probes to the methylated regression line were calculated as formula (7). The log likelihood ratio (LR) of all the microarray probes was estimated by formula (8) for further analysis.</p> <p></p> <p></p> <title>Determination of log likelihood ratio cut-offs</title> <p>Two inconsistency rates (IR<sub>meth </sub>and IR<sub>unmeth</sub>) and two consistency rates (CR<sub>meth </sub>and CR<sub>unmeth</sub>) between dye-swap arrays were used to determine the log like likelihood ratio threshold. IR<sub>meth </sub>(formula 9) represents the rate of the microarray probes identified as methylated in one array but as unmethylated in the other one, while IR<sub>unmeth </sub>(formula 10) is the rate of the microarray probes identified as unmethylated in one array but as methylated in the other one. CR<sub>meth </sub>(formula 11) and CR<sub>unmeth </sub>(formula 12) are the rates for the spots identified as methylated (CR<sub>meth</sub>) and unmethylated (CR<sub>unmeth</sub>) in both dye-swap arrays (Figure ##FIG##6##7##).</p> <p></p> <p></p> <p></p> <p></p> <p>b The log likelihood ratio thresholds (LR<sub>meth </sub>and LR<sub>unmeth</sub>) for methylated and unmethylated microarray probes, which kept the IR rates low (at or close to 1%) and the CR rates high (at or close to 140%), were used as the cut-offs for methylated and unmethylated loci. IR tends to rise with the increase of CR slowly, but starts to increase dramatically when the CR goes above 140%, at which point the inconsistency rate is generally about 1%. We have therefore used CR &gt; 140% and IR &lt; 1% as the criteria for determining the methylation cut-offs.</p> <title>Identification of robust regression outliers</title> <p>Each microarray probe was scored based on the cut-offs of likelihood ratios for methylation and unmethylation on dye-swap arrays using the weighted methylation scoring scheme shown in Figure ##FIG##2##3##. The microarray probes consistently identified as methylated candidates on dye-swap arrays were scored of 1; similarly unmethylated microarray probes were scored of -1. The rest of the microarray probes were assigned a weighted score based on their location on the plot.</p> <p>A robust regression model [##UREF##2##28##] was constructed with the averaged scores in one class of samples as the explanatory variable, and the corresponding scores in the other class of samples as the dependent variable. The degree of trimming was determined according to Barnett <italic>et al</italic>. [##UREF##3##29##] when estimating the variance of residuals to the robust linear regression model.</p> <p>It was assumed that the standardised residuals (SRs) from the robust regression line followed a normal distribution <italic>N</italic>(<italic>μ</italic>, <italic>σ</italic><sup>2</sup>). <italic>μ </italic>and <italic>σ </italic>were estimated excluding outliers using the MAD-Median Rule [##UREF##4##30##]. The p value for each SR cut-off was calculated as described by Simon <italic>et al </italic>[##UREF##1##25##]. This p-value reflects the probability of observing a group of more extreme residuals from the fitted normal distribution. Microarray probes were identified as outliers if their SRs were larger than the cut-off for which the p-value was less than 0.01.</p> <title>Estimation of misclassification rate</title> <p>The misclassification rate was estimated by drawing bootstrap samples 500 times with replacement from the two classes (sensitive and resistant) and carrying out hierarchical clustering based on the loci identified as differentially methylated using weighted scores for MLDA and log ratios without between-group normalisation for SAM and PAM, respectively. Clustering was carried out using Euclidean distance as the distance metric, and clusters were agglomerated using the average linkage criterion. The clustering tree was cut into two groups and the number of misclassified cell lines was counted. The misclassification rate was obtained from the averaged number of misclassified samples in 500 bootstraps divided by the total number of samples.</p> <title>SAM and PAM analysis</title> <p>The raw signal intensities of each channel were subtracted by the median signal intensities of corresponding channel of controls on HCGI12K array. After this correction, SAM in samr package and PAM in pamr package were applied using log ratios (digested/undigested) in R version 2.7.0. Between-group normalisation was not used in SAM and PAM to avoid over-correction masking the differential methylation.</p>
<title>Results</title> <title>Outline of MLDA</title> <p>In this study, we have developed a novel approach, named MLDA, for analysing CpG island microarray hybridisation data that allows the identification of differentially methylated loci. MLDA was programmed in R (version 2.7.0) and the package is available at CRAN [##UREF##0##1##]. This approach uses three relatively simple linear regression models. The first one is constructed by the log-transformed signal intensities of unmethylated features and used as the reference for unmethylation (Figure ##FIG##1##2b##). The second one is the intermediate model constructed through the point corresponding to the 97.5-quantiles residual below the first linear regression line (Figure ##FIG##1##2c##). The features with a standardised residual less than 2 from this intermediate model are used to generate the third model which is used as the reference for methylation (Figure ##FIG##1##2d##). The log likelihood ratio of a locus being methylated is then proportional to the difference between the squared standardised residual from the methylated line and that from the unmethylated line. The log likelihood threshold of zero then provides a more rational basis for distinguishing between methylated and unmethylated loci than a robust undigested/digested ratio of 1.5, as it takes into account the observed variability in the experiment.</p> <p>In our approach the consistency and inconsistency rates of log likelihood ratios on dye-swapped/duplicate arrays are used to determine methylation and unmethylation cut-offs, which keep the consistency rate (CR) relatively high (about 140%) and the inconsistency rate (IR) low (about 1%). Each loci is assigned a score based on the cut-offs using the weighted methylation scoring scheme. The feature consistently identified as methylated candidates on dye-swapped/duplicate arrays are scored as 1; similarly unmethylated features are scored as -1; the rest of the feature are assigned a weighted score corresponding to their location on the plot of log-likelihood ratios (Figure ##FIG##2##3##).</p> <p>The averaged score for each locus is calculated in each sample class (e.g. resistant or sensitive) and plotted against each other. A robust regression model is then fitted to these data. The standardised residuals from the robust regression model are assumed to follow a normal distribution <italic>N</italic>(<italic>μ</italic>, <italic>σ</italic><sup>2</sup>). The outliers of the standardised residuals are identified as the differentially methylated loci between the class samples.</p> <title>DMH dataset</title> <p>MLDA was applied to identify the CGIs differentially methylated from DMH data derived from sensitive A2780 derivatives (A2780, A2780p3, A2780p5, A2780p6, A2780p13, A2780p14) and isogenically matched, resistant lines [##REF##8640828##23##] derived by multiple exposures to cytotoxic levels of cisplatin and which are 2–5 fold resistant to cisplatin in clonogenic assays (A2780cp70, A2780/MCP1, A2780/MCP2, A2780/MCP3, A2780/MCP4, A2780/MCP5, A2780/MCP6, A2780/MCP7, A2780/MCP8, A2780/MCP9). After background correction, the log-transformed digested and undigested intensities of the 13056 microarray probes show three approximately parallel linear patterns (Figure ##FIG##1##2a##). The first pattern (digested/undigested is close to 1) represents the unmethylated sequences. The second pattern represents either hemi-methylated sequences or the unmethylated sequences cross-hybridised with the methylated ones on the panel. The third pattern represents the methylated sequences in target DNA. The methylated and unmethylated loci in target DNA can be characterised by a linear regression model for each pattern. As previously mentioned, normalisation may not be appropriate for DMH data, so the log ratios of signal intensities in two classes of samples are not at the same level (Figure ##FIG##3##4##). Normalisation is not required for MLDA as the determination of the methylation score is based on the data within each experiment.</p> <p>Mitochondrial DNA is unmethylated [##REF##15277367##24##], therefore, the signal intensities of both channels of microarray probes for mitochondrial sequences are expected to be equal. However, a bi-modal distribution is observed in the log-transformed fluorescence ratios (digested/undigested) of 121 mitochondrial sequences. The first peak represents the unmethylated mitochondrial sequences and the second lower peak is assumed to be the mitochondrial sequences cross-hybridised with other methylated sequences on the panel. Thus, we selected 94 of 121 mitochondrial sequences that were consistently unmethylated through all the cell lines and used them as the unmethylation reference in target DNA.</p> <p>The parameters of those two models in all 16 cell lines were estimated (Table ##TAB##0##1##). The slope of the unmethylated regression line constructed by 94 mitochondrial sequences is indeed close to 1. After computing the log-likelihood ratios, the methylation and unmethylation cut-offs and associated IRs and CRs were determined from the dye-swapped array pairs (details in Method section). As shown in Figure ##FIG##4##5##, IR tends to rise with the increase of CR slowly, but starts to increase dramatically when the CR goes above 140%, at which point IR is generally about 1%. We have therefore used CR &gt; 140% and IR &lt; 1% as the criteria for determining the methylation and unmethylation cut-offs. Each locus was scored using the weighted scoring scheme based on those cut-offs. The averaged scores in 6 cisplatin-sensitive cell lines and 10 cisplatin-resistant cell lines were used to construct a robust regression model. Figure ##FIG##5##6a## shows that the standardised residuals (residual/<italic>σ</italic>) from the robust regression model roughly follow a normal distribution. The positive and negative outliers are determined as described in Method section.</p> <p>Finally, 115 loci were identified as candidates differentially methylated between A2780 sensitive and these resistant cell lines (additional file ##SUPPL##0##1##). Noticeably, 113 of 115 loci (<italic>p </italic>= 8.8 × 10<sup>-3</sup>, outlier detection test [##UREF##1##25##]) were hypermethylated, but only 2 loci (<italic>p </italic>&lt; 0.001, outlier detection test) lost methylation in the resistant cell lines (Figure ##FIG##5##6b##). This is consistent with the unbalanced shift in DMH data and indicates cisplatin treatment of cells selects preferentially for hypermethylation of loci, rather than hypomethylation in these tumor cells.</p> <title>Validation of differential methylation</title> <p>To confirm the differential methylation of loci identified in this study, we experimentally tested the methylation of 26 loci by methylation-specific PCR (MSP) and/or pyrosequencing of bisulphite modified DNA [##REF##8790415##26##] in sensitive A2780 derivatives and cisplatin resistant derivatives. Twenty-three out of the 26 loci have been confirmed as differentially methylated (additional file ##SUPPL##0##1##). It should be noted that MSP and pyrosequencing only examine methylation at a limited number of CpG sites of the sequence present on the DMH analysis. It is possible that the loci which were not confirmed as differentially methylated are methylated at other CpG sites which are detected by DMH but not targeted by MSP and/or pyrosequencing primers and so 23 out of 26 loci confirmed as differentially methylated is a minimum estimation.</p> <p>To compare the results from MLDA, SAM and PAM, we analysed the DMH dataset by all three methods. MLDA identified 115 loci (113 hypermethylated and 2 hypomethylated loci, misclassification error &lt; 0.001), SAM identified 152 loci (149 hypermethylated and 3 hypomethylated loci, misclassification error = 0.227, FDR = 6.17 × 10<sup>-3</sup>), and PAM found 24 hypermethylated loci (misclassification error = 0.084, FDR &lt; 0.001) in the resistant cell lines. Twenty-four loci identified by all three methods are listed in Table ##TAB##1##2##.</p>
<title>Discussion</title> <p>Hypermethylation of promoter CpG islands is strongly correlated to transcriptional gene silencing and epigenetic maintenance of the silenced state and is a potential rich source of biomarkers of cancer. Differential Methylation Hybridisation (DMH) is one technique used for genome-wide DNA methylation analysis. The study of such microarray data sets should ideally account for the specific biological features of DNA methylation and the non-symmetrical distribution of the ratios of unmethylated and methylated sequences hybridised on the array. We have therefore developed a novel algorithm tailored to this type of data, Methylation Linear Discriminant Analysis (MLDA). MLDA utilises log likelihood ratios representing the relative probability that loci are methylated instead of log ratios of signal intensities used in previous studies [##REF##11731411##6##, ####REF##12114427##7##, ##REF##12163706##8##, ##REF##12095276##9##, ##REF##11951657##10####11951657##10##,##REF##12724229##27##]. Validation of 23/26 identified loci using independent methods of methylation analysis shows that MLDA can robustly identify differential methylated loci between ovarian cancer sensitive and resistant cell lines without requiring the data to be normalised.</p> <p>Although a log likelihood ratio above zero means that the locus tends to be methylated, we did not use zero as the cut-off to determine the number of methylated and unmethylated sequences, as the existence of cross-hybridisation and measurement errors in the DMH assay makes this unreliable. To increase the precision of the methylation classification, we used the inconsistency (IR) and consistency (CR) rates between the dye-swap arrays to determine likelihood ratio cut-offs for methylation and unmethylation and assigned each locus a methylation score based on the position relative to these cut-offs. As shown in Figure ##FIG##4##5##, not all cell lines can reach the point that CR is around 140% and IR is about 1%. IR and CR need to be carefully selected as the methylation scores of loci are consequently influenced by the change of IR and CR. We also observed a lower CR (about 120%) and a higher IR (about 2%) in another CpG island array using DMH (data not shown), therefore, further examination of what factors influence the achievable CR and IR rates may improve the utility of the MLDA approach.</p> <p>Data on methylation status for 121 mitochondrial derived sequences were available in this study. Mitochondrial sequences would be expected to be unmethylated. We used 94 mitochondrial sequences to construct unmethylated linear model at the beginning of the study, and indeed, 93 of 121 mitochondrial loci were defined as unmethylated and 25 loci being of uncertain methylation status by MLDA. However, three mitochondrial loci were identified as hypermethylated candidates in the resistant ovarian carcinoma cell lines by both MLDA and SAM. One explanation of this discrepancy is that all these three loci have more than one BLAT hit indicating the existence of homology with nuclear DNA sequences, raising the possibility of hybridisation with these nuclear DNA sequences which may be differentially methylated. As shown in Figure ##FIG##1##2a##, the loci in the middle pattern represent either hemi-methylated sequences or the unmethylated sequences cross-hybridised with the methylated ones on the panel. No specific allowance is made for these intermediate points in analysis by SAM and PAM, whereas MLDA attempts specifically to down-weight these points in the identification of the methylation regression line. By giving a lower weighted score (close to 0) (Figure ##FIG##2##3##) to those loci, MLDA reduces the influence of cross-hybridisation among this group of sequences. Of course cross-hybridisation may also occur in the loci in the other two patterns (methylated and unmethylated patterns), but it is not possible for any mathematical approach to identify this.</p> <p>The misclassification error of MLDA based on the methylation score is much lower than that for either SAM or PAM based on the log ratios, indicating the potential of MLDA methylation scores to be used as a reliable discriminator between classes of samples.</p>
<title>Conclusion</title> <p>We have developed a novel method, named MLDA, for genome-wide DNA methylation studies. MLDA can transform the signal intensities to log-likelihood ratios through three linear regression models. Using this approach MLDA allows determination of the methylation status of a locus based on dye-swapped/duplicate arrays. The method has been applied to assess the methylation status of each locus and identified 115 loci that exhibit differential methylation between A2780 sensitive and resistant cell lines. A minimum of 23 out of 26 loci have been confirmed by independent methods as differentially methylated.</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>Hypermethylation of promoter CpG islands is strongly correlated to transcriptional gene silencing and epigenetic maintenance of the silenced state. As well as its role in tumor development, CpG island methylation contributes to the acquisition of resistance to chemotherapy. Differential Methylation Hybridisation (DMH) is one technique used for genome-wide DNA methylation analysis. The study of such microarray data sets should ideally account for the specific biological features of DNA methylation and the non-symmetrical distribution of the ratios of unmethylated and methylated sequences hybridised on the array. We have therefore developed a novel algorithm tailored to this type of data, Methylation Linear Discriminant Analysis (MLDA).</p> <title>Results</title> <p>MLDA was programmed in R (version 2.7.0) and the package is available at CRAN [##UREF##0##1##]. This approach utilizes linear regression models of non-normalised hybridisation data to define methylation status. Log-transformed signal intensities of unmethylated controls on the microarray are used as a reference. The signal intensities of DNA samples digested with methylation sensitive restriction enzymes and mock digested are then transformed to the likelihood of a locus being methylated using this reference. We tested the ability of MLDA to identify loci differentially methylated as analysed by DMH between cisplatin sensitive and resistant ovarian cancer cell lines. MLDA identified 115 differentially methylated loci and 23 out of 26 of these loci have been independently validated by Methylation Specific PCR and/or bisulphite pyrosequencing.</p> <title>Conclusion</title> <p>MLDA has advantages for analyzing methylation data from CpG island microarrays, since there is a clear rational for the definition of methylation status, it uses DMH data without between-group normalisation and is less influenced by cross-hybridisation of loci. The MLDA algorithm successfully identified differentially methylated loci between two classes of samples analysed by DMH using CpG island microarrays.</p>
<title>Authors' contributions</title> <p>WD conducted the statistical analysis and algorithm development supervised by JP and RB. The DMH data was produced by JMT in collaboration with TH and PY. RB, JP and KV conceived the study. JMT, JG and CZ conducted validation by MSP or pyrosequencing in RB's lab. WD, JP and RB prepared the manuscript with review by all authors. Funding was obtained by RB and TH.</p> <title>Supplementary Material</title>
<title>Acknowledgements</title> <p>This work was supported by a CR-UK programme grant to RB, NIH R21 grant CA110475 to RB and TH, and an Ovarian Cancer Action project grant to RB, JMT and JG.</p>
<fig position="float" id="F1"><label>Figure 1</label><caption><p><bold>Distribution of log-transformed ratio of gene expression data in breast cancer and DMH data in A2780 cell line</bold>. The left histogram shows the distribution of log-transformed ratios (cy3/cy5) in gene expression profiling data from a previous study of breast cancer 36] which is symmetric, while the right histogram shows the log-transformed ratios (undigested/digested) of DMH data from the present study which is skewed.</p></caption></fig> <fig position="float" id="F2"><label>Figure 2</label><caption><p><bold>An illustration of unmethylated and methylated model construction in MLDA in A2780 cell line</bold>. a: Three patterns can be observed on the scatter plot of log-transformed Cy3 (undigested) against log-transformed Cy5 (digested) intensities. b: The unmethylated model constructed using 94 mitochondrial sequences as a unmethylation reference. c: The intermediate model constructed through the 97.5 quantile residual. The point X is the 97.5 quantile residual. The microarray probes colored in blue (standardised residual to the intermediate model is less than 2) are selected to construct the methylated model. d: Methylated (in blue) and unmethylated (in red) models in A2780 cell line.</p></caption></fig> <fig position="float" id="F3"><label>Figure 3</label><caption><p><bold>Weighted scoring scheme</bold>. The microarray probes consistently identified as methylated candidates on dye-swap arrays were scored 1; similarly unmethylated microarray probes were scored -1. The rest of the microarray probes were assigned a weighted score based on their location on the plot. LRmeth: log likelihood ratio cut-off for methylated loci; LRunmeth: log likelihood ratio cut-off for unmethylated loci. LR: log likelihood ratio on dye-swapped arrays.</p></caption></fig> <fig position="float" id="F4"><label>Figure 4</label><caption><p><bold>Box plot of log ratios of undigested signal intensities against digested signal intensities in 16 cell lines (dye-swapped arrays)</bold>. The boxes colored in red are the A2780 sensitive cell lines; in blue are the A2780 resistant cell lines. As normalisation is not applied, the center and scale of log ratios for the 16 cell lines are not at the same level.</p></caption></fig> <fig position="float" id="F5"><label>Figure 5</label><caption><p><bold>CR against IR in 16 cell lines</bold>. X axis is the consistency rate (CR) and y axis is the inconsistency rate (IR). IR tends to rise with the increase of CR slowly, but starts to increase dramatically when the CR goes above 140%, at which point the inconsistency rate is generally about 1%. Not all cell lines could reach this point e.g. MCP3.</p></caption></fig> <fig position="float" id="F6"><label>Figure 6</label><caption><p><bold>Outliers identifications</bold>. a: Distribution of the observed (histogram) standardised residuals and the theoretical distribution based on the fitted model (dashed smooth line in red). The red and blue solid line are the positive and negative cut-offs, respectively. b: Scatter plot of sensitive scores against resistant scores in A2780 series cell lines. The hypermethylated loci are colored in red and hypomethylated loci are in blue. The robust regression model is Y = 0.9956X + 0.0019.</p></caption></fig> <fig position="float" id="F7"><label>Figure 7</label><caption><p><bold>Determination of methylation and unmethylation cut-offs of likelihood ratios on dye-swapped arrays</bold>. LRmeth: log likelihood ratio cut-off for methylated spots; LRunmeth: log likelihood ratio cut-off for unmethylated spots.</p></caption></fig>
<table-wrap position="float" id="T1"><label>Table 1</label><caption><p>Parameters of linear models in MLDA for 16 cell lines in DMH dataset I</p></caption><table frame="hsides" rules="groups"><thead><tr><td align="left" colspan="9"><bold>Unmethylation linear regression model</bold></td></tr><tr><td align="left">cell line</td><td align="left">intercetp (<italic>α</italic>)</td><td align="left">slope (<italic>β</italic>)</td><td align="left"><italic>σ</italic></td><td align="left">R2</td><td align="left">intercept<sub>ds </sub>(<italic>α</italic>)</td><td align="left">slope<sub>ds </sub>(<italic>β</italic>)</td><td align="left"><italic>σ</italic><sub>ds</sub></td><td align="left">R2<sub>ds</sub></td></tr></thead><tbody><tr><td align="left">A2780</td><td align="left">-0.0003</td><td align="left">1.0122</td><td align="left">0.1727</td><td align="left">0.9829</td><td align="left">0.0005</td><td align="left">1.0574</td><td align="left">0.1897</td><td align="left">0.978</td></tr><tr><td align="left">A2780p3</td><td align="left">-0.0003</td><td align="left">1.0343</td><td align="left">0.1212</td><td align="left">0.9897</td><td align="left">0.0018</td><td align="left">1.1065</td><td align="left">0.1425</td><td align="left">0.9882</td></tr><tr><td align="left">A2780p5</td><td align="left">-0.0003</td><td align="left">1.0138</td><td align="left">0.1684</td><td align="left">0.984</td><td align="left">-0.0002</td><td align="left">1.0728</td><td align="left">0.1425</td><td align="left">0.9883</td></tr><tr><td align="left">A2780p6</td><td align="left">0.0002</td><td align="left">0.9914</td><td align="left">0.1605</td><td align="left">0.9778</td><td align="left">-0.0001</td><td align="left">1.0012</td><td align="left">0.1638</td><td align="left">0.9747</td></tr><tr><td align="left">A2780p13</td><td align="left">-0.0012</td><td align="left">1.024</td><td align="left">0.1628</td><td align="left">0.9786</td><td align="left">-0.0005</td><td align="left">1.0744</td><td align="left">0.1436</td><td align="left">0.9852</td></tr><tr><td align="left">A2780p14</td><td align="left">-0.0022</td><td align="left">1.0499</td><td align="left">0.1523</td><td align="left">0.9809</td><td align="left">-0.0009</td><td align="left">1.034</td><td align="left">0.2069</td><td align="left">0.9691</td></tr><tr><td align="left">A2780cp70</td><td align="left">0.0013</td><td align="left">0.9604</td><td align="left">0.2532</td><td align="left">0.9524</td><td align="left">-0.0002</td><td align="left">1.0119</td><td align="left">0.2402</td><td align="left">0.9479</td></tr><tr><td align="left">MCP1</td><td align="left">0.0002</td><td align="left">0.9946</td><td align="left">0.145</td><td align="left">0.9803</td><td align="left">-0.0023</td><td align="left">1.112</td><td align="left">0.1452</td><td align="left">0.9836</td></tr><tr><td align="left">MCP2</td><td align="left">0.0000</td><td align="left">0.998</td><td align="left">0.137</td><td align="left">0.9727</td><td align="left">-0.0028</td><td align="left">1.0475</td><td align="left">0.1719</td><td align="left">0.9653</td></tr><tr><td align="left">MCP3</td><td align="left">0.0004</td><td align="left">0.9932</td><td align="left">0.2253</td><td align="left">0.9183</td><td align="left">-0.0023</td><td align="left">1.0517</td><td align="left">0.2795</td><td align="left">0.8978</td></tr><tr><td align="left">MCP4</td><td align="left">0.0006</td><td align="left">0.9838</td><td align="left">0.1739</td><td align="left">0.9718</td><td align="left">-0.0028</td><td align="left">1.077</td><td align="left">0.1947</td><td align="left">0.9751</td></tr><tr><td align="left">MCP5</td><td align="left">0.0009</td><td align="left">0.9857</td><td align="left">0.2464</td><td align="left">0.9639</td><td align="left">-0.0008</td><td align="left">1.017</td><td align="left">0.2166</td><td align="left">0.9692</td></tr><tr><td align="left">MCP6</td><td align="left">-0.0022</td><td align="left">1.0352</td><td align="left">0.122</td><td align="left">0.9751</td><td align="left">-0.0068</td><td align="left">1.1283</td><td align="left">0.154</td><td align="left">0.9752</td></tr><tr><td align="left">MCP7</td><td align="left">-0.0005</td><td align="left">1.0079</td><td align="left">0.1379</td><td align="left">0.9791</td><td align="left">-0.0045</td><td align="left">1.1529</td><td align="left">0.1588</td><td align="left">0.9764</td></tr><tr><td align="left">MCP8</td><td align="left">-0.0028</td><td align="left">1.0578</td><td align="left">0.1903</td><td align="left">0.9431</td><td align="left">-0.0068</td><td align="left">1.1193</td><td align="left">0.1885</td><td align="left">0.9575</td></tr><tr><td align="left">MCP9</td><td align="left">-0.0017</td><td align="left">1.0331</td><td align="left">0.1834</td><td align="left">0.9614</td><td align="left">-0.0091</td><td align="left">1.1538</td><td align="left">0.1691</td><td align="left">0.9674</td></tr><tr><td colspan="9"><hr/></td></tr><tr><td align="left" colspan="9"><bold>Methylation linear regression model</bold></td></tr><tr><td align="left">cell line</td><td align="left">intercetp (<italic>α</italic>)</td><td align="left">slope (<italic>β</italic>)</td><td align="left"><italic>σ</italic></td><td align="left">R2</td><td align="left">intercept<sub>ds </sub>(<italic>α</italic>)</td><td align="left">slope<sub>ds </sub>(<italic>β</italic>)</td><td align="left"><italic>σ</italic><sub>ds</sub></td><td align="left">R2<sub>ds</sub></td></tr><tr><td colspan="9"><hr/></td></tr><tr><td align="left">A2780</td><td align="left">-0.8839</td><td align="left">0.8917</td><td align="left">0.1491</td><td align="left">0.9438</td><td align="left">-0.8055</td><td align="left">0.9086</td><td align="left">0.1706</td><td align="left">0.926</td></tr><tr><td align="left">A2780p3</td><td align="left">-1.1672</td><td align="left">0.9797</td><td align="left">0.1518</td><td align="left">0.9553</td><td align="left">-0.7414</td><td align="left">0.9774</td><td align="left">0.179</td><td align="left">0.9438</td></tr><tr><td align="left">A2780p5</td><td align="left">-0.8991</td><td align="left">0.8978</td><td align="left">0.1515</td><td align="left">0.9476</td><td align="left">-0.9246</td><td align="left">0.9766</td><td align="left">0.1673</td><td align="left">0.9523</td></tr><tr><td align="left">A2780p6</td><td align="left">-0.9455</td><td align="left">0.9562</td><td align="left">0.1838</td><td align="left">0.9378</td><td align="left">-1.1995</td><td align="left">0.9641</td><td align="left">0.1788</td><td align="left">0.9324</td></tr><tr><td align="left">A2780p13</td><td align="left">-1.8918</td><td align="left">0.9807</td><td align="left">0.2535</td><td align="left">0.8962</td><td align="left">-1.8049</td><td align="left">0.9652</td><td align="left">0.2936</td><td align="left">0.8512</td></tr><tr><td align="left">A2780p14</td><td align="left">-1.5637</td><td align="left">0.9142</td><td align="left">0.2549</td><td align="left">0.8857</td><td align="left">-1.4468</td><td align="left">0.9066</td><td align="left">0.2128</td><td align="left">0.8837</td></tr><tr><td align="left">A2780cp70</td><td align="left">-1.0317</td><td align="left">0.8501</td><td align="left">0.1581</td><td align="left">0.9115</td><td align="left">-1.3074</td><td align="left">0.8967</td><td align="left">0.1541</td><td align="left">0.9265</td></tr><tr><td align="left">MCP1</td><td align="left">-1.199</td><td align="left">0.9781</td><td align="left">0.1692</td><td align="left">0.9467</td><td align="left">-1.0935</td><td align="left">1.0384</td><td align="left">0.1775</td><td align="left">0.9525</td></tr><tr><td align="left">MCP2</td><td align="left">-0.8037</td><td align="left">0.9292</td><td align="left">0.1486</td><td align="left">0.9557</td><td align="left">-0.9738</td><td align="left">0.9381</td><td align="left">0.2176</td><td align="left">0.8848</td></tr><tr><td align="left">MCP3</td><td align="left">-1.1244</td><td align="left">0.9151</td><td align="left">0.1755</td><td align="left">0.9482</td><td align="left">-0.9303</td><td align="left">0.9205</td><td align="left">0.2599</td><td align="left">0.8098</td></tr><tr><td align="left">MCP4</td><td align="left">-1.4326</td><td align="left">0.9171</td><td align="left">0.1418</td><td align="left">0.966</td><td align="left">-1.6205</td><td align="left">0.961</td><td align="left">0.2348</td><td align="left">0.8323</td></tr><tr><td align="left">MCP5</td><td align="left">-1.1187</td><td align="left">0.9425</td><td align="left">0.1839</td><td align="left">0.9404</td><td align="left">-1.2007</td><td align="left">0.9546</td><td align="left">0.1757</td><td align="left">0.9295</td></tr><tr><td align="left">MCP6</td><td align="left">-1.246</td><td align="left">0.925</td><td align="left">0.1966</td><td align="left">0.9294</td><td align="left">-1.2182</td><td align="left">0.9826</td><td align="left">0.2248</td><td align="left">0.8989</td></tr><tr><td align="left">MCP7</td><td align="left">-1.8972</td><td align="left">0.9909</td><td align="left">0.1977</td><td align="left">0.9442</td><td align="left">-1.4894</td><td align="left">1.0139</td><td align="left">0.2458</td><td align="left">0.8886</td></tr><tr><td align="left">MCP8</td><td align="left">-1.0219</td><td align="left">0.9905</td><td align="left">0.1975</td><td align="left">0.9468</td><td align="left">-0.4735</td><td align="left">0.9421</td><td align="left">0.228</td><td align="left">0.8761</td></tr><tr><td align="left">MCP9</td><td align="left">-1.3399</td><td align="left">0.9837</td><td align="left">0.2073</td><td align="left">0.9352</td><td align="left">-1.1497</td><td align="left">1.0078</td><td align="left">0.1967</td><td align="left">0.9115</td></tr></tbody></table></table-wrap> <table-wrap position="float" id="T2"><label>Table 2</label><caption><p>24 loci identified by MLDA, PAM and SAM as differentially methylated candidates in the comparison between A2780 cisplatin sensitive and cisplatin multiple-selected resistant cell lines.</p></caption><table frame="hsides" rules="groups"><thead><tr><td align="center"><bold>microarray ID</bold></td><td align="center"><bold>status</bold></td><td align="center"><bold>validation</bold></td><td align="center"><bold>MLDA rank*</bold></td><td align="center"><bold>PAM rank</bold></td><td align="center"><bold>SAM rank</bold></td><td align="center"><bold>CGI***</bold></td><td align="center"><bold>gene symbol**</bold></td><td align="center"><bold>GenBank Accession</bold></td><td align="center"><bold>Chromosome</bold></td></tr></thead><tbody><tr><td align="center">66_G_6</td><td align="center">hypermethylated</td><td align="center">Yes</td><td align="center">1</td><td align="center">1</td><td align="center">1</td><td align="center">Yes</td><td/><td/><td/></tr><tr><td align="center">121_D_9</td><td align="center">hypermethylated</td><td align="center">Yes</td><td align="center">2</td><td align="center">5</td><td align="center">6</td><td align="center">Yes</td><td align="center">CRABP1</td><td align="center"><ext-link ext-link-type="gen" xlink:href="NM_004378">NM_004378</ext-link></td><td align="center">15</td></tr><tr><td align="center">39_E_1</td><td align="center">hypermethylated</td><td align="center">ND</td><td align="center">3</td><td align="center">2</td><td align="center">2</td><td align="center">No</td><td/><td/><td/></tr><tr><td align="center">122_D_9</td><td align="center">hypermethylated</td><td align="center">No</td><td align="center">4</td><td align="center">11</td><td align="center">11</td><td align="center">Yes</td><td align="center">SOX12</td><td align="center"><ext-link ext-link-type="gen" xlink:href="NM_006943">NM_006943</ext-link></td><td align="center">20</td></tr><tr><td align="center">123_D_9</td><td align="center">hypermethylated</td><td align="center">No</td><td align="center">5</td><td align="center">10</td><td align="center">10</td><td align="center">Yes</td><td align="center">SOX12</td><td align="center"><ext-link ext-link-type="gen" xlink:href="NM_006943">NM_006943</ext-link></td><td align="center">20</td></tr><tr><td align="center">51_H_8</td><td align="center">hypermethylated</td><td align="center">Yes</td><td align="center">6</td><td align="center">18</td><td align="center">19</td><td align="center">No</td><td align="center">FEZF2</td><td align="center"><ext-link ext-link-type="gen" xlink:href="NM_018008">NM_018008</ext-link></td><td align="center">3</td></tr><tr><td align="center">58_A_1</td><td align="center">hypermethylated</td><td align="center">ND</td><td align="center">7</td><td align="center">22</td><td align="center">22</td><td align="center">Yes</td><td/><td/><td/></tr><tr><td align="center">80_H_5</td><td align="center">hypermethylated</td><td align="center">ND</td><td align="center">8</td><td align="center">14</td><td align="center">16</td><td align="center">No</td><td/><td/><td/></tr><tr><td align="center">21_A_11</td><td align="center">hypermethylated</td><td align="center">Yes</td><td align="center">9</td><td align="center">17</td><td align="center">18</td><td align="center">Yes</td><td align="center">NTN4</td><td align="center"><ext-link ext-link-type="gen" xlink:href="NM_021229">NM_021229</ext-link></td><td align="center">12</td></tr><tr><td align="center">38_D_7</td><td align="center">hypermethylated</td><td align="center">Yes</td><td align="center">11</td><td align="center">23</td><td align="center">24</td><td align="center">Yes</td><td align="center">AGBL2</td><td align="center"><ext-link ext-link-type="gen" xlink:href="NM_024783">NM_024783</ext-link></td><td align="center">11</td></tr><tr><td align="center">40_E_1</td><td align="center">hypermethylated</td><td align="center">ND</td><td align="center">12</td><td align="center">19</td><td align="center">17</td><td align="center">No</td><td/><td/><td/></tr><tr><td align="center">18_A_7</td><td align="center">hypermethylated</td><td align="center">ND</td><td align="center">13</td><td align="center">20</td><td align="center">20</td><td align="center">No</td><td align="center">EDIL3</td><td align="center"><ext-link ext-link-type="gen" xlink:href="NM_005711">NM_005711</ext-link></td><td align="center">5</td></tr><tr><td align="center">55_F_8</td><td align="center">hypermethylated</td><td align="center">Yes</td><td align="center">14</td><td align="center">9</td><td align="center">9</td><td align="center">Yes</td><td align="center">BC127881</td><td align="center"><ext-link ext-link-type="gen" xlink:href="BC127881">BC127881</ext-link></td><td align="center">7</td></tr><tr><td align="center">122_B_1</td><td align="center">hypermethylated</td><td align="center">ND</td><td align="center">15</td><td align="center">16</td><td align="center">14</td><td align="center">No</td><td/><td/><td/></tr><tr><td align="center">109_A_6</td><td align="center">hypermethylated</td><td align="center">No</td><td align="center">18</td><td align="center">12</td><td align="center">13</td><td align="center">Yes</td><td/><td/><td/></tr><tr><td align="center">41_D_9</td><td align="center">hypermethylated</td><td align="center">Yes</td><td align="center">22</td><td align="center">8</td><td align="center">8</td><td align="center">Yes</td><td align="center">WNT1</td><td align="center"><ext-link ext-link-type="gen" xlink:href="NM_005430">NM_005430</ext-link></td><td align="center">12</td></tr><tr><td align="center">42_D_9</td><td align="center">hypermethylated</td><td align="center">ND</td><td align="center">23</td><td align="center">4</td><td align="center">4</td><td align="center">No</td><td/><td/><td/></tr><tr><td align="center">119_A_6</td><td align="center">hypermethylated</td><td align="center">Yes</td><td align="center">24</td><td align="center">6</td><td align="center">5</td><td align="center">Yes</td><td align="center">NR2E1</td><td align="center"><ext-link ext-link-type="gen" xlink:href="NM_003269">NM_003269</ext-link></td><td align="center">6</td></tr><tr><td align="center">63_A_8</td><td align="center">hypermethylated</td><td align="center">ND</td><td align="center">26</td><td align="center">15</td><td align="center">15</td><td align="center">No</td><td/><td/><td/></tr><tr><td align="center">6_D_4</td><td align="center">hypermethylated</td><td align="center">Yes</td><td align="center">31</td><td align="center">3</td><td align="center">3</td><td align="center">Yes</td><td align="center">LMX1A</td><td align="center"><ext-link ext-link-type="gen" xlink:href="NM_177398">NM_177398</ext-link></td><td align="center">1</td></tr><tr><td align="center">17_H_9</td><td align="center">hypermethylated</td><td align="center">Yes</td><td align="center">34</td><td align="center">13</td><td align="center">12</td><td align="center">Yes</td><td align="center">HRASLS3</td><td align="center"><ext-link ext-link-type="gen" xlink:href="NM_006290">NM_006290</ext-link></td><td align="center">6</td></tr><tr><td align="center">5_D_4</td><td align="center">hypermethylated</td><td align="center">Yes</td><td align="center">35</td><td align="center">7</td><td align="center">7</td><td align="center">Yes</td><td align="center">LMX1A</td><td align="center"><ext-link ext-link-type="gen" xlink:href="NM_177398">NM_177398</ext-link></td><td align="center">1</td></tr><tr><td align="center">24_D_3</td><td align="center">hypermethylated</td><td align="center">Yes</td><td align="center">75</td><td align="center">24</td><td align="center">23</td><td align="center">Yes</td><td align="center">SP5</td><td align="center"><ext-link ext-link-type="gen" xlink:href="NM_001003845">NM_001003845</ext-link></td><td align="center">2</td></tr><tr><td align="center">122_G_1</td><td align="center">hypermethylated</td><td align="center">ND</td><td align="center">101</td><td align="center">21</td><td align="center">21</td><td align="center">Yes</td><td/><td/><td/></tr></tbody></table></table-wrap> <table-wrap position="float" id="T3"><label>Table 3</label><caption><p>Classification of loci</p></caption><table frame="hsides" rules="groups"><thead><tr><td/><td align="center">Unmethylated</td><td align="center">uncertain</td><td align="center">methylated</td></tr></thead><tbody><tr><td align="left">Unmethylated</td><td align="right">a</td><td align="right">b</td><td align="right">c</td></tr><tr><td align="left">Uncertain</td><td align="right">d</td><td align="right">e</td><td align="right">f</td></tr><tr><td align="left">Methylated</td><td align="right">g</td><td align="right">h</td><td align="right">i</td></tr></tbody></table></table-wrap>
<disp-formula id="bmcM1"><label>(1)</label><italic>G</italic><sub><italic>i </italic></sub>= <italic>α </italic>+ <italic>βR</italic><sub><italic>i </italic></sub>+ <italic>ξ</italic><sub><italic>i </italic></sub><italic>i </italic>= 1,2,3.......<italic>k</italic></disp-formula> <disp-formula id="bmcM2"><label>(2)</label><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M1" name="1471-2105-9-337-i1" overflow="scroll"><mml:semantics><mml:mrow><mml:mtable><mml:mtr><mml:mtd><mml:mrow><mml:mover accent="true"><mml:mi>β</mml:mi><mml:mo>^</mml:mo></mml:mover><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:mstyle displaystyle="true"><mml:mo>∑</mml:mo><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:msub><mml:mi>R</mml:mi><mml:mi>i</mml:mi></mml:msub><mml:mo>−</mml:mo><mml:mover accent="true"><mml:mi>R</mml:mi><mml:mo>¯</mml:mo></mml:mover><mml:mo stretchy="false">)</mml:mo><mml:mo stretchy="false">(</mml:mo><mml:msub><mml:mi>G</mml:mi><mml:mi>i</mml:mi></mml:msub><mml:mo>−</mml:mo><mml:mover accent="true"><mml:mi>G</mml:mi><mml:mo>¯</mml:mo></mml:mover><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:mstyle></mml:mrow><mml:mrow><mml:mstyle displaystyle="true"><mml:mo>∑</mml:mo><mml:mrow><mml:msup><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:msub><mml:mi>R</mml:mi><mml:mi>i</mml:mi></mml:msub><mml:mo>−</mml:mo><mml:mover accent="true"><mml:mi>R</mml:mi><mml:mo>¯</mml:mo></mml:mover><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mn>2</mml:mn></mml:msup></mml:mrow></mml:mstyle></mml:mrow></mml:mfrac></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:mi>i</mml:mi><mml:mo>=</mml:mo><mml:mn>1</mml:mn><mml:mo>,</mml:mo><mml:mn>2</mml:mn><mml:mo>,</mml:mo><mml:mn>3.......</mml:mn><mml:mi>k</mml:mi></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:mrow></mml:semantics></mml:math></disp-formula> <disp-formula id="bmcM3"><label>(3)</label><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M2" name="1471-2105-9-337-i2" overflow="scroll"><mml:semantics><mml:mrow><mml:mover accent="true"><mml:mi>α</mml:mi><mml:mo>^</mml:mo></mml:mover><mml:mo>=</mml:mo><mml:mover accent="true"><mml:mi>β</mml:mi><mml:mo>^</mml:mo></mml:mover><mml:mover accent="true"><mml:mi>R</mml:mi><mml:mo>¯</mml:mo></mml:mover><mml:mo>−</mml:mo><mml:mover accent="true"><mml:mi>G</mml:mi><mml:mo>¯</mml:mo></mml:mover></mml:mrow></mml:semantics></mml:math></disp-formula> <inline-formula><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M3" name="1471-2105-9-337-i3" overflow="scroll"><mml:semantics><mml:mover accent="true"><mml:mi>G</mml:mi><mml:mo>¯</mml:mo></mml:mover></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-337-i4" overflow="scroll"><mml:semantics><mml:mover accent="true"><mml:mi>R</mml:mi><mml:mo>¯</mml:mo></mml:mover></mml:semantics></mml:math></inline-formula> <disp-formula id="bmcM4"><label>(4)</label><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M5" name="1471-2105-9-337-i5" overflow="scroll"><mml:semantics><mml:mrow><mml:mi>S</mml:mi><mml:msub><mml:mi>R</mml:mi><mml:mrow><mml:mi>m</mml:mi><mml:mi>i</mml:mi><mml:mi>t</mml:mi><mml:mi>o</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:mi>r</mml:mi><mml:mi>e</mml:mi><mml:mi>s</mml:mi><mml:mi>i</mml:mi><mml:mi>d</mml:mi><mml:mi>u</mml:mi><mml:mi>a</mml:mi><mml:mi>l</mml:mi><mml:msub><mml:mi>s</mml:mi><mml:mrow><mml:mi>m</mml:mi><mml:mi>i</mml:mi><mml:mi>t</mml:mi><mml:mi>o</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mrow><mml:msub><mml:mi>σ</mml:mi><mml:mrow><mml:mi>m</mml:mi><mml:mi>i</mml:mi><mml:mi>t</mml:mi><mml:mi>o</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:mfrac></mml:mrow></mml:semantics></mml:math></disp-formula> <disp-formula id="bmcM5"><label>(5)</label><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M6" name="1471-2105-9-337-i6" overflow="scroll"><mml:semantics><mml:mrow><mml:mover accent="true"><mml:mi>α</mml:mi><mml:mo>^</mml:mo></mml:mover><mml:mo>=</mml:mo><mml:msub><mml:mi>G</mml:mi><mml:mrow><mml:mn>.975</mml:mn></mml:mrow></mml:msub><mml:mo>−</mml:mo><mml:msub><mml:mi>R</mml:mi><mml:mrow><mml:mn>.975</mml:mn></mml:mrow></mml:msub><mml:mo>+</mml:mo><mml:mn>1.96</mml:mn><mml:msub><mml:mi>σ</mml:mi><mml:mrow><mml:mi>m</mml:mi><mml:mi>i</mml:mi><mml:mi>t</mml:mi><mml:mi>o</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:semantics></mml:math></disp-formula> <disp-formula id="bmcM6"><label>(6)</label><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M7" name="1471-2105-9-337-i7" overflow="scroll"><mml:semantics><mml:mrow><mml:mi>S</mml:mi><mml:msub><mml:mi>R</mml:mi><mml:mrow><mml:mn>.975</mml:mn></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:mi>r</mml:mi><mml:mi>e</mml:mi><mml:mi>s</mml:mi><mml:mi>i</mml:mi><mml:mi>d</mml:mi><mml:mi>u</mml:mi><mml:mi>a</mml:mi><mml:mi>l</mml:mi><mml:msub><mml:mi>s</mml:mi><mml:mrow><mml:mn>.975</mml:mn></mml:mrow></mml:msub></mml:mrow><mml:mrow><mml:msub><mml:mi>σ</mml:mi><mml:mrow><mml:mi>m</mml:mi><mml:mi>i</mml:mi><mml:mi>t</mml:mi><mml:mi>o</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:mfrac></mml:mrow></mml:semantics></mml:math></disp-formula> <disp-formula id="bmcM7"><label>(7)</label><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M8" name="1471-2105-9-337-i8" overflow="scroll"><mml:semantics><mml:mrow><mml:mi>S</mml:mi><mml:msub><mml:mi>R</mml:mi><mml:mrow><mml:mi>m</mml:mi><mml:mi>e</mml:mi><mml:mi>t</mml:mi><mml:mi>h</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:mi>r</mml:mi><mml:mi>e</mml:mi><mml:mi>s</mml:mi><mml:mi>i</mml:mi><mml:mi>d</mml:mi><mml:mi>u</mml:mi><mml:mi>a</mml:mi><mml:mi>l</mml:mi><mml:msub><mml:mi>s</mml:mi><mml:mrow><mml:mi>m</mml:mi><mml:mi>e</mml:mi><mml:mi>t</mml:mi><mml:mi>h</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mrow><mml:msub><mml:mi>σ</mml:mi><mml:mrow><mml:mi>m</mml:mi><mml:mi>e</mml:mi><mml:mi>t</mml:mi><mml:mi>h</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:mfrac></mml:mrow></mml:semantics></mml:math></disp-formula> <disp-formula id="bmcM8"><label>(8)</label><italic>LR </italic>= <italic>SR</italic><sup>2</sup><sub><italic>mito </italic></sub>- <italic>SR</italic><sup>2</sup><sub><italic>meth</italic></sub></disp-formula> <disp-formula id="bmcM9"><label>(9)</label><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M9" name="1471-2105-9-337-i9" overflow="scroll"><mml:semantics><mml:mrow><mml:mi>I</mml:mi><mml:msub><mml:mi>R</mml:mi><mml:mrow><mml:mi>m</mml:mi><mml:mi>e</mml:mi><mml:mi>t</mml:mi><mml:mi>h</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mfrac><mml:mi>c</mml:mi><mml:mrow><mml:mi>c</mml:mi><mml:mo>+</mml:mo><mml:mi>f</mml:mi><mml:mo>+</mml:mo><mml:mi>i</mml:mi></mml:mrow></mml:mfrac><mml:mo>+</mml:mo><mml:mfrac><mml:mi>g</mml:mi><mml:mrow><mml:mi>g</mml:mi><mml:mo>+</mml:mo><mml:mi>h</mml:mi><mml:mo>+</mml:mo><mml:mi>i</mml:mi></mml:mrow></mml:mfrac></mml:mrow></mml:semantics></mml:math></disp-formula> <disp-formula id="bmcM10"><label>(10)</label><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M10" name="1471-2105-9-337-i10" overflow="scroll"><mml:semantics><mml:mrow><mml:mi>I</mml:mi><mml:msub><mml:mi>R</mml:mi><mml:mrow><mml:mi>u</mml:mi><mml:mi>n</mml:mi><mml:mtext> </mml:mtext><mml:mi>m</mml:mi><mml:mi>e</mml:mi><mml:mi>t</mml:mi><mml:mi>h</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mfrac><mml:mi>g</mml:mi><mml:mrow><mml:mi>a</mml:mi><mml:mo>+</mml:mo><mml:mi>d</mml:mi><mml:mo>+</mml:mo><mml:mi>g</mml:mi></mml:mrow></mml:mfrac><mml:mo>+</mml:mo><mml:mfrac><mml:mi>c</mml:mi><mml:mrow><mml:mi>a</mml:mi><mml:mo>+</mml:mo><mml:mi>b</mml:mi><mml:mo>+</mml:mo><mml:mi>c</mml:mi></mml:mrow></mml:mfrac></mml:mrow></mml:semantics></mml:math></disp-formula> <disp-formula id="bmcM11"><label>(11)</label><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M11" name="1471-2105-9-337-i11" overflow="scroll"><mml:semantics><mml:mrow><mml:mi>C</mml:mi><mml:msub><mml:mi>R</mml:mi><mml:mrow><mml:mi>m</mml:mi><mml:mi>e</mml:mi><mml:mi>t</mml:mi><mml:mi>h</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mfrac><mml:mi>i</mml:mi><mml:mrow><mml:mi>c</mml:mi><mml:mo>+</mml:mo><mml:mi>f</mml:mi><mml:mo>+</mml:mo><mml:mi>i</mml:mi></mml:mrow></mml:mfrac><mml:mo>+</mml:mo><mml:mfrac><mml:mi>i</mml:mi><mml:mrow><mml:mi>g</mml:mi><mml:mo>+</mml:mo><mml:mi>h</mml:mi><mml:mo>+</mml:mo><mml:mi>i</mml:mi></mml:mrow></mml:mfrac></mml:mrow></mml:semantics></mml:math></disp-formula> <disp-formula id="bmcM12"><label>(12)</label><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M12" name="1471-2105-9-337-i12" overflow="scroll"><mml:semantics><mml:mrow><mml:mi>C</mml:mi><mml:msub><mml:mi>R</mml:mi><mml:mrow><mml:mi>m</mml:mi><mml:mi>e</mml:mi><mml:mi>t</mml:mi><mml:mi>h</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mfrac><mml:mi>a</mml:mi><mml:mrow><mml:mi>a</mml:mi><mml:mo>+</mml:mo><mml:mi>b</mml:mi><mml:mo>+</mml:mo><mml:mi>c</mml:mi></mml:mrow></mml:mfrac><mml:mo>+</mml:mo><mml:mfrac><mml:mi>a</mml:mi><mml:mrow><mml:mi>a</mml:mi><mml:mo>+</mml:mo><mml:mi>d</mml:mi><mml:mo>+</mml:mo><mml:mi>g</mml:mi></mml:mrow></mml:mfrac></mml:mrow></mml:semantics></mml:math></disp-formula>
<supplementary-material content-type="local-data" id="S1"><caption><title>Additional file 1</title><p>115 differential methylated candidates identified by MLDA in A2780 series cell lines.</p></caption></supplementary-material>
<table-wrap-foot><p>ds: dye swap</p><p><italic>σ</italic>: standard deviation</p><p>R<sup>2</sup>: coefficient of determination</p></table-wrap-foot> <table-wrap-foot><p>ND: not done. Yes: validated. No: not validated</p><p>MLDA rank*: the rank of standardised residuals to the robust regression line constructed by the averaged sensitive scores against averaged resistant scores</p><p>Gene symbol**: only the gene of which transcription start site (TSS) is within 5 kb span of the loci</p><p>CGI***: CpG island defined by Gardiner-Garden and Frommer [##REF##3656447##31##].</p></table-wrap-foot>
<graphic xlink:href="1471-2105-9-337-1"/> <graphic xlink:href="1471-2105-9-337-2"/> <graphic xlink:href="1471-2105-9-337-3"/> <graphic xlink:href="1471-2105-9-337-4"/> <graphic xlink:href="1471-2105-9-337-5"/> <graphic xlink:href="1471-2105-9-337-6"/> <graphic xlink:href="1471-2105-9-337-7"/>
<media xlink:href="1471-2105-9-337-S1.xls" mimetype="application" mime-subtype="vnd.ms-excel"><caption><p>Click here for file</p></caption></media>
[{}, {"surname": ["Simon", "Korn", "McShane", "Radmacher", "Wright", "Zhao"], "given-names": ["RM", "EL", "LM", "MD", "GW", "Y"], "article-title": ["Design and Analysis of DNA Microarray Investigations"], "source": ["Statistics for Biology and Health"], "year": ["2003"], "publisher-name": ["New York , Springer"]}, {"surname": ["Tukey"], "given-names": ["JW"], "source": ["Exploratory Data Analysis"], "year": ["1977"], "publisher-name": ["Rading Massachusetts , Addison-Wesley"]}, {"surname": ["Wilcox", "Holland BA"], "given-names": ["RR"], "article-title": ["Robust Regression"], "source": ["Robust Esitmation and Hypothesis Testing"], "year": ["2005"], "publisher-name": ["London , Elsevier Academic Press"]}, {"surname": ["Wilcox", "Holland BA"], "given-names": ["RR"], "article-title": ["Estimating Measures of Location and Scale"], "source": ["Robust Esitmation and Hypothesis Testing"], "year": ["2005"], "publisher-name": ["London , Elsevier Academic Press"]}]
BMC Bioinformatics. 2008 Aug 8; 9:337