{r add_column, warning=FALSE, message=FALSE} # model[which(model$type=="Gene"),] #
{r add_column2, warning=FALSE, message=FALSE} # model[which(model$type=="Disease"),] #
InterMine is a powerful open source data warehouse system integrating diverse biological data sets (e.g. genomic, expression and protein data) for various organisms. Integrating data makes it possible to run sophisticated data mining queries that span domains of biological knowledge. A selected list of databases powered by InterMine is shown in Table 1:
Database | Organism | Data |
---|---|---|
FlyMine | Drosophila | Genes, homology, proteins, interactions, gene ontology, expression, regulation, phenotypes, pathways, diseases, resources, publications |
HumanMine | H. sapiens | Genomics, SNPs, GWAS, proteins, gene ontology, pathways, gene expression, interactions, publications, disease, orthologues, alleles |
MouseMine | M. musculus | Genomics, proteins, gene ontology, expression, interactions, pathways, phenotypes, diseases, homology, publications |
RatMine | R. norvegicus | Disease, gene ontology, genomics, interactions, phenotype, pathway, proteins, publication QTL, SNP |
WormMine | C. elegans | Genes, alleles, homology, go annotation, phenotypes, strains |
YeastMine | S. cerevisiae | Genomics, proteins, gene ontology, comparative genomics, phenotypes, interactions, literature, pathways, gene expression |
ZebrafishMine | D. rerio | Genes, constructs, disease, gene ontology, genotypes, homology, morpholinos, phenotypes |
TargetMine | H. sapiens, M. musculus | Genes, protein structures, chemical compounds, protein domains, gene function, pathways, interactions, disease, drug targets |
MitoMiner | H. sapiens, M. musculus, R. norvegicus, D. rerio, S. cerevisiae, S. pombe | Genes, homology, localisation evidence, Mitochondrial reference gene lists, phenotypes, diseases, expression, interactions, pathways, exome |
IndigoMine | Archae | Genomics |
ThaleMine | A. thaliana | Genomics, proteins, domains, homology, gene ontology, interactions, expression, publications, pathways, GeneRIF, stocks, phenotypes, alleles, insertions, TAIR |
MedicMine | Medicago truncatula | Genomics, pathways, gene ontology, publications, proteins, homology |
PhytoMine | over 50 plant genomes | Genes, proteins, expression, transcripts, homology |
Please see the InterMine home page for a full list of available InterMines.
InterMine includes an attractive, user-friendly web interface that works ‘out of the box’ and a powerful, scriptable web-service API to allow programmatic access to your data. This R package provides an interface with the InterMine-powered databases through Web services.
Let’s start with a simple task - find the pathways of the gene ABO.
First, we look at what databases are available.
## BMAP
## "https://bmap.jgi.doe.gov/bmapmine/"
## BeanMine
## "https://mines.legumeinfo.org/beanmine"
## BovineMine
## "http://genomes.missouri.edu/bovinemine"
## CHOmine
## "https://chomine.boku.ac.at/chomine"
## ChickpeaMine
## "https://mines.legumeinfo.org/chickpeamine"
## CovidMine
## "https://test.intermine.org/covidmine/"
## CowpeaMine
## "https://mines.legumeinfo.org/cowpeamine"
## FawMine
## "http://insectmine.org:8080/fawmine"
## FlyMine
## "https://www.flymine.org/flymine"
## GrapeMine
## "http://urgi.versailles.inra.fr/GrapeMine"
## HumanMine
## "https://www.humanmine.org/humanmine"
## HymenopteraMine
## "http://128.206.116.3:8080/hymenopteramine"
## IndigoMine
## "http://www.cbrc.kaust.edu.sa/indigo"
## JointvetchMine
## "https://mines.legumeinfo.org/jointvetchmine"
## LegumeMine
## "https://mines.legumeinfo.org/legumemine"
## LocustMine
## "http://locustmine.org:8080/locustmine"
## LupinMine
## "https://mines.legumeinfo.org/lupinmine"
## MaizeMine
## "http://maizemine.rnet.missouri.edu:8080/maizemine"
## MedicMine
## "https://mines.legumeinfo.org/medicmine"
## MitoMiner
## "http://mitominer.mrc-mbu.cam.ac.uk/release-4.0"
## MouseMine
## "http://www.mousemine.org/mousemine"
## OakMine
## "https://urgi.versailles.inra.fr/OakMine_PM1N"
## PeanutMine
## "https://mines.legumeinfo.org/peanutmine"
## PhytoMine
## "https://phytozome.jgi.doe.gov/phytomine/"
## PlanMine
## "http://planmine.mpi-cbg.de/planmine"
## RatMine
## "http://ratmine.mcw.edu/ratmine"
## RepetDB
## "http://urgi.versailles.inra.fr/repetdb"
## SoyMine
## "https://mines.legumeinfo.org/soymine"
## TargetMine
## "https://targetmine.mizuguchilab.org/targetmine"
## ThaleMine
## "https://bar.utoronto.ca/thalemine"
## WheatMine
## "https://urgi.versailles.inra.fr/WheatMine"
## WormMine
## "http://intermine.wormbase.org/tools/wormmine/"
## YeastMine
## "https://yeastmine.yeastgenome.org/yeastmine"
## ZebrafishMine
## "http://zebrafishmine.org"
Since we would like to query human genes, we select HumanMine.
## An object of class "Service"
## Slot "mine":
## HumanMine
## "https://www.humanmine.org/humanmine"
##
## Slot "token":
## [1] ""
Both in InterMine database website and in InterMineR, you are able to build custom queries. However, to facilitate the retrieval of information from InterMine databases, a variety of pre-built queries, called templates, have also been made available. Templates are queries that have already been created with a fixed set of output columns and one or more constraints.
## name
## 1 Tissue_Expression_illumina
## 2 humDisGeneOrthol2
## 3 PhenotypeGene
## 4 Protein_complex_details
## 5 disExprGene
## 6 protein_interactions
## title
## 1 Tissue --> Gene Expression (Illumina body map)
## 2 Human Disease --> Human Gene + Orthologue Gene(s)
## 3 Mouse Phenotype --> Mouse Genes + Orthologous genes
## 4 Protein --> Protein Complex
## 5 Disease Expression --> Genes
## 6 Protein --> Interactions
We would like to find templates involving genes.
## name
## 2 humDisGeneOrthol2
## 3 PhenotypeGene
## 5 disExprGene
## 7 Gene_Interactions2
## 8 Protein_Gene_Ortho
## 9 GOterm_Gene
## 12 Disease_gene_RNAseq
## 13 Gene_Alleles_Disease2
## 15 ChromRegion_GenesTransExon
## 18 GeneExpress
## 19 Disease_Genes2
## 20 Gene_Location
## 21 Protein_GeneChromosomeLength
## 22 Gene_Identifiers
## 24 Gene_Pathway
## 25 gene_complex_details
## 28 Gene_protein_sequence
## 29 PathwayGenes
## 31 Gene_Protein
## 32 Gene_OverlapppingGenes
## 35 Gene_To_Publications
## 36 Gene_Interactions_forReportPage
## 38 Gene_Disease_HPO
## 39 testPathwayGenes
## 40 Gene_GO
## 42 GeneInteractorsExpression
## 43 Gene_particularGoannotation
## 44 Gene_TissueExpressionIllumina
## 46 Gene_HPOphenotype_2
## 47 domain_protein_gene
## 48 Gene_Expression_GTex
## 49 Gene_ExpressionProteinAtlas
## 50 Pathway_ProteinGene
## 51 Gene_description
## 55 Gene_Interact_disease
## 56 GeneHPOparent_Genes_2
## 57 Gene_proteindomain
## 58 HPO_Gene
## 59 Gene_SigSNP
## 60 Gene_inGWAS
## 61 geneGWAS_reportPg
## 62 geneInteractiongene
## 63 Gene_Disease2
## 64 Term_inGWASoptionalGene
## 65 Gene_proteinAtlasExpression2
## 66 GeneOrthAllele
## 68 Gene_Orth2
## 71 Gene_Orth
## 72 ChromRegion_Genes
## 74 GenePathway_interactions2
## 75 Gene_AllelePhen
## title
## 2 Human Disease --> Human Gene + Orthologue Gene(s)
## 3 Mouse Phenotype --> Mouse Genes + Orthologous genes
## 5 Disease Expression --> Genes
## 7 Gene --> Interactions
## 8 Protein --> Gene and Orthologues
## 9 GO term --> Genes
## 12 Disease -> Genes + RNA-seq Expression
## 13 Gene --> Alleles and Disease (clinVar data)
## 15 Chromosomal Location --> All Genes + Transcripts + Exons
## 18 Gene --> Gene Expression (Tissue, Disease; Array Express, E-MTAB-62)
## 19 Disease --> Gene(s)
## 20 Gene --> Chromosomal location.
## 21 Protein --> Gene.
## 22 Gene --> All identifiers.
## 24 Gene --> Pathway
## 25 Gene --> protein complex
## 28 Gene -> Protein + protein sequence
## 29 Pathway --> Genes
## 31 Gene --> Proteins.
## 32 Gene --> Overlapping genes.
## 35 Gene --> Publications.
## 36 Gene --> Physical and Genetic Interactions
## 38 Gene --> Disease + HPO annotations (Human Phenotype Ontology
## 39 testPathway --> Genes
## 40 Gene --> GO terms.
## 42 Gene + Tissue Expression --> Interactors that are expressed in that tissue
## 43 Gene + GO term --> Genes by GO term
## 44 Gene --> Tissue Expression (Illumina body map)
## 46 Gene -> HPO annotation (Human Phenotype Ontology)
## 47 Protein Domain --> Protein and Genes
## 48 Gene --> Tissue Expression (GTex data)
## 49 Gene(s) --> Tissue Expression (Protein Atlas RNA-seq)
## 50 Pathway --> Protein and Gene
## 51 Gene -> Description
## 55 Gene -> Interactions + diseases
## 56 Gene + HPO Phenotype parent term -> Genes
## 57 Gene --> Protein + Domains
## 58 HPO term --> Genes
## 59 Gene(s) --> Significant SNPs (GTex data)
## 60 Gene --> GWAS hit
## 61 Gene Report --> GWAS hit
## 62 Gene A --> Interaction <-- Gene B
## 63 Gene --> Disease (OMIM)
## 64 GWAS term --> SNP + Associated gene if available
## 65 Gene --> Protein tissue Localisation
## 66 Gene (Hum OR Rat) --> Mouse Allele (Phenotype)
## 68 Gene --> Orthologues2
## 71 Gene --> Orthologues
## 72 Region --> Genes
## 74 Gene + Pathway --> Interactions
## 75 Mouse Gene --> Allele [Phenotype]
The template Gene_Pathway seems to be what we want. Let’s look at this template in more detail.
# Query for gene pathways
queryGenePath = getTemplateQuery(
im = im,
name = "Gene_Pathway"
)
queryGenePath
## $model
## name
## "genomic"
##
## $title
## [1] "Gene --> Pathway"
##
## $description
## [1] "For a given Gene (or List of Genes) show any associated Pathway(s) (Data Source: KEGG or REACTOME). Keywords: pathways, metabolism, cascade "
##
## $select
## [1] "Gene.primaryIdentifier" "Gene.symbol"
## [3] "Gene.name" "Gene.pathways.name"
## [5] "Gene.pathways.dataSets.name" "Gene.pathways.identifier"
## [7] "Gene.organism.shortName"
##
## $constraintLogic
## [1] "B and A"
##
## $name
## [1] "Gene_Pathway"
##
## $comment
## [1] "Added 26OCT2010: ML"
##
## $tags
## [1] "im:public" "im:report" "im:favourite"
##
## $rank
## [1] "1"
##
## $orderBy
## $orderBy[[1]]
## Gene.primaryIdentifier
## "ASC"
##
## $orderBy[[2]]
## Gene.primaryIdentifier
## "ASC"
##
## $orderBy[[3]]
## Gene.primaryIdentifier
## "ASC"
##
##
## $where
## $where[[1]]
## $where[[1]]$path
## [1] "Gene"
##
## $where[[1]]$op
## [1] "LOOKUP"
##
## $where[[1]]$code
## [1] "A"
##
## $where[[1]]$editable
## [1] TRUE
##
## $where[[1]]$switchable
## [1] FALSE
##
## $where[[1]]$switched
## [1] "LOCKED"
##
## $where[[1]]$value
## [1] "pparg"
##
##
## $where[[2]]
## $where[[2]]$path
## [1] "Gene.organism.name"
##
## $where[[2]]$op
## [1] "="
##
## $where[[2]]$code
## [1] "B"
##
## $where[[2]]$editable
## [1] TRUE
##
## $where[[2]]$switchable
## [1] FALSE
##
## $where[[2]]$switched
## [1] "LOCKED"
##
## $where[[2]]$value
## [1] "Homo sapiens"
There are three essential members in a query - SELECT, WHERE and constraintLogic.
What does ‘Gene.symbol’ mean? What is ‘Gene.pathway.identifier’?
Let’s take a look at the data model. NOTE: Section temporarily removed due to errors
Let's look at the children of the Gene data type.
# ```{r gene_type_data, warning=FALSE, message=FALSE}
# model[which(model$type=="Gene"),]
Gene has a field called ‘symbol’ (hence the column Gene.symbol). Gene also references the Pathways class, which is of the Pathway data type.
```
Let’s now run our template.
## Gene.primaryIdentifier Gene.symbol
## 1 5468 PPARG
## 2 5468 PPARG
## 3 5468 PPARG
## 4 5468 PPARG
## 5 5468 PPARG
## 6 5468 PPARG
## Gene.name
## 1 peroxisome proliferator activated receptor gamma
## 2 peroxisome proliferator activated receptor gamma
## 3 peroxisome proliferator activated receptor gamma
## 4 peroxisome proliferator activated receptor gamma
## 5 peroxisome proliferator activated receptor gamma
## 6 peroxisome proliferator activated receptor gamma
## Gene.pathways.name Gene.pathways.dataSets.name
## 1 Developmental Biology Reactome pathways data set
## 2 Gene expression (Transcription) Reactome pathways data set
## 3 Generic Transcription Pathway Reactome pathways data set
## 4 Intracellular signaling by second messengers Reactome pathways data set
## 5 MECP2 regulates transcription factors Reactome pathways data set
## 6 Metabolism Reactome pathways data set
## Gene.pathways.identifier Gene.organism.shortName
## 1 R-HSA-1266738 H. sapiens
## 2 R-HSA-74160 H. sapiens
## 3 R-HSA-212436 H. sapiens
## 4 R-HSA-9006925 H. sapiens
## 5 R-HSA-9022707 H. sapiens
## 6 R-HSA-1430728 H. sapiens
Let’s modify the query to find the pathways of the gene ABO. We want to change the ‘value’ attribute from PPARG to ABO.
There are two ways to build a query in InterMineR.
We can either build a query as a list object with newQuery
function, and assign all input values (selection of retrieved data type, constraints, etc.) as items of that list,
Or we can build the query as an InterMineR-class
object with the functions setConstraint
, which allows us to generate a new or modify an existing list of constraints, and setQuery
, which generates the query as a InterMineR-class
object.
setConstraints
and setQuery
functions are designed to facilitate the generation of queries for InterMine instances and avoid using multiple iterative loops, especially when it is required to include multiple constraints or constraint values (e.g. genes, organisms) in your query.
# modify directly the value of the first constraint from the list query
queryGenePath$where[[1]][["value"]] <- "ABO"
# or modify the value of the first constraint from the list query with setConstraints
queryGenePath$where = setConstraints(
modifyQueryConstraints = queryGenePath,
m.index = 1,
values = list("ABO")
)
queryGenePath$where
## [[1]]
## [[1]]$path
## [1] "Gene"
##
## [[1]]$op
## [1] "LOOKUP"
##
## [[1]]$code
## [1] "A"
##
## [[1]]$editable
## [1] TRUE
##
## [[1]]$switchable
## [1] FALSE
##
## [[1]]$switched
## [1] "LOCKED"
##
## [[1]]$value
## [1] "ABO"
##
##
## [[2]]
## [[2]]$path
## [1] "Gene.organism.name"
##
## [[2]]$op
## [1] "="
##
## [[2]]$code
## [1] "B"
##
## [[2]]$editable
## [1] TRUE
##
## [[2]]$switchable
## [1] FALSE
##
## [[2]]$switched
## [1] "LOCKED"
##
## [[2]]$value
## [1] "Homo sapiens"
Note the value is now equal to ‘ABO’. Let’s rerun our query with the new constraint.
## Gene.primaryIdentifier Gene.symbol
## 1 28 ABO
## 2 28 ABO
## 3 28 ABO
## 4 28 ABO
## Gene.name
## 1 ABO, alpha 1-3-N-acetylgalactosaminyltransferase and alpha 1-3-galactosyltransferase
## 2 ABO, alpha 1-3-N-acetylgalactosaminyltransferase and alpha 1-3-galactosyltransferase
## 3 ABO, alpha 1-3-N-acetylgalactosaminyltransferase and alpha 1-3-galactosyltransferase
## 4 ABO, alpha 1-3-N-acetylgalactosaminyltransferase and alpha 1-3-galactosyltransferase
## Gene.pathways.name Gene.pathways.dataSets.name
## 1 ABO blood group biosynthesis Reactome pathways data set
## 2 Blood group systems biosynthesis Reactome pathways data set
## 3 Metabolism Reactome pathways data set
## 4 Metabolism of carbohydrates Reactome pathways data set
## Gene.pathways.identifier Gene.organism.shortName
## 1 R-HSA-9033807 H. sapiens
## 2 R-HSA-9033658 H. sapiens
## 3 R-HSA-1430728 H. sapiens
## 4 R-HSA-71387 H. sapiens
Now we are seeing pathways for the ABO gene.
You can also add additional filters. Let’s look for a specifc pathway.
There are four parts of a constraint to add:
newConstraint <- list(
path=c("Gene.pathways.name"),
op=c("="),
value=c("ABO blood group biosynthesis"),
code=c("B")
)
queryGenePath$where[[2]] <- newConstraint
queryGenePath$where
## [[1]]
## [[1]]$path
## [1] "Gene"
##
## [[1]]$op
## [1] "LOOKUP"
##
## [[1]]$code
## [1] "A"
##
## [[1]]$editable
## [1] TRUE
##
## [[1]]$switchable
## [1] FALSE
##
## [[1]]$switched
## [1] "LOCKED"
##
## [[1]]$value
## [1] "ABO"
##
##
## [[2]]
## [[2]]$path
## [1] "Gene.pathways.name"
##
## [[2]]$op
## [1] "="
##
## [[2]]$value
## [1] "ABO blood group biosynthesis"
##
## [[2]]$code
## [1] "B"
Our new filter has been added successfully. Rerun the query and you see you only have one pathway,ABO blood group biosynthesis, returned.
## Gene.primaryIdentifier Gene.symbol
## 1 28 ABO
## 2 MGI:2135738 Abo
## Gene.name
## 1 ABO, alpha 1-3-N-acetylgalactosaminyltransferase and alpha 1-3-galactosyltransferase
## 2 ABO blood group (transferase A, alpha 1-3-N-acetylgalactosaminyltransferase, transferase B, alpha 1-3-galactosyltransferase)
## Gene.pathways.name Gene.pathways.dataSets.name
## 1 ABO blood group biosynthesis Reactome pathways data set
## 2 ABO blood group biosynthesis Reactome pathways data set
## Gene.pathways.identifier Gene.organism.shortName
## 1 R-HSA-9033807 H. sapiens
## 2 R-MMU-9033807 M. musculus
You can also add additional columns to the output. For instance, is the Gene also involved in any disease? Let’s add this information.
Let’s see what we know about diseases.
{r add_column, warning=FALSE, message=FALSE} # model[which(model$type=="Gene"),] #
The Gene data type has an ‘Diseases’ reference of type ‘Disease’.
{r add_column2, warning=FALSE, message=FALSE} # model[which(model$type=="Disease"),] #
Disease has an attribute called “name”. Add Gene.diseases.name to the view. We’ll add it as the last column, we can see from above there are 7 other columns already so we’ll put it as #8:
# use setQuery function which will create an InterMineR-class query
queryGenePath.InterMineR = setQuery(
inheritQuery = queryGenePath,
select = c(queryGenePath$select,
"Gene.diseases.name")
)
getSelect(queryGenePath.InterMineR)
## [1] "Gene.primaryIdentifier" "Gene.symbol"
## [3] "Gene.name" "Gene.pathways.name"
## [5] "Gene.pathways.dataSets.name" "Gene.pathways.identifier"
## [7] "Gene.organism.shortName" "Gene.diseases.name"
#queryGenePath.InterMineR@select
# or assign new column directly to the existing list query
queryGenePath$select[[8]] <- "Gene.diseases.name"
queryGenePath$select
## [1] "Gene.primaryIdentifier" "Gene.symbol"
## [3] "Gene.name" "Gene.pathways.name"
## [5] "Gene.pathways.dataSets.name" "Gene.pathways.identifier"
## [7] "Gene.organism.shortName" "Gene.diseases.name"
# run queries
resGenePath.InterMineR <- runQuery(im, queryGenePath.InterMineR)
resGenePath <- runQuery(im, queryGenePath)
all(resGenePath == resGenePath.InterMineR)
## [1] TRUE
## Gene.primaryIdentifier Gene.symbol
## 1 28 ABO
## Gene.name
## 1 ABO, alpha 1-3-N-acetylgalactosaminyltransferase and alpha 1-3-galactosyltransferase
## Gene.pathways.name Gene.pathways.dataSets.name
## 1 ABO blood group biosynthesis Reactome pathways data set
## Gene.pathways.identifier Gene.organism.shortName Gene.diseases.name
## 1 R-HSA-9033807 H. sapiens BLOOD GROUP, ABO SYSTEM
NB: adding columns can result in changing the row count.
The constraintLogic, if not given, is ‘A and B’. We would now try to explicitly specify the constraintLogic. A and B corresponds to the ‘code’ for each constraint.
## [1] "A and B"
Run the query again and see no change:
## Gene.primaryIdentifier Gene.symbol
## 1 28 ABO
## Gene.name
## 1 ABO, alpha 1-3-N-acetylgalactosaminyltransferase and alpha 1-3-galactosyltransferase
## Gene.pathways.name Gene.pathways.dataSets.name
## 1 ABO blood group biosynthesis Reactome pathways data set
## Gene.pathways.identifier Gene.organism.shortName Gene.diseases.name
## 1 R-HSA-9033807 H. sapiens BLOOD GROUP, ABO SYSTEM
Change to be ‘A or B’ and see how the results change.
- Start with the template Gene GO
## $model
## name
## "genomic"
##
## $title
## [1] "Gene --> GO terms."
##
## $description
## [1] "Search for GO annotations for a particular gene (or List of Genes)."
##
## $select
## [1] "Gene.primaryIdentifier"
## [2] "Gene.symbol"
## [3] "Gene.goAnnotation.ontologyTerm.identifier"
## [4] "Gene.goAnnotation.ontologyTerm.name"
## [5] "Gene.goAnnotation.ontologyTerm.namespace"
## [6] "Gene.goAnnotation.evidence.code.code"
## [7] "Gene.goAnnotation.ontologyTerm.parents.identifier"
## [8] "Gene.goAnnotation.ontologyTerm.parents.name"
## [9] "Gene.goAnnotation.qualifier"
##
## $name
## [1] "Gene_GO"
##
## $comment
## [1] "Added 15NOV2010: ML"
##
## $tags
## [1] "im:aspect:Function" "im:aspect:Gene Ontology"
## [3] "im:aspect:Genomics" "im:frontpage"
## [5] "im:public" "im:report"
##
## $rank
## [1] "4"
##
## $orderBy
## $orderBy[[1]]
## Gene.primaryIdentifier
## "ASC"
##
##
## $where
## $where[[1]]
## $where[[1]]$path
## [1] "Gene"
##
## $where[[1]]$op
## [1] "LOOKUP"
##
## $where[[1]]$code
## [1] "A"
##
## $where[[1]]$editable
## [1] TRUE
##
## $where[[1]]$switchable
## [1] FALSE
##
## $where[[1]]$switched
## [1] "LOCKED"
##
## $where[[1]]$value
## [1] "PPARG"
##
## $where[[1]]$extraValue
## [1] "H. sapiens"
- Modify the view to display a compact view
## [1] "Gene.symbol"
## [2] "Gene.goAnnotation.ontologyTerm.identifier"
## [3] "Gene.goAnnotation.ontologyTerm.name"
## [4] "Gene.goAnnotation.ontologyTerm.namespace"
- Modify the constraints to look for gene ABO.
## [[1]]
## [[1]]$path
## [1] "Gene"
##
## [[1]]$op
## [1] "LOOKUP"
##
## [[1]]$code
## [1] "A"
##
## [[1]]$editable
## [1] TRUE
##
## [[1]]$switchable
## [1] FALSE
##
## [[1]]$switched
## [1] "LOCKED"
##
## [[1]]$value
## [1] "ABO"
##
## [[1]]$extraValue
## [1] "H. sapiens"
- Run the query
## Gene.symbol Gene.goAnnotation.ontologyTerm.identifier
## 1 ABO GO:0000166
## 2 ABO GO:0003823
## 3 ABO GO:0004380
## 4 ABO GO:0004381
## 5 ABO GO:0005576
## 6 ABO GO:0005794
## Gene.goAnnotation.ontologyTerm.name
## 1 nucleotide binding
## 2 antigen binding
## 3 glycoprotein-fucosylgalactoside alpha-N-acetylgalactosaminyltransferase activity
## 4 fucosylgalactoside 3-alpha-galactosyltransferase activity
## 5 extracellular region
## 6 Golgi apparatus
## Gene.goAnnotation.ontologyTerm.namespace
## 1 molecular_function
## 2 molecular_function
## 3 molecular_function
## 4 molecular_function
## 5 cellular_component
## 6 cellular_component
- Start with the template Gene GO
## $model
## name
## "genomic"
##
## $title
## [1] "GO term --> Genes"
##
## $description
## [1] "Search for Genes in a specified organism that are associated with a particular Gene Ontology (GO) annotation."
##
## $select
## [1] "Gene.primaryIdentifier"
## [2] "Gene.symbol"
## [3] "Gene.name"
## [4] "Gene.goAnnotation.ontologyTerm.identifier"
## [5] "Gene.goAnnotation.ontologyTerm.name"
## [6] "Gene.organism.shortName"
##
## $constraintLogic
## [1] "A and B"
##
## $name
## [1] "GOterm_Gene"
##
## $comment
## [1] "Added 26OCT2010: ML"
##
## $tags
## [1] "im:aspect:Function" "im:aspect:Gene Ontology"
## [3] "im:aspect:Genomics" "im:public"
## [5] "im:report"
##
## $rank
## [1] "2"
##
## $orderBy
## $orderBy[[1]]
## Gene.symbol
## "ASC"
##
##
## $where
## $where[[1]]
## $where[[1]]$path
## [1] "Gene.goAnnotation.ontologyTerm.name"
##
## $where[[1]]$op
## [1] "LIKE"
##
## $where[[1]]$code
## [1] "A"
##
## $where[[1]]$editable
## [1] TRUE
##
## $where[[1]]$switchable
## [1] FALSE
##
## $where[[1]]$switched
## [1] "LOCKED"
##
## $where[[1]]$value
## [1] "DNA binding"
##
##
## $where[[2]]
## $where[[2]]$path
## [1] "Gene.organism.shortName"
##
## $where[[2]]$op
## [1] "="
##
## $where[[2]]$code
## [1] "B"
##
## $where[[2]]$editable
## [1] TRUE
##
## $where[[2]]$switchable
## [1] FALSE
##
## $where[[2]]$switched
## [1] "LOCKED"
##
## $where[[2]]$value
## [1] "H. sapiens"
- Modify the view to display a compact view
## [1] "Gene.symbol"
## [2] "Gene.name"
## [3] "Gene.goAnnotation.ontologyTerm.identifier"
## [4] "Gene.goAnnotation.ontologyTerm.name"
- Modify the constraints to look for GO term ‘metal ion binding’
## [[1]]
## [[1]]$path
## [1] "Gene.goAnnotation.ontologyTerm.name"
##
## [[1]]$op
## [1] "LIKE"
##
## [[1]]$code
## [1] "A"
##
## [[1]]$editable
## [1] TRUE
##
## [[1]]$switchable
## [1] FALSE
##
## [[1]]$switched
## [1] "LOCKED"
##
## [[1]]$value
## [1] "metal ion binding"
##
##
## [[2]]
## [[2]]$path
## [1] "Gene.organism.shortName"
##
## [[2]]$op
## [1] "="
##
## [[2]]$code
## [1] "B"
##
## [[2]]$editable
## [1] TRUE
##
## [[2]]$switchable
## [1] FALSE
##
## [[2]]$switched
## [1] "LOCKED"
##
## [[2]]$value
## [1] "H. sapiens"
- Run the query
## Gene.symbol Gene.name
## 1 A3GALT2 alpha 1,3-galactosyltransferase 2
## 2 AARSD1 alanyl-tRNA synthetase domain containing 1
## 3 ABAT 4-aminobutyrate aminotransferase
## 4 ABCG5 ATP binding cassette subfamily G member 5
## 5 ABCG8 ATP binding cassette subfamily G member 8
## 6 ABLIM1 actin binding LIM protein 1
## Gene.goAnnotation.ontologyTerm.identifier Gene.goAnnotation.ontologyTerm.name
## 1 GO:0046872 metal ion binding
## 2 GO:0046872 metal ion binding
## 3 GO:0046872 metal ion binding
## 4 GO:0046872 metal ion binding
## 5 GO:0046872 metal ion binding
## 6 GO:0046872 metal ion binding
- Start with the Gene_Location template, update to search for ABCA6 gene.
queryGeneLoc = getTemplateQuery(im, "Gene_Location")
queryGeneLoc$where[[2]][["value"]] = "ABCA6"
resGeneLoc= runQuery(im, queryGeneLoc)
resGeneLoc
## Gene.primaryIdentifier Gene.secondaryIdentifier Gene.symbol
## 1 23460 ENSG00000154262 ABCA6
## Gene.name Gene.chromosome.primaryIdentifier
## 1 ATP binding cassette subfamily A member 6 17
## Gene.locations.start Gene.locations.end Gene.locations.strand
## 1 69062044 69141927 -1
We’re going to use the output (gene location) as input for the next query.
- Define a new query
# set constraints
constraints = setConstraints(
paths = c(
"Gene.chromosome.primaryIdentifier",
"Gene.locations.start",
"Gene.locations.end",
"Gene.organism.name"
),
operators = c(
"=",
">=",
"<=",
"="
),
values = list(
resGeneLoc[1, "Gene.chromosome.primaryIdentifier"],
as.character(as.numeric(resGeneLoc[1, "Gene.locations.start"])-50000),
as.character(as.numeric(resGeneLoc[1, "Gene.locations.end"])+50000),
"Homo sapiens"
)
)
# set InterMineR-class query
queryNeighborGene = setQuery(
select = c("Gene.primaryIdentifier",
"Gene.symbol",
"Gene.chromosome.primaryIdentifier",
"Gene.locations.start",
"Gene.locations.end",
"Gene.locations.strand"),
where = constraints
)
summary(queryNeighborGene)
## path op value code
## 1 Gene.chromosome.primaryIdentifier = 17 A
## 2 Gene.locations.start >= 69012044 B
## 3 Gene.locations.end <= 69191927 C
## 4 Gene.organism.name = Homo sapiens D
- Run the query
## Gene.primaryIdentifier Gene.symbol Gene.chromosome.primaryIdentifier
## 1 100616316 MIR4524A 17
## 2 100847008 MIR4524B 17
## 3 23460 ABCA6 17
## Gene.locations.start Gene.locations.end Gene.locations.strand
## 1 69099564 69099632 -1
## 2 69099542 69099656 1
## 3 69062044 69141927 -1
- Plot the genes
resNeighborGene$Gene.locations.strand[which(resNeighborGene$Gene.locations.strand==1)]="+"
resNeighborGene$Gene.locations.strand[which(resNeighborGene$Gene.locations.strand==-1)]="-"
gene.idx = which(nchar(resNeighborGene$Gene.symbol)==0)
resNeighborGene$Gene.symbol[gene.idx]=resNeighborGene$Gene.primaryIdentifier[gene.idx]
annTrack = AnnotationTrack(
start=resNeighborGene$Gene.locations.start,
end=resNeighborGene$Gene.locations.end,
strand=resNeighborGene$Gene.locations.strand,
chromosome=resNeighborGene$Gene.chromosome.primaryIdentifier[1],
genome="GRCh38",
name="around ABCA6",
id=resNeighborGene$Gene.symbol)
gtr <- GenomeAxisTrack()
itr <- IdeogramTrack(genome="hg38", chromosome="chr17")
plotTracks(list(gtr, itr, annTrack), shape="box", showFeatureId=TRUE, fontcolor="black")
## R version 4.1.0 (2021-05-18)
## Platform: x86_64-pc-linux-gnu (64-bit)
## Running under: Ubuntu 20.04.2 LTS
##
## Matrix products: default
## BLAS: /home/biocbuild/bbs-3.13-bioc/R/lib/libRblas.so
## LAPACK: /home/biocbuild/bbs-3.13-bioc/R/lib/libRlapack.so
##
## locale:
## [1] LC_CTYPE=en_US.UTF-8 LC_NUMERIC=C
## [3] LC_TIME=C LC_COLLATE=C
## [5] LC_MONETARY=en_US.UTF-8 LC_MESSAGES=en_US.UTF-8
## [7] LC_PAPER=en_US.UTF-8 LC_NAME=C
## [9] LC_ADDRESS=C LC_TELEPHONE=C
## [11] LC_MEASUREMENT=en_US.UTF-8 LC_IDENTIFICATION=C
##
## attached base packages:
## [1] grid parallel stats4 stats graphics grDevices utils
## [8] datasets methods base
##
## other attached packages:
## [1] Gviz_1.36.0 GenomicRanges_1.44.0 GenomeInfoDb_1.28.0
## [4] IRanges_2.26.0 S4Vectors_0.30.0 BiocGenerics_0.38.0
## [7] InterMineR_1.14.1 BiocStyle_2.20.0
##
## loaded via a namespace (and not attached):
## [1] colorspace_2.0-1 rjson_0.2.20
## [3] ellipsis_0.3.2 biovizBase_1.40.0
## [5] htmlTable_2.2.1 XVector_0.32.0
## [7] base64enc_0.1-3 dichromat_2.0-0
## [9] rstudioapi_0.13 bit64_4.0.5
## [11] AnnotationDbi_1.54.0 fansi_0.5.0
## [13] sqldf_0.4-11 xml2_1.3.2
## [15] splines_4.1.0 cachem_1.0.5
## [17] knitr_1.33 Formula_1.2-4
## [19] jsonlite_1.7.2 Rsamtools_2.8.0
## [21] cluster_2.1.2 dbplyr_2.1.1
## [23] png_0.1-7 BiocManager_1.30.15
## [25] compiler_4.1.0 httr_1.4.2
## [27] backports_1.2.1 lazyeval_0.2.2
## [29] assertthat_0.2.1 Matrix_1.3-3
## [31] fastmap_1.1.0 htmltools_0.5.1.1
## [33] prettyunits_1.1.1 tools_4.1.0
## [35] igraph_1.2.6 gtable_0.3.0
## [37] glue_1.4.2 GenomeInfoDbData_1.2.6
## [39] dplyr_1.0.6 rappdirs_0.3.3
## [41] Rcpp_1.0.6 Biobase_2.52.0
## [43] jquerylib_0.1.4 vctrs_0.3.8
## [45] Biostrings_2.60.0 RJSONIO_1.3-1.4
## [47] rtracklayer_1.52.0 xfun_0.23
## [49] stringr_1.4.0 proto_1.0.0
## [51] lifecycle_1.0.0 ensembldb_2.16.0
## [53] restfulr_0.0.13 XML_3.99-0.6
## [55] zlibbioc_1.38.0 scales_1.1.1
## [57] BSgenome_1.60.0 VariantAnnotation_1.38.0
## [59] ProtGenerics_1.24.0 hms_1.1.0
## [61] MatrixGenerics_1.4.0 SummarizedExperiment_1.22.0
## [63] AnnotationFilter_1.16.0 RColorBrewer_1.1-2
## [65] yaml_2.2.1 curl_4.3.1
## [67] memoise_2.0.0 gridExtra_2.3
## [69] ggplot2_3.3.3 sass_0.4.0
## [71] biomaRt_2.48.0 rpart_4.1-15
## [73] latticeExtra_0.6-29 stringi_1.6.2
## [75] RSQLite_2.2.7 highr_0.9
## [77] BiocIO_1.2.0 checkmate_2.0.0
## [79] GenomicFeatures_1.44.0 filelock_1.0.2
## [81] BiocParallel_1.26.0 chron_2.3-56
## [83] rlang_0.4.11 pkgconfig_2.0.3
## [85] matrixStats_0.58.0 bitops_1.0-7
## [87] evaluate_0.14 lattice_0.20-44
## [89] purrr_0.3.4 htmlwidgets_1.5.3
## [91] GenomicAlignments_1.28.0 bit_4.0.4
## [93] tidyselect_1.1.1 magrittr_2.0.1
## [95] bookdown_0.22 R6_2.5.0
## [97] magick_2.7.2 generics_0.1.0
## [99] Hmisc_4.5-0 DelayedArray_0.18.0
## [101] DBI_1.1.1 gsubfn_0.7
## [103] pillar_1.6.1 foreign_0.8-81
## [105] survival_3.2-11 KEGGREST_1.32.0
## [107] RCurl_1.98-1.3 nnet_7.3-16
## [109] tibble_3.1.2 crayon_1.4.1
## [111] utf8_1.2.1 BiocFileCache_2.0.0
## [113] rmarkdown_2.8 jpeg_0.1-8.1
## [115] progress_1.2.2 data.table_1.14.0
## [117] blob_1.2.1 digest_0.6.27
## [119] munsell_0.5.0 bslib_0.2.5.1
The InterMine model could be accessed from the mine homepage by clicking the tab “QueryBuilder” and selecting the appropriate data type under “Select a Data Type to Begin a Query”:
Here we select Gene as the data type:
We could select Symbol and Chromosome->Primary Identifier by clicking Show on the right of them. Then click “Export XML” at the bottom right corner of the webpage: