Abstract:Motivation: The negative binomial distribution has been shown to be a good model for counts data from both bulk and single-cell RNA-sequencing (RNA-seq). Gaussian process (GP) regression provides a useful non-parametric approach for modelling temporal or spatial changes in gene expression. However, currently available GP regression methods that implement negative binomial likelihood models do not scale to the increasingly large datasets being produced by single-cell and spatial transcriptomics. Results: The GPcounts package implements GP regression methods for modelling counts data using a negative binomial likelihood function. Computational efficiency is achieved through the use of variational Bayesian inference. The GP function models changes in the mean of the negative binomial likelihood through a logarithmic link function and the dispersion parameter is fitted by maximum likelihood. We validate the method on simulated time course data, showing better performance to identify changes in over-dispersed counts data than methods based on Gaussian or Poisson likelihoods. To demonstrate temporal inference, we apply GPcounts to single-cell RNA-seq datasets after pseudotime and branching inference. To demonstrate spatial inference, we apply GPcounts to data from the mouse olfactory bulb to identify spatially variable genes and compare to two published GP methods. We also provide the option of modelling additional dropout using a zero-inflated negative binomial. Our results show that GPcounts can be used to model temporal and spatial counts data in cases where simpler Gaussian and Poisson likelihoods are unrealistic. Availability and implementation: GPcounts is implemented using the GPflow library in Python and is available at https://github.com/ManchesterBioinference/GPcounts along with the data, code and notebooks required to reproduce the results presented here. The version used for this paper is archived at https://doi.org/10.5281/zenodo.5027066. Supplementary information: Supplementary data are available at Bioinformatics online.

Bayesian modeling of spatial molecular profiling data via Gaussian process

A bayesian multivariate mixture model for high throughput spatial transcriptomics

Bayesian hidden mark interaction model for detecting spatially variable genes in imaging-based spatially resolved transcriptomics data

Bayesian statistics in spatial epidemiology]

A Bayesian modified Ising model for identifying spatially variable genes from spatial transcriptomics data

High-dimensional Bayesian Model for Disease-Specific Gene Detection in Spatial Transcriptomics

Bayesian Flexible Modelling of Spatially Resolved Transcriptomic Data

An interpretable Bayesian clustering approach with feature selection for analyzing spatially resolved transcriptomics data

Non-parametric modelling of temporal and spatial counts data from RNA-seq experiments

Joint identification of spatially variable genes via a network-assisted Bayesian regularization approach

Spanve: an Statistical Method to Detect Clustering-friendly Spatially Variable Genes in Large-scale Spatial Transcriptomics Data

Semi-Supervised Non-Parametric Bayesian Modelling of Spatial Proteomics

Bayesian Nonparametric Clustering with Feature Selection for Spatially Resolved Transcriptomics Data

Generalized Bayesian nonparametric clustering framework for high-dimensional spatial omics data

Latent Space Inference For Spatial Transcriptomics

A Bayesian mark interaction model for analysis of tumor pathology images

Spatial transcriptomics at subspot resolution with BayesSpace

A Bayesian spatial model for imaging genetics

A robust statistical approach for finding informative spatially associated pathways

Bayesian Integrative Region Segmentation in Spatially Resolved Transcriptomic Studies

Joint Bayesian estimation of cell dependence and gene associations in spatially resolved transcriptomic data