Skip to main content

AT excursion: a new approach to predict replication origins in viral genomes by locating AT-rich regions

Abstract

Background

Replication origins are considered important sites for understanding the molecular mechanisms involved in DNA replication. Many computational methods have been developed for predicting their locations in archaeal, bacterial and eukaryotic genomes. However, a prediction method designed for a particular kind of genomes might not work well for another. In this paper, we propose the AT excursion method, which is a score-based approach, to quantify local AT abundance in genomic sequences and use the identified high scoring segments for predicting replication origins. This method has the advantages of requiring no preset window size and having rigorous criteria to evaluate statistical significance of high scoring segments.

Results

We have evaluated the AT excursion method by checking its predictions against known replication origins in herpesviruses and comparing its performance with an existing base weighted score method (BWS1). Out of 43 known origins, 39 are predicted by either one or the other method and 26 origins are predicted by both. The excursion method identifies six origins not predicted by BWS1, showing that the AT excursion method is a valuable complement to BWS1. We have also applied the AT excursion method to two other families of double stranded DNA viruses, the poxviruses and iridoviruses, of which very few replication origins are documented in the public domain. The prediction results are made available as supplementary materials at [1]. Preliminary investigation shows that the proposed method works well on some larger genomes too.

Conclusion

The AT excursion method will be a useful computational tool for identifying replication origins in a variety of genomic sequences.

Background

Recent advances in biotechnology have rendered sequencing a complete genome routine. With the increasing availability of DNA sequences, computational methods to predict likely locations of important functional sites before experimental search are highly valuable because the computational predictions can often help design finely tuned experiments to find these functional sites in shorter time with less labor and fewer resources. Replication origins, which are places on the DNA molecules where replication processes are initiated, are considered important sites for understanding the molecular mechanisms involved in DNA replication. For some viruses with double stranded DNA (dsDNA) genomes in particular, detailed knowledge of their replication processes have had significant impact in developing effective strategies to control the growth and spread of viruses (see, for example, [2]).

A number of computational methods have been developed for predicting replication origins in bacterial, archaeal, and eukaryotic genomes. All these algorithms exploit certain characteristic sequence features found around the replication origins. For example, Lobry [3] employs the GC skew plot to predict replication origins and terminus in bacterial genomes. The skew (G-C)/(G+C), where G and C respectively stand for the percentages of guanine and cytosine bases in a sliding window, switches polarity in the vicinity of the replication origin and terminus, with the leading strand manifesting a positive skew. Salzberg et al. [4] predict the replication origins for a number of bacterial and archaeal genomes by identifying some 7-mers and/or 8-mers whose orientation is preferentially skewed around the replication origins. Zhang and Zhang [5] use the Z-curve method successfully to identify several replication origins in bacterial and archaeal genomes. The Z-curve of any given DNA sequence is a three-dimensional curve which uniquely represents the sequence so that unusual sequence compositional features, such as those around a replication origin, can sometimes be visually recognized. Mackiewicz et al. [6] propose three methods, based on DNA asymmetry, the distribution of DnaA boxes and dnaA gene location, were applied to identify the putative replication origins in 112 bacterial chromosomes. They find that DNA asymmetry is the most universal method of putative oriC identification and better prediction can be achieved when the method is applied together with others.

For eukaryotic DNA, Breier et al. [7] develop the Oriscan algorithm to predict replication origins in the S. cerevisiae genome by searching for sequences similar to a training set of 26 known yeast origins pinpointed by site-directed mutagenesis. Oriscan uses both the origin recognition complex binding site and its flanking regions to identify candidates, and then ranks potential origins by their likelihood of activity. More recently, wavelet based multi-scale analysis of DNA strand asymmetries have also been developed [8, 9] for detecting mammalian DNA replication origins.

It is important to note that a prediction method designed for one kind of genomes may not necessarily work well on others because the differences in DNA replication mechanisms in different organisms naturally lead to differences in sequence features around their replication origins. One would not expect that the prediction methods designed for bacterial, archaeal, and eukaryotic genomes can be applied directly to viral genomes and produce accurate results. Indeed, when we attempted to use the above algorithms on some herpesviruses genomes with known replication origins like those listed in Table 5 of [10], a variety of difficulties were encountered. For instance, no clear cut switches of polarity were observed in the GC skew plot. No definitive peaks can be visually identified from the Z-curves as potential replication origins of the viruses. When we mined for DnaA boxes [6] in the herpesviruses, just one cluster of DnaA boxes was observed, but it is not near to any known replication origins. Information about origin recognition complex binding sites for herpesvirus genomes, needed for applying Oriscan, are not readily available. While the method based on oligomers skew [4] is designed to work for genomes with single replication origins, the herpesviruses and many other dsDNA viruses contain multiple replication origins in their genomes.

