This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Capped expression values in single cell RNAseq (Tabula Sapiens/scanpy/CellXGene)

Hi all,

If anyone familiar with using Tabula Sapiens or scanpy could address a question I'd be grateful. I'm looking into using Tabula sapiens (10x genomics data across several organs) to check out some cell-specific markers but finding that some of them appear to have their expression values capped. Looks like a small proportion of cells overall have expression values of 10.00 which leads to some funny looking distributions e.g. the below violin plot of normalised expression values for ANKRD1 - a gene which has these values capped at 10.00 in Tabula sapiens:

image: violin plot

Hard to tell from the manuscript but seems likely the data was normalised with scanpy, which I haven't used before myself. Would this be the source of the capped data? Not seeing any info anywhere on why this data looks like this. It seems likely to me to lead to some quite skewed diff. expression results and logFC values...

scanpy cellxgene scrna-seq 10x normalisation
• 883 views • link
written by matt.a.bennett25890 • matt.a.bennett25890

1 answer

For anyone falling down a similar rabbit hole in the future: my working hypothesis is that this was done to facilitate easier visualisation on the UMAP. To prevent extreme expression outliers from obscuring the color scale for the vast majority of remaining values.

But probably means the values should be re-computed if doing any type of differential expression analysis on the normalised counts for Tabula Sapiens data

• 0 views • link
written by matt.a.bennett25890 • matt.a.bennett25890

Log in to answer this question.

More posts like this