Question: How to search for gene isoforms in FASTQ file
gravatar for nkabo
15 days ago by
nkabo50 wrote:

Dear all,

I am a very beginner of RNA-seq analysis. I had completed the RNA-seq process by following the workflow in the link below; I had used Salmon and tximport:

After this part, I had used Ensembl for transcript annotation, but I only get the names for genes (not for transcripts) as shown below:

  1. ENSG00000000003
  2. ENSG00000000005
  3. ENSG00000000419
  4. ENSG00000000457 etc.

However, I was asked to detect whether a specific isoform of a given gene is present or not in the samples. The gene is present in the final result table (as an "average" of the present isoforms I assume). However, the result table provided me only with list of genes and not the transcripts, isoforms are not annotated.

It might be a silly question, but is it possible to obtain FASTA format of this desired specific isoform of the given gene and search it into the FASTQ file of the samples? Is it logical? Or could you suggest me a better way? I have the quantification and abundances, maybe I should use another annotation method to obtain transcript IDs but I am not sure about which part I should change. Thanks in advance

ADD COMMENTlink modified 15 days ago by ATpoint42k • written 15 days ago by nkabo50
gravatar for ATpoint
15 days ago by
ATpoint42k wrote:

tximport aggregates or "summarizes" the isoforms to the gene level, so it is the sum rather than the average. That is the whole point of tximport so take the transcript level counts and make it ready for gene level analysis.

If you want to check the expression level of the isoforms you can simply take the TPMs per transcript produced by salmon in the quant.sf file in every of the quantification directories per sample.

ADD COMMENTlink written 15 days ago by ATpoint42k

Thanks a lot, it worked.

ADD REPLYlink written 14 days ago by nkabo50
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1205 users visited in the last hour