Cited by 2,235Bioinformatics and Genomic NetworksBiomedical Text Mining and OntologiesGenomics and Phylogenetic Studies

The DOI matches its Crossref record; 58 of the 58 references it cites were found in a registry or an index.

Enrichment Map: A Network-Based Method for Gene-Set Enrichment Visualization and Interpretation

Open access

Abstract

BACKGROUND: Gene-set enrichment analysis is a useful technique to help functionally characterize large gene lists, such as the results of gene expression experiments. This technique finds functionally coherent gene-sets, such as pathways, that are statistically over-represented in a given gene list. Ideally, the number of resulting sets is smaller than the number of genes in the list, thus simplifying interpretation. However, the increasing number and redundancy of gene-sets used by many current enrichment analysis software works against this ideal. PRINCIPAL FINDINGS: To overcome gene-set redundancy and help in the interpretation of large gene lists, we developed "Enrichment Map", a network-based visualization method for gene-set enrichment results. Gene-sets are organized in a network, where each set is a node and edges represent gene overlap between sets. Automated network layout groups related gene-sets into network clusters, enabling the user to quickly identify the major enriched functional themes and more easily interpret the enrichment results. CONCLUSIONS: Enrichment Map is a significant advance in the interpretation of enrichment analysis. Any research project that generates a list of genes can take advantage of this visualization framework. Enrichment Map is implemented as a freely available and user friendly plug-in for the Cytoscape network visualization software (http://baderlab.org/Software/EnrichmentMap/).

Overview

The claims and statistics of this paper haven’t been checked yet.

Run the check — free, no account needed. It reads the open-access text and checks the sources it cites.

A few free checks a day · papers already verified open instantly · results are public

Integrity signals

Registry

Crossref · 10.1371/journal.pone.0013984

Cites

58 of 58 references verified0 citation errors · 0 not checkable

Statistics

Figures

Figure 1 From ranked gene lists to the enrichment map. High-throughput genomic experiments often output large gene lists, which are typically ranked according to a statistic measuring difference in one experimental condition versus another. Ranked lists are analyzed for enrichment in known sets of functionally related genes (e.g. pathways) from publicly accessible databases. An enrichment map is drawn, representing the enrichment results as a network of gene-sets (nodes) related by their similarity (edges), with enrichment significance encoded by the node color gradient, where color intensity represents significance and color hue (red / blue) represents the class (i.e. biological condition) of interest. Node size represents the gene-set size and edge thickness represents the degree of overlap between two gene-sets.
Figure 1 From ranked gene lists to the enrichment map. High-throughput genomic experiments often output large gene lists, which are typically ranked according to a statistic measuring difference in one experimental condition versus another. Ranked lists are analyzed for enrichment in known sets of functionally related genes (e.g. pathways) from publicly accessible databases. An enrichment map is drawn, representing the enrichment results as a network of gene-sets (nodes) related by their similarity (edges), with enrichment significance encoded by the node color gradient, where color intensity represents significance and color hue (red / blue) represents the class (i.e. biological condition) of interest. Node size represents the gene-set size and edge thickness represents the degree of overlap between two gene-sets.CC BY · PubMed Central
Figure 2 Hierarchical visualization of enrichment results for estrogen treatment of breast cancer cells. Hierarchical organization of GO gene-set enrichment results for estrogen-treated compared to untreated breast cancer cells at 24 hours of culture. Nodes represent gene-sets and edges represent GO defined relations (Is-a, Part-of, Regulates). Gene-sets that did not pass the enrichment significance threshold are not shown. Nodes are colored according to enrichment results: red represents enrichment in estrogen-treated cells (i.e. up-regulation after estrogen treatment), whereas blue represents enrichment in untreated cells (i.e. down-regulation after estrogen treatment). Color intensity is proportional to enrichment significance. Since conservative thresholds were used to select gene-sets, most of the node colors are intense (corresponding to highly significant gene-sets). Subnetworks (i.e. clusters) are annotated according to the corresponding function. The acronym in brackets represents the specific GO ontology the gene-sets belong to: Molecular Function (MF), Cellular Component (CC), Biological Process (BP). Microtubule cytoskeleton (purple labels) and tRNA processing (green labels) were highlighted to show absence of connections between related gene-sets.
Figure 2 Hierarchical visualization of enrichment results for estrogen treatment of breast cancer cells. Hierarchical organization of GO gene-set enrichment results for estrogen-treated compared to untreated breast cancer cells at 24 hours of culture. Nodes represent gene-sets and edges represent GO defined relations (Is-a, Part-of, Regulates). Gene-sets that did not pass the enrichment significance threshold are not shown. Nodes are colored according to enrichment results: red represents enrichment in estrogen-treated cells (i.e. up-regulation after estrogen treatment), whereas blue represents enrichment in untreated cells (i.e. down-regulation after estrogen treatment). Color intensity is proportional to enrichment significance. Since conservative thresholds were used to select gene-sets, most of the node colors are intense (corresponding to highly significant gene-sets). Subnetworks (i.e. clusters) are annotated according to the corresponding function. The acronym in brackets represents the specific GO ontology the gene-sets belong to: Molecular Function (MF), Cellular Component (CC), Biological Process (BP). Microtubule cytoskeleton (purple labels) and tRNA processing (green labels) were highlighted to show absence of connections between related gene-sets.CC BY · PubMed Central
Figure 3 Enrichment map for estrogen treatment of breast cancer cells at 24 hours of culture. The map displays the enriched gene-sets in estrogen-treated vs. untreated breast cancer cells at 24 hours of culture. As in Figure 2, red node color represents enrichment in estrogen-treated cells (i.e. up-regulation after estrogen treatment), whereas blue represents enrichment in untreated cells (i.e. down-regulation after estrogen treatment); color intensity is proportional to enrichment significance. Clusters of functionally related gene-sets were manually circled and assigned a label.
Figure 3 Enrichment map for estrogen treatment of breast cancer cells at 24 hours of culture. The map displays the enriched gene-sets in estrogen-treated vs. untreated breast cancer cells at 24 hours of culture. As in Figure 2, red node color represents enrichment in estrogen-treated cells (i.e. up-regulation after estrogen treatment), whereas blue represents enrichment in untreated cells (i.e. down-regulation after estrogen treatment); color intensity is proportional to enrichment significance. Clusters of functionally related gene-sets were manually circled and assigned a label.CC BY · PubMed Central
Figure 4 Zoom in of the microtubule cytoskeleton cluster in the 24 hours estrogen treatment enrichment map. Cytoskeleton-related gene-sets from different GO partitions, such as Spindle (CC) and Microtubule-based process (BP) have been grouped together, unlike in the purely hierarchical visualization in Figure 2. As in the previous figures, red node color represents enrichment in estrogen-treated cells (i.e. up-regulation after estrogen treatment), whereas blue represents enrichment in untreated cells (i.e. down-regulation after estrogen treatment); color intensity is proportional to enrichment significance.
Figure 4 Zoom in of the microtubule cytoskeleton cluster in the 24 hours estrogen treatment enrichment map. Cytoskeleton-related gene-sets from different GO partitions, such as Spindle (CC) and Microtubule-based process (BP) have been grouped together, unlike in the purely hierarchical visualization in Figure 2. As in the previous figures, red node color represents enrichment in estrogen-treated cells (i.e. up-regulation after estrogen treatment), whereas blue represents enrichment in untreated cells (i.e. down-regulation after estrogen treatment); color intensity is proportional to enrichment significance.CC BY · PubMed Central
Figure 5 Enrichment map for estrogen treatment of breast cancer cells at 12 and 24 hours of culture. The map displays the enriched gene-sets in estrogen-treated vs. untreated breast cancer cells at 12 and 24 hours of culture. Enrichments were mapped to the inner node area and to the node borders, respectively. As in the previous figures, red represents enrichment in estrogen-treated cells (i.e. up-regulation after estrogen treatment), whereas blue represents enrichment in untreated cells (i.e. down-regulation after estrogen treatment); color intensity is proportional to enrichment significance. Clusters of functionally related gene-sets were manually circled and assigned a label.
Figure 5 Enrichment map for estrogen treatment of breast cancer cells at 12 and 24 hours of culture. The map displays the enriched gene-sets in estrogen-treated vs. untreated breast cancer cells at 12 and 24 hours of culture. Enrichments were mapped to the inner node area and to the node borders, respectively. As in the previous figures, red represents enrichment in estrogen-treated cells (i.e. up-regulation after estrogen treatment), whereas blue represents enrichment in untreated cells (i.e. down-regulation after estrogen treatment); color intensity is proportional to enrichment significance. Clusters of functionally related gene-sets were manually circled and assigned a label.CC BY · PubMed Central
Figure 6 Heat-maps displaying gene-set expression patterns in the estrogen treatment experiment. Two gene-sets displaying different enrichment patterns at 12 and 24 hours of the estrogen treatment experiment were selected from the enrichment map in Figure 5 and their expression patterns were explored using heat maps within the Enrichment Map software. For APC-dependent protein degradation (left), there is a marked increase of gene expression in estrogen treated cells at 24 hours compared to 12 hours, whereas the gene levels for untreated cells are substantially the same at the two time points; the pattern observed is consistent with the presence of significant enrichment only at 24 hours. On the other hand, for Replication fork (right), gene expression in estrogen treated cells at 12 and 24 hours is globally at similar levels, whereas there is an increase of gene levels in untreated cells. This suggests an explanation of why Replication fork is enriched only at 12 hours.
Figure 6 Heat-maps displaying gene-set expression patterns in the estrogen treatment experiment. Two gene-sets displaying different enrichment patterns at 12 and 24 hours of the estrogen treatment experiment were selected from the enrichment map in Figure 5 and their expression patterns were explored using heat maps within the Enrichment Map software. For APC-dependent protein degradation (left), there is a marked increase of gene expression in estrogen treated cells at 24 hours compared to 12 hours, whereas the gene levels for untreated cells are substantially the same at the two time points; the pattern observed is consistent with the presence of significant enrichment only at 24 hours. On the other hand, for Replication fork (right), gene expression in estrogen treated cells at 12 and 24 hours is globally at similar levels, whereas there is an increase of gene levels in untreated cells. This suggests an explanation of why Replication fork is enriched only at 12 hours.CC BY · PubMed Central
Figure 7 Enrichment map for early-onset colon cancer and overlap with known disease genes. The map displays the enriched gene-sets in early onset colon cancer patients vs. normal controls. The yellow triangle represents the set of known colon cancer genes from the DiseaseHub database, which integrates disease gene lists from several genotype-phenotype association resources. Purple edges indicate overlap between the disease signature and enriched gene-sets; thickness represents significance. Only edges with a Fisher's Exact Test nominal p-value smaller than 10−4 were visualized.
Figure 7 Enrichment map for early-onset colon cancer and overlap with known disease genes. The map displays the enriched gene-sets in early onset colon cancer patients vs. normal controls. The yellow triangle represents the set of known colon cancer genes from the DiseaseHub database, which integrates disease gene lists from several genotype-phenotype association resources. Purple edges indicate overlap between the disease signature and enriched gene-sets; thickness represents significance. Only edges with a Fisher's Exact Test nominal p-value smaller than 10−4 were visualized.CC BY · PubMed Central

Verification

Claims the paper makes

The claims and statistics of this paper haven’t been checked yet.

Run the check — free, no account needed. It reads the open-access text and checks the sources it cites.

A few free checks a day · papers already verified open instantly · results are public

Claims other papers make about this one

Each is checked against this paper’s own text: does it say what the citing paper says it does?

Self-check

Percentages, sums and reported differences re-computed from the paper’s own numbers.

The self-check needs the paper’s full text.

Citations

Cites

Show all 58 references

Research works printed but not linked: 1

  1. S Draghici, 2003, Global functional profiling of gene expression., Genomics, 81, 98

Cited by

Sorted by citations

Items 1,001–1,005 of 1,005

Metadata

Authors
Merico, D., Isserlin, R., Stueker, O., Emili, A., & Bader, G. D.
Affiliations
University of Toronto (CA)
Corresponding author
  • Merico, Daniele
  • Bader, Gary D.
Venue
PLoS ONE
Published
2010-11-15
Volume
5
Issue
11
Pages
e13984
ISSN
1932-6203 (ISSN-L)
OpenAlex
W2097948046
Cited by
2,235
References
58
Citations by year
2026: 89 · 2025: 145 · 2024: 176 · 2023: 216 · 2022: 238 · 2021: 233 · 2020: 237 · 2019: 164 · 2018: 158 · 2017: 145 · 2016: 127 · 2015: 101 · 2014: 79 · 2013: 69 · 2012: 43
FWCI
16.30
Type
journal article
Language
en
Open access
Locations
  • PLoS ONE (publishedVersion, cc-by)
  • PubMed Central (publishedVersion, cc-by)
  • PubMed (publishedVersion)
  • DOAJ (DOAJ: Directory of Open Access Journals) (submittedVersion, cc-by-sa)
  • Figshare (submittedVersion, cc-by)
Funders
  • Heart and Stroke Foundation of Canada
  • Genome Canada
  • Ontario Genomics
  • Ontario Genomics Institute
  • Hospital for Sick Children
  • National Institutes of Health
  • National Human Genome Research Institute (P41 P41HG04118, P41HG04118, P41HG004118)
Topics
  • Life Sciences › Biochemistry, Genetics and Molecular Biology › Molecular Biology › Bioinformatics and Genomic Networks
  • Health Sciences › Medicine › Health Informatics › Biomedical Text Mining and Ontologies
  • Life Sciences › Biochemistry, Genetics and Molecular Biology › Molecular Biology › Genomics and Phylogenetic Studies
Keywords
gene set enrichment analysis, network visualization, gene set overlap, network clustering, pathway over-representation, gene expression data interpretation, Cytoscape plugin
MeSH terms
Algorithms, Breast Neoplasms, Colonic Neoplasms, Estrogens, Female, Humans, Software, Reproducibility of Results, Gene Expression Regulation, Neoplastic, Cluster Analysis, Computational Biology, Internet, Gene Expression Profiling, Gene Regulatory Networks

About these checks

Registry, retraction, Cites and statistics results are computed, not judged, with one exception: a printed reference that no index lists is read once by a model before it counts as a citation error. Claims are extracted by a model and kept only when their quote appears verbatim in the text; each cited source is then checked, with a model’s help, for whether it says what the claim relies on. Those verdicts can be wrong — read the claim’s quote and the cited source before relying on one.

Statistics and arithmetic are recomputed, not judged.

Confidence bands are shown only for fields calibrated on our annotated test set. None is calibrated yet, so no band is shown.

What counts as a claim ›How we measure confidence ›How we check this ›

Are you an author, or do you see a problem with this page? Request a review