This is a test version of Biostars. For the public version, visit https://www.biostars.org.
get gene ids running sleuth kallisto for differential gene expression

Hi,

I'm trying to run sleuth to look for differential gene expression between two conditions, but I find the manual a bit hard to interpret.

After running kallisto I have extensions in the transcript IDs (like "ENST00000632859.1") and I assume that's why I cannot get the output data mapped to gene IDs while running biomaRt as described in the "getting started" webpage, anyone knows how to deal with this? Can I just remove the extensions in the .tsv file or do I need to change something in the bootstrap file as well?

kallisto sleuth rna-seq
• 4,707 views • link
updated by jocelyn.petitto • jocelyn.petitto • written by frida.danielsson • frida.danielsson

Did you use their transcriptome index? Because that one gives me nicely formatted name, causing no problem to map them to the gene names.

• 1 views • link
written by WouterDeCoster WouterDeCoster

I used a fasta file for the human genome looking like this:

>ENST00000448914.1 havana_ig_gene:known chromosome:GRCh38:14:22449113:22449125:1 gene:ENSG00000228985.1 gene_biotype:TR_D_gene transcript_biotype:TR_D_gene
ACTGGGGGATACG...

Maybe I should try rebuild my index and rerun kallisto...

• 1 views • link
updated by Ram Ram • written by frida.danielsson • frida.danielsson

1 answer

Removing the .1 shouldn't be a problem. Whether it's there or not is due to whether your files came from Ensembl or Gencode (or maybe UCSC) originally.

• 1 views • link
written by Devon Ryan Devon Ryan

That means you can remove it if it's causing problems.

• 1 views • link
written by Devon Ryan Devon Ryan

Log in to answer this question.

More posts like this