- Methodology article
- Open Access
Biological sequences as pictures – a generic two dimensional solution for iterated maps
BMC Bioinformatics volume 10, Article number: 100 (2009)
Representing symbolic sequences graphically using iterated maps has enjoyed an enduring popularity since it was first proposed in Jeffrey 1990 as chaos game representation (CGR). The usefulness of this representation goes beyond the convenience of a scale independent representation, it provides a variable memory length representation of transition. This includes the representation of succession with non-integer order, which comes with the promise of generalizing Markovian formalisms. The original proposal targeted genomic sequences only but since then several generalizations have been proposed, many specifically designed to handle protein data.
The challenge of a general solution is that of deriving a bijective transformation of symbolic sequences into bi-dimensional planes. More specifically, it requires the regular fractal nesting of polygons. A first attempt at a general solution was proposed by Fiser 1994 by using non-overlapping circles that contain the polygons. This was used as a starting point to identify a more efficient solution where the encapsulating circles can overlap without the same happening for the sequence maps which are circumscribed to fractal polygon domains.
We identified the optimal inscribed packing solution for iterated maps of any Biological sequence, indeed of any symbolic sequence. The new solution maintains the prized bijective mapping property and includes the Sierpinski triangle and the CGR square as particular solutions of the more encompassing formulation.
The use of iterative functions to represent nucleotide sequence was originally proposed in 1990 with the designation of Chaos Game Representation, CGR . Since then, the CGR technique evolved from being a graphic representation technique into becoming a platform for pattern recognition [2–9], screening entropic properties [2, 10–12], and finally into a generalization of Markov Transition tables [13, 14]. That conclusion was further generalized for any symbolic sequence  where a dissimilarity metric was also proposed and routes for efficient implementation were established . The emerging interest in alignment-free sequence analysis techniques , in particular for application to proteomic sequences, raises the prospect of a wider use of CGR and CGR-related techniques [18, 19] capable of exploratory sequence analysis of whole genomes [20–22]. We have subsequently examined the advantage of using non-genomic word-statistics (integer order) for Biological sequence analysis  with an application to the SCOP protein database . The strengths of that approach suggest that even more interesting results would be achieved by order-free analysis. However, the potential of hyper-dimensional generalization of CGR such as the use of Universal Sequence Maps, USM , are hard to convey and realize in the absence of a 2D projection that retains the bijective mapping property of this technique. This conclusion is particularly clear in a recent report exploring CGR application to aminoacid sequences . In that work, in order to benefit from the scale independency of Chaos Game applied to proteins, the reverse encoding of the aminoacid sequences back to the 4 unit alphabet nucleotide sequences is the solution of choice. Furthermore, as detailed in the derivation of sequence similarity metrics and density kernel functions [13–15], CGR is most useful as a bijective mapping of a sequence of symbols into a numeric vector of the same length, rather than as a technique to compress sequences into individual points. Accordingly, the work reported here is driven by the need to advance the identification of a 2D projection of the vector of CGR positions that is still applicable to longer, often non-genomic, alphabets.
Either as a dimension reduction technique to represent hyper-dimensional CGR  or an extension of 2D representation beyond nucleotide sequences, the generalization of bijective iterated maps is desirable for scale-free visualization and analysis of Biological sequences. Three different approaches have been proposed to represent the iterative function results for sequences with more than 4 unique type of units. Those efforts have often been specifically driven towards extending to protein sequences the advantages of CGR of nucleotide sequences.
A possible solution is to keep the quadrangular representation for each unit and arrange them in a tabular format . However, this solution sacrifices the regularity of the CGR solution where changing the position of the elements of the alphabet (ACGT) corresponds to a simple translation of the coordinates of the map. In order to preserve equivalence between units of the alphabet, the iterated dividing ratios would have to have the same value for all elements of the larger protein alphabet.
To satisfy the desirable equivalence between all elements of the alphabet (all 20 aminoacids in the case of proteins) one could simply keep the 1/2 ratio of the original CGR  independently from the length of the alphabet. This solution has in fact been explored  by morphing the unit square into a polygon with as many vertices as units in the alphabet. However, this second option sacrifices the bijective mapping property for alphabets with more than 4 units. Without this property, the graphical projection becomes just a visualization technique from which it is no longer possible to completely recover the sequence composition.
Finally, the third possibility, chronologically the first extension to be proposed, is to adjust the iterative step such that the prized bijective property is preserved , which enables recovery of sequence from map position. However, in that work, the iterative step was adjusted to inscribe the circles where the polygons are inscribed, not the polygons themselves. As a consequence, neither the CGR nor the iconic three alphabet unit map, the Spierpinsky triangle, appear as particular solutions of this formulation.
The study reported here takes the third of these approaches  as the starting point and then seeks to identify a new rule that maximizes the packing density of the projection to the point where the polygons are inscribed directly within each other. In addition to providing a more efficient visualization of the sequence, the specific goal of this study if to obtain a solution that also generalizes the 3 alphabet unit solution by the Spierpinsky triangle and the 4 alphabet unit solution of genomic CGR.
The algorithms were identified using MathCad (MathSoft Inc) and Matlab (Mathworks Inc) programming environment was used to deploy them. Matlab was also used to produce the graphic displays presented in this report. The pseudo-random number generator of this programming environment was used to produce the 10,000-long random sequences used for the simulations. This manuscript is also accompanied by a graphic user interface (Figure 1) to a Matlab application where arbitrary symbolic sequences can be processed using different dividing ratios, including the generic solution identified by this study.
The main goal of the work described here is achieved with the identification of the dividing ratio in Equation 2 for which maximum packing of nexted, non-ovelapping polygons is obtained. The identification of this equation was achieved empirically with the assistance of MathCad's symbolic processing engine. Specifically, different n-polygons were inscribed in a circle of unitary radius and, using basic trigonometry: the dividing ratio was obtained as the fraction of the distance between any two edges that would identify the perimeter of the target polygon. The observation that the dividing ratio between any two edges is a constant value for each n-polygon was itself an empirical observation. The collection of solutions obtained was symbolically processed using MathCad's equation simplification tool to generate the solution in Equation 2.
The starting point for the attempt to find a generic formulation of CGR that is applicable to all alphabets is twofold. Firstly and foremost there is the original CGR formulation for nucleotide sequences , using a dividing ratio of 1/2 to fill a unit square domain homogeneously. Secondly, there is the proposition that there is a function of the number of edges in a regular polygon, n, that will yield a maximum dividing ratio, s, where the bijective property is still preserved. Fiser et. al  proposed Equation 1 as a suitable description of that dependency.s = (1 + sin (π/n))-1
That solution was designed to preserve the bijective mapping property by avoiding overlap between the circles that enclose the individual projections. However, that solution does actually not produce the full packing of the original CGR for n = 4, as noted by the proponents themselves in the original report (color plate 1 in ), and graphically recalled in Figure 2. It also sets a looser packing for n = 3, where the ratio, s, is 0.54, instead of 1/2 as necessary to produce the Sierpinki's triangle.
Here we seek to identify an alternative general formulation for the chaos game dividing ratio parameter, s, that will produce the reference structures as special solutions. The relevant reference structures include three maps that use dividing rations of 1/2, for n = 2, 3 and 4. The first is the 1D CGR, i.e. n = 2, used to repeatedly generate hyper-dimensional representations in the Unified Sequence Map technique . The second reference representation is the emblematic Sierpinski triangle , for n = 3. Finally and most important, the original formulation of CGR for n = 4 , where the unit square is covered homogeneously for a dividing ration of 1/2. We have achieved this goal by following a deductive procedure that, while fastidious, was derived using basic trigonometry. The proposed solution was found by noting that a tighter packing would be achieved if instead of inscribing the circles defined by the polygonic projections areas, the inscription was determined for the polygons themselves. The final solution (see Methods for its derivation) is presented in Equation 2:
The claim that Equation 2 is an accurate generalization of Chaos Game for any alphabet is best first illustrated graphically. In the right hand column of Figure 3, the projection areas defined by the proposed solution are produced for simulated sequences made using alphabets of increasing length. This includes the reference lengths of 2, 3 and 4 where the desired special solutions are correctly defined. For comparison, the formulation that served as starting point for this exercise, Equation 1, is illustrated on the left hand column of the same figure.
It is important to note (see Methods for description of computational derivation) that similarly to the reference circle packing solution by Fiser et. al (Equation 1), the polygon packing achieved by Equation 2 is also an heuristic solution. This is not a rare situation when it comes to fractal geometry  where the numerical solutions for a novel graphical configuration are often the starting point, rather than the conclusion, of the deductive process . It nevertheless suggests that further analysis of this solution is needed to uncover simpler and more meaningful patterns.
The relative gains of the proposed formulation for different alphabets are assessed in Figure 4. It is noted that the two formulations converge as the alphabet length increases but not monotonically. Instead Figure 4 shows that equivalence between the two formulations occur at regular increments of 4 in the number of vertices of the projection, starting with binary vocabularies (n = 2, 6, 10, 14, ...). The explanation for this observation becomes immediately clear when the two solutions are compared graphically, in Figure 3. That figure shows that for the equivalent solutions the vertices shared between distinct polygon areas occur at the exact position where the two involving circles intercept.
As overviewed in the Background section, the enduring popularity of iterated maps to represent Biological sequences has gradually expanded from a convenient graphical illustration into an efficient computational framework for computation. Two routes appear now possible for furthering the CGR transformation. One possibility if to relinquish the requirement to project sequences in two dimensions and instead use unitary hypercube maps . The other possibility assisted by the formulation reported here is to revisit CGR Kernel functions  to enable feature detection through entropic profiling  for non-genomics sequences. Interestingly, the 4 unit alphabet of genomic code is the exact point where both routes are equivalent, an observation that has been noted, and used, by other authors analyzing non-genomic sequences .
A generic formulation of Chaos Game Representation dividing ratios for 2-dimensional displaying was identified. The new formulation determined ratios for the iterative function that produce the reference representations for alphabets length 2 to 4. Specifically, it produces Sierpinsky's triangle for n = 3 and the original, homogeneously covered Chaos Game representation for n = 4. Given the fact that the density of CGR correspond to a order-free transition matrix – each consecutively nested polygon corresponds to an additional Markov chain order – the value of consistent graphical representation techniques is, potentially, enormous. Furthermore, the growing interest on alignment-free sequence dissimilarity metrics suggests a new role for Chaos Game iterative functions as a scale-free approach to word-statistics. The formulation identified to optimize the dividing ratio is, as is often the case with fractal processes, an empirical result that should be object of further analysis to uncover simpler and more meaningful patterns. Nevertheless, and regardless of the actual CGR computation being performed in two or more dimensions for sequence with longer alphabets, such as proteins, a generic graphical visualization technique is now at hand.
The m-code (Matlab, Mathworks Inc) used to generate the figures as well as the application depicted in Figure one are made available with open source at http://genechaos.org.
Jeffrey HJ: Chaos game representation of gene structure. Nucleic Acids Res. 1990, 18 (8): 2163-2170. 10.1093/nar/18.8.2163.
Oliver JL, Bernaola-Galvan P, Guerrero-Garcia J, Roman-Roldan R: Entropic profiles of DNA sequences through chaos-game-derived images. J Theor Biol. 1993, 160 (4): 457-470. 10.1006/jtbi.1993.1030.
Pleissner KP, Wernisch L, Oswald H, Fleck E: Representation of amino acid sequences as two-dimensional point patterns. Electrophoresis. 1997, 18 (15): 2709-2713. 10.1002/elps.1150181504.
Hill KA, Schisler NJ, Singh SM: Chaos game representation of coding regions of human globin genes and alcohol dehydrogenase genes of phylogenetically divergent species. J Mol Evol. 1992, 35 (3): 261-269. 10.1007/BF00178602.
Dutta C, Das J: Mathematical characterization of Chaos Game Representation. New algorithms for nucleotide sequence analysis. J Mol Biol. 1992, 228 (3): 715-719. 10.1016/0022-2836(92)90857-G.
Deschavanne PJ, Giron A, Vilain J, Fagot G, Fertil B: Genomic signature: characterization and classification of species assessed by chaos game representation of sequences. Mol Biol Evol. 1999, 16 (10): 1391-1399.
Choi JH, Cho HG: Analysis of common k-mers for whole genome sequences using SSB-tree. Genome Inform. 2002, 13: 30-41.
Basu S, Pan A, Dutta C, Das J: Chaos game representation of proteins. J Mol Graph Model. 1997, 15 (5): 279-289. 10.1016/S1093-3263(97)00106-X.
Wu ZB: Metric representation of DNA sequences. Electrophoresis. 2000, 21 (12): 2321-2326. 10.1002/1522-2683(20000701)21:12<2321::AID-ELPS2321>3.0.CO;2-S.
Tino P: Spacial Representation of Symbolic Sequences Through Iterative Function Systems. IEEE Transationa on Systems, Man, and Cybernetics – Part A. 1999, 29 (4): 386-393. 10.1109/3468.769757.
Vinga S, Almeida JS: Local Renyi entropic profiles of DNA sequences. BMC Bioinformatics. 2007, 8 (1): 393-10.1186/1471-2105-8-393.
Vinga S, Almeida JS: Renyi continuous entropy of DNA sequences. J Theor Biol. 2004, 231 (3): 377-388. 10.1016/j.jtbi.2004.06.030.
Almeida JS, Carrico JA, Maretzek A, Noble PA, Fletcher M: Analysis of genomic sequences by Chaos Game Representation. Bioinformatics. 2001, 17 (5): 429-437. 10.1093/bioinformatics/17.5.429.
Almeida JS, Vinga S: Computing distribution of scale independent motifs in biological sequences. Algorithms Mol Biol. 2006, 1: 18-10.1186/1748-7188-1-18.
Almeida JS, Vinga S: Universal sequence map (USM) of arbitrary discrete sequences. BMC Bioinformatics. 2002, 3 (1): 6-10.1186/1471-2105-3-6.
Schwacke J, Almeida JS: Efficient Boolean implementation of universal sequence maps (bUSM). BMC Bioinformatics. 2002, 3 (1): 28-10.1186/1471-2105-3-28.
Vinga S, Almeida J: Alignment-free sequence comparison-a review. Bioinformatics. 2003, 19 (4): 513-523. 10.1093/bioinformatics/btg005.
Yu ZG, Anh V, Lau KS: Chaos game representation of protein sequences based on the detailed HP model and their multifractal and correlation analyses. J Theor Biol. 2004, 226 (3): 341-348. 10.1016/j.jtbi.2003.09.009.
Cenac P: Test on the structure of biological sequences via Chaos Game Representation. Stat Appl Genet Mol Biol. 2005, 4: Article27-
Fu W, Wang Y, Lu D: Multifractal Analysis of Genomic Sequences CGR Images. Conf Proc IEEE Eng Med Biol Soc. 2005, 5: 4783-4786.
Fu W, Wang Y, Lu D: Multifractal analysis of genomes sequences' CGR graph. Journal of Biomedical Engineering. 2007, 24 (3): 522-525.
Joseph J, Sasikumar R: Chaos game representation for comparison of whole genomes. BMC Bioinformatics. 2006, 7: 243-10.1186/1471-2105-7-243.
Vinga S, Gouveia-Oliveira R, Almeida JS: Comparative evaluation of word composition distances for the recognition of SCOP relationships. Bioinformatics. 2004, 20: 206-215. 10.1093/bioinformatics/btg392.
Andreeva A, Howorth D, Brenner SE, Hubbard TJ, Chothia C, Murzin AG: SCOP database in 2004: refinements integrate structure and sequence family data. Nucleic Acids Res. 2004, 32 (1): D226-229. 10.1093/nar/gkh039.
Almeida JS, Vinga S: Universal sequence map (USM) of arbitrary discrete sequences. BMC bioinformatics. 2002, 3: 6-10.1186/1471-2105-3-6.
Deschavanne P, Tuffery P: Exploring an alignment free approach for protein classification and structural class prediction. Biochimie. 2008, 90 (4): 615-625. 10.1016/j.biochi.2007.11.004.
Fiser A, Tusnady GE, Simon I: Chaos game representation of protein structures. J Mol Graph. 1994, 12 (4): 302-304. 10.1016/0263-7855(94)80109-6.
Jeffrey HJ: Chaos Game Visualization of Sequences. Computation and Graphics. 1992, 16 (1): 25-33. 10.1016/0097-8493(92)90067-6.
Giles J: Benoit Mandelbrot: father of fractals. Nature. 2004, 432 (7015): 266-267. 10.1038/432266a.
Gallos LK, Song C, Havlin S, Makse HA: Scaling theory of transport in complex biological networks. Proc Natl Acad Sci USA. 2007, 104 (19): 7746-7751. 10.1073/pnas.0700250104.
This work was supported by the NHLBI Proteomics Initiative through contract N01-HV-28181.
JSA identified the optimal packing ration (Equation 2) and wrote the report in close consultation with SV. JSA and SV have jointly performed the exploratory work leading to the three point conceptualization of optimal increment ratios in iterative maps, described at the end of the Background section.
About this article
Cite this article
Almeida, J.S., Vinga, S. Biological sequences as pictures – a generic two dimensional solution for iterated maps. BMC Bioinformatics 10, 100 (2009). https://doi.org/10.1186/1471-2105-10-100
- Biological Sequence
- Symbolic Sequence
- Circle Packing
- Aminoacid Sequence
- Markov Chain Order