Skip to contents

Disease Ontology over-representation analysis (ORA)

gly_enrich_do.Rd

[Deprecated]

This function was deprecated because we decided to move all enrichment analysis functions to the separate glyfun package, which has more features and better API design.

Perform Disease Ontology ORA for protein UniProt accessions using DOSE::enrichDO().

  • gly_enrich_do() accepts a glyexp::experiment() and extracts protein information from the "protein" column in the variable information tibble.

  • gly_enrich_do_() accepts a character vector of UniProt IDs.

As DOSE::enrichDO() only accepts Entrez IDs, the UniProt IDs will be first transformed into Entrez IDs with clusterProfiler::bitr().

Usage

gly_enrich_do(
  exp,
  add_info = TRUE,
  ont = "HDO",
  organism = "hsa",
  orgdb = "org.Hs.eg.db",
  universe = NULL,
  p_adj_method = "BH",
  p_cutoff = 0.05,
  q_cutoff = 0.2
)

gly_enrich_do_(
  proteins,
  ont = "HDO",
  organism = "hsa",
  orgdb = "org.Hs.eg.db",
  universe = NULL,
  p_adj_method = "BH",
  p_cutoff = 0.05,
  q_cutoff = 0.2
)

Arguments

exp

(Only for gly_enrich_do()) A glyexp::experiment() object.

add_info

A logical value. This parameter is included for API consistency but has no effect since enrichment results do not contain variable or sample columns. Only applicable to top-level APIs.

ont

Passed to ont of DOSE::enrichDO(). Disease Ontology type. Options: "HDO" (Human Disease Ontology), "HPO" (Human Phenotype Ontology), "MPO" (Mammalian Phenotype Ontology). Defaults to "HDO".

organism

Passed to organism of DOSE::enrichDO(). Organism code. "hsa" for human (Homo sapiens), "mmu" for mouse (Mus musculus). Defaults to "hsa".

orgdb

Passed to OrgDb of clusterProfiler::bitr(). Organism database name (e.g., "org.Hs.eg.db" for human). Defaults to "org.Hs.eg.db".

universe

Background genes. If a character vector, it is expected to contain UniProt accession IDs; these will be converted to Entrez Gene IDs and then passed to universe of DOSE::enrichDO(). You can also provide a glyexp::experiment() object with "glycoproteomics" type. In this case all detected proteins in this experiment will be extracted as UniProt IDs, converted to Entrez IDs, and then passed to universe of DOSE::enrichDO().

p_adj_method

Passed to pAdjustMethod of DOSE::enrichDO().

p_cutoff

Passed to pvalueCutoff of DOSE::enrichDO().

q_cutoff

Passed to qvalueCutoff of DOSE::enrichDO().

proteins

(Only for gly_enrich_do_()) A character vector of UniProt accession IDs.

Value

A list with three elements:

  • tidy_result: A tibble with enrichment results containing the following columns:

    • id: Term ID (e.g., DOID:XXXX)

    • description: Term description

    • gene_ratio: Ratio of genes in the term to total genes in the input

    • bg_ratio: Ratio of genes in the term to total genes in the background

    • p_val: Raw p-value from hypergeometric test

    • p_adj: Adjusted p-value

    • q_value: Q-value (FDR)

    • gene_id: Gene IDs in the term (separated by "/")

    • count: Number of genes in the term

  • raw_result: The raw DOSE enrichResult object

  • meta_data (only for gly_enrich_do()): A list containing metadata from the input experiment The list has classes glystats_do_ora_res and glystats_res.

Required packages

These functions require the following packages to be installed:

  • clusterProfiler for ID conversion

  • DOSE for enrichment analysis

  • org.Hs.eg.db for human gene annotation or other OrgDb packages

See also

DOSE::enrichDO()