Comparison of novel and existing methods for detecting differentially methylated regions

Epigenome-wide association study of short-term temperature fluctuations based on within-sibship analyses in Australian females

BACKGROUND

Temperature fluctuations can affect human health independent of the effect of mean temperature. However, no study has evaluated whether short-term temperature fluctuations could affect DNA methylation.

METHODS

Peripheral blood DNA methylation for 479 female siblings of 130 families were analysed. Gridded daily temperatures data were obtained, linked to each participant's home address, and used to calculate nine different metrics of short-term temperature fluctuations: temperature variabilities (TVs) within the day of blood draw and preceding one to seven days (TV 0-1 to TV 0-7), diurnal temperature range (DTR), and temperature change between neighbouring days (TCN). Within-sibship design was used to perform epigenome-wide association analyses, adjusting for daily mean temperatures, and other important covariates (e.g., smoking, alcohol use, cell-type proportions). Differentially methylated regions (DMRs) were further identified. Multiple-testing comparisons with a significant threshold of 0.01 for cytosine-guanine dinucleotides (CpGs) and 0.05 for DMRs were applied.

RESULTS

Among 479 participants (mean age ± SD, 56.4 ± 7.9 years), we identified significant changes in methylation levels in 14 CpGs and 70 DMRs associated with temperature fluctuations. Almost all identified CpGs were associated with exposure to temperature fluctuations within three days. Differentially methylated signals were mapped to 68 genes that were linked to human diseases such as cancer (e.g., colorectal carcinoma, breast carcinoma, and metastatic neoplasms) and mental disorder (e.g., schizophrenia, mental depression, and bipolar disorder). The top three most significantly enriched gene ontology terms were Response to bacterium (TV 0-3), followed by Hydrolase activity, acting on ester bonds (TCN), and Oxidoreductase activity (TV 0-3).

CONCLUSIONS

Short-term temperature fluctuations were associated with differentially methylated signals across the human genome, which provides evidence on the potential biological mechanisms underlying the health impact of temperature fluctuations. Future studies are needed to further clarify the roles of DNA methylation in diseases associated with temperature fluctuations.

An evaluation of the genome-wide false positive rates of common methods for identifying differentially methylated regions using illumina methylation arrays

Differentially methylated regions (DMRs) are genomic regions with specific methylation patterns across multiple loci that are associated with a phenotype. We examined the genome-wide false positive (GFP) rates of five widely used DMR methods: comb-p, Bumphunter, DMRcate, mCSEA and coMethDMR using both Illumina HumanMethylation450 (450 K) and MethylationEPIC (EPIC) data and simulated continuous and dichotomous null phenotypes (i.e., generated independently of methylation data). coMethDMR provided well-controlled GFP rates (~5%) except when analysing skewed continuous phenotypes. DMRcate generally had well-controlled GFP rates when applied to 450 K data except for the skewed continuous phenotype and EPIC data only for the normally distributed continuous phenotype. GFP rates for mCSEA were at least 0.096 and comb-p yielded GFP rates above 0.34. Bumphunter had high GFP rates of at least 0.35 across conditions, reaching as high as 0.95. Analysis of the performance of these methods in specific regions of the genome found that regions with higher correlation across loci had higher regional false positive rates on average across methods. Based on the false positive rates, coMethDMR is the most recommended analysis method, and DMRcate had acceptable performance when analysing 450 K data. However, as both could display higher levels of FPs for skewed continuous distributions, a normalizing transformation of skewed continuous phenotypes is suggested. This study highlights the importance of genome-wide simulations when evaluating the performance of DMR-analysis methods.

Obesity-Associated Differentially Methylated Regions in Colon Cancer

Obesity with adiposity is a common disorder in modern days, influenced by environmental factors such as eating and lifestyle habits and affecting the epigenetics of adipose-based gene regulations and metabolic pathways in colorectal cancer (CRC). We compared epigenetic changes of differentially methylated regions (DMR) of genes in colon tissues of 225 colon cancer cases (154 non-obese and 71 obese) and 15 healthy non-obese controls by accessing The Cancer Genome Atlas (TCGA) data. We applied machine-learning-based analytics including generalized regression (GR) as a confirmatory validation model to identify the factors that could contribute to DMRs impacting colon cancer to enhance prediction accuracy. We found that age was a significant predictor in obese cancer patients, both alone (p = 0.003) and interacting with hypomethylated DMRs of ZBTB46, a tumor suppressor gene (p = 0.008). DMRs of three additional genes: HIST1H3I (p = 0.001), an oncogene with a hypomethylated DMR in the promoter region; SRGAP2C (p = 0.006), a tumor suppressor gene with a hypermethylated DMR in the promoter region; and NFATC4 (p = 0.006), an adipocyte differentiating oncogene with a hypermethylated DMR in an intron region, are also significant predictors of cancer in obese patients, independent of age. The genes affected by these DMR could be potential novel biomarkers of colon cancer in obese patients for cancer prevention and progression.

A Mixture of Endocrine Disrupting Chemicals Associated with Lower Birth Weight in Children Induces Adipogenesis and DNA Methylation Changes in Human Mesenchymal Stem Cells

