 Research article
 Open Access
 Published:
Hierarchical structure of cascade of primary and secondary periodicities in Fourier power spectrum of alphoid higher order repeats
BMC Bioinformatics volume 9, Article number: 466 (2008)
Abstract
Background
Identification of approximate tandem repeats is an important task of broad significance and still remains a challenging problem of computational genomics. Often there is no single best approach to periodicity detection and a combination of different methods may improve the prediction accuracy. Discrete Fourier transform (DFT) has been extensively used to study primary periodicities in DNA sequences. Here we investigate the application of DFT method to identify and study alphoid higher order repeats.
Results
We used method based on DFT with mapping of symbolic into numerical sequence to identify and study alphoid higher order repeats (HOR). For HORs the power spectrum shows equidistant frequency pattern, with characteristic twolevel hierarchical organization as signature of HOR. Our case study was the 16 mer HOR tandem in AC017075.8 from human chromosome 7. Very long array of equidistant peaks at multiple frequencies (more than a thousand higher harmonics) is based on fundamental frequency of 16 mer HOR. Pronounced subset of equidistant peaks is based on multiples of the fundamental HOR frequency (multiplication factor n for n mer) and higher harmonics. In general, n mer HORpattern contains equidistant secondary periodicity peaks, having a pronounced subset of equidistant primary periodicity peaks. This hierarchical pattern as signature for HOR detection is robust with respect to monomer insertions and deletions, random sequence insertions etc. For a monomeric alphoid sequence only primary periodicity peaks are present. The 1/f^{β}– noise and periodicity three pattern are missing from power spectra in alphoid regions, in accordance with expectations.
Conclusion
DFT provides a robust detection method for higher order periodicity. Easily recognizable HOR power spectrum is characterized by hierarchical twolevel equidistant pattern: higher harmonics of the fundamental HORfrequency (secondary periodicity) and a subset of pronounced peaks corresponding to constituent monomers (primary periodicity). The number of lower frequency peaks (secondary periodicity) below the frequency of the first primary periodicity peak reveals the size of n mer HOR, i.e., the number n of monomers contained in consensus HOR.
Background
Introduction
Repeat sequences are a common feature of genomes [1–3]. The detection and study of periodicity in genomic sequences has been an area of increasing interest. Signal processing approaches to periodicity detection methods are attracting significant attention in genomic DNA investigations of approximate repeats because they are rather robust in the presence of substitutions, insertions and deletions and may identify approximate periodicities in DNA sequences. Different computational techniques have been used: Fourier spectral analysis [4–20], wavelet transform [21], DNA walk analysis [22–25], information theory measures [26–28], informational decomposition [29, 30], quaternionic periodicity transform [31], exactly periodic subspace decomposition [32, 33], portrait method [34], enhance algorithm for distance frequency distribution [35], etc.
Discrete Fourier transformation (DFT) based methods
Spectral analysis employing Discrete Fourier transform is used to reveal periodicity in symbolic sequences, like genomic and protein sequences [7, 9, 14, 16, 17, 20, 36–53], to investigate longrange correlations [4, 5, 54, 55] and to study the problem of sequence similarity [14, 56–62].
DFT identification of approximate repeats
A peak at a frequency f in Fourier power spectrum of base correlations of a given genomic sequence shows a kind of l = 1/f – base periodicity, exact or approximate [14–16, 63]. In the ideal case of perfect periodicity, where a fragment of the length l is exactly repeated N times, periodicity generates a series of l1 equidistant peaks in the power spectrum, at frequencies [14, 16]:
f_{1} = 1/l, f_{2} = 2/l, f_{3} = 3/l, ... f_{l1}= (l1)/l
Approximate repeats, modified by random insertions and/or deletions with respect to perfect repeats, typical for genomic sequences of higher organisms, can often be identified using Fourier transform [14, 16, 17]. This procedure results in a characteristic system of equidistant peaks. However, it was noted that a disadvantage of methods based on Fourier transform may be that in cases of more pronounced deletions or insertions the periodicity cannot be detected, while deletions and insertions are frequent mutational events in genomic sequences [29, 50].
DFT identification of period three hidden periodicity
A sharp peak of period three was found in a search for periodic regularities on a sample set of human exons [5, 9, 10, 22, 54, 60, 64]. The threebase periodicity in exons is caused by unbalanced nucleotide distributions in the three coding positions, while in intron sequences the nucleotides distribute uniformly. The relative height of the corresponding peak in Fourier spectrum is a good discriminator of coding potential and has been used to detect coding regions [9, 14, 37, 45, 49, 65–75].
DFT identification of longrange correlations
Statistical studies of DNA sequences have been instigated by finding of the 1/f^{β}longrange powerlaw correlations in human genomic sequences, indicating the presence of scale invariant structure [4, 5, 22], implying that the underlying system shows fractal properties [25, 76, 77]. The lack of long enough sequences and the use of different methods of estimating the correlations, leading to some results not strictly comparable to each other contributed to controversies regarding findings on longrange correlations, like the presence of these correlations only in noncoding or in all human genomic sequences, and their presence in other organisms [5, 6, 23, 36, 78–83]. Nonstationary analysis of DNA sequences has shown that both coding and noncoding sequences exhibit longrange correlations, with the average spectral exponent of noncoding segments being higher than its counterpart for coding segments [84]. With the availability of large sequences and extended statistical computations, showing powerlaw correlations over four or five orders of magnitude, with exponents which are consistent with previous results obtained analyzing short sequences, such correlations in human DNA, with fractallike scaling, are now commonly accepted [27, 28, 45]. It has been pointed out that the mosaic structure of genome is presumably responsible for longrange correlations [79, 85, 86]. At very low frequencies (for example, f < 10^{6}) the power spectrum flattens out [87–89]. It should be noted that the attribution of longrang correlations exclusively to largescale variations of nucleotide density responsible for 1/f^{β}spectra is not quite correct. Generally, even largescale variations of nucleotide density may produce patterns different from 1/f^{β}spectra.
DFT identification of alphoid higher order repeats (HOR)
Here we investigate the application of Fourier analysis to human alpha satellite tandem repeats and the associated higher order repeats (HORs). Alphoid arrays consist of tandem repeats of alpha satellite monomer unit of approximately 171 bp, which form chromosomespecific higher order repeats (HOR) or monomeric organization consisting of diverged monomers [90–104]. Alpha satellite monomers within HOR exhibit substantial mutual sequence divergence (20–40%), while HORs exhibit much lower mutual divergence (< 5%) [98]. Such a case is interesting for Fourier analysis because it has a twolevel hierarchy of approximate homology.
Alpha satellite DNA is characterized by many levels of hierarchical organization in genomes, from suprachromosomal families to chromosomespecific subsets, to polymorphic variation within these subsets [90–103]. The higher order repeat organization is consistent with linear sets of diverged monomers becoming the unit of crossingover during the process of sequence homogenization. The HOR units of alpha satellite monomers are organized in largely chromosome specific manner. The centromere of each human chromosome is characterized by one or more subsets of distinct alpha satellite HOR units. Analyses have revealed the presence of up to several thousand repeat units arranged in an apparently uninterrupted fashion in the centromere and forming arrays of several million base pairs. Alpha satellite HORs have been studied using restriction enzymes that cut higher order repeat unit [98, 101]. Recently, HORs and monomeric alpha satellites have been studied by computational analysis of genomic sequences from the NCBI genome assembly [104–108].
As a case study we consider a 16 mer HOR at the loci D7Z2 and D7Z1 in human chromosome 7 [94–97]. In [104] 16 mers were identified by DOTTER analysis; the presence of 16 mer was reported, but detailed HOR structure was not presented. In detailed computational studies of genomic sequence of the 193277bp clone AC017075.8 (contig NT_023603.5), the 46 complete and 14 incomplete copies of 16 mer alphoid HOR were identified in the central domain (positions 31338 to 177434, total length 148147 bp) [105–107]. Preliminary study of power spectra discussed the general pattern and the signaltonoise ratio [17]. These HOR copies are highly homologous (divergence from consensus less than 0.6% on the average), while divergence among monomers within each HOR copy is sizeable (20% on the average). (In accordance with common practice, monomer deletions or insertions, which appear in some HOR copies, are not taken into account in calculating divergence among HOR copies.) Such a long genomic sequence enables a highly precise determination of higher order periodicities. In the front domain of genomic sequence (31337 bp) and in the back domain (15843 bp), 199 alpha satellite monomers are present which are not organized into HORs and therefore are all mutually divergent by 20% or more. Only 29% of this bordering domain is not of alpha satellite type.
Our goal is to investigate the periodicities in the short, medium, and longrange order, related both to less homologous alphoid monomeric pattern and to more homologous alphoid higherorder repeats and to correlate the two levels of periodicity, primary (basic monomer periodicity) and secondary (HOR periodicity).
Results and discussion
DFT identification of HOR in AC017075.8 based on quartic mapping
The genomic sequence AC017075.8 (193277 bp) from chromosome 7 was transformed into numerical sequence using quartic mapping (Eq. 5) with parameters (Eq. 6) (see section Methods). The AC017075.8 sequence is used as a case study for the use of DFT method for interplay of monomeric and HOR repeats. In general, regions containing higher order repeat sequences can be located through the sliding window analysis, similarly as used in [16] for primary periodicity sequences. Analyzing complete nucleotide sequence we found domains having different repeat pattern, the central HOR domain and the bordering domains (front and back domains), in accordance with identifications obtained using other methods [94–97, 105–108].
The lowfrequency part f < 0.01 bp^{1} of the power spectrum of HOR domain in AC017075.8 is displayed in Figure 1a) and the power spectrum up to the frequency 0.15 bp^{1} in Figure 1b). The computed power spectrum at low frequencies shows equidistant peaks at frequencies
f_{ n }= n·f_{1}, n = 1, 2, 3, ...
The fundamental frequency f_{1} corresponds to the 2734bp HOR. (Due to truncation of data set and the associated precision limit of 7.6·10^{6}, a more precise value can be deduced from the systematic of higher multiples).
These equidistant peaks are identified over a very broad interval, up to n ≈ 1000. In fact, all prominent peaks above the white noise background in the power spectrum are multiples of the fundamental frequency f_{1}. We note that such an extremely regular pattern can be rarely found even in the most regular dynamical systems in physics and engineering.
In addition to the standard spectral density S_{ f }(f_{ n }) at frequency f_{ n }(square of absolute value of Fourier amplitude), we define an effective spectral density
renormalized in order to increase the relative weight of lowfrequency with respect to high frequency peaks. The effective values S_{ eff }corresponding to frequencies f_{1}, f_{2} = 2f_{1}, f_{6} = 6f_{1}, and f_{16} = 16f_{1} are 2.025, 0.973, 0.895, and 6.584, respectively.
The prominent peak at the frequency f_{16} = 0.005852 bp^{1} corresponds to approximately 171 bp length. More precisely, 1/f_{16} = 170.88 bp. It corresponds to a set of alpha satellite monomers which constitute consensus HOR (nine 171bp, five 170bp, one 172bp, and one 173bp copy variants). Alternatively, the HOR period 1/f_{1} = 2734 bp could be also expressed as multiple of monomer period 1/f_{16} = 171 bp. The lowfrequency peaks at f_{1}, f_{2}, ..., f_{15} are subharmonics of the monomer frequency f_{16}.
In the frequency region above the monomer frequency f_{16} (Figure 1b), within the set of multiple frequencies nf_{1} (n > 16) we find a prominent subset of higher harmonics at frequencies that are multiples of the monomer frequency f_{16} : 2f_{16} , 3f_{16}, 4f_{16}, ... This subset with band head at the frequency f_{16} will be referred to as monomeric band.
Fourier analysis works well enough for studying relatively short periodicities while the statistical significance of longer periodicities will be decreased by the presence of shorter periodicities [29]. Thus, the statistical significance of longer periods was predicted to be a sort of smeared through statistical significance of shorter periods, i.e., for harmonics with longer periods the damping effect may be more pronounced [29]. We show here that the DFT method is applicable to alphoid HORs up to very long periodicities (up to several kilobases).
Generally, the fundamental frequency for equidistant pattern in the power spectrum corresponds to the periodicity of highest order in a given sequence, i.e., to the period of HOR (secondary periodicity). Specifically, in the HOR domain of AC017075.8 the fundamental frequency in power spectrum corresponds to the period of HOR consensus unit, 2734 bp.
Although the HOR copies are much more homologous to each other than the constituent alpha satellite monomers among themselves, the number of monomers corresponding to primary periodicity (at frequencies f_{16}, 2f_{16}, 3f_{16}, ...) is much higher than the number of HOR copies corresponding to secondary periodicity (at frequencies f_{1}, f_{2}, f_{3}, ...), and therefore the peaks of primary periodicity have higher spectral strengths.
Robustness of DFT results for hierarchical structure of cascade of primary and secondary periodicities using different genomic into numerical sequence mapping
The difficulty with DFT approach may be dependence on a particular labeling adopted. For example, some of the relevant harmonic structure can be hidden (or exposed) by the symbolictonumeric mapping [111]. To check the required mapping invariance, we investigate whether the hierarchical periodic pattern shown in this paper is robust with respect to a particular choice of procedure for transforming symbolic to numerical sequence.
Test computation for quartic mapping deduced from systematic of purine/pyrimidine and strong/weak bond characteristics
To test robustness of hierarchical structure obtained in Figures 1a) and 1b), we have first computed the power spectrum of HOR domain in AC017075.8, using the quartic mapping (Eq. 5) with parameters (Eq. 8) (Figures 2a) and 2b)). The quartic parametrization (Eq. 8) was based on combined consideration of purine/pyrimidine and strong/weak characteristics of nucleotides [113]. Up to an overall normalization, this spectrum is practically identical as obtained in the computation in Figure 1. This is understandable because of linear relation between mapping parameters (Eq. 6) and (Eq. 8).
Test computation for quartic mapping deduced from reduced dimensionality of frequency spectrum in symmetric manner
A further test was performed using quartic mapping (Eq. 5) with reduced dimensionality of the frequency spectrum representation from four to three with parameters (Eq. 9)–(Eq. 11) from [37]. The computed spectrum in Figures 3a) and 3b) shows a similar pattern of hierarchy of primary and secondary periodicity peaks as in Figure 1, confirming robustness of the method.
Test computation by summing the squares of Fourier transforms of indicator sequences
Finally, we test the robustness of hierarchical primary and secondary periodicity pattern by computing total power spectrum obtained by summing squares of Fourier transform of indicator sequences (Eq. 4). The resulting power spectrum in Figure 4 shows a similar hierarchical pattern as in Figure 1, confirming robustness of the method.
Hierarchical primary and secondary periodicity pattern in perfect HOR sequence
The HOR sequence from AC017075.8 in chromosome 7, studied in Figures 1, 2, 3, 4, is characterized by a low divergence among 54 HOR copies in the sequence of only a few percent [105–108]. Here we construct an exact HOR sequence, with divergence among copies equal to zero, by forming a sequence of 54 identical HOR copies, equal to the 2734bp consensus HOR corresponding to AC017075.8 in chromosome 7 [108]. The resulting power spectrum (Figure 5) shows a much more pronounced hierarchical secondary periodicity pattern than obtained for realistic HOR sequence in Figure 1. In this way analysis was extended to give some feel of how this perfect case appears when the periodicity is disrupted for a realistic genomic sequence.
Robustness of power spectrum for hierarchical structure of cascade of primary and secondary periodicities for imperfect HORs
In the next step we have investigated the robustness of the hierarchical periodicity pattern with an increase of imperfection in the HOR sequence. This is shown by random insertion of increased length inserted into the HOR sequence AC017075.8. Power spectra are presented for: 10000bp random insertion (Figure 6a) and 6b)), 30000bp random insertion (Figure 6c) and 6d)), and 60000bp random insertion (Figure 6e) and 6f)). It is seen that the level of noise increases with increase of insertion length, but even in the case of 60000bp random insertion (which is 40% of the total length of HOR copies) the hierarchical structure of primary and secondary periodicity can be identified (Figures 6e) and 6f)).
Noisy power spectrum of a random artificial sequence
In order to test that the hierarchical primary and secondary periodicity pattern is not a numerical artifact, we have computed the power spectrum corresponding to a random sequence generated by random number generator, having the same length as the HOR sequence in AC017075.8 (148147 bp). Computation is performed using quartic mapping with paramerization (Eq. 8). From this power spectrum, shown in Figure 7, it is seen that the computational method does not generate hierarchical periodicity.
Power spectrum pattern of monomeric alphoid domain in AC017075.8
The lowfrequency power spectrum (f < 0.01 bp^{1}) of combined front and backdomains of AC017075.8 is displaced in Figure 8a) and the higherfrequency section in Figure 8b). A significant difference with respect to the central HOR domain is seen in the lowfrequency region (Figure 8a): there are no prominent peaks below the frequency f_{16} (1/171 bp^{1}). In the front and backdomains there is no peak corresponding to 16 mer HOR (at frequency f_{1} in Figure 1a), as well as to the multiples of f_{1}: f_{2} = 2f_{1}, f_{3} = 3f_{1}, ..., f_{15} = 15f_{1}, at frequencies below f_{16} which corresponds to the 171bp monomer. In that case the frequency of the lowest peak in the power spectrum corresponds to the period 171 bp of consensus alpha monomer and the power spectrum contains only the monomeric band (primary periodicity). This reveals that HOR is absent in the front and back domains. A tandem of alpha satellite monomers, not organized into HORs, is referred to as monomeric [96, 97]. As seen from our results, the noise is stronger in monomeric domains (front and backdomains) than in the central HOR domain. On the other hand, the equidistant pattern associated with prominent peaks above the frequency f_{16} (Figure 8b), for monomeric domain is rather similar as for HOR domain (Figure 1b).
Absence of lowfrequency 1/f^{β}noise in DFT power spectrum
The 1/f^{β}– noise is absent in the lowfrequency region of power spectrum of AC017075.8, both in the central HOR domain (Figure 1a) and in the monomeric front and backdomains (Figure 8a). This result is in accordance with expectations, because the sequence mainly consists of approximate repeats, without sizeable sequencewide base composition fluctuations. Previously, some cases of absence of long range correlations in repeat sequences have been found. For example, in a sequence for beta globin on human chromosome 11 (HUMHBB, 73326 bp) two 6kb segments without longrange correlations were identified, both including stretches of repetitive DNA [76].
The A+T fraction in the AC017075.8 sequence is almost constant along the sequence (Figure 9). Using bins of 1 kb, the calculated fraction of A+T nucleotides is 0.625 ± 0.008%, with small fluctuations around the average value. This homogeneity of nucleotide density is in accordance with expectations that the 1/f^{β}noise is related to varying ratio of pyrimidines and purines, or other nucleotide combinations, at base positions along DNA sequence [24, 81].
For comparison, to show that the DFT power spectrum method used here identifies the lowfrequency 1/f^{β}– noise, if present, we display the power spectrum computed for contig NT_004434.18 from chromosome 1 (Figure 10). This contig of about 1 Mb lies outside of (peri)centromeric region and is characterized by the presence of genes and absence of HORs. The power spectrum computed using quartic mapping at parametrization (Eq. 7) clearly shows the presence of lowfrequency 1/f^{β}– noise.
Rank ordering for harmonics in alpha monomeric spectrum
In the power spectrum of our case study for genomic sequence in HOR domain the equidistant peaks corresponding to multiples of monomer frequency f_{16} are sizably stronger than the other peaks (Figure 1b). Among the lowfrequency peaks in Figure 1b) the most pronounced peaks are 10f_{16}, 12f_{16}, and 14f_{16}. The corresponding lengths are approximately 17 bp, 14 bp, and 12 bp, respectively. The 14bp length may be related to the highest frequency of appearance of the 6bp key string TTTTGA at the distance of 14 bp between two neighboring key strings. However, in general, the chosen mapping may influence the rank ordering of harmonics, as seen by comparing their relative heights in Figures 1, 2, 3, 4, 5. Thus, the effect of parameter choice for symbolictonumeric transformation may overshadow the effect of hidden genomic substructure.
Absence of periodicity three in power spectrum of HORs
In previous investigations of Fourier power spectra of coding DNA sequences a major peak was found at the frequency f = 1/3 bp^{1}, related to the codon structure [5, 14, 37]. In the present case of a segment with entirely noncoding sequence, no peak appears at f = 1/3 bp^{1} (Figure 11). This is in accordance with previous conclusions that the period3 feature is usually lacking or is weak in noncoding regions [7, 9, 37, 39, 41, 66]. For comparison, using quartic mapping we computed the power spectrum of CDS from the gene DNAH11 in chromosome 7 (Figure 12). In this power spectrum obtained by computation using quartic mapping (Eq. 8), at the frequency f = 1/3 bp^{1} a pronounced peak is present.
Use of power spectrum for identification of hierarchical primary and secondary periodicity pattern in chromosome 1
The computation of power spectrum, shown here for the test case of 16 mer HOR in chromosome 7, can be extended for HOR identification and study in other chromosomes as well. As an example, we present in Figure 13 the power spectrum computed for contig NT_077389.3 in chromosome 1, using quartic mapping with parametrization (Eq. 8). Here we find a hierarchical pattern of primary and secondary periodicity (11 mers).
Conclusion
We have demonstrated that DFT is a robust and efficient method to identify alphoid HORs in alpha satellite domain of genomic sequence. In the case of n mer HOR the lowest peak is at the fundamental frequency 1/(171n bp), which will be referred to as HORfrequency. It is a head of band of equidistant peaks at frequencies equal to consecutive multiples of HORfrequency, i.e., 1/(171n bp), 2/(171n bp), 3/(171n bp), ... This band is referred to as the HORband. Some peaks within the HORband form a strongspectralpower subset with band head at monomerfrequency 1/(171 bp). This subset forms a band of equidistant peaks at frequencies which are multiples of monomerfrequency 1/(171 bp), i.e., it corresponds to the peaks at frequencies 1/(171 bp), 2/(171 bp), 3/(171 bp), .... This subband is referred to as the monomericband.
In the case of monomeric alpha satellites (not organized into HOR) the lowest peak is at the monomerfrequency 1/(171 bp). It is a head of monomericband built from peaks at frequencies 1/(171 bp), 2/(171 bp), 3/(171 bp), ....
DFT was applied here in the case study of genomic sequence AC017075.8 (193277 bp) from centromeric region in human chromosome 7. The central domain of AC017075.8 consists of 16mer alphoid HOR copies. Thus the frequency of the lowest peak in the power spectrum (HORfrequency) is 1/(171·16 bp). We identified in the power spectrum as many as one thousand peaks at frequencies equal to multiples of HORfrequency, forming a HORband. Among these peaks in the HORband a subset of peaks at frequencies 1/(171 bp), 2/(171 bp), 3/(171 bp), .... is characterized by pronounced spectral power and represents the monomericband. This reveals hidden periodicities in the 171bp monomer, i.e., a hierarchy of periodicities within the monomer sequence. Power spectra of both the HOR region and of the monomeric region show this pattern of hidden higher frequencies.
The case study shows that DFT is robust in detecting approximate HORs, even in the presence of substantial sequence insertions and deletions.
Additionally, the applicability of DFT method was shown for chromosome 1, where a hierarchical pattern of 11 mer HOR is present.
Computing DFT power spectra for anonymous genomic sequence using sliding windows for bins of about 50 kb and step size of about 10 kb provides an easily recognizable hierarchical twolevel equidistant pattern in the power spectrum as signature of presence of HOR and gives a simple method to determine the size of HOR.
Methods
Discrete Fourier transform of genomic sequence
To apply DFT, one should first represent genomic sequence, a symbol sequence over the alphabet {A, T, G, C}, as a numerical sequence reflecting the characteristics of the symbol sequence. Several approaches have been used for solving the problem of transformation of a symbol sequence to numerical sequence.
A common mapping scheme is to decompose genomic sequence into four component indicator sequences. These binary indicator sequences, u_{A(m)}, u_{T(m)}, u_{C(m)}, and u_{G(m)}, take the value of either 1 or 0 at position m depending on whether the corresponding character is present or absent at that location, respectively. These indicator sequences were analyzed by respective Fourier transforms [5, 9, 39, 40, 55]. For pure DNA character strings (i.e., without assigning numerical values), to the binary indicator sequences u_{A(m)}, u_{T(m)}, u_{C(m)}, and u_{G(m)} correspond the DFT sequencesu_{A(k)}, u_{T(k)}, u_{C(k)}, and u_{G(k)}
respectively, providing a fourdimensional representation of the frequency spectrum of the character string. The quantity obtained by summing the squares of the Fourier transform of indicator sequences:S(k) = u_{ A }(k)^{2} + u_{ T }(k)^{2} + u_{ C }(k)^{2} + u_{ G }(k)^{2}
is used as a measure of the total spectral content of DNA character string at frequency k [9, 37, 39, 40, 111].
Fourier transform of a nucleotide sequence was represented also by sum of pure sequences (Eq. 3) or by their product [15, 109]. A single binary sequence was used by mapping genomic sequence into purine/pyrimidine representation [22], or into weak bond/strong bond representation [109]. Alternatively, mapping of DNA symbolic sequence into a set of quaternions could be utilized via the use of quaternionic Fourier transform [31].
A quartic mapping of genomic into numerical sequence of length N was performed by mapping each symbol to a number [17, 37, 111]:
x(m) = au_{A(m)} + tu_{T(m)} + cu_{C(m)} + gu_{G(m)}, m = 0, 1, 2, ..., N1 (5)
where a, t, c, and g are numerical values assigned to the characters A, T, C, and G, respectively.
We define the quartic map by ordering numbers of nucleotides with increasing frequency in the sequence AC017075.8 (corresponding to the orientation ) in chromosome 7, which is used for our case study:
a = 4, t = 3, c = 2, g = 1 (6)
These values are in accordance with ordering of nucleotides with decreasing frequencies in the HOR region, and therefore they are biased in favor of A and T.
When using sequences from the Build 36.2 assembly (corresponding to the orientation +), the corresponding parametrization for quartic mapping is complement to (Eq. 6):
t = 4, a = 3, g = 2, c = 1 (7)
In [113] the purine/pyrimidine and strong/weak bond properties of the four kinds of nucleotides were considered. The point (1,1) was used to represent nucleotide C corresponding to its pyrimidine and strong bond properties; the point (1,1) to represent nucleotide G corresponding to its purine and strong bond properties; the point (1,1) to represent nucleotide A corresponding to its purine and weak bond properties; and the point (1,1) to represent nucleotide T corresponding to its pyrimidine and weak bond properties. Then the vectors connecting the origin to the four points (1,1), (1,1), (1,1) and (1,1) have the rotational angles π/4, 3π/4, 5π/4,7π/4 with the xaxis and correspondingly the map defined [113]:
a = 7, t = 5, c = 3, g = 1 (8)
These quartic map parameters are linearly related to parameters in (Eq. 6), b = 2b'1 (b and b' stand for the corresponding nucleotides in (Eq. 8) and (Eq. 6), respectively).
In [37] the dimensionality of the frequency spectrum representation was reduced from four to three in a symmetric manner with respect to all four components. Three numerical sequences ξ_{r}, ξ_{g}, ξ_{b} were defined from the corresponding coefficients (a_{r}, t_{r}, c_{r}, g_{r}), (a_{g}, t_{g}, c_{g}, g_{g}), (a_{b}, t_{b}, c_{b}, g_{b}) by considering the four threedimensional vectors having magnitude equal to 1 and pointing to the four directions from the center to the vertices of regular tetrahedron:
from which the DFTs are calculated.
In all computations the DFT was computed using Fast Fourier Transform (FFT) computer program [115] with the 1/$\sqrt{N}$ normalization.
A search for regions of higher order repeats in anonymous sequence, without prior knowledge on its structure, can be performed by sliding window analysis, similarly as used in Spectral Repeat Finder [16]. Once a region of HOR structure is detected, a more precise edge detection of HOR region can be determined by performing more precise local search using smaller step size.
Abbreviations
 HOR:

