Question: which file used for GO
gravatar for fatimarasool135
5 months ago by
fatimarasool13540 wrote:

I have created the following output files by using the stringtie tool 1) Output GTF file containing the assembled transcripts by stringtie

2) A merged GTF file from all sample GTF

I am not sure which of the above file is used further for Gene Ontology by using tool DAVID or agriGO. The output of my file is as:-

Gene ID     Gene Name   Reference   Strand  Start   End     Coverage    FPKM        TPM
  NM_000451   SHOX        chrX        +       624344  646823  0.000000    0.000000    0.000000
  NM_006883   SHOX        chrX        +       624344  659411  0.000000    0.000000    0.000000

Am I supposed to pick the gene Id and expression value to find gene ontology results. Or there is something else I need to do.

rna-seq R • 201 views
ADD COMMENTlink written 5 months ago by fatimarasool13540

You need to perform differential transcript or differential gene expression to then potentially see if there are over-represented gene ontology terms or whatever your goal is. There is somewhat of a guide here if you want to use the developers' of StringTie tool called Ballgown.

ADD REPLYlink modified 5 months ago • written 5 months ago by jean.elbers1.3k

Hi Jean elbers I did differential gene expression by hisat ,string-tie and ballgown. next, i have to do GO ,Singular Enrichment Analysis and KEGG. by using the tool agriGo.I do not know which gene list I have to provide as input to this agriGo tool.

ADD REPLYlink written 5 months ago by fatimarasool13540

list of genes is as

enter code here

feature id fc pval qval gene MSTRG.28632 0.341220148 1.95E-05 0.176761218 gene MSTRG.3615 5.155727198 2.21E-05 0.176761218 gene MSTRG.7507 0.251906716 2.22E-05 0.176761218

I have to use the database AgriGO....According to this tool, gene ids should be in format GenBank ID: AAP50233.1 DDBJ ID: BAB11514.1 EMBL ID: CAA18188.1 UniProt ID: Q9LYA9 RefSeq Peptide ID: NP_564434 PDB ID: 2zsi_B Wheat Affymetrix Genome Array: Ta.10183.2.S1_at but as you have seen mine data set gene list contain different gene ids ,how can I chang it

ADD REPLYlink modified 5 months ago • written 5 months ago by fatimarasool13540
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: 1175 users visited in the last hour