No base count matrix(.tsv) files in Step 2 output but bam files of those cell types are there in Step1 Output.
lipikakalson opened this issue · 16 comments
I am trying to detect somatic mutations from my single-cell RNA-seq data for lung cancer samples.
Step 1 command used:
Output files:
Step 2:
I am not sure why there is no output files for all other cell types. Is there a specific reason that 10 out of 12 files are not there.
But all the files are processed (attaching the ss of subset of output of Step2)
Note: I also tried using the Example data and its the same issue:
Only here, one file is missing.
Hi, I am finding the same problem here, The step 1 runs with no problems and I have the bams files, however the step 2 code is generating the temporary folders (I have removed in the code the rm -rf $temp line to check what was going on on those folders so I am keeping them, but they are empty).
However, it never adds anything inside the temp folders and the output file is not produced. It doesn't throw an error anywhere in the process either.
Yes, I tried it with other sample also, no file is produced now in Step 2 output even though it runs without any error. Also, all my temp files are empty.
Dear users,
Could you check couple of things?
- Could you show the size of the bam files. For instance using
ls -lh *bam
. - Could you paste here the output (*report.txt) of the SplitBamCellTypes.py ?
- Could you run the BaseCellCounter.py with the parameter
--nprocs 1
and print the error? - Finally, could you confirm that the reference genome version used in the SComatic computation is the same as the one used for the alignment (CellRanger)?
Thanks,
Fran
Hi,
It seems that the first step (SplitBamCellTypes.py) filtered a huge number of reads (66M out of 70M) because of CBs not being found in the metadata file. As consequence, step2 files are empty because there are almost no PASS reads to analyse (see column Pass_reads in the report.txt file). The CBs provided in the medatada file do rarely match with the ones observed in the bam file.
Cheers,
Fran
Hi, thank you for the help.
I figure out the error. As you said it was a problem with the reads being filtered... in my case it was that I left the nM filter and my BAM file don't have it. Now, I am running the STEP 2 and looks like is working... I'll update with the results.
Hi, It seems that the first step (SplitBamCellTypes.py) filtered a huge number of reads (66M out of 70M) because of CBs not being found in the metadata file. As consequence, step2 files are empty because there are almost no PASS reads to analyse (see column Pass_reads in the report.txt file). The CBs provided in the medatada file do rarely match with the ones observed in the bam file.
Cheers, Fran
I see the issue now, thank you for your reply and help.
Is there any way I can see which CB are found, which are mismatched or not found?
Best,
Lipika
Hi,
You could grep
the CB columns in the bam file using this command samtools view scRNAseq.bam | grep -o 'CB:Z:[ACTGN]*-1' | sort | uniq > unique.cb.txt
and then compare the output file with your metadata.
You might need to slightly change the previous command depending on how your CBs look like.
Thanks,
Fran
Hi, It seems that the first step (SplitBamCellTypes.py) filtered a huge number of reads (66M out of 70M) because of CBs not being found in the metadata file. As consequence, step2 files are empty because there are almost no PASS reads to analyse (see column Pass_reads in the report.txt file). The CBs provided in the medatada file do rarely match with the ones observed in the bam file.
Cheers, Fran
Hi Fran, I tried with my other sample, in this case there is almost 85% of passed reads, still there is no file in Step 2. I am attaching the ss here. Also I could not see the unique CB of these bam files.
Could you paste the command you used to run for instance the endothelial cells?
Hi Fran,
I used the same command for Endothelial cells.
I managed to get past this error by using meta file from annotating raw feature matrix without doing any QC or any filtering on it instead of doing the annotation later, so that no reads are removed.
I am guessing that the 4 cell types did not have enough passed reads, that's why those 4 cell types are missing in Step 2.
If I do Step 3 now, there is no file in this output even it completes without any error.
Could you paste the command you used to run the step 3? Why does it say that there are 15 tsv files in this folder? There should be only the ones you want to be merged for this sample/individual.
Have you check the hidden files in your folder? (Ctrl + H). It happened to me that I forgot to define the Sample object in the terminal and the script generated a file that started with a period, which in Linux it makes the file hidden
Yes, thanks, it worked now. It was hidden.
Have you check the hidden files in your folder? (Ctrl + H). It happened to me that I forgot to define the Sample object in the terminal and the script generated a file that started with a period, which in Linux it makes the file hidden
Could you paste the command you used to run the step 3? Why does it say that there are 15 tsv files in this folder? There should be only the ones you want to be merged for this sample/individual.
Yes, i was guessing 15 files are the unique cell types (base count matrix) tsv files, that will be merged into Step3.
I successfully ran the SComatic without any errors. Thankyou so much for all the help!
One more question, how to visualise these results graphically? (Sorry, I am very new to scRNA analysis :( )
Thanks in advance!
We used to perform the downstream analysis using R or Python depending on your priorities. However, this is not included in SComatic.
Cheers,
Fran
Thank you so much for your help :)