Disagreement between two common biomarkers of global DNA methylation
© Knothe et al. 2016
Received: 21 March 2016
Accepted: 10 May 2016
Published: 23 May 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
Samples, conditions, and concentrations used for the assessment of methylation levels by means of LINE-1 and LUMA, of which the cell lines had been obtained previously 
Tissue (cell type)
Age (mean ± SD) (years)
n (total number of replicates)
Human breast cancer cell line
5-Aza-CdR + SAM
0.3 + 50/100 (μM)
5-Aza-CdR + DCP
0.3 (μM) + 10 (mg/l)
5-Aza-CdR + Methadone
0.3 + 10/25/50/75/100 (μM)
5-Aza-CdR + SAM + DCP
0.3 (μM) + 100 (μM) + 10 (mg/l)
RG108 + DCP
100 (μM) + 10 (mg/l)
SAM + DCP
100 (μM) + 10 (mg/l)
Human neuronal cell line
3 days/7 days
3 days/7 days
3 days/7 days
5-Aza-CdR + SAM
0.3 + 30 (μM)
5-Aza-CdR + Methadone
0.3 + 10/100 (μM)
39.57 ± 7.02
Pain patients, opioid treated
47.83 ± 7.36
52.69 mg ± 22.11 OMEa
Pain patients non-opioid treated
45.68 ± 11.63
Descriptive and inference statistics of the two assays
Method of 95 % limits of agreement between measurements (“Bland-Altman”)
Data subset (cell type)
LINE-1 (mean ± SD and (range)) (%)
LUMA (mean ± SD and (range)) (%)
Wilcoxon signed rank test
Spearman Correlation (95 % CI)
Mean difference (fixed bias) (95 % CI)
Significance of mean deviation from 0
KS-test of normal distribution of differences
95 % confidence limits of agreement
Slope (proportional bias) (CI of estimate)
Significance of deviation of slope from 0
60.96 ± 11.99 (34, 76.18)
67.8 ± 11.99 (34.36, 78.34)
V = 2476, p = 6.22 × 10−10
ρ = 0.58 (0.38, 0.74), p = 5.53 × 10−8
−6.8 (−8.3, −5.3)
1.93 × 10−13
D = 0.1183, p = 0.24
−0.00046 (−0.13, 0.13)
59.01 ± 11.54 (31.91, 71.65)
55.66 ± 11.99 (24.89, 67.65)
V = 74, p = 4.66 × 10−5
ρ = 0.8 (0.56, 0.92), p = 1.31 × 10−8
3.3 (1.9, 4.8)
3.16 × 10−5
D = 0.068241, p = 0.99
−0.04 (−0.16, 0.08)
92.4 ± 2.75 (77.9, 98.8)
72.3 ± 3.1 (61.9, 78.5)
V = 0, p <2.2 × 10−16
ρ = 0.35 (0.18, 0.51), p = 3.25 × 10−5
20.1 (19.5, 20.7)
3.41 × 10−103
D = 0.086144, p = 0.29
−0.2 (−0.44, 0.045)
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
θ 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.
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.
Parameters and estimated values of the final linear model of the data-subset specific relation of LINE-1 and LUMA assay-based measurements of global DNA methylation. The full model was given as MethylationLINE1 = (Intersection + θ Int, Subset1.. Subset3) + (Slope × θ Slope, Subset1.. Subset3) × MethylationLUMA + ε, where Intersection and Slope are structural parameters of the linear model denoted during the fitting as θ 1 and θ 1, respectively, θ Int,Subset1.. Subset3 and θ Slope, Subset1.. Subset3 are data subset specific modulators of the structural parameter values, denoted during the fitting as θ 3..12, and η 1 and ε accounts for the additive error in the fit of the percent methylation data acquired by means of two different assays. The final model was the result of the model building favoring the best but sparsest model based on goodness-of-fit statistics
Value (and % SEE)
95 % bootstrap CI
Y-intersection (% methylation) = θ 1
θ 3 = θ Int,Subset1
θ 4 = θ Int,Subset2
θ 5 = θ Int,Subset3
Slope = θ 2
θ 6 = θ Slope,Subset1
θ 6 = θ7 = 2.68 (0.7)
θ 7 = θ Slope,Subset2
SlopeSubset1,2 = 2.68 · 0.38 = 0.9
θ 8 = θ Slope,Subset3
SlopeSubset3 = 0.88 · 0.38 = 0.29
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.
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.
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.
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.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- 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.View ArticlePubMedGoogle Scholar
- Diorio J, Meaney MJ. Maternal programming of defensive responses through sustained effects on gene expression. J Psychiatry Neurosci. 2007;32(4):275–84.PubMedPubMed CentralGoogle Scholar
- Meaney MJ. Maternal care, gene expression, and the transmission of individual differences in stress reactivity across generations. Annu Rev Neurosci. 2001;24:1161–92.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- Kamakura M. Royalactin induces queen differentiation in honeybees. Nature. 2011;473(7348):478–83.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Smith SS. DNA methylation in eukaryotic chromosome stability. Mol Carcinog. 1991;4(2):91–2.View ArticlePubMedGoogle Scholar
- Jones PA. Functions of DNA methylation: islands, start sites, gene bodies and beyond. Nat Rev Genet. 2012;13(7):484–92.View ArticlePubMedGoogle Scholar
- Ivanov M, Barragan I, Ingelman-Sundberg M. Epigenetic mechanisms of importance for drug treatment. Trends Pharmacol Sci. 2014;35(8):384–96.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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<.View ArticlePubMedPubMed CentralGoogle Scholar
- Wolffe AP, Matzke MA. Epigenetics: regulation through repression. Science. 1999;286(5439):481–6.View ArticlePubMedGoogle Scholar
- Nelson HH, Marsit CJ, Kelsey KT. Global methylation in exposure biology and translational medical science. Environ Health Perspect. 2011;119(11):1528–33.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- Kazazian Jr HH, Goodier JL. LINE drive. Retrotransposition and genome instability. Cell. 2002;110(3):277–80.View ArticlePubMedGoogle Scholar
- Ehrlich M. DNA methylation in cancer: too much, but also too little. Oncogene. 2002;21(35):5400–13.View ArticlePubMedGoogle Scholar
- Deininger PL, Moran JV, Batzer MA, Kazazian Jr HH. Mobile elements and mammalian genome evolution. Curr Opin Genet Dev. 2003;13(6):651–8.View ArticlePubMedGoogle Scholar
- Mohsen K, Johansson S, Ekström T. Using LUMA: a luminometric-based assay for global DNA-methylation. Epigenetics. 2006;1(1):46–9.View ArticleGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- Fazzari MJ, Greally JM. Epigenomics: beyond CpG islands. Nat Rev Genet. 2004;5(6):446–55.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.Google Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Knothe C, Doehring A, Ultsch A, Lötsch J. Methadone induces hypermethylation of human DNA. Epigenomics. 2015.Google Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Striebel W. Therapie chronischer Schmerzen. Ein praktischer Leitfaden. Stuttgart: Schattauer; 2002.Google Scholar
- Twycross R, Wilcock A. Palliative care formulary. In., 3 edn. Oxford: palliativebooks.com; 2007: 269.Google Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.PubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Magari RT. Statistics for laboratory method comparison studies. Biopharm Int. 2002;15(1):28–32.Google Scholar
- Wilcoxon F. Individual comparisons by ranking methods. Biometrics. 1945;1:80–3.View ArticleGoogle Scholar
- Hochberg Y. A sharper Bonferroni procedure for multiple tests of significance. Biometrika. 1988;75(4):800–2.View ArticleGoogle Scholar
- Spearman C. The proof and measurement of association between two things. Am J Psychol. 1904;15:72–101.View ArticleGoogle Scholar
- Efron B, Tibshirani RJ. An introduction to the bootstrap. San Francisco: Chapman and Hall; 1995.Google Scholar
- Bland JM, Altman DG. Applying the right statistics: analyses of measurement studies. Ultrasound Obstet Gynecol. 2003;22(1):85–93.View ArticlePubMedGoogle Scholar
- Bland JM, Altman DG. Statistical methods for assessing agreement between two methods of clinical measurement. Lancet. 1986;1:307–10.View ArticlePubMedGoogle Scholar
- Smirnov N. Table for estimating the goodness of fit of empirical distributions. Ann Math Statist. 1948;19(2):279-281.Google Scholar
- 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.Google Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- Stork DG, Duda RO, Hart PE, Stork D. Pattern classification. IO/Il7 2001.Google Scholar
- 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.Google Scholar
- Breimann L, Friedman JH, Olshen RA, Stone CJ. Classification and regression trees. Boca Raton: Chapman and Hall; 1993.Google Scholar
- Hill T, Lewicki P. STATISTICS: methods and applications. Tulsa: StatSoft; 2007.Google Scholar
- Boeckmann AJ, Sheiner LB, Beal SL. NONMEM user’s guide. San Francisco: University of California San Francisco; 1994.Google Scholar
- 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.View ArticlePubMedGoogle Scholar
- Wigler MH. The inheritance of methylation patterns in vertebrates. Cell. 1981;24(2):285–6.View ArticlePubMedGoogle Scholar
- Clark SJ, Harrison J, Frommer M. CpNpG methylation in mammalian cells. Nat Genet. 1995;10(1):20–7.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Ioshikhes IP, Zhang MQ. Large-scale human promoter mapping using CpG islands. Nat Genet. 2000;26(1):61–3.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- Wong JMCaNC. Epigenetics: a reference manual: Caister Academic Press; 2011.Google Scholar
- Federico C, Andreozzi L, Saccone S, Bernardi G. Gene density in the Giemsa bands of human chromosomes. Chromosome Res. 2000;8(8):737–46.View ArticlePubMedGoogle Scholar
- Arhondakis S, Auletta F, Torelli G, D’Onofrio G. Base composition and expression level of human genes. Gene. 2004;325:165–9.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticleGoogle Scholar
- Cohen J. A power primer. Psych Bull. 1992;112:155–9.View ArticleGoogle Scholar
- Kampstra P. Beanplot: a boxplot alternative for visual comparison of distributions. J Stat Software. 2008;28:Code Snippet 1.View ArticleGoogle Scholar
- Williams GJ. Data mining with rattle and R: the art of excavating data for knowledge discovery. Springer; 2011.Google Scholar