This is a test version of Biostars. For the public version, visit https://www.biostars.org.
How to remove a gene before making a cellranger reference?

Hello,

I was asked to make a new reference where I remove one of the genes from a reference and add some isoforms. I have found several instances of this gene in the gtf file and I will remove them using regex. Do I need to remove anything from the fasta file too?

Thank you

fasta gtf cellranger scrnaseq
• 1,755 views • link
updated by ATpoint ATpoint • written by bioinfo • bioinfo

1 answer

No, the fasta is just the genome so you don't have to do anything.

• 0 views • link
written by ATpoint ATpoint

Hello,

Thank you very much for replying. For the fasta file then it is ok to just add the isoform sequences?

Thank you

• 0 views • link
written by bioinfo • bioinfo

Why would you add isoform sequences to a reference genome file? A reference genome contains the sequence for all possible isoforms in all cells.

• 0 views • link
written by Ram Ram

I think you need some background on the basics. A genome is what is encoded as DNA on the chromosomes as plain sequence so any isoform is already there. It is the GTF that guides the aligner and the quantification at which position we hsve features such as intron gaps, exons etc. I am pretty sure 10x has a guide towards creation of custom CellRanger references.

Edit: https://support.10xgenomics.com/single-cell-gene-expression/software/pipelines/latest/using/tutorial_mr

• 0 views • link
written by ATpoint ATpoint

Log in to answer this question.

More posts like this