Higher Order Repeat
 KSA:

Key String Algorithm
 DFT:

Discrete Fourier Transform.
References
 1.
Haubold B, Wiehe T: How repetitive are genomes? BMC Bioinformatics 2006, 7: 541.
 2.
Rocha EPC, Danchin A, Viari A: Functional and evolutionary rules of long repeats in prokaryotes. Res Microbiol 1999, 150: 725–733.
 3.
Gregory TR: Synergy between sequence and size in largescale genomics. Nature Rev Genet 2005, 6: 699–708.
 4.
Li W, Kaneko K: Longrange correlation and partial 1/ f^{α}spectrum in a noncoding DNA sequence. Europhys Lett 1992, 17: 655–660.
 5.
Voss RF: Evolution of longrange correlations and 1/f noise in DNA base sequences. Phys Rev Lett 1992, 68: 3805–3808.
 6.
Borštnik B, Pumpernik D, Lukman D: Analysis of apparent 1/ f^{α}spectrum in DNA sequences. Europhys Lett 1993, 23: 389–394.
 7.
Chechetkin VR, Turygin AY: Search of hidden periodicities in DNA sequences. J Theor Biol 1995, 175: 477–494.
 8.
Coward E: Equivalence of two Fourier methods for biological sequences. J Math Biol 1997, 36: 64–70.
 9.
