splatPop: simulating population scale single-cell RNA sequencing data

Abstract

With improving technology and decreasing costs, single-cell RNA sequencing (scRNA-seq) at the population scale has become more viable, opening up the doors to study functional genomics at the single-cell level. This development has lead to a rush to adapt bulk methods and develop new single-cell-specific methods and tools for computational analysis of these studies. Many single-cell methods have been tested, developed, and benchmarked using simulated data. However, current scRNA-seq simulation frameworks do not allow for the simulation of population-scale scRNA-seq data. Here, we present splatPop, a new Splatter model, for flexible, reproducible, and well documented simulation of population-scale scRNA-seq data with known expression quantitative trait loci (eQTL) effects. The splatPop model also allows for the simulation of complex batch effects, cell group effects, and conditional effects between individuals from different cohorts.