Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

MultiVelo for multipo samples #17

Open
Lei-Guo opened this issue Jan 13, 2023 · 1 comment
Open

MultiVelo for multipo samples #17

Lei-Guo opened this issue Jan 13, 2023 · 1 comment

Comments

@Lei-Guo
Copy link

Lei-Guo commented Jan 13, 2023

Hi,

Thanks for the awesome package. I've got some questions.

  1. How do I use MultiVelo on multiple samples? The demo seems to be for 1 sample only.

  2. The peaks in my scATACseq dataset were called using MACS2. How do I deal with the "feature_linkage.bedpe" files? Should I make my own linkage file using "LinkPeaks"?

@danielee0707
Copy link
Collaborator

Hi, thank you for using the package.
To use multiple samples, one needs to first remove batch effects if there are any. The procedure can be a little complicated, and I wrote the steps we took in discussion #10.
I think LinkPeaks should be able to output an object that describes peak-gene correlations. I have never used this function before, but you may need to manually change the format of the output to mimic the structure of a bedpe file from 10X. Either CellRanger ARC V1 or V2 would work. See https://support.10xgenomics.com/single-cell-multiome-atac-gex/software/pipelines/latest/output/analysis#fl_output for column descriptions. Alternatively, you can start a test run by using an empty bedpe file, so only promoters are considered.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants