Email updates

Keep up to date with the latest news and content from BMC Genomics and BioMed Central.

This article is part of the supplement: The International Conference on Intelligent Biology and Medicine (ICIBM) – Genomics

Open Access Research

Detecting the borders between coding and non-coding DNA regions in prokaryotes based on recursive segmentation and nucleotide doublets statistics

Suping Deng12, Yixiang Shi3, Liyun Yuan34, Yixue Li23 and Guohui Ding2*

Author Affiliations

1 School of Life Sciences and Technology, Tongji University, 1239 Siping Road, Shanghai, 200092, P.R. China

2 Key lab of systems biology, Shanghai Institutes for Biological Sciences, Chinese Academy of Sciences, 320 Yueyang Road, Shanghai, 200031, P.R.China

3 Shanghai Center for Bioinformation Technology, 100 Qinzhou Road, Shanghai, 200235, P. R. China

4 Shanghai Jiaotong University, Shanghai, 200240, P. R. China

For all author emails, please log on.

BMC Genomics 2012, 13(Suppl 8):S19  doi:10.1186/1471-2164-13-S8-S19

The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1471-2164/13/S8/S19


Published:17 December 2012

© 2012 Deng et al.; licensee 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.

Abstract

Background

Detecting the borders between coding and non-coding regions is an essential step in the genome annotation. And information entropy measures are useful for describing the signals in genome sequence. However, the accuracies of previous methods of finding borders based on entropy segmentation method still need to be improved.

Methods

In this study, we first applied a new recursive entropic segmentation method on DNA sequences to get preliminary significant cuts. A 22-symbol alphabet is used to capture the differential composition of nucleotide doublets and stop codon patterns along three phases in both DNA strands. This process requires no prior training datasets.

Results

Comparing with the previous segmentation methods, the experimental results on three bacteria genomes, Rickettsia prowazekii, Borrelia burgdorferi and E.coli, show that our approach improves the accuracy for finding the borders between coding and non-coding regions in DNA sequences.

Conclusions

This paper presents a new segmentation method in prokaryotes based on Jensen-Rényi divergence with a 22-symbol alphabet. For three bacteria genomes, comparing to A12_JR method, our method raised the accuracy of finding the borders between protein coding and non-coding regions in DNA sequences.

Background

The prediction of protein coding regions in DNA sequences is a major goal and a long-lasting topic in molecular biology, especially for the genome projects [1-6]. Lots of methods for finding probable borders are based on strong signals between the coding regions and the non-coding ones [7,8]. Staden [9] used the intersection method to detect the borders between coding and non-coding regions. The information entropy measures for signals are useful for identifying the homogeneous regions and evaluating the genomic complexity [10-12]. The entropy-based segmentation methods can be used to identify the borders between coding and non-coding regions [10,13,14]. The Jensen-Shannon divergence measure has provided an impelling tool in doing this [8,9,15]. Bernaola-Galvan et al. presented an entropic segmentation method to search the borders [12]. The accuracy of their results was higher than those obtained with the intersection method [9,12]. The segmentation method presented by Nicorici et al. [8] was based on the Jensen-Rényi divergence measure in both DNA strands. In 2007, Zhang et al. [16] introduced a segmentation method based on a R14 alphabet and the β-KL divergence. However, its accuracy is not higher than Nicorici's method [8].

In this study, we constructed a 22-symbol alphabet to represent DNA sequences. Based on the entropy theory, we used recursive segmentation to detect the borders between coding and non-coding DNA regions. Comparing to previous methods, it is shown that our accuracy was well improved.

Materials and methods

The data set

