Disagreement between two common biomarkers of global DNA methylation
Clinical Epigeneticsvolume 8, Article number: 60 (2016)
The quantification of global DNA methylation has been established in epigenetic screening. As more practicable alternatives to the HPLC-based gold standard, the methylation analysis of CpG islands in repeatable elements (LINE-1) and the luminometric methylation assay (LUMA) of overall 5-methylcytosine content in “CCGG” recognition sites are most widely used. Both methods are applied as virtually equivalent, despite the hints that their results only partly agree. This triggered the present agreement assessments.
Three different human cell types (cultured MCF7 and SHSY5Y cell lines treated with different chemical modulators of DNA methylation and whole blood drawn from pain patients and healthy volunteers) were submitted to the global DNA methylation assays employing LINE-1 or LUMA-based pyrosequencing measurements. The agreement between the two bioassays was assessed using generally accepted approaches to the statistics for laboratory method comparison studies. Although global DNA methylation levels measured by the two methods correlated, five different lines of statistical evidence consistently rejected the assumption of complete agreement. Specifically, a bias was observed between the two methods. In addition, both the magnitude and direction of bias were tissue-dependent. Interassay differences could be grouped based on Bayesian statistics, and these groups allowed in turn to re-identify the originating tissue.
Although providing partly correlated measurements of DNA methylation, interchangeability of the quantitative results obtained with LINE-1 and LUMA was jeopardized by a consistent bias between the results. Moreover, the present analyses strongly indicate a tissue specificity of the differences between the two methods.
Epigenetic effects are exerted by various factors such as early social experiences [1–3], physical training , age , nutritional or chemical factors such as royal jelly , benzene , asbestos, smoking , and drugs . For example, drugs may interfere with epigenetics  including all classical mechanisms such as histone modifications, DNA methylation [9, 10], and further regulatory processes of gene expression [11, 12]. This interference can be intended as with epigenetic therapeutics or unintended as common drugs may also exert epigenetic (side) effects . The consequences reach from none, via modulating the disease or introducing disease independent symptoms, up to a possible hereditability of epigenetic fingerprints provided that epigenetic transmission, so far only shown in plants , extends to animals. Moreover, the influences between drugs and epigenetics are mutual. Not only can drugs modulate epigenetics, but epigenetics is also increasingly recognized as a source of interindividual variability in drug effects . Quantification of epigenetic modulations has therefore manifold applications .
Assessing global DNA methylation is a frequent used marker for epigenetic screening. It captures the DNA methylation also at unknown genetic locations while the results of average DNA methylation correlate with the methylation of some trait-relevant genes . The gold standard in this field is HPLC quantification of the 5-methylcytosine content (methyl group added to the 5-carbon position of a cytosine) within the whole genome that occur mostly at CpG sites . However, due to its high demand in DNA amount and the difficulty to apply the method in high throughput approaches , alternatives have been developed . Widely accepted are (i) the cumulative analysis of methylation at CpG sites in repeatable elements such as long interspersed nuclear element 1 (LINE-1)  dispersed in more than 500,000 copies across the whole human genome [19–21] and (ii) the luminometric methylation assay (LUMA) of overall 5-mC content in “CmCGG” recognition sites across the whole genome .
The utility of both biomarkers is supported by the correlation of their results with the HPLC gold standard [23, 24]. However, both methods address different recognition sites within the genome. Specifically, while LINE-1 is overrepresented in A+T rich regions, LUMA recognition sites are more dense in C+G rich regions . Nevertheless, both methods are widely used as markers of global methylation [26–36] as if they were equivalent. Consistent with their biological differences, however, doubts have been raised about such equivalence. Indeed, the two methods quantified the global DNA methylation in colon biopsies, cell lines , and human blood cells  differently. However, this was based on limited sample sizes [23, 24] and DNA methylation ranges . Considering the increasing importance of the assessment of unknown epigenetic effects such as of drugs  or alimentary materials , where a definite set of target genes for such epigenetic effects can often not be predefined, the present study aimed at systematic assessment of the agreement between the two bioassays. With a clinical focus analogously to a previous investigation , the analyses were performed in human blood cells that are frequently the only easily available biological material in human studies [26–36]. However, to increase the range of global DNA methylation, human-derived cell populations (MCF7 and SHSY5Y cell lines) were added following treatment with methylation modifying substances. This provided a total of 238 samples. The underlying hypothesis of the present method comparison was the non-agreement between the two assays, as suggested by the biological differences of their recognition sites.
Study design and subjects
The two different global DNA methylation markers (LINE-1 pyrosequencing, LUMA) were assessed in three independent sample sets that were generated (i) in vitro from human cell lines (MCF7 and SHSY5Y) or (ii) whole blood samples acquired from healthy volunteers (one set) or (iii) pain patients (two sets). The in vitro approach served to induce a broad variation of DNA methylation as the required suitable basis for correlation analysis using known modulators of DNA methylation to the cell culture under highly controlled laboratory conditions. The LINE-1-based data had been analyzed previously in a non-redundant context addressing the effect of drug exposure on DNA methylation . To picture the clinical setting, whole blood-derived DNA collected from either healthy subjects or chronic pain patients was analyzed. The in vivo assessments followed the Declaration of Helsinki and were approved by the Ethics Committee of the Goethe University, Frankfurt am Main, Germany. Informed written consent from each participating subject had been obtained. The actual health status of the healthy volunteers was ascertained by medical history and physical examination including vital signs. Exclusion criteria were a current clinical condition, any other actual diseases, and drug intake within a week except oral contraceptives. Patient’s samples were available from a previous assessment of chronic pain patients treated with either opioid or non-opioid analgesics .
DNA sample acquisition
Human blood samples
DNA from whole blood samples and opioid-related phenotypes was available from previous studies [40–42]. Cohort 1 consisted of 83 (26 men, 57 women, aged 39.6 ± 7.02) healthy subjects that were drawn randomly from a control cohort. Cohorts 2 and 3 comprised pain patients in tertiary care. Cohort 2 contained 29 (14 men, 15 women, aged 47.8 ± 7.36) pain patients with an opioid treatment duration of ≥1 year and an average daily opioid dose of 52.69 ± 22.11 mg of oral morphine equivalents [43, 44]. Cohort 3 consisted of 19 (2 men, 17 women, aged 45.7 ± 11.63) pain patients who had received no opioids during their analgesic therapy (Tables 1 and 2). Sex and age matching was not possible due to insufficient human material for LINE-1 and LUMA assessments each in two independent measurements.
Cell culture and exposure to known modulators of DNA methylation
Since the human blood samples did not provide sufficient variability in DNA methylation for the present analyses, human cell lines were added in which the desired broad range of methylation could be induced by treating the cells with known modulators of DNA methylation. Two human cell lines were chosen, i.e., MCF7 cells that is a breast cancer-derived cell line and SHSY5Y cells that possess a neuronal character. The choice was based on the criteria (i) human origin consistent with the present clinical focus, (ii) reported ability to respond with decreased DNA methylation to the stimulation with the demethylating agent 5-Aza-2′-deoxycytidine (5-Aza-CdR), which can be inhibited by simultaneous treatment with S-adenosyl methionine (SAM).
MCF7 cells were cultured in Dulbecco’s modified Eagle’s medium (DMEM) + GlutaMax™ (Gibco, Darmstadt, Germany) supplemented with 10 % (v/v) fetal calf serum (FCS) and 1 % penicillin/streptomycin (PAA, Cölbe, Germany). SHSY5Y cells were obtained from the DZMS Collection of Microorganisms and Cell Cultures (Braunschweig, Germany) and were grown in 1:1 mixed Ham F12 and Minimum Essential Medium (MEM) (Gibco, Darmstadt, Germany) supplemented with 15 % (v/v) FCS, 1 % Minimum Essential Medium Non-Essential Amino Acids (MEM NEAA, Gibco, Darmstadt, Germany), 2 mM L-glutamine (Gibco, Darmstadt, Germany), and 1 % penicillin/streptomycin (PAA, Cölbe, Germany) at 37 °C in humidified atmosphere containing 5 % CO2. In prior to drug treatment, cells were seeded at a density of 3 × 105/10 cm2 and allowed to settle for 24 h in the complete media. Subsequently, cells were incubated for 72 h (MCF7, SHSY5Y) or 7 days (SHSY5Y) in the presence of the known or potential modulators of DNA methylation at various concentrations Table 1.
Drug treatment conditions have been reported elsewhere . In brief, 5-Aza-2′-deoxycytidine (5-Aza-CdR), a cytidine analogue, covalently traps DNMTs, and RG108, a specific DNMT inhibitor, directly blocks the active site of the enzyme which are expected to lead to global DNA hypomethylation [45–48]. S-adenosyl methionine (SAM) is a methyl donor that is catalyzed by DNMTs to form 5-methyl cytosine at CpG sites , thereby, it is expected to increase DNA methylation or at least inhibit global hypomethylation induced by 5-Aza-CdR . 2,4-Dichlorophenol (DCP) is an environmental pollutant reported to increase global methylation . Methadone was chosen as opioid because it had been involved in the largest group of patients (heroin addicts) in whom the clinical association of opioid-induced hypermethylation had been observed . On every day, media were replaced and compounds were added freshly. Methadone hydrochloride (Fagron, Barsbüttel, Germany) was dissolved in Dulbecco’s phosphate-buffered saline (DPBS) without CaCl2 and MgCl2 (Gibco, Germany, Darmstadt; 14190-094). 5-Aza-CdR, SAM, DCP (Sigma-Aldrich, Taufkirchen, Germany), and RG108 (Biomol, Hamburg, Germany) were dissolved in DMSO and mixed with solvent to obtain a final concentration of 0.1 % DMSO (0.25 % for RG108) to the cell media during incubation. Cells incubated with 0.1 % solvent alone or without any substance addition (i.e., the control condition) served as controls.
Quantification of global DNA methylation
Genomic DNA was extracted from cell line materials and whole blood samples with the DNeasy Blood and Tissue Kit (Qiagen, Hilden, Germany) according to the manufacturer’s protocol and eluted in water. Genomic DNA obtained from blood samples was concentrated using vacuum rotation (45 °C for 25 min) to reach at least a final concentration of 50 ng/μl.
Methylation analysis of retrotransposon LINE-1
The analysis of LINE-1 DNA methylation was performed identically as described previously in full detail . Bisulfite treatment was performed using the EZ DNA Methylation-Gold Kit (Zymo Research, Freiburg, Germany) with 0.5–1 μg genomic DNA as instructed by the manufacturer.
The analyzed region of a CpG island located in the promoter region (L1Hs) DNA (PubMed GenBank X58075.1; lower strand) has the bisulfite-converted sequence 5′-TTTTG AGTTAGGTGTGGGATATA GTTTYGTGGTGYGTYGTTTTTTAAGTYGGTTTGAAAAGCTAATATTCGGGTGGGAGTGATTCGATTTTTTAGGTGCGTTCGTTATTTTTTTTTTTGATTCGGAAAGGGAATTTTTTGATTTT-3′ where the 146-bp PCR product contains four analyzed CpG methylation sites (bold) and annealing sites for the PCR primers (underlined) and the sequencing primer (italic), respectively [7, 52]. PCR reactions were run on a Mastercycler nexus gradient flexlid device (Eppendorf, Hamburg, Germany) in a 50-μl reaction volume including 5-μl bisulfite-treated DNA, mixed with 0.5 μl MyTaq™ HS DNA Polymerase (5 U/μl) (Bioline, Luckenwalde, Germany), 10 μl 5× MyTaq Reaction Buffer, 0.2 μl of each PCR primer (100 μM), and 34.1 μl HPLC-purified water. The following PCR program was used: 95 °C for 1 min, 40 amplification cycles at 95 °C for 15 s, 56 °C for 15 s, 72 °C for 15 s, and a final elongation step at 72 °C for 5 min.
The analysis of the global methylation marker LINE-1 was done by means of Pyrosquencing™ (Qiagen, Hilden, Germany) as described previously [7, 38, 52]. In brief, 50 μl of the PCR templates were processed and purified with the PyroMark Vacuum Prep Worktable (Biotage, Uppsala Sweden) and subsequently annealed to the sequencing primer (5′-AGTTAGGTGTGGGATATAGT-3′) at 80 °C for 2 min as instructed by the manufacturer.
Sequence analysis took place on a PSQ 96 MA System using the PyroMark Gold Q96 Reagents (Qiagen, Hilden, Germany) with the sequence to analyze TTYGTGGTGYGTYGTTTTTTAAGTYGGTTT. Pyro Q-CpG methylation software (version 1.0.9) had been used to determine the nucleotide dispensation order (ATCAGTGTGTCAGTCAGTCTAGTCTG). LINE-1 methylation values represent the mean percentage methylation across all four CpG sites, which were measured in duplicate samples within one run. In addition, each sample was measured in two independent runs, which were subsequently averaged.
The accuracy of the analyses was verified by adding positive and negative control samples. Specifically, each run included control DNA from the EpiTect PCR Control DNA Set (Qiagen, Hilden, Germany) that contained both bisulfite-converted 100 % methylated and completely unmethylated DNA as positive controls and unconverted unmethylated DNA as negative control. The bisulfite-converted methylated control DNA reached on average 75.08 ± 0.68 % methylation while the bisulfite-converted unmethylated control DNA reached only 3.37 ± 0.21 % methylation, which agrees with published values . The negative PCR control did not show specific spikes for any injected nucleotide, which demonstrated assay specificity. All absolute methylation values were subsequently calibrated to the methylated and unmethylated control DNA to cover a range from 0 to 100 %. Non-CpG cytosine residues were used as built-in controls to verify the bisulfite conversion. The acceptable percentages for passed and checked quality were adjusted to the complete bisulfite-converted controls supplied by Qiagen. Samples not meeting the criteria for complete bisulfite conversion or pyrosequencing™ quality control checks were excluded. The interassay coefficients of variation for duplicates were 2.38 % for cell-line samples and 1.18 % for blood samples.
Luminometric methylation assay
The luminometric methylation assay (LUMA) was performed as described previously  with modifications previously proposed . A common used isoschizomer pair to investigate global DNA methylation pattern is HpaII and MspI; HpaII digestion is inhibited if the internal cytosine is methylated (CmCGG) at recognition site whereas MspI is insensitive to CpG methylation within this sequence . DNA methylation level is defined as the HpaII/MspI ratio that would be 1.0 if the DNA is completely unmethylated and would approach zero if the DNA is completely methylated [22, 53]. Because of reported star activity of EcoRI , we used MfeI, a methylation-insensitive restriction enzyme, as normalization reference. Four hundred to five hundred nanograms of genomic DNA was cleaved with either HapII + MfeI or MspI + MfeI in two separate 20 μl reactions containing 2 μl of 10× Tango Buffer (330 mM Tris-acetate, 100 mM Mg-acetate, 660 mM K-acetate, 0.1 mg/ml BSA, Thermo Scientific, Schwerte, Germany), 5 U of HpaII (10 U/μl; NEB, Frankfurt, Germany) or MspI (20 U/μl; NEB, Frankfurt, Germany), and 2.5 U of MfeI (10 U/μl; NEB, Frankfurt, Germany) at 37 °C for 16 h using a PSQ 96 Plate Low (Qiagen, Hilden, Germany). The incubation time could be reduced to 4 h without impacting the completion of the enzymatic reaction. Subsequent to digestion, 20 μl of annealing buffer (Qiagen, Hilden, Germany) was added to the cleavage reactions and samples were assayed in duplicate using the PSQ 96 MA System (Biotage AB, Uppsala, Sweden) and PyroMArk Gold Q96 reagents (Qiagen, Hilden, Germany). The sequence AC/TCGA was analyzed in SNP mode with ACTCGA nucleotide dispensation order. The dispensation order of dNTPs were dATPαS (step 1); mixture of dGTP + dCTP (step 2); dTTP (step 3); mixture of dGTP + dCTP (step 4); water (step 5); and dATPαS (step 6). Peak heights were calculated using the PyroMark™ ID software and HpaII/MfeI and MspI/MfeI ratios were determined as (dGTP + dCTP)/mean(dATP,dTTP) for each reaction. The HpaII/MspI ratio was then calculated as (HpaII/MfeI)/(MspI/MfeI), and methylation level was obtained as Methylation (%) = (1-HpaII/MspI) × 100. Samples with peak heights <2 (blood samples) or <1 (cell samples), MspI/MfeI ratio >4.2, and peak heights at dispensation peak 6 of more than 25 % relative to dispensation peak 1 were excluded from the analysis (modified ). The interassay coefficients of variation for duplicates were 3.28 % for cell samples and 2.17 % for blood samples.
The accuracy of the analysis was verified by including in each run an unmethylated probe of lambda phage DNA as 0 % control and a completely methylated probe of lambda phage DNA as 100 % control. All absolute methylation values measured in the three subsets of DNA samples, respectively, in human cell types were calibrated to the methylated and unmethylated lambda phage DNA control to cover the range of 0–100 %.
The data analysis employed several bioinformatics methods to assess the agreement between the percentages of DNA methylation quantified by the LINE-1- or LUMA-based method. The analyses cover and extend previously proposed approaches to the statistics for laboratory method comparison studies . It included (i) standard analysis of variance and correlation, (ii) visual inspection, (iii) the method of 95 % limits of agreement between measurements by two methods, (iv) Gauss mixture modeling, and (v) linear regression.
Analysis of variance and correlation, visual inspection
The first four analyses were performed using the R (version 3.2.1 for Linux; http://CRAN.R-project.org/) and SPSS (version 23 for Linux, IBM SPSS Statistics, Chicago, USA) software environments on an Intel Xeon® computer running on Ubuntu Linux 14.04 64-bit. In a first analytical approach, differences between DNA methylation assessed either by LINE-1 methylation or by LUMA were explored by submitting the data to analysis of variance for repeated measures (rm-ANOVA). “LINE-1/LUMA” was taken as within-subject factor and “data subset” (n = 3, Table 1) as between-subject factor, with post hoc Wilcoxon signed rank test-based  exploration of single differences. The α level was set at 0.05 and corrected for multiple testing according to Bonferroni . Additional statistics included nonparametric correlation analyses calculating Spearman’s ρ , for which 95 % confidence intervals (CI) were obtained using 1000 bootstrap resamplings . This was followed by the second approach, visual inspection of the scatter plot of the data and its placement relative to the line of equality.
Assessment of method agreement and bias
Absent correlation would discourage an agreement between the two assays. However, as pointed out previously, correlation analysis assesses the degree of association rather than the agreement between the methods and is insensitive to a possible bias . Therefore, a third analytical approach employed the method of 95 % limits of agreement between measurements by two methods proposed by Bland and Altman . For each data subset, differences in DNA methylation between the LINE-1- and LUMA-obtained magnitudes of DNA methylation were plotted against the mean of the two measurements (Fig. 3). The mean difference was an estimate of the fixed bias and tested for significant deviation from 0 on the basis of a one-sample t test. The 95 % confidence interval of the differences marked the limits of agreement for the two methods. A linear regression of the difference between the methods against their average indicated a relationship of the discrepancies between the measurements and the true value, which in the case of a slope significantly differing from 0 denoted the proportional bias. Calculations were performed using the R packages “BlandAltmanLeh” (B. Lehnert, https://cran.r-project.org/web/packages/BlandAltmanLeh/index.html) and “epade” (A. Schulz, https://cran.r-project.org/web/packages/epade/). Normality of the distribution of the differences between the two methods was assessed by means of Kolmogorov-Smirnov tests .
Pattern analysis of interassay differences
A data-subset specificity of the above differences was explored by fitting a mixture of Gaussian distributions (Gaussian mixture model (GMM) to their empirical distribution (Pareto density estimation (PDE)) as given by the equation
where N(x|Mean i , SD i ) denotes Gaussian probability densities with means, Means i , and standard deviations, SD i , while the w i is the mixture weights indicating the relative contribution of each component Gaussian to the overall distribution, and M denotes the number of components in the mixture. GMM fitting was performed with our R package “AdaptGauss” (M. Thrun, https://cran.r-project.org/web/packages/AdaptGauss/index.html; ), using the root mean square error between PDE and GMM as the fit criterion. The limits between the different Gaussians are defined by Bayes decision boundaries , i.e., the probability a data point being assigned to a specific Gaussian was calculated by an application of Bayes’ theorem , and the resulting grouping of the data was subsequently explored for association with data subsets, respectively, tissues, by applying a decision-tree algorithm  that used the information index, f(p) = − p ⋅ ln(p), to find optimal (local) dichotomic decisions. The method is invariant under transformations of the variables, robust with respect to outliers, and allows estimation of the misclassification rate . The resulting tree model was cross-validated using a leave-k-out approach, where k was a randomly picked tenth of the total sample and the tree models were built 100 times with the respective remaining data. Calculations were done using the “rpart” function of the similarly named R package (B. Ripley; https://cran.r-project.org/web/packages/rpart/index.html).
The tissue-dependent relation between LINE-1 and LUMA measurements was further explored in a fifth analytical approach that employed linear modeling performed with the non-linear mixed effects modeling (NONMEM) software (version 7.3, Icon, Dublin, Ireland ). The analysis searched data-subset specific deviations from a y-intersection of zero and a slope of the value of one of the linear model expressed by extending its reduced form of
θ Int,Subset1.. Subset3 and θ Slope,Subset1.. Subset3
The θs were allowed to take values differing from 0 or 1, for intersections and slopes, respectively, for each specific data subset while the values of the θs were fixed at values of 0 or 1, respectively, for all other subsets. For example, for the MCF7 cell line, θ Int,Subset1and θ Slope,Subset1 described the deviations of the linear relationship from the other samples. The parameters were fit only for the MCF7 data while they remained fixed at 0 or 1 when other data was analyzed. The full linear model thus consisted of eight structural parameters θ, of which θ 1 and θ 2 denoted the global intersection and slope of the linear relationship, and θ 3..8 accounted for the set-specific deviations from this global relationship and an additive residual error modeled as MethylationObserved = MethylationPredicted + ε, in which ε is a parameter with a mean of zero and a variance of σ 2.
During the fitting process, parameters θ Int,Subset1.. Subset3 and θ Slope,Subset1.. Subset3 were introduced into the model in a stepwise fashion. Whether or not a specific θ remained part of the final model was established based on goodness-of-fit assessments, i.e., Occam’s razor or the principle of parsimony was applied. The simpler model was preferred to the more complex model as long as an additional parameter did not provide a significantly better fit. The main test was a likelihood ratio test, and therefore, the indicator of improvement of the fit was a change in minus twofold the log likelihood (Δ−2LL), and the χ 2 approximation with the number of degrees of freedom equal to the difference in the number of parameters between two models was applied to judge statistical significance. Thus, the full model included an additional term and the reduced model involved the fixing of the respective term to a neutral value, i.e., 1 for factors and exponents and zero for summands. The α level was set at 0.05, which implies a significance criterion of Δ−2LL <-3.84; for further details of the fitting process refer to . Calculations were performed using “first order conditional estimation” .
Confidence intervals (95 %) of parameter values were calculated from 1000 runs of the final model with data sets that were obtained by bootstrap resampling  from the original data set , using PDxPop (version 5.10, Icon, Dublin, Ireland) for NONMEM. The limits of the 95 % confidence intervals of the parameter estimates were obtained as the 2.5th and 97.5th percentiles of the results of the 1000 model runs.
Five applied methods were consistent with finding a disagreement between the magnitudes of DNA methylation measured with the LINE-1 or the LUMA method (Fig. 1). The first analytical approach, standard analysis of variance, identified differences between the assays that additionally differed among the three data subsets. A second approach was visual in nature. By looking at the scatterplot in Fig. 2, it is clear that the majority of DNA methylation values measured by the LINE-1 sequence lay above or below the line of equality. The third approach, the Bland and Altman method (Fig. 3), verified this observation by finding a fixed bias in all three data subsets. This in its extent additionally appeared to differ among data subsets. The fourth analytical approach (Fig. 4), Gauss mixture modeling (GMM), found groups of interassay differences that partially allowed identifying the underlying tissue. The fifth analytical approach, linear regression (Fig. 2), substantiated the observation of a tissue dependency of the relationship between the measurements taken by the two methods. The results of the different analyses will be reported in detail in the following.
Visual inspection and analyses of variance and correlation
The distribution of the global DNA methylation of human blood, MCF7 cells, and SHSY5Y cells indicated differences between the two methylation markers LINE-1 and LUMA (Fig. 1). In MCF7 cells, the global DNA methylation appeared to be slightly smaller when assessed using LINE-1 than when assessed using LUMA. In the other data subsets, the opposite was observed, i.e., DNA methylation appeared to be slightly larger when assessed using LINE-1 than when assessed using LUMA. Data-subset specific interassay differences were substantiated by significant main effects of the rm-ANOVA factors “data subset” (df = 2,235, F = 200.11, p = 1.81 × 10−51) and “LINE1/LUMA” (df = 1,235, F = 252.57, p = 4.11 × 10−39) and by a significant interaction “LINE1/LUMA” by “data subset” (df = 2,235, F = 823.72, p = 6.62 × 10−107). The differences were statistically significant for all data subsets (Wilcoxon signed rank tests: all p value <0.001, Table 2). Additional visual inspection of the scatterplot of the LINE-1 versus the LUMA measurements (Fig. 2) indicated that only data acquired in SHSY5Y cells was scattered around the line of equality while samples acquired in MCF7 cells were located below and those acquired in human blood samples above that line. Finally, a statistically significant correlation of the DNA methylation between the two assays was observed in all data subsets (Table 2), however, only weakly in blood samples.
Method agreement and bias
Applying the method of 95 % limits of agreement identified significant bias between the results obtained with the two bioassays. This was observed as a deviation from zero of the differences between the DNA methylation measured using LINE-1 and that measured in the same sample using LUMA (Fig. 3). The deviation of that difference from zero was statistically significant in all data subsets (one-sample t tests: p < 0.001, Table 2). Differences between assay results were normally distributed as indicated by non-significant Kolmogorov-Smirnov tests. A fixed bias between the results obtained with the two assays was observed, that is, the deviation of that difference from zero was consistent across the observed range of DNA methylation, i.e., the slope of a regression line through this difference did not significantly differ from zero (Table 2) indicating that the bias between both measurements was independent form the degree of methylation.
Pattern of interassay differences
Interassay differences among human blood, MCF7 cells, or SHSY5Y were large enough to render the LINE-1-LUMA difference in DNA methylation as a good predictor of tissue origin. The multimodal distribution of the differences could be described by a mixture model with M = 3 Gaussians (Fig. 4). Bayesian decision limits of −6.3 and 12.3 % DNA methylation were observed. Different data subsets were unequally represented among the Gaussians (χ 2 = 299.67, df = 8, p < 2.2 × 10−16). This provided a basis to build a decision-tree algorithm (Fig. 4) that was able to predict from the originating tissues (data subsets) in which Gaussian an interassay difference will be placed at a cross-validated accuracy of 83.6 %.
Finally, linear regression analysis was used to further characterize deviations of the results obtained from both assays from the line of identity (Fig. 2). The approach employed goodness-of-fit-based statistics to substantiate tissue-dependent deviations of the linear model from a y-intersection of zero and a slope of a value of one (Eq. 3). The goodness of fit was greatly improved when allowing separate parameter values for each data subset (Δ−2LL = −532.877, p < 0.0001). However, this was not statistically supported for every data subset, and when allowing certain subsets to share parameter value, the fit was not always worsened.
The final model indicated the following results (Table 3). Firstly, the two in vitro cell lines differed with respect to their y-intercepts, i.e., the y-intercept for the MCF7 cell line was zero while that for the SHSY5Y cell line was located at 9 % and that for the blood cells was located at 71.1 % DNA methylation. This was statistically supported by a non-significant increase in −2LL when the intersection parameter was fixed at a value of 0 for the MFC7 cell line (θ Int,Subset1; Δ−2LL = +1.178) but a significant increase when the same was done for the SHSY5Y cell line (θ Int,Subset2; Δ−2LL = +97.882; p < 0.0001) or for the blood cells (θ Int,Subset3; Δ−2LL = + 445.892; p < 0.0001). Secondly, the final linear regression model indicated that the two cell lines shared the same slope of 0.8978 as −2LL raised only non-significantly by 0.02 (p > 0.05) when setting θ Slope,Subset1 equal to θSlope,Subset2 while associating the same slope also to the blood cells worsened the fit. Thus, the results of the regression analysis indicated differences among tissue types consisting of (i) the cell lines shared the same slope but had a significantly different y-intercept and (ii) the blood sample differed from the relationship observed in the cell lines with respect to both y-intercept and slope of the linear relationship between LINE-1- and LUMA-derived measurements of global DNA methylation.
Different approaches applied to the agreement of global DNA methylation measured by LINE-1 and LUMA in three different DNA sample subsets consistently rejected the assumption of complete agreement between the two bioassays (Fig. 1). Moreover, the differences between the two assays were tissue-dependent.
The disagreement of the two assays seems biologically plausible as the two assays pursue different basic approaches not necessarily leading to the same picture of global DNA methylation. LINE-1 and LUMA differ with respect to their CpG recognition sites at the DNA. Specifically, DNA methylation occurs to 70–80 % of cytosines that locate within CpG dinucleotides . This corresponds to 3–5 % of all cytosines of the human genome . CpG dinucleotides are enriched in CpG islands, repetitive sequences, and CpG island shores  and in approximately 60 % of all gene promoters [74, 75]. The LINE-1-based assay selectively measures the methylation of CpG islands located within long interspersed nucleotide elements (LINE). These have a length of up to 6 kb and with >500,000 copies account for approximately 17–20 % of the human genome [19–21]. However, LINE-1 is unevenly distributed throughout the genome , and in addition, most of them are excluded from genomic regions containing housekeeping genes . LINE-1 is most frequently methylated in somatic tissues, where an estimated one third of DNA methylation occurs in these repetitive sequences  and particularly dense in X chromosomes [17, 25]. By including four CpG positions within the LINE-1 sequence, the pyrosequencing assay recognized 2,000,000 CpG sites when 500,000 copy numbers are estimated through the genome, i.e., roughly 7 % of the whole CpG dinucleotide contents of the human genome. However, not all LINE-1 elements are of full length; most of them are truncated and just about 10,000 LINE-1 elements contain a 5′UTR. Therefore, the effective recognized CpG dinucleotide content should be lower than 7 % .
In contrast to the LINE-1-based assay, the LUMA method measures the DNA methylation also outside repetitive elements . However, the target sequence CCGG of its restriction enzyme HpaII does also not cover all CpGs. Of the 28,000,000 CpG dinucleotides in the human genome, 4.14 % are within HpaII target sites (CCGG) located in repetitive elements and 3.90 % in unique sequences . HpaII covers 11.7 and 12.9 % of CpGs in promoter and CpG islands, respectively . C+G-rich regions of the genome have been associated with increased gene numbers [81, 82], higher amounts of CpG islands , and enhanced transcriptional activity . HpaII target sites are 15-fold enriched in CpG islands so that analysis of HpaII digested sites may over-represent potentially important regulatory sequences [25, 84]. This is a further contrast to LINE-1 sequences, which are enriched in A+T-rich gene regions  associated with fewer gene numbers and a lower transcription rate. Furthermore, from their target sequences across the genome, LINE-1 and LUMA may indeed measure different DNA methylation facets of epigenetic regulation of gene expression. This makes the consistently observed disagreement between them as biomarkers of DNA methylation biologically plausible.
The biological differences between LINE-1 and LUMA may add technical differences of the assays. Specifically, the CpG sequence targeted by LINE-1 pyrosequencing in the 5′ region tends to be deleted at unknown frequency. Approximately 2000 of the LINE-1 elements are active  that can reintegrate into the human genome results in generation of new LINE-1 sequences. Therefore, the count of the analyzed elements is unstable and may vary among different samples and individuals . Moreover, primers should amplify the region of interest regardless of its methylation status, but in practice, complete independence of the methylation is often not achievable due to a PCR bias favoring amplification of unmethylated templates . This may explain why the methylation of the completely methylated controls was not quantified as 100 % by the LINE-1-based assay (data not shown). This resembles observations with this assay in other laboratories . Although linearity of the calibration curve between non-methylated and completely methylated controls allowed for a valid recalibration of the results, the difference to the LUMA assay that quantified the completely methylated control closer to 100 % (95.7 %) is a factor contributing to the dissimilarity between the methods.
The analyses also pointed at a tissue dependency of the degree of agreement between the LINE-1- and the LUMA-based measurements of DNA methylation. This agrees with independent evidence of tissue-specific differential DNA methylation among 17 human somatic tissues . To further strengthen the observation of a cell population-dependent effect, additional measurements were performed assessing the DNA methylation level in further seven cell lines using the two assays. This additional set of samples comprised a mixture of HEK 293 (human embryonic kidney), KELLY (human neuroblastoma cells), Jurkat (human T lymphocytes), MDA-MB-468 (human mammary gland/breast cancer), HeLa (human cervical cancer), HT29 (human colon carcinoma), and THP1 (human monocytes from acute monocytic leukemia) cells. Following laboratory assays, the DNA methylation data obtained with the LINE-1 were plotted against those obtained with the LUMA method (Fig. 5). This scatterplot showed that the location of most data points was not on or close to the line of equality. Hence, differences between LINE-1- and LUMA-based measurements of global DNA methylation seem to be a consistent observation in various cell or tissue types, not restricted to the tissue types chosen for the present main analyses.
As in the present analysis, human blood was used as a frequent material in clinical epigenetic studies [26–36]; it is important to note that the heterogeneous composition of different blood cells within whole blood samples has an influence on the determined global methylation [37, 88, 89]. Therefore, associations of global methylation patterns with certain health-associated conditions, such as with inflammatory diseases, can be the result of a cell composition effect. To avoid such influences, the use of better-defined cell population should be considered. Moreover, the DNA extraction method can influence the measurement of global DNA methylation. In the present case, however, both assays were run on the same DNA samples, which reduces the probability that the differences between LINE-1- and LUMA-based readouts were caused by cell composition effects within blood samples or by an extraction bias.
Finally, while the statistical analyses and biological reasoning provided support for rejecting an agreement between the two common biomarkers of global DNA methylation, the present method comparison cannot provide a choice of the better method as a test against the gold standard as provided by an HPLC analysis. Since the gold standard HPLC specifies the DNA methylation level as the percentage of 5-mC relative to the whole cytosine amount of genomic DNA, the surrogate markers LINE-1 and LUMA specify the DNA methylation as the percentage of 5-mC in CpG dinucleotides within the recognition side. Therefore, a comparison of absolute values with the true value from HPLC analysis is not possible. Moreover, the statistically significant correlation between the results of the two assays in all data subsets (Table 2) supported the utility of both methods as global methylation markers. Indeed, an exploration of the effects sizes, estimated as Cohen’s d , produced by the various treatments to which the two cell lines were exposed supported the suitability of both methods to assess changes in DNA methylation. Specifically, a permutation approach provided a total of 72 paired comparisons between all different treatments (66 in MCF7 cells and six in SHSY5Y cells). The values of Cohen’s d  calculated for the effects resulting when the LINE-1 based assay was used, and again, when the LUMA method was used, it indicated comparable effect sizes (Fig. 6) and were significantly correlated (Spearman’s ρ: 0.79, p < 0.0001). As the sample sizes were often very small, the numerical results of this accessory analysis have, however, to be interpreted with caution.
Different approaches to the agreement of LINE-1- and LUMA-based measurements of global DNA methylation were applied to three human-derived cell types, and the assumption of complete agreement between the two bioassays were consistently rejected. Although providing partly correlated measurements of DNA methylation, interchangeability of the quantitative results obtained from the two methods was jeopardized by a consistent bias between the results. Moreover, present analyses strongly indicate a tissue specificity of the differences between the two methods.
Liu D, Diorio J, Tannenbaum B, Caldji C, Francis D, Freedman A, Sharma S, Pearson D, Plotsky PM, Meaney MJ. Maternal care, hippocampal glucocorticoid receptors, and hypothalamic-pituitary-adrenal responses to stress. Science. 1997;277(5332):1659–62.
Diorio J, Meaney MJ. Maternal programming of defensive responses through sustained effects on gene expression. J Psychiatry Neurosci. 2007;32(4):275–84.
Meaney MJ. Maternal care, gene expression, and the transmission of individual differences in stress reactivity across generations. Annu Rev Neurosci. 2001;24:1161–92.
Zhang FF, Cardarelli R, Carroll J, Zhang S, Fulda KG, Gonzalez K, Vishwanatha JK, Morabia A, Santella RM. Physical activity and global genomic DNA methylation in a cancer-free population. Epigenetics. 2011;6(3):293–9.
Christensen BC, Houseman EA, Marsit CJ, Zheng S, Wrensch MR, Wiemels JL, Nelson HH, Karagas MR, Padbury JF, Bueno R, et al. Aging and environmental exposures alter tissue-specific DNA methylation dependent upon CpG island context. PLoS Genet. 2009;5(8):e1000602.
Kamakura M. Royalactin induces queen differentiation in honeybees. Nature. 2011;473(7348):478–83.
Bollati V, Baccarelli A, Hou L, Bonzini M, Fustinoni S, Cavallo D, Byun HM, Jiang J, Marinelli B, Pesatori AC, et al. Changes in DNA methylation patterns in subjects exposed to low-dose benzene. Cancer Res. 2007;67(3):876–80.
Lötsch J, Schneider G, Reker D, Parnham MJ, Schneider P, Geisslinger G, Doehring A. Common non-epigenetic drugs as epigenetic modulators. Trends Mol Med. 2013;19(12):742–53.
Smith SS. DNA methylation in eukaryotic chromosome stability. Mol Carcinog. 1991;4(2):91–2.
Jones PA. Functions of DNA methylation: islands, start sites, gene bodies and beyond. Nat Rev Genet. 2012;13(7):484–92.
Ivanov M, Barragan I, Ingelman-Sundberg M. Epigenetic mechanisms of importance for drug treatment. Trends Pharmacol Sci. 2014;35(8):384–96.
Sasai N, Defossez PA. Many paths to one goal? The proteins that recognize methylated DNA in eukaryotes. Int J Dev Biol. 2009;53(2-3):323–34.
Ou X, Zhang Y, Xu C, Lin X, Zang Q, Zhuang T, Jiang L, von Wettstein D, Liu B. Transgenerational inheritance of modified DNA methylation patterns and enhanced tolerance induced by heavy metal stress in rice (Oryza sativa L.). PLoS One. 2012;7(9):e41143.
Rasoulpour RJ, LeBaron MJ, Ellis-Hutchings RG, Klapacz J, Gollapudi BB. Epigenetic screening in product safety assessment: are we there yet? Toxicol Mech Methods. 2011;21(4):298–311.
Ohka F, Natsume A, Motomura K, Kishida Y, Kondo Y, Abe T, Nakasu Y, Namba H, Wakai K, Fukui T, et al. The global DNA methylation surrogate LINE-1 methylation is correlated with MGMT promoter methylation and is a better prognostic factor for glioma. PLoS One. 2011;6(8):e23332<.
Wolffe AP, Matzke MA. Epigenetics: regulation through repression. Science. 1999;286(5439):481–6.
Nelson HH, Marsit CJ, Kelsey KT. Global methylation in exposure biology and translational medical science. Environ Health Perspect. 2011;119(11):1528–33.
Florea AM. DNA methylation pyrosequencing assay is applicable for the assessment of epigenetic active environmental or clinical relevant chemicals. BioMed Res Int. 2013;2013:486072.
Kazazian Jr HH, Goodier JL. LINE drive. Retrotransposition and genome instability. Cell. 2002;110(3):277–80.
Ehrlich M. DNA methylation in cancer: too much, but also too little. Oncogene. 2002;21(35):5400–13.
Deininger PL, Moran JV, Batzer MA, Kazazian Jr HH. Mobile elements and mammalian genome evolution. Curr Opin Genet Dev. 2003;13(6):651–8.
Mohsen K, Johansson S, Ekström T. Using LUMA: a luminometric-based assay for global DNA-methylation. Epigenetics. 2006;1(1):46–9.
Weisenberger DJ, Campan M, Long TI, Kim M, Woods C, Fiala E, Ehrlich M, Laird PW. Analysis of repetitive element DNA methylation by MethyLight. Nucleic Acids Res. 2005;33(21):6823–36.
Lisanti S, Omar WA, Tomaszewski B, De Prins S, Jacobs G, Koppen G, Mathers JC, Langie SA. Comparison of methods for quantification of global DNA methylation in human cells and tissues. PLoS One. 2013;8(11):e79044.
Fazzari MJ, Greally JM. Epigenomics: beyond CpG islands. Nat Rev Genet. 2004;5(6):446–55.
Di Francesco A, Arosio B, Falconi A, Micioni Di Bonaventura MV, Karimi M, Mari D, Casati M, Maccarrone M, D’Addario C. Global changes in DNA methylation in Alzheimer’s disease peripheral blood mononuclear cells. Brain Behav Immun. 2015;45:139–44.
Xu X, Gammon MD, Hernandez-Vargas H, Herceg Z, Wetmur JG, Teitelbaum SL, Bradshaw PT, Neugut AI, Santella RM, Chen J. DNA methylation in peripheral blood measured by LUMA is associated with breast cancer in a population-based study. FASEB J. 2012;26(6):2657–66.
Akers SN, Moysich K, Zhang W, Collamat Lai G, Miller A, Lele S, Odunsi K, Karpf AR. LINE1 and Alu repetitive element DNA methylation in tumors and white blood cells from epithelial ovarian cancer patients. Gynecol Oncol. 2014;132(2):462–7.
Andreotti G, Karami S, Pfeiffer RM, Hurwitz L, Liao LM, Weinstein SJ, Albanes D, Virtamo J, Silverman DT, Rothman N, et al. LINE1 methylation levels associated with increased bladder cancer risk in pre-diagnostic blood DNA among US (PLCO) and European (ATBC) cohort study participants. Epigenetics : official journal of the DNA Methylation Society. 2014;9(3):404–15.
Marques-Rocha JL, Milagro FI, Mansego ML, Mourao DM, Martinez JA, Bressan J. LINE-1 methylation is positively associated with healthier lifestyle but inversely related to body fat mass in healthy young individuals. Epigenetics : official journal of the DNA Methylation Society. 2016;11(1):49–60.
Bozovic IB, Stankovic A, Zivkovic M, Vranekovic J, Kapovic M, Brajenovic-Milic B. Altered LINE-1 methylation in mothers of children with down syndrome. PLoS One. 2015;10(5):e0127423.
Martin-Nunez GM, Rubio-Martin E, Cabrera-Mulero R, Rojo-Martinez G, Olveira G, Valdes S, Soriguer F, Castano L, Morcillo S. Type 2 diabetes mellitus in relation to global LINE-1 DNA methylation in peripheral blood: a cohort study. Epigenetics : official journal of the DNA Methylation Society. 2014;9(10):1322–8.
McCullough LE, Chen J, Cho YH, Khankari NK, Bradshaw PT, White AJ, et al. DNA methylation modifies the association between obesity and survival after breast cancer diagnosis. Breast Cancer Res Treat. 2016.
Ding R, Jin Y, Liu X, Zhu Z, Zhang Y, Wang T, Xu Y. Characteristics of DNA methylation changes induced by traffic-related air pollution. Mutat Res Genet Toxicol Environ Mutagen. 2016;796:46–53.
Zheng P, Yang T, Ju L, Jiang B, Lou Y. Epigenetics in Legg-Calve-Perthes disease: a study of global DNA methylation. J Int Med Res. 2015;43(6):758–64.
Brown TA, Lee JW, Holian A, Porter V, Fredriksen H, Kim M, Cho YH. Alterations in DNA methylation corresponding with lung inflammation and as a biomarker for disease development after MWCNT exposure. Nanotoxicology. 2016;10(4):453–61.
Wu HC, Delgado-Cruzata L, Flom JD, Kappil M, Ferris JS, Liao Y, Santella RM, Terry MB. Global methylation profiles in DNA from different blood cell types. Epigenetics : official journal of the DNA Methylation Society. 2011;6(1):76–85.
Doehring A, Oertel BG, Sittl R, Lötsch J. Chronic opioid use is associated with increased DNA methylation correlating with increased clinical pain. Pain. 2013;154(1):15–23.
Knothe C, Doehring A, Ultsch A, Lötsch J. Methadone induces hypermethylation of human DNA. Epigenomics. 2015.
Doehring A, Freynhagen R, Griessinger N, Zimmermann M, Sittl R, Hentig N, Geisslinger G, Lötsch J. Cross-sectional assessment of the consequences of a GTP cyclohydrolase 1 haplotype for specialized tertiary outpatient pain care. Clin J Pain. 2009;25(9):781–5.
Doehring A, Hentig N, Graff J, Salamat S, Schmidt M, Geisslinger G, Harder S, Lötsch J. Genetic variants altering dopamine D2 receptor expression or function modulate the risk of opiate addiction and the dosage requirements of methadone substitution. Pharmacogenet Genomics. 2009;19(6):407–14.
Lötsch J, Freynhagen R, von Hentig N, Griessinger N, Zimmermann M, Sittl R, Geisslinger G. Higher pain scores, similar opioid doses and side effects associated with antipyretic analgesics in specialised tertiary pain care. Inflamm Res. 2010;59(11):989–95.
Striebel W. Therapie chronischer Schmerzen. Ein praktischer Leitfaden. Stuttgart: Schattauer; 2002.
Twycross R, Wilcock A. Palliative care formulary. In., 3 edn. Oxford: palliativebooks.com; 2007: 269.
Stresemann C, Brueckner B, Musch T, Stopper H, Lyko F. Functional diversity of DNA methyltransferase inhibitors in human cancer cell lines. Cancer Res. 2006;66(5):2794–800.
Momparler RL, Momparler LF, Samson J. Comparison of the antileukemic activity of 5-AZA-2′-deoxycytidine, 1-beta-D-arabinofuranosylcytosine and 5-azacytidine against L1210 leukemia. Leuk Res. 1984;8(6):1043–9.
Brueckner B, Garcia Boy R, Siedlecki P, Musch T, Kliem HC, Zielenkiewicz P, Suhai S, Wiessler M, Lyko F. Epigenetic reactivation of tumor suppressor genes by a novel small-molecule inhibitor of human DNA methyltransferases. Cancer Res. 2005;65(14):6305–11.
Hollenbach PW, Nguyen AN, Brady H, Williams M, Ning Y, Richard N, Krushel L, Aukerman SL, Heise C, MacBeth KJ. A comparison of azacitidine and decitabine activities in acute myeloid leukemia cell lines. PLoS ONE. 2010;5(2):e9001.
Ham M-S, Lee J-K, Kim K-C. S-adenosyl methionine specifically protects the anticancer effect of 5-FU via DNMTs expression in human A549 lung cancer cells. Mol Clin Oncol. 2013;1(2):373–8.
Chik F, Machnes Z, Szyf M. Synergistic anti-breast cancer effect of a combined treatment with the methyl donor S-adenosyl methionine and the DNA methylation inhibitor 5-aza-2′-deoxycytidine. Carcinogenesis. 2014;35(1):138–44.
Zhang X, Li H, Qiu Q, Qi Y, Huang D, Zhang Y. 2,4-Dichlorophenol induces global DNA hypermethylation through the increase of S-adenosylmethionine and the upregulation of DNMTs mRNA in the liver of goldfish Carassius auratus. Comp Biochem Physiol Toxicol Pharmacol. 2014;160:54–9.
Kile ML, Baccarelli A, Tarantini L, Hoffman E, Wright RO, Christiani DC. Correlation of global and gene-specific DNA methylation in maternal-infant pairs. PLoS One. 2010;5(10):e13730.
Karimi M, Johansson S, Stach D, Corcoran M, Grander D, Schalling M, Bakalkin G, Lyko F, Larsson C, Ekstrom TJ. LUMA (LUminometric Methylation Assay)—a high throughput method to the analysis of genomic DNA methylation. Exp Cell Res. 2006;312(11):1989–95.
Keller M, Kralisch S, Rohde K, Schleinitz D, Dietrich A, Schon MR, Gartner D, Lohmann T, Dressler M, Tonjes A et al. Global DNA methylation levels in human adipose tissue are related to fat distribution and glucose homeostasis. Diabetologia. 2014;57(11):2374–83.
Magari RT. Statistics for laboratory method comparison studies. Biopharm Int. 2002;15(1):28–32.
Wilcoxon F. Individual comparisons by ranking methods. Biometrics. 1945;1:80–3.
Hochberg Y. A sharper Bonferroni procedure for multiple tests of significance. Biometrika. 1988;75(4):800–2.
Spearman C. The proof and measurement of association between two things. Am J Psychol. 1904;15:72–101.
Efron B, Tibshirani RJ. An introduction to the bootstrap. San Francisco: Chapman and Hall; 1995.
Bland JM, Altman DG. Applying the right statistics: analyses of measurement studies. Ultrasound Obstet Gynecol. 2003;22(1):85–93.
Bland JM, Altman DG. Statistical methods for assessing agreement between two methods of clinical measurement. Lancet. 1986;1:307–10.
Smirnov N. Table for estimating the goodness of fit of empirical distributions. Ann Math Statist. 1948;19(2):279-281.
Ultsch A. Pareto density estimation: a density estimation for knowledge discovery. In: Innovations in classification, data science, and information systems—proceedings 27th annual conference of the German classification society (GfKL): 2003. Berlin: Springer; 2003.
Ultsch A, Thrun MC, Hansen-Goos O, Lötsch J. Identification of molecular fingerprints in human heat pain thresholds by use of an interactive mixture model R toolbox (AdaptGauss). Int J Mol Sci. 2015;16(10):25897–911.
Stork DG, Duda RO, Hart PE, Stork D. Pattern classification. IO/Il7 2001.
McGrayne SB. The theory that would not die: how Bayes’ rule cracked the enigma code, hunted down Russian submarines and emerged triumphant from two centuries of controversy. New Haven: Yale University Press; 2011.
Breimann L, Friedman JH, Olshen RA, Stone CJ. Classification and regression trees. Boca Raton: Chapman and Hall; 1993.
Hill T, Lewicki P. STATISTICS: methods and applications. Tulsa: StatSoft; 2007.
Boeckmann AJ, Sheiner LB, Beal SL. NONMEM user’s guide. San Francisco: University of California San Francisco; 1994.
Yafune A, Ishiguro M. Bootstrap approach for constructing confidence intervals for population pharmacokinetic parameters. I: A use of bootstrap standard error. Stat Med. 1999;18(5):581–99.
Wigler MH. The inheritance of methylation patterns in vertebrates. Cell. 1981;24(2):285–6.
Clark SJ, Harrison J, Frommer M. CpNpG methylation in mammalian cells. Nat Genet. 1995;10(1):20–7.
Weber M, Hellmann I, Stadler MB, Ramos L, Paabo S, Rebhan M, Schubeler D. Distribution, silencing potential and evolutionary impact of promoter DNA methylation in the human genome. Nat Genet. 2007;39(4):457–66.
Ioshikhes IP, Zhang MQ. Large-scale human promoter mapping using CpG islands. Nat Genet. 2000;26(1):61–3.
Fahy J, Jeltsch A, Arimondo PB. DNA methyltransferase inhibitors in cancer: a chemical and therapeutic patent overview and selected clinical studies. Expert Opin Ther Pat. 2012;22(12):1427–42.
Aporntewan C, Phokaew C, Piriyapongsa J, Ngamphiw C, Ittiwut C, Tongsima S, Mutirangura A. Hypomethylation of intragenic LINE-1 represses transcription in cancer cells through AGO2. PLoS One. 2011;6(3):e17934.
Eller CD, Regelson M, Merriman B, Nelson S, Horvath S, Marahrens Y. Repetitive sequence environment distinguishes housekeeping genes. Gene. 2007;390(1-2):153–65.
Yang AS, Estecio MR, Doshi K, Kondo Y, Tajara EH, Issa JP. A simple method for estimating global DNA methylation using bisulfite PCR of repetitive DNA elements. Nucleic Acids Res. 2004;32(3):e38.
Penzkofer T, Dandekar T, Zemojtel T. L1Base: from functional annotation to prediction of active LINE-1 elements. Nucleic Acids Res. 2005;33(Database issue):D498–500.
Wong JMCaNC. Epigenetics: a reference manual: Caister Academic Press; 2011.
Federico C, Andreozzi L, Saccone S, Bernardi G. Gene density in the Giemsa bands of human chromosomes. Chromosome Res. 2000;8(8):737–46.
Arhondakis S, Auletta F, Torelli G, D’Onofrio G. Base composition and expression level of human genes. Gene. 2004;325:165–9.
Jabbari K, Bernardi G. CpG doublets, CpG islands and Alu repeats in long human DNA sequences from different isochore families. Gene. 1998;224(1-2):123–7.
Bjornsson HT, Sigurdsson MI, Fallin MD, Irizarry RA, Aspelund T, Cui H, Yu W, Rongione MA, Ekstrom TJ, Harris TB et al. Intra-individual change over time in DNA methylation with familial clustering. Jama. 2008;299(24):2877–83.
Pavlicek A, Jabbari K, Paces J, Paces V, Hejnar JV, Bernardi G. Similar integration but different stability of Alus and LINEs in the human genome. Gene. 2001;276(1-2):39–45.
Wojdacz TK, Hansen LL, Dobrovic A. A new approach to primer design for the control of PCR bias in methylation studies. BMC Res Notes. 2008;1:54.
Lokk K, Modhukur V, Rajashekar B, Martens K, Magi R, Kolde R, Koltsina M, Nilsson TK, Vilo J, Salumets A, et al. DNA methylome profiling of human tissues identifies global and tissue-specific methylation patterns. Genome Biol. 2014;15(4):r54.
Houseman EA, Kim S, Kelsey KT, Wiencke JK. DNA methylation in whole blood: uses and challenges. Curr Environ Health Rep. 2015;2(2):145–54.
Delgado-Cruzata L, Vin-Raviv N, Tehranifar P, Flom J, Reynolds D, Gonzalez K, Santella RM, Terry MB. Correlations in global DNA methylation measures in peripheral blood mononuclear cells and granulocytes. Epigenetics : official journal of the DNA Methylation Society. 2014;9(11):1504–10.
Cohen J. A power primer. Psych Bull. 1992;112:155–9.
Kampstra P. Beanplot: a boxplot alternative for visual comparison of distributions. J Stat Software. 2008;28:Code Snippet 1.
Williams GJ. Data mining with rattle and R: the art of excavating data for knowledge discovery. Springer; 2011.
The authors thank S. Grösch, Goethe University, Frankfurt am Main, Germany, for kindly providing the MFC7 cell line and C. Feinweber for kindly providing the extracted DNA from further cell lines (MDA-MB-468; HeLa and HT29 cells).
The work has been supported by the Else Kröner-Fresenius Foundation (EKFS), Research Training Group Translational Research Innovation—Pharma (TRIP, GG, JL), and the Landesoffensive zur Entwicklung wissenschaftlich-ökonomischer Exzellenz (LOEWE), LOEWE-Zentrum für Translationale Medizin und Pharmakologie (GG, JL). Additional support of the analytical environment was gained from the European Union Seventh Framework Programme (FP7/2007–2013) under grant agreement no. 602919 (“GLORIA”, JL,) The funders had no role in the method design, data selection and analysis, decision to publish, or preparation of the manuscript.
CK, JL, and AD conceived and designed the experiments. CK, HS, and ER performed the experiments. JL and AU analyzed the data. JL and CK wrote the paper. GG contributed regencies/materials. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.