Question: Extrapolate transcripts that dictate difference between samples in PCA
1
gravatar for lucatucciarone
5 months ago by
lucatucciarone50 wrote:

I'm interested in extrapolating what transcripts dictates differences between sample in a PCA from RNA-seq data. This is the PCA of RNA-seq on different samples (A,B,C,D,E,F,G):

Samples are in duplicate and are circled for better understanding. Bold samples are controls, samples connected with arrows are treated samples.

What I want to understand is what are the genes that "moves" A to B (Red Arrow) or A to C (Blue arrow). Given the partially similar trajectory of A to C and D to E the second question would be to check whether some genes in A to C are also responsible of the movement of D to E. This same comparison can be made from comparing A to B and F to G.

Do you happen to know a tool that can allow me to better address this questions?

rna-seq R • 234 views
ADD COMMENTlink modified 5 months ago • written 5 months ago by lucatucciarone50

I have thought of a strategy to partially address my question:

My observation starts with A to C and D to E going into the same direction.
Now, one way to answer that question would be to: Extrapolate from both PCs (1 and 2) genes that are most variable (and thus dictates the PC1 and 2) Check the expression of those gene in: A Vs C D Vs E and select for genes that has a similar expression pattern.

Do you think that this could be a valid strategy? What tools can you suggest me?

ADD REPLYlink written 5 months ago by lucatucciarone50

Someone please!!?!?!?!

ADD REPLYlink written 5 months ago by lucatucciarone50
2
gravatar for geek_y
5 months ago by
geek_y9.8k
Barcelona
geek_y9.8k wrote:

First of all, there are no trajectories in PCA. It's not a graph/network.

You can get the 'n' list of genes that are most variable between any two groups by rowVars function that is part of several R packages. These highly variable genes are the drivers of PCA. Usually top 500 genes are considered for PCA.

ADD COMMENTlink written 5 months ago by geek_y9.8k

This is very interesting reply, thank you. I'd have one more question on the subject:

Given that there are no trajectory, do you think is legitimate to think that " A to C and D to E" may have some genes that go in the same direction and for that reason their "path/differences" may be similar? In other words, can I at least say that the effect that C has on A may be have some genes overlapping what E does on D? If not, what can I say about the similarity of A to C and D to E? If it is possible to see any

ADD REPLYlink written 5 months ago by lucatucciarone50

You may need to do differential analysis on all pairs and try to build a network from differential genes.

ADD REPLYlink written 5 months ago by geek_y9.8k

Would you kindly suggest me some papers/tools that explains this pipeline that you are suggesting?

ADD REPLYlink written 5 months ago by lucatucciarone50

First you need to accept the existing answer if it clarifies your question.

ADD REPLYlink written 5 months ago by geek_y9.8k

Oki! Sorry about that, I am new to the website : )

ADD REPLYlink written 5 months ago by lucatucciarone50

any news on this subject?

ADD REPLYlink written 5 months ago by lucatucciarone50
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1316 users visited in the last hour