Question about preprocessing ATAC-Seq

Question

Question about preprocessing ATAC-Seq

Opened this issue 9 months ago · 4 comments

Hello,

I'm interested in using scglue to integrate my scRNAseq with scATACseq data, comprising four paired samples. These datasets are not multiome; they're only paired samples.

After processing with Cellranger, I used episcanpy's epi.ct.bld_mtx_fly to construct a count matrix from the TSV/TBI files. Subsequently, I filtered the data using epi.pp.filter_cells and epi.pp.filter_features.

My current concern pertains to the input matrix for the scglue.data.lsi function.

Could you please advise on whether I should use the Raw Matrix, Binarized Matrix, or normalizedMatrix for scglue.data.lsi?

Thank you. Your assistance is appreciated.

Answer 1 · 2024-04-03T11:10:08.000Z

Hi @luglilab. Thanks for your interest in GLUE!

I would recommend using the raw count matrix, as was also recently confirmed in this paper.

Answer 2 · 2024-04-03T15:18:55.000Z

Thanks a lot!

Answer 3 · 2024-04-22T21:40:35.000Z

Hello @Jeff1995 ,

could you suggest a good preprocessing pipeline for 10X data compatible with your tool? I'm trying to figure it out if is better snapatac, episcanpy or signac but i'd like to have your comment.

Thank you.

Answer 4 · 2024-04-24T08:16:04.000Z

Hi @luglilab, I haven't really fiddled with different preprocessing pipelines. As long as it outputs count matrices, it would work well with scglue. scglue has builtin lsi implementation and does not need further processing than the count matrices.