eQTL studies are essential for understanding genomic regulation. The effects of genetic variations on gene regulation are cell-type-specific and cellular-context-related, so studying eQTLs at a single-cell level is crucial. The ideal solution is to use both mutation and expression data from the same cells. However, the current technology of such paired data in single cells is still immature. We present a new method, eQTLsingle, to discover eQTLs only with single-cell RNA-seq (scRNA-seq) data, without genomic data. It detects mutations from scRNA-seq data and models gene expression of different genotypes with the zero-inflated negative binomial (ZINB) model to find associations between genotypes and phenotypes at the single-cell level. On a glioblastoma and gliomasphere scRNA-seq dataset, eQTLsingle discovered hundreds of cell-type-specific tumor-related eQTLs, most of which cannot be found in bulk eQTL studies. Detailed analyses on examples of the discovered eQTLs revealed important underlying regulatory mechanisms. eQTLsingle is a uniquely powerful tool for utilizing the vast scRNA-seq resources for single-cell eQTL studies, and it is available for free academic use at https://github.com/horsedayday/eQTLsingle.