Tiwari S, Ramachandran S, Bhattacharya S, Ramaswami R: Prediction of probable genes by Fourier analysis of genomic sequences. Comp Appl Biosci 1997, 13: 263–270.
 10.
Kutuzova GI, Frank GK, Makeev VY, Esipova NG, Polozov RV: Analysis of nucleotide sequences – periodicities in Ecoli promoter sequences. Biofizika 1999, 42: 354–362.
 11.
Guharay S, Hunt BR, Yorke JA, White OR: Correlations in DNA sequences across the three domains of life. Physica D 2000, 146: 388–396.
 12.
Nagai N, Kuwata K, Hayashi T, Kuwata H, Era S: Evolution of the periodicity and the selfsimilarity in DNA sequence: A Fourier analysis. Jap J Physiol 2001, 51: 159–168.
 13.
Yu ZG, Anh V, Lau KS: Measure representation and multifractal analysis of complete genomes. Phys Rev E 2001, 64: 031903.
 14.
Lobzin VV, Chechetkin VR: Order and correlations in genomic DNA sequences. The spectral approach. Uspekhi Fizicheskih Nauk 2000, 170: 57–81.
 15.
Tran TT, Emanuele VA II, Zhou GT: Techniques for detecting approximate tandem repeats in DNA. In Proceedings of International Conference for Acoustics, Speech and Signal Processing (ICASSP). Volume 5. Montreal, Canada; 2004:449–452.
 16.
