This is a test version of Biostars. For the public version, visit https://www.biostars.org.
raw matrices or filtered matrices

Hi everyone,

I am working on a small pilot study using single-nucleus RNA-seq in an Activity-Based Anorexia (ABA) mouse model. We have a limited number of samples (5 hypothalamic punches pooled into an ABA sample and 5 pooled into a control sample).

We initially started from the raw Cell Ranger output and applied cell calling, QC filtering, and doublet removal. After these steps, the number of nuclei decreased substantially (approximately 22,000 nuclei per sample down to about 8,000).

I have a few questions:

  1. In a pilot study with limited biological replication, how would you evaluate whether the QC pipeline is removing mostly low-quality nuclei versus potentially removing biologically relevant signal?

  2. How much cell loss after QC and doublet removal would you consider typical or acceptable for snRNA-seq data from hypothalamic tissue?

  3. Would it be reasonable to perform a sensitivity analysis by running the same downstream workflow on both:

    • the raw matrices
    • the filtered matrices,

    and then compare cell type composition, clustering, and major biological findings? or would you pick one over the other

  4. If the major conclusions differ substantially between the two approaches, how would you determine which preprocessing strategy is more appropriate?

Any thoughts or recommendations would be greatly appreciated.

cellranger snrnaseq

The raw matrix is mostly empty droplets and noise. I would not even consider using it. Results will be different for sure, but not out of genuine biology, but because you are analyzing noise. Stick with the filtered one. Probably not what you want to hear, but if you do not have enough cells or depth, then including noise won't change that.

The raw matrix contains all barcodes, including cells/nuclei and empty droplets. In most cases, I would recommend using the filtered matrix, unless you suspect that biologically relevant cells were removed during filtering.

The number of recovered cells/nuclei depends on factors such as cell loading, tissue type, sample quality, and library preparation, ... Likewise, cell loss after QC varies widely between datasets, so it is difficult to give a general estimate.

Hope this helps. If you have any other questions, I'd be happy to help. My GitHub & email: https://github.com/hongthai-lenguyen

0 answers

No answers yet.

Log in to answer this question.