This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Extract exon Intron reads from RNA seq data

Hello everyone,

I am working with a plant fungi of which I have a DNA sequence. I am trying to make primers for some differentially expressed genes for which I need to know exon-exon junctions. Is there any way I can know the exon intron information of genes I am interested in so that I can easily design the primers.

Many thanks in advance,

Ambika

rna-seq
• 2,734 views • link
written by AP • AP

Have you checked to see if there is a gene annotation file available? You should be able to figure the junctions out from that.

• 0 views • link
written by GenoMax GenoMax

Genomax, I have a gtf file from augustus. But I am not sure how can I infer information out of that. Can you please tell me if thats the file I need and how can I interpret that file

• 0 views • link
written by AP • AP

Does the GTF file have entries for exons (and perhaps introns)?

• 0 views • link
written by GenoMax GenoMax

It has CDS instead of exons

• 0 views • link
written by AP • AP 
tig00007990     AUGUSTUS        CDS     869652  869790  0.90    +       0       transcript_id "g17308.t1"; gene_id "g17308";
tig00007990     AUGUSTUS        CDS     869843  870931  0.99    +       2       transcript_id "g17308.t1"; gene_id "g17308";
tig00007990     AUGUSTUS        CDS     870982  871333  1.00    +       2       transcript_id "g17308.t1"; gene_id "g17308";
tig00007990     AUGUSTUS        CDS     871385  871442  0.99    +       1       transcript_id "g17308.t1"; gene_id "g17308";

The format looks like this

• 0 views • link
written by AP • AP

Those 4 are exons of that one gene transcript (g17308.t1) so presumably there are more than one transcript for each gene (g17308)?

• 0 views • link
written by GenoMax GenoMax

So these numbers represent the starting and ending position of sequence which are exons. I tried to search for this gene in the file and I just have one transcript (g17308.t1) . So can I say this gene have 4 exons?

• 0 views • link
written by AP • AP

If that is all blocks then yes. Check what the number following the strand signifies in specifications. I don't recollect off the top of my head.

• 0 views • link
written by GenoMax GenoMax

0 answers

No answers yet.

Log in to answer this question.

More posts like this