I am a beginner in R I did Differential expressed gene analysis and I have two questions:

1- depending on the graphs is the differentially expressed genes list reliable or not For further analysis?

I am disappointed because I didn't get the segregation of the samples in the heatmap and PCA1 and PCA2 are low

I noticed that the sample S1 varies from other samples resulting in low PCA1 but also different from the 10 controls. is there any way to handle it?!

2- when I am exporting the list of DEGs I get (some genes appears many times) For example, if the gene id is df3t_00100 I got records as following:

df3t_00100

df3t_00100.1

df3t_00100.2

df3t_00100.3

what are those and how can I handle them?

Thank you in advance!

Hi Mohammad,

Regarding the PCA bi-plot, I see no major issue, assuming that you have generated this PCA bi-plot in an unbiased ('unsupervised') way using all genes. Can you share the code that you used? Your 2 groups (Control + Sample) are almost exclusively segregated along PC1. The sample at the bottom-right is behaving differently, but it is still not grouping with Control.

Then, in your second figure generated with pheatmap(), it seems that —yes— your groups are segregated perfectly via hierarchical clustering, and the heatmap colour shade also indicates this.

Regarding the gene naming issue, which species is this? Can you confirm how the read count quantification was performed and with which reference GTF? Generally, to help, please explain your broader analysis pipeline so that we can begin to try to solve this.

Kevin