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

Merging mixcr scRNA-seq output with Seurat data #1020

Open
meghatron21 opened this issue Feb 28, 2023 · 1 comment
Open

Merging mixcr scRNA-seq output with Seurat data #1020

meghatron21 opened this issue Feb 28, 2023 · 1 comment

Comments

@meghatron21
Copy link

Hi. I have 10X genomic 5' scRNA-seq fastq files. I ran mixcr like so:
~/mixcr analyze 10x-5gex-full-length --species hsa $sample_R1_001.fastq.gz $sample_R2_001.fastq.gz $output/$sample

I exported clones like so:
~/Programs/mixcr exportClones --split-by-tags Cell -tag cell -cellGroup -uniqueTagCount Molecule $sample.contigs.clns $sample.scRNA.clones.tsv

I see that everything ran without errors, but I'm uncertain on how to map this back to my Seurat dataframe. I tried going through the docs but I'm having trouble. Is there a good resource for this? I essentially want to map each of my cells to likely clonotype.

Thanks,
Meghana

@mehdiborji
Copy link

@meghanasp21
For now, you have to go through the clones.tsv file, group the file by barcodes, then within each barcode write some clean up function where a dominant chain is picked as the most likely present chain in the cell. In that sense this is more complicated than cellranger's output, but assigning confident chains to cells is not super easy even with cellranger's newest version where you can input both RNA and VDJ data.

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