# On the Existence of Wavelet Symmetries in Archaea DNA

Authors: Carlo Cattani
Journal: Computational and Mathematical Methods in Medicine
 (2012)

---

## Abstract

This paper deals with the complex unit roots representation
of archea DNA sequences and the analysis of symmetries in
the wavelet coefficients of the digitalized sequence. It is shown that
even for extremophile archaea, the distribution of nucleotides
has to fulfill some (mathematical) constraints in such a way that the
wavelet coefficients are symmetrically distributed, with respect to the
nucleotides distribution.


---

## Body

## 1. Introduction

In some recent papers the existence of symmetries in nucleotide distribution has been studied for several living organisms [1–6] including mammals, fungi [1–4], and viruses [5, 6]. Thus showing that any (investigated) DNA sequence, when converted into a digital sequence, features some fractal shape of its DNA walk and an apparently random-like distribution. However, when the short wavelet transform maps the digital sequence into the space of wavelet coefficients, and these coefficients are clustered then they are located along some symmetrical shapes.One of the main tasks of this paper is to show that although the distribution of nucleotide, in any DNA sequence, can be considered as randomly given, when we compare a random sequence (and the corresponding random walk) with a DNA sequence (and walk) it can be seen that there exists some distinctions. So that the nucleotides distribution seems to side with a random distribution with some constraints. These constraints (rules) are singled out in the following, by showing the existence of hidden geometry which underlies the structure of a DNA sequence.In other words, nucleotides are distributed along any DNA sequence at first apparently randomly but at second analysis according to some (statistical) mathematical constraints which does not allow a given nucleotide to be arbitrarily followed by any other remaining nucleotides.It is interesting to notice that even in the primitives organisms which billions of years ago have been colonizing the earth under extreme conditions of life, their DNA has to fulfill the same constraints of the more evolved DNAs.In order to achieve this goal some fundamental steps have to be taken into consideration and discussed.(1)
Since DNA is a sequence of symbols, a map of these symbols into numbers has to be defined. In the following we will consider the complex unit roots map, which has the advantage of being unitary and distributed along the unit circle.(2)
The indicator matrix is defined on the the indicator map. This matrix is important in order to draw the dot plot of the DNA sequence and from this plot we can see that apparently nucleotides seem to be randomly distributed. However, we will show by wavelet analysis that they look randomly distributed, while they are not.(3)
The Ulam spiral adapted to DNA sequences is defined in order to single out some geometrical patterns.(4)
Random walks on DNA, or short DNA walks, show that the random walks look like fractals.(5)
The analysis of clusters of wavelet coefficients show that DNA walks have to fulfill some geometrical constraints.In all DNA sequences, analyzed so far, for different kinds of living organisms, this geometrical symmetry [1–6] has been detected. In the following this analysis is extended also to archaea, since they might be considered at the early stage of life and their DNA is compared with more evolved microorganisms as bacteria.It will be shown that, inspite of the many similarities with random sequences, only the wavelet analysis makes it possible to single out some distinctions. In particular, the wavelet coefficients of all (analyzed) organisms tend to fulfill a minimum principle for the energy of the signal. Also the archaea which often live in extreme environments have to fulfill the same geometrical rule of any other living organism.The analysis of DNA by wavelets [7–9], as seen in [8–12], helps to single out local behavior and singularities [7, 13] or to express the scale invariance of coefficients [14]. Also multifractal nature of the time series [15–17] can be easily detected by wavelet analysis.Some previous paper have studied various sequences of DNA such as leukemia tet variants, influenza viruses such as the A (H1N1) variant, mammalian, and a fungus (see [1–3, 14]) provided by the National Center for Biotechnology Information [18–21]. In all these papers it was observed that DNA has to fulfill not only some chemical steady state given by the chemical ligands but also some symmetrical distribution of nucleotide along the sequence. In other words, base pairs have to be placed exactly in some positions.According to previous results, it will be shown that as any other living organisms also these elementary organisms have DNA walks with fractal shape and wavelet coefficients bounded on a short-range wavelet transform. In other words, also anaerobic organism which should be understood as the most elementary at the first step of life have the same symmetries on wavelet coefficients as for more evolved organism, so that life has to fulfill some constrained distribution of nucleotides in order to give rise to some organism even at the most elementary step.In particular, in Section2, some remarks about the analysed data are given. Section 3 deals with some elementary plots which can easily visualize the distribution of nucleotides. The Ulam spiral plot is also proposed for the first time and it is observed a different distribution of weak/strong hydrogen bonds. Section 4 provides some definitions about parameters of complexity. We will notice that all these parameters give rise to the same classification of organism. Section 4 proposes a complex numerical representation of DNA chains and random walks, while in final Section 6 the short wavelet trasform is given in order to single out some symmetries at the lower order of transform.


## 2. Materials and Methods

In the following we will take into consideration some genome, complete sequences of DNA, concerning the following archaea:h1:
Aeropyrum pernix K1, complete genome. DNA, circular, 1669696 bp, [18–21], accession BA000002.3. Lineage: Archaea; Crenarchaeota; Thermoprotei; Desulfurococcales; Desulfurococcaceae; Aeropyrum; Aeropyrum pernix; Aeropyrum pernix K1. This organism, which was the first strictly aerobic hyperthermophilic archaeon sequenced, was isolated from sulfuric gases in Kodakara-Jima Island, Japan in 1993.h2:
Acidianus hospitalis W1, complete genome. DNA, circular, 2137654 bp, [18–21], accession CP002535. Lineage: Archaea; Crenarchaeota; Thermoprotei; Sulfolobales; Sulfolobaceae; Acidianus; Acidianus hospitalis; Acidianus hospitalis W1h3:
Acidilobus saccharovorans 345-15. complete genome. DNA, circular, 2137654 bp, [18–21], accession CP001742.1. Lineage: Archaea; Crenarchaeota; Thermoprotei; Acidilobales; Acidilobaceae; Acidilobus; Acidilobus saccharovorans; Acidilobus saccharovorans 345-15. Anaerobic bacteria found in hot springs.to be compared with the following (aerobic/anaerobic) bacteria/fungi:b1:
Mycoplasma putrefaciens KS1 chromosome, complete genome. DNA, circular, length 832603 bp, [18–21], accession NC 015946,. Lineage: Bacteria; Tenericutes; Mollicutes; Mycoplasmatales; Mycoplasmataceae; Mycoplasma; Mycoplasma putrefaciens; Mycoplasma putrefaciens KS1.b2:
Mortierella verticillata mitochondrion, complete genome. dsDNA, circular, length 58745 bp, [18–21], accession NC 006838. Lineage: Eukaryota; Opisthokonta; Fungi; Fungi incertae sedis; Basal fungal lineages; Mucoromycotina; Mortierellales; Mortierellaceae; Mortierella; Mortierella verticillata.b3:
Blattabacterium sp. (Periplaneta Americana) str. BPLAN, complete genome. DNA, circular, length 636994 nt, [18–21], accession NC 013418. Lineage: Bacteria; Bacteroidetes/Chlorobi group; Bacteroidetes; Flavobacteria; Flavobacteriales; Blattabacteriaceae; Blattabacterium; Blattabacterium sp. (Periplaneta Americana); Blattabacterium sp. (Periplaneta Americana) str. BPLAN.Moreover we will compare DNA sequences with artificial sequences of nucleotides randomly taken (see Section4).


