1 / 39

TTh  11:00-12:15 in Clark S361 Profs: Serafim Batzoglou, Gill Bejerano

TTh  11:00-12:15 in Clark S361 Profs: Serafim Batzoglou, Gill Bejerano TAs: George Asimenos, Cory McLean. Lecture 18. Chains & Nets Non-coding Transcripts. Chaining Alignments. Chaining bridges the gulf between syntenic blocks and base-by-base alignments.

hchris
Download Presentation

TTh  11:00-12:15 in Clark S361 Profs: Serafim Batzoglou, Gill Bejerano

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. TTh  11:00-12:15 in Clark S361 Profs: Serafim Batzoglou, Gill Bejerano TAs: George Asimenos, Cory McLean http://cs273a.stanford.edu [Bejerano Spr06/07]

  2. Lecture 18 • Chains & Nets • Non-coding Transcripts http://cs273a.stanford.edu [Bejerano Spr06/07]

  3. Chaining Alignments • Chaining bridges the gulf between syntenic blocks and base-by-base alignments. • Local alignments tend to break at transposon insertions, inversions, duplications, etc. • Global alignments tend to force non-homologous bases to align. • Chaining is a rigorous way of joining together local alignments into larger structures. [Jim Kent’s slides]

  4. Chains join together related local alignments Protease Regulatory Subunit 3

  5. Chains • a chain is a sequence of gapless aligned blocks, where there must be no overlaps of blocks' target or query coords within the chain. • Within a chain, target and query coords are monotonically non-decreasing. (i.e. always increasing or flat) • double-sided gaps are a new capability (blastz can't do that) that allow extremely long chains to be constructed. • not just orthologs, but paralogs too, can result in good chains. but that's useful! • chains should be symmetrical -- e.g. swap human-mouse -> mouse-human chains, and you should get approx. the same chains as if you chain swapped mouse-human blastz alignments. • chained blastz alignments are not single-coverage in either target or query unless some subsequent filtering (like netting) is done. • chain tracks can contain massive pileups when a piece of the target aligns well to many places in the query. Common causes of this include insufficient masking of repeats and high-copy-number genes (or paralogs). [Angie Hinrichs, UCSC wiki] http://cs273a.stanford.edu [Bejerano Spr06/07]

  6. Affine penalties are too harsh for long gaps Log count of gaps vs. size of gaps in mouse/human alignment correlated with sizes of transposon relics. Affine gap scores model red/blue plots as straight lines.

  7. Before and After Chaining

  8. Chaining Algorithm • Input - blocks of gapless alignments from blastz • Dynamic program based on the recurrence relationship:score(Bi) = max(score(Bj) + match(Bi) - gap(Bi, Bj)) • Uses Miller’s KD-tree algorithm to minimize which parts of dynamic programming graph to traverse. Timing is O(N logN), where N is number of blocks (which is in hundreds of thousands) j<i

  9. Netting Alignments • Commonly multiple mouse alignments can be found for a particular human region, particularly for coding regions. • Net finds best match mouse match for each human region. • Highest scoring chains are used first. • Lower scoring chains fill in gaps within chains inducing a natural hierarchy.

  10. Net Focuses on Ortholog

  11. Nets • a net is a hierarchical collection of chains, with the highest-scoring non-overlapping chains on top, and their gaps filled in where possible by lower-scoring chains, for several levels. • a net is single-coverage for target but not for query. • because it's single-coverage in the target, it's no longer symmetrical. • the netter has two outputs, one of which we usually ignore: the target-centric net in query coordinates. The reciprocal best process uses that output: the query-referenced (but target-centric / target single-cov) net is turned back into component chains, and then those are netted to get single coverage in the query too; the two outputs of that netting are reciprocal-best in query and target coords. Reciprocal-best nets are symmetrical again. • nets do a good job of filtering out massive pileups by collapsing them down to (usually) a single level. [Angie Hinrichs, UCSC wiki] http://cs273a.stanford.edu [Bejerano Spr06/07]

  12. "LiftOver chains" are actually chains extracted from nets, or chains filtered by the netting process. Same-species liftOver chains are generated by a series of scripts that use blat -fastMap as the alignment method. [Angie Hinrichs, UCSC wiki] http://cs273a.stanford.edu [Bejerano Spr06/07]

  13. Before and After Chaining

  14. Net highlights rearrangements A large gap in the top level of the net is filled by an inversion containing two genes. Numerous smaller gaps are filled in by local duplications and processed pseudo-genes.

  15. Useful in finding pseudogenes Ensembl and Fgenesh++ automatic gene predictions confounded by numerous processed pseudogenes. Domain structure of resulting predicted protein must be interesting!

  16. Mouse/HumanRearrangement Statistics Number of rearrangements of given type per megabase excluding known transposons.

  17. A Rearrangement Hot Spot Rearrangements are not evenly distributed. Roughly 5% of the genome is in hot spots of rearrangements such as this one. This 350,000 base region is between two very long chains on chromosome 7.

  18. Cautionary Note 1 http://cs273a.stanford.edu [Bejerano Spr06/07]

  19. Cautionary Note 2 http://cs273a.stanford.edu [Bejerano Spr06/07]

  20. Same Region… same in all the other fish http://cs273a.stanford.edu [Bejerano Spr06/07]

  21. Orthology vs. Paralogy http://cs273a.stanford.edu [Bejerano Spr06/07]

  22. non coding transcripts http://cs273a.stanford.edu [Bejerano Spr06/07]

  23. http://cs273a.stanford.edu [Bejerano Spr06/07]

  24. http://cs273a.stanford.edu [Bejerano Spr06/07]

  25. http://cs273a.stanford.edu [Bejerano Spr06/07]

  26. http://cs273a.stanford.edu [Bejerano Spr06/07]

  27. http://cs273a.stanford.edu [Bejerano Spr06/07]

  28. http://cs273a.stanford.edu [Bejerano Spr06/07]

  29. http://cs273a.stanford.edu [Bejerano Spr06/07]

  30. Human Specific Rapid Evolution maximally changed m c h r m h r 100%id 100%id http://cs273a.stanford.edu [Bejerano Spr06/07]

  31. Nearest Neighbor Model for RNA Secondary Structure Free Energy at 37 OC: Mathews, Disney, Childs, Schroeder, Zuker, & Turner. 2004. PNAS 101: 7287.

  32. http://cs273a.stanford.edu [Bejerano Spr06/07]

  33. http://cs273a.stanford.edu [Bejerano Spr06/07]

  34. Transcripts, transcripts everywhere Human Genome Leaky tx? Functional? Transcribed (Tx) Tx from both strands http://cs273a.stanford.edu [Bejerano Spr06/07]

  35. http://cs273a.stanford.edu [Bejerano Spr06/07]

  36. http://cs273a.stanford.edu [Bejerano Spr06/07]

  37. http://cs273a.stanford.edu [Bejerano Spr06/07]

  38. http://cs273a.stanford.edu [Bejerano Spr06/07]

  39. http://cs273a.stanford.edu [Bejerano Spr06/07]

More Related