| Title: | R-side access to BioPlex protein-protein interaction data |
|---|---|
| Description: | The BioPlex package implements access to the BioPlex protein-protein interaction networks and related resources from within R. Besides protein-protein interaction networks for HEK293 and HCT116 cells, this includes access to CORUM protein complex data, and transcriptome and proteome data for the two cell lines. Functionality focuses on importing the various data resources and storing them in dedicated Bioconductor data structures, as a foundation for integrative downstream analysis of the data. |
| Authors: | Ludwig Geistlinger [aut, cre], Robert Gentleman [aut] |
| Maintainer: | Ludwig Geistlinger <[email protected]> |
| License: | Artistic-2.0 |
| Version: | 1.18.0 |
| Built: | 2026-10-10 11:18:33 UTC |
| Source: | https://github.com/bioc/BioPlex |
This function adds PFAM domain annotations to the node metadata of the BioPlex PPI graph.
annotatePFAM(bp.gr, orgdb)annotatePFAM(bp.gr, orgdb)
bp.gr |
an object of class |
orgdb |
an |
An object of class graphNEL containing PFAM domain annotations
in the nodeData.
BioPlex: https://bioplex.hms.harvard.edu/interactions.php
PFAM: http://pfam.xfam.org
# (1) Obtain the latest version of the 293T PPI network bp.293t <- getBioPlex(cell.line = "293T", version = "3.0") # (2) Turn the data into a graph bp.gr <- bioplex2graph(bp.293t) # (3) Obtain orgdb package from AnnotationHub ah <- AnnotationHub::AnnotationHub() orgdb <- AnnotationHub::query(ah, c("orgDb", "Homo sapiens")) orgdb <- orgdb[[1]] # (4) Annotate PFAM domains bp.gr <- annotatePFAM(bp.gr, orgdb)# (1) Obtain the latest version of the 293T PPI network bp.293t <- getBioPlex(cell.line = "293T", version = "3.0") # (2) Turn the data into a graph bp.gr <- bioplex2graph(bp.293t) # (3) Obtain orgdb package from AnnotationHub ah <- AnnotationHub::AnnotationHub() orgdb <- AnnotationHub::query(ah, c("orgDb", "Homo sapiens")) orgdb <- orgdb[[1]] # (4) Annotate PFAM domains bp.gr <- annotatePFAM(bp.gr, orgdb)
Representation of BioPlex PPIs in a graphNEL object
from the graph package.
bioplex2graph(bioplex.df)bioplex2graph(bioplex.df)
bioplex.df |
a |
An object of class graphNEL.
BioPlex: https://bioplex.hms.harvard.edu/interactions.php
# (1) Obtain the latest version of the 293T PPI network bp.293t <- getBioPlex(cell.line = "293T", version = "3.0") # (2) Turn the data into a graph bp.gr <- bioplex2graph(bp.293t)# (1) Obtain the latest version of the 293T PPI network bp.293t <- getBioPlex(cell.line = "293T", version = "3.0") # (2) Turn the data into a graph bp.gr <- bioplex2graph(bp.293t)
Functionality for storing the protein expression data
from the Cancer Cell Line Encyclopedia (CCLE) in a
SummarizedExperiment.
ccleProteome2SummarizedExperiment(df, cell.line = "HCT116")ccleProteome2SummarizedExperiment(df, cell.line = "HCT116")
df |
a |
cell.line |
character. One or more cell line IDs such as |
A SummarizedExperiment storing protein expression
data for the specified cell line(s).
CCLE proteomics: https://gygi.hms.harvard.edu/publications/ccle.html
# Connect to ExperimentHub eh <- ExperimentHub::ExperimentHub() # Obtain CCLE proteome data frame AnnotationHub::query(eh, c("gygi", "depmap")) ccle.prot <- eh[["EH3459"]] ccle.prot <- as.data.frame(ccle.prot) # Turn into a SummarizedExperiment se <- ccleProteome2SummarizedExperiment(ccle.prot)# Connect to ExperimentHub eh <- ExperimentHub::ExperimentHub() # Obtain CCLE proteome data frame AnnotationHub::query(eh, c("gygi", "depmap")) ccle.prot <- eh[["EH3459"]] ccle.prot <- as.data.frame(ccle.prot) # Turn into a SummarizedExperiment se <- ccleProteome2SummarizedExperiment(ccle.prot)
Functionality for storing CORUM protein complex data in a
list of graph instances.
corum2graphlist(corum.df, subunit.id.type = c("UNIPROT", "ENTREZID"))corum2graphlist(corum.df, subunit.id.type = c("UNIPROT", "ENTREZID"))
corum.df |
A |
subunit.id.type |
character. Supported options include |
A list with an entry for each complex. Each entry is an
object of class graphNEL connecting all subunit IDs with each other
by undirected edges.
CORUM: http://mips.helmholtz-muenchen.de/corum/#download
# (1) Obtain the core set of CORUM complexes ... core <- getCorum(set = "core") # (2) ... turn into a list of graphs core.glist <- corum2graphlist(core)# (1) Obtain the core set of CORUM complexes ... core <- getCorum(set = "core") # (2) ... turn into a list of graphs core.glist <- corum2graphlist(core)
Functionality for storing CORUM protein complex data in a
list.
corum2list(corum.df, subunit.id.type = c("UNIPROT", "ENTREZID"))corum2list(corum.df, subunit.id.type = c("UNIPROT", "ENTREZID"))
corum.df |
A |
subunit.id.type |
character. Supported options include |
A list with an entry for each complex. Each entry is a
character vector of subunit IDs.
CORUM: http://mips.helmholtz-muenchen.de/corum/#download
# (1) Obtain the core set of CORUM complexes ... core <- getCorum(set = "core") # (2) ... turn into a list core.list <- corum2list(core)# (1) Obtain the core set of CORUM complexes ... core <- getCorum(set = "core") # (2) ... turn into a list core.list <- corum2list(core)
Functionality for retrieving the BioPlex protein-protein interaction data. Available networks include:
BioPlex 293T cells (versions 1.0, 2.0, and 3.0),
BioPlex HCT116 cells (version 1.0).
See references.
getBioPlex( cell.line = c("293T", "HCT116"), version = c("3.0", "1.0", "2.0"), remap.uniprot.ids = FALSE, cache = TRUE )getBioPlex( cell.line = c("293T", "HCT116"), version = c("3.0", "1.0", "2.0"), remap.uniprot.ids = FALSE, cache = TRUE )
cell.line |
character. Valid options include:
Defaults to |
version |
character. Valid options include |
remap.uniprot.ids |
logical. Should the protein-to-gene mappings from BioPlex
(i.e. UNIPROT-to-SYMBOL and UNIPROT-to-ENTREZID) be updated using Bioc annotation
functionality?
Defaults to |
cache |
logical. Should a locally cached version used if available?
Defaults to |
A data.frame.
BioPlex: https://bioplex.hms.harvard.edu/interactions.php
# (1) Obtain the latest version of the 293T PPI network bp.293t <- getBioPlex(cell.line = "293T", version = "3.0") # (2) Obtain the latest version of the HCT116 PPI network bp.hct116 <- getBioPlex(cell.line = "HCT116", version = "1.0")# (1) Obtain the latest version of the 293T PPI network bp.293t <- getBioPlex(cell.line = "293T", version = "3.0") # (2) Obtain the latest version of the HCT116 PPI network bp.hct116 <- getBioPlex(cell.line = "HCT116", version = "1.0")
Functionality for retrieving the BioPlex3 protein expression data comparing expression in the HCT116 and the 293T cell lines.
getBioplexProteome(cache = TRUE)getBioplexProteome(cache = TRUE)
cache |
logical. Should a locally cached version used if available?
Defaults to |
A SummarizedExperiment storing protein expression
data for the both cell line(s) with 5 replicates each.
BioPlex: https://bioplex.hms.harvard.edu
se <- getBioplexProteome()se <- getBioplexProteome()
Functionality for retrieving the CORUM protein complex data. Available complex collections include:
complete set of complexes,
core set of complexes,
complexes with splice variants.
See references.
getCorum( set = c("all", "core", "splice"), organism = "Human", remap.uniprot.ids = FALSE, cache = TRUE, mode = c("ehub", "web") )getCorum( set = c("all", "core", "splice"), organism = "Human", remap.uniprot.ids = FALSE, cache = TRUE, mode = c("ehub", "web") )
set |
character. Valid options include:
Defaults to |
organism |
character. Use |
remap.uniprot.ids |
logical. Should the protein-to-gene mappings from CORUM
(i.e. UNIPROT-to-SYMBOL and UNIPROT-to-ENTREZID) be updated using Bioc annotation
functionality? Currently only supported in combination with |
cache |
logical. Should a locally cached version used if available?
Defaults to |
mode |
character. Should CORUM complexes be obtained from ExperimentHub
or via a web download from the CORUM homepage? Defaults to |
A data.frame.
CORUM: http://mips.helmholtz-muenchen.de/corum/#download
# Obtain the core set of CORUM complexes core <- getCorum(set = "core")# Obtain the core set of CORUM complexes core <- getCorum(set = "core")
Functionality for storing the 293T RNA-seq data
from GSE122425 in a
SummarizedExperiment. The dataset includes three
wild type samples and three NSUN2 knockout samples.
getGSE122425(cache = TRUE)getGSE122425(cache = TRUE)
cache |
logical. Should a locally cached version used if available?
Defaults to |
A SummarizedExperiment storing RNA-seq
data for the 293T cell line.
GSE122425: https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE122425
# Obtain the data as a SummarizedExperiment se <- getGSE122425()# Obtain the data as a SummarizedExperiment se <- getGSE122425()
Screens a list of graph instances storing
CORUM protein complex data for a subunit of choice.
hasSubunit(glist, subunit, id.type = "SYMBOL")hasSubunit(glist, subunit, id.type = "SYMBOL")
glist |
A |
subunit |
character. A gene ID corresponding to the subunit of interest. |
id.type |
character. Gene ID type of the given subunit. Defaults to |
A logical vector indicating which graphs have a node with the given subunit.
# (1) Obtain the core set of CORUM complexes ... core <- getCorum(set = "core") # (2) ... turn into a list of graphs ... core.glist <- corum2graphlist(core) # (3) .. check for a particular subunit of interest has.cdk2 <- hasSubunit(core.glist, subunit = "CDK2")# (1) Obtain the core set of CORUM complexes ... core <- getCorum(set = "core") # (2) ... turn into a list of graphs ... core.glist <- corum2graphlist(core) # (3) .. check for a particular subunit of interest has.cdk2 <- hasSubunit(core.glist, subunit = "CDK2")
Functionality for mapping experimental data stored in
a SummarizedExperiment onto a
graph object.
mapSummarizedExperimentOntoGraph( gr, se, col.names = NULL, rowdata.cols = NULL, prefix = "" )mapSummarizedExperimentOntoGraph( gr, se, col.names = NULL, rowdata.cols = NULL, prefix = "" )
gr |
an object of class |
se |
an object of class |
col.names |
character. Column names of |
rowdata.cols |
character. Column names of |
prefix |
character. Informative prefix that should be pasted together
with the selected |
An object of class graph.
# (1) Obtain the latest version of the 293T PPI network ... bp.293t <- getBioPlex(cell.line = "293T", version = "3.0") # (2) ... and turn into a graph bp.gr <- bioplex2graph(bp.293t) # (3) Obtain the BioPlex3 proteome data ... se <- getBioplexProteome() # (4) ... and map onto the graph bp.gr <- mapSummarizedExperimentOntoGraph(bp.gr, se)# (1) Obtain the latest version of the 293T PPI network ... bp.293t <- getBioPlex(cell.line = "293T", version = "3.0") # (2) ... and turn into a graph bp.gr <- bioplex2graph(bp.293t) # (3) Obtain the BioPlex3 proteome data ... se <- getBioplexProteome() # (4) ... and map onto the graph bp.gr <- mapSummarizedExperimentOntoGraph(bp.gr, se)