Sharma D, Isaac B, Raghava GPS, Ramaswamy R: Spectral Repeat Finder (SRF): identification of repetitive sequences using Fourier transformation. Bioinformatics 2004, 20: 1405–1412.
 17.
Paar V, Pavin N, Basar I, Rosandić M, Luketin I, Durajlija Žinić S: Spectral densities and frequencies in the power spectrum of higher order repeat alpha satellite in human DNA molecule. Cro Chem Acta 2004, 77: 73–81.
 18.
Pop PG: Spectral techniques in finding DNA approximate tandem repeats. Automation, Quality and Testing, Robotics, IEEE International Conference 2006, 2: 441–444.
 19.
Du L, Zhou H, Yan H: OMWSA: detection of DNA repeats using moving window spectral analysis. Bioinformatics 2007, 23: 631–633.
 20.
Ramakrishna R, Srinivasan R: Gene identification in bacterial and organellar genomes using GeneScan. Computers Chem 1999, 23: 165–174.
 21.
Arneodo A, Bacry E, Graves PV, Muzy JF: Characterizing longrange correlations in DNA sequences from wavelet analysis. Phys Rev Lett 1995, 74: 3293–3296.
 22.
Peng CK, Buldyrev SV, Goldberger AL, Havlin S, Sciortino F, Simons M, Stanley HE: Longrange correlations in nucleotide sequences. Nature 1992, 356: 168–170.
 23.