Computational prediction of replication origins, based on the observation of a high concentration of palindromes around the origins, for dsDNA viral genomes was first attempted by Masse et al. [11] on the human cytomegalovirus. Leung et al. [10] formalize the procedure by laying down the mathematical foundation to justify the use of scan statistics for identifying statistically significant palindrome clusters. The location of such palindrome clusters are then taken to be the likely locations of replication origins in herpesviruses. Viewing the scan statistics approach as equivalent to counting the palindromes in sliding windows, Chew et al. [12] offer two more refined schemes of quantifying palindrome concentration to improve the sensitivity of the prediction. One of these schemes, namely the base weighted scheme (BWS1), which scores each palindrome according to how rarely it is expected to occur in a nucleotide sequence generated randomly as a first order Markov chain, is found to be the most sensitive for the herpesviruses.

Because of the lack of strong family-wide sequence similarities around the origins, the above prediction methods designed for relatively large and complex dsDNA viruses like the herpesviruses with over 100,000 base pairs in the genomes are based on various sequence statistics rather than the actual nucleotide sequences around replication origins.

Herpesviruses utilize two different types of replication origins during lytic and latent infections. For each type of origins, the count and locations in the genome vary from one kind of herpesvirus to another. Most herpesviruses have one to two copies of latent and lytic origins. It has been documented in various studies (e.g. [11, 13, 14]) that the nucleotide sequences around the replication origins are specific to the individual viruses. Yet the presence of clusters of direct or inverted repetitive sequences, including palindromes, is quite common in both types of origins in many members of the herpesvirus family (see [12] and references therein).

Lin et al [15] have observed that in some herpesvirus genomes, the nucleotide sequences around replication origins are richer in A and T bases. This is not surprising because DNA replication typically requires the binding of an assembly of enzymes (e.g., helicases) to locally unwind the DNA helical structure, and pull apart the two complementary strands (see Chapter 1 in [16, 17]). Higher AT content around the origins makes the two complementary DNA strands bond less strongly to each other. This facilitates the two strands to be pulled apart and initiate the replication process. Indeed, Segurado et al. [18] have used a sliding window approach to find "islands" within the Schizosaccharomyces pombe genome that have high AT content. They measure base composition using sliding windows of different sizes and find that AT content of windows in regions containing replication origins are significantly higher than those that do not.

Chew et al. [12] have also reported using sliding windows of AT percentages on herpesviruses. Using windows with top AT percentages they are able to predict 65% of replication origins in their dataset. Moreover, this method has successfully identified four origins not predicted by BWS1, suggesting that the AT percentages may be a useful sequence feature to be incorporated into the set of replication origin prediction tools for dsDNA viruses. This motivates us to seek a means to better quantify the AT content variation in genome sequences. We find that the general score based excursion approach first proposed by Karlin and Altschul in [19] fits our purpose very well when it is applied appropriately to quantify local AT abundance. The excursion approach has the advantages of not requiring a preset sliding window size and having rigorous criteria to evaluate statistical significance of high scoring segments [20–22].

There are three main objectives in this paper. First, we shall develop the AT excursion method as a possible alternative to existing approaches for replication origin prediction in DNA sequences. Second, we shall assess the performance of AT excursion in comparison with the prediction results of BWS1 on a data set of currently known origins of the herpesviruses. The herpes family is chosen as it is one of the bigger families of viruses with known replication origins so that the performance of our prediction method can be assessed. Our results demonstrate that the AT excursion method not only can compare with but can also complement the BWS1 predictions very well. Having established that AT excursion method is a credible prediction tool, our third objective is to use it for predicting likely replication origin locations for two other families of dsDNA viruses, namely the poxviruses and iridoviruses of which very few replication origins are documented in the public domain. To demonstrate the generality of the AT excursion approach, we also apply it to several larger genomes.

Methods

We adopt the score-based excursion approach [19] to identify segments of a genome having high AT concentration. This, in turn, forms the basis of our proposed method to predict replication origins for the herpesviruses. Table 1 presents the viruses to be analyzed. The data set comprises all complete genome sequences of the herpesvirus family downloaded from GenBank at the NCBI web site in March 2006. For each virus, we list its abbreviation, accession number, sequence length, and AT percentages.

Table 1 The list of herpesviruses to be analyzed.

Score-based sequence analysis

Score-based sequence analysis is a powerful and yet flexible tool to identify segments of a biological (DNA, RNA or amino acids) sequence containing high concentration of residues of interest according to the users' objectives. One assigns high positive scores to residues of interest, high negative scores to contrasting residues and low or zero scores for the rest. Using various score schemes, Karlin and his collaborators applied this approach with success to gene finding, identification of transmembrane protein segments, and DNA-binding domains. For details and other applications, see, for example, [20–22] and the references therein.

