Abstract:Background: RNA methylome has been discovered as an important layer of gene regulation and can be profiled directly with count-based measurements from high-throughput sequencing data. Although the detailed regulatory circuit of the epitranscriptome remains uncharted, clustering effect in methylation status among different RNA methylation sites can be identified from transcriptome-wide RNA methylation profiles and may reflect the epitranscriptomic regulation. Count-based RNA methylation sequencing data has unique features, such as low reads coverage, which calls for novel clustering approaches. Objective: Besides the low reads coverage, it is also necessary to keep the integer property to approach clustering analysis of count-based RNA methylation sequencing data. Method: We proposed a nonparametric generative model together with its Gibbs sampling solution for clustering analysis. The proposed approach implements a beta-binomial mixture model to capture the clustering effect in methylation level with the original count-based measurements rather than an estimated continuous methylation level. Besides, it adopts a nonparametric Dirichlet process to automatically determine an optimal number of clusters so as to avoid the common model selection problem in clustering analysis. Results: When tested on the simulated system, the method demonstrated improved clustering performance over hierarchical clustering, K-means, MClust, NMF and EMclust. It also revealed on real dataset two novel RNA N6-methyladenosine (m6A) co-methylation patterns that may be induced directly by METTL14 and WTAP, which are two known regulatory components of the RNA m6A methyltransferase complex. Conclusion: Our proposed DPBBM method not only properly handles the count-based measurements of RNA methylation data from sites of very low reads coverage, but also learns an optimal number of clusters adaptively from the data analyzed. Availability: The source code and documents of DPBBM R package are freely available through the Comprehensive R Archive Network (CRAN): https://cran.r-project.org/web/packages/DPBBM/.

A Multi-Poisson Dynamic Mixture Model to Cluster Developmental Patterns of Gene Expression by RNA-seq

A Bi-Poisson Model for Clustering Gene Expression Profiles by RNA-seq.

How to Cluster Gene Expression Dynamics in Response to Environmental Signals

A sparse negative binomial mixture model for clustering RNA-seq count data

Bayesian clustering of replicated time-course gene expression data with weak signals

A Computational Algorithm For Functional Clustering Of Proteome Dynamics During Development

A Multivariate Poisson-Log Normal Mixture Model for Clustering Transcriptome Sequencing Data

Clustering of time-course gene expression profiles using normal mixture models with AR(1) random effects

Differential Expression Analysis for RNAseq Using Poisson Mixed Models.

A bayesian multivariate mixture model for high throughput spatial transcriptomics

A Skellam Model to Identify Differential Patterns of Gene Expression Induced by Environmental Signals

Time-variant Clustering Model for Understanding Cell Fate Decisions.

Clustering gene expression time series data using an infinite Gaussian process mixture model

Clustering gene expression data based on predicted differential effects of GV interaction.

Bi-clustering interpretation and prediction of correlation between gene expression and protein abundance

Visualization, benchmarking and characterization of nested single-cell heterogeneity as dynamic forest mixtures

Single-Cell Transcriptome Data Clustering via Multinomial Modeling and Adaptive Fuzzy K-Means Algorithm

Integrative Model-based clustering of microarray methylation and expression data

ZINBMM: a general mixture model for simultaneous clustering and gene selection using single-cell transcriptomic data

Clustering Count-based RNA Methylation Data Using a Nonparametric Generative Model

A probabilistic model-based bi-clustering method for single-cell transcriptomic data analysis