Gene Prediction: Similarity-Based Approaches in Bioinformatics
Gene prediction in bioinformatics involves predicting gene locations in a genome using different approaches like statistical methods and similarity-based approaches. The similarity-based approach uses known genes as a template to predict unknown genes in newly sequenced DNA fragments. This method involves finding local similarities between genomic and target protein sequences and selecting non-overlapping substrings with high similarity as putative exon structures. Other techniques discussed include the Exon Chaining Problem and Reverse Translation for inferring coding DNA from known proteins.
Download Presentation
![](/assets/img/so-down.gif)
Please find below an Image/Link to download the presentation.
The content on the website is provided AS IS for your information and personal use only. It may not be sold, licensed, or shared on other websites without obtaining consent from the author. Download presentation by click this link. If you encounter any issues during the download, it is possible that the publisher has removed the file from their server.
E N D
Presentation Transcript
Dynamic Programming II Dynamic Programming II Gene Prediction: Similarity Gene Prediction: Similarity- -Based Approaches The idea of similarity-based approach to gene prediction Exon Chaining Problem Spliced Alignment Problem Based Approaches
An Introduction to Bioinformatics Algorithms www.bioalgorithms.info Gene Prediction Computational problem of predicting the locations of genes in a genome given only the genomic DNA sequence (gene is broken into pieces called as exons that are separated by junk NDA/introns). Two approaches 1. Statistical Approaches : limited success since it tends to match frequently in the genome at non-splice sites using a profile describing the propensities of different nucleotides to occur at different positions. use stop codons, TAA, TAG, TGA, or start codon ATG AG and GT on the left- and right-hand sides of an exon are highly conserved. 2. Similarity-Based Approaches: using previously sequenced genes and their protein products as a template for the recognition of unknown genes in newly sequenced DNA fragments. Find all local similarities between the genomic sequence and the target protein sequence (using local alignment algorithm) Shorten or extend the substrings with high similarity such that they start at AG and end at GT. The resulting set may contain overlapping substrings. Choose the best subset of non-overlapping substrings as a putative exon structure.
An Introduction to Bioinformatics Algorithms www.bioalgorithms.info Similarity-Based Approach to Gene Prediction Genes in different organisms are similar The similarity-based approach uses known genes in one genome to predict (unknown) genes in another genome Problem: Given a known gene and an unannotated genome sequence, find a set of substrings of the genomic sequence whose concatenation best fits the gene
An Introduction to Bioinformatics Algorithms Comparing Genes in Two Genomes www.bioalgorithms.info Small islands of similarity corresponding to similarities between exons
An Introduction to Bioinformatics Algorithms www.bioalgorithms.info Reverse Translation Given a known protein, find a gene in the genome which codes for it One might infer the coding DNA of the given protein by reversing the translation process Inexact: amino acids map to > 1 codon UUA, UUG, CUU, CUC, CUA, CUG: leucine This problem is essentially reduced to an alignment problem
An Introduction to Bioinformatics Algorithms Comparing Genomic DNA Against mRNA www.bioalgorithms.info (codon sequence) mRNA { { { { { exon1 intron1 exon2 intron2 exon3 Portion of genome
An Introduction to Bioinformatics Algorithms www.bioalgorithms.info Using Similarities to Find the Exon Structure The known frog gene is aligned to different locations in the human genome Find the best path to reveal the exon structure of human gene Frog Gene (known) Human Genome
An Introduction to Bioinformatics Algorithms Finding Local Alignments www.bioalgorithms.info Use local alignments to find all islands of similarity Frog Genes (known) Human Genome
An Introduction to Bioinformatics Algorithms Chaining Local Alignments www.bioalgorithms.info Find substrings that match a given gene sequence (candidate exons) Define a candidate exons as (l, r, w) (left, right, weight defined as score of local alignment) Look for a maximum chain of substrings Chain: a set of non-overlapping nonadjacent intervals.
An Introduction to Bioinformatics Algorithms Exon Chaining Problem www.bioalgorithms.info 5 5 15 9 11 4 3 0 2 3 5 6 11 13 16 20 25 27 28 30 32 Locate the beginning and end of each interval (2n points) Find the best path
An Introduction to Bioinformatics Algorithms Exon Chaining ProblemP Given a set of putative exons, where each exon is represented by (l,r,w), l and r are the left- and right-hand positions, and w is the weight reflecting the likelihood that this interval is an exon (e.g., the local alignment score), find a maximum set of non- overlapping putative exons. Input: A set of weighted intervals (putative exons) Output: A maximum chain of intervals from this set. 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 www.bioalgorithms.info 3 1 0 4 5 12 6 7 10 5 10 12 3 6 4 7 1 0 0 3 3 5 5 5 9 9 10 10 15 15 15 17 17 17 21
An Introduction to Bioinformatics Algorithms ExonChaining(G,n) for i := 1 to 2n s[i] :=0; l[i] = 0; for i : = 1 to 2n if vertex v[i] in G corresponds to the right end of an interval I j := index of vertex for left end of the interval I w := weight of the interval I s[i]:= max{s[j] + w, s[i-1]} if (s[i] = s[j]+w) then l[i]= j; else s[i] := s[i-1] v www.bioalgorithms.info PrintChain(l, m) If m = 0 return Else if ( l[m] 0 ) PrintChain(l, l[m]) Print ( l[m], m ) else PrintChain (l, m-1) 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 3 1 0 4 5 13 I 6 7 10 5 10 13 3 6 4 7 1 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 0 0 3 3 5 5 5 9 9 10 10 15 15 15 17 17 18 21 s l 0 0 2 0 1 0 0 4 0 9 0 5 0 0 11 0 7 16
An Introduction to Bioinformatics Algorithms Exon Chaining: Deficiencies www.bioalgorithms.info Poor definition of the putative exon endpoints Optimal chain of intervals may not correspond to any valid alignment First interval may correspond to a suffix, whereas second interval may correspond to a prefix Combination of such intervals is not a valid alignment
An Introduction to Bioinformatics Algorithms Infeasible Chains www.bioalgorithms.info Red local similarities form two non -overlapping intervals but do not form a valid global alignment Frog Genes (known) Human Genome
An Introduction to Bioinformatics Algorithms www.bioalgorithms.info Gene Prediction Analogy: Selecting Putative Exons The cell carries DNA as a blueprint for producing proteins, like a manufacturer carries a blueprint for producing a car.
An Introduction to Bioinformatics Algorithms Using Blueprint www.bioalgorithms.info
An Introduction to Bioinformatics Algorithms Assembling Putative Exons www.bioalgorithms.info
An Introduction to Bioinformatics Algorithms Still Assembling Putative Exons www.bioalgorithms.info
An Introduction to Bioinformatics Algorithms Spliced Alignment www.bioalgorithms.info Mikhail Gelfand and colleagues proposed a spliced alignment approach of using a protein within one genome to reconstruct the exon-intron structure of a (related) gene in another genome. Begins by selecting either all putative exons between potential acceptor and donor sites or by finding all substrings similar to the target protein (as in the Exon Chaining Problem). This set is further filtered in a such a way that attempt to retain all true exons, with some false ones.
An Introduction to Bioinformatics Algorithms www.bioalgorithms.info Spliced Alignment Problem: Formulation Goal: Find a chain of blocks in a genomic sequence that best fits a target sequence Input: Genomic sequences G, target sequence T, and a set of candidate exons B. Output: A chain of exons such that the global alignment score between * and T is maximum among all chains of blocks from B. * - concatenation of all exons from chain
An Introduction to Bioinformatics Algorithms Lewis Carroll Example www.bioalgorithms.info 1 2 3 4 1 2 4 3
An Introduction to Bioinformatics Algorithms Spliced Alignment: Idea www.bioalgorithms.info Compute the best alignment between i-prefix of genomic sequence G and j-prefix of target T: S(i,j) But what is i-prefix of G? There may be a few i-prefixesof G depending on which block B we are in. Compute the best alignment between i-prefix of genomic sequence G and j-prefix of target T under the assumption that the alignment uses the block B at position i S(i,j,B)
An Introduction to Bioinformatics Algorithms Spliced Alignment Recurrence If i is not the starting vertex of block B: S(i, j, B) = max { S(i 1, j, B) indel penalty S(i, j 1, B) indel penalty S(i 1, j 1, B) + (gi, tj) } If i is the starting vertex of block B: S(i, j, B) = max { S(i, j 1, B) indel penalty maxall blocks B preceding block B S(end(B ), j, B ) indel penalty maxall blocks B preceding block B S(end(B ), j 1, B ) + (gi, tj) } T www.bioalgorithms.info j j-1 G End(B ) i-1 i S(i,j,B)
An Introduction to Bioinformatics Algorithms Spliced Alignment Solution www.bioalgorithms.info After computing the three-dimensional table S(i, j, B), the score of the optimal spliced alignment is: maxall blocksBS(end(B), length(T), B)
An Introduction to Bioinformatics Algorithms Spliced Alignment: Complications www.bioalgorithms.info Considering multiple i-prefixes leads to slow down. running time: O(mn2|B|) where m is the target length, n is the genomic sequence length and|B|is the number of blocks. A mosaic effect: short exons are easily combined to fit any target protein
An Introduction to Bioinformatics Algorithms Spliced Alignment: Speedup www.bioalgorithms.info
An Introduction to Bioinformatics Algorithms Spliced Alignment: Speedup www.bioalgorithms.info
An Introduction to Bioinformatics Algorithms Spliced Alignment: Speedup www.bioalgorithms.info P(i,j)=maxall blocks B preceding position i S(end(B), j, B)
An Introduction to Bioinformatics Algorithms Exon Chaining vs Spliced Alignment www.bioalgorithms.info In Spliced Alignment, every path spells out string obtained by concatenation of labels of its edges. The weight of the path is defined as optimal alignment score between concatenated labels (blocks) and target sequence Defines weight of entire path in graph, but not the weights for individual edges. Exon Chaining assumes the positions and weights of exons are pre-defined
An Introduction to Bioinformatics Algorithms www.bioalgorithms.info Gene Prediction Tools GENSCAN/Genome Scan TwinScan Glimmer GenMark
An Introduction to Bioinformatics Algorithms www.bioalgorithms.info Gene Prediction: Aligning Genome vs. Genome Align entire human and mouse genomes Predict genes in both sequences simultaneously as chains of aligned blocks (exons) This approach does not assume any annotation of either human or mouse genes.
An Introduction to Bioinformatics Algorithms The GENSCAN Algorithm www.bioalgorithms.info Algorithm is based on probabilistic model of gene structure similar to Hidden Markov Models (HMMs). GENSCAN uses a training set in order to estimate the HMM parameters, then the algorithm returns the exon structure using maximum likelihood approach standard to many HMM algorithms (Viterbi algorithm). Biological input: Codon bias in coding regions, gene structure (start and stop codons, typical exon and intron length, presence of promoters, presence of genes on both strands, etc) Covers cases where input sequence contains no gene, partial gene, complete gene, multiple genes.
An Introduction to Bioinformatics Algorithms GENSCAN Limitations Does not use similarity search to predict genes. Does not address alternative splicing. Could combine two exons from consecutive genes together www.bioalgorithms.info
An Introduction to Bioinformatics Algorithms GenomeScan www.bioalgorithms.info Incorporates similarity information into GENSCAN: predicts gene structure which corresponds to maximum probability conditional on similarity information Algorithm is a combination of two sources of information Probabilistic models of exons-introns Sequence similarity information
An Introduction to Bioinformatics Algorithms TwinScan www.bioalgorithms.info Aligns two sequences and marks each base as gap ( - ), mismatch (:), match (|), resulting in a new alphabet of 12 letters: {A-, A:, A |, C-, C:, C |, G-, G:, G |, T-, T:, T|}. Run Viterbi algorithm using emissions ek(b) where b {A-, A:, A|, , T|}. http://www.standford.edu/class/cs262/ Spring2003/Notes/ln10.pdf
An Introduction to Bioinformatics Algorithms TwinScan(cont d) www.bioalgorithms.info The emission probabilities are estimated from from human/mouse gene pairs. Ex. eI(x|) < eE(x|) since matches are favored in exons, and eI(x-) > eE(x-) since gaps (as well as mismatches) are favored in introns. Compensates for dominant occurrence of poly-A region in introns
An Introduction to Bioinformatics Algorithms Glimmer www.bioalgorithms.info Gene Locator and Interpolated Markov ModelER Finds genes in bacterial DNA Uses interpolated Markov Models
An Introduction to Bioinformatics Algorithms The Glimmer Algorithm www.bioalgorithms.info Made of 2 programs BuildIMM Takes sequences as input and outputs the Interpolated Markov Models (IMMs) Glimmer Takes IMMs and outputs all candidate genes Automatically resolves overlapping genes by choosing one, hence limited Marks suspected to truly overlap genes for closer inspection by user
An Introduction to Bioinformatics Algorithms GenMark www.bioalgorithms.info Based on non-stationary Markov chain models Results displayed graphically with coding vs. noncoding probability dependent on position in nucleotide sequence
An Introduction to Bioinformatics Algorithms www.bioalgorithms.info Homework Assignment (1) Solve Exon Chaining Problem in slide 11 when the weight for interval (6,12) change to 18 and weight for interval change to (13,14) = 4 (2) Implement Exon Chaining Algorithm in slide 12. Print out the tables s, l and exon chain for the problem in slide 11.