1 / 44

Bioinformatics I

Bioinformatics I. Date: October 4, 2014. To have an understanding of: Sequence analysis Genome assembly and annotation Publicly available molecular biology and genetic databases Identification of sequence similarity Defining function of biological sequences

zarola
Download Presentation

Bioinformatics I

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. Bioinformatics I Date: October 4, 2014

  2. To have an understanding of: Sequence analysis Genome assembly and annotation Publicly available molecular biology and genetic databases Identification of sequence similarity Defining function of biological sequences Using sequence information to hypothesize function Learning Objectives

  3. Bioinformatics is the intersection of computer science, statistics and one of the following life science disciplines: biology, biochemistry or medicine Bioinformatics Introduction Bioinformatics Computer Science Statistics

  4. Clinical Informatics Systems used to deal with patient data Clinical trial management systems, electronic health records, etc. Laboratory Informatics Systems to deal with scientific instruments and data management Connecting instruments together, managing laboratory flow, etc. Bioinformatics Systems to deal with basic research data DNA, proteins, ‘molecular’ things What is bioinformatics?

  5. DNA sequencing is a common activity, so analysis of biological sequences has become critical to modern molecular biology and genetics. What kind of sequences? DNA Protein RNA Today we will learn about how these sequences are managed and used by biologists Sequence Analysis

  6. FASTA Format: >TITLE SEQUENCE 1 SEQUENCE1VSLLSNNLNKFDEGLALAHFVWIAPLQVALLMGLIWELLQASAFCGLGFLIVLALFQAGL >TITLE SEQUENCE 2 SEQUENCE2ASKKNPKLINALRRCFFWRFMFYGIFLYLGEVTKAVQPLLLGRIIASYDPDNKEERSIAIYLGIGLCLLFIVRTLLLHPAIFGLHHIGMQMRIAMFSLIYKKTLKLSSRVLDKISIGQLVSLLSNNLNKFDEGLALAHFVWIAPLQVALLMGLIWELLQASAFCGLGFLIVLALFQAG >TITLE SEQUENCE 3 AINKIQDFLQKQEYKTLEYNLTTTEVVMENVTAFWEEGFGELFEKAKQNNNNRKTSNGDDSLFFSNFSLLGTPVLKDINFKIERGQL Storing A Sequence On A Computer

  7. >gi|90421312|ref|NM_000492.3| Homo sapiens cystic fibrosis transmembrane conductance regulator (ATP-binding cassette sub-family C, member 7) (CFTR), mRNA AATTGGAAGCAAATGACATCACAGCAGGTCAGAGAAAAAGGGTTGAGCGGCAGGCACCCAGAGTAGTAGG TCTTTGGCATTAGGAGCTTGAGCCCAGACGGCCCTAGCAGGGACCCCAGCGCCCGAGAGACCATGCAGAG GTCGCCTCTGGAAAAGGCCAGCGTTGTCTCCAAACTTTTTTTCAGCTGGACCAGACCAATTTTGAGGAAA GGATACAGACAGCGCCTGGAATTGTCAGACATATACCAAATCCCTTCTGTTGATTCTGCTGACAATCTAT CTGAAAAATTGGAAAGAGAATGGGATAGAGAGCTGGCTTCAAAGAAAAATCCTAAACTCATTAATGCCCT TCGGCGATGTTTTTTCTGGAGATTTATGTTCTATGGAATCTTTTTATATTTAGGGGAAGTCACCAAAGCA GTACAGCCTCTCTTACTGGGAAGAATCATAGCTTCCTATGACCCGGATAACAAGGAGGAACGCTCTATCG CGATTTATCTAGGCATAGGCTTATGCCTTCTCTTTATTGTGAGGACACTGCTCCTACACCCAGCCATTTT TGGCCTTCATCACATTGGAATGCAGATGAGAATAGCTATGTTTAGTTTGATTTATAAGAAGACTTTAAAG CTGTCAAGCCGTGTTCTAGATAAAATAAGTATTGGACAACTTGTTAGTCTCCTTTCCAACAACCTGAACA AATTTGATGAAGGACTTGCATTGGCACATTTCGTGTGGATCGCTCCTTTGCAAGTGGCACTCCTCATGGG GCTAATCTGGGAGTTGTTACAGGCGTCTGCCTTCTGTGGACTTGGTTTCCTGATAGTCCTTGCCCTTTTT Actual FASTA File

  8. GenBank is the public domain database of biological (RNA, DNA, protein) sequences. If you sequence a novel sequence, you can submit it! Started at Los Alamos National Laboratory Moved to the National Center for Biotechnology Information in 1993 Stores almost every type of sequence possible GenBank

  9. You can access GenBank at http://www.ncbi.nlm.nih.gov Using GenBank

  10. Now bigger than 100,000,000,000 bases 240,000 named organisms More than 60 million records GenBank http://woldlab.caltech.edu/biohub/scipy2006/genbankgrowth.jpg

  11. GenBank is an uncurated mess… RefSeq – the reference sequence project. This is a subset of GenBank, and a curated set of sequences and annotations RefSeq IDs have the form of “Letter-Letter-Underscore-Number.” The letter-letter-underscore prefixes are: “NM_” - mRNA “NP_” - protein “NT_” – genomic (automated) “XP_” – genomic protein (automated), etc.. The need for quality: RefSeq

  12. It will come back to you in something like the GenBank file format. When you look for a sequence…

  13. LOCUS NM_000492 6132 bp mRNA linear PRI 07-OCT-2007 DEFINITION Homo sapiens cystic fibrosis transmembrane conductance regulator (ATP-binding cassette sub-family C, member 7) (CFTR), mRNA. ACCESSION NM_000492 VERSION NM_000492.3 GI:90421312 KEYWORDS . SOURCE Homo sapiens (human) ORGANISM Homo sapiens Eukaryota; Metazoa; Chordata; Craniata; Vertebrata; Euteleostomi; Mammalia; Eutheria; Euarchontoglires; Primates; Haplorrhini; Catarrhini; Hominidae; Homo. REFERENCE 1 (bases 1 to 6132) AUTHORS Pall,H., Zielenski,J., Jonas,M.M., DaSilva,D.A., Potvin,K.M., Yuan,X.W., Huang,Q. and Freedman,S.D. TITLE Primary sclerosing cholangitis in childhood is associated with abnormalities in cystic fibrosis-mediated chloride channel function JOURNAL J. Pediatr. 151 (3), 255-259 (2007) PUBMED 17719933 REMARK GeneRIF: There is a high prevalence of CFTR-mediated ion transport dysfunction in subjects with childhood primary sclerosing cholangitis. ......... What are the important elements of GenBank format?

  14. LOCUS NM_000492 6132 bp mRNA linear PRI 07-OCT-2007 DEFINITION Homo sapiens cystic fibrosis transmembrane conductance regulator (ATP-binding cassette sub-family C, member 7) (CFTR), mRNA. ACCESSION NM_000492 VERSION NM_000492.3 GI:90421312 KEYWORDS . SOURCE Homo sapiens (human) ORGANISM Homo sapiens Eukaryota; Metazoa; Chordata; Craniata; Vertebrata; Euteleostomi; Mammalia; Eutheria; Euarchontoglires; Primates; Haplorrhini; Catarrhini; Hominidae; Homo. REFERENCE 1 (bases 1 to 6132) AUTHORS Pall,H., Zielenski,J., Jonas,M.M., DaSilva,D.A., Potvin,K.M., Yuan,X.W., Huang,Q. and Freedman,S.D. TITLE Primary sclerosing cholangitis in childhood is associated with abnormalities in cystic fibrosis-mediated chloride channel function JOURNAL J. Pediatr. 151 (3), 255-259 (2007) PUBMED 17719933 REMARK GeneRIF: There is a high prevalence of CFTR-mediated ion transport dysfunction in subjects with childhood primary sclerosing cholangitis. ......... What are the important elements of GenBank format? Accession: NM_000492 Version: 3 GenBank ID:90421312 Symbol: CFTR

  15. The gene ID or symbol is not a cite-able sequence! Without the version, the RefSeq ID is not unique! The GI is always unique, if a change occurs to the entry, a new GI is issued. When in doubt, use the GI! What do you need to know?

  16. FEATURES Location/Qualifiers CDS 133..4575 /gene="CFTR" /protein_id="NP_000483.3" /db_xref="GI:90421313" /db_xref="CCDS:CCDS5773.1" /db_xref="GeneID:1080" /db_xref="HGNC:1884" /db_xref="HPRD:03883" /db_xref="MIM:602421" /translation="MQRSPLEKASVVSKLFFSWTRPILRKGYRQRLELSDIYQIPSVD SADNLSEKLEREWDRELASKKNPKLINALRRCFFWRFMFYGIFLYLGEVTKAVQPLLL GRIIASYDPDNKEERSIAIYLGIGLCLLFIVRTLLLHPAIFGLHHIGMQMRIAMFSLI YKKTLKLSSRVLDKISIGQLVSLLSNNLNKFDEGLALAHFVWIAPLQVALLMGLIWEL LQASAFCGLGFLIVLALFQAGLGRMMMKYRDQRAGKISERLVITSEMIENIQSVKAYC WEEAMEKMIENLRQTELKL...... GenBank files contain annotation

  17. FEATURES Location/Qualifiers CDS 133..4575 /gene="CFTR" /protein_id="NP_000483.3" /db_xref="GI:90421313" /db_xref="CCDS:CCDS5773.1" /db_xref="GeneID:1080" /db_xref="HGNC:1884" /db_xref="HPRD:03883" /db_xref="MIM:602421" /translation="MQRSPLEKASVVSKLFFSWTRPILRKGYRQRLELSDIYQIPSVD SADNLSEKLEREWDRELASKKNPKLINALRRCFFWRFMFYGIFLYLGEVTKAVQPLLL GRIIASYDPDNKEERSIAIYLGIGLCLLFIVRTLLLHPAIFGLHHIGMQMRIAMFSLI YKKTLKLSSRVLDKISIGQLVSLLSNNLNKFDEGLALAHFVWIAPLQVALLMGLIWEL LQASAFCGLGFLIVLALFQAGLGRMMMKYRDQRAGKISERLVITSEMIENIQSVKAYC WEEAMEKMIENLRQTELKL...... Note protein translation is present! There is also another protein record for this entry (NP_000483.3)

  18. Why so many IDs? Biology! Gene Symbol (CFTR) GeneID – Unique to organism (1080) • Transcript ID -2 etc • Transcript 2 version1 • Transcript 2 version2 • Transcript ID -1 • Transcript 1 version1 • Transcript 1 version2 • Transcript 1 version3

  19. Gene symbol/name is the official alphanumeric name of a gene It is defined by Human Genome Organization (HUGO) It generally does not reliably tell you the species or describe a transcript Example: CFTR, ER, AR, BRCA1 Gene Identifier A gene identifier, generally a number, allows you to uniquely link to the gene and species Common identifiers: Locus Link or Gene – NIH identifiers Mouse Genome Informatics IDs (MGI) – Mouse identification Ensembl IDs – European gene identifiers Transcript (the mRNA product of gene) RefSeq GenBank Protein Next topic! The heirarchy of names

  20. The most popular protein database is called Uniprot Uniprot is a superset of Swiss-Prot, at http://uniprot.org Name is CFTR_HUMAN, Accession is Letter+Number (P13569) Protein databases

  21. When publishing genes or proteins: When publishing information about a gene, use the HUGO approved symbol. If publishing information specifically about the gene, put the gene ID (preferably the locus link or entrez gene id). When publishing information about a specific transcript of a gene, use the RefSeq ID. If including specifics about the sequence, include the version Uniprot is a great way to reference proteins. Some simple rules for publication

  22. We now know how to find and cite a sequence from a database We know that genomes, genes, transcripts and proteins are treated differently How do we use sequences to do things like predict function? Sequence analysis! What next?

  23. Given two sequences of letters, and a scoring scheme for evaluating matching letters, find the optimal pairing of letters from one sequence to letters of the other sequence. (The sequence analysis slides were adapted from Russ Altman’s bioinformatics course, Stanford University) Sequence analysis is based on the sequence alignment! Align: THIS IS A RATHER LONGER SENTENCE THAN THE NEXT.THIS IS A SHORT SENTENCE. THIS IS A RATHER LONGER SENTENCE THAN THE NEXT. THIS IS A ------SHORT-- SENTENCE--------------. OR THIS IS A SHORT---------SENTENCE--------------.

  24. DNA (4 letter alphabet + gap) TTGACAC TTTACAC Proteins (20 letter alphabet + gap) RKVA--GMAKPNM RKIAVAAASKPAV Aligning biological sequences

  25. Given: 2 sequences scoring system for evaluating match(or mismatch) of two characters penalty function for gaps in sequences Produce: Optimal pairing of sequences that retains the order of characters introduces gaps maximizes total score Statement of Problem

  26. Lots of sequences with unknown structure and function. A few sequences with known structure and function. If they align, they are similar, maybe due to common descent. If they are similar, then they might have same structure or function. If one of them has known structure/function, then alignment to the other yields insight about how the structure or function works. Why align sequences?

  27. Alignments have parameters Exact Matches OK, Inexact Costly, Gaps cheap. This is a rather longer sentence than the next. This is a ------------- sentence--------------. OR This is a *rather longer sentence than the next. This is a s---h----o---rtsentence--------------. Exact Matches OK, Inexact Moderate, Gaps cheap. This is a rather longer sentence than the next. This is a --short------ sentence--------------. Exact Matches cheap, Inexact cheap, Gaps expensive. This is a rather longer sentence than the next. This is a short sentence.----------------------

  28. Pairwise alignment (two at a time) is much easier than multiple alignment (N at a time). To be discussed later. Multiple alignment This is a rather longer sentence than the next. This is a short sentence. This is the next sentence. Rather long is the next concept. Rather longer than what is the next concept.

  29. Multiple Alignment This is a rather longer sentence than --------the next---------. This is a short-------- sentence-------------------------------. This is --------------------------------------the next sentence. ----------Rather long is --------the------------- next concept-. ----------Rather longer ---------than what is the next concept-.

  30. Put one sequence along the top row of a matrix. Put the other sequence along the left column of the matrix. Plot a dot everytime there is a match between an element of row sequence and an element of the column sequence. Diagonal lines indicate areas of match. Dot Plots To Visualize Sequence Similarity

  31. Rely on visual analysis Difficult to find optimal alignments Need scoring schemes more sophisticated that “identical match” Difficult to estimate significance of alignments Problems with dot matrices

  32. Sequence alignment generally uses an algorithm called dynamic programming Dynamic programming allows for fast, optimal alignment of sequences Very informally, dynamic programming finds the best path through a dot plot. The Dynamic Programming Algorithm

  33. The degree of match between two letters can be represented in a matrix. Area of active research Changing matrix changes alignments context-specific matching information theoretic interpretation of scores modeling evolution with different matrices Substitution Matrices

  34. A Sample Match Matrix for the amino acids (PAM-250).

  35. 1. Manually align protein structures (or, more risky, sequences) 2. Look for frequency of amino acid substitutions at structurally nearly constant sites. 3. Entry ~ log(freq(observed)/freq(expected)) + —> More likely than random 0 —> At random base rate - —> Less likely than random Where do matrices come from?

  36. Global alignment: find alignment in which total score is highest, perhaps at expense of areas of great local similarity. Local alignment: find alignment in which the highest scoring subsequences are identified, at the expense of the overall score. Global vs. Local Alignment

  37. Two types of alignment approaches for DNA and Protein Global alignments Originally described with the Needleman-Wunsch algorithm [1970] using dynamic programming Optimally finds an alignment between two sequences ClustalW is software for global alignment Local alignments Originally described with the Smith-Waterman algorithm [1981], again using dynamic programming Optimally finds subsequence alignments Sequence alignment

  38. Searching databases Assembling genome sequences Annotation transfer: similar sequences implies similar function Constructing multiple sequence alignments of groups of sequences Understanding evolutionary relationships between sequences DNA/RNA hybridization probe design Peptide and protein identification in proteomic experiments For more detail on sequence alignment algorithms and their statistics see: Biological Sequence Analysis: Probabilistic models of proteins and nucleic acids by Durbin, Eddy, Krogh and Mitchison (Cambridge Press) Uses of sequence alignment

  39. TTGACACCCTCCCAATTGTA ACCCCAGGCTTTACACAT TTGACACCCTCCCAATTGTA--- |||| || | -----ACCCCAGGCTTTACACAT ---------TTGACACCCTCCCAATTGTA || |||| ACCCCAGGCTTTACACAT----------- (from Gribskov/Devereaux page 133) Global vs. Local Alignment

  40. NOT: Are these two sequences generally the same? (Global question) BUT: Do these two sequences contain high scoring subsequences? (Local question) Local similarities may occur in sequences with different structure or function that share common substructure/subfunction. MOTIFS… Local alignment asks different question than Global alignment.

  41. ClustalW, http://www.ebi.ac.uk/Tools/clustalw/index.html Enter sequences in FASTA format! Parameters: penalties for gaps, substitution matrix, etc. Obtaining a Global Alignment

  42. We will learn about local alignment tomorrow. Local Alignment?

  43. Biological sequences are stored in databases with annotations relevant to their function The primary sequence database is GenBank, hosted at the National Institutes of Health The providers of these databases have come up with database IDs that we should use when writing papers for genes, transcripts and proteins Sequence analysis allows us to compare two sequences through sequence alignments which come in two flavors, local and global Summary

More Related