### 2.1. Archaea

Archaea are a group of elementary single-cell microorganisms, having no cell nucleus or any other membrane-bound organelles within their cells. They are similar to bacteria, since they have the same size and shape (apart few exceptions) and the generally similar cell structure. However, the evolutionary history of archaea and their biochemistry has significant differences with regard to other forms of life. Therefore they are considered as members of a phylogenetic group distinct from bacteria and eukaryota.Archaea during their evolution have been spreading all over the Earth in almost all habitats [22, 23] existing in a broad range of habitats, being one of the major contribution (20%) to earth’s biomass. The most peculiar feature of archaea is that they can live in some environments with extreme life conditions (thus being considered as extremophiles [22, 24]). Indeed, some archaea survive to high temperatures, over 100°C, while others can live in very cold habitats or highly saline, acidic, or alkaline water. Nevertheless some archaea are living in mild conditions.It has been also recognized that the archaea may be the most ancient organisms on the Earth, so that archaea, and eukaryotes are probably diverged early from an ancestral colony of organisms.We will see, in the following, that archaea DNA it looks very close to random sequences so that we can assume that the ancestral organism were evolving by random permutations from a primitive assembly of nucleotides.  So that the evolution can be seen as a tendency to a steady state far from the randomness. Therefore, the bacteria’s DNA (and other eukaryotes’ [1–6]), as a result of the evolution, shows the existence of some hidden stability.


## 2.1. Archaea

Archaea are a group of elementary single-cell microorganisms, having no cell nucleus or any other membrane-bound organelles within their cells. They are similar to bacteria, since they have the same size and shape (apart few exceptions) and the generally similar cell structure. However, the evolutionary history of archaea and their biochemistry has significant differences with regard to other forms of life. Therefore they are considered as members of a phylogenetic group distinct from bacteria and eukaryota.Archaea during their evolution have been spreading all over the Earth in almost all habitats [22, 23] existing in a broad range of habitats, being one of the major contribution (20%) to earth’s biomass. The most peculiar feature of archaea is that they can live in some environments with extreme life conditions (thus being considered as extremophiles [22, 24]). Indeed, some archaea survive to high temperatures, over 100°C, while others can live in very cold habitats or highly saline, acidic, or alkaline water. Nevertheless some archaea are living in mild conditions.It has been also recognized that the archaea may be the most ancient organisms on the Earth, so that archaea, and eukaryotes are probably diverged early from an ancestral colony of organisms.We will see, in the following, that archaea DNA it looks very close to random sequences so that we can assume that the ancestral organism were evolving by random permutations from a primitive assembly of nucleotides.  So that the evolution can be seen as a tendency to a steady state far from the randomness. Therefore, the bacteria’s DNA (and other eukaryotes’ [1–6]), as a result of the evolution, shows the existence of some hidden stability.


## 3. Correlation Plots

In this section we will consider some elementary plots from where it is possible to visualize autocorrelation, distribution law of nucleotides and to measure some fundamental parameters by using frequency count.Let(1)A=def{A,C,G,T}
be the finite set (alphabet) of nucleotides (nucleic acids): adenine (A), cytosine (C), guanine (G), thymine (T), and x∈𝒜 any member of the alphabet. Nucleic acids are further grouped according to their ligand properties as(a)
purine{A,G}, pyrimidine {C,T},(b)
amino{A,C}, keto {G,T},(c)
weak hydrogen bonds{A,T}, strong hydrogen bond {G,C}.A DNA sequence is the finite symbolic sequence(2)S=N×A
so that(3)S≝{xh}h=1,…,N,N<∞
with(4)xh≝(h,x)=x(h),(h=1,2,…,N;x∈A)
being the nucleotide x at the position h.In general we can define anℓ-length alphabet as follows: let the ℓ-length DNA word be defined by the ℓ-combination of the 4 nucleotides (1). For each fixed length ℓ there are 4ℓ words, however not all of them can be considered, from biological point of view, as independent instances (see, e.g., Table 1), for this we define the ℓ-length alphabet as the set of ℓ-length independent words:(5)Al≝{a1,a2,…,aMl},Ml≝|Al|≤4l
with |⋯| cardinality of the set and(6)l≝length(aj),(j=1,…,Ml).
For instance with ℓ=1, the alphabet is 𝒜1=𝒜={A,C,G,T}, with ℓ=3 the alphabet is given by the 20 amino acids(7)A3={M,E,Q,D,R,T,N,H,V,G,L,S,P,F,I,C,A,K,Y,W}
each amino acid being represented by a 3-length word of Table 1.Table 1
Correspondence codons to amino acids.
Amino acidCodon1MMethionineATG2EGlutamic acidGAA, GAG3QGlutamineCAA, CAG4DAspartic acidGAT, GAC5RArginineCGT, CGC, CGA, CGG, AGA, AGG6TThreonineACT, ACC, ACA, ACG7NAsparagineAAT, AAC8HHistidineCAT, CAC9VValineGTT, GTC, GTA, GTG10GGlycineGGT, GGC, GGA, GGG11LLeucineTTA, TTG, CTT, CTC, CTA, CTG12SSerineTCT, TCC, TCA, TCG, AGT, AGC13PProlineCCT,CCC,CCA,CCG14FPhenylalanineTTT,TTC15IIsoleucineATT,ATC,ATA16CCysteineTGT,TGC17AAlanineGCT,GCC,GCA,GCG18KLysineAAA,AAG19YThyroxineTAT,TAC20WTryptophanTGGStopTAA,TAG,TGALet𝒮N be an N-length ordered sequence of nucleotides {A,C,G,T} and 𝒜ℓ the chosen alphabet, a DNA sequence of words is the finite symbolic sequence(8)Dl(SN)=N×Al
so that(9)Dl(SN)=def{xh}h=1,…,N,(x∈Al;N<∞)
with(10)xh=def(h,x),(h=1,2,…,N;x∈Al)
being the word x at the position h.


