1 / 31

Motif instance identification using comparative genomics

Motif instance identification using comparative genomics. Pouya Kheradpour Joint work with: Alexander Stark, Sushmita Roy and Manolis Kellis. Background and goal. Regulators bind to short (5 to 20bp) sequence specific patterns (motifs)

olisa
Download Presentation

Motif instance identification using comparative genomics

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript


  1. Motif instance identification using comparative genomics Pouya Kheradpour Joint work with: Alexander Stark, Sushmita Roy and Manolis Kellis

  2. Background and goal • Regulators bind to short (5 to 20bp) sequence specific patterns (motifs) • Genes are largely controlled through the binding of regulators • Transcription factors (TFs) are proteins that bind near the transcription start site (TSS) of genes and either activate or repress transcription • miRNAs bind to the 3’ un-translated region (UTR) of mRNAs to repress translation • The goal of our work is to identify these binding sites (motif instances) TF1 TF2 microRNA1

  3. Motivation Network: Davidson and Erwin, Science (2006) Mouse: Pennacchio, et al., Nature (2006) Fly: Tomancak, et al., Genome Biology (2002) • In all animals, genes are both temporally and spatially regulated to produce complex expression patterns • Identifying the targets of regulators is vital to understanding this expression • Conservation allows for identifying targets that are evolutionarily meaningful

  4. Previous work • Single genome approaches • Generally use positional clustering of motif matches to increase signal (e.g. Berman, et al. 2002; Schroeder, et al. 2004; Philippakis, et al. 2006) • A single 5mer match occurs on average 3 million times in mammalian genome • Requires set of specific factors that act together • Miss instances of motifs that may occur alone • Multi-genome approaches (phylogentic footprinting) • Blanchette and Tompa 2002 use an alignment free phylogenetic approach to find k-mers that are unusually well conserved • Moses, et al. 2004 use a strict phylogenetic model to find regions that evolve according to the motif and not the background • Etwiller, et al. 2005 use both nearby species and distant species (fish) to identify motif instances • Lewis, et al. 2005 finds putative microRNA binding sites requiring full conservation in five species

  5. Approach outline • Produce a raw conservation score for each motif match (branch length score or BLS) • For each motif and region, produce a mapping from BLS to confidence Advantages • Now we have many, complete, closely related genomes • Gives enough power to identify binding sites (Eddy, 2005) • Do not have to worry about dramatic divergence • Account for non-motif conservation using globally derived statistics • Robust against errors and evolutionary turnover • Computationally feasible to run genome wide for all available motifs

  6. Don’t require perfect conservation:  Branch length score Don’t require exact alignment:  Search within a window Large phylogeny challenges in instance identification Motif instance movement missing sequence • Sequencing / assembly / alignment artifacts • Low coverage sequencing, mis-alignments • Evolutionary variation • Individual binding sites can move / mutate • Some instances found only in subset of species

  7. Computing Branch Length Score (BLS) CTCF mutations movement missing short branches Does not over count redundant branch length Allows for: Mutations permitted by motif degeneracy Misalignment/movement of motifs within window (up to hundreds of nucleotides) Missing motif matches in dense species tree BLS = 2.23sps (78%)

  8. Branch Length Score  Confidence • Evaluate non-motif probability of a given score • Sequence could also be conserved due to overlap with un-annotated element (e.g. non-coding RNA) • Account for differences in motif composition and length • For example, short motif more likely to be conserved by chance

  9. Control motifs • Control motifs are the basis of our estimation of the background level of conservation and for evaluating enrichment • Each motif has its own set of controls • They are chosen to: • Have the same composition as the original motif • Match the target regions (e.g. promoters) with approximately the same frequency (+/- 20%) • Not too similar to each other (to preserve diversity) • Not be similar to known motifs (including the one being shuffled) • Background level is estimated separately in each region type (e.g. Promoters or 3’ UTRs)

  10. Branch Length Score  Confidence Use motif-specific shuffled control motifs determine the expected number of instances at each BLS by chance alone or due to non-motif conservation Compute Confidence Score as fraction of instances over noise at a given BLS(=1 – false discovery rate) Select movement window that leads to the most instances at each confidence

  11. Confidence selects for functional instances Transcription factor motifs MicroRNA motifs 3’UTR 3’UTR Intron Intron CDS CDS 5’UTR 5’UTR Promoter Promoter Confidence selects for transcription factor motif instances in promoters and miRNA motifs in 3’ UTRs

  12. Confidence selects for functional instances Strand Bias Confidence selects for transcription factor motif instances in promoters and miRNA motifs in 3’ UTRs miRNA motifs are found preferentially on the plus strand, whereas no such preference is found for TF motifs

  13. Experimental identification of binding sites • Chromatin immunoprecipitation (ChIP) combined with either sequencing (seq) or with microarrays (chip) are experimental procedures that are used to identify binding sites • Not all binding is functional, can have high false positive rate • Only binding that is active in the surveyed conditions is found ChIP-seq Maridis 2007

  14. Intersection with CTCF ChIP-Seq regions ChIP data from Barski, et al., Cell (2007) • Conserved CTCF motif instances highly enriched in ChIP-Seq sites • High enrichment does not require low sensitivity • Many motif instances are verified CTCF 50% motifs verified ≥ 50% of regions with a motif 50% confidence

  15. Enrichment found for other factors in mammals and flies Mammals Flies Zeitlinger, et al., Genes & Devel (2007) Sandmann, et al., Genes & Devel (2007) Robertson, et al., Nature Methods (2006) Odom, et al., Nature Genetics (2007) Barski, et al., Cell (2007) Abrams and Andrew, Devel (2005) (Not ChIP) Lin, et al., PLoS Genetics (2007) Lim, et al., Molecular Cell (2007) Sandmann, et al., Devel Cell (2006) Wei, et al., Cell (2006) Zeller, et al., PNAS (2006)

  16. Enrichment increases in conserved bound regions ChIP bound regions may not be conserved (Odom, et al. 2007) For CTCF we also have binding data in mouse Enrichment in intersection is dramatically higher Human: Barski, et al., Cell (2007) Mouse: Bernstein, unpublished

  17. Enrichment increases in conserved bound regions Odom, et al., Nature Genetics (2007) Human: Barski, et al., Cell (2007) Mouse: Bernstein, unpublished ChIP bound regions may not be conserved (Odom, et al. 2007) For CTCF we also have binding data in mouse Enrichment in intersection is dramatically higher Trend persists for other factors where we have multi-species ChIP data

  18. Enrichment of instancesin fly muscle genes • Motifs at 60% confidence and ChIP have similar enrichments (depletion for the repressor Snail) in the functional promoters • Enrichments persist even when you look at non-overlapping subsets • Intersection of two has strongest signal • Evolutionary and experimental evidence is complementary • ChIP includes species specific regions and differentiates tissues • Conserved instances include binding sites not seen in tissues surveyed ChIP data from: Zeitlinger, et al., G&D (2007); Sandmann, et al,. G&D (2007); Sandmann, et al., Dev Cell (2006)

  19. Fly regulatory network at 60% confidence TFs: 67 of 83 (81%) 46k instances miRNAs: 49 of 67 (86%) 4k instances • Several connections confirmed by literature (either directly or indirectly) • Global view of instances allows us to make network level observations: • TFs were more targeted by TFs (P < 10-20) and by miRNAs (P < 5 x 10-5) • TF in-degree associated with miRNA in-degree (high-high: P < 10-4; low-low P < 10-6)

  20. Contributions • A general methodology for regulatory motif instance identification using many, closely related genomes • Robust against errors from sequencing, assembly and alignment • Allows limited functional turnover and motif movement • Provides statistical measurement of confidence for each instance, correcting for length, composition and overlap with other functional elements • Validation and comparison to experimental data • High enrichment of binding sites in ChIP regions for a variety of factors • Functional enrichments suggest comparable ability to identify functional instances as ChIP

  21. Future directions • Our predicted network was static, but real regulatory networks are dynamic • They change throughout development and in different conditions • They can vary greatly in different species • We want to expand this work to learn about this network dynamics • ChIP data is becoming increasingly available in a variety of conditions – we can use this to learn what causes changes in binding • Multi-species data is also becoming more available • Can match motif binding to cross-species expression changes • We can train on this data to find motifs that act together or compensate for each other

  22. Acknowledgments • Alexander Stark • Sushmita Roy • Manolis Kellis Mouse CTCF ChIP-Seq • Tarjei Mikkelsen • Brad Bernstein Funding • William C.H. Chao Fellowship • NSF Graduate Research Fellowship • MIT CSAIL • Matt Rasmussen • Mike Lin • Issao Fujiwara • Rogerio Candeias • Broad Institute • Or Zuk • Michele Clamp • Manuel Garber • Mitch Guttman • Eric Lander

  23. The End

  24. Implementation details • Table lookup on the next 8 bases of the genome are used to find potential matches to the target genome • Results in an order-of-magnitude increase in speed over scanning through all motifs • In a first run, 100 shuffles of each motif are evaluated and up to 10 that fulfill the requirements are selected • All motifs and their selected shuffles are matched to the target genome and their BLS scores are computed • The matches are evaluated at each branch length cutoff and a mapping is produced for each motif from branch length score to confidence • All code is designed to run on BROAD cluster (often with parallelization) and is written in C

  25. Performance on mammalian TRANSFAC motifs 2.5x increase 3.5x 6.5x • Most motifs have confident instances into 90% confidence with 18 mammals • Substantial increase in the number of instances compared to only human, mouse rat and dog.

  26. The promise of many genomes • Eddy showed that with many genomes, resolving binding sites using conservation is possible • The goal of our work is to make this practical • Integrate evidence from multiple informant species • Determine which of the thousands of motif matches are functional using conservation

  27. Slides on motif discovery

  28. Related problem: computational motif discovery • Discovery of the regulatory motifs (as opposed to their binding sites) has also been an active area of research for several years • Single species work has generally required sequences thought to have similar regulation (for comparison, see Tompa, et al. 2005; Elemento, et al. 2007) • Looked for patterns that were enriched in target sequences • Use of conservation has been generally successful in re-identifying known binding affinities for TFs and miRNAs (e.g. Kellis, et al. 2003; Xie, et al. 2005; Etwiller, et al. 2005) • Requires fewer species (i.e. less branch length) than instance identification because signal can be integrated over thousands of instances found genome-wide

  29. Motif discovery pipeline • Enumerate motif seeds • Six non-degenerate characters with variable size gap in the middle • Score seed motifs • Use a conservation ratio corrected for composition and small counts to rank seed motifs • Expand seed motifs • Use expanded nucleotide IUPAC alphabet to fill unspecified bases around seed using hill climbing • Cluster to remove redundancy • Using sequence similarity gap gap S R T T G G C C Y W T T A A G G R

  30. Top 30 discovered fly motifs Many of the top discovered motifs match known motifs Motifs are associated with genes that are preferentially expressed in tissues

  31. Discovered motifs have functional enrichments Enrichment or depletion of a motif in the promoters of genes expressed in a tissue Tissues Motifs • Most motifs avoided in ubiquitously expressed genes 2. Functional clusters emerge

More Related