Cluster analysis Cluster analysis Function Function Places genes with similar expression patterns in groups. Sometimes genes of unknown function will be grouped with genes of known function. The functions that are known allow the investigator to hypothesize regarding the functions of genes not yet characterized. Examples: Identify genes important in cell cycle regulation Identify genes that participate in a biosynthetic pathway Identify genes involved in a drug response Identify genes involved in a disease response
42
Embed
Cluster analysis Function Places genes with similar expression patterns in groups. Sometimes genes of unknown function will be grouped with genes.
This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Transcript
Cluster analysis Cluster analysis Function Function Places genes with similar expression patterns in
groups. Sometimes genes of unknown function will be
grouped with genes of known function. The functions that are known allow the investigator
to hypothesize regarding the functions of genes not yet characterized.
Examples: Identify genes important in cell cycle regulation Identify genes that participate in a biosynthetic pathway Identify genes involved in a drug response Identify genes involved in a disease response
Cluster analysis software is developed to group genes with similar patterns of expression.
In this example, the columns represent different timepoints, and the rows represent the results for a single gene.
The products of the genes expressed in a single cluster may have related or similar functions.
FUNCTIONAL GENOMICSCLUSTER ANALYSIS OF MICROARRAY DATA
(11)
CLUSTER ANALYSISCLUSTER ANALYSIS
OUTLINE OF TALK
Clustering - GoalClustering - Goal
Partition of the genes in the dataset into distinct sets (clusters), according to similarity in their expression profiles across the probed conditions
Euclidian distance – takes into account the magnitude of the expression
Correlation distance - insensitive to the amplitude of expression, takes into account the trends of the change.
Common trends are considered biologically relevant, the magnitude is considered less important → correlation
Gene X
Gene Y
What correlation distance sees
What euclidean distance sees
PCAPCA
A technique for projecting the expression data set onto a reduced (2 or 3 dimensional) easily visualized space
Dataset: Thousands of genes probed in 10 conditions. The expression profile of each gene is presented by the
vector of its expression levels: X = (X1, X2, X3, X4, X5) Imagine each gene X as a point in a 5-dimentional
space. Each direction/axis corresponds to a specific condition Genes with similar profiles are close to each other in
this space PCA- Project this dataset to 2 dimensions, preserving
as much information as possible
PCA transformation of a microarray PCA transformation of a microarray datasetdataset
Visual estimation of the number of clusters in the data
Clustering AlgorithmsClustering Algorithms
K–meansSOMsHierarchical clustering
K-MEANSK-MEANS1. The user sets the number of clusters- k2. Initialization: each gene is randomly assigned
to one of the k clusters3. Average expression vector is calculated for
each cluster (cluster’s profile) 4. Iterate over the genes:
• For each gene- compute its similarity to the cluster profiles.
• Move the gene to the cluster it is most similar to.• Recalculated cluster profiles.
5. Score current partition: sum of distances between genes and the profile of the cluster they are assigned to (homogeneity of the solution).
6. Stop criteria: further shuffling of genes results in minor improvement in the clustering score
How to choose the number of clusters How to choose the number of clusters needed to informatively partition the data needed to informatively partition the data
Try several parameters (number of desired clusters,distance metric) and compare the clustering
User sets the number of clusters in a form of a rectangular grid (e.g., 3x2) – ‘map nodesmap nodes’
Imagine genes as points in (M-dimensional) space
Initialization: map nodes are randomly placed in the data space
Genes – data points
Clusters – map nodes
SOM - SchemeSOM - Scheme
• Randomly choose a data point (gene).
• Find its closest map node
• Move this map node towards the data point
• Move the neighbor map nodes towards this point, but to lesser extent (thinner arrows show weaker shift)
• Iterate over data points
• Each successive gene profile (black dot) has less of an influence on the displacement of the nodes.
• Iterate through all profiles several times (10-100)
• When positions of the cluster nodes have stabilized, assign each gene to its closest map node (cluster)
Hierarchical ClusteringHierarchical Clustering Goal#1: Organize the genes in a
structure of a hierarchical tree 1) Initial step: each gene is
regarded as a cluster with one item 2) Find the 2 most similar clusters
and merge them into a common node (red dot)
3) Merge successive nodes until all genes are contained in a single cluster
Goal#2: Collapse branches to group genes into distinct clusters g1 g2 g3 g4 g5
{1,2}
{4,5}
{1,2,3}
{1,2,3,4,5}
Mathematical evaluation of Mathematical evaluation of clustering solutionclustering solution
Merits of a ‘good’ clustering solution: Homogeneity:
Genes inside a cluster are highly similar to each other. Average similarity between a gene and the center
(average profile) of its cluster.
Separation: Genes from different clusters have low similarity to each
other. Weighted average similarity between centers of clusters.
These are conflicting features: increasing the number of clusters tends to improve with-in cluster Homogeneity on the expense of between-cluster Separation
“True”CAST*
GeneCluster
K-means
CLICK
Homogeneity
Separa
tion
Performance on Yeast Cell Cycle Data
*Ben-Dor, Shamir, Yakhini
1999
698 genes, 72 conditions (Spellman et al. 1998). Each algorithm was run by its authors in a “blind” test.
Overall strategy:
PCA-transformationClustering and evaluation of clusteringCheck for bio-significance
Which genes to cluster? Which genes to cluster? Apply filtering prior to clustering – focus the
analysis on the ‘responding genes’ Applying controlled statistical tests to identify
‘responding genes’ usually ends up with too few genes that do not allow for a global characterization of the response.
Fold change: choose genes that changed by at least M-folds in at least L conditions
Variance: filter out genes that do not vary greatly among the conditions of the experiment.
Try various filtering schemes to find the setting that gives the best biological results
Ascribing Biological Meaning to Ascribing Biological Meaning to ClustersClusters
Identify over-represented functional categories in the clusters (i.e., cluster contains much more genes of specific biological process than expected by chance)
Requirements for systematic analysis: Controlled vocabulary for describing biological