* Your assessment is very important for improving the workof artificial intelligence, which forms the content of this project
Download Inferring Function From Known Genes
RNA interference wikipedia , lookup
Heritability of IQ wikipedia , lookup
Gene therapy wikipedia , lookup
Vectors in gene therapy wikipedia , lookup
Transposable element wikipedia , lookup
Pharmacogenomics wikipedia , lookup
X-inactivation wikipedia , lookup
Epigenetics in learning and memory wikipedia , lookup
Gene nomenclature wikipedia , lookup
Genetically modified crops wikipedia , lookup
Epigenetics of diabetes Type 2 wikipedia , lookup
Long non-coding RNA wikipedia , lookup
Oncogenomics wikipedia , lookup
Therapeutic gene modulation wikipedia , lookup
Genetic engineering wikipedia , lookup
Gene desert wikipedia , lookup
Epigenetics of neurodegenerative diseases wikipedia , lookup
Public health genomics wikipedia , lookup
Pathogenomics wikipedia , lookup
Quantitative trait locus wikipedia , lookup
Polycomb Group Proteins and Cancer wikipedia , lookup
Site-specific recombinase technology wikipedia , lookup
Essential gene wikipedia , lookup
Nutriepigenomics wikipedia , lookup
History of genetic engineering wikipedia , lookup
Gene expression programming wikipedia , lookup
Genome evolution wikipedia , lookup
Genomic imprinting wikipedia , lookup
Artificial gene synthesis wikipedia , lookup
Microevolution wikipedia , lookup
Designer baby wikipedia , lookup
Minimal genome wikipedia , lookup
Ridge (biology) wikipedia , lookup
Epigenetics of human development wikipedia , lookup
Genome (book) wikipedia , lookup
Inferring Function From Known Genes Naomi Altman Nov. 06 Objective There are 3 major objectives for microarray studies: 1) Understand the function of genes. 2) Understand a biological process. 3) Classify samples. 3 is "easy" since it is purely observational. 1&2 start from genetic methods (i.e. 1 gene at a time) Genetic Methods Tools: "knock-out" genotypes in which the gene is mutated to be non-functional "tagged" genotypes in which a fluorescent tag is activated when the gene is activated "in situ" hybridization in which the gene product is labeled in thin sections which can be viewed under a microscope "genetically modified" strains in which foreign genes have been added "genetically modified" strains in which selected genes are over- or under-expressed Using Known Genes There are several ways in which known genes can be used to infer the function of unknown genes in a microarray experiment. 1) Seeded clustering Assume that the nearest neighbors (according to some expression metric) of the known gene have similar function. 2) Unsupervised clustering majority rule - in a cluster consisting of both known and unknown genes, assume that the cluster has the majority function Using Known Genes There are several ways in which known genes can be used to infer the function of unknown genes in a microarray experiment. 3) Pathway analysis If the genes are sufficiently well understood, they may be assembled into networks showing which genes regulate other genes. Unknown genes that have expression patterns similar to those in the network can be placed in the network. BioPixie (for yeast) will be demonstrated by 2 project groups. PathAssist ? Understanding the Biological Process We also use the known genes to infer the biological processes underlying our experimental treatments. The primary tools are gene classification methods based on function and/or sequence. The most used tool is probably the Gene Ontology Project. Gene Ontology Project • http://en.wikipedia.org/wiki/Gene_Ontology • http://www.geneontology.org/ (try entering "DNA repair" "ribosome" "protease" or the term of your choice under gene or protein: try "ap2" How does GO "work" The terms are in a tree-like structure but some nodes have multiple "parents". The annotations are assigned by a team of collaborators and may also be submitted by other biologists (i.e. somewhat like Wikipedia, but with a more formalized central group) Each annotation also has an "evidence" annotation, which can be used to assess reliability. Using GO with Microarray Data 1) Compile a list of differentially expressed genes. 2) Obtain the GO annotations of all genes on the array. 3) Extract the GO annotations of the DE genes. Determine which GO annotations are over- or underrepresented among the DE genes. Descriptive Use of GO Testing For Enrichment DE not DE In N11 Out N12 N•1 N21 N1• N22 N2• N•2 N• • H0: The percentage of DE genes In the GO category is proportional to the number on the array in the category percentage on array: N1•/N•• expected: N•1N•1/N•• observed: N11 test (O-E)2/E (Chi-squared test) Testing For Enrichment DE not DE In N11 Out N12 N•1 N21 N1• N22 N2• N•2 N• • Problems: 1) Multiple testing: e.g. the Ontology for humans contains over 1000 terms 2) The ontology categories are nested. 3) The test statistic assumes that the genes are selected independently (but they are generally dependent). 4) Large N•• can lead to spuriously small p-values. GO in R The GOstats package in R will: take a genelist with Entrez IDs take a reference genelist with Entrez IDs take an annotation package (you can make your own) Find all the "significantly enriched" or "significantly depleted" GO categories. (The documentation was not very readable, but it was simple once I found an example.) Using GOstats with the Human vs Chimp Brain Data The limma output was saved in efit.contrast. I selected all genes with highly significant differential expression among the treatments. select=efit.contrast$F.p.value<.0000001 genes=efit.contrast$genes[select,] #affy probe id library(hgu95av2) w1<-as.list(hgu95av2ENTREZID) #Entrez id for every probeset on the array Using GOstats with the Human vs Chimp Brain Data entrez=w1[genes] entrezID=unlist(entrez) params <- new("GOHyperGParams", geneIds = entrezID, universeGeneIds = unlist(w1), annotation = "hgu95av2", ontology = "BP", pvalueCutoff = .001, conditional = FALSE, testDirection = "over") HGOver=hyperGTest(params) #tests all nodes htmlReport(HG0ver,"h.html")#report on sig nodes summary(HG0ver) #report in R Results Gene to GO BP test for over-representation GOBPID Pvalue OddsRatio GO:0008065 0.00 -29.22 GO:0007399 0.00 2.69 GO:0048731 0.00 2.65 GO:0019226 0.00 3.15 GO:0007417 0.00 4.36 GO:0007420 0.00 8.60 GO:0007268 0.00 2.97 GO:0050877 0.00 2.22 GO:0007154 0.00 1.49 GO:0030072 0.00 58.68 GO:0030073 0.00 58.68 GO:0007215 0.00 11.04 GO:0007267 0.00 1.80 GO:0006813 0.00 3.16 GO:0009187 0.00 5.42 ExpCount 0 21 21 12 6 2 11 25 133 0 0 1 27 5 2 Count 2 48 48 31 20 11 29 48 171 4 4 6 45 13 7 N11 Size 1.00 325.00 329.00 181.00 89.00 30.00 177.00 381.00 2066.00 5.00 5.00 14.00 425.00 74.00 26.00 N21 Term establishment of blood-nerve barri nervous system development system development transmission of nerve impulse central nervous system developme brain development synaptic transmission neurophysiological process cell communication peptide hormone secretion insulin secretion glutamate signaling pathway cell-cell signaling potassium ion transport cyclic nucleotide metabolism