/hmdbQuery

utilities for querying Human Metabolome Database

Primary LanguageR

title author date vignette output
hmdbQuery: working with Human Metabolome Database (hmdb.ca)
Vincent J. Carey, stvjc at channing.harvard.edu
`r format(Sys.time(), '%B %d, %Y')`
%\VignetteEngine{knitr::rmarkdown} %\VignetteIndexEntry{hmdbQuery: working with Human Metabolome Database (hmdb.ca)}
html_document pdf_document
highlight number_sections theme toc
pygments
true
united
true
toc number_sections
true
true

Initial remarks

The human metabolomics database (HMDB, http://www.hmdb.ca) includes XML documents describing 114000 metabolites. We will show how to manipulate the metadata on metabolites fairly flexibly.

suppressMessages({
suppressPackageStartupMessages({
library(hmdbQuery)
library(gwascat)
})
})

Key utilities of the package

The hmdbQuery package includes a function for querying HMDB directly over HTTP:

library(hmdbQuery)
lk1 = HmdbEntry(prefix = "http://www.hmdb.ca/metabolites/", 
       id = "HMDB00001")

The result is parsed and encapsulated in an S4 object

lk1

The size of the complete import of information about a single metabolite suggests that it would not be too convenient to have comprehensive information about all HMDB constituents in memory. The most effective approach to managing the metadata will depend upon use cases to be developed over the long run.

Note however that this package does provide snapshots of certain direct associations derived from all available information as of Sept. 23 2017. Information about direct associations reported in the database is present in tables hmdb_disease, hmdb_gene, hmdb_protein, hmdb_omim. For example

data(hmdb_disease)
hmdb_disease

Working with the metadata

Disease associations

Some HMDB metabolites have been mapped to diseases.

diseases(lk1)
pmids = unlist(diseases(lk1)[1,]$pmids)
library(annotate)
pm = pubmed(pmids[1])
ab = buildPubMedAbst(xmlRoot(pm)[[1]])
ab

Biofluid and tissue location metadata

There are arbitrarily many biofluid and tissue associations provided for each HMDB entry. We have direct accessors, and by default we capture all metadata, available through the store method.

biofluids(lk1)
tissues(lk1)
st = store(lk1)
head(names(st))
length(names(st))
sapply(store(lk1)$protein_associations, "[[", "name")
sapply(store(lk1)$protein_associations, "[[", "gene_name")