### 3.1. Indicator Matrix

The 2D indicator function, based on the 1D definition given in [25], is the map(11)u:S×S⟶{0,1}
such that
(12)u(xh,xk)=def{1ifxh=xk,0ifxh≠xk,(xh∈S,xk∈S),
with
(13)u(xh,xk)=u(xk,xh),u(xh,xh)=1
and, where for short, we have assumed(14)S=defD1(SN).
According to (12), the indicator of an N-length sequence can be easily represented by the N×N sparse symmetric matrix of binary values {0,1} which results from the indicator matrix (see also [3–5])(15)uhk=defu(xh,xk),(xh∈S,xk∈S;h,k=1,…,N),
being, explicitly(16)⋮⋮⋮⋮⋮⋮⋮⋮⋮⋮⋮⋰G0100000001⋯C0001000010⋯A1000101100⋯A1000101100⋯T0010010000⋯A0000100100⋯C0001000010⋯T0010010000⋯G0100000001⋯A1000100100⋯uhkAGTCATAACG⋯This squared matrix can be plotted in 2 dimensions by putting a black dot whereuhk=1 and white spot when uhk=0 (Figure 1) thus giving rise to the two-dimensional dot plot, which is a special case of the recurrence plot [26].Indicator matrix for: (a1) pseudorandom 70-length sequence; (a2) pseudo-periodic 70-length sequence with period π=35; (b1) 70-length DNA sequence of Mycoplasma KS1 bacter; (h3) 70-length DNA sequence of Acidilobus Archaea.A simple generalization of this matrix can be considered for the alphabets𝒜ℓ, as follows. By choosing the 3 alphabet of amino acids, the 2D indicator function is the map(17)u:D3(SN)×D3(SN)⟶{0,1}
such that(18)u(xh,xk)=def{1ifxh=xk,0ifxh≠xk,(xh∈D3(SN),xk∈D3(SN)),
with(19)u(xh,xk)=u(xk,xh),u(xh,xh)=1.According to (12), the indicator, on the 3-alphabet of amino acids of an N-length sequence can be easily represented by the N×N sparse symmetric matrix of binary values {0,1}:(20)uhk=defu(xh,xk),(xh∈D3(SN),xk∈D3(SN);h,k=1,…,N),
being, explicitly(21)⋮⋮⋮⋮⋮⋮⋮⋮⋮⋮⋮⋰M1000000001⋯Q0000000010⋯R0011000100⋯T0000011000⋯T0000011000⋯E0000100100⋯R0011000000⋯R0011000100⋯K0100000000⋯M1000000001⋯uhkMKRRETTRQM⋯
With the graphical representation of this matrix we can also show the correlation of amino acids.


### 3.2. Test Sequences

In the following, in order to single out the main features of biological sequences, we will compare the DNA sequence with some test sequences.(1)
PseudorandomN-length sequence of nucleotides is the sequence {ℛi}i=1,…,Nℓ where ri is a symbol randomly chosen in the alphabet 𝒜ℓ, like for example, (ℓ=1):
(22){A,C,A,G,T,A,T,G,G,A,T,T,A,C,C,G,…}.(2)
PseudoperiodicN-sequence of nucleotides with period π is the direct sum of a given π-length pseudorandom sequence, such that N=kπ,(k∈ℕ) and ℛi=ℛi+π, for example,
(23){A,C,A,G,A,C,A,G,A,C,A,G,A,C,A,G,…},(π=4).
When π=1 we have a pseudorandom sequence.If we plot the indicator matrix of some bacteria and compare it with a pseudorandom and periodic sequence, we can see that (Figure1)(1)
the main diagonal is a symmetry axis for the plot;(2)
there are some motifs which are repeated at different scales like in a fractal;(3)
periodicity is detected by parallel lines to the main diagonal (Figure1(a2));(4)
empty spaces are more distributed than filled spaces, in the sense that the matrixuhk is a sparse matrix (having more 0’s than 1’s);(5)
it seems that there are some square-like islands where black spots are more concentrated; these islands show the persistence of a nucleotide (Figures1(a2) and 1(b1));(6)
the dot plot of archaea is very similar to the dot plot of a random sequence (Figures1(a1) and 1(h3)).It can be noticed that DNA sequences of a living organism resemble (Figure1) random sequences, with some short range influence, built on the same alphabet. This has been taken as an axiom of nucleotides distribution, so that DNA sequences are often considered as Markov chain [27]. However, there are some hidden rules in combining the nucleotides and these rules lead, during the evolution, to a steady distribution. In fact, the more primitive the sequence is, the more randomly distributed the nucleotides are. It seems that as a consequence of the evolution, nucleotides move from a disordered aggregation toward a more organized structure, shown by the growing islands in the dot plot. The biological evolution is such that the challenge for the self-organization might follow from random permutations of a primitive disordered sequence so that the organization, that is, the complexity, is only the result of many arbitrary permutations of randomness. During the challenge for complexity, DNA sequence becomes “less random” and it loses some kind of energy.From the graphical representation of the indicator matrix for bacteria and amino acids we can see a more sparse matrix, but with some typical plots (Figure2).Indicator matrix for the first 100 amino acids of (h1) Aeropyrum pernix K1, (h2) Acidianus hospitalis W1, (h3) Acidilobus saccharovorans 345-15 (b1) Mycoplasma putrefaciens KS1, (b2) Mortierella verticillata, and (b3) Blattabacterium sp.


### 3.3. Spiral Plot

In this section we consider a 2D distribution of nucleotides, following the idea given by Ulam for the distribution of primes, along an Ulam-like spiral [28]. In order to find some patterns in their distribution, nucleotides are arranged along a rectangular spiral. This is equivalent to mapping the 1D sequence of integers into a 2D sequence as follows:(24)X11{0,0}X22{1,0}X33{1,1}X44{0,1}X55{-1,1}X66{-1,0}X77{-1,-1}X88{0,-1}X99{1,-1}X1010{2,-1}X1111{2,0}⋮⋮⋮For instance the sequence(25){A,T,G,G,A,A,G,A,T,A,A,G,…}
distributed along the spiral looks like Figure 3.Figure 3
Distribution of nucleotides on a rectangular spiral.For each nucleotide we can draw a spiral containing the distribution of only one acid nucleic. To each organism there correspond four plots, forA,C,G,T, respectively.Let us first note that on a random sequence (Figure4) the four distribution are equivalent.Spiral distribution of the first 3752 nucleotides for the random sequence.By comparing the spirals of bacteria, random and archaea (Figures4, 5, 6, 7, 8, 9, 10) we can see that there is a different distribution of each nucleotide. However the more evolved organism tends to have a higher percentage of weak hydrogen bonds (Figures 5, 6 and 7), so that we can assume the following.Spiral distribution of the first 3752 nucleotides forMycoplasma putrefaciens KS1.Spiral distribution of the first 3752 nucleotides forMortierella verticillata.Spiral distribution of the first 3752 nucleotides forBlattabacterium sp.Spiral distribution of the first 3752 nucleotides forAeropyrum pernix K1..Spiral distribution of the first 3752 nucleotides forAcidianus hospitalis W1.Spiral distribution of the first 3752 nucleotides forAcidilobus saccharovorans 345-15.Conjecture 1.
During the evolution, the distribution of nucleotides changes in a such way that strong hydrogen bonds tend to become weak.It should be noticed that along these spirals, there is a one-to-one mapλ between ℕ and the points of the spiral (with integer coordinates) in ℜ2(26)λ:N⟼γ⊂R×R
so that(27)λ(n)=(a,b),(n∈N;(a,b)∈γ⊂R×R;a∈Z,b∈Z),λ-1(a,b)=n.
This bijective map can be considered also between ℕ and the complex space ℂ so that each natural number corresponds to a complex number (with integer coefficients)(28)λ(n)=z=defa+ib,(n∈N;a,b∈Z;z∈C).Since these spirals seem to fill in a finite region of the plane we can evaluate the complexity of each curve by typical fractal measures.


## 3.1. Indicator Matrix

The 2D indicator function, based on the 1D definition given in [25], is the map(11)u:S×S⟶{0,1}
such that
(12)u(xh,xk)=def{1ifxh=xk,0ifxh≠xk,(xh∈S,xk∈S),
with
(13)u(xh,xk)=u(xk,xh),u(xh,xh)=1
and, where for short, we have assumed(14)S=defD1(SN).
According to (12), the indicator of an N-length sequence can be easily represented by the N×N sparse symmetric matrix of binary values {0,1} which results from the indicator matrix (see also [3–5])(15)uhk=defu(xh,xk),(xh∈S,xk∈S;h,k=1,…,N),
being, explicitly(16)⋮⋮⋮⋮⋮⋮⋮⋮⋮⋮⋮⋰G0100000001⋯C0001000010⋯A1000101100⋯A1000101100⋯T0010010000⋯A0000100100⋯C0001000010⋯T0010010000⋯G0100000001⋯A1000100100⋯uhkAGTCATAACG⋯This squared matrix can be plotted in 2 dimensions by putting a black dot whereuhk=1 and white spot when uhk=0 (Figure 1) thus giving rise to the two-dimensional dot plot, which is a special case of the recurrence plot [26].Indicator matrix for: (a1) pseudorandom 70-length sequence; (a2) pseudo-periodic 70-length sequence with period π=35; (b1) 70-length DNA sequence of Mycoplasma KS1 bacter; (h3) 70-length DNA sequence of Acidilobus Archaea.A simple generalization of this matrix can be considered for the alphabets𝒜ℓ, as follows. By choosing the 3 alphabet of amino acids, the 2D indicator function is the map(17)u:D3(SN)×D3(SN)⟶{0,1}
such that(18)u(xh,xk)=def{1ifxh=xk,0ifxh≠xk,(xh∈D3(SN),xk∈D3(SN)),
with(19)u(xh,xk)=u(xk,xh),u(xh,xh)=1.According to (12), the indicator, on the 3-alphabet of amino acids of an N-length sequence can be easily represented by the N×N sparse symmetric matrix of binary values {0,1}:(20)uhk=defu(xh,xk),(xh∈D3(SN),xk∈D3(SN);h,k=1,…,N),
being, explicitly(21)⋮⋮⋮⋮⋮⋮⋮⋮⋮⋮⋮⋰M1000000001⋯Q0000000010⋯R0011000100⋯T0000011000⋯T0000011000⋯E0000100100⋯R0011000000⋯R0011000100⋯K0100000000⋯M1000000001⋯uhkMKRRETTRQM⋯
With the graphical representation of this matrix we can also show the correlation of amino acids.


## 3.2. Test Sequences

In the following, in order to single out the main features of biological sequences, we will compare the DNA sequence with some test sequences.(1)
PseudorandomN-length sequence of nucleotides is the sequence {ℛi}i=1,…,Nℓ where ri is a symbol randomly chosen in the alphabet 𝒜ℓ, like for example, (ℓ=1):
(22){A,C,A,G,T,A,T,G,G,A,T,T,A,C,C,G,…}.(2)
PseudoperiodicN-sequence of nucleotides with period π is the direct sum of a given π-length pseudorandom sequence, such that N=kπ,(k∈ℕ) and ℛi=ℛi+π, for example,
(23){A,C,A,G,A,C,A,G,A,C,A,G,A,C,A,G,…},(π=4).
When π=1 we have a pseudorandom sequence.If we plot the indicator matrix of some bacteria and compare it with a pseudorandom and periodic sequence, we can see that (Figure1)(1)
the main diagonal is a symmetry axis for the plot;(2)
there are some motifs which are repeated at different scales like in a fractal;(3)
periodicity is detected by parallel lines to the main diagonal (Figure1(a2));(4)
empty spaces are more distributed than filled spaces, in the sense that the matrixuhk is a sparse matrix (having more 0’s than 1’s);(5)
it seems that there are some square-like islands where black spots are more concentrated; these islands show the persistence of a nucleotide (Figures1(a2) and 1(b1));(6)
the dot plot of archaea is very similar to the dot plot of a random sequence (Figures1(a1) and 1(h3)).It can be noticed that DNA sequences of a living organism resemble (Figure1) random sequences, with some short range influence, built on the same alphabet. This has been taken as an axiom of nucleotides distribution, so that DNA sequences are often considered as Markov chain [27]. However, there are some hidden rules in combining the nucleotides and these rules lead, during the evolution, to a steady distribution. In fact, the more primitive the sequence is, the more randomly distributed the nucleotides are. It seems that as a consequence of the evolution, nucleotides move from a disordered aggregation toward a more organized structure, shown by the growing islands in the dot plot. The biological evolution is such that the challenge for the self-organization might follow from random permutations of a primitive disordered sequence so that the organization, that is, the complexity, is only the result of many arbitrary permutations of randomness. During the challenge for complexity, DNA sequence becomes “less random” and it loses some kind of energy.From the graphical representation of the indicator matrix for bacteria and amino acids we can see a more sparse matrix, but with some typical plots (Figure2).Indicator matrix for the first 100 amino acids of (h1) Aeropyrum pernix K1, (h2) Acidianus hospitalis W1, (h3) Acidilobus saccharovorans 345-15 (b1) Mycoplasma putrefaciens KS1, (b2) Mortierella verticillata, and (b3) Blattabacterium sp.


## 3.3. Spiral Plot

In this section we consider a 2D distribution of nucleotides, following the idea given by Ulam for the distribution of primes, along an Ulam-like spiral [28]. In order to find some patterns in their distribution, nucleotides are arranged along a rectangular spiral. This is equivalent to mapping the 1D sequence of integers into a 2D sequence as follows:(24)X11{0,0}X22{1,0}X33{1,1}X44{0,1}X55{-1,1}X66{-1,0}X77{-1,-1}X88{0,-1}X99{1,-1}X1010{2,-1}X1111{2,0}⋮⋮⋮For instance the sequence(25){A,T,G,G,A,A,G,A,T,A,A,G,…}
distributed along the spiral looks like Figure 3.Figure 3
Distribution of nucleotides on a rectangular spiral.For each nucleotide we can draw a spiral containing the distribution of only one acid nucleic. To each organism there correspond four plots, forA,C,G,T, respectively.Let us first note that on a random sequence (Figure4) the four distribution are equivalent.Spiral distribution of the first 3752 nucleotides for the random sequence.By comparing the spirals of bacteria, random and archaea (Figures4, 5, 6, 7, 8, 9, 10) we can see that there is a different distribution of each nucleotide. However the more evolved organism tends to have a higher percentage of weak hydrogen bonds (Figures 5, 6 and 7), so that we can assume the following.Spiral distribution of the first 3752 nucleotides forMycoplasma putrefaciens KS1.Spiral distribution of the first 3752 nucleotides forMortierella verticillata.Spiral distribution of the first 3752 nucleotides forBlattabacterium sp.Spiral distribution of the first 3752 nucleotides forAeropyrum pernix K1..Spiral distribution of the first 3752 nucleotides forAcidianus hospitalis W1.Spiral distribution of the first 3752 nucleotides forAcidilobus saccharovorans 345-15.Conjecture 1.
During the evolution, the distribution of nucleotides changes in a such way that strong hydrogen bonds tend to become weak.It should be noticed that along these spirals, there is a one-to-one mapλ between ℕ and the points of the spiral (with integer coordinates) in ℜ2(26)λ:N⟼γ⊂R×R
so that(27)λ(n)=(a,b),(n∈N;(a,b)∈γ⊂R×R;a∈Z,b∈Z),λ-1(a,b)=n.
This bijective map can be considered also between ℕ and the complex space ℂ so that each natural number corresponds to a complex number (with integer coefficients)(28)λ(n)=z=defa+ib,(n∈N;a,b∈Z;z∈C).Since these spirals seem to fill in a finite region of the plane we can evaluate the complexity of each curve by typical fractal measures.


## 4. Parameters of Complexity

In this section we define some parameters, based on frequency distribution, which can measure the complexity of a DNA by computing the complexity of its representation in the complex plane (for a more detailed analysis see [29] and references therein).Let𝒮N be an N-length-ordered sequence of nucleotides, and(29)px(h),x∈A1={A,C,G,T}
be the probability to find the nucleotide x at the position h,1≤h≤N. According to (12) we define(30)ah=def∑j=1huAj,ch=def∑j=1huCj,gh=def∑j=1huGj,th=def∑j=1huTj,(1≤h≤N)
as the number of nucleotides in the h-length segment of 𝒮N, so that(31)ah+ch+gh+th=h.
The corresponding frequencies are(32)vx(h)=def1h∑j=1huxj,x∈A1,(1≤h≤N),
so that(33)vA(h)=ahh,vC(h)=chh,vG(h)=ghh,vT(h)=thh.We can assume that for large sequences(34)px(h)≅vx(h).


### 4.1. Randomness

Since for a random sequence the frequencies of nucleotides coincide for largen,(35)vA(n)≅vC(n)≅vG(n)≅vT(n)
we can define as randomness index the following:(36)R=def1-σ(vA(n),vC(n),vG(n),vT(n))
with σ being the variance, so that ℛ=1 for random sequence and ℛ=0 for a nonrandom sequence. Over the first 10000 nucleotides we have the randomness value of Table 2.Table 2
Randomness.
Mycoplasma putrefaciens0.696Mortierella verticillata0.779Blattabacterium0.743Aeropyrum pernix0.982Acidianus hospitalis0.828Acidilobus saccharouorans0.934pseudorandom0.999However, if we compute the randomness index over the frequencies of amino acids in the𝒜3 alphabet then we can observe a different distribution of values. Over the first 30000 nucleotides corresponding to 10000 amino acids, we have the randomness value of Table 3.Table 3
Randomness of amino acids distribution.
Mycoplasma putrefaciens0.946Mortierella verticillata0.938Blattabacterium0.953Aeropyrum pernix0.962Acidianus hospitalis0.916Acidilobus saccharouorans0.950pseudorandom0.963So that we can comment that the arising complexity of the words and alphabets shows a different randomness in each alphabet.


### 4.2. Complexity

As a simple measure of complexity [30–32], for an n-length sequence, the following has been proposed [33]: (37)K=1nlog⁡n!an!cn!gn!tn!.In Table4 the complexity of the first 100-length segment of the DNA sequences is computed. It is interesting to notice the more similarities between the archaea Acidilobus with the pseudorandom sequence than with the pseudoperiodic. Nucleotide distribution in primitive biosequences is more likely random than pseudodeterministic. Moreover, the evolution reduces the complexity of the sequence.Table 4
Complexity.
Mycoplasma putrefaciens1.151Mortierella verticillata1.285Blattabacterium1.197Aeropyrum pernix1.212Acidianus hospitalis1.231Acidilobus saccharouorans1.296Pseudorandom1.295


### 4.3. Fractal Dimension

The fractal dimension is computed on the dot plot, by the box counting algorithm [34, 35], as the average of the number p(n) of 1’s in the randomly taken n×n minors of the N×N indicator matrix uhk or equivalently the number p(n) of black dots in the randomly taken n×n squares over the dot plot(38)D=12N∑n=2Nlog⁡⁡p(n)log⁡⁡n.The explicit computation enables us to compare the fractal dimension on the first 100-length segments of DNA chains, with an approximation up to10-3 (see Table 5).Table 5
Fractal dimensions.
Mycoplasma putrefaciens1.283Mortierella verticillata1.296Blattabacterium1.287Aeropyrum pernix1.288Acidianus hospitalis1.290Acidilobus saccharouorans1.297pseudorandom1.298pseudoperiodic1.285If we compare the fractal dimensions of the bacteria with pseudorandom and pseudoperiodic we can see that the fractal dimension of nucleotide distribution ranges, for all variants, in the interval[1.28–1.30]. As expected, the more “random” sequences have higher fractal dimension.


### 4.4. Entropy

Another fundamental parameter, related to the information content of a sequence which measures the heterogeneity of data, is the information entropy (or Shannon entropy) [36–42]. Based on the axiom that less information implies a larger uncertainty and vice versa that more information leads us to a more deterministic model, the entropy concept has been recently offering some interesting interpretations about uncertainty in DNA. In fact, DNA as any other signal has been considered as a sequence of symbols carrying chemical-functional information.The normalized Shannon entropy [39, 40, 42] is defined, over the alphabet 𝒜ℓ, as(39)H(n)=-1log⁡l∑x∈Alpx(n)×{log⁡⁡px(n)ifpx(n)≠0,0ifpx(n)=0,
where px(n) should be computed for large sequences. According to (32), (34), we will approximate its value with(40)px(n)≅1n∑i=1nuxi,(x∈Al,1≤n≤N).However, the entropy is a parameter very similar to the complexity. In fact, it can be easily seen that (for the proof see [29]) the entropy H and the measure of complexity K differ for a factor. There follows that the entropy does not give any new information comparing with the previous parameters. As expected also the table of entropies classifies bacteria and archaea in the same way (Table 6).Table 6
Shannon entropy.
Mycoplasma putrefaciens0.877Mortierella verticillata0.976Blattabacterium0.911Aeropyrum pernix0.922Acidianus hospitalis0.937Acidilobus saccharovorans0.984pseudorandom0.984


## 4.1. Randomness

Since for a random sequence the frequencies of nucleotides coincide for largen,(35)vA(n)≅vC(n)≅vG(n)≅vT(n)
we can define as randomness index the following:(36)R=def1-σ(vA(n),vC(n),vG(n),vT(n))
with σ being the variance, so that ℛ=1 for random sequence and ℛ=0 for a nonrandom sequence. Over the first 10000 nucleotides we have the randomness value of Table 2.Table 2
Randomness.
Mycoplasma putrefaciens0.696Mortierella verticillata0.779Blattabacterium0.743Aeropyrum pernix0.982Acidianus hospitalis0.828Acidilobus saccharouorans0.934pseudorandom0.999However, if we compute the randomness index over the frequencies of amino acids in the𝒜3 alphabet then we can observe a different distribution of values. Over the first 30000 nucleotides corresponding to 10000 amino acids, we have the randomness value of Table 3.Table 3
Randomness of amino acids distribution.
Mycoplasma putrefaciens0.946Mortierella verticillata0.938Blattabacterium0.953Aeropyrum pernix0.962Acidianus hospitalis0.916Acidilobus saccharouorans0.950pseudorandom0.963So that we can comment that the arising complexity of the words and alphabets shows a different randomness in each alphabet.


## 4.2. Complexity

As a simple measure of complexity [30–32], for an n-length sequence, the following has been proposed [33]: (37)K=1nlog⁡n!an!cn!gn!tn!.In Table4 the complexity of the first 100-length segment of the DNA sequences is computed. It is interesting to notice the more similarities between the archaea Acidilobus with the pseudorandom sequence than with the pseudoperiodic. Nucleotide distribution in primitive biosequences is more likely random than pseudodeterministic. Moreover, the evolution reduces the complexity of the sequence.Table 4
Complexity.
Mycoplasma putrefaciens1.151Mortierella verticillata1.285Blattabacterium1.197Aeropyrum pernix1.212Acidianus hospitalis1.231Acidilobus saccharouorans1.296Pseudorandom1.295


## 4.3. Fractal Dimension

The fractal dimension is computed on the dot plot, by the box counting algorithm [34, 35], as the average of the number p(n) of 1’s in the randomly taken n×n minors of the N×N indicator matrix uhk or equivalently the number p(n) of black dots in the randomly taken n×n squares over the dot plot(38)D=12N∑n=2Nlog⁡⁡p(n)log⁡⁡n.The explicit computation enables us to compare the fractal dimension on the first 100-length segments of DNA chains, with an approximation up to10-3 (see Table 5).Table 5
Fractal dimensions.
Mycoplasma putrefaciens1.283Mortierella verticillata1.296Blattabacterium1.287Aeropyrum pernix1.288Acidianus hospitalis1.290Acidilobus saccharouorans1.297pseudorandom1.298pseudoperiodic1.285If we compare the fractal dimensions of the bacteria with pseudorandom and pseudoperiodic we can see that the fractal dimension of nucleotide distribution ranges, for all variants, in the interval[1.28–1.30]. As expected, the more “random” sequences have higher fractal dimension.


## 4.4. Entropy

Another fundamental parameter, related to the information content of a sequence which measures the heterogeneity of data, is the information entropy (or Shannon entropy) [36–42]. Based on the axiom that less information implies a larger uncertainty and vice versa that more information leads us to a more deterministic model, the entropy concept has been recently offering some interesting interpretations about uncertainty in DNA. In fact, DNA as any other signal has been considered as a sequence of symbols carrying chemical-functional information.The normalized Shannon entropy [39, 40, 42] is defined, over the alphabet 𝒜ℓ, as(39)H(n)=-1log⁡l∑x∈Alpx(n)×{log⁡⁡px(n)ifpx(n)≠0,0ifpx(n)=0,
where px(n) should be computed for large sequences. According to (32), (34), we will approximate its value with(40)px(n)≅1n∑i=1nuxi,(x∈Al,1≤n≤N).However, the entropy is a parameter very similar to the complexity. In fact, it can be easily seen that (for the proof see [29]) the entropy H and the measure of complexity K differ for a factor. There follows that the entropy does not give any new information comparing with the previous parameters. As expected also the table of entropies classifies bacteria and archaea in the same way (Table 6).Table 6
Shannon entropy.
Mycoplasma putrefaciens0.877Mortierella verticillata0.976Blattabacterium0.911Aeropyrum pernix0.922Acidianus hospitalis0.937Acidilobus saccharovorans0.984pseudorandom0.984


## 5. Complex Root Representation of DNA Words

The complex (digital) representation of a DNA sequence of words is the map of the symbolic sequence of words into a set of complex numbers and it is defined as(41)Dl(SN)⟶ρC
such that for each xh∈𝒟ℓ(SN) it is ρ(xh)∈ℂ.The complex root representation of the sequenceSN is the sequence 𝒟ℓ(SN) of complex numbers {yh}h=1,…,N defined as(42)yh=ρ(xh)=defe2πi(j-1)/|Al|,(j=1,…,|Al|,h=1,…,N)
with i=-1 being the imaginary unit. There follows that, independently on the alphabet, it is(43)|yh|=|e2πi(j-1)/|Al||=1,(∀l;h=1,…,N)
being all complex roots, of the unit, located on the unit circle of the complex plane ℂ1.For instance, with𝒜1={A,C,G,T}, the cardinality of the alphabet is |𝒜1|=4 and(44)ρ(A)=e0/4=1,j=1,ρ(C)=eπi/2=i,j=2,ρ(G)=eπi=-1,j=3,ρ(T)=eπi3/2=-i,j=4.Analogously, with𝒜3={M,E,…,W} it is |𝒜3|=20 and the 20 complex roots of unit(45)ρ(xn)=e2πi(n-1)/20,(n=1,…,20;xn∈A3)
so that explicitly is(46)ρ(M)=e2πi0/20=1,j=1,ρ(E)=eπi/10=14[2(5+5)+i(5-1)],j=2,ρ(Q)=eπi/5=14[1+5+i2(5-5)],j=3,⋮⋮ρ(W)=eπi19/10=14[2(5+5)-i(5-1)],j=20.
Therefore the complex representation of a DNA sequence is a sequence of complex numbers(47)yh=ξh+ηhi,ξh=R(yh),ηh=I(yh)
with yh given by (42).Ann-length pseudorandom (white noise) complex sequence belonging to the unit circle can be defined directly by using some random exponents(48)Rn=def(-1)rnisn,|Rn|=1,
with rn,sn being random values in the set {0,ℕ}.


### 5.1. Random Walks

Random walk on the complex sequenceYN is defined as the series ZN={zn}n=1,…,N(49)zn=def∑k=1,…,nyk,n=1,…,N
which is the cumulative sum(50){y1,y1+y2,…,∑s=1nys…,∑s=1Nys}.
When yk=ρ(xk) with xk∈𝒜ℓ and Xk∈SN we will properly call these walks as DNA walk. When the yk are randomly generated we will call them random walks.By remembering the definition of frequencies, DNA walk is the complex value signal{Zn}n=0,…,N-1 with(51)zn=(R[zn],I[zn])=(an-gn)+(tn-cn)i,zn∈C1,
where the coefficients an,gn,tn,cn given by (12) fulfill the condition (31).If we compare the DNA walks (Figure11) some primitive archaea such as h3  are very similar to a random walk (Figure 13). In particular archaea seem to grow less than other bacteria (with the exception of b2).Walks on the first 200 nucleotides: (b1) Mycoplasma putrefaciens, (b2) Mortierella verticillata, (b3) Blattabacterium, (h1) Aeropyrum pernix, (h2) Acidianus hospitalis, and (h3) Acidilobus saccharovorans.It is interesting also to notice that the random walks on amino acids (Figure12) show that more evolved organisms have some “periodic” behavior, while the absolute value of walks on archaea is growing fast.Absolute value of walks on the first 100 amino acids: (b1) Mycoplasma putrefaciens, (b2) Mortierella verticillata, (b3) Blattabacterium, (h1) Aeropyrum pernix, (h2) Acidianus hospitalis, (h3) Acidilobus saccharovorans.Cluster analysis of the 4th short Haar wavelet transform of a 4000-length random sequence (left) and its 2000-length random walk (right): (a)(α,α*); (b) (β00,β*00); (c) (β01,β*01); (d) (β11,β*11).
(a)(b)(c)(d)


## 5.1. Random Walks

Random walk on the complex sequenceYN is defined as the series ZN={zn}n=1,…,N(49)zn=def∑k=1,…,nyk,n=1,…,N
which is the cumulative sum(50){y1,y1+y2,…,∑s=1nys…,∑s=1Nys}.
When yk=ρ(xk) with xk∈𝒜ℓ and Xk∈SN we will properly call these walks as DNA walk. When the yk are randomly generated we will call them random walks.By remembering the definition of frequencies, DNA walk is the complex value signal{Zn}n=0,…,N-1 with(51)zn=(R[zn],I[zn])=(an-gn)+(tn-cn)i,zn∈C1,
where the coefficients an,gn,tn,cn given by (12) fulfill the condition (31).If we compare the DNA walks (Figure11) some primitive archaea such as h3  are very similar to a random walk (Figure 13). In particular archaea seem to grow less than other bacteria (with the exception of b2).Walks on the first 200 nucleotides: (b1) Mycoplasma putrefaciens, (b2) Mortierella verticillata, (b3) Blattabacterium, (h1) Aeropyrum pernix, (h2) Acidianus hospitalis, and (h3) Acidilobus saccharovorans.It is interesting also to notice that the random walks on amino acids (Figure12) show that more evolved organisms have some “periodic” behavior, while the absolute value of walks on archaea is growing fast.Absolute value of walks on the first 100 amino acids: (b1) Mycoplasma putrefaciens, (b2) Mortierella verticillata, (b3) Blattabacterium, (h1) Aeropyrum pernix, (h2) Acidianus hospitalis, (h3) Acidilobus saccharovorans.Cluster analysis of the 4th short Haar wavelet transform of a 4000-length random sequence (left) and its 2000-length random walk (right): (a)(α,α*); (b) (β00,β*00); (c) (β01,β*01); (d) (β11,β*11).
(a)(b)(c)(d)


## 6. Wavelet Analysis

Wavelet analysis is a powerful method extensively applied to the analysis of biological signals [12, 43–45] aiming to single out the most significant parameters of complexity and heterogeneity in a time series and, in particular, in a DNA sequence. This method is based on the analysis of wavelet coefficients which are obtained by the wavelet transform.We will consider in the following the Haar wavelet basis (see, e.g., [3, 4, 29]) made by scaling functions:(52)φkn(x)=def2n/2φ(2nx-k),(0≤n,0≤k≤2n-1),φ(2nx-k)={1,x∈Ωkn,Ωkn=def[k2n,k+12n),0,x∉Ωkn,
and the Haar wavelets:(53)ψkn(x)=def2n/2ψ(2nx-k),‖ψkn(x)‖L2=1,ψ(2nx-k)=def{-1,x∈[k2n,k+1/22n),1,x∈[k+1/22n,k+12n),(0≤n,0≤k≤2n-1),0,elsewhere.Thediscrete Haar wavelet transform is the N×N matrix 𝒲N:𝕂N⊂ℓ2→𝕂N⊂ℓ2 which maps the vector(54)Y≡{Yi},(i=0,…,2M-1,2M=N<∞,M∈N)
into the vector of wavelet coefficientsβN={α,βkn}:(55)WNY=βN,βN=def{α,β00,…,β2M-1-1M-1},Y=def{Y0,Y1,…,YN-1},(2M=N).The matrix𝒲N can be easily computed by some recursive product [3, 4, 13, 29, 46] so that with N=4,M=2, we have [3, 4, 29](56)W4=(12121212-12-121212-12120000-1212).From (55) with M=2,N=4, by explicit computation, we have(57)α=14(Y0+Y1+Y2+Y3)
and [1–3, 14](58)β00=12(Y2-Y0+Y3-Y1),β01=12(Y0-Y1),β11=12(Y3-Y2).Thus the first wavelet coefficientα represents the average value of the sequence and the other coefficients β the finite differences. The wavelet coefficients β’s, also called details coefficients, are strictly connected with the first-order properties of the discrete time series.In the following we will consider the short wavelet transform which consists in the subdivision of the DNA sequence into 4-length segments and apply the wavelet transform to each segment. As a result, from theN=2M-length complex vector Y, which is subdivided into 2M-2 segments, the 4-parameter short Haar wavelet transform gives the cluster of points(59)(WpR(Ys),WpI(Ys)),s=0,…,σ=Np,p=4
in the 8-dimensional space ℝ4×ℝ4, that is,(60)(α,α*),(β00,β*00),…,(β2p-1-1p-1,β*2p-1-1p-1),p=4.This algorithm enables us to construct clusters of wavelet coefficients and to study the correlation between the real and imaginary coefficients of the DNA representation and DNA walk. It has been observed [3, 4, 29] that some symmetry arises from the plots of wavelet coefficients of DNA walks.


### 6.1. Cluster Analysis of the Wavelet Coefficients of the Complex DNA Representation

Let us first compute the clusters of wavelet coefficients for the random sequence (48). As can be seen the wavelet coefficients both for the sequence and for its series range in some discrete set of values (see Figure 13).The cluster algorithm applied to the complex representation sequence shows that the values of the wavelet coefficients belong to some discrete finite sets (Figure14).Cluster analysis of the 4th short Haar wavelet transform of the complex representation for a DNA walk on the first 2000 nucleotides of (h1) Aeropyrum, (h2) Acidianus, (h3) Acidilobus saccharovorans in the planes: (a) (α,α*); (b) (β00,β*00); (c) (β01,β*01); (d) (β11,β*11).
(a)(b)(c)(d)It should be noticed that this symmetry on detail coefficients is lost for wavelet transform on longer segments (Figures15, 16 and 17).Cluster analysis of the 8th (left), 16th (middle column), 32th (right) short Haar wavelet transform of the DNA walk on the first 1000 nucleotides of h1 (Aeropyrum) in the planes: (a)(α,α*); (b) (β00,β*00); (c) (β01,β*01); (d) (β11,β*11).
(a)(b)(c)(d)Cluster analysis of the 8th (left), 16th (middle column), 32th (right) short Haar wavelet transform of the DNA walk on the first 1000 nucleotides of h2 (Acidianus) in the planes: (a)(α,α*); (b) (β00,β*00); (c) (β01,β*01); (d) (β11,β*11).
(a)(b)(c)(d)Cluster analysis of the 8th (left), 16th (middle column), 32th (right) short Haar wavelet transform of the DNA walk on the first 1000 nucleotides of h3 (Acidilobus saccharovorans) in the planes: (a)(α,α*); (b) (β00,β*00); (c) (β01,β*01); (d) (β11,β*11).
(a)(b)(c)(d)There follows that DNA sequences have to be considered as Markov chain with short range dependence; in other words any acid nucleic is attached to the chain on the base of a correlation of the previous acid nucleic. In other words, if we look for a dependence rule on the DNA nucleotides this dependence might be summarized by a function as(61)xn+1=f(xn),(n=1,…,N).


## 6.1. Cluster Analysis of the Wavelet Coefficients of the Complex DNA Representation

Let us first compute the clusters of wavelet coefficients for the random sequence (48). As can be seen the wavelet coefficients both for the sequence and for its series range in some discrete set of values (see Figure 13).The cluster algorithm applied to the complex representation sequence shows that the values of the wavelet coefficients belong to some discrete finite sets (Figure14).Cluster analysis of the 4th short Haar wavelet transform of the complex representation for a DNA walk on the first 2000 nucleotides of (h1) Aeropyrum, (h2) Acidianus, (h3) Acidilobus saccharovorans in the planes: (a) (α,α*); (b) (β00,β*00); (c) (β01,β*01); (d) (β11,β*11).
(a)(b)(c)(d)It should be noticed that this symmetry on detail coefficients is lost for wavelet transform on longer segments (Figures15, 16 and 17).Cluster analysis of the 8th (left), 16th (middle column), 32th (right) short Haar wavelet transform of the DNA walk on the first 1000 nucleotides of h1 (Aeropyrum) in the planes: (a)(α,α*); (b) (β00,β*00); (c) (β01,β*01); (d) (β11,β*11).
(a)(b)(c)(d)Cluster analysis of the 8th (left), 16th (middle column), 32th (right) short Haar wavelet transform of the DNA walk on the first 1000 nucleotides of h2 (Acidianus) in the planes: (a)(α,α*); (b) (β00,β*00); (c) (β01,β*01); (d) (β11,β*11).
(a)(b)(c)(d)Cluster analysis of the 8th (left), 16th (middle column), 32th (right) short Haar wavelet transform of the DNA walk on the first 1000 nucleotides of h3 (Acidilobus saccharovorans) in the planes: (a)(α,α*); (b) (β00,β*00); (c) (β01,β*01); (d) (β11,β*11).
(a)(b)(c)(d)There follows that DNA sequences have to be considered as Markov chain with short range dependence; in other words any acid nucleic is attached to the chain on the base of a correlation of the previous acid nucleic. In other words, if we look for a dependence rule on the DNA nucleotides this dependence might be summarized by a function as(61)xn+1=f(xn),(n=1,…,N).


## 7. Conclusions

In this paper archaea DNAs have been studied by focussing on the main parameters for complexity. It has been shown that more or less the main indices for complexity and heterogeneity, such as entropy, fractal dimension, and complexity do not differ too much when we have to classify the complexity of the sequence. However, some DNA sequences look more close to random sequences than others, thus suggesting that the evolution involves a process of complexity reduction: the more evolved a sequence is, the more far from a random distribution it is. In any case seems to be apparently impossible to distinguish between a random sequence and a DNA chain. By using the short wavelet transform instead we have shown that on short range (4-nucleotides) a DNA sequence shows some symmetries that slowly disappear by increasing the length of the analysed segment. Moreover, more evolved organisms have a more symmetrical distribution of wavelet coefficients.



---

*Source: 673934-2012-03-11.xml*