[BioC] Pathway Information

John Zhang jzhang at jimmy.harvard.edu
Thu Jun 16 16:54:52 CEST 2005

Suppose you are looking for gene that are involved in human "Biosynthesis of 
steriods" pathway, the following will give you a vector of Entrez Gene ids:

>get(paste("hsa", get("Biosynthesis of steriods", KEGGPATHNAME2ID), sep = ""), 

>X-Original-To: jzhang at jimmy.harvard.edu
>Delivered-To: jzhang at jimmy.harvard.edu
>From: Ramsi Haddad <rhaddad at genetics.wayne.edu>
>To: michael.auer at meduniwien.ac.at
>Mime-Version: 1.0
>Date: 16 Jun 2005 09:58:13 -0400
>Received-SPF: pass (hypatia.math.ethz.ch: is authenticated by a 
trusted mechanism)
>Received-SPF: none (hypatia: domain of rhaddad at genetics.wayne.edu does not 
designate permitted sender hosts)
>X-Virus-Scanned: by amavisd-new at stat.math.ethz.ch
>Cc: bioconductor mailing list <bioconductor at stat.math.ethz.ch>
>Subject: Re: [BioC] Pathway Information
>X-BeenThere: bioconductor at stat.math.ethz.ch
>X-Mailman-Version: 2.1.6
>List-Id: The Bioconductor Project Mailing List <bioconductor.stat.math.ethz.ch>
>List-Unsubscribe: <https://stat.ethz.ch/mailman/listinfo/bioconductor>, 
<mailto:bioconductor-request at stat.math.ethz.ch?subject=unsubscribe>
>List-Archive: <https://stat.ethz.ch/pipermail/bioconductor>
>List-Post: <mailto:bioconductor at stat.math.ethz.ch>
>List-Help: <mailto:bioconductor-request at stat.math.ethz.ch?subject=help>
>List-Subscribe: <https://stat.ethz.ch/mailman/listinfo/bioconductor>, 
<mailto:bioconductor-request at stat.math.ethz.ch?subject=subscribe>
>Content-Transfer-Encoding: 7bit
>X-Spam-Checker-Version: SpamAssassin 3.0.1 (2004-10-22) on  
>X-Spam-Status: No, score=-2.2 required=3.0 tests=AWL,BAYES_00 autolearn=ham  
>Dear Michael,
>	There is a pretty good description in the globaltest vignette.  Most of
>the stuff shown below is based on it:
>### this next line generates a list of kegg pathways and all
>### the probe sets that belong in a given pathway.
>mouse.kegg.sets <- as.list(mgu74av2PATH2PROBE)
>### the next 3 lines give the names of the pathways from their
>### kegg ID numbers
>all.kegg.ids <- ls(KEGGPATHID2NAME)
>all.kegg.names <- mget(all.kegg.ids, KEGGPATHID2NAME)
>table.of.keggs <- as.matrix(unlist(all.kegg.names))
>write.table(table.of.keggs, "keggs.csv", sep=",",col.names=NA,
>### here's how to get the name of a given kegg pathway:
>what.is.00010 <- get("00010", KEGGPATHID2NAME) ## this is glycolysis.
>### here's how to get the probesets in the Glycolysis Path from the
>mouse array:
>glycolysis.affyIDs <- get("00010", mgu74av2PATH2PROBE) ## there are 68
>such probe sets.
>another.glycolysis.affyIDs <- mouse.kegg.sets[["00010"]] ## more than
>one way to skin a cat
>### here's how to get the kegg pathway numbers for anything with
>colera.kegg.paths.index <- grep("Cholera", table.of.keggs) ## this
>provides the index
>colera.keggs <- table.of.keggs[colera.kegg.paths]
>### here it is in one step:
>### table.of.keggs[grep("Cholera", table.of.keggs)]
>hope that gets you started.
>> Hi
>> Can anybody tell me the following
>> I want to identify genes which are located on a certain pathway NF Kappa.
>> How can I obtain such an information. Which package is appropriate and
>> which commands are needed. If any body ever encounterd the same problem
>> please let me know.
>> Michael
>Ramsi Haddad, Ph.D.
>Center for Molecular Medicine and Genetics.
>Functional Genomics Laboratory,
>Perinatology Research Branch, NICHD, NIH.
>259 Mack Avenue,
>Room 3146 Applebaum Bldg.
>Detroit, MI 48201, USA.
>phone:(313) 577-2569 / fax:(313) 577-7736
>Bioconductor mailing list
>Bioconductor at stat.math.ethz.ch

Jianhua Zhang
Department of Medical Oncology
Dana-Farber Cancer Institute
44 Binney Street
Boston, MA 02115-6084

More information about the Bioconductor mailing list