How to choose the important clusters for RNAseq data of different groups for downstream analysis/annotation
1
0
Entering edit mode
18 months ago
alwayshope ▴ 40

Dear experts,

Could you shed some light on the choice of the clusters of genes in bulk RNAseq analysis?

The mfuzz package can do cluster analysis for groups, while how to choose which cluster is the most important one for the downstream analysis?

Some clusters have a lot of genes that have constant expression but also include a small number of genes that have increased expression. Is there a need to check all the clusters or only focus on the cluster with constantly increased/decreased genes?

Thank you very much!

enter image description here

RNA-seq cluster • 893 views
ADD COMMENT
1
Entering edit mode
18 months ago
rfran010 ★ 1.3k

I don't really consider myself an expert, but these clusters look strange to me. Did you normalize the data specifically for the clustering Looks like mfuzz may offer this function with

 standardise()

Besides that consideration, you would normally find a cluster that just looksinteresting, and that might include all of them. In a way, as the researcher, you will annotate the clusters. For example, I had 5 clusters and after observing them, I determined there to be 1) early upregulated, 2) late upregulated, 3) transiently upregulted, 4) downregulated, and 5) no change. After that though, what is interesting depends on what your experimental question is.

ADD COMMENT
0
Entering edit mode

Thanks a lot for your guidance!

I add the standardise() to the eset for the control sample group(average the values of the replicates for each time point). The results are as below. As you can see, some increased genes can be shadowed by the constantly expressed genes, not so obvious as in cluster3, obviously shadowed as in cluster1, I saw the process when generating the cluster plot, the genes increased at 24h actually have very low changes at 5 and 10h, then the constantly expressed lines ploted and covered the lines of the few increased genes.

enter image description here

ADD REPLY

Login before adding your answer.

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