Homolog Database: Mouse homologs of human nicotinic acetylcholine receptors subunits

Libraries

library(hom.Hs.inp.db)

## Loading required package: AnnotationDbi

## Loading required package: stats4

## Loading required package: BiocGenerics

## Loading required package: parallel

## 
## Attaching package: 'BiocGenerics'

## The following objects are masked from 'package:parallel':
## 
##     clusterApply, clusterApplyLB, clusterCall, clusterEvalQ,
##     clusterExport, clusterMap, parApply, parCapply, parLapply,
##     parLapplyLB, parRapply, parSapply, parSapplyLB

## The following objects are masked from 'package:stats':
## 
##     IQR, mad, sd, var, xtabs

## The following objects are masked from 'package:base':
## 
##     anyDuplicated, append, as.data.frame, basename, cbind,
##     colMeans, colnames, colSums, dirname, do.call, duplicated,
##     eval, evalq, Filter, Find, get, grep, grepl, intersect,
##     is.unsorted, lapply, lengths, Map, mapply, match, mget, order,
##     paste, pmax, pmax.int, pmin, pmin.int, Position, rank, rbind,
##     Reduce, rowMeans, rownames, rowSums, sapply, setdiff, sort,
##     table, tapply, union, unique, unsplit, which, which.max,
##     which.min

## Loading required package: Biobase

## Welcome to Bioconductor
## 
##     Vignettes contain introductory material; view with
##     'browseVignettes()'. To cite Bioconductor, see
##     'citation("Biobase")', and for packages 'citation("pkgname")'.

## Loading required package: IRanges

## Loading required package: S4Vectors

## 
## Attaching package: 'S4Vectors'

## The following object is masked from 'package:base':
## 
##     expand.grid

## 
## Attaching package: 'IRanges'

## The following object is masked from 'package:grDevices':
## 
##     windows

## 

List of unique ENSMBL protein ids for nAChR subunits obtained earlier

[1] “ENSP00000276410” “ENSP00000433871” “ENSP00000434659” “ENSP00000405989” “ENSP00000488147” [6] “ENSP00000494974” “ENSP00000312663” “ENSP00000250699” “ENSP00000437107” “ENSP00000432757” [11] “ENSP00000492221” “ENSP00000491344” “ENSP00000491163” “ENSP00000491113” “ENSP00000491764” [16] “ENSP00000492078” “ENSP00000357461” “ENSP00000489703” “ENSP00000490474” “ENSP00000490203” [21] “ENSP00000433913” “ENSP00000289957” “ENSP00000408819” “ENSP00000404950” “ENSP00000258385” [26] “ENSP00000398143” “ENSP00000410801” “ENSP00000438380” “ENSP00000293780” “ENSP00000497829” [31] “ENSP00000496907”

Homologous proteins

head(keytypes(hom.Hs.inp.db))#lists species, there are 100 of them, 51=Hs, 60=Mm

## [1] "ACYRTHOSIPHON_PISUM"   "AEDES_AEGYPTI"         "ANOPHELES_GAMBIAE"    
## [4] "APIS_MELLIFERA"        "ARABIDOPSIS_THALIANA"  "ASPERGILLUS_FUMIGATUS"

head(columns(hom.Hs.inp.db))##same list of species 

## [1] "ACYRTHOSIPHON_PISUM"   "AEDES_AEGYPTI"         "ANOPHELES_GAMBIAE"    
## [4] "APIS_MELLIFERA"        "ARABIDOPSIS_THALIANA"  "ASPERGILLUS_FUMIGATUS"

Homolog of human CHRNA6 (Entrenz id =8973, ENSEMBLPROT=ENSP00000276410) in mouse

select(hom.Hs.inp.db, 
       keys="ENSP00000276410", 
       columns="MUS_MUSCULUS", 
       keytype="HOMO_SAPIENS")

## 'select()' returned 1:1 mapping between keys and columns

##      HOMO_SAPIENS       MUS_MUSCULUS
## 1 ENSP00000276410 ENSMUSP00000033882

Find homologs of human CHRN proteins in mouse

Make a vector of the ENSEMBLPROT ids

HsCHRNensemblID = c("ENSP00000276410", "ENSP00000433871", "ENSP00000434659", "ENSP00000405989", "ENSP00000488147", "ENSP00000494974", "ENSP00000312663", "ENSP00000250699", "ENSP00000437107", "ENSP00000432757", "ENSP00000492221", "ENSP00000491344", "ENSP00000491163", "ENSP00000491113", "ENSP00000491764", "ENSP00000492078", "ENSP00000357461", "ENSP00000489703", "ENSP00000490474", "ENSP00000490203", "ENSP00000433913", "ENSP00000289957", "ENSP00000408819", "ENSP00000404950", "ENSP00000258385", "ENSP00000398143", "ENSP00000410801", "ENSP00000438380", "ENSP00000293780",
"ENSP00000497829", "ENSP00000496907")

MmHomologsforHsCHRN <- select(hom.Hs.inp.db,
                              keys=HsCHRNensemblID, 
                              columns="MUS_MUSCULUS", 
                              keytype="HOMO_SAPIENS")

## Warning in `[<-.factor`(`*tmp*`, iseq, value = c("ENSP00000433871",
## "ENSP00000434659", : invalid factor level, NA generated

## 'select()' returned 1:many mapping between keys and columns

head(MmHomologsforHsCHRN)

##      HOMO_SAPIENS       MUS_MUSCULUS
## 1 ENSP00000276410 ENSMUSP00000033882
## 2 ENSP00000312663 ENSMUSP00000031108
## 3 ENSP00000250699 ENSMUSP00000081891
## 4 ENSP00000357461 ENSMUSP00000029562
## 5 ENSP00000289957 ENSMUSP00000052297
## 6 ENSP00000258385 ENSMUSP00000084162