I have a RNAseq dataset with genes as rows and samples as columns. I have managed to do hierarchical clustering of the samples, and I want to look for subgroups in the disease. I additionally have metadata with sample names as one column and lesion type in the other column. I want to add the metadata to the hierarchical cluster dendrogram so that I can display the "lesional", "control", or "non lesional" category of the samples. How do I do that in R?