RNAseq analysis for aspergillus niger ?
1
0
Entering edit mode
5 weeks ago
najibveto ▴ 50

hello, I am working on analysis of RNAseq from three strains of aspergillus niger with three replicates for each condition (treatment and control). I have the count file for each strain as table limited file and the sequences for each transcripted are stored as text file with their transcript id. I want to use R for analysis the differential expressed gene in each strain, also combine the same transcripts from each strain in file to compare their expression through use of heatmap. further, I would like to do gene ontology analysis for the different genes for each strain as the organism is non-model organism. could someone help me in this task? thank you for your help.

aspergillus niger RNAseq • 474 views
ADD COMMENT
1
Entering edit mode
5 weeks ago
Mensur Dlakic ★ 12k

There is a series of Perl script that come bundled with Trinity installation that should be suitable for your task. For an explanation on how to prepare the files and what each script can do:

https://github.com/trinityrnaseq/trinityrnaseq/wiki/Trinity-Differential-Expression

ADD COMMENT
0
Entering edit mode

thank you a lot for your help, Just I don't have raw data and the gene ontology for aspergillus niger is lacking. so is it possible to use trinity for that purpose?

ADD REPLY
1
Entering edit mode

Hi

Can you show us a bunch of transcript IDs from each strain?

edit: Is important to know if for transcript quantification a reference strain of aspergillus niger was used.

ADD REPLY
0
Entering edit mode

thank you for your answer, the reference strain was N402 and its transcript and count as fellow enter image description here

for the second strain KCN5:

enter image description here

for the third strain KJC3:

enter image description here

and for each strain, there is separated file with gene id and their sequences.

ADD REPLY
1
Entering edit mode

If the reference strain was N402, then KCN5 and KJC3 transcript IDs should look like those in Table 1 (N402).

also combine the same transcripts from each strain in file to compare their expression through use of heatmap.

This is only possible if the N402, KCN5 and KJC3 count tables share the same transcript ID.

for each strain, there is separated file with gene id and their sequences.

This is not very helpful since you have transcript level counts

ADD REPLY
0
Entering edit mode

thank you for your reply, how is it possible to combine the same transcripts based on sequences? there is a tool or script for it?

ADD REPLY
0
Entering edit mode

any homology-based approach would do the job but it is a subpar solution in my opinion. I would contact the guy who did the mapping for transcript quantification and ask him to keep the transcript IDs of the reference genome.

ADD REPLY
0
Entering edit mode

thank you, just the person who generated this data isn’t around anymore. is there a tutorial for similar analysis?

ADD REPLY
0
Entering edit mode

thank you, just the person who generated this data isn’t around anymore.

I don't get this. It can't possibly be easier to ask on the internet than to get in touch with a person who worked in the same lab. Did they not leave any notes or contact address?

ADD REPLY
0
Entering edit mode

Yes, that is the problem, there is no notes plus the person who generate the data didn't left any contact and didn't left any note on how he generate the data counts. can you please let me know which format of table can be made to be used for R? like should i make table with the three strains and their counts number? for gene ontology, should have same processes? thank you and sorry for the trouble.

ADD REPLY
1
Entering edit mode

You said that you have counts - that should work with these scripts. The only other thing you need is a feature file (GFF or GTF) and you are good to go for DGE analysis. You will need something for gene ontology - Trinity is a de novo RNAseq assembler.

ADD REPLY
0
Entering edit mode

thank you for the reply, except count and genes sequences, I don't have the feature file. for gene ontology, is it possible to use the reference strain N402?

ADD REPLY

Login before adding your answer.

Traffic: 2174 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6