Guidelines for cell-type heterogeneity quantification based on a comparative analysis of reference-free DNA methylation deconvolution software

BMC Bioinformatics

Aleksandar Milosavljevic¹  Eugene Lurie¹  Michael Scherer²  Pavlo Lutsik³  Eugene Andres Houseman⁴  Arthur Waguet⁵  Raphael Bacher⁵  Michael G. B. Blum⁵  Daniel Jost⁵  Clémentine Decamps⁵  Magali Richard⁵  Florian Privé⁵ 

[1] Bioinformatics Research Laboratory, Molecular and Human Genetics Department, Baylor College of Medicine, Houston, TX, USA;Department of Genetics/Epigenetics, Saarland University, 66123, Saarbruecken, Germany;Division of Cancer Epigenomics, German Cancer Research Center (DKFZ), Heidelberg, Germany;Independent Statistical Consultant, La Center, WA, USA;Laboratory TIMC-IMAG, UMR 5525, Univ. Grenoble Alpes, CNRS, F-38700, Grenoble, France;

关键词: Cell heterogeneity;    Deconvolution;    DNA methylation;    Epigenetics;    Matrix factorization;    R package/pipeline;   

DOI  :  10.1186/s12859-019-3307-2

来源: Springer

【 摘 要 】

BackgroundCell-type heterogeneity of tumors is a key factor in tumor progression and response to chemotherapy. Tumor cell-type heterogeneity, defined as the proportion of the various cell-types in a tumor, can be inferred from DNA methylation of surgical specimens. However, confounding factors known to associate with methylation values, such as age and sex, complicate accurate inference of cell-type proportions. While reference-free algorithms have been developed to infer cell-type proportions from DNA methylation, a comparative evaluation of the performance of these methods is still lacking.ResultsHere we use simulations to evaluate several computational pipelines based on the software packages MeDeCom, EDec, and RefFreeEWAS. We identify that accounting for confounders, feature selection, and the choice of the number of estimated cell types are critical steps for inferring cell-type proportions. We find that removal of methylation probes which are correlated with confounder variables reduces the error of inference by 30–35%, and that selection of cell-type informative probes has similar effect. We show that Cattell’s rule based on the scree plot is a powerful tool to determine the number of cell-types. Once the pre-processing steps are achieved, the three deconvolution methods provide comparable results. We observe that all the algorithms’ performance improves when inter-sample variation of cell-type proportions is large or when the number of available samples is large. We find that under specific circumstances the methods are sensitive to the initialization method, suggesting that averaging different solutions or optimizing initialization is an avenue for future research.ConclusionBased on the lessons learned, to facilitate pipeline validation and catalyze further pipeline improvement by the community, we develop a benchmark pipeline for inference of cell-type proportions and implement it in the R package medepir.

【 授权许可】

CC BY