Hello,
I relatively new to programming with R and I was wondering if someone could help me in providing suggestions on annotating proteins in a heatmap. Here I have tried the following code, but the heatmap generated is not color coded correctly. I am also unsure how to add multiple annotations to a single protein. E.g. if I wanted to annotate Protein 4 with all 3 paths or with 2 specific paths how could I do this? Any help would be greatly appreciated.
# Read the CSV file
cultureconditions <- read.csv("Sample proteomic data3.csv", row.names = 1)
# These are the names of the columns
new_column_names <- c("condition1", "condition2", "condition3", "condition4")
colnames(cultureconditions) <- new_column_names
#Generate annotation for rows
annotation_row = data.frame(
GenePaths = factor(rep(c("Path1", "Path2", "Path3"), c(17, 5, 18)))
)
#Color code annotations
ann_colors = list(GenePaths = c(Path1 = "#7570B3", Path2 = "#E7298A", Path3 = "#66A61E")
)
# Create the heatmap
# scale=row allows for normalization with z scores
#cluster_cols prevents the columns from clustering
heatmap <- pheatmap(cultureconditions, scale = "row", cluster_row = FALSE, cluster_cols = FALSE, annotation_row = annotation_row,
annotation_colors = ann_colors, angle_col = 45, fontsize = 16)
This is the result I get:
1 answer
Your annotation_row dataframe needs to have the rows named. This is how pheatmap matches which rows to annotate. Note your row names in your original dataframe need to match exactly the row names in the annotation_row. So Protein1 and Protein 1 are different.
annotation_row = data.frame(
GeneClass = factor(rep(c("Path1", "Path2", "Path3"), c(10, 4, 6)))
)
rownames(annotation_row) = paste("Protein ", 1:20, sep = "")
I also noticed your annotation_row has the wrong dimensions. It should be 20 rows like your first dataframe.
To improve your code I suggest these edits:
Instead of:
new_column_names <- c("condition1", "condition2", "condition3", "condition4")
colnames(cultureconditions) <- new_column_names
Just do:
colnames(cultureconditions) <- c("condition1", "condition2", "condition3", "condition4")
Thank you very much again! I will give it a try. Also, my apologies but could clarify which pheatmap doc you are referring to? This has been my reference for pheatmaps currently: https://jokergoo.github.io/ComplexHeatmap-reference/book/integrate-with-other-packages.html
The reference manual has a huge tutorial. See page 5 https://cran.r-project.org/web/packages/pheatmap/pheatmap.pdf
Thank you very much! I did look through this doc. I think the trouble that I am having is that I would like to use the same data frame (GenePaths) to annotate multiple proteins. In other words, I would like to annotate select proteins to suggest that they are a part of Path1 and Path2 etc. Would you happen to know if I can do this with the same data frame?
There's two ways to do this. I would recommend the what I suggested above about creating another item in the dataframe. Another way is to add another factor to your items to indicate shared paths. eg: factor(rep(c("Path1", "Path2", "Path3", "Path1.Path3"), c(17, 5, 18, 3))), here "Path1.Path3" indicate a shared path.
Don't forget to create a color for the extra items you create.
Thank you very much! Apologies for the repeated questions. I tried the second option that you suggested but it gave the following error shown below.
Error in .rowNamesDF<-(x, value = value) : invalid 'row.names' length
Here is my code below:
colnames(cultureconditions) <- c("condition1", "condition2", "condition3", "condition4")
#Generate annotation for rows
annotation_row = data.frame(
GenePaths = factor(rep(c("Path1", "Path2", "Path3","Path1.Path3"), c(17, 5, 18,3)))
)
rownames(annotation_row) = paste("Protein ", 1:40, sep = "")
#Color code annotations
ann_colors = list(GenePaths = c(Path1 = "#7570B3", Path2 = "#E7298A", Path3 = "#66A61E", Path1.Path2 = "#38B41E")
)
# Create the heatmap
heatmap <- pheatmap(cultureconditions, scale = "row", cluster_row = FALSE, cluster_cols = FALSE, annotation_row = annotation_row,
annotation_colors = ann_colors, angle_col = 45, fontsize = 16)
It seems that my rows for my "culturecondition" and annotation.rows are not matching as it provides me with this:
> dim(cultureconditions)
[1] 40 4
> dim(annotation_row)
[1] 43 1
I am not sure how to provide these additional row names given that it is the same row that I am trying annotate a second time.
Log in to answer this question.
When I run your code with a dummy dataset, it produces a figure as expected with Genepaths colored. There's something with your dataset. Could you post the first 20 rows of your dataset? An easy way to do this is to
deputlike this:dput(cultureconditions[1:20,]). If you can not post your dataset, then create a dummy one (see pheatmap docs for an example)Thank you for your response! I have attached the code results from the code below.
Appreciate your time in effort in helping me with this!
Thanks for posting the data. See my answer.