Our interest in this paper is to identify segments of genomic sequences with high AT content. Towards this end, we label bases C or G as "strongly bonding" base S; and bases A or T as "weakly bonding" base W. Under this label, S bases (i.e., C or G) are given a score of s and W bases (i.e., A or T) a score of w. The scores s and w will be specified below. We next model the genomic sequence as a realization of a sequence of independent and identically distributed random variables, X1, X2, ..., X n (where n is the genome length), taking values in {s, w}. If the i th base is labeled as W, X i is given a score w otherwise X i = s. We let p := P(X i = s) and P (X i = w) = 1 - p (denoted by q). The parameter p is naturally estimated by the CG percentage in the genome. An additional constraint needed to be imposed on the choice of s and w is that the expected score per base μ = ps + qw has to be negative. This condition prevents favoring long segments to be high scoring segments. A moment's reflection shows that we can always standardize one of the scores to be 1. Here we let w = 1 and choose s to be a negative integer (integer-value choice due to a technical reason as pointed out after equation (3)) so that the expected score per base, μ = ps + qw is close to the value of -0.5 (where we adopt Karlin's choice of expected value as in [21]). In other words, w := 1 and

s : = ⌊ μ − q w p ⌋ , MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaacqWGZbWCcqGG6aGocqGH9aqpdaGbdaqaamaalaaabaacciGae8hVd0MaeyOeI0IaemyCaeNaem4DaChabaGaemiCaahaaaGaayj84laawUp+aiabcYcaSaaa@3D1A@
(1)

where μ = -0.5 and ⌊·⌋ denotes the integer floor function.

Excursions and their values

We next compute the cumulative scores and seek to identify segments of the genome that have significantly high scores. As we are only interested in segments with positive additive scores, we reset our cumulative scores to zero whenever it becomes negative.

The excursion scores E i 's are defined recursively as

E0 = 0, E i = max{Ei-1+ X i , 0}, for 1 ≤ i ≤ n.

Using this recursive definition, we are able to construct "excursions" for each of the genomes. An excursion starts at a point i where E i is zero and ends at j > i where E j is the very next zero. The score then stays at zero until it first becomes positive again for the start of the next excursion. The value of an excursion is defined to be the peak score during the course of that particular excursion.

Distribution of the Maximal Aggregate Score

For each value of x, the maximal aggregate score

M n = max 1 ≤ k ≤ n E k MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaacqWGnbqtdaWgaaWcbaGaemOBa4gabeaakiabg2da9maaxababaGagiyBa0MaeiyyaeMaeiiEaGhaleaacqaIXaqmcqGHKjYOcqWGRbWAcqGHKjYOcqWGUbGBaeqaaOGaemyrau0aaSbaaSqaaiabdUgaRbqabaaaaa@3E95@

satisfies

P ( M n > ln n λ ∗ + x ) ≈ 1 − exp { − K ∗ e − λ * x } , MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaacqWGqbaudaqadaqaaiabd2eannaaBaaaleaacqWGUbGBaeqaaOGaeyOpa4ZaaSaaaeaacyGGSbaBcqGGUbGBcqWGUbGBaeaaiiGacqWF7oaBdGaGKYbaaSqajaiPbGaGKIGaaiad0HSFxiIkaaaaaOGaey4kaSIaemiEaGhacaGLOaGaayzkaaGaeyisISRaeGymaeJaeyOeI0IagiyzauMaeiiEaGNaeiiCaaNaei4EaSNaeyOeI0Iaem4saS0aaWbaaSqabeaacWaDK63fIOcaaOGaemyzau2aaWbaaSqabeaacqGHsislcqWF7oaBdaahaaadbeqaaiadaYUGQaGkaaWccqWG4baEaaGccqGG9bqFcqGGSaalaaa@5979@
(2)

where λ* is the unique positive solution to the equation E ( e λ X 1 ) MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaacqWGfbqrdaqadaqaaiabdwgaLnaaCaaaleqabaacciGae83UdWMaemiwaG1aaSbaaWqaaiabigdaXaqabaaaaaGccaGLOaGaayzkaaaaaa@34E3@ = peλs+ qeλw= 1 and K* is a parameter given by an explicit series expansion (See [23]).

When X is an integer-valued variable of span δ, we have a simpler expression for K* ([23]):