Nee S: Uncorrelated DNA walks. Nature 1992, 357: 450.
 24.
ChatzidimitriouDreismann CA, Larhammar D: Longrange correlations in DNA. Nature 1993, 361: 212–213.
 25.
Peng CK, Buldyrev SV, Goldberger AL, Havlin S, Mantegna RN, Simons M, Stanley HE: Statistical properties of DNA sequences. Physica A 1995, 221: 180–192.
 26.
Herzel H, Grosse I: Measuring correlations in symbol sequences. Physica A 1995, 216: 518–542.
 27.
Holste D, Grosse I, Herzel H: Statistical analysis of the DNA sequence of human chromosome 22. Phys Rev E 2001, 64: 041917.
 28.
BernaolaGalvan P, Carpena P, RomanRoldan R, Oliver JL: Study of statistical correlations in DNA sequences. Gene 2002, 300: 105–115.
 29.
Korotkov EV, Korotkova MA, Frenkel FE, Kudryashov NA: The informational concept of searching for periodicity in symbol sequences. Mol Biol 2003, 37: 372–376.
 30.
Korotkov EV, Korotkova MA, Kudryashov NA: Information decomposition method to analyze symbolical sequences. Phys Lett A 2003, 312: 198–210.
 31.
Brodzik AK: Quaternionic periodicity transform: an algebraic solution to the tandem repeat detection problem. Bioinformatics 2007, 23: 694–700.
 32.
