1 / 31

Introduction to Bioinformatics - Tutorial no. 2

Introduction to Bioinformatics - Tutorial no. 2. BLAST. BLAST. BLAST – Outline. Sequence Alignment Complexity and indexing BLASTN and BLASTP Basic parameters PAM and BLOSUM matrices Affine gap model E Values (once again). Advanced BLAST. Databases BLAST options BLAST output

goro
Download Presentation

Introduction to Bioinformatics - Tutorial no. 2

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. Introduction to Bioinformatics - Tutorial no. 2 BLAST

  2. BLAST

  3. BLAST – Outline • Sequence Alignment • Complexity and indexing • BLASTN and BLASTP • Basic parameters • PAM and BLOSUM matrices • Affine gap model • E Values (once again)

  4. Advanced BLAST • Databases • BLAST options • BLAST output • Taxonomic BLAST • Pairwise BLAST

  5. BLAST Variations • Genomic translations test all 6 possibilities: 3x for codon frames, 2x for reverse complement

  6. BLASTN Databases

  7. BLASTP Databases

  8. BLASTN/P Options (1) Only search part of database using NCBI Entrez query format Search specific organism Remove low information content, e.g. short repeats or rich in only 2 nucleotides Remove known human repeats (LINEs, SINEs)

  9. BLASTN/P Options (2) Threshold for results significance Allowed gap scores: 10/1, 10/2, 11/1, 8/2, 9/2 Use index based on words of 7, 11 or 15 nucleotides Costs to open and extend gap, score for nucleotide match or mismatch.

  10. BLASTP Options Scoring matrix: PAM, etc… Costs to open and extend gap Search for a motif (PSI-BLAST)

  11. BLASTN/P Formatting (1) Show colored bar chart Other (less important) options on what to show Number of sequences listed Number of alignments shown

  12. BLASTN/P Formatting (2) How to display alignments Only show results which match Entrez search or are from specific organism Only show results with E values in this range

  13. BLASTN Results Query sequence representation Matched areas of database sequences

  14. BLAST Output Header Request ID for later retrieval Query sequence details Database details Tax BLAST

  15. BLAST Alignments (1) Sequence Identifier Sequence description Score andE value

  16. BLAST Alignments (2) Normalized score of alignment Expected number of such hits (2e-11 = 2  10-11) Number of insertion / deletions Number of exact matches Number of matches with positive score

  17. BLAST Alignments (3) Insertion / deletion Exact match Query sequence Mismatch with positive score Matched sequence Position within sequence Masked low complexity region

  18. Expectation Values Increases linearly with length of query sequence Decreases exponentially with score of alignment Increases linearly with length of database

  19. Tax BLAST Lineage of organism with strongest hit Score of organism’s strongest hit Shared ancestry in taxonomic tree Number of organism hits

  20. Type of program BLAST2SEQ This tool produces the alignment of two given sequences using BLAST engine for local alignment. Scoring matrix Scoring scheme Gap model, Expect Value, Advanced options Sequences Sequences GO !

  21. Questions You have two query sequences: query1 and query2: >query1 CCGTCCGTCCGTCGTCCTCCTCGCTTGCGGGGCGCCGGGCCCGTCCTCGAGCCCCCNNNNNCCGTCCGGC CGCGTCGGGGCCTCGCCGCGCTCTACCTACCTACCTGGTTGATCCTGCCAGTAGCATATGCTTGTCTCAA AGATTAAGCCATGCATGTCTAAGTACGCACGGCCGGTACAGTGAAACTGCGAATGGCTCATTAAATCAGT TATGGTTCCTTTGGTCGCTCGCTCCTCTCCTACTTGGATAACTGTGGTAATTCTAGAGCTAATACATGCC GACGGGCGCTGACCCCCTTCGCGGGGGGGATGCGTGCATTTATCAGATCAAAACCAACCCGGTCAGCCCC TCTCCGGCCCCGGCCGGGGGGCGGGCCGCGGCGGCTTTGGTGACTCTAGATAACCTCGGGCCGATCGCAC GCCCCCCGTGGCGGCGACGACCCATTCGAACGTCTGCCCTATCAACTTTCGATGGTAGTCGCCGTGCCTA CCATGGTGACCACGGGTGACGGGGAATCAGGGTTCGATTCCGGAGAGGGAGCCTGAGAAACGGCTACCAC ATCCAAGGAAGGCAGCAGGCGCGCAAATTACCCACTCCCGACCCGGGGAGGTAGTGACGAAAAATAACAA TACAGGACTCTTTCGAGGCCCTGTAATTGGAATGAGTCCACTTTAAATCCTTTAACGAGGATCCATTGGA GGGCAAGTCTGGTGCCAGCAGCCGCGGTAATTCCAGCTCCAATAGCGTATATTAAAGTTGCTGCAGTTAA AAAGCTCGTAGTTGGATCTTGGGAGCGGGCGGGCGGTCCGCCGCGAGGCGAGCCACCGCCCGTCCCCGCC CCTTGCCTCTCGGCGCCCCCTCGATGCTCTTAGCTGAGTGTCCCGCGGGGCCCGAAGCGTTTACTTTGAA AAAATTAGAGTGTTCAAAGCAGGCCCGAGCCGCCTGGATACCGCAGCTAGGAATAATGGAATAGGACCGC GGTTCTATTTTGTTGGTTTTCGGAACTGAGGCCATGATTAAGAGGGACGGCCGGGGGCATTCGTATTGCG CCGCTAGAGGTGAAATTCTTGGACCGGCGCAAGACGGACCAGAGCGAAAGCATTTGCCAAGAATGTTTTC ATTAATCAAGAACGAAAGTCGGAGGTTCGAAGACGATCAGATACCGTCGTAGTTCCGACCATAAACGATG CCGACCGGCGATGCGGCGGCGTTATTCCCATGACCCGCCGGGCAGCTTCCGGGAAACCAAAGTCTTTGGG TTCCGGGGGGAGTATGGTTGCAAAGCTGAAACTTAAAGGAATTGACGGAAGGGCACCACCAGGAGTGGAG CCTGCGGCTTAATTTGACTCAACACGGGAAACCTCACCCGGCCCGGACACGGACAGGATTGACAGATTGA TAGCTCTTTCTCGATTCCGTGGGTGGTGGTGCATGGCCGTTCTTAGTTGGTGGAGCGATTTGTCTGGTTA ATTCCGATAACGAACGAGACTCTGGCATGCTAACTAGTTACGCGACCCCCGAGCGGTCGGCGTCCCCCAA CTTCTTAGAGGGACAAGTGGCGTTCAGCCACCCGAGATTGAGCAATAACAGGTCTGTGATGCCCTTAGAT GTCCGGGGCTGCACGCGCGCTACACTGACTGGCTCAGCGTGTGCCTACCCTACGCCGGCAGGCGCGGGTA ACCCGTTGAACCCCATTCGTGATGGGGATCGGGGATTGCAATTATTCCCCATGAACGAGGAATTCCCAGT AAGTGCGGGTCATAAGCTTGCGTTGATTAAGTCCCTGCCCTTTGTACACACCGCCCGTCGCTACTACCGA TTGGATGGTTTAGTGAGGCCCTCGGATCGGCCCCGCCGGGGTCGGCCCACGGCCTGGCGGAGCGCTGAGA AGACGGTCGAA

  22. Questions >query2 TACGAACGCTGGCGGCATGCTAATACATGCAAGTCGAACGAGACCTTCGGGTCTAGTGGCGCACGGGTGG CTAACGCGTGGGAATCTGCCCTTGGGTTCGGAATAACTTCGGGAAACTGAAGCTAATACCGGATGATGAC GAAAGTCCAAAGATTTATCGCCCAGGGATGAGCCCGCGTAGGATTAGCTAGTTGGTGGGGTAAAGGCTCA CCAAGGCAACGATCCTTAGCTGGTCTGAGAGGATGATCAGCCACACTGGGACTGAGACACGGCCCAGACT CCTACGGGAGGCAGCAGTAGGGAATATTGGACAATGGGCGAAAGCCTGATCCAGCAATGCCGCGTGAGTG ATGAAGGCCTTAGGGTTGTAAAGCTCTTTTACCCGAGATGATAATGACAGTATCGGGAGAATAAGCTCCG GCTAACTCCGTGCCAGCAGCCGCGGTAATACGGAGGGAGCTAGCGTTGTTCGGAATTACTGGGCGTAAAG CGCACGTAGGCGGCGATTTAAGTCAGAGGTGAAAGCCCGGGCTCAACCCCGAACTGCCTTTGAGACTGGA TTGCTAGAATCTTGGAGAGGCGAGTGGAATTCCGAGTGTAGAGGTGAAATTCGTAGATATTCGGAAGAAC ACCAGTGCGAAGGCGGCTCGCTGGACAAGTATTGACGCTGAGGTGCGAAAGCGTGGGGAGCAAACAGGAT TAGATACCCTGGTAGTCCACGCCGTAAACGATGATAACTAGCTGCCGGGGCACATGGTGTTTCGGTGGCG CACGTAACGCATTAAGTTATCCGCCTGGGGAGTACGGTCGCAAGATTAAAACTCAAAGGAATTGACGGGG GCCTGCACAAGCGGTGGAGCATGTGGTTTAATTCGAAGCAACGCGCAGAACCTTACCAGCGTTTGACATC CTCATCGCGGATTTCAGAGATGATTTCCTTCAGTTCGGCTGGATGAGTGACAGGTGCTGCATGGCTGTCG TCAGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTCGCCTTTAGTTGCCAGCAT TTAGTTGGGTACTCTAAAGGAACCGCCGGTGATAAGCCGGAGAAGGTGGGGATGACGTCAAGTCCTCATG GCCCTTACGCGCTGGGCTACACACGTGCTACAATGGCGACTACAGTGGGCTGCAACCGTGCGAGCGGTAG CTAATCTCCAAAAGTCGTCTCAGTTCGGATTGTTCTCTGCAACTCGAGAGCATGAAGGCGGAATCGCTAG TAATCGCGGATCAGCATGCCGCGGTGAATACGTTCCCAGGCCTTGTACACACCGCCCGTCACACCATGGG ATTTGGATTCACCCGAAGGCACTGCGTTAACCCGCAAGGGAGACAGGTGACCACGGTGGGTTTAGAGACT GGGGTGAA

  23. Questions • Using BLASTN • Find what do each one of these sequences code for.

  24. Questions

  25. Questions • To which organism each sequence is related? • Do these sequences code for proteins? Pretend the information for answering previous questions is not available to you could you suggest a way to answer these questions anyway? BLASTX

  26. Questions • Look carefully at the e-value column of the first 50 results of each query. What can you learn about these sequences? Are these sequences generally conserved between other organisms? 5 last answers

  27. Questions • Use bl2seq to align the two query sequences. What can you say about the relation between them? Based does this last result make sense?

  28. Questions You have two query sequences. >query3 ATGTCTGCTCCACAAGCCAAGATTTTGTCTCAAGCTCCAACTGAATTGGAATTACAAGTT GCTCAAGCTTTCGTTGAATTGGAAAATTCTTCTCCAGAATTGAAAGCTGAGTTGAGACCT TTGCAATTCAAGTCCATCAGAGAAGT >query4 GTATGTTATTAATTTGAATCTAAACTTAAGAATAATGGAGAGTAACAAAGGAAAAAAGTG TGAACGGGACGATACCAGAATGTTTCAATCTAGAAAAGTATAAAAGATAAGGACTAGGAC TCAAATGTATTTGGCTGACTATCGCCTGAACCTTGATGCTAAGCAAATACCATATCTTCA AGAAAAAGCCTACTCCAGTGTTTAAGAAGAAGGGAACGATTTACTAGATCATGCTATACG CAGTAAGGTTCTGATAGTTAATTACAATCGGTCCAAGTTCTAAGCGGTGTCGTCCATGCA TATATCATTTACAAGTTACTGGCGTCAACTCTTCAAATATTCAAAATATCACCTAATCAA ACTTACTAACATTTTCCTTTTTTGTTTTCCTTCTTTTATAG • Now use BlastX • To what protein does these sequences code for? • are these proteins conserved in other organisms?

  29. Questions • Now use BlastX • To what protein does these sequences code for? • are these proteins conserved in other organisms? Query 3 Query 4 No protein – e-value 3.2 A conserved protein component of the small (40S) subunit of S. cerevisiae.

  30. Questions • You are told that the sequences were extracted from the same gene. How could you explain the above results? • Answer: query4 is extracted from a non-coding region (intron) and thus doesn’t code for any protein.

More Related