exp { − K + e − λ x } ≤ lim inf n → ∞ P ( M n − ln n λ * < x ) ≤ lim sup n → ∞ P ( M n − ln n λ * < x ) ≤ exp { − K − e − λ * x } , MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakqaaeeqaaiGbcwgaLjabcIha4jabcchaWjabcUha7jabgkHiTiabdUealnaaBaaaleaacqGHRaWkaeqaaOGaemyzau2aaWbaaSqabeaacqGHsisliiGacqWF7oaBcWaDez4fIOIaemiEaGhaaOGaeiyFa0NaeyizIm6aaCbeaeaacyGGSbaBcqGGPbqAcqGGTbqBcyGGPbqAcqGGUbGBcqGGMbGzaSqaaiabd6gaUjabgkziUkabg6HiLcqabaGccqWGqbaudaqadaqaaiabd2eannaaBaaaleaacqWGUbGBaeqaaOGaeyOeI0YaaSaaaeaacyGGSbaBcqGGUbGBcqWGUbGBaeaacqWF7oaBcWal0jOkaOcaaiabgYda8iabdIha4bGaayjkaiaawMcaaaqaaiabgsMiJoaaxababaGagiiBaWMaeiyAaKMaeiyBa0Magi4CamNaeiyDauNaeiiCaahaleaacqWGUbGBcqGHsgIRcqGHEisPaeqaaOGaemiuaa1aaeWaaeaacqWGnbqtdaWgaaWcbaGaemOBa4gabeaakiabgkHiTmaalaaabaGagiiBaWMaeiOBa4MaemOBa4gabaGae83UdWMamGdGcQcaQaaacqGH8aapcqWG4baEaiaawIcacaGLPaaaaeaacqGHKjYOcyGGLbqzcqGG4baEcqGGWbaCcqGG7bWEcqGHsislcqWGlbWsdaWgaaWcbaGaeyOeI0cabeaakiabdwgaLnaaCaaaleqabaGaeyOeI0Iae83UdWMamqhAcQcaQiabdIha4baakiabc2ha9jabcYcaSaaaaa@9178@

where

K − = λ * δ e λ * δ − 1 K * , K + = λ * δ 1 − e − λ * δ K * . MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaacqWGlbWsdaWgaaWcbaGaeyOeI0cabeaakiabg2da9maalaaabaacciGae83UdW2aaWbaaSqabeaacWaGulOkaOcaaOGae8hTdqgabaGaemyzau2aaWbaaSqabeaacqWF7oaBdaahaaadbeqaaiadacRGQaGkaaWccqWF0oazaaGccqGHsislcqaIXaqmaaGaem4saS0aaWbaaSqabeaacWaGilOkaOcaaOGaeiilaWIaeeiiaaIaem4saS0aaSbaaSqaaiabgUcaRaqabaGccqGH9aqpdaWcaaqaaiab=T7aSnacmskhaaWcbKaJKgacmsQamai1cQcaQaaakiab=r7aKbqaaiabigdaXiabgkHiTiabdwgaLnaaCaaaleqabaGaeyOeI0Iae83UdW2aiGjGCaaameqcycyaiGjGcWaGWkOkaOcaaSGae8hTdqgaaaaakiabdUealnaaCaaaleqabaGamaiYcQcaQaaakiabc6caUaaa@6263@
(3)

For the simple score scheme with values {-m, ..., -1, 0, 1} occurring with probabilities {p-m, ..., p-1, p0, p1} we have,

K- = (e-λ*- e-2λ*) E (Xeλ*X).

We can set the left hand side of Equation (2) to some predetermined significance level, say P = 0.05 or 0.01, and solve for x. A segment with score exceeding M P = ln n λ * + x MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaacqWGnbqtdaWgaaWcbaGaemiuaafabeaakiabg2da9maalaaabaGagiiBaWMaeiOBa4MaemOBa4gabaacciGae83UdW2aaWbaaSqabeaacWaAulOkaOcaaaaakiabgUcaRiabdIha4baa@3ADC@ is then said to be significant at the 100P% level.

In this paper, we use K- in place of K* in Equation (2) for a "conservative" estimate of the probability and K+ for a "generous" one.

We use Equation (2) with P = 0.05 and P = 0.01 to get M0.05 and M0.01 respectively. If the value of an excursion exceeds the critical value M0.05 (or M0.01), then the segment from the beginning of the excursion up to the base where the peak value is realized is said to be a high-scoring segment (HSS) significant at the 5% (or 1%) level.

HSS Selection

For each of the genomic sequences listed in Table 1, we obtain a set of HSS, significant at the 5% (or 1%) level. In each set of HSS, it is common to find several of them located close to one another. We thus apply a filtering procedure so that, if this happens, we shall only select one of several neighboring excursions as a representative for that part of the genome. In fact, we first sort all the HSS according to their aggregate scores. Starting with the one with the highest value, say segment A, we 'discard' neighboring HSS that are within 2 map units of it. After that, we pick among the rest (not including segment A and the discarded HSS), the HSS with the next highest value, say segment B, and repeat the process. Only the representative segments A, B, and so forth, will be used in replication origin prediction.

Results and Discussion

HSS Tables and Excursion Plots

Table 2 lists the HSS for each herpesvirus in Table 1. We have also tried locating high-scoring segments by running the excursions from the 3' end to 5' end of the genome. The results obtained are not much different from the "vanilla" version (i.e., from 5' to 3').

Table 2 Herpesviruses : HSS at 5% level using the conservative bound.