Mauresan DD, Parks TW: Orthogonal, exactly periodic subspace decomposition. IEEE Transactions on Signal Processing 2003, 51: 2270–2279.
 33.
Gupta R, Sarthi D, Mittal A, Singh K: A novel signal processing measure to identify exact and inexact tandem repeat patterns in DNA sequences. EURASIP J Bioinform Syst Biol 2007, 43596.
 34.
Yu ZG, Jiang P: Distance, correlation and mutual information among portraits of organisms based on complete genomes. Phys Lett A 2001, 286: 34–46.
 35.
Pizzi E, Liuni S, Frontali C: Detection of latent sequence periodicities. Nucleic Acids Res 1990, 18: 3745–3752.
 36.
Buldyrev SV, Goldberger AL, Havlin S, Mantegna RN, Matsa ME, Peng CK, Simons M, Stanley HE: Longrange correlation properties of coding and noncoding DNA sequences: GenBank analysis. Phys Rev E 1995, 51: 5084–5091.
 37.
Anastassiou D: Frequencydomain analysis of biomolecular sequences. Bioinformatics 2000, 16: 1073–1081.
 38.
Widom J: Shortrange order in two eukaryotic genomes: relation to chromosome structure. J Mol Biol 1996, 259: 579–588.
 39.
Silverman BD, Linsker R: A measure of DNA periodicity. J Theor Biol 1986, 118: 295–300.
 40.
Li W, Marr TG, Kaneko K: Understanding longrange correlations in DNA sequences. Physica D 1994, 75: 392–416.
 41.
Trifonov EN: 3, 10.5, 200 and 400base periodicities in genome sequences. Physica A 1998, 249: 511–516.
 42.
Zhang MQ: Computational prediction of eucaryotic proteincoding genes. Nature Genet 2002, 3: 698–710.
 43.
Chechetkin VR, Lobzin VV: Nucleosome units and hidden periodicities in DNA sequences. J Biomol Structure Dynamics 1998, 15: 937–947.
 44.
Kim BR, Littel RC, Wu R: Clustering periodic patterns of gene expression based on Fourier approximations. Curr Genomics 2006, 7: 197–203.
 45.
Som A, Sahoo S, Mukhopadhyay I, Chakrabarti J, Chaudhury R: Scaling violation in coding DNA. Europhys Lett 2003, 62: 271–277.
 46.
Yin C, Yau SST: A Fourier characteristic of coding sequences: origins and a new nonFourier approximation. J Comput Biol 2005, 12: 1153–1165.
 47.
Fuentes AR, Ginori JVL, Abalo RG: Detection of coding regions in large DNA sequences using the short time Fourier transform with reduced computational load. Progress in pattern recognition, image analysis and applications, Proceedings of the Lecture Notes in Computer Science 2006, 4225: 902–909.
 48.
Lynn AM, Jain CK, Kosalai K, Barman P, Thakur N, Batra H, Bhattacharaya A: An automated annotation tool for genomic DNA sequences using GeneScan and Blast. J Genet 2001, 80: 9–16.
 49.
Hall R, Stern L: A rapid method for illustrating features in both coding and noncoding regions of a genome. Bioinformatics 2004, 20: 982–983.
 50.
Turutina VP, Laskin AA, Kudryashov NA, Skryabin KG, Korotkov EV: Identification of latent periodicity in amino acid sequences of protein families. Biochemistry (Moscow) 2006, 71: 18–31.
 51.
Dodin G, Vandergheynst P, Levoir P, Cordier C, Marcourt L: Fourier and wavelet transform analysis, a tool for visualizing regular patterns in DNA sequences. J Theor Biol 2000, 206: 323–326.
 52.
Jackson JH, George R, Herring PA: Vectors of Shannon information from Fourier signals characterizing base periodicity in genes and genomes. Biochem Biophys Res Commun 2000, 268: 289–292.
 53.
