Dot plot (bioinformatics)
One way to visualize the similarity between two protein or nucleic acid sequences is to use a similarity matrix, known as a dot plot. These were introduced by Gibbs and McIntyre in 1970 and are two-dimensional matrices that have the sequences of the proteins being compared along the vertical and horizontal axes. For a simple visual representation of the similarity between two sequences, individual cells in the matrix can be shaded black if residues are identical, so that matching sequence segments appear as runs of diagonal lines across the matrix.
Some idea of the similarity of the two sequences can be gleaned from the number and length of matching segments shown in the matrix. Identical proteins will obviously have a diagonal line in the center of the matrix. Insertions and deletions between sequences give rise to disruptions in this diagonal. Regions of local similarity or repetitive sequences give rise to further diagonal matches in addition to the central diagonal. Because of the limited protein alphabet, many matching sequence segments may simply have arisen by chance. One way of reducing this noise is to only shade runs or 'tuples' of residues, e.g. a tuple of 3 corresponds to three residues in a row. This is effective because the probability of matching three residues in a row by chance is much lower than single-residue matches. It can be seen from Figures 3.3h,c that the number of diagonal runs in the matrix has been considerably reduced by looking for 2-tuples or 3-tuples.
Dot plots compare two sequences by organizing one sequence on the x-axis, and another on the y-axis, of a plot. When the residues of both sequences match at the same location on the plot, a dot is drawn at the corresponding position. Note, that the sequences can be written backwards or forwards, however the sequences on both axes must be written in the same direction. Also note, that the direction of the sequences on the axes will determine the direction of the line on the dot plot. Once the dots have been plotted, they will combine to form lines. The closeness of the sequences in similarity will determine how close the diagonal line is to what a graph showing a curve demonstrating a direct relationship is. This relationship is affected by certain sequence features such as frame shifts, direct repeats, and inverted repeats. Frame shifts include insertions, deletions, and mutations. The presence of one of these features, or the presence of multiple features, will cause for multiple lines to be plotted in a various possibility of configurations, depending on the features present in the sequences. A feature that will cause a very different result on the dot plot is the presence of low-complexity region/regions. Low-complexity regions are regions in the sequence with only a few amino acids, which in turn, causes redundancy within that small or limited region. These regions are typically found around the diagonal, and may or may not have a square in the middle of the dot plot.
- Gibbs, Adrian J.; McIntyre, George A. (1970). "The Diagram, a Method for Comparing Sequences. Its Use with Amino Acid and Nucleotide Sequences". Eur. J. Biochem 16: 1–11. doi:10.1111/j.1432-1033.1970.tb01046.x.
Software to create plots
- Genomdiff – An open source Java dot plot program for viruses.
- Gepard - Dot plot tool suitable for even genome scale.
- ANACON – Contact analysis of dot plots.
- General introduction to dot plots with example algorithms and a software tool to create small and medium size dot plots.
- Dotlet – Provides a program allowing you to construct a dot plot with your own sequences.
- UGENE Dot Plot viewer – Opensource dot plot visualizer.
- seqinr - R package to generate dot plots.
- dotplot - R package to rapidly generate dot plots as either traditional or ggplot graphics.
- dotmatcher - Web tool to generate dot plots.
- Dotter - Stand alone program to generate dot plots.