For visualizing the locations of the selected HSS relative to the entire genome, the excursion plot is a convenient tool. The excursion plot of the Human Herpesvirus 3 (vzv) is presented in Figure 1, where the AT excursion values are plotted against the bases along the genome. The general appearance of Figure 1 is typical of the excursion plots for all the herpesviruses analyzed. In the case of vzv, three peaks with excursion values exceeding the 5% significance level are observed. Two of these peaks are close to the centers of the only two known replication origins of vzv (see Table 3).

Figure 1
figure 1

The Excursion Plot of the vzv virus. The horizontal line corresponds to the 5% significant level. The two triangles denote the locations of known replication origins of the vzv.

Table 3 Prediction results at 5% level using the conservative bound.

Prediction Performance

The high-scoring segments are checked against known replication origins in herpesviruses to evaluate their performance as a prediction tool. Table 3 lists all the known replication origins for the herpesviruses in Table 1. These origins are reported either in published literature or GenBank annotations. For each replication origin, we list the HSS (at 5% level) closest to it. For this table we had used the "conservative" estimate for the value of K* (See Equations (2) and (3)). When the peak of an HSS is less than 2 map units (one map unit is one percent of the genome length) away from the center of a replication origin, we say that our method has correctly predicted that particular replication origin. From Table 3, we see that of the 43 replication origins known, compiled from literature or annotations, 32 of them are close to HSS that have been identified.

We had also tried using the "generous" estimate for K* at the 5% and 1% level of significance. Table 4 gives a summary of the performance of our prediction scheme when those bounds were used. The first two columns of the table gives the sensitivity level and positive prediction value of our scheme. Sensitivity refers to the percentage of replication origins predicted by our method, and PPV (positive predictive value) the proportion of HSS that correctly predict replication origins. APD (average predictive distance), given in map units (± one standard deviation), shows the average of the distances (in map units) between the center of each replication origin and the HSS that predicts it. Note that the APD values say that on average, when a prediction by an HSS is successful, the replication origin is about 0.35 map units away from it. We have also done some simple analysis of the location of the center of each replication origin with respect to the HSS closest to it. We count the number of times the center of replication origin falls within the left, right or center of the HSS. The columns %L, %R, and %C in Table 4 give these proportions. Our results show that the origin falls within the center of the HSS half the time.

Table 4 Prediction Performance Summary.

Comparison with Other Approaches

How does the AT excursion method compare with the sliding window approach using palindrome based scoring schemes previously presented in [12]? Since the BWS1 scheme has been shown to perform best among the various palindrome based schemes, we have examined the numbers of replication origins correctly predicted by AT excursion and by BWS1. The results are summarized in Figure 2.

Figure 2
figure 2

Predictions of AT excursion and BWS 1 . In this figure, the set A consists of origin replications predicted by the AT excursion method and B consists of those predicted by the BWS1 method. A ⋂ BC= {cehv71, cehv72, ehv41, hsv21, hsv22, hsv23}, AC⋂ B = {cehv162, cehv163, ebv1, ebv3, hhv6, hhv6b, rcmv}, (A ⋃ B)C= {bohv4, ehv42, ehv43, hhv7}. The rest of the replication origins (26 of them) are predicted by both methods. Note that for viruses with several known replication origins, such as the hsv2, which has three (see Table 3), we denote the replication origins as hsv21, hsv22, hsv23, etc.

The majority of the 43 known origins in the herpesviruses listed in Table 1 are predicted by both methods and most of the remaining ones are predicted by one method or the other. Only four of the origins fail to be predicted by either method. This suggests that the AT excursion method and the BWS1 scheme complement each other very well.

There are certain advantages in the AT excursion approach over BWS1. First, AT excursion does not require any sequence specific parameters to be prescribed by the user. It is window size free because it does not require any sliding window to measure AT concentration. Moreover, while the palindrome based methods require the specification of a minimal palindrome length before the analysis can be carried out, no such parameter is needed for AT excursion. Second, the AT excursion method is statistically based, as the probabilistic distribution has already been established [20–22]. This allows the statistical significance for HSS be evaluated easily.

We also note that the more elaborate AT excursion approach performs better than the simpler procedure of measuring the percentage of A and T bases on a sliding window in terms of number of correct predictions and the proximity of these predictions to the true origins. Out of the 43 known replication origins for the herpesviruses in Table 1, 32 are correctly predicted by AT excursion but only 28 by AT sliding window plot. Furthermore, the boxplots of the predictive distances (Figure 3) of the AT excursion approach suggests that the predictions given by the AT excursion approach are much closer to known replication origins as compared to those of the AT sliding window plot approach. (In fact, the predictive distances of the AT excursion approach compared to that of the PLS and BWS1 approaches mentioned in [12] are observably shorter. See Figure 3.) This suggests that the excursion values might more correctly capture the essence of A/T abundance variation along genomic sequences.

Figure 3
figure 3

