270 likes | 446 Views
Clustering Method for Repeat Analysis in DNA sequences. Gnana Sundar Rajendiran Joyesh Mishra Rishi Mishra FALL 2008 Bioinformatics. Abstract. Implement a Proposed Clustering Technique by Volfovsky et al. for Repeat Analysis
E N D
Clustering Method for Repeat Analysis in DNA sequences Gnana Sundar RajendiranJoyeshMishraRishiMishra FALL 2008 Bioinformatics
Abstract Implement a Proposed Clustering Technique by Volfovsky et al. for Repeat Analysis Distribute Merging Repeats into Classes/Clusters based on Similarity Measures Why analyze repeats?
Algorithm Description • Selection using RepeatMatch • Steps • Preprocessing • Merging and Repeat Map Generation • Classification • BLAST searches & further merging of Clusters • Results
RepeatMatch Uses Suffix Tree algorithm to determine all the exact repeats in a given sequence It is a part of the MUMmer package used for the rapid alignment of very large DNA and amino acid sequences.
Example… • Forward & Reverse Complement Repeats
Definitions • Exact Repeat: Subsequence that occurs in DNA at least twice. • Exact repeat is represented by pair of co-ordinates (A1,A2) delimiting its location in the genome sequence and by the repeat length. • Maximal Repeat: Repeat that can not be extended in either direction without incurring a mismatch • Initial Repeat Set: Set of repeats chosen initially, from which the repeat classes will be constructed
Preprocessing Output from RepeatMatch is used to partition the original genome sequence Each partition point has a reference to the pair co-ordinates (A1, A2) and repeat length l For each repeat starting at co-ordinates A1 and A2, with length l, this list will include both (A1, A2, l) and (A2, A1,l)
Merging & Repeat Map Generation • This procedure works by repeatedly working together two exact repeats that either overlap or that occur within a limited distance(a gap) of each other. • Significant subsequences of the new merging repeats appear at least twice in the genome sequence. • Merging with Gap • Merging with Overlap
Merging with Gap Given 2 partition points: p1 = (A1, A2, lA) and p2 = (B1, B2, lB) [A1 < B1] Compute the distance between the non-overlapping repeats as d(p1,p2) = max(0, B1 - A1 - lA + 1) Given a maximum Gap size G > 0 Sequences corresponding to p1 and p2 are merged if d(p1, p2) < G
Merging with Overlap Merges sequences which are partially identical Overlap of 2 sequences is denoted as: o(p1, p2) = max(o, A1 + lA – B1 + 1) for A1<B1 Given a minimum overlap proportion op, where 0 <= op <= 1, repeat points (A1, A2, lA) and (B1, B2, lB) are merged if at least one of the four repeats has overlap satisfying o(p1, p2) > op min(lA, lB) op : Interpreted as a fraction of the shorter of the 2 repeats. E.g. if op = 0.75, the two overlapped sequences are merged if the length of their overlap is at least 75% of the length of the shorter sequence.
Output: Merging Repeat The new sequence is defined as merging repeat with starting position M = A1 and with length lM = max(A1 + lA, B1 + lB) – A1 A data structure stores with each merging repeat its start co-ordinate, the length(nM), and a list of references to itself and to other repeats.
Classification If a merging repeat has at least one reference in common with one another, then they belong to the same class If a merging repeat has references that belong to multiple distinct classes, then those classes are combined into one. If a merging repeat contains no reference to an existing class, then the merging repeat forms a new class.
BLAST searches and further merging • To merge similar but non-exact repeats • Search all merging repeats against all others • A local BLAST database is created with all repeat sequences using formatdb • Classes are merged if any of their underlying sequences have a BLAST E-value less than a User-Specified Threshold when compared to any sequence in another class • If a class appears in multiple similarity pairs, all these similar classes are merged with the original class
What’s Next? The results collected help to cluster similar repeats and generate a database. This DB, in future, can be used to find Similar Repeats faster. Also, if on future classification, we find common traits among repeats, they can be analyzed only for that cluster alone.
Resources • Softwares Used: • RepeatMatch (MUMmer 3.0) • BLAST (formatdb, blastall) • MS Excel 2007 • Programmming Language Used: C++, UNIX Shell Script • Source for Project & Presentation • www.cise.ufl.edu/~jmishra/bioinformatics.html
References • Volfovsky N., Haas Brian J. and Salzberg Steven L. A clustering method for repeat analysis in DNA sequences, 2001 • Mummer software : http://mummer.sourceforge.net/ • http://www.ncbi.nlm.nih.gov/staff/tao/URLAPI/ • blastall, formatdb • www.ncbi.nlm.nih.gov/BLAST/
References http://www.tigr.org http://www.animalgenome.org/blast/docs/blast_faq.html Of course, Wikipedia … Data obtained from: ftp://ftp.ncbi.nih.gov/genomes/Bacteria/