1 / 54

Searching Molecular Databases with BLAST

Searching Molecular Databases with BLAST. Searching Molecular Databases with BLAST. B asic L ocal A lignment S earch T ool How BLAST works Interpreting search results The NCBI Web BLAST interface Demonstration and exercises. Why learn sequence database searching?. What have I cloned ?

MikeCarlo
Télécharger la présentation

Searching Molecular Databases with BLAST

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. Searching Molecular Databases with BLAST

  2. Searching Molecular Databases with BLAST • Basic Local Alignment Search Tool • How BLAST works • Interpreting search results • The NCBI Web BLAST interface • Demonstration and exercises

  3. Why learn sequence database searching? • What have I cloned ? • Is this really “my gene” ? • Has someone else already found it ? • What is this protein’s function ? • What is it related to ? • Can I get more sequence easily ?

  4. Search programs are sequence alignment programs • They try to find the best alignment between your probe sequence and every target sequence in the database • Finding optimal alignments is computationally a very resource intensive process • It is usually not necessary to find optimal alignments, particularly for large databases • Alignments are ranked and only top scores are reported

  5. Practical database search methods incorporate shortcuts • The fastest sequence database searching programs use heuristic algorithms • The basic concept is to break the search and alignment process down into several steps • At each step, only a best scoring subset is retained for further analysis

  6. What does ‘HEURISTIC’ mean? • “a commonsense rule (or set of rules) intended to increase the probability of solving some problem” • Why consider every possible alignment once a reasonably good alignment is found?

  7. Heuristic programs find approximate alignments • They are less sensitive than “dynamic programming” algorithms such as Smith-Waterman for detecting weak similarity • In practice, they run much faster and are usually adequate • The BLAST program developed by Stephen Altschul and coworkers at the NCBI is the most widely used heuristic program

  8. BLAST is a collection of five programs for different combinations of query and database sequences

  9. Why BLAST is great • Very fast and can be used to search extremely large databases • Sufficiently sensitive and selective for most purposes • Robust - the default parameters can usually be used

  10. BLAST scores are reported in two columns • Raw values based on the specific scoring matrix employed • As bits, which are matrix independent normalized values (bigger = better) • Significance is represented by E values (smaller = better)

  11. Sorted by E value Typical BLAST Output

  12. The EXPECT (E) threshold is used to control score reporting • A match will only be reported if its E value falls below the threshold set • The default value for E is 10, which means that 10 matches with scores this high are expected to be found by chance • Lower EXPECT thresholds are more stringent, and report fewer matches

  13. Interpreting BLAST scores • Score interpretation is based on context • What is the question? • What else do you know about the sequences? • Scoring is highly dependent on probe length • Exact matches will usually have the highest scores (and lowest E values) • Short exact matches may score lower than longer partial matches

  14. Interpreting BLAST scores • Short exact matches are expected to occur at random. • Partial matches over the entire length of a query are stronger evidence for homology than are short exact matches.

  15. Homology vs Identity • Homologous sequences are descended from a common ancestral sequence. • Homology is either true or false. It can never be partial! Saying two sequences are 45% homologous is a misuse of the term. • Sequence identity and similarity can be described as a percentage and are used as evidence of homology.

  16. BLAST Example Is this sequence known? What does it encode?

  17. Search Strategy • Choose the BLAST program: • nucleotide query vs. nucleotide db • megabalst: optimized to find identical sequences • blastn: will find identical and similar sequences • Choose the Database • nr (non-redundant) – everything • genome specific

  18. blastn Options Paste Query Sequence HERE Choose Database HERE Choose search program HERE

  19. Each line is a hit in the database sorted vertically by E value Colored rectangles along the X axis show where in the query sequence a similarity in the database has been found. Color indicates degree of similarity

  20. Output sorted by E value

  21. Link to GenBank file

  22. Link to alignment

  23. Link to Entrez Gene

  24. blastn Alignment

  25. BLASTP Example

  26. blastp input

  27. blastp Databases

  28. blastp Databases • nr - All non-redundant GenBank CDS translations + PDB + SwissProt+PIR • swissprot - the last major release of the SWISS-PROT protein sequence database • pat - patented sequences • pdb - Sequences derived from the 3-dimensional structure Protein Data Bank • env_nr - Non-redundant environmental samples

  29. BLASTP Output Conserved Domain Search Conserved domains are shown graphically. Link to explanation of the domain.

  30. blastp Output

  31. blastp Alignment

  32. Protein Scoring Matrices Blosom 62 is the default BLASTP scoring matrix

  33. Different Matrices Produce slightly different alignments

  34. Other BLAST Programs:Psi-BLAST 4.6 PSI-BLAST is designed for more sensitive protein-protein similarity searches. Position-Specific Iterated (PSI)-BLAST is the most sensitive BLAST program, making it useful for finding very distantly related proteins or new members of a protein family. Use PSI-BLAST when your standard protein-protein BLAST search either failed to find significant hits, or returned hits with descriptions such as "hypothetical protein" or "similar to...".

  35. Other BLAST Programs:Phi-BLAST 4.7 PHI-BLAST can do a restricted protein pattern search. Pattern-Hit Initiated (PHI)-BLAST is designed to search for proteins that contain a pattern specified by the user AND are similar to the query sequence in the vicinity of the pattern. This dual requirement is intended to reduce the number of database hits that contain the pattern, but are likely to have no true homology to the query.

  36. Sequence filters • Since only a limited number of matches are reported, hits to simple repeats and other low complexity sequences can obscure other more biologically meaningful similarities • Filters are used to remove low complexity sequences from the probe • Low Complexity, human repeats (blastn)

  37. Low Complexity Sequences are Filtered Out

  38. BLASTN vs BLASTP • Protein sequences have much higher information content than nucleotide sequence • To find evidence for sequence homology, use BLASTP and search protein sequences • Is my sequence already in the database? • To find identical sequences, search nucleotide databases

  39. Translated BLAST Searches • translations use all 6 frames • computationally intensive • tblastx searches can be very slow with some large databases • must specify genetic code

  40. Alternate Genetic Codes

  41. Translated BLAST Searches

  42. Taxonomy Reports

  43. Taxonomy Reports

  44. BLAST Genomes

More Related