Coward E, Drablos F: Detecting periodic patterns in biological sequences. Bioinformatics 1998, 14: 498–507.
 54.
Chechetkin VR, Knizhnikova LA, Turygin AY: Threequasiperiodicity, mutual correlations, ordering and longrange modulations in genomic nucleotide sequences for viruses. J Biomol Structure Dynamics 1994, 12: 271–299.
 55.
Stanley HE, Buldyrev SV, Goldberger AL, Goldberger ZD, Havlin S, Mantegna RN, Ossadnik SM, Peng CK, Simons M: Statistical mechanics in biology: how ubiquitous are longrange correlations? Physica A 1994, 205: 214–253.
 56.
Felsenstein J, Sawyer S, Kochin R: An efficient method for matching nucleicacid sequences. Nucleic Acids Res 1982, 10: 133–139.
 57.
Benson DC: Fourier methods for biosequence analysis. Nucleic Acids Res 1990, 18: 6305–6310.
 58.
Cleever EA, Overton GC, Searls DB: Fast Fourier transformbased correlation of DNA sequences using complexplane encoding. Comp Appl Biosci 1991, 7: 143–154.
 59.
Arques DG, Michel CJ, Orieux K: Analysis of gene evolution – the software age. Computer Appl Biosci 1992, 8: 5–14.
 60.
Chechetkin VR, Turygin AY: Study of correlations in DNA sequences. J Theor Biol 1996, 178: 205–217.
 61.
Chechetkin VR, Lobzin VV: Study of correlations in segmented DNA sequences – applications to structural coupling between exons and introns. J Theor Biol 1998, 190: 69–83.
 62.
Aghili SA, Agrawal D, El Abbadi A: Sequence similarity search using discrete Fourier and wavelet transformation techniques. Int J Artificial Intelligence Tools 2005, 14: 733–754.
 63.
Li W: The study of correlation structures of DNA sequences: a critical review. Comput Chem 1997, 21: 257–271.
 64.
Trifonov EN, Sussman JL: The pitch of chromatin DNA is reflected in its nucleotide sequence. Proc Natl Acad Sci USA 1980, 77: 3816–3820.
 65.
Herzel H, Trifonov EN, Weiss O, Grosse I: Interpreting correlations in biosequences. Physica A 1998, 249: 449–459.
 66.
Fickett JW: Recognition of protein coding regions in RNA sequences. Nucleic Acids Res 1982, 10: 5303–5318.
 67.
Fickett JW: The gene identification problem: an overview for developers. Comput Chem 1996, 20: 103–118.
 68.
Yin CC, Yau SST: Prediction of protein coding regions by the 3base periodicity analysis of a DNA sequence. J Theor Biol 2007, 247: 687–694.
 69.
Tian YX, Chen C, Zou XY, Qiu JD, Cai PX, Mo JY: Study on period3 behavior of exons. Acta Chimica Sinica 2005, 63: 1215–1219.
 70.
Kotlar D, Lavner Y: Gene prediction by spectral rotation measure: a new method for identifying protein coding regions. Genome Res 2003, 13: 1930–1937.
 71.
Jin J: Identification of protein coding regions of rice genes using alternative spectral rotation measure and linear discriminant analysis. Genomics, Proteomics & Bioinformatics 2004, 2: 167–173.
 72.
Gao J, Qi Y, Cao Y, Tung WW: Protein coding sequence identification by simultaneously characterizing the periodic and random features of DNA sequences. J Biomed Biotechnol 2005, 2: 139–146.
 73.
Lee W, Luo L: Periodicity of base correlation in nucleotide sequence. Phys Rev E 1997, 56: 848–851.
 74.
Dunham I, Shimizu N, Roe BA, Chissoe S: The DNA sequence of human chromosome 22. Nature 1999, 402: 489–495.
 75.
Weiss A, McDonough D, Wertman B, Acakposatchivi L, Montgomery K, Kucherpalati R, Leinwand L, Krauter K: Organization of human and mouse skeletal myosin heavy chain gene clusters is highly conserved. Proc Natl Acad Sci USA 1999, 96: 2958–2963.
 76.
BernaolaGalvan P, RomanRoldan R, Oliver JL: Compositional segmentation and longrange fractal correlations in DNA sequences. Phys Rev E 1996, 53: 5181–5189.
 77.
Stanley HE, Buldyrev SV, Goldberger AL, Goldberger ZD, Havlin S, Mantegna RN, Ossadnik SM, Peng CK, Simons M: Scaling features of noncoding DNA. Physica A 1999, 273: 1–18.
 78.
Buldyrev SV, Goldberger AL, Havlin S, Peng CK, Simons M, Sciortino F, Stanley HE: Long range fractal correlations in DNA. Phys Rev Lett 1993, 71: 1776.
 79.
Karlin S, Brendel V: Patchiness and correlations in DNA sequences. Science 1993, 259: 677–680.
 80.
Pradbu VV, Claverie JM: Correlations in intronless DNA. Nature 1992, 359: 782.
 81.
Larhammar D, ChatzidimitriouDreismann CA: Biological origins of longrange correlations and compositional variations in DNA. Nucleic Acids Res 1993, 21: 5167–5170.
 82.
Arneodo A: What can we learn with wawelets about DNA sequences? Physica A 1998, 249: 439–448.
 83.
ChatzidimitriouDreismann CA, Streffer RM, Larhammar D: A quantitative test of longrange correlations and compositional fluctuations in DNA sequences. Eur J Biochem 1994, 224: 365–371.
 84.
Bouayanaya N, Schonfeld D: Nonstationary analysis of DNA sequences. Proceedings of the 14th Workshop on Statistical Signal Processing 2007, 2007: 200–204.
 85.
Bernardi G, Olofsson B, Filipski J, Zerial M, Salinas J, Cuny G, Meunierrotival M, Rodier F: The mosaic genome of warm blooded vertebrates. Science 1985, 228: 953–958.
 86.
Herzel H, Schmitt AO, Ebeling W: Finitesample effects in sequence analysis. Chaos, Solitons & Fractals 1994, 4: 97–113.
 87.
Vaidyanathan PP, Yoon BJ: The role of signalprocessing concepts in genomics and proteomics. J Franklin Inst 2004, 341: 111–135.
 88.
Sousa Vieira de M: Statistics of DNA sequences: a low frequency analysis. Phys Rev E 1999, 60: 5932–5937.
 89.
