English

Multivariate regression with missing response data for modelling regional DNA methylation QTLs

Methodology 2025-07-09 v1

Abstract

Identifying genetic regulators of DNA methylation (mQTLs) with multivariate models enhances statistical power, but is challenged by missing data from bisulfite sequencing. Standard imputation-based methods can introduce bias, limiting reliable inference. We propose \texttt{missoNet}, a novel convex estimation framework that jointly estimates regression coefficients and the precision matrix from data with missing responses. By using unbiased surrogate estimators, our three-stage procedure avoids imputation while simultaneously performing variable selection and learning the conditional dependence structure among responses. We establish theoretical error bounds, and our simulations demonstrate that \texttt{missoNet} consistently outperforms existing methods in both prediction and sparsity recovery. In a real-world mQTL analysis of the CARTaGENE cohort, \texttt{missoNet} achieved superior predictive accuracy and false-discovery control on a held-out validation set, identifying known and credible novel genetic associations. The method offers a robust, efficient, and theoretically grounded tool for genomic analyses, and is available as an R package.

Keywords

Cite

@article{arxiv.2507.05990,
  title  = {Multivariate regression with missing response data for modelling regional DNA methylation QTLs},
  author = {Shomoita Alam and Yixiao Zeng and Sasha Bernatsky and Marie Hudson and Inés Colmegna and David A. Stephens and Celia M. T. Greenwood and Archer Y. Yang},
  journal= {arXiv preprint arXiv:2507.05990},
  year   = {2025}
}
R2 v1 2026-07-01T03:51:25.523Z