This is a test version of Biostars. For the public version, visit https://www.biostars.org.
How can I get Transcript ID from the gene ID?

Dear all

I have a list of gene IDs in a tabular format. How I can extract the transcript IDs for the list of genes IDs mentioned above, from a gff file?

Thank you so much.

Carlos

rna-seq
• 5,842 views • link
updated by Jeffin Rockey Jeffin Rockey • written by Carlos Caicedo Carlos Caicedo

Depends on what genome this is but you could try BioMart tool from Ensembl.

• 0 views • link
written by GenoMax GenoMax

I have a data from a bacterium specie, so I think BioMart does not function in this case.

• 0 views • link
written by Carlos Caicedo Carlos Caicedo

If this is a bacterium then you should have a single transcript from each gene since there is no alternate splicing, right?

• 0 views • link
written by GenoMax GenoMax

Of course, you are right. I going to try to do a better explanation of my question.

A gff file is something like this:

chromosome  ena gene    661 1041    .   -   .   ID=gene:SCLAV_0001;biotype=protein_coding;description=Hypothetical protein;gene_id=SCLAV_0001;logic_name=ena;version=1
chromosome  ena transcript  661 1041    .   -   .   ID=transcript:EFG05077;Parent=gene:SCLAV_0001;biotype=protein_coding;transcript_id=EFG05077;version=1
chromosome  ena exon    661 1041    .   -   .   Parent=transcript:EFG05077;Name=EFG05077-1;constitutive=1;ensembl_end_phase=0;ensembl_phase=0;exon_id=EFG05077-1;rank=1;version=1
chromosome  ena CDS 661 1041    .   -   0   ID=CDS:EFG05077;Parent=transcript:EFG05077;protein_id=EFG05077

I have a list with ID:gene 

SCLAV_0001
SCLAV_0002

And I need to get for each gene in the list the transcript ID

For instance: 

SCLAV_0001   EFG05077
SCLAV_0002  EFG0XXY

And so on.

• 0 views • link
updated by GenoMax GenoMax • written by Carlos Caicedo Carlos Caicedo

1 answer

Hope the below one liner helps or at least, indicates the way to go ..

awk<yourGeneModel.gff3 -F'\t' '$3=="transcript" {print$9}'  | sed -e 's|ID=transcript:\([^;]*\)\(.*\)Parent=gene:\([^;]*\)\(.*\)|\2\t\1|g'

If it is one gene one transcript for the genemodel, this should do. Else one more script to combine multiple transcripts per gene would be required.

• 0 views • link
written by Jeffin Rockey Jeffin Rockey

Log in to answer this question.

More posts like this