RNA seq - TMM values
0
0
Entering edit mode
2.9 years ago

Hello!!

I did an RNA seq and then I did the normalization of TMM. Afterwards I chose:

  1. 1 contig which was over-expressed
  2. my 4 plants populations
  3. the TMM values regarding this specific contig and the 4 populations.

So I thought to do an ANOVA test to check for any st.sign. diff. among the populations for the specific contig. In the end I did a t test with Bonferroni adjustment.

Do you think that this whole thing I did make sense? If not what you could recommend to me?

Btw I tried the heatmap but I think my excel set up is somehow wrong...

Thank you very much in advance!

RNA-Seq tmm • 1.1k views
ADD COMMENT
1
Entering edit mode

Any reason why you did not use an algorithm such as DESeq2 or edgeR for your analysis?

ADD REPLY
0
Entering edit mode

I did. I used edgeR...

ADD REPLY
1
Entering edit mode

Then why the ANOVA test? I'm afraid I don't understand your question.

ADD REPLY
0
Entering edit mode

Sorry. It s my fault. I am beginner in Bioinformatics. I did it with edgeR. And then I did an exact test analysis in Blast2GoPro based on some threshold for logFC, FDR, p value. And I took the contigs with FDR<0.05 as upregulated. Now I want to demostrate them as a figure. So now I can plot the TMM as a boxplot for my different populations and write on the top of my plot the FDR value. But I want also to show statistically significant differences btw the populations. The FDR is only about the contig. I dont have a unique FDR for each population. So my question is: How to show st.sign.differences btw the populations regarding the same contig? Sorry for the confusion....

ADD REPLY
1
Entering edit mode

I think the question was why didn't you go ahead and use edgeR/limma for the statistics part? Most people will use the linear modeling in limma/DESeq2 for statistical testing.

ADD REPLY
0
Entering edit mode

I dont know what are these stuff...sorry...

ADD REPLY
1
Entering edit mode

Hi evlampia123x,

there are a few things that need clarification: What exactly do you mean by contig? Is contig = gene? So far, I understand that you did RNA-seq but which groups did you compare with edgeR? Do you have 4 populations under different conditions or compare 2vs2? What is the experimental setup? Please add some details to help understand your problem.

ADD REPLY
0
Entering edit mode

I did RNA sequence using untreated plants coming from 4 different populations. I had 6 plants from each population. So overall 24 untreated plants. Then I took the results of the sequence and aligned them against a reference transcriptome. Then I did an BWA-MEM mapping. Then I took the raw data results and then I took the TMM (normalization) values (all these took place online in GeneData). When I say contig I mean contig which corresponds to a gene (let s say contig=gene). And I compared everything with edgeR (GeneData). Like this I took the TMM values. Are these infos enough? Sorry for the confusion. I am a beginner....if you cannot help me I will understand it...it s my fault bcs I cannot give you enough details because I dont know many things. I did the analysis with an other person...My question actually is: Is it correct to an ANOVA for TMM values or this is absolutely bulshit? Thank you guys for your help!!

ADD REPLY

Login before adding your answer.

Traffic: 1588 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