Critical evaluation of linear regression models for cell-subtype specific methylation signal from mixed blood cell DNA

PLoS One. 2018 Dec 20;13(12):e0208915. doi: 10.1371/journal.pone.0208915. eCollection 2018.

Abstract

Epigenome-wide association studies seek to identify DNA methylation sites associated with clinical outcomes. Difference in observed methylation between specific cell-subtypes is often of interest; however, available samples often comprise a mixture of cells. To date, cell-subtype estimates have been obtained from mixed-cell DNA data using linear regression models, but the accuracy of such estimates has not been critically assessed. We evaluated linear regression performance for cell-subtype specific methylation estimation using a 450K methylation array dataset of both mixed-cell and cell-subtype sorted samples from six healthy males. CpGs associated with each cell-subtype were first identified using t-tests between groups of cell-subtype sorted samples. Subsequent reduced panels of reliably accurate CpGs were identified from mixed-cell samples using an accuracy heuristic (D). Performance was assessed by comparing cell-subtype specific estimates from mixed-cells with corresponding cell-sorted mean using the mean absolute error (MAE) and the Coefficient of Determination (R2). At the cell-subtype level, methylation levels at 3272 CpGs could be estimated to within a MAE of 5% of the expected value. The cell-subtypes with the highest accuracy were CD56+ NK (R2 = 0.56) and CD8+T (R2 = 0.48), where 23% of sites were accurately estimated. Hierarchical clustering and pathways enrichment analysis confirmed the biological relevance of the panels. Our results suggest that linear regression for cell-subtype specific methylation estimation is accurate only for some cell-subtypes at a small fraction of cell-associated sites but may be applicable to EWASs of disease traits with a blood-based pathology. Although sample size was a limitation in this study, we suggest that alternative statistical methods will provide the greatest performance improvements.

Publication types

  • Research Support, Non-U.S. Gov't

MeSH terms

  • Blood Cells
  • Cell Lineage / genetics*
  • CpG Islands / genetics
  • DNA / blood*
  • DNA / genetics
  • DNA Methylation / genetics*
  • Epigenesis, Genetic
  • Epigenomics*
  • Female
  • Humans
  • Linear Models
  • Male
  • Signal Transduction / genetics
  • Treatment Outcome

Substances

  • DNA

Grants and funding

RL is partially supported by Multiple Sclerosis Research Australia funding for Bioinformatics. This research is supported by an Australian Government Research Training Program (RTP) Scholarship, and the Australian Research Council Australian Centre for Mathematical and Statistical Frontiers (ACEMS).