Comparative Genomics Comparative Genomics Common Themes Gene and - - PowerPoint PPT Presentation

comparative genomics comparative genomics
SMART_READER_LITE
LIVE PREVIEW

Comparative Genomics Comparative Genomics Common Themes Gene and - - PowerPoint PPT Presentation

Comparative Genomics Comparative Genomics Common Themes Gene and functional pathway presence/absence Identification, alignment of orthologs Structural analysis and comparison Comparative Genomics Basic Local Alignment Search Tool


slide-1
SLIDE 1

Comparative Genomics

slide-2
SLIDE 2

Comparative Genomics

Common Themes

  • Gene and functional pathway presence/absence
  • Identification, alignment of orthologs
  • Structural analysis and comparison
slide-3
SLIDE 3

Comparative Genomics

Basic Local Alignment Search Tool (BLAST)

slide-4
SLIDE 4

Comparative Genomics

The Core BLAST Programs

Program Query Sequences Database Sequences Alignment blastn nucleotide nucleotide nucleotide blastp amino acid amino acid amino acid blastx nucleotide amino acid amino acid tblastn amino acid nucleotide amino acid tblastx nucleotide nucleotide amino acid

slide-5
SLIDE 5

Comparative Genomics

Downstream Applications

  • Identify orthologs
  • Reciprocal BLAST
  • OrthoMCL
  • Annotation and Pathway Analysis
  • e.g., BLAST2GO
  • Many others
slide-6
SLIDE 6

Comparative Genomics

Advantages of Running BLAST Locally

  • Customize databases with unpublished sequences
  • Run thousands of simultaneous searches
  • Integrate into custom scripts and pipelines
slide-7
SLIDE 7

Comparative Genomics

Parsing BLAST Output

slide-8
SLIDE 8

Comparative Genomics

Organization of BLAST Output

  • RESULTS – One results block for each sequence in the
  • riginal query file
  • HITS – Within each results block, one hit block for

each sequence in the database that produced a significant alignment

  • HSPs – Within each hit block, one high-scoring-

pair block for each significant local alignment to that hit

slide-9
SLIDE 9

Comparative Genomics

Parsing BLAST Output

http://bioperl.org/howtos/SearchIO_HOWTO.html

slide-10
SLIDE 10

Comparative Genomics

Parsing BLAST Output – BioPerl Results Information

slide-11
SLIDE 11

Comparative Genomics

Parsing BLAST Output – BioPerl Hit Information

slide-12
SLIDE 12

Comparative Genomics

Parsing BLAST Output – BioPerl HSP Information

slide-13
SLIDE 13

Comparative Genomics

Parsing BLAST Output – BioPerl HSP Information (Cont.)

slide-14
SLIDE 14

Exercise

  • Make BLAST databases
  • Run local BLAST searches
  • Parse BLAST output with BioPerl
  • Make dot-plot comparing two bacterial genomes in R using parsed BLAST
  • utput

https://dbsloan.github.io/TS2019/exercises/local_blast.html