Endocrine Disrupting Chemicals (EDCs) are man-made compounds that alter functions of the endocrine system. Environmental mixtures of EDCs might have adverse effects on human health, even though their individual concentrations are below regulatory levels of concerns. However, studies identifying and experimentally testing adverse effects of real-life mixtures are scarce. In this study, we aimed at evaluating an epidemiologically identified EDC mixture in an experimental setting to delineate its cellular and epigenetic effects. The mixture was established using data from the Swedish Environmental Longitudinal Mother and child Asthma and allergy (SELMA) study where it was associated with lower birth weight, an early marker for prenatal metabolic programming. This mixture was then tested for its ability to change metabolic programming of human mesenchymal stem cells. In these cells, we assessed if the mixture induced adipogenesis and genome-wide DNA methylation changes. The mixture increased lipid droplet accumulation already at concentrations corresponding to levels measured in the pregnant women of the SELMA study. Furthermore, we identified differentially methylated regions in genes important for adipogenesis and thermogenesis. This study shows that a mixture reflecting human real-life exposure can induce molecular and cellular changes during development that could underlie adverse outcomes.

Identification of two methylated fragments of an SDC2 CpG island using a sliding window technique for early detection of colorectal cancer

Colorectal cancer (CRC) is one of the most common cancer types globally with a 5‐year survival rate of < 50% in China. Aberrant DNA methylation is one of the hallmarks of tumor initiation, progression, and metastasis. Here, we investigated the clinical performance of two differentially methylated regions (DMRs) in SDC2 CpG islands for the detection of CRC. A sliding window technique was used to identify the DMRs, and methylation‐specific PCR assay was used to assess the DMRs in 198 CRC samples and 54 normal controls. Two DMRs (DMR2 and DMR5) were identified using The Cancer Genome Atlas (TCGA) data, and the hypermethylation of DMR2 and DMR5 was detected in 90.91% (180/198) and 89.90% (178/198) of CRC samples, respectively. When combining DMR2 and DMR5, the sensitivity for CRC detection was 94.4% higher than that of DMR2 or DMR5 alone. Based on the above results, we propose using DMR2 and DMR5 as a sensitive biomarker to detect CRC.

Detecting differentially methylated regions with multiple distinct associations

Aim: We evaluated five methods for detecting differentially methylated regions (DMRs): DMRcate, comb-p, seqlm, GlobalP and dmrff. Materials & methods: We used a simulation study and real data analysis to evaluate performance. Additionally, we evaluated the use of an ancestry-matched reference cohort to estimate correlations between CpG sites in cord blood. Results: Several methods had inflated Type I error, which increased at more stringent significant levels. In power simulations with 1-2 causal CpG sites with the same direction of effect, dmrff was consistently among the most powerful methods. Conclusion: This study illustrates the need for more thorough simulation studies when evaluating novel methods. More work must be done to develop methods with well-controlled Type I error that do not require individual-level data.

Fibrate pharmacogenomics: expanding past the genome

Fibrates are a medication class prescribed for decades as 'broad-spectrum' lipid-modifying agents used to lower blood triglyceride levels and raise high-density lipoprotein cholesterol levels. Such lipid changes are associated with a decrease in cardiovascular disease, and fibrates are commonly used to reduce risk of dangerous cardiovascular outcomes. As with most drugs, it is well established that response to fibrate treatment is variable, and this variation is heritable. This has motivated the investigation of pharmacogenomic determinants of response, and multiple studies have discovered a number of genes associated with fibrate response. Similar to other complex traits, the interrogation of single nucleotide polymorphisms using candidate gene or genome-wide approaches has not revealed a substantial portion of response variation. However, recent innovations in technological platforms and advances in statistical methodologies are revolutionizing the use and integration of other 'omes' in pharmacogenomics studies. Here, we detail successes, challenges, and recent advances in fibrate pharmacogenomics.

Application of novel and existing methods to identify genes with evidence of epigenetic association: results from GAW20

BACKGROUND

The rise in popularity and accessibility of DNA methylation data to evaluate epigenetic associations with disease has led to numerous methodological questions. As part of GAW20, our working group of 8 research groups focused on gene searching methods.

RESULTS

Although the methods were varied, we identified 3 main themes within our group. First, many groups tackled the question of how best to use pedigree information in downstream analyses, finding that (a) the use of kinship matrices is common practice, (b) ascertainment corrections may be necessary, and (c) pedigree information may be useful for identifying parent-of-origin effects. Second, many groups also considered multimarker versus single-marker tests. Multimarker tests had modestly improved power versus single-marker methods on simulated data, and on real data identified additional associations that were not identified with single-marker methods, including identification of a gene with a strong biological interpretation. Finally, some of the groups explored methods to combine single-nucleotide polymorphism (SNP) and DNA methylation into a single association analysis.

CONCLUSIONS

A causal inference method showed promise at discovering new mechanisms of SNP activity; gene-based methods of summarizing SNP and DNA methylation data also showed promise. Even though numerous questions still remain in the analysis of DNA methylation data, our discussions at GAW20 suggest some emerging best practices.