Hi Fernando, In addition to Vince's suggestions, here is an approach where you actually compute the mapping yourself based on the gene locations. It uses the GenomicFeatures package to retrieve the transcript and exon locations for each gene, and then findOverlaps() to find the overlapping genes. All the information can be downloaded from the "RefSeq Genes" and "Ensembl Genes" UCSC tracks. The first track is associated with Entrez Gene IDs (NCBI) and the second track with Ensembl Gene IDs. library(GenomicFeatures) txdb1 <- makeTranscriptDbFromUCSC("bosTau4", tablename="refGene") txdb2 <- makeTranscriptDbFromUCSC("bosTau4", tablename="ensGene") genes1 <- range(transcriptsBy(txdb1, by="gene")) # 11668 Entrez genes genes2 <- range(transcriptsBy(txdb2, by="gene")) # 25669 Ensembl genes mm <- findOverlaps(genes1, genes2) ## Entrez genes that don't overlap: unmapped_ezids <- names(genes1)[-unique(queryHits(mm))] # 101 Entrez genes ## Ensembl genes that don't overlap: unmapped_ensids <- names(genes2)[-unique(subjectHits(mm))] # 12876 Ensembl genes ## Mapping between Entrez genes and Ensembl genes: map <- data.frame(ezid=names(genes1)[queryHits(mm)], ensid=names(genes2)[subjectHits(mm)]) As you can see, it's not a 1-to-1 mapping: > head(map) ezid ensid 1 100034674 ENSBTAG00000006026 2 100036590 ENSBTAG00000038843 3 100048947 ENSBTAG00000009091 4 100048949 ENSBTAG00000033312 5 100101492 ENSBTAG00000032159 6 100101492 ENSBTAG00000032234 Note that alternatively makeTranscriptDbFromBiomart() could be used to make 'txdb2' but that means you would first need to make sure that the current version of Ensembl (Ensembl 61) is also mapping the genes, transcripts and exons of Bos taurus to bosTau4. Hope this helps, H. On 03/08/2011 11:49 AM, Vincent Carey wrote: > There are various approaches using Bioconductor. The fundamental > resource is the package org.Bt.eg.db, which you can acquire using > biocLite(). > > You can find associations between ENSEMBL ids and Entrez ids using > mappings in that package. > > You may find GSEABase useful also. For example > >> dput(ensex) > c("ENSBTAG00000000005", "ENSBTAG00000000008", "ENSBTAG00000000009", > "ENSBTAG00000000010", "ENSBTAG00000000011", "ENSBTAG00000000012", > "ENSBTAG00000000013", "ENSBTAG00000000014", "ENSBTAG00000000015", > "ENSBTAG00000000016", "ENSBTAG00000000018", "ENSBTAG00000000019", > "ENSBTAG00000000020", "ENSBTAG00000000021", "ENSBTAG00000000022", > "ENSBTAG00000000023", "ENSBTAG00000000024", "ENSBTAG00000000025", > "ENSBTAG00000000026", "ENSBTAG00000000027") >> e1 = GeneSet(ensex, geneIdType=ENSEMBLIdentifier("org.Bt.eg.db")) >> e1 > setName: NA > geneIds: ENSBTAG00000000005, ENSBTAG00000000008, ..., > ENSBTAG00000000027 (total: 20) > geneIdType: ENSEMBL (org.Bt.eg.db) > collectionType: Null > details: use 'details(object)' >> g1 = e1 >> geneIdType(g1) = EntrezIdentifier("org.Bt.eg.db") >> g1 > setName: NA > geneIds: 282136, 539250, ..., 512788 (total: 21) > geneIdType: EntrezId (org.Bt.eg.db) > collectionType: Null > details: use 'details(object)' > > This shows that there are 21 Entrez IDs associated with the 20 ENSEMBL > ids given above. > After converting sets to common ID type, you can use intersect, > setdiff methods to answer some of the > questions you pose. > > >> sessionInfo() > R version 2.13.0 Under development (unstable) (2011-03-01 r54628) > Platform: x86_64-apple-darwin10.4.0/x86_64 (64-bit) > > locale: > [1] C > > attached base packages: > [1] stats graphics grDevices datasets tools utils methods > [8] base > > other attached packages: > [1] GSEABase_1.13.3 graph_1.29.3 annotate_1.29.3 > [4] org.Bt.eg.db_2.4.6 RSQLite_0.9-4 DBI_0.2-5 > [7] AnnotationDbi_1.13.15 Biobase_2.11.9 weaver_1.17.0 > [10] codetools_0.2-8 digest_0.4.2 > > loaded via a namespace (and not attached): > [1] Matrix_0.999375-47 XML_3.2-0 grid_2.13.0 lattice_0.19-17 > [5] xtable_1.5-6 > > > On Tue, Mar 8, 2011 at 12:13 PM, Biase, Fernando<biase at="" illinois.edu=""> wrote: >> Hi everyone, >> >> I have a list of ENSEMBL gene _IDS and a list with NCBI gene_IDs. I need to find which ids correspond to genes in both list (overlapping genes) and each genes are in each one of them but not present in the other list (non-overlapping genes). >> Can anyone give me some advice on this task? Or indicate a material do read? >> In case it is relevant, the organism is Bos taurus. >> >> Thanks in advance, >> Fernando >> >> _______________________________________________ >> Bioconductor mailing list >> Bioconductor at r-project.org >> https://stat.ethz.ch/mailman/listinfo/bioconductor >> Search the archives: http://news.gmane.org/gmane.science.biology.informatics.conductor >> > > _______________________________________________ > Bioconductor mailing list > Bioconductor at r-project.org > https://stat.ethz.ch/mailman/listinfo/bioconductor > Search the archives: http://news.gmane.org/gmane.science.biology.informatics.conductor -- Hervé Pagès Program in Computational Biology Division of Public Health Sciences Fred Hutchinson Cancer Research Center 1100 Fairview Ave. N, M2-B876 P.O. Box 19024 Seattle, WA 98109-1024 E-mail: hpages at fhcrc.org Phone: (206) 667-5791 Fax: (206) 667-1319