 Research article
 Open access
 Published:
Pervasive properties of the genomic signature
BMC Genomics volumeÂ 3, ArticleÂ number:Â 23 (2002)
Abstract
Background
The dinucleotide relative abundance profile can be regarded as a genomic signature because, despite diversity between species, it varies little between 50 kilobase or longer windows on a given genome. Both the causes and the functional significance of this phenomenon could be illuminated by determining if it persists on smaller scales. The profile is computed from the base step "odds ratios" that compare dinucleotide frequencies to those expected under the assumption of stochastic equilibrium (thorough shuffling). Analysis is carried out on 22 sequences, representing 19 species and comprised of about 53 million bases all together, to assess stability of the signature in windows ranging in size from 50 kilobases down to 125 bases.
Results
Dinucleotide relative abundance distance from the global signature is computed locally for all nonoverlapping windows on each sequence. These distances are lognormally distributed with nearly constant variance and with means that tend to zero slower than reciprocal square root of window size. The mean distance within genomes is larger for protist, plant, and human chromosomes, and smaller for archaea, bacteria, and yeast, for any window size.
Conclusions
The imprint of the global signature is locally pervasive on all scales considered in the sequences (either genomes or chromosomes) that were scanned.
Background
Compositional heterogeneity pervades the genome on all scales [1] and attempts to partition DNA sequences into homogeneous segments give different results depending on the method chosen, even for genomes as small as the 49 kilobases (kb) of bacteriophage lambda [2]. It therefore seems remarkable that dinucleotide relative abundance profiles exhibit local stability in the sense that, when computed for any 50 kb window on a given microbial genome, the profile is about the same as if computed globally from the bulk genomic DNA of the organism [3, 4]. The dinucleotide relative abundance profile, previously called the "general design [5]," can be viewed as a "genomic signature" that reflects a "total net response to selective pressure [6]." Yet the 50 kb window spans the complete genome of bacteriophage lambda and hence it would not be surprising to discover local instability in the profile seen through smaller windows.
Bacterial genomes commonly exhibit an approximate balance between purine (A+G) and pyrimidine (C+T) fractions, and between amino (A+C) and keto (G+T) fractions, when the whole leading strand is examined. Locally, however, these fractions fluctuate from their global averages almost as dramatically as the strong (C+G) and weak (A+T) fractions that do not show global balance. Indeed, these compositional biases, extending over hundreds of kilobases, are strongly correlated with the direction of replication in some species [7]. Perhaps the stability of the dinucleotide relative abundance profile in 50 kb and longer windows conceals highly variable behavior on a shorter scale. By analogy, the variability of the profile within 50 kb windows could have as much functional significance as invariance between such windows.
The stability of the dinucleotide relative abundance profile could result from constraints on dinucleotide stacking energy and DNA helicity, contextdependent mutation pressures, and replication and repair mechanisms [3, 4, 6, 8]. To gauge the influence of these or other factors in contributing to total net response, it again seems pertinent to ask whether stability persists on smaller scales than 50 kb. If it is reasonable to suppose that signature variation between species is due to differences in replication machinery, operating through local DNA structures and base step conformational tendencies [9], then we may expect to find stability on the scale of the machinery. For example, the proximity of dinucleotide relative abundance profiles in two samples could reflect the similarity of enzymes [10] that engage them in their replication processes. When these samples come from the same genome, which encodes the same enzymes, a close proximity should be expected.
Replicational constraints that act on all parts of the genome could also manifest themselves in codon usage preference. The genomic signature may be associated with codon usage to the extent that it embodies these constraints. Given that the signature "pervades both coding and noncoding DNA," its uniformity throughout the genome "cannot be explained by preferential codon usage [9]." The converse proposition, that codon usage is explained by dinucleotide relative abundance, is hardly more defensible, since dinucleotide frequencies do not determine trinucleotide frequencies [11]. Yet the replicational constraints could be a common factor underlying both phenomena. If the global signature pervades an open reading frame (ORF), the trinucleotide frequencies will to some extent express the signature as modulated by the local composition. (Codon usage may exhibit more diversity within genomes by virtue of its sensitivity to composition.) This kind of association would be supported by discovery of signature stability in windows smaller than the average ORF and negated by the breakdown of signature stability below some window size that far exceeds the typical ORF length, unless protein coding sequences are less sensitive to local replicational constraints than dinucleotide relative abundance [12].
Genomes are identifiable by their signatures and dissimilarity between signatures is used to estimate the evolutionary distance between species [4, 6]. If the imprint of the global signature is locally pervasive, down to the scale of the single gene or coding sequence, large deviations on that scale could highlight segments introduced by recent horizontal transfer from another species [13]. Socalled filtration methods, based on dissimilarity measures computed from dinucleotide counts, have been employed for the alignmentfree computation of evolutionary distances between homologous sequences [14, 15]. The "transition matrix method" was a similar technique involving raw counts of amino acid pairs in protein primary sequences [16]. Phylogeny construction based on dinucleotide relative abundance distance ("deltadistance") would seem to have a similar intent although its application has been to whole genomes [9]. Mathematically, however, the use of relative abundance instead of raw abundance (frequencies or counts) is a subtle innovation that compensates for compositional variation. This subtle change has a profound effect and is essential to achieving logical results since whole genome comparisons based on raw abundance often fail to find a close proximity between closely related species.
Can the same improvement be achieved by dinucleotide relative abundance distance calculations in smaller windows? If so, a deltascan of one bacterial genome could detect outliers bearing the imprint of a host or foreign signature. Whether this is practical depends on how the variability of the intragenomic deltadistance grows with shrinking window size. Purely statistical considerations will limit the detectability of small scale deviations. From the investigations of Karlin et al [3, 4, 6, 8, 9, 12, 17] it is clear that deltadistances between 50 kb and larger windows show fluctuations that are small compared to distances between closely related species.
Hypothesis formulation
The assessment of bias in dinucleotide relative abundance begins with the "odds ratios" r _{ xy } = f _{ xy }/f _{ x } f _{ y } where f _{ x } denotes the (normalized) frequency of nucleotide (base) x and f _{ xy } is the frequency of dinucleotide (base step) xy in the leading strand. These ratios compare observed dinucleotide frequencies to those expected from the base composition alone under the assumption of statistical independence (i.e., thorough shuffling of the sequence). When r _{ xy } is greater (less) than one, xy is over(under)represented. The symmetrized version r*_{ xy } is computed from frequencies of the sequence concatenated with its inverted complement. The numbers {r*_{ xy }} comprise the dinucleotide relative abundance profile [3â€“6].
The statistical problem is to test the hypotheses that patterns of dinucleotide over and underrepresentation in a given genome are invariant. Using symbol f for frequencies in windows, let g be used to represent the global frequencies computed for a complete sequence (a genome or chromosome). The hypothesis asserts that r*_{ xy } in any window is approximately equal to a constant. This constant is the global signature c*_{ xy } = g*_{ xy }/g*_{ x } g*_{ y }.
Karlin, Landunga, and Blaisdell [8] assessed homogeneity of the dinucleotide relative abundance profile through the deltadistance Î´* = (1/16) Î£ r*_{ xy }  c*_{ xy }. They provided standards for classifying Î´* in 100 kb windows as follows: "random" (0 < 1000Î´* < 15), "very close" (15 < 1000Î´* < 30), "close" (30 < 1000Î´* < 45), "moderately related" (45 < 1000Î´* < 65), and "distantly related" (65 < 1000Î´* < 95). The upper limit of the "random" range, which typifies thoroughly shuffled sequences, scales as 1/âˆšn, where n is the number of bases in the window.
The local stability of r _{ xy } would be an ancillary result under the assumption of a stationary stochastic process, since then the frequencies converge in probability to fixed limits, f _{ x } â†’ p _{ x } and f _{ xy }/f _{ x } â†’ p _{ xy } as n â†’ âˆž. The simplest case would be a homogeneous Markov chain with base step transition probabilities p _{ xy } and stationary base composition p _{ x }. In this case the differences f _{ x }  p _{ x } and f _{ xy }/f _{ x }  p _{ xy } tend to zero as 1/âˆšn and it is clear that r _{ xy }  p _{ xy }/p _{ y } will vanish at the same rate [18]. Moreover, the globally computed quotient c _{ xy } = g _{ xy }/g _{ x } g _{ y } would be a consistent estimate of p _{ xy }/p _{ y }. Thus the separate terms of Î£ r _{ xy } c _{ xy } tend to zero as 1/âˆšn and the same must obviously apply to Î´*.
We start however with the understanding that the sequence is fundamentally nonstationary, exhibiting statistically significant variations in base frequencies between nonoverlapping windows [1, 2]. Locally or globally estimated Markov models may describe it better than assuming that the bases are independent and identically distributed [19] but they fail to reflect the salient features of natural sequences [20]. For example, Robin and Daudin [21] compared the frequency of a specific motif in the genome sequence of Haemophilus influenzae to Markovian predictions and found that the observed frequencies were everywhere higher than predicted.
This point aside, the stationary Markov analogy provides a useful benchmark in assessing local stability of the genomic signature. If nucleotide sequences behaved like Markov chains, then Î´*âˆšn would not depend on n. We will examine this scaled deltadistance Î´*âˆšn to see if it exhibits any trend. A decreasing ("superMarkov") trend would imply that signature stability emerges as the scale increases and could indicate the breakdown of stability for some window size below 50 kb.
Scope of the investigation
This survey examines 22 sequences from 19 species and 17 genera. The sequences are listed in Table 1 along with serial numbers (SN) and 4letter abbreviations (Abbr). Most of them have been previously studied and found to show stability of the genomic signature in 50 kb windows [9]. The shortest is the 580 kb complete genome of Mycoplasma genitalium. The longest complete sequence is the 7657 kb human chromosome XXII.
The selected sequences, which are not always complete genomes, fall into two main classes, being typically (1) the chromosome that constitutes the largest single element in a prokaryotic genome or (2) one of the chromosomes in a eukaryotic genome. In the first case, for example, is the Borrelia burgdorferi sequence that excludes 21 identified plasmids. The second case is exemplified by Plasmodium falciparum where chromosomes II and III are selected to the exclusion of the other 12. Since our investigation focuses on scaling properties of the genomic signature, it is appropriate to consider sequences comprised of many 50 kb contigs, assuming that variation between such contigs is small compared to variation between species.
The present sample, which spans a wide range of G+C proportion, is hoped diverse enough that any consistent trends and features in the statistical picture it produces cannot be easily attributed to chance. A broader range of sequences, including mitochondrial and large viral genomes, has been surveyed by Karlin et al [3, 4, 6, 8, 9, 12, 17] using 50 kb and larger windows. This investigation, which applies similar methodology to smaller window sizes, concerns the intragenomic homogeneity of dinucleotide relative abundance, and intersequence distance calculations are beyond its scope.
Our use of nonoverlapping windows is consistent with the methodology employed in prior studies using 50 kb and longer windows. (Overlapping windows with a high percentage overlap would be required to localize and sort out the significance of nonconforming segments that could possibly reflect a foreign signature.) Overlap would introduce statistical dependence between successive windows and such dependence could only complicate the analysis of variance within sequences. Window size was varied by factors of approximately two, the specific values being 50 kb, 25 kb, 10 kb, 5 kb, 2 kb, 1 kb, 500 b, 250 b, and 125 b.
Results and discussion
Increasing trend in mean scaled deltadistance
Table 2 shows mean scaled deltadistance by sequence (Abbr) and window size (n). Scaled deltadistance, defined above as Î´*âˆšn, is a statistical invariant for any benchmark sequence generated by a Markov chain exhibiting the same signature as the given sequence. Except for the Borrelia burgdorferi sequence, the mean scaled deltadistance is increasing in window size for every sequence. The lone exception is a drop, for B. burgdorferi, from 3.580 to 3.553, as window size increases from 10 kb to 25 kb.
The essentially increasing trend in Î´*âˆšn has an obvious implication for the scalability of standard binning levels used in classifying intragenomic deltadistance [8]. These levels cannot be rescaled by the reciprocal square root of window size without admitting that the profiles seen through smaller windows are statistically closer to the global signature. The profiles seen through larger windows obviously tend toward the signature but local fluctuations tend to zero slower than 1/âˆšn (i.e., the convergence rate is "subMarkov").
Quasistable hierarchy of mean scaled deltadistances
The rows of Table 2 are grouped as archaea (Aful, Mjan), protist (Pfa2, Pfa3), yeast (Sc11, Sc15), plant (Ath4), human (Hs22), or bacteria (all 14 others). Rows are averaged in groups and the results are plotted in Figure 1. Granted that the curves for archaea and bacteria are almost indistinguishable, we see that ranking mean scaled delta distances by groups produces about the same result for every window size. The one clear exception is the crossing of the curves for plant (Arabidopsis thaliana, chromosome IV) and human (chromosome XXII) between n = 2 kb and n = 5 kb. Thus the view through 50 kb windows shows that the imprint of the global signature is weaker in the latter; but 2 kb and smaller windows show relatively clearer imprints in the human chromosome.
Normality and homoscedasticity of log delta
The increasing trend in all the curves of Figure 1 suggests that local deviations from the global signature are better described by a multiplicative error process (instead of additive). Therefore we examine the (natural) logarithms of the deltadistances in windows and find that they are close to normally distributed with nearly constant variance over the range of window size. Mean values of log(Î´*) are shown in Table 3 with corresponding standard deviations in Table 4.
Each pair of corresponding cells in Tables 3 and 4 defines a normal distribution with mean a variance determined by a sample of observed deltadistances. Each sample is subjected to the KolmogorovSmirnov test of fit to the corresponding normal distribution [22]. (We generate a random sample of the same size as the observed sample from a normal distribution with the same parameters. The KolmogorovSmirnov twosample test is then invoked. The null hypothesis is that the observed and simulated samples share a common underlying distribution.) Pvalues from the test, which measure support for the null hypothesis of normality, are shown in Table 5. The bottom row of Table 5 counts (by window size) the number of times the test accepts normality of log(Î´*) at the 5% significance level.
Deltadistances tend to fit the lognormal distribution better for window sizes 1 kb and larger. Rejection of lognormality for the smallest window sizes is to some extent a consequence of increasing sample size. The apparent goodnessoffit actually gets better as window size decreases as illustrated in Figure 2 for three indicated window sizes on the Haemophilus influenzae sequence. (The sample CDF approaches a smooth curve as window size declines and the number of windows increases. The heavy plotted points are from lognormal distribution functions with means and variances of the corresponding samples.)
SubMarkov convergence rate of mean log delta
Since Î´*âˆšn is essentially increasing in n, we infer that the locally computed dinucleotide relative abundance converges to the global signature more slowly than the reciprocal square root of window size. The lognormality of the intragenomic delta distance suggests a loglinear model of the form log(Î´*) = Î±  Î²log(n) for intercept (Î±) and slope (Î²) coefficients that can be estimated for each sequence (or group) by simple linear regression. Table 3, columns 10 and 11, give intercept and slope coefficients of least squares fits to the scatter plots of mean log(Î´*) versus log(n). Column 12 gives the mean absolute residual (MAR) to indicate the closeness of the straight line fit. All estimated slope coefficients are greater than the benchmark value 1/2 implied by the Markovian analogy. Figure 3 shows least squares fits to loglog scatter plots obtained by grouping rows of Table 3 (in the same as way described above for going from Table 2 to Figure 1).
Weakly consistent patterns of intragenomic variability in deltadistance
Residual accumulated deltadistance (RADD) is defined as window size times the cumulative sum of terms Î´*(t)  mean(Î´*), t = 1, 2, ..., T where t is the position index (counted in windows from the 5' end) and T is the total number of windows. Here mean(Î´*) is just the average of unscaled deltadistances in windows of size n. Since Î´*(1) +...+ Î´*(T) = T mean(Î´*), a plot of RADD versus t always returns to zero. The RADD plot will superficially resemble random walks obtained by integrating counts of purine minus pyrimidine bases. Such "walking plots" are useful in depicting long range compositional biases concomitant to replication [7] and their selfsimilarity with respect to rescaling has been said to instantiate the "fractal geometry of nature [23, 24]."
Here we use the RADD plot to examine how well deltadistances in smaller windows are smoothed by corresponding larger windows. For example, take a single 50 kb window on the M. genitalium sequence, compute its deltadistance from the global signature, then divide it into 50 contiguous 1 kb windows and repeat the calculation for each piece. Do the 50 deltadistances exhibit random fluctuations about the level indicated by the single larger window? If so, the patterns would be judged consistent, and that result would typify a stationary process. Such consistency is only weakly evidenced In Figure 4, however, as larger scales give rise to emergent trends. The same effect is seen in Figure 5 where the 125 b and 1 kb traces track each other closely but diverge from the other traces in the central region.
Conclusions
The imprint of the global signature is locally pervasive on all scales considered in the sequences that were scanned. No lower bound can yet be placed on the local scale on which the global signature is reflected. The intergenomic hierarchy of mean intragenomic deltadistances is essentially preserved across the range of window sizes. Intragenomic deltadistance is approximately lognormally distributed (in windows down to 1 kb) and the variance of logdelta is fairly uniform across the set of sequences. Deltadistance tends to zero with increasing window size but the rate of this convergence is significantly slower than for simple random processes.
Methods
The sequences listed in Table 1 were downloaded from the (National Center for Biotechnology Information) GenBank [25] in FASTA format and saved as plain text files. The last two columns of Table 1 provide GenBank accession numbers and approximate dates of the revisions used in this analysis. The sequences, as text files, were processed by routines written in SPlus, Version 4.5. The texts were read in blocks (contigs) of 50 kb (when computing deltadistances in 50 kb and 25 kb windows) or 20 kb (for smaller windows). Counts of overlapping base steps in windows were computed with the last base of one window serving as the first base of the next. Thus the number of base steps in a window is equal to the number of bases (not one less). However, one base step was uncounted at the start of every block. Incomplete windows in the last block were discarded and blocks past endoffile were ignored. Sometimes complete blocks near endoffile were omitted. With the exception of Arabidopsis thaliana, chromosome IV, the total sequence lengths in kilobases are listed in the fourth column (kb) of Table 1. All corresponding sample lengths are at least 96% of total length, and most are close to 99%, when texts were read in 20 kb blocks. For A. thaliana, however, the global signature was computed from a scan of the 99% complete sequence; but local deltadistances from the global signature were computed only for the first 50% of the complete sequence due to computing difficulties. The 8750 kb length for this sequence in Table 1 is therefore about half the total bases in the chromosome.
Abbreviations
 CDF:

Cumulative Distribution Function
 MAR:

Mean Absolute Residual
 RADD:

Residual Accumulated DeltaDistance, Sequence abbreviations are as shown in Table 1.
References
Karlin S, Brendel V: Patchiness and correlations in DNA sequences. Science. 1993, 259: 667679.
Braun JV, MÃ¼ller HG: Statistical methods for DNA sequence segmentation. Statistical Science. 1998, 13: 142162. 10.1214/ss/1028905933. [http://projecteuclid.org/Dienst/UI/1.0/Display/euclid.ss/1028905933?abstract]
MrÃ¡zek J, Karlin S: Strand compositional asymmetry in bacterial and large viral genomes. Proc Natl Acad Sci USA. 1998, 95: 37203725. 10.1073/pnas.95.7.3720.
Karlin S, MrÃ¡zek J, Campbell AM: Compositional biases of bacterial genomes and evolutionary implications. J Bacteriology. 1997, 179: 38993913.
Russell GJ, SubakSharpe JH: Similarity of the general designs of protochordates and invertebrates. Nature. 1977, 266: 533535.
Karlin S, Burge C: Dinucleotide relative abundance extremes: a genomic signature. Trends in Genetics. 1995, 11: 283290. 10.1016/S01689525(00)890769.
Freeman JM, Plasterer TN, Smith TF, Mohr SC: Patterns of genome organization in bacteria. Science. 1996, 279: 18271829. 10.1126/science.279.5358.1827a.
Karlin S, Landunga I, Blaisdell BE: Heterogeneity of genomes: measures and values. Proc Natl Acad Sci USA. 1994, 91: 1283712841.
Campbell A, MrÃ¡zek J, Karlin S: Genome signature comparisons among prokaryote, plasmid, and mitochondiral DNA. Proc Natl Acad Sci USA. 1999, 96: 91849189. 10.1073/pnas.96.16.9184.
Frick DN, Richardson CC: DNA primases. Annu Rev Biochem. 2001, 70: 3980. 10.1146/annurev.biochem.70.1.39.
Arnold J, Cuticchia AJ, Newsome DA, Jennings WW, Ivarie R: Monothrough hexanucleotide composition of the sense strand of yeast DNA: a Markov chain analysis. Nucleic Acids Res. 1988, 18: 71457158.
Karlin S, Landunga I: Comparisons of eukaryotic genome sequences. Proc Natl Acad Sci USA. 1994, 91: 1283212836.
Ochman H, Lawrence JG, Groisman EA: Lateral gene transfer and the nature of bacterial innovation. Nature. 2000, 405: 299304. 10.1038/35012500.
Blaisdell BE: Effectiveness of measures requiring and not requiring prior sequence alignment for estimating the dissimilarity of natural sequences. J Molec Evol. 1989, 29: 526537.
Pevzner PA: Statistical distance between texts and filtration methods in sequence comparison. ABIOS. 1992, 8: 121127.
Gibbs AJ, Dale MB, Kinns HR, MacKenzie HG: The transition matrix method for comparing sequences. Systematic Zoology. 1971, 20: 417425.
Cardon LR, Burge C, Cayton DA, Karlin S: Pervasive CpG suppression in animal and mitochondrial genomes. Proc Natl Acad Sci USA. 1994, 91: 37993803.
Billingsley P: Statistical methods in Markov chains. Ann Math Stat. 1961, 12: 488497.
Avery PJ, Henderson DA: Fitting Markov chain models to discrete state series such as DNA sequences. Applied Statistics. 1999, 48: 5361. 10.1111/14679876.00139.
Pevzner PA: Nucleotide sequences versus Markov models. Computers Chem. 1992, 16: 103106. 10.1016/00978485(92)80036Y.
Robin S, Daudin JJ: Exact distribution of the distances between any occurrences of a set of words. Ann Inst Statist Math. 2001, 4: 895905. 10.1023/A:1014633825822.
Daniel WW: Applied Nonparametric Statistics,. Boston, PWSKent Pub Co. 1990, 2
Peng CK, Buldyrev SV, Goldberger AL, Havlin S, Mantegna RN, Simon M, Stanley HE: Finitesize effects on longrange correlations: implications for analyzing DNA sequences. Phys Rev E. 1993, 47: 37303733. 10.1103/PhysRevE.47.3730.
Peng CK, Buldyrev SV, Goldberger AL, Havlin S, Mantegna RN, Simon M, Stanley HE: Statistical properties of DNA sequences. Physica A. 1995, 221: 180192.
Benson DA, I KarschMizrachi, Lipman DJ, Ostell J, Rapp BA, Wheeler DL: GenBank. Nucleic Acids Res. 2000, 28: 1518. 10.1093/nar/28.1.15.
Author information
Authors and Affiliations
Corresponding author
Additional information
Authors' contributions
RHB conceived the study and carried out the computations. RWJ participated in experimental design and data analysis. Both authors read 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: Verbatim copying and redistribution of this article are permitted in all media for any purpose, provided this notice is preserved along with the article's original URL.
About this article
Cite this article
Jernigan, R.W., Baran, R.H. Pervasive properties of the genomic signature. BMC Genomics 3, 23 (2002). https://doi.org/10.1186/14712164323
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/14712164323