Three tested genomes were downloaded from the National Center for Biotechnology Information (https://www.ncbi.nlm.nih.gov/ webcite): Rickettsia prowazekii (GenBank: AJ235269), Borrelia burgdorferi (GenBank: NC_000948 and AE000783) and E.coli (GenBank: NC_009837, NC_008563 and NC_010468).

A22 alphabet

The statistical properties of DNA sequences were commonly used to recognize protein coding regions [11,12,14,17-20]. The statistical properties of doublets of nucleotides (called di-nucleotide for short) in coding regions are also different from those in non-coding regions. This may be used to predict the coding DNA regions. Each di-nucleotide of the DNA sequence is substituted by the symbols from A16 = {AA, AC, AG, AT, CA, CC, CG, CT, GA, GC, GG, GT, TA, TC, TG, TT} (Table 1).

Table 1. Di-nucleotides mapping in 22-symbol alphabet.

The distribution of stop codon patterns (SCPs for short) in DNA coding regions differs from that in the non-coding regions [3,21]. It is well known that the SCPs are strong signals in DNA sequences, so that we can effectively use these signals to detect borders between coding and non-coding DNA regions [22]. The SCPs, TGA, TAG and TAA appear in one given DNA strand, and the three SCPs corresponding to TCA, CTA and TTA appear on the reverse strand. In this way, the SCPs statistics on both DNA strands is the same as the statistics of the six codons TAA, TAG, TGA, TCA, CTA, and TTA on a single DNA strand.

In our study, we introduced a 22-symbol alphabet (called A22 for short) that took into account the non-uniform distribution of di-nucleotides and SCPs in both DNA strands (Table 1 and Table 2). Thus the di-nucleotides and the SCPs are substituted by the symbols from A22 = {AA, AC, AG, AT, CA, CC, CG, CT, GA, GC, GG, GT, TA, TC, TG, TT, S1, S2, S3, S'1, S'2, S'3}(Table 1 and Table 2). The phase of the nucleotide is defined as m = ((n-1)mod 3)+1, where m ∈{1,2,3}, and n is the position of the nucleotide in the DNA sequence. The phase of a SCP is defined in the same way with the exception that n represents the position of the first nucleotide of the given codon. For example, the DNA sequence ACGTAATC is converted using the A22 alphabet as AC, CG, GT, TA, S1, AA, S'2, AT, TC.

Table 2. SCPs mapping in 22-symbol alphabet.

Detecting borders between coding and non-coding DNA regions

In order to partition a DNA sequence, we used the approach proposed by Nicorici et al [8]. and Li [10,13]. A sliding pointer is moving along the sequence. At each position, the pointer divided the sequence into two subsequences and we computed the Jensen-Rényi divergence <a onClick="popup('http://www.biomedcentral.com/1471-2164/13/S8/S19/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2164/13/S8/S19/mathml/M1">View MathML</a>. Then, we found the maximum <a onClick="popup('http://www.biomedcentral.com/1471-2164/13/S8/S19/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2164/13/S8/S19/mathml/M1">View MathML</a> and computed its segmentation strength s (see below). If this segmentation strength s exceeded a given threshold s0, the position was identified as a significant cut (or a probable border) between coding and non-coding DNA regions. The procedure continued recursively for each of the two resulting subsequences created by each cut until none of the cuts had a segmentation strength level exceeding the s0. Then such a sequence was segmented at the segmentation strength level s0.

In this study, the Jensen-Rényi divergence [23,24] is defined as follows:

<a onClick="popup('http://www.biomedcentral.com/1471-2164/13/S8/S19/mathml/M2','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2164/13/S8/S19/mathml/M2">View MathML</a>

(1)

Where Ra, Ra, l and Ra, r are the Jensen-Rényi entropies of the whole, left, and right subsequences, respectively.

Stopping criterion

To decide when the segmentation process has to be stopped, we adopted the method proposed in references [8,25,26] and introduced a segmentation strength, derived empirically, as

<a onClick="popup('http://www.biomedcentral.com/1471-2164/13/S8/S19/mathml/M3','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2164/13/S8/S19/mathml/M3">View MathML</a>

(2)

The recursive segmentation continues as long as s s0 and the segmented sequence has SCPs in all three phases, where s0 can be set by the user. K is a constant, which was set as 16 [8].

Actually, the probable borders (the significant cuts) predicted by the recursive segmentation method is generally not the actual borders but are close to them. Since it is well known that the codons at the real borders between coding and non-coding DNA regions must be one kind of start or stop codons, we could use start or stop codons like nuclear acid pattern around the border as border cut. Then, we filter the segment region less than 20 bp. The procedure for finding borders between coding and non-coding DNA regions can be described by the flow chart (Figure 1).

thumbnailFigure 1. The flow chart of the procedure for finding borders between coding and non-coding DNA regions.

Evaluation

In order to evaluate how well the predicted borders matching the actual borders between coding and non-coding regions, we use the following measure introduced by Bernaola-Galvan et al. [12].

<a onClick="popup('http://www.biomedcentral.com/1471-2164/13/S8/S19/mathml/M4','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2164/13/S8/S19/mathml/M4">View MathML</a>

(3)

Where {bi} is the set of all known borders (called KBs for short) between coding and non-coding regions, and { cj} is the set of all predicted borders (called PBs for short), and NT is the total length of the DNA sequence. The first summation measures the discrepancy between PBs and KBs by adding the distance from each KB to the closest PB, and the second summation performs the same operation, but includes the distance from each PB to the closest KB. Both are required to take into account not only the correctness for the cutting position (CBC would be zero only when the PBs overlap the KBs), but also the difference between the number of PBs and KBs. CBC can be viewed as an average of the error in determining the correct boundaries between coding and non-coding regions, so (1-CBC) is a reasonable measure of the accuracy of the method.

Results and discussion

In Figure 2, we plotted the Jensen-Rényi divergence (α = 0.5 and used in the following experiments as the prediction results have no change when α is adjusted from 0 to 1) with A12 [12] and A22 alphabets along a DNA segment. The DNA segment was randomly chosen from the bacterium genome Borrelia burgdorferi and Rickettsia prowazekii. In Figure 2(a), the analyzed DNA segment was chosen from bacterium Rickettsia prowazekii (AJ235269, 3757-6226 bp). The left part (length 2121 bp) belongs to a coding region and the right part (length 350 bp) belongs to a non-coding region. In Figure 2(b), the analyzed DNA segment was chosen from bacterium Rickettsia prowazekii (AJ235269, 10683-11820 bp). The left part (length 1074 bp) belongs to a coding region and the right one belongs to a non-coding region. From Figure 2, the cuts predicted by A22-JR (the method with A22 alphabet, Jensen-Rényi divergence) are closer to the real borders than those by A12-JR.

thumbnailFigure 2. Jensen-Rényi divergence versus cutting position for a DNA sequence. The DNA sequence contains a coding region followed by a non-coding region. The maximum values for the divergences are circled on the graph. (a) The analyzed DNA segment was chosen from bacterium Rickettsia prowazekii (AJ235269, 3757-6226 bp). (b) The analyzed DNA segment was chosen from bacterium Rickettsia prowazekii (AJ235269, 10683-11820 bp).

We also applied the two methods to whole genome respectively. There are multiple coding and non-coding regions in those sequences. The results are summarized in Table 3. The accuracy of A22_JR is better than that of A12_JR for each DNA sequence (p = 0.0015, Table 3).

Table 3. The maximum accuracy of different methods applied to different data sets.

For visualizing the borders predicted by our proposed method, we plotted the known coding regions in the first 22000 bp of the bacterium genome Borrelia burgdorferi (AE000783) and the unmodified predicted borders from our results (Figure 3).

thumbnailFigure 3. Comparison between the known coding regions and the predicted borders of a DNA sequence. The known coding regions are gray regions with solid lines as borders. The predicted borders (vertical dotted lines) is obtained through recursive segmentation using A22_JR (a = 0.5). The DNA sequence is from bacteria Rickettsia prowazekii and the borders. The coding regions shown downwards are on the opposite DNA strand.

Finally, we described how to choose an appropriate threshold s0 of segmentation strength. After having gotten the cuts and their corresponding segmentation strength, s0 ranged from 0.30 to 1.00 stepping by 0.01. For each s0, the accuracy was computed. From Figure 4, we can find that the accuracy is much higher when s0 is about -0.50. Thus the appropriate threshold s0 of segmentation strength can be set as -0.50.

thumbnailFigure 4. Accuracies of recursive segmentation for different thresholds of segmentation strength. The DNA sequences are from the genomes of Rickettsia prowazekiiand and Borrelia burgdorferi (the first 30000 bp).

In this study, we introduced a new segmentation method for finding the borders between coding and non-coding regions. It is based on the Jensen-Rényi divergence, a 22-symbol alphabet, and a new stopping criterion. Tested on three bacteria genomes, our method improved the accuracies of the borders detection over the previously reported A12-JR segmentation approach. Most of the existing segmentation algorithms [10,12,13] rely heavily on statistical properties of the coding, non-coding or other interested regions in DNA sequences. Moreover, since the gene-finding systems [24,26,27] use biological knowledge regarding functional sites, together with statistics for finding genes, they require extensive training on known datasets. The recursive segmentation needs no prior training. It should be noted that the value of the segmentation strength threshold s0 is generally set as -0.50 for bacterium and may be adjusted accordingly in different species. For a new unknown genomic sequence, the optimal threshold s0 of segmentation strength or significance level can be computed using the genomic sequence of the same or the closest organism.

Conclusions

The borders between coding and non-coding regions are found more efficient and accurate will raise the vital effect for DNA sequences annotation. This paper presents a new segmentation method based on Jensen-Rényi divergence with a 22-symbol alphabet, new stopping criterion for finding the borders between coding and non-coding DNA regions in prokaryotes. For three bacteria genomes, comparing to A12_JR method, our method raised the accuracy of finding the borders between coding and non-coding regions in DNA sequences. The success comes from the utilization of the di-nucleotides and SCPs statistics in all three phases along the DNA sequence, and use of Jensen-Rényi divergence.

Competing interests

The authors declare that they have no competing interests.

Acknowledgements

This research was supported by grants from from National High-Tech R&D Program (863) (2009AA02Z304, 2012AA020404), State key basic research program (973) (2006CB910705, 2010CB529206, 2011CBA00801), Research Program of CAS (KSCX2-YW-R-112, KSCX2-YW-R-190, 2011KIP204), National Natural Science Foundation of China (30900272) and SA-SIBS Scholarship Program.

This article has been published as part of BMC Genomics Volume 13 Supplement 8, 2012: Proceedings of The International Conference on Intelligent Biology and Medicine (ICIBM): Genomics. The full contents of the supplement are available online at http://www.biomedcentral.com/bmcgenomics/supplements/13/S8.

References

  1. Li W: The complexity of DNA.

    Complexity 1997, 3:33-37. Publisher Full Text OpenURL

  2. Zhang CT, W J: Recognition of protein coding genes in the yeast genome at better than 95% accuracy based on Z curve.

    Nucleic Acids Res 2000, 28:2804-2814. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  3. Stanke M, W S: Gene prediction with a hidden Markov model and a new intron submodel.

    Bioinformatics 2003, 19:ii215-ii225. PubMed Abstract | Publisher Full Text OpenURL

  4. Haimovich AD, B B, Ramaswamy R, Welsh WJ: Wavelet analysis of DNA walks.

    J Comput Biol 2006, 13:1289-1298. PubMed Abstract | Publisher Full Text OpenURL

  5. Orlov YL, T R, Abnizova I: Statistical measures of the structure of genomic sequences: entropy, complexity and position information.

    JBioinform Comput Biol 2006, 4:523-526. PubMed Abstract | Publisher Full Text OpenURL

  6. TeBoekhorst R, A I, Nehaniv C: Discriminating coding, non-coding and regulatory regions using rescaled range and detrended fluctuation analysis.

    BioSystems 2008, 91:183-194. PubMed Abstract | Publisher Full Text OpenURL

  7. Bennetzen JL, H BD: Codon selection in yeast.

    J Biol Chem 1982, 257:3026-3031. PubMed Abstract | Publisher Full Text OpenURL

  8. Nicorici Daniel, A J: Segmentation of DNA into Coding and Noncoding Regions Based on Recursive Entropic Segmentation and Stop-Codon Statistics.

    EURASIP Journal on Applied Signal Processing 2004, 1:81-91. OpenURL

  9. Staden R: Measurements of the effects that coding for a protein has on a DNA sequence and their use for finding genes.

    Nucleic Acid Res 1984, 12:551-567. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  10. Li Wentian , P B-G, Haghighi F, Grosse I: Applications of recursive segmentation to the analysis of DNA sequences.

    Computers and Chemistry 2002, 26(5):491-510. PubMed Abstract | Publisher Full Text OpenURL

  11. Nicorici D, B JA, Astola J, Mitra SK: Finding borders between coding and noncoding DNA regions using recursive segmentation and statistics of stop codons.

    Proceedings of the 2003 Finnish Signal Processing Symposium: May 2003 2003; Tampere, Finland231-235. OpenURL

  12. Bernaola-Galvan P, G I, Carpena P, Oliver JL, Roman-Roldan R, Stanley HE: Finding borders between coding and noncoding DNA regions by an entropic segmentation method.

    Phys Rev Lett 2000, 85(6):1342-1345. PubMed Abstract | Publisher Full Text OpenURL

  13. Li W: New stopping criteria for segmenting DNA sequences.

    PhysRevLett 2001, 86(25):5815-5818. PubMed Abstract | Publisher Full Text OpenURL

  14. Bernaola-Galvan P, R-R R, Oliver JL: Compositional segmentation and long-range fractal correlations in DNA sequences.

    PhysRevE 1996, 53(5):5181-5189. PubMed Abstract | Publisher Full Text OpenURL

  15. Ramaswamy R: Prediction of probable genes by Fourier analysis of genomic sequences.

    CABIOS 1997, 13(3):263-270. PubMed Abstract OpenURL

  16. Zhang Jingxiang , X Z: Finding Borders Between Coding and Noncoding DNA Regions By β-KL Divergence.

    ICBBE 2007 2007, 77:286-289. OpenURL

  17. Fickett JW: Recognition of protein coding regions in DNA sequences.

    Nucleic Acids Research 1982, 10(17):5303-5318. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  18. Staden R, M AD: Codon preference and its use in identifying protein coding regions in long DNA sequences.

    Nucleic Acids Research 1982, 10:141-156. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  19. Shepherd JCW: Method to determine the reading frame of a protein from the purine/pyrimidine genome sequence and its possible evolutionary justification.

    Proc Natl Acad Sci 1981, 78:1596-1600. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  20. Herzel H, G I: Measuring correlations in symbolic sequences.

    Physica A 1995, 216:518-542. Publisher Full Text OpenURL

  21. Grantham R, G C, Gouy M, Jacobzone M, Mercier R: Codon catalog usage is a genome strategy modulated for gene expressivity.

    Nucleic Acids Res 1981, 9(1):R43-74. PubMed Abstract | PubMed Central Full Text OpenURL

  22. Grosse I, H H, Buldyrev SV, Stanley HE: Species independence of mutual information in coding and noncoding DNA.

    Phys Rev E 2000, 61(5):5624-5629. PubMed Abstract | Publisher Full Text OpenURL

  23. Voss RF: Evolution of long-range fractal correlations and 1/f noise in DNA base sequences.

    Phys Rev Lett 1992, 68(1):3805-3808. PubMed Abstract | Publisher Full Text OpenURL

  24. Nicorici D, A J, Tabus I: Computational identification of exons in DNA with a hidden Markov model. In Work shop on Genomic Signal Processing and Statistics. Raleigh, NC, USA; 2002. OpenURL

  25. He Y, H AB, Krim H: A generalized divergence measure for robust image registration.

    IEEE Trans Signal Process 2003, 51(5):1211-1220. Publisher Full Text OpenURL

  26. Henderson J, S S, Fasman KH: Finding genes in DNA with a hidden Markov model.

    Journal of Computational Biology 1997, 4(2):127-141. PubMed Abstract | Publisher Full Text OpenURL

  27. Salzberg S, D A, Fasman K, Henderson J: A decision tree system for finding genes in DNA.

    Journal of Computational Biology 1998, 5(4):667-680. PubMed Abstract | Publisher Full Text OpenURL