590 likes | 765 Views
Tools for understanding the sequence, evolution, and function of the human genome. Jim Kent and the Genome Bioinformatics Group University of California Santa Cruz. The Goal. Make the human genome understandable by humans. Step 1. Sequence the human genome.
E N D
Tools for understanding the sequence, evolution, and function of the human genome. Jim Kent and the Genome Bioinformatics Group University of California Santa Cruz
The Goal Make the human genome understandable by humans.
Step 1 Sequence the human genome
Mapping 300,000 BAC Clones Were Digested and Run on Agarose Gels Cari Soderlund’s FPC and Wash U Pathfinders Made Fingerprint Map Contigs Bob Waterston escaping management Genetic and radiation hybrid maps placed contigs on chromsomes
Sequence and Assembly • BAC Clones shotgun sequenced at high throughput to 4x ‘draft’. • Assembled with Phil Green’s Phrap
GigAssembler Jim Kent David Haussler (meanwhile Celera working on whole genome shotgun version)
The Truth + - ? + + ? + - ? - - ? + - ? + light - darkness Keeping strands straight is the hard part
“Finishing” Sequence • Using primers to end of contigs close gaps. • Checking automatic assembly especially near tandem repeats. • Checking in-silico restriction digest of BAC matches actual digest. • Time consuming - 1 year to ‘draft’ genome, 2 years to ‘finish’. • Human finished. Mouse will be finished (currently half finished). Other genomes may stay at draft stage, though draft stage can be very good these days.
Now What? TGGCTTTTGAAGGGAGTTCTGTTTATATATACGTCAACATCCAGTTGGAGGTGAAAAGGTTAGCACTTGACCCAGGAAGTATCCATGTTTGTTTCAAAAATAAATCTGCTTCATAAATTTCTTCATCAGTCTTTTTTTCCATTATGAGCTTTGATTATAATAAAGGAGCTGTTATTAACTTTTATTCAAGAAAAGGCCCATCTCTTTGAAAATATTTACCACCCTTCTCCCTTTCCCCTCATGAAATGTGCCAACTTCATAGGAATTAACAAATTGTAGCCCAGCCAAATACACGGATGCTTAAGCATACCTGAAACTTGAGTATATTTATTTATTACAGACATCCTAAGACCCGTAAACTCTGCTCTGGATCATATCACTCCAGGATCTCAGAGCTGTTCATGATTGTACAGGAAATGGGGAATATCATAGGCTCACAAAGGATAACTGATAGAACTCAGTGTGGTACTTTGGGGACATCAAACATTGTGCGACATGCAAAAGACTATTCACGAATAACACAAAATATACATTCATTGTGCCATCCATCACATTAACAATTGAGCTGAAAATACATTATATCCAGCTAAGATAACTGTGGAAGGAAGAAATTGGTTTGAATAATACTTTTAGGTTCTGAATAACCCAGCACAAATTTTAAACAGAGGGTGGCCCGAGAAGAAAGGGGTAGAGATTGGGAAAGACTTAGCACAGGAAGCCGGGTTTCTGAAGTTTGTGCTCTGCAGGGCTTCTTAACTGTAAGAACAAATCAAGGCTACCCTCTGAGGCATCTGATTGGGTTTAAATGAGGGAATTTTTTCTTTCACCTATAAAATTGTACCAGTTTAGAGAGTTTGCCCACCCTGTTTTAGTAACCTAAACATTTCTAGAAAATCTGTATAAAGATAAATCTCTTAGGACAAAGTATTTACAACCAGCAAACTCACACACATGAAAATGACTTAAATTAAGGGATGAATTAATTGTGTAAACATATAGTGCATCTCTTCTTCCTGAGCTCCTGGACTCGCCTTTCGCTATATCCTACTTTCAAGGACAAGGGAGGGGAGAGCTGTACATATAGTTAGATAAAAGATGAGAAGATTCCTTCTGGCATGTTTCTGTTGGCAAAGGGAACTATTTTCCAAAAGGTCATCTGAAAGGAACAGTAGGTTCTGTGAATTCTCCTAAAAGCAGGAGGGATGTTAAGGCCCACCAGAAAATGTATGCTGGCACCCAATCTGGATGAAGGTGTTAACCCCGCACCAAGTCTCTGGTCCAGAATTATCTGCAAATATATTATCCTGGCCAGGAGCTCCCCAGATAGGATTAGAAAGGAAGAAAGAGACTGTAAATGGAAAGAAAGATAAGCTAAGCATGTGCTTTGGGTAAGAAGTCCCAGCCCAAGGAGATGCCTGGGCTGTTGTCTGGGGCTGGAGCCGCCTCAGTGGGAGGTAGTCAGAGTGTCTGAGGTAGAAGACCCCGGGGAAGGAACGCAGGGCGAAGAGCTGGACTTCTCTGAGGATTCCTCGGCCTTCTCGTCGTTTCCTGGCGGGGTGGCCGGAGAGATGGGCAAGAGACCCTCCTTCTCACGTTTCTTTTGCTTCATTCGGCGGTTCTGGAACCAGATCTTCACTTGGGTCTCGTTGAGCTGCAGGGATGCAGCGATCTCCACCCTGCGGGCGCGCGTCAGGTACTTGTTGAAGTGGAACTCCTTCTCCAGTTCCGTGAGCTGCTTGGTAGTGAAGTTGGTGCGCACCGCGTTGGGTTGACCCAGGTAGCCGTACTCTCCAACTTTCCCTGGGGCAAAGTGGGAAGCCATGAGACGGAAATGTAAAAATTTTTAAATCGACTTGAGATTCCCCACACGCTTCATGGCAACACTCAGGTAAAGAAAAGATCAAGAACTCAGCACAAATCGGGCTGTGGAGGGTGAGTGATGAGGTGTAAAGTGTTAACCTGATGTAAACCATTAGCATGGTCAGACCGGTGATTAATGGAGCCTCAAGATATTAACAGAACACTACCGTCACAATAACCACCCCCACATACTTCCTATTTCCCAAATGTATAAAATCCTTGAAAACACACCAATCCCTGAGACTTCTTTGCCCCAACACCTCTGGGCACCCTCTCCATGCACTACAACACTAGTCTGATACAAAAGCCTTTTAAAAAAAAGATCATTATTAATTTCCTTGGAAATTAAGCATACCAGCTCCTTCCAGAATAATCAAGGAGCATCCACCAACCAGCAGGACTGACCTGTTTTGGGAGGGTTTCTTTTGACTTTCATCCAGTCAAAAGTCTGCGCTGGAGAAGATGTCTCCGATGCGGGGGAGCGACAGGCTTCTTGGTGGCTGGCGTGGAGAGGGGACAAGGAGTTATTATACGTAGCCAGGGCCAGGCTCTGGTGCTCCTGTCCATATGAGTGGTGAATGTATTGAGGCGAGCCCACCGCGCCCCCAGCATAACCCTGGTGGTGGTGGTGATGCTGGACCATGGGAGATGAGAGATTTCCAGAGTAAACAGCGGGAGCGCACTGGGGGTACCCACCACTTACGTCTGCTTCCTGATTTAACGCGTAGGGGCTGTAAGGCGCACTGAAGTTCTGTGAGCCATAGCTTGGACCACAACTTGAGTGGGAGTAGGACACCCCCAGGTTCCCGGAAGTCTGGTAGGTAGCCGGCTGGGGGTGGCGATGGTGGTGGTGGTGGTGGTGGTGGGGCGAACCGATCTGCACCCCCCTGCCCACTAGGAAGCGGTCGTCGCCGCCGCAACTGTTGGCGCTGACCGCGCACGACTGGAAAGTTGTAATCCTATGGTCCGAGGGGTAGGCTCGGGCTGAGCAGGTCCCCGAGTCGCCACTGCTAAGTATGGGGTATTCCAGGAAGGAGTTCATTCTTGCATTGTCCATCTGTCACTGAGTGACCTGGTCCTGCGAAGCCCGGCGTGACTGTGCCAACTTTCTCACTTCCTC
Finding the Genes Dr. Blat helping a gene find itself.
SIGLEC7 - a gene with some transcriptional complexity. Sialic Acid Binding/Ig-like Lectin 7 displayed in UCSC Genome Browser
Genes: Lines of Evidence • Full length human mRNA (the best!) • Protein homology with other species. • EST evidence - 1st step for much mRNA. • Evidence from genome/genome alignments • HMM based gene finders
Transferrin Clicking on a “known gene” brings up a large page of information on the gene.
Current state of human genome • ~99% of human genome sequenced. Last 1% will still be a challenge. • ~85% of human genes located. Substantial resources are being devoted to last 15%. • ~20% of human genes with any depth of functional annotation. Curation and integrated database are key to progress. • <1% of human regulatory regions located.
Transferrin Receptor Note peaks of conservation in 3’ UTR. These include iron response elements which regulate translation of this gene.
Comparative Genomics Webb Miller
Conservation of Gene Features Conservation pattern across 3165 mappings of human RefSeq mRNAs to the genome. A program sampled 200 evenly spaced bases across 500 bases upstream of transcription, the 5’ UTR, the first coding exon, introns, middle coding exons, introns, the 3’ UTR and 500 bases after polyadenylatoin. There are peaks of conservation at the transition from one region to another.
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.
Chains join together related local alignments Protease Regulatory Subunit 3
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.
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
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.
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.
Useful in finding pseudogenes Ensembl and Fgenesh++ automatic gene predictions confounded by numerous processed pseudogenes. Domain structure of resulting predicted protein must be interesting!
Mouse/HumanRearrangement Statistics Number of rearrangements of given type per megabase excluding known transposons.
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.
Reconstructed ancestral (boreutherian) genome for one chromosome
Finding Function • We’ve located 85% of the genes, on track for 95% in a year or two. • We have SOME idea of what 30% of the genes do. • We have virtually NO idea of what the rest do.
How to Find Function • Homology - guilt by association. Orthologs very valuable. • Genetics/knockouts - what happens when a gene gets broken? • RNAi is speeding this up amazingly in worms and other model organisms. • Expression - when and where is gene used? • Microarrays, in situs, GFP fusions. • Interactions - what molecules are touching? • Yeast 2 hybrid, Immunoprecipitations • Literature - finding out what we already know.
VisiGene • Image browser for in-situ and other gene- oriented pictures • Hopefully in the long run will have a million images covering almost all vertebrate genes. • Currently has 6000 images covering 1000 mouse transcription factors courtesy of Paul Gray et al.
Gene Browser Staff • Programming: Hiram Clawson, Mark Diekhans, Rachel Harte, Angie Hinrichs, Fan Hsu, Andy Pohl, Kate Rosenbloom, Chuck Sugnet, • Docs, quality, support: Gill Barber, Ron Chao, Jennifer Jackson, Donna Karolchik, Bob Kuhn, Crystal Lynch, Ali Sultan-Qurraie, Heather Trumbower • Computer systems: Jorge Garcia, Patrick Gavin, Paul Tatarsky