Hello, I have 8 samples (4 WT and 4 mutants). I have performed bulk RNA seq on these (Kallisto, and then DESeq2) I am attaching my PCA plot, I was expecting the wild type to cluster together and all the mutants to be together, but this was not the result. Does it indicate that my samples have a lot of variance?
Can anyone suggest to me any other method/plot that can represent my WT vs mutant samples?
I would really appreciate it.
Right now, the WT samples are clustering together, but you have one mutant sample off by itself, which could be an outlier. How did you scale the data? If you used Z-scores, I would suggest trying center log ratios.
Besides PCA, you could also try hierarchical clustering.
Everything that Jeremy said, plus that there may be no data manipulation that will give you the clusters you expected. Sometimes there are outliers in the data, and that won't change no matter how much we torture the data.
Log in to answer this question.
How much variance were explained by PC1 and PC2? It doesn't look like that the largest source of variation came from the genotype difference.
Cross posted here (including the code used to make up the counts):
https://support.bioconductor.org/p/9157911/#9157959