Predictive Distances for PLS, BWS 1 , AT-swp and AT excursion. These boxplots show the predictive distances for PLS, BWS1, AT-swp and AT excursion.

Herpesvirus Replication Origins Alignment and Motif Finding

One might ask whether or not the nucleotide sequences around replication origins in various viruses of the same family share sufficient similarities so that the origins can be identified by sequence alignments and motif finding techniques. We therefore extracted the nucleotide sequences of the known herpesvirus origins according to their documented locations for closer examination. These sequences are available as supplementary materials on the companion website. A multiple alignment using CLUSTAL W [24] and motif searches using MEME and MAST [25, 26] have been conducted for the herpesvirus origin sequences. No significant sequence similarity or common motif pattern across all the origin sequences has been found, agreeing with the findings of [11, 13, 14].

What if we first classify these nucleotide sequences according to some classification schemes, will the members within each class share noticeable sequence similarities? We classified the origins according to (i) the sub-family of the virus (herpesviruses are classified into the alpha, beta, and gamma sub-families by their biological properties [27]), (ii) the type of origin (i.e., whether the origin is a oriL, oriLyt or oriS). We ran MEME and MAST separately on the sequences in each sub-family/type of origins to detect common motif patterns. From the outputs under classification (i), we note that the origins from the alpha sub-family can be further divided into two groups. Each group has a common motif pattern across its members. For the beta and gamma sub-families, no distinct patterns can be found. However, the rcmv and ebv origins contain many repeat patterns. For classification (ii), we find that both the oriL and oriLyt origins contain sequence motifs common to a number of their members. No motif was found for oriS sequences. The results of our motif search are made available in the supplementary materials.

Although our investigations are preliminary, the motifs found in these subsets of herpesvirus genomes may suggest new information that can be incorporated into the replication origin prediction procedures.

Other Families of Viruses

Aside from the herpesviruses, we have also applied the AT excursion method to search for HSS in the poxviruses and iridoviruses. These two viral families are chosen because, like the herpesviruses, they are large, complex dsDNA viruses with no RNA stage. Their genome lengths are also similar in magnitude to those of the herpesviruses.

Poxviruses infect a large variety of animal species that gather in swarms and herds (e.g., mosquitoes, cows). Smallpox is a major disease caused by the variola virus, a member of the poxvirus family. Smallpox was eradicated in 1977 by preventive inoculations with cowpox or vaccinia viruses through the dedicated efforts of the World Health Organization and many individuals. In the recent few years, as the threat of the variola virus being used as a biological weapon is raised, there is growing interest in further studying poxviruses for biodefense purposes [28, 29]. Iridoviruses are found in a variety of fish, amphibians, and reptiles. Some iridoviruses have been associated with serious diseases (e.g., viral erythrocytic necrosis of salmonids), while others have only been found in apparently healthy animals (e.g., goldfish iridovirus). Iridovirus infection is considered a serious concern in modern aquaculture, fish farming, and wildlife conservation [30].

Amongst these two families, only one genome, namely the Chilo iridescent virus, has documented replication origin locations [31]. Our method has correctly predicted one of these locations. Due to the lack of confirmed origin locations, prediction accuracy cannot be tested on these families. Nevertheless, our predictions may assist researchers to investigate these viruses experimentally to identify and confirm the exact locations of replication origins in their genomes. We have, therefore, made our prediction results available at [1].

AT excursion applied to larger genomes

To gauge whether the AT excursion approach can potentially be generalized to predict replication origins for non-viral genomes, we apply it to several archaeal and bacterial genomes which have been previously analyzed. From [4, 5, 32] we are able to compile a list of 15 known or suggested replication origins (11 known, 4 suggested). Using the AT excursion method, we manage to correctly predict 9 of the replication origins (6 known, 3 suggested). Although our studies are preliminary, the results show that the AT excursion method can work reasonably well even on larger genomes.

Conclusion

This paper introduces the AT excursion method to quantify local AT abundance in genomic sequences. The simple and intuitive idea of locating regions with high AT content as potential replication origin sites proves to be effective in identifying several replication origins not previously predicted. This shows that the AT excursion approach is a valuable addition to existing prediction tools. However, we have also observed that quite a number of the statistically significant HSS found by AT excursions are not close to replication origins. Whether these HSS correspond to other important functional sites in the genomic sequences remains an interesting question to be investigated.

The availability of statistical significance criteria and the independence of ad hoc parameters like the minimal palindrome length and sliding window size make the AT excursion method particularly easy to apply to those viral genomes where no replication origin information in similar and related genomes is available. On the other hand, if such information is available, the AT excursion method is not capable of taking advantage of it. To address this issue, machine learning approaches (e.g., neural networks and support vector machines), which better allow us to use knowledge in related genomes, are currently being explored. We anticipate that a combination of score based statistics with machine learning approaches will provide a highly accurate prediction tool set for replication origins.

