m 6 Aexpress-enet: Predicting the regulatory expression m 6 A sites by an enet-regularization negative binomial regression model.

As the most abundant mRNA modification, m6 A controls and influences many aspects of mRNA metabolism including the mRNA stability and degradation. However, the role of specific m6 A sites in regulating gene expression still remains unclear. In additional, the multicollinearity problem caused by the correlation of methylation level of multiple m6 A sites in each gene could influence the prediction performance. To address the above challenges, we propose an elastic-net regularized negative binomial regression model (called m6 Aexpress-enet) to predict which m6 A site could potentially regulate its gene expression. Comprehensive evaluations on simulated datasets demonstrate that m6 Aexpress-enet could achieve the top prediction performance. Applying m6 Aexpress-enet on real MeRIP-seq data from human lymphoblastoid cell lines, we have uncovered the complex regulatory pattern of predicted m6 A sites and their unique enrichment pathway of the constructed co-methylation modules. m6 Aexpress-enet proves itself as a powerful tool to enable biologists to discover the mechanism of m6 A regulatory gene expression. Furthermore, the source code and the step-by-step implementation of m6 Aexpress-enet is freely accessed at https://github.com/tengzhangs/m6Aexpress-enet.