BGGN 213 Pathway Analysis and the Interpretation of Gene Lists Lecture 15 Barry Grant http://thegrantlab.org/bggn213 Inputs Steps Reads R1 Reads R2 [optional] Control FastQ FastQ 1. Quality Reads R1 Reads R2 Control [optional] FastQC FastQ FastQ Treatment 2. Reference Alignment Genome (Mapping) UCSC Fasta TopHat2 Count Last day’s step Table 4. 3. data.frame Annotation Read Differential Counting expression UCSC GTF Count analysis! CuffLinks Table DESeq2 data.frame Volcano Plot Fold change vs P-value Significant (P < 0.01 & log2 > 2) My high-throughput experiment generated a long list of genes/proteins… What do I do now? Pathway analysis! (a.k.a. geneset enrichment) Use bioinformatics methods to help extract biological meaning from such lists… Inputs Steps Reads R1 Reads R2 [optional] Control FastQ FastQ 1. Quality Reads R1 Reads R2 Control [optional] FastQC 5. Gene set FastQ FastQ enrichment Treatment analysis 2. nalysis! Reference KEGG, GO, … Alignment Genome (Mapping) UCSC Fasta TopHat2 Count Table 4. 3. data.frame Annotation Read Differential Counting expression UCSC GTF Count analysis! CuffLinks Table DESeq2 data.frame Basic idea Differentially Expressed Genes (DEGs) Gene-sets (Pathways, annotations, etc...) Annotate... Basic idea Differentially Expressed Genes (DEGs) Gene-sets (Pathways, annotations, etc...) Annotate... Differentially Pathways Expressed Overlap... Genes (DEGs) Pathway analysis (geneset enrichment) Pathway analysis (a.k.a. geneset enrichment) Principle DEGs Pathway Differentially Pathway Expressed (geneset) Genes (DEGs) Enriched Not enriched • DEGs come from your experiment ➢ Critical, needs to be as clean as possible • Pathway genes (“geneset”) come from annotations ➢ Important, but typically not a competitive advantage • Variations of the math: overlap, ranking, networks... ➢ Not critical, different algorithms show similar performances Side-note: Pathway analysis (a.k.a. geneset enrichment) Limitations • Geneset annotation bias: can only discover what is already known • Non-model organisms: no high-quality genesets available • Post-transcriptional regulation is neglected • Tissue-specific variations of pathways are not annotated • e.g. NF-κB regulates metabolism, not inflammation, in adipocytes • Size bias: stats are influenced by the size of the pathway • Many pathways/receptors converge to few regulators e.g. Tens of innate immune receptors activate four TFs: NF-kB, AP-1, IRF3/7, NFAT Starting point for pathway analysis: Your gene list • You have a list of genes/proteins of interest • You have quantitative data for each gene/protein 228018_at F26A1.8NP_000192 C20orf58 • Fold change 226925_atENSG00000090339 F26A1.8 NP_057219 ACPL2 207113_s_atENSG00000010030 3383F53F8.4 NP_055029 TNF • p-value 225078_atENSG00000141401 51513C09H5.2NP_000585 EMP2 221577_x_atENSG00000206439 3613C09H5.2NP_006125 GDF15 1553997_a_atENSG00000204490 7124F22E10.3NP_689495 ASPHD1 • Spectral counts 218433_atENSG00000206328 757F09F7.8NP_001032249 PANK3 204975_atENSG00000142188 92370C34F11.8NP_078870 EMP2 • Presence/absence 202638_s_atENSG00000155893 79646F07F6.5NP_064515 ICAM1 230954_atENSG00000120137 56892F35B3.4NP_733839 C20orf112 228017_s_atENSG00000176907 124540W08E12.3NP_620412 C20orf58 1554018_atENSG00000153944 253982B0222.10NP_859069 GPNMB 203126_atENSG00000174939 140688C03A7.14NP_542183 IMPA2 225182_atENSG00000197183 NP_00250110457Y40B10A.6 TMEM50B 225079_atENSG00000136235 NP_0010053409518M176.6 EMP2 NP_004855 243010_atENSG00000130513 2013H04M03.2 MSI2 NP_001415 230668_atENSG00000213853 4050C03A7.4 C20orf58 218541_s_atENSG00000204487 NP_033666C03A7.7 C8orf4 224225_s_atENSG00000206437 NP_002332F44F4.4 ETV7 207339_s_atENSG00000206327 F32A5.2 LTB 202637_s_at W03F8.6 ICAM1 Translating between identifiers • Many different identifiers exist for genes and proteins, e.g. UniProt, Entrez, etc. • Often you will have to translate one set of ids into another • A program might only accept certain types of ids • You might have a list of genes with one type of id and info for genes with another type of id Translating between identifiers • Many different identifiers exist for genes and proteins, e.g. UniProt, Entrez, etc. • Often you will have to translate one set of ids into another • A program might only accept certain types of ids • You might have a list of genes with one type of id and info for genes with another type of id • Various web sites translate ids -> best for small lists • UniProt < www.uniprot.org>; IDConverter < idconverter.bioinfo.cnio.es > Translating between identifiers: UniProt < www.uniprot.org > Translating between identifiers • Many different identifiers exist for genes and proteins, e.g. UniProt, Entrez, etc. • Often you will have to translate one set of ids into another • A program might only accept certain types of ids • You might have a list of genes with one type of id and info for genes with another type of id • Various web sites translate ids -> best for small lists • UniProt < www.uniprot.org>; IDConverter < idconverter.bioinfo.cnio.es > • VLOOKUP in Excel - good if you are an excel whizz - I am not! • Download flat file from Entrez, Uniprot, etc; Open in Excel; Find columns that correspond to the 2 IDs you want to convert between; Sort by ID; Use vlookup to translate your list Translating between identifiers: Excel VLOOKUP VLOOKUP(lookup_value, table_array, col_index_num) Translating between identifiers • Many different identifiers exist for genes and proteins, e.g. UniProt, Entrez, etc. • Often you will have to translate one set of ids into another • A program might only accept certain types of ids • You might have a list of genes with one type of id and info for genes with another type of id • Various web sites translate ids -> best for small lists • UniProt < www.uniprot.org >; IDConverter < idconverter.bioinfo.cnio.es > • VLOOKUP in Excel -> good if you are an excel whizz - I am not! • Download flat file from Entrez, Uniprot, etc; Open in Excel; Find columns that correspond to the two ids you want to convert between; Use vlookup to translate your list • Use the merge() or mapIDs() functions in R - fast, versatile & reproducible! • Also clusterProfiler::bitr() function and many others… [Link to clusterProfiler vignette] Reminder Using the merge() function anno <- read.csv("data/annotables_grch38.csv") This is an annotation file merge(mygenes, anno, by.x="row.names", by.y= "ensgene") Using mapIDs() Thisfunction is our from differential bioconductor expressed genes library("AnnotationDbi") library("org.Hs.eg.db") mygenes$symbol <- mapIds( org.Hs.eg.db, keys=row.names(mygenes), column="SYMBOL", keytype="ENSEMBL" ) Reminder Using the merge() function anno <- read.csv("data/annotables_grch38.csv") merge(mygenes, anno, by.x="row.names", by.y= "ensgene") Using mapIDs() function from bioconductor library("AnnotationDbi") library("org.Hs.eg.db") Load the required Bioconductor packages mygenes$symbol <- mapIds( org.Hs.eg.db, Annotation we want to add column="SYMBOL", keys=row.names(mygenes), keytype="ENSEMBL" ) Our vector of gene names & their format Alternative... See package vignette: https://bioconductor.org/packages/release/bioc/html/clusterProfiler.html What functional set databases do you want? Pathway • Most commonly used: DEGs GO KEGG • Gene Ontology (GO) IPA etc.... • KEGG Pathways (mostly metabolic) • GeneGO MetaBase • Ingenuity Pathway Analysis (IPA) • Many others... • Enzyme Classification, PFAM, Reactome, • Disease Ontology, MSigDB, Chemical Entities of Biological Interest, Network of Cancer Genes etc… • See: Open Biomedical Ontologies ( www.obofoundry.org ) GO < www.geneontology.org > • What function does HSF1 perform? • response to heat; sequence-specific DNA binding; transcription; etc • Ontology => a structured and controlled vocabulary that allows us to annotate gene products consistently, interpret the relationships among annotations, and can easily be handled by a computer • GO database consists of 3 ontologies that describe gene products in terms of their associated biological processes, cellular components and molecular functions GO Annotations • GO is not a stand-alone database of genes/proteins or sequences • Rather gene products get annotated with GO terms by UniProt and other organism specific databases, such as Flybase, Wormbase, MGI, ZFIN, etc. • Annotations are available through AmiGO < amigo.geneontology.org > GO is structured as a “directed graph” Relationships (edges) Parent terms are more general & child terms more specific GO Terms (nodes) GO evidence codes Use and misuse of the gene ontology annotations Seung Yon Rhee, Valerie Wood, Kara Dolinski & Sorin Draghici Nature Reviews Genetics 9, 509-515 (2008) • See AmiGO for details: http://amigo.geneontology.org/amigo/base_statistics Can now do gene list analysis with GeneGO online! Another popular online tool: DAVID at NIAID < david.abcc.ncifcrf.gov > DAVID • Functional Annotation Chart Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources Da Wei Huang, Brad T Sherman & Richard A Lempicki Nature Protocols 4, 44 - 57 (2009) Overlapping functional sets • Many functional sets overlap • In particular those from databases that are hierarchical in nature (e.g. GO) • Hierarchy enables: • Annotation flexibility (e.g. allow different degrees of annotation completeness based on what is known) • Computational methods to “understand” function relationships (e.g. ATPase function is a subset of enzyme function) • Unfortunately, this also makes functional profiling trickier • Clustering of functional sets can be helpful in these cases DAVID