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

empty result when using cellsnp-lite #12

Open
onmyojiyys opened this issue Nov 11, 2024 · 3 comments
Open

empty result when using cellsnp-lite #12

onmyojiyys opened this issue Nov 11, 2024 · 3 comments

Comments

@onmyojiyys
Copy link

Thank you very much for your research in the field of CNA analysis. I have encountered some issues while using CalicoST. When I executed cellsnp-lite with three samples (10x Genomics Visium v2), the output files were all empty, and no SNPs were identified. Could you provide a download link for the bam files in the examples data? I need to confirm whether the issue is with my usage of CalicoST or the data itself. Best regards!

image

image

cellsnp-lite -s /disk/pipeline/STtools/CNV/calicost/MU_4174_F1-0829/possorted_genome_bam.bam -b /disk/pipeline/STtools/CNV/calicost/MU_4174_F1-0829/barcodes_change.txt -O /disk/pipeline/STtools/CNV/calicost/MU_4174_F1-0829/genotyping/ -p 20 --minMAF 0 --minCOUNT 1 --UMItag Auto --cellTAG CB --gzip

@Congm12
Copy link
Contributor

Congm12 commented Nov 11, 2024

The command for cellsnp-lite lacks -R argument for specifying the VCF file of the known SNPs for genotyping. The VCF of known human SNPs can be obtained from the https://cellsnp-lite.readthedocs.io/en/latest/main/data.html#data-list-of-common-snps.

In addition, can you check whether your SRT samples are probe-based? As probe-based transcriptomics doesn't contain SNP information, we recommend running STARCH (https://github.com/raphael-group/STARCH/tree/develop) or other total copy number inference methods.

@onmyojiyys
Copy link
Author

Thank you. After adding candidate SNPs with the -R parameter, the results I obtained are still empty. My data is Visium FFPE, captured with probes.
After downloading the data from the 10x database and completing the analysis with spaceranger count, the calicoST pipeline is running normally.

@onmyojiyys
Copy link
Author

The data used were human lung cancer samples, and the data summary report is as follows.
image

When running the estimate_tumor_proportion.py script, some errors were encountered and a response is hoped for.

Issue with the MIN_TOTAL parameter in the identify_normal_spots function
The configuration_purity configuration file is as follows:
image

When using the default parameters, the following error occurred. Is it because the MIN_TOTAL parameter is set too high, resulting in this_bafs being empty? If so, does this mean the data is not suitable for analysis with the CalicoST method, or can the MIN_TOTAL threshold be appropriately lowered?
image

The default MIN_TOTAL is set to 10. When setting MIN_TOTAL to 2 and rerunning, the same error still occurred.
image

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