search

linnarsson-lab / loompy

Python implementation of the Loom file format - http://loompy.org
BSD 2-Clause "Simplified" License
139 stars 37 forks source link

Unable to retrieve spliced/unspliced RNA from single nuclei isolation #189

Closed hyjforesight closed 14 hours ago

hyjforesight commented 1 day ago

Hello loompy, Thank you for making such an amazing packge!

We’re unable to retrieve spliced/unspliced RNA from single nuclei isolation procedure of scRNA-seq Flex kit. Please find the attached figure. Comparison

ACI sample is used as control. We do single cell isolation of tumor, live cell sorting, then use 3’ chemistry v3.1 kit to do the scRNA-seq. We use CellRanger v8.0.1 to convert Fastq to matrices. Then, use velocyto.py v0.17 to retrieve spliced/unspliced RNA (generate loom file) from CellRanger outputs. We can see 7000 n_genes_by_counts for matrices and 5000 n_genes_by_counts for loom file, which is OK for our RNA velocity analysis. These results validate that these procedures are good to retrieve spliced/unspliced RNA from single cell isolation.

1503095E, 2113752H, 2300648F, 2121949E samples are FFPE samples. We do hydration, dissociation, single nuclei isolation, add probes embedding for multiplexing, then use Flex kit to do the snRNA-seq. We use CellRanger v8.0.1 to convert Fastq to matrices. Then, use velocyto.py v0.17 to retrieve spliced/unspliced RNA (generate loom file) from CellRanger per_sample_outs. We can see n_genes_by_counts of loom file is extreme lower than matrices, which cannot be used for our RNA velocity analysis.

For this situation (unable to retrieve spliced/unspliced RNA for our FFPE samples), could you please help with our questions:

  1. Is velocyto.py not compatible with the single nuclei isolation? Both Fastq files from v3.1 kit and Flex kit are compatible for Cell Ranger, so I don't think it is the kit problem. The splicing of mRNA happens in nucleus, so single nuclei isolation should retrieve good spliced/unspliced RNA, isn’t it?
  2. Did your amazing loompypackage met similar issue or can it retrieve the right spliced/unspliced RNAs?

Thank you! Best,

hyjforesight commented 14 hours ago

It's the problem of Flex kit. See this: https://kb.10xgenomics.com/hc/en-us/articles/25938615598477-Is-trajectory-analysis-compatible-with-the-Fixed-RNA-Profiling-Flex-assay