O'Neil PV: Advanced Engineering Mathematics. Belmont, California: Wadsworth Publishing Company; 1991.
 90.
Manuelidis L: Chromosomal location of complex and simple repeated human DNAs. Chromosoma 1978, 66: 23–32.
 91.
Jorgensen AL, Bostock CJ, Bak AL: Chromosomespecific subfamilies within human alphoid repetitive DNA. J Mol Biol 1986, 187: 185–196.
 92.
Waye JS, Willard HF: Nucleotide sequence heterogeneity of alpha satellite DNA: a survey of alphoid sequences from different human chromosomes. Nucleic Acids Res 1987, 15: 7549–7572.
 93.
TylerSmith C, Brown WRA: Structure of the major block of alphoid satellite DNA on the human Y chromosome. J Mol Biol 1987, 195: 457–470.
 94.
Waye JS, England SB, Willard HF: Genomic organization of alpha satellite DNA on human chromosome 7: evidence of two distinct alphoid domains on a single chromosome. Mol Cell Biol 1987, 7: 349–356.
 95.
Wevrick R, Willard HF: Physical map of centromeric region of human chromosome 7: relationship between two distinct alpha satellite arrays. Nucleic Acids Res 1991, 19: 2295–2301.
 96.
Wevrick R, Willard VP, Willard HF: Structure of DNA near long tandem arrays of alpha satellite DNA at the centromere of human chromosome 7. Genomics 1992, 14: 912–923.
 97.
Puente de la A, Velasco E, Perez Jurado LA, Hernandez Chico C, Rijke FM, Scherer SW, Raap AK, Cruces J: Analysis of the monomeric alphoid sequences in the pericentromeric region of human chromosome 7. Cytogenet Cell Genet 1998, 83: 176–181.
 98.
Warburton PE, Willard HF: Evolution of centromeric alpha satellite DNA: molecular organization within and between human and primate chromosomes. In Human Genome Evolution. Edited by: Jackson M, Strachan T, Dover G. Oxford: BIOS Scientific; 1996:121–145.
 99.
Lee C, Wevrick R, Fisher RB, FergusonSmith MA, Lin CC: Human centromeric DNAs. Hum Genet 1997, 100: 291–304.
 100.
Alexandrov IA, Kazakov A, Tumeneva I, Shepelev V, Yurov Y: Alphasatellite DNA of primates: old and new families. Chromosoma 2001, 110: 253–266.
 101.
Cho KHA: The Centromere. Oxford: Oxford University Press; 1997.
 102.
Haaf T, Willard HF: Organization, polymorphism, and molecular cytogenetics of chromosomespecific alphasatellite DNA from the centromere of chromosome 2. Genomics 1992, 13: 122–128.
 103.
Romanova LY, Deriagin GV, Mashkova TD, Tumeneva IG, Mushegian AR, Kisselev LL, Alexandrov IA: Evidence for selection in evolution of alpha satellite DNA: the central role of CENPB/pJ α binding region. J Mol Biol 1996, 261: 334–340.
 104.
Rudd MK, Willard HF: Analysis of he centromeric regions of the human genome assembly. Trends Genet 2004, 20: 529–533.
 105.
Rosandić M, Paar V, Basar I: Keystring segmentation algorithm and higherorder repeat 16 mer (54 copies) in human alpha satellite DNA in chromosome 7. J Theor Biol 2003, 221: 29–37.
 106.
Paar V, Pavin N, Rosandić M, Glunčić M, Basar I, Pezer R, Durajlija Žinić S: ColorHOR – novel graphical algorithm for fast scan of alpha satellite higherorder repeats and HOR annotation for GenBank sequence of human genome. Bioinformatics 2005, 21: 846–852.
 107.
Rosandić M, Paar V, Basar I, Glunčić M, Pavin N, Pilaš I: CENPB box and pJα sequence distribution in human alpha satellite higherorder repeats (HOR). Chromosome Res 2006, 14: 735–753.
 108.
Paar V, Basar I, Rosandić M, Glunčić M: Consensus higher order repeats and frequency of string distributions in human genome. Curr Genomics 2007, 8: 93–111.
 109.
Emanuele VA, Tran TT, Zhou GT: A Fourier product method for detecting approximate tandem repeats in DNA. Proceedings of the 13th Workshop on Statistical Signal Processing 2005 IEEE/SP 2005, 1390–1395.
 110.
Azbel MY: Random twocomponent onedimensional Ising model for heteropolymer melting. Phys Rev Lett 1973, 31: 589–592.
 111.
Afreixo V, Ferreira PJSG, Santos D: Fourier analysis of symbolic data: A brief review. Digital Signal Processing 2004, 14: 523–530.
 112.
Zhang R, Zhang CT: Zcurves, an intuitive tool for visualizing and analyzing the DNAsequences. J Biomol Structure Dynamics 1994, 4: 767–782.
 113.
Zhou LQ, Yu ZG, Deng JQ, Anh V, Long SC: A fractal method to distinguish coding and noncoding sequences in a complete genome based on a number sequence representation. J Theor Biol 2005, 232: 559–567.
 114.
Gutierez JM, Rodriguez MA, Abramson G: Multifractal analysis of DNA sequences using novel chaosgame representation. Physica 2001, A300: 271–284.
 115.
Wolfram Mathematica 6: Version Number 6.0.1.0 (32bit).
Acknowledgements
This work was supported by Ministry of Science, Education and Sports of Croatia. The authors express special thanks to anonymous reviewers for their valuable discussions and comments on this manuscript.
Author information
Additional information
Authors' contributions
VP initiated the project and guided the whole work. VP and MR drafted the manuscript. NP, VP and IB implemented and developed the program for power spectrum with quartic mapping. NP, IB, MR, MG and NP were involved in computations and analysis of results. All authors participated in discussions and approved the final manuscript.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
Rights and permissions
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
About this article
Cite this article
Paar, V., Pavin, N., Basar, I. et al. Hierarchical structure of cascade of primary and secondary periodicities in Fourier power spectrum of alphoid higher order repeats. BMC Bioinformatics 9, 466 (2008). https://doi.org/10.1186/147121059466
Received:
Accepted:
Published:
Keywords
 Power Spectrum
 Discrete Fourier Transformation
 Alpha Satellite
 Fourier Power Spectrum
 High Order Repeat