Cargando…

Pathway size matters: the influence of pathway granularity on over-representation (enrichment analysis) statistics

BACKGROUND: Enrichment or over-representation analysis is a common method used in bioinformatics studies of transcriptomics, metabolomics, and microbiome datasets. The key idea behind enrichment analysis is: given a set of significantly expressed genes (or metabolites), use that set to infer a small...

Descripción completa

Detalles Bibliográficos
Autores principales: Karp, Peter D., Midford, Peter E., Caspi, Ron, Khodursky, Arkady
Formato: Online Artículo Texto
Lenguaje:English
Publicado: BioMed Central 2021
Materias:
Acceso en línea:https://www.ncbi.nlm.nih.gov/pmc/articles/PMC7967953/
https://www.ncbi.nlm.nih.gov/pubmed/33726670
http://dx.doi.org/10.1186/s12864-021-07502-8
_version_ 1783665976419549184
author Karp, Peter D.
Midford, Peter E.
Caspi, Ron
Khodursky, Arkady
author_facet Karp, Peter D.
Midford, Peter E.
Caspi, Ron
Khodursky, Arkady
author_sort Karp, Peter D.
collection PubMed
description BACKGROUND: Enrichment or over-representation analysis is a common method used in bioinformatics studies of transcriptomics, metabolomics, and microbiome datasets. The key idea behind enrichment analysis is: given a set of significantly expressed genes (or metabolites), use that set to infer a smaller set of perturbed biological pathways or processes, in which those genes (or metabolites) play a role. Enrichment computations rely on collections of defined biological pathways and/or processes, which are usually drawn from pathway databases. Although practitioners of enrichment analysis take great care to employ statistical corrections (e.g., for multiple testing), they appear unaware that enrichment results are quite sensitive to the pathway definitions that the calculation uses. RESULTS: We show that alternative pathway definitions can alter enrichment p-values by up to nine orders of magnitude, whereas statistical corrections typically alter enrichment p-values by only two orders of magnitude. We present multiple examples where the smaller pathway definitions used in the EcoCyc database produces stronger enrichment p-values than the much larger pathway definitions used in the KEGG database; we demonstrate that to attain a given enrichment p-value, KEGG-based enrichment analyses require 1.3–2.0 times as many significantly expressed genes as does EcoCyc-based enrichment analyses. The large pathways in KEGG are problematic for another reason: they blur together multiple (as many as 21) biological processes. When such a KEGG pathway receives a high enrichment p-value, which of its component processes is perturbed is unclear, and thus the biological conclusions drawn from enrichment of large pathways are also in question. CONCLUSIONS: The choice of pathway database used in enrichment analyses can have a much stronger effect on the enrichment results than the statistical corrections used in these analyses. SUPPLEMENTARY INFORMATION: The online version contains supplementary material available at (10.1186/s12864-021-07502-8).
format Online
Article
Text
id pubmed-7967953
institution National Center for Biotechnology Information
language English
publishDate 2021
publisher BioMed Central
record_format MEDLINE/PubMed
spelling pubmed-79679532021-03-22 Pathway size matters: the influence of pathway granularity on over-representation (enrichment analysis) statistics Karp, Peter D. Midford, Peter E. Caspi, Ron Khodursky, Arkady BMC Genomics Research Article BACKGROUND: Enrichment or over-representation analysis is a common method used in bioinformatics studies of transcriptomics, metabolomics, and microbiome datasets. The key idea behind enrichment analysis is: given a set of significantly expressed genes (or metabolites), use that set to infer a smaller set of perturbed biological pathways or processes, in which those genes (or metabolites) play a role. Enrichment computations rely on collections of defined biological pathways and/or processes, which are usually drawn from pathway databases. Although practitioners of enrichment analysis take great care to employ statistical corrections (e.g., for multiple testing), they appear unaware that enrichment results are quite sensitive to the pathway definitions that the calculation uses. RESULTS: We show that alternative pathway definitions can alter enrichment p-values by up to nine orders of magnitude, whereas statistical corrections typically alter enrichment p-values by only two orders of magnitude. We present multiple examples where the smaller pathway definitions used in the EcoCyc database produces stronger enrichment p-values than the much larger pathway definitions used in the KEGG database; we demonstrate that to attain a given enrichment p-value, KEGG-based enrichment analyses require 1.3–2.0 times as many significantly expressed genes as does EcoCyc-based enrichment analyses. The large pathways in KEGG are problematic for another reason: they blur together multiple (as many as 21) biological processes. When such a KEGG pathway receives a high enrichment p-value, which of its component processes is perturbed is unclear, and thus the biological conclusions drawn from enrichment of large pathways are also in question. CONCLUSIONS: The choice of pathway database used in enrichment analyses can have a much stronger effect on the enrichment results than the statistical corrections used in these analyses. SUPPLEMENTARY INFORMATION: The online version contains supplementary material available at (10.1186/s12864-021-07502-8). BioMed Central 2021-03-16 /pmc/articles/PMC7967953/ /pubmed/33726670 http://dx.doi.org/10.1186/s12864-021-07502-8 Text en © The Author(s) 2021 Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
spellingShingle Research Article
Karp, Peter D.
Midford, Peter E.
Caspi, Ron
Khodursky, Arkady
Pathway size matters: the influence of pathway granularity on over-representation (enrichment analysis) statistics
title Pathway size matters: the influence of pathway granularity on over-representation (enrichment analysis) statistics
title_full Pathway size matters: the influence of pathway granularity on over-representation (enrichment analysis) statistics
title_fullStr Pathway size matters: the influence of pathway granularity on over-representation (enrichment analysis) statistics
title_full_unstemmed Pathway size matters: the influence of pathway granularity on over-representation (enrichment analysis) statistics
title_short Pathway size matters: the influence of pathway granularity on over-representation (enrichment analysis) statistics
title_sort pathway size matters: the influence of pathway granularity on over-representation (enrichment analysis) statistics
topic Research Article
url https://www.ncbi.nlm.nih.gov/pmc/articles/PMC7967953/
https://www.ncbi.nlm.nih.gov/pubmed/33726670
http://dx.doi.org/10.1186/s12864-021-07502-8
work_keys_str_mv AT karppeterd pathwaysizematterstheinfluenceofpathwaygranularityonoverrepresentationenrichmentanalysisstatistics
AT midfordpetere pathwaysizematterstheinfluenceofpathwaygranularityonoverrepresentationenrichmentanalysisstatistics
AT caspiron pathwaysizematterstheinfluenceofpathwaygranularityonoverrepresentationenrichmentanalysisstatistics
AT khodurskyarkady pathwaysizematterstheinfluenceofpathwaygranularityonoverrepresentationenrichmentanalysisstatistics