References

  1. Supplementary Materials[http://www.stat.nus.edu.sg/~stadc/papers/ATexcursion]

  2. Delecluse HJ, Hammerschmidt W: The genetic approach to the Epstein-Barr virus: from basic virology to gene therapy. Mol Pathol 2000, 53(5):270–279. 10.1136/mp.53.5.270

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  3. Lobry JR: Asymmetric substitution patterns in the two DNA strands of bacteria. Mol Biol Evol 1996, 13(5):660–665.

    Article  CAS  PubMed  Google Scholar 

  4. Salzberg SL, Salzberg AJ, Kerlavage AR, Tomb JF: Skewed oligomers and origins of replication. Gene 1998, 217(1–2):57–67. 10.1016/S0378-1119(98)00374-6

    Article  CAS  PubMed  Google Scholar 

  5. Zhang R, Zhang CT: Identification of replication origins in archaeal genomes based on the Z-curve method. Archaea 2005, 1(5):335–346.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  6. Mackiewicz P, Zakrzewska-Czerwinska J, Zawilak A, Dudek MR, Cebrat S: Where does bacterial replication start? Rules for predicting the oriC region. Nucleic Acids Res 2004, 32(13):3781–3791. [http://0-dx-doi-org.brum.beds.ac.uk/10.1093/nar/gkh699] 10.1093/nar/gkh699

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  7. Breier AM, Chatterji S, Cozzarelli NR: Prediction of Saccharomyces cerevisiae replication origins. Genome Biol 2004, 5(4):R22. [http://0-dx-doi-org.brum.beds.ac.uk/10.1186/gb-2004–5-4-r22] 10.1186/gb-2004-5-4-r22

    Article  PubMed Central  PubMed  Google Scholar 

  8. Brodie of Brodie EB, Nicolay S, Touchon M, Audit B, d'Aubenton Carafa Y, Thermes C, Arneodo A: From DNA sequence analysis to modeling replication in the human genome. Phys Rev Lett 2005, 94(24):248103. 10.1103/PhysRevLett.94.248103

    Article  CAS  PubMed  Google Scholar 

  9. Touchon M, Nicolay S, Audit B, Brodie of Brodie EB, d'Aubenton Carafa Y, Arneodo A, Thermes C: Replication-associated strand asymmetries in mammalian genomes: toward detection of replication origins. Proc Natl Acad Sci USA 2005, 102(28):9836–9841. [http://0-dx-doi-org.brum.beds.ac.uk/10.1073/pnas.0500577102] 10.1073/pnas.0500577102

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  10. Leung MY, Choi KP, Xia A, Chen LHY: Nonrandom clusters of palindromes in herpesvirus genomes. J Comput Biol 2005, 12(3):331–354. [http://0-dx-doi-org.brum.beds.ac.uk/10.1089/cmb.2005.12.331] 10.1089/cmb.2005.12.331

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  11. Masse MJ, Karlin S, Schachtel GA, Mocarski ES: Human cytomegalovirus origin of DNA replication (oriLyt) resides within a highly complex repetitive region. Proc Natl Acad Sci USA 1992, 89(12):5246–5250. 10.1073/pnas.89.12.5246

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  12. Chew DSH, Choi KP, Leung MY: Scoring schemes of palindrome clusters for more sensitive prediction of replication origins in herpesviruses. Nucleic Acids Res 2005, 33(15):e134. [http://0-dx-doi-org.brum.beds.ac.uk/10.1093/nar/gni135] 10.1093/nar/gni135

    Article  PubMed Central  PubMed  Google Scholar 

  13. Hamzeh FM, Lietman PS, Gibson W, Hayward GS: Identification of the lytic origin of DNA replication in human cytomegalovirus by a novel approach utilizing ganciclovir-induced chain termination. J Virol 1990, 64(12):6184–6195.

    PubMed Central  CAS  PubMed  Google Scholar 

  14. Dykes C, Chan H, Krenitsky DM, Dewhurst S: Stringent structural and sequence requirements of the human herpesvirus 6B lytic-phase origin of DNA replication. J Gen Virol 1997, 78(Pt 5):1125–1129.

    Article  CAS  PubMed  Google Scholar 

  15. Lin CL, Li H, Wang Y, Zhu FX, Kudchodkar S, Yuan Y: Kaposi's sarcoma-associated herpesvirus lytic origin (ori-Lyt)-dependent DNA replication: identification of the ori-Lyt and association of K8 bZip protein with the origin. J Virol 2003, 77(10):5578–5588. 10.1128/JVI.77.10.5578-5588.2003

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  16. Kornberg A, Baker TA: DNA Replication. 2nd edition. New York, New York: WH Freeman and Company; 1992.

    Google Scholar 

  17. Bramhill D, Kornberg A: A model for initiation at origins of DNA replication. Cell 1988, 54(7):915–918. 10.1016/0092-8674(88)90102-X

    Article  CAS  PubMed  Google Scholar 

  18. Segurado M, de Luis A, Antequera F: Genome-wide distribution of DNA replication origins at A+T-rich islands in Schizosaccharomyces pombe. EMBO Rep 2003, 4(11):1048–1053. [http://0-dx-doi-org.brum.beds.ac.uk/10.1038/sj.embor.embor7400008] 10.1038/sj.embor.7400008

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  19. Karlin S, Altschul SF: Methods for assessing the statistical significance of molecular sequence features by using general scoring schemes. Proc Natl Acad Sci USA 1990, 87(6):2264–2268. 10.1073/pnas.87.6.2264

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  20. Karlin S, Brendel V: Chance and statistical significance in protein and DNA sequence analysis. Science 1992, 257(5066):39–49. 10.1126/science.1621093

    Article  CAS  PubMed  Google Scholar 

  21. Karlin S: Statistical studies of biomolecular sequences: score-based methods. Phil Trans R Soc Lond B 1994, 344(1310):391–402. 10.1098/rstb.1994.0078

    Article  CAS  Google Scholar 

  22. Karlin S: Statistical signals in bioinformatics. Proc Natl Acad Sci USA 2005, 102(38):13355–13362. [http://0-dx-doi-org.brum.beds.ac.uk/10.1073/pnas.0501804102] 10.1073/pnas.0501804102

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  23. Karlin S, Dembo A, Kawabata T: Statistical composition of high-scoring segments from molecular sequences. Ann Statist 1990, 18(2):571–581.

    Article  Google Scholar 

  24. Thompson JD, Higgins DG, Gibson TJ: CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res 1994, 22(22):4673–4680. 10.1093/nar/22.22.4673

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  25. Bailey TL, Elkan C: Fitting a mixture model by expectation maximization to discover motifs in biopolymers. In Proceedings of the Second International Conference on Intelligent Systems for Molecular Biology. Volume 2. Menlo Park, California: AAAI Press; 1994:28–36.

    Google Scholar 

  26. Bailey TL, Gribskov M: Combining evidence using p-values: application to sequence homology searches. Bioinformatics 1998, 14: 48–54. 10.1093/bioinformatics/14.1.48

    Article  CAS  PubMed  Google Scholar 

  27. Cann AJ: Principles of Molecular Virology. 4th edition. San Diego: Elsevier Academic Press; 2005.

    Google Scholar 

  28. Henderson DA: The looming threat of bioterrorism. Science 1999, 283(5406):1279–1282. 10.1126/science.283.5406.1279

    Article  CAS  PubMed  Google Scholar 

  29. Miller SE: Bioterrorism and electron microscopic differentiation of poxviruses from herpesviruses: dos and don'ts. Ultrastruct Pathol 2003, 27(3):133–140. 10.1080/01913120309932

    Article  PubMed  Google Scholar 

  30. Tsai CT, Ting JW, Wu MH, Wu MF, Guo IC, Chang CY: Complete genome sequence of the grouper iridovirus and comparison of genomic organization with those of other iridoviruses. J Virol 2005, 79(4):2010–2023. [http://0-dx-doi-org.brum.beds.ac.uk/10.1128/JVI.79.4.2010–2023.2005] 10.1128/JVI.79.4.2010-2023.2005

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  31. Handermann M, Schnitzler P, Rösen-Wolff A, Raab K, Sonntag KC, Darai G: Identification and mapping of origins of DNA replication within the DNA sequences of the genome of insect iridescent virus type 6. Virus Genes 1992, 6: 19–32. 10.1007/BF01703754

    Article  CAS  PubMed  Google Scholar 

  32. Freeman JM, Plasterer TN, Smith TF, Mohr SC: Patterns of genome organization in bacteria. Science 1998, 279(5358):1827a. [http://www.sciencemag.org] 10.1126/science.279.5358.1827a

    Article  Google Scholar 

Download references

Acknowledgements

We would like to thank the editor and the anonymous reviewers for helpful comments and suggestions which improve this paper. Kwok Pui Choi was supported by National University of Singapore ARF Research grant R-155-000-051-112; and Ming-Ying Leung by NIH grants S06GM08012-35, 5G12RR008124-11 and 3T34GM008048-20S1.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to David SH Chew.

Additional information

Authors' contributions

DC participated in the design of the study and performed the data and statistical analysis. KPC and MYL conceived the study, and participated in its design and coordination. All authors contributed to writing, reading and approving the final manuscript.

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.

Reprints and permissions

About this article

Cite this article

Chew, D.S., Leung, MY. & Choi, K.P. AT excursion: a new approach to predict replication origins in viral genomes by locating AT-rich regions. BMC Bioinformatics 8, 163 (2007). https://0-doi-org.brum.beds.ac.uk/10.1186/1471-2105-8-163

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://0-doi-org.brum.beds.ac.uk/10.1186/1471-2105-8-163

Keywords