230 likes | 243 Views
This overview introduces the concept of clustering, including its mechanics, parameter choices, and examples. It discusses gene and sample clustering, as well as different clustering methods and their effects. The text is in English.
E N D
An Overview of Clustering Methods Michael D. Kane, Ph.D.
Topics • What is clustering? • Clustering mechanics (how the computer does it). • Parameter choices and their effect. • Examples.
What is clustering? Grouping by similarity.
Samples Gene clustering Similar genes.Group genes that have similar expression profiles when observed over multiple samples. Genes
Samples Sample clustering Similar samples.Group samples that are similar when observed over multiple genes. Genes
Why cluster? • Similar gene expression infers common biology. Function of uncharacterized genes may be deduced from co- expression with known genes. • Associate expression patterns with: Response to environmental change. Disease pathology/progression.
+ E2 E1 E2 E2 c Gene a Gene b a e d - + Gene c E1 b Gene d Gene e f Gene f - Clustering Mechanics For gene clustering, we must measure similarity between genes.
+ E2 c a e d - + E1 b f - Distance (similarity) measure Euclidean distance (1.0, 1.7) dbe (4.6, 0.5)
Distance Measure Pearson Correlation S=(-1 . . . +1) Used in “Eisen” clustering
+ E2 c a e d - + E1 b f - Hierarchical Clustering a b c d e f
Measuring distance between clusters Single linkage The minimum distance between clusters. May form loose clusters. Produces “chained” clusters. Complete linkage The maximum distance between clusters. Tends to form compact clusters.
Methods for joining clusters UPGMA unweighted pair group method (Average linkage) The average distance between clusters. Weighted pair group method Same as UPGMA but the distance is weighted by cluster size. Use when clusters are expected to be significantly uneven in size!
Effect of distance measure Euclidean Single Linkage Euclidean Complete Linkage
Effect of distance measure Euclidean UPGMA Euclidean Ward’s Method
k-means Alternatives to hierarchical clustering • Number of clusters specified by user. • Good when prior knowledge available.
1. Number of clusters specified by user. 2. Genes randomly assigned to clusters. + E2 c a e d - + E1 3. Assess inter and intra-cluster similarity. 3. Assess inter and intra-cluster similarity. b f 4. Move genes to alternative cluster if distance is reduced. 4. Move genes to alternative cluster if distance is reduced. - k-means clustering
SOM Self-organizing maps Alternatives to hierarchical clustering • Number of clusters specified by user. • Good when prior knowledge available.
E1 E2 E2 E1 E1 E1 E1 E1 E2 E2 E2 E2 E2 + 0 - + 0 - + 0 - + 0 - + 0 - cluster 1 cluster 1 cluster 1 cluster 1 cluster 1 E1 E2 + 0 - Gene a E1 E1 E1 E1 E1 E2 E2 E2 E2 E2 + 0 - + 0 - + 0 - + 0 - + 0 - cluster 2 cluster 2 cluster 2 cluster 2 cluster 2 + 0 - Gene b + 0 - Gene c E1 E1 E1 E1 E1 E2 E2 E2 E2 E2 + 0 - + 0 - + 0 - + 0 - + 0 - + 0 - Gene d cluster 3 cluster 3 cluster 3 cluster 3 cluster 3 + 0 - Gene e + 0 - After training, assign each gene to the most similar cluster. Increase the similarity by adjusting the cluster representation. “Training” For a gene, find the most similar cluster representation. User specified number of clusters. Each initially given a random expression representation. Iteratively train the cluster representations. Gene f SOM
Gene clustering Eisen et al., Cluster analysis and display of genome-wide expression patterns. PNAS v95,14863-14868, 1998 cholesterol biosynthesis 24 hour time course after re-introduction of serum to serum-deprived human fibroblasts. Pearson correlation, average linkage. cell cycle immediate-early response signaling wound healing
Sample clustering Ross et al., Systematic variation in gene expression patterns in human cancer cell lines. Nature Genetics v24, 227-235, 2000 Note breast cancer cell lines, derived from the same patient. 64 cancer cell lines clustered. 8,000 genes. Clustering performed with 2 different subsets of genes. Similar results. Pearson correlation, average linkage.
Summary • Different methods often provide different clusters. • No overall “best” clustering method. • Clustering applied to unrelated data will still provide clusters. • Use biological insight in method selection and interpretation.
+ E2 c a e d - + E1 b f - Clustering a b c d e f
E1 E2 + 0 - cluster 1 E1 E2 + 0 - Gene a E1 E2 + 0 - cluster 2 + 0 - Gene b + 0 - Gene c E1 E2 + 0 - + 0 - Gene d cluster 3 + 0 - Gene e + 0 - After training, assign each gene to the most similar cluster. Gene f SOM E2