For Friday:

• read the introduction and the caveat section of the PSI Blast tutorial
• For those interested to learn about a popular alternative to PSI Blast:
  Read through http://en.wikipedia.org/wiki/HMMER and http://en.wikipedia.org/wiki/Hidden_Markov_model
  In case the principle of hidden profiles remains too cryptic, you can try the section in the HMMER manual on "What profile HMMs are".
• Questions due on Friday 5pm

For Monday:

• Think about project for which one could use PSI-blast or Hmmer searches.

For Wednesday:

• Take Exam #9 due
• Submit questions for study guide till 5pm .

For Today

• Bring one question to class regarding selection and the modern sysnthesis, or selection and evolution of modern humans.
• Ponder the possibility of complexity arising as a consequences of neutral processes (constructive neutral evolution). Do you know an example?

Examples for jumps in evolution (as opposed to gradual changes).

Constructive neutral evolution: Are mutations directed (tangle of cables as example)? Spliceeosomal introns, split inteins (discussion section 5 this article), mRNA editing, rearrangement of gene fragments in the transistion from micro to macronucleus in some ciliates (image here), emerging interactions between proteins (here, here, and here).

Slides for today on selection and PSI blast

PSI BLAST example PSI-blast provides an enormous advantage over normal blast in the detection of distantly related sequences.  It only works, if some closely related sequences are already available, but if this is the case it finds a lot of other distantly related sequences.  The NCBI page describes PSI blast as follows: Position-Specific Iterated BLAST (PSI-BLAST) provides an automated, easy-to-use version of a "profile" search, which is a sensitive way to look for sequence homologues. The program first performs a gapped BLAST database search. The PSI-BLAST program uses the information from any significant alignments returned to construct a position-specific score matrix, which replaces the query sequence for the next round of database searching. PSI-BLAST may be iterated until no new significant alignments are found. At this time PSI-BLAST may be used only for comparing protein queries with protein databases.  A diagram giving an overview on the PSI-blast procedure is here. The results of a normal blast search are aligned and a pattern of conserved residues is extracted from the alignment.  This pattern (the Position Specific Scoring Matrix) is used as query for the next iteration.  An important parameter to adjust is the E-value threshold up to which matches are included in the alignment and pattern extraction.  At higher iterations a PSI blast profile can be corrupted and false positives are identified with significant E-values.   I.e., in a traditional blast search one can be quite certain that a match with an E-value of 10^-13 represents a homologue; this is not clear with PSI blast.  Test studies indicate that profile corruptions are likely after more than 5 iterations. On the positive side: there are many fewer false negatives with PSI blast than with normal blast. False negatives: Homologous sequences that are not detected -- less in PSI blast than in normal blast False positives: Non-homologous sequences that are listed as matches -- very few in normal blast - possible in PSI blast after profile corruption The "problem" is that the E-value reported in a PSI-blast search represents the match with the profile, not with the original sequence!! PSI BLAST Example Query sequence: >gi|18892111|gb|AAL80306.1|[18892111] Pfu VMA intein (including -1 and +1 extein residues) KCVDGDTLILTKEFGLIKIKDLYEKLDGKGRKTVEGNEEWTELEEPITVYGYKNGKIVEIKATHVYKG ASSGMIEIKTRTGRKIKVTPIHKLFTGRVTKDGLVLEEVMAMHIKPGDRIAVVKKIDGGEYVKLDTSS VTKIKVPEVLNEELAEFLGYVIGDGTLKPRTVAIYNNDESLLKRANFlAMKLFGVSGKIVQERTVKAL LIHSKYLVDFLKKLGIPGNKKARTWKVPKElLLSPPSVVKAFINAYIACDGYYNKEKGEIEIVTASEE GAYGLTYLLAKLGIYATIRRKTINGREYYRVVISGKANLEKLGVKREARGYTSIDVVPVDVESIYEAL GRPYSELKKEGIEIHNYLSGENMSYETFRKFAKVVGLEEIAENHLQHILFDEVVEVNYISEPQEVYDI TTETHNFVGGNMPTLLHNT Blast page Use SWISSPROT Check format for PSI-blast Check inclusion limit Turn on Filter for low complexity Increase number of reported matches Point out manual selection of sequences to be included in profile Use PSSM to search something else (e.g., individual genomes (here))

Goals Class 24

• Understand the different approaches to detect genes under positive selection.
• Know that HGT helped create biochemical pathways that changed the face of this planet (acetoclastic methanogenesis, oxygen producing photosynthesis).
• Understand that neutral processes lead to more complex intermolecular relations that at least initially do not provide a selective advantage to the organism.
• Understand the purifying selection in small populations is less effective, and that even slightly deleterious changes in gene and genome architecture (gene and mRNA editing, self splicing and spliceosomal introns, split inteins) can be fixed by genetic drift.