BACKGROUND
a promoter is a signal element on a dna molecule that specifies a controlling region of a gene where rna polymerase binds to initiate the transcription of the gene.
rna polymerase ii  in eukaryotic cell binds the promoter signals of all protein-coding sequences.
as there is no common collective occurrence of these signals it is extremely difficult to predict the promoter efficiently using an arbitrary or small set of predetermined signals.
the approach adopted in this research considers a much larger set of signals that are not mutually exclusive and increases the success rate of finding a promoter in a given sequence significantly.
a number of methods for the prediction of promoters, tss  and tf  binding sites in eukaryotic dna sequence presently exist  <cit> .
although contemporary algorithms have improved on predictive ability, a lot more research is required to achieve satisfactory levels of performance.
many general-purpose promoter prediction implementations can typically recognize only ~50% of the promoters with a false positive  rate of ~ <dig> per 700– <dig> bp  <cit> .
a stepwise strategy method corepromoter based on localization of functional promoter and tss region was suggested by zhang.
<cit>  the use of markov chain models in promoter prediction tools by ohler et al.
<cit>  improved the results slightly but they identified the same 50% of promoters from the dataset analyzed by fickett and hatzigeorgiou  <cit> , while having a false positive prediction rate of 1/ <dig> bp.
later ohler  <cit>  improved the result by incorporating some prominent motifs such as tata, dre, inr, dpe, mte  <cit>  to discriminate promoter and non-promoter using hidden markov model.
also, another promoter identification program, promoter  <dig> , designed by knudsen  <cit>  applied a combination of neural networks and genetic algorithms.
promoter  <dig>  was experienced on recognizing promoters in the complete adenovirus genome .
the program predicted some known promoter sites such as tata box, cap site, ccaat box and gc box on the positive strand and  <dig> false promoters.
pandey et.
al.
<cit>  analyzed plant promoters in terms of their structural and sequence dependent properties like curvature and periodicity.
promoters with tata box and without tata box have also been analyzed in this regard.
since completion of the sequencing of the human genome, the efficiency of promoter prediction tools became critical and still poses a major challenge.
promoterinspector program  <cit>  was the first software tool used to identify the promoters in human chromosome  <dig>  this method identifies ~50% of known promoters as genomic regions up to  <dig> kb in length by discriminating them from the exon, intron and 3'-untranslated region  sequences.
davuluri et al.
proposed firstef program which included a decision tree based on structural and compositional discriminating features such as cpg islands, promoter regioins and first splice-donor sites  <cit> .
in  <dig>  bajic et al.
<cit>  reported the dragon promoter finder  program  <cit> , which used sensors for three functional regions: promoters, exons and introns.
their findings claimed that dbf has higher accuracy than three other promoter finding programs which they had investigated namely: nnpp  <dig>   <cit>  promoter <dig>   <cit> and promoterinspector  <cit> .
also in  <dig>  down and hubbard  <cit>  reported a novel hybrid machine-learning method capable of predicting > 50% of human tss with a specificity of > 70%.
xiao et.
al.
<cit>  predicted pol ii promoter sequence by cooperating transcription factor binding site  with the promoter sequence.
they have identified that nearly 71% of the promoter sequences contain transcription factor with known cooperation.
rajeev et.
al.
<cit>  proposed another promoter identification tool where they identified the promoter sequence from non promoter sequence using non linear time series descriptor using non linear machine learning algorithm, such as svm.
their approach showed 87% accuracy in  <dig> fold cross-validation test with an independent test having accuracy of near 85% in identifying promoter and non-promoter.
yet another promoter prediction tool, tssp-tcm has been trained and adapted for plants  <cit> .
in the test set of tata promoters, the program correctly predicted tss for  <dig> out of  <dig>  genes.
for  <dig> tata-less promoters, tsss were predicted for  <dig> out of  <dig>  genes, including  <dig> cases of  <dig> bp distance between annotated and predicted tsss.
a different approach was undertook by gershenzon and ioshikhes, besides identifying core promoter elements in human they tried to demonstrate the synergistic effect of the combination of several core promoter elements for the initiation of the transcription  <cit> .
while jin et al.
later in their study revealed a three-way synergistic effect using a degree of conservation between orthologous mouse and human sequence  <cit> .
the proposed method in this paper demonstrates that with the aid of prommachine- a machine-learning tool, promoters are distinguished from the non-promoter sequences on the basis of abundance of some characteristic 4-mer motifs.
the prommachine is trained with  <dig> distinguishing 4-mers that can discriminate between promoter and non promoter.
using this knowledge the machine-learning tool can efficiently decide whether a given sequence contains a promoter or not.
with high sensitivity, specificity and accuracy the proposed approach shows high efficiency in promoter prediction in eukaryotic genomic sequences.
being applicable for any reasonable length  of given sequences this approach can become a dynamic tool for finding promoters of other eukaryotes.
<dig>
RESULTS
<dig>  prediction accuracy
different statistical measures have been used to analyze the performance of the proposed model.
in order to test the prediction accuracy,  <dig> known promoter sequences from each species i.e.
human, mouse, plant, rat and drosophila as well as  <dig> known non-promoter sequences for each of those five species are initially taken for the test.
then taking one species at a time i.e.
<dig> known promoter and  <dig> known non-promoter for that species, the previously trained model is applied and tested upon these test data.
it is to be mentioned that these test dataset are completely independent from the training set.
in most of the case the model correctly predicts the promoters and non-promoter in these test data.
sensitivity and specificity are two widely used techniques for performance evaluation, which are defined as following, where tp is true positive, fn is false negative, fp is false positive and tn is truly negative.
sensitivity = tp/
specificity = tn/
the results of sensitivity and specificity for all the species are recorded in table  <dig>  as it can be seen from table  <dig> a sensitivity of 86% and specificity of 90% for plant; 96% and 92% for drosophila; 88% and 92% for human; 78% and 84% for mouse and finally 82% and 80% for rat are achieved, which demonstrates the high performance of the proposed model for correctly identifying a promoter as a promoter and a non-promoter as a non-promoter.
also it demonstrates that, the model is less prone to false positive and false negative.
besides sensitivity and specificity, a number of other statistical procedures such as correlation coefficient , false positive rate , false negative rate , positive predictive value , negative predictive value , positive likelihood  and negative likelihood  are also calculated for each of the species.
the fpr values for each species show that the proportions of negative instances that are erroneously reported as being positive are low.
the ppv values in table  <dig> shows that about 80% times a promoter is correctly identified as promoter in rat, whereas this accuracy increased to 83% in mouse, 89% in plant and 92% for both drosophila and human.
however, because the ppv depends on the prevalence of the dataset, a same proportion of positive and negative data are taken for the test.
the likelihood ratios are also calculated, as it is independent of the prevalence .
cross validation is another statistical process for estimating the predictive accuracy of a classifier on data.
in this process, the whole dataset i.e.
both positive and negative data for a given species is partitioned into subsets in such a way that the analysis is made on one subset while all the other subsets are used for training.
to analyze the performance of some other n-mers; first,  <dig>   <dig> and 5-mers are generated and tested against the same dataset.
a cross-validation accuracy of  <dig> %,  <dig> % and  <dig> % are obtained, which justifies the use of 4-mer motif for proposed promoter prediction methodology.
a differential hexamer technique was used earlier by hutchinson for identifying vertebrate promoter; on  <dig> test sequences he correctly distinguished  <dig> promoters as true positive whereas  <dig> were false positive, which gave him a sensitivity of  <dig> %.
the result improved up to  <dig> % when only sequences of length above  <dig>  were considered  <cit> .
also chan and kibler have used 6-mer distribution for identifying cis-regulatory motifs in drosophila and obtained a sensitivity and specificity of  <dig> %,  <dig> %, and the ppv was  <dig> %  <cit> , which is significantly improved upon by the method proposed in this paper .
this supports the significance of 4-mer distribution used here compared to other n-mers.
the rational of 4-mer giving better results is, for lower or higher mers the frequency distribution of the feature set are not descriptive enough for the svm classifier to distinguish promoters from the non-promoters.
moreover a significant overrepresentation of 4-mer motif such as ataa, tata or tcag was confirmed earlier by ohler et al.
<cit> , it may also speculate the justification of 4-mer performing better than some other n-mers.
when the model presented in the paper is trained with the  <dig> discriminating 4-mer features of a given species and tested; the results from the model are promising.
a 7-fold cross validation for all five species displays accuracy of  <dig> %,  <dig> %,  <dig> %,  <dig> % and  <dig> % for plant, drosophila, human, mouse and rat respectively .
the high-percentage of correct value, correlation coefficient for this proposed model clearly indicates that calculated frequencies of 4-mer sequences are capable of discriminating with high accuracy between promoter and non-promoter regions.
the top five 4-mers which are most frequent in the promoter dataset of all five species are listed in table  <dig>  a union of these motifs across the species results in a total of  <dig> motifs all together.
the positional distributions of these  <dig> motifs for each of the five species are shown in figure  <dig>  although measuring the conservation of these 4-mer motifs across different species is outside the scope of this research but positional conserveness of the  <dig> motifs in each five species i.e.
where in the sequence dataset for a given species these  <dig> motifs mostly exists, can be observed from the figures presented.
also the top  <dig> discriminating features  between plant promoter and plant non-promoter are listed in table  <dig>
in addition, the proposed method identifies the location of the tata box using tata weighted matrix  <cit> .
for this a sliding window technique is used with a window size of  <dig> and cutoff score of the matrix is selected as  <dig> .
each time the window is shifted by  <dig> bp across the sequences.
the windows with a tata box score exceeding the cutoff score are accepted as tata box.
this indicates the possible location of promoter.
as tata box is usually located  <dig> base pairs upstream to the transcription start site , using this information tss can also be identified in the tata box containing promoter sequences.
<dig>  comparison with existing methods
similarity based methods such as blast is normally the first choice to distinguish between two different set of data.
but no significant distinguishable features are identified using promoter-promoter, non-promoter – non-promoter and promoter-non-promoter blast search.
to date various algorithms and methods are used for promoter prediction.
using these algorithms, some widely used promoter prediction tools have been developed, e.g.
soft berry, dragon promoter finder, neural network promoter prediction, promoter  <dig>  prediction server and promoter scan.
however, the model proposed here demonstrates even better performance than these tools.
the results shown in table  <dig> clearly indicate that the prediction accuracy of the model is much higher in comparison to other tools.
there are a number of key features in this model which contribute to the improved performance such as the use of larger sets  of discriminating features between promoters and non promoters and the utilisation of the supervised learning system- svm.
other prominent promoter prediction tools use either a statistical approach or neural network such as softberry .
however, svm has out-performed both of these approaches in pattern matching and supervised learning.
also most of the afore-mentioned approaches are limited in the number of features used as the promoter signals such as tata box or initiator etc.
in contrast, the model described in this paper, uses svm to detect a greater set of signals than previously attempted to date to decide whether a sequence is a promoter or non promoter.
the additional features developed in the course of this research show that a significant increase in promoter prediction can be achieved using this method.
<dig>
CONCLUSIONS
the successful prediction of promoters with high accuracy using frequency distribution of 4-mer sequences shows that this novel method has significant merit as an approach for successful eukaryotic promoter prediction.
the principal objective of this work was to develop an efficient tool that can discriminate between promoter and non-promoter in an unknown sequence with high accuracy.
the high accuracy of promoter prediction in plant, human, drosophila, mouse and rat  using our approach has validated the use of frequency distribution of 4-mers in discrimination between promoter and non promoter regions in eukaryotic sequences.
however, though the approach is very efficient in predicting the presence of promoter in a given sequence, it cannot locate the position of the tss when tata box is not present.
this challenge should be resolved in future if other signals can be characterized for specific position in the sequence.
the improvement of identification of all tss site will be considered in future using various published databases such as mpromdb, tred, dbtss etc.
though the epd is known to have a tata bias  <cit> , the result presented in this paper is not likely to vary much using other database, cause the proposed method does not rely on the identification of the tata box for discriminating promoters from non-promoters.
incorporating comparative genomics certainly can extend the method.
it would reduce the false discovery rate further, as at a moderate phylogenetic distance the functional elements are known to be more conserved for orthologous gene  <cit> .
nevertheless, it is suggested that the approach proposed here would be a extremely useful and an efficient tool to meet the demands of the molecular biologists.
<dig> methods
<dig>  promoter sequence database
to accomplish the task of rna polymerase ii promoter prediction, the plant promoter dataset is taken from plantprom database  <cit> .
a total of  <dig> entries of plant promoter sequence are obtained from plantprom db.
also to make the model applicable to predict promoters in drosophila, human, mouse and rat, a number of  <dig>   <dig>   <dig> and  <dig> promoter sequences respectively are collected from epd   <cit> .
these data are used for training the proposed model and tested against mutually exclusive data from the same data source.
sequence segments from - <dig> to + <dig> relative to tss are considered from the epd for the experiment.
on the other hand - <dig> to + <dig> segments are taken from plantprom db.
<dig>  non-promoter sequence database
the non-promoter sequences are extracted from unigene database that belongs to embl.
cdss are the best non promoter sequence because it contain no promoter sequences and so for every organism specific model an equal number of cdss are usaed in the training process as the non promoter sequence.
<dig>  support vector machine
svm, a supervised machine-learning technique has been used for discriminating between promoter and non-promoter sequences.
svm classifiers solve multi-class classification problems using the structural minimization principle.
given a training set in a vector space, svms can find the best decision hyper plane, which separates two classes.
the quality of the decision hyper plane depends on the difference margin between the two hyperplanes defined by the svm  <cit> .
for a typical learning task p = pp, an inductive svm learner aims to build a decision function.
here xi ∈ x are objects of the training set and yi ∈ {- <dig>  1} are their known classes.
fn: x → {- <dig> } based on a training set s which is fn = fn
where strain = , ,..., .
it is necessary to select a kernel function and the regularization parameter in each binary classifier and in this instance, the radial basis function   <cit>  is selected as the kernel function.
because the rbf kernel nonlinearly maps samples into a higher dimensional space, it can handle cases where the relations are non linear.
the svm classification problem can be formulated in terms of a convex quadratic optimization problem as
max
in the above equation n is the total number of input vectors, αi is any real value that maximizes the function, xi is any real number as the input vectors and yi be their corresponding target class, which is either - <dig> or  <dig> in binary classifier.
libsvm – a library for svm developed by chih-chung chang and chih-jen lin  is used to train in a supervised manner on a collection of promoter and non-promoter training sequences.
the model developed from this training set is then used to predict the promoter sequences from a test sequence.
<dig>  procedure
<dig> .
<dig> selection of features
taking all four nucleotides  a  <dig> =  <dig> different 4-mers combinations are generated.
in order to select most discriminating features the frequency of each of these features are calculated.
in this procedure the fij and fnij are calculated, where fij is the frequency of each ith 4-mer motif  in jth known promoter sequence and fnij is the frequency of each ith 4-mer in jth known non-promoter sequence.
then this fij and fnij are used for calculating pi and npi for which the equation is given below.
pi=∑j=1nfij
npi=∑j=1nfnij
where, pi is the summation of frequencies in n promoter sequences for each ith 4-mer  and npi is the summation of frequencies in n non-promoter sequences for each ith 4-mer.
the absolute difference between the numbers of occurrences of all  <dig> possible combinations of 4-mers in known promoter and non-promoter sequences are calculated using diffi = |pi - npi |, where diffi is the absolute difference of the occurrence of ith 4-mer between known promoter and non promoter.
the value of these  <dig> diffi are then sorted in descending order and the top  <dig> 4-mers for which the diffi is maximum are selected as features of the proposed model.
all the selected n known promoter and n known non-promoter sequence are of length  <dig> bp.
a single strand search of  <dig> bp window has been conducted during the frequency calculation.
<dig> .
<dig> training with svm
the frequency fij of the top  <dig> 4-mer motifs are used to find the promoter  in a given sequence.
here i range from  <dig> to  <dig> and j range from  <dig> to n. a scaling factor is introduced to bring the numerical value of fij in between  <dig> and  <dig> using equation  <dig>
dij=fij−minmax
here, the numerical value of dij ranges from  <dig> to  <dig>
here, the value dij is the scaled value of ith 4-mer in jth sequence.
a complete set of dij values of all the promoter and non-promoter sequences is then used to train the svm for a given species.
libsvm  <cit> , which is an efficient software for svm classification and regression, is trained with the frequency patterns of the  <dig> characteristic 4-mer motifs in known promoters and non-promoters.
a model is then built that can be used to distinguish between promoter and non promoter in test sequences.
<dig> .
<dig> testing
to substantiate the machine-learning model, jackknife  <cit>  validation could be the ideal process.
jackknife is a method that creates a series of statistics from a single data set by generating that statistic repeatedly on the data set leaving one data value out each time which produces a mean estimate of the parameter and a standard deviation of the estimates of the parameter.
this validation method is time consuming.
so to speed up the process, a 7-fold cross-validation is performed where the whole dataset is divided into  <dig> random parts and  <dig> iterations are performed to predict the accuracy.
in every iteration  <dig> parts are used for training and the other part is used for testing.
in the following iteration the other  <dig> parts are used for training and the remaining part is used for testing.
authors' contributions
fa was responsible for implementation and reviewing the research and manuscript.
smb was responsible for coding, implementation and manuscript preparation.
tj was responsible for implementation, coding and testing.
mmh was responsible for biological information analysis and interpretation.
ms proposed the idea for the research and generated the project specification.
hk reviewed the research manuscript and inspired overall analysis.
rw reviewed the research manuscript, revised previous versions, edited and formatted the final document.

 a promoter is a signal element on a dna molecule that specifies a controlling region of a gene where rna polymerase binds to initiate the transcription of the gene. rna polymerase ii in eukaryotic cell binds the promoter signals of all protein-coding sequences. as there is no common collective occurrence of these signals it is extremely difficult to predict the promoter efficiently using an arbitrary or small set of predetermined signals. the approach adopted in this research considers a much larger set of signals that are not mutually exclusive and increases the success rate of finding a promoter in a given sequence significantly. a number of methods for the prediction of promoters, tss and tf binding sites in eukaryotic dna sequence presently exist. although contemporary algorithms have improved on predictive ability, a lot more research is required to achieve satisfactory levels of performance. many general-purpose promoter prediction implementations can typically recognize only ~50% of the promoters with a false positive rate of ~ per 700– bp. a stepwise strategy method corepromoter based on localization of functional promoter and tss region was suggested by zhang. the use of markov chain models in promoter prediction tools by ohler et al. improved the results slightly but they identified the same 50% of promoters from the dataset analyzed by fickett and hatzigeorgiou, while having a false positive prediction rate of 1/ bp. later ohler improved the result by incorporating some prominent motifs such as tata, dre, inr, dpe, mte to discriminate promoter and non-promoter using hidden markov model. also, another promoter identification program, promoter, designed by knudsen applied a combination of neural networks and genetic algorithms. promoter was experienced on recognizing promoters in the complete adenovirus genome. the program predicted some known promoter sites such as tata box, cap site, ccaat box and gc box on the positive strand and false promoters. pandey et. al. analyzed plant promoters in terms of their structural and sequence dependent properties like curvature and periodicity. promoters with tata box and without tata box have also been analyzed in this regard. since completion of the sequencing of the human genome, the efficiency of promoter prediction tools became critical and still poses a major challenge. promoterinspector program was the first software tool used to identify the promoters in human chromosome this method identifies ~50% of known promoters as genomic regions up to kb in length by discriminating them from the exon, intron and 3'-untranslated region sequences. davuluri et al. proposed firstef program which included a decision tree based on structural and compositional discriminating features such as cpg islands, promoter regioins and first splice-donor sites. in bajic et al. reported the dragon promoter finder program, which used sensors for three functional regions: promoters, exons and introns. their findings claimed that dbf has higher accuracy than three other promoter finding programs which they had investigated namely: nnpp promoter and promoterinspector. also in down and hubbard reported a novel hybrid machine-learning method capable of predicting > 50% of human tss with a specificity of > 70%. xiao et. al. predicted pol ii promoter sequence by cooperating transcription factor binding site with the promoter sequence. they have identified that nearly 71% of the promoter sequences contain transcription factor with known cooperation. rajeev et. al. proposed another promoter identification tool where they identified the promoter sequence from non promoter sequence using non linear time series descriptor using non linear machine learning algorithm, such as svm. their approach showed 87% accuracy in fold cross-validation test with an independent test having accuracy of near 85% in identifying promoter and non-promoter. yet another promoter prediction tool, tssp-tcm has been trained and adapted for plants. in the test set of tata promoters, the program correctly predicted tss for out of genes. for tata-less promoters, tsss were predicted for out of genes, including cases of bp distance between annotated and predicted tsss. a different approach was undertook by gershenzon and ioshikhes, besides identifying core promoter elements in human they tried to demonstrate the synergistic effect of the combination of several core promoter elements for the initiation of the transcription. while jin et al. later in their study revealed a three-way synergistic effect using a degree of conservation between orthologous mouse and human sequence. the proposed method in this paper demonstrates that with the aid of prommachine- a machine-learning tool, promoters are distinguished from the non-promoter sequences on the basis of abundance of some characteristic 4-mer motifs. the prommachine is trained with distinguishing 4-mers that can discriminate between promoter and non promoter. using this knowledge the machine-learning tool can efficiently decide whether a given sequence contains a promoter or not. with high sensitivity, specificity and accuracy the proposed approach shows high efficiency in promoter prediction in eukaryotic genomic sequences. being applicable for any reasonable length of given sequences this approach can become a dynamic tool for finding promoters of other eukaryotes.  prediction accuracy different statistical measures have been used to analyze the performance of the proposed model. in order to test the prediction accuracy, known promoter sequences from each species i.e. human, mouse, plant, rat and drosophila as well as known non-promoter sequences for each of those five species are initially taken for the test. then taking one species at a time i.e. known promoter and known non-promoter for that species, the previously trained model is applied and tested upon these test data. it is to be mentioned that these test dataset are completely independent from the training set. in most of the case the model correctly predicts the promoters and non-promoter in these test data. sensitivity and specificity are two widely used techniques for performance evaluation, which are defined as following, where tp is true positive, fn is false negative, fp is false positive and tn is truly negative. sensitivity = tp/ specificity = tn/ the results of sensitivity and specificity for all the species are recorded in table as it can be seen from table a sensitivity of 86% and specificity of 90% for plant; 96% and 92% for drosophila; 88% and 92% for human; 78% and 84% for mouse and finally 82% and 80% for rat are achieved, which demonstrates the high performance of the proposed model for correctly identifying a promoter as a promoter and a non-promoter as a non-promoter. also it demonstrates that, the model is less prone to false positive and false negative. besides sensitivity and specificity, a number of other statistical procedures such as correlation coefficient, false positive rate, false negative rate, positive predictive value, negative predictive value, positive likelihood and negative likelihood are also calculated for each of the species. the fpr values for each species show that the proportions of negative instances that are erroneously reported as being positive are low. the ppv values in table shows that about 80% times a promoter is correctly identified as promoter in rat, whereas this accuracy increased to 83% in mouse, 89% in plant and 92% for both drosophila and human. however, because the ppv depends on the prevalence of the dataset, a same proportion of positive and negative data are taken for the test. the likelihood ratios are also calculated, as it is independent of the prevalence. cross validation is another statistical process for estimating the predictive accuracy of a classifier on data. in this process, the whole dataset i.e. both positive and negative data for a given species is partitioned into subsets in such a way that the analysis is made on one subset while all the other subsets are used for training. to analyze the performance of some other n-mers; first, and 5-mers are generated and tested against the same dataset. a cross-validation accuracy of %, % and % are obtained, which justifies the use of 4-mer motif for proposed promoter prediction methodology. a differential hexamer technique was used earlier by hutchinson for identifying vertebrate promoter; on test sequences he correctly distinguished promoters as true positive whereas were false positive, which gave him a sensitivity of %. the result improved up to % when only sequences of length above were considered. also chan and kibler have used 6-mer distribution for identifying cis-regulatory motifs in drosophila and obtained a sensitivity and specificity of %, %, and the ppv was %, which is significantly improved upon by the method proposed in this paper. this supports the significance of 4-mer distribution used here compared to other n-mers. the rational of 4-mer giving better results is, for lower or higher mers the frequency distribution of the feature set are not descriptive enough for the svm classifier to distinguish promoters from the non-promoters. moreover a significant overrepresentation of 4-mer motif such as ataa, tata or tcag was confirmed earlier by ohler et al., it may also speculate the justification of 4-mer performing better than some other n-mers. when the model presented in the paper is trained with the discriminating 4-mer features of a given species and tested; the results from the model are promising. a 7-fold cross validation for all five species displays accuracy of %, %, %, % and % for plant, drosophila, human, mouse and rat respectively. the high-percentage of correct value, correlation coefficient for this proposed model clearly indicates that calculated frequencies of 4-mer sequences are capable of discriminating with high accuracy between promoter and non-promoter regions. the top five 4-mers which are most frequent in the promoter dataset of all five species are listed in table a union of these motifs across the species results in a total of motifs all together. the positional distributions of these motifs for each of the five species are shown in figure although measuring the conservation of these 4-mer motifs across different species is outside the scope of this research but positional conserveness of the motifs in each five species i.e. where in the sequence dataset for a given species these motifs mostly exists, can be observed from the figures presented. also the top discriminating features between plant promoter and plant non-promoter are listed in table in addition, the proposed method identifies the location of the tata box using tata weighted matrix. for this a sliding window technique is used with a window size of and cutoff score of the matrix is selected as. each time the window is shifted by bp across the sequences. the windows with a tata box score exceeding the cutoff score are accepted as tata box. this indicates the possible location of promoter. as tata box is usually located base pairs upstream to the transcription start site, using this information tss can also be identified in the tata box containing promoter sequences. comparison with existing methods similarity based methods such as blast is normally the first choice to distinguish between two different set of data. but no significant distinguishable features are identified using promoter-promoter, non-promoter – non-promoter and promoter-non-promoter blast search. to date various algorithms and methods are used for promoter prediction. using these algorithms, some widely used promoter prediction tools have been developed, e.g. soft berry, dragon promoter finder, neural network promoter prediction, promoter prediction server and promoter scan. however, the model proposed here demonstrates even better performance than these tools. the results shown in table clearly indicate that the prediction accuracy of the model is much higher in comparison to other tools. there are a number of key features in this model which contribute to the improved performance such as the use of larger sets of discriminating features between promoters and non promoters and the utilisation of the supervised learning system- svm. other prominent promoter prediction tools use either a statistical approach or neural network such as softberry. however, svm has out-performed both of these approaches in pattern matching and supervised learning. also most of the afore-mentioned approaches are limited in the number of features used as the promoter signals such as tata box or initiator etc. in contrast, the model described in this paper, uses svm to detect a greater set of signals than previously attempted to date to decide whether a sequence is a promoter or non promoter. the additional features developed in the course of this research show that a significant increase in promoter prediction can be achieved using this method.  the successful prediction of promoters with high accuracy using frequency distribution of 4-mer sequences shows that this novel method has significant merit as an approach for successful eukaryotic promoter prediction. the principal objective of this work was to develop an efficient tool that can discriminate between promoter and non-promoter in an unknown sequence with high accuracy. the high accuracy of promoter prediction in plant, human, drosophila, mouse and rat using our approach has validated the use of frequency distribution of 4-mers in discrimination between promoter and non promoter regions in eukaryotic sequences. however, though the approach is very efficient in predicting the presence of promoter in a given sequence, it cannot locate the position of the tss when tata box is not present. this challenge should be resolved in future if other signals can be characterized for specific position in the sequence. the improvement of identification of all tss site will be considered in future using various published databases such as mpromdb, tred, dbtss etc. though the epd is known to have a tata bias, the result presented in this paper is not likely to vary much using other database, cause the proposed method does not rely on the identification of the tata box for discriminating promoters from non-promoters. incorporating comparative genomics certainly can extend the method. it would reduce the false discovery rate further, as at a moderate phylogenetic distance the functional elements are known to be more conserved for orthologous gene. nevertheless, it is suggested that the approach proposed here would be a extremely useful and an efficient tool to meet the demands of the molecular biologists. methods promoter sequence database to accomplish the task of rna polymerase ii promoter prediction, the plant promoter dataset is taken from plantprom database. a total of entries of plant promoter sequence are obtained from plantprom db. also to make the model applicable to predict promoters in drosophila, human, mouse and rat, a number of and promoter sequences respectively are collected from epd. these data are used for training the proposed model and tested against mutually exclusive data from the same data source. sequence segments from - to + relative to tss are considered from the epd for the experiment. on the other hand - to + segments are taken from plantprom db. non-promoter sequence database the non-promoter sequences are extracted from unigene database that belongs to embl. cdss are the best non promoter sequence because it contain no promoter sequences and so for every organism specific model an equal number of cdss are usaed in the training process as the non promoter sequence. support vector machine svm, a supervised machine-learning technique has been used for discriminating between promoter and non-promoter sequences. svm classifiers solve multi-class classification problems using the structural minimization principle. given a training set in a vector space, svms can find the best decision hyper plane, which separates two classes. the quality of the decision hyper plane depends on the difference margin between the two hyperplanes defined by the svm. for a typical learning task p = pp, an inductive svm learner aims to build a decision function. here xi ∈ x are objects of the training set and yi ∈ {- 1} are their known classes. fn: x → {- } based on a training set s which is fn = fn where strain =,,...,. it is necessary to select a kernel function and the regularization parameter in each binary classifier and in this instance, the radial basis function is selected as the kernel function. because the rbf kernel nonlinearly maps samples into a higher dimensional space, it can handle cases where the relations are non linear. the svm classification problem can be formulated in terms of a convex quadratic optimization problem as max in the above equation n is the total number of input vectors, αi is any real value that maximizes the function, xi is any real number as the input vectors and yi be their corresponding target class, which is either - or in binary classifier. libsvm – a library for svm developed by chih-chung chang and chih-jen lin is used to train in a supervised manner on a collection of promoter and non-promoter training sequences. the model developed from this training set is then used to predict the promoter sequences from a test sequence. procedure. selection of features taking all four nucleotides a = different 4-mers combinations are generated. in order to select most discriminating features the frequency of each of these features are calculated. in this procedure the fij and fnij are calculated, where fij is the frequency of each ith 4-mer motif in jth known promoter sequence and fnij is the frequency of each ith 4-mer in jth known non-promoter sequence. then this fij and fnij are used for calculating pi and npi for which the equation is given below. pi=∑j=1nfij npi=∑j=1nfnij where, pi is the summation of frequencies in n promoter sequences for each ith 4-mer and npi is the summation of frequencies in n non-promoter sequences for each ith 4-mer. the absolute difference between the numbers of occurrences of all possible combinations of 4-mers in known promoter and non-promoter sequences are calculated using diffi = |pi - npi |, where diffi is the absolute difference of the occurrence of ith 4-mer between known promoter and non promoter. the value of these diffi are then sorted in descending order and the top 4-mers for which the diffi is maximum are selected as features of the proposed model. all the selected n known promoter and n known non-promoter sequence are of length bp. a single strand search of bp window has been conducted during the frequency calculation.. training with svm the frequency fij of the top 4-mer motifs are used to find the promoter in a given sequence. here i range from to and j range from to n. a scaling factor is introduced to bring the numerical value of fij in between and using equation dij=fij−minmax here, the numerical value of dij ranges from to here, the value dij is the scaled value of ith 4-mer in jth sequence. a complete set of dij values of all the promoter and non-promoter sequences is then used to train the svm for a given species. libsvm, which is an efficient software for svm classification and regression, is trained with the frequency patterns of the characteristic 4-mer motifs in known promoters and non-promoters. a model is then built that can be used to distinguish between promoter and non promoter in test sequences.. testing to substantiate the machine-learning model, jackknife validation could be the ideal process. jackknife is a method that creates a series of statistics from a single data set by generating that statistic repeatedly on the data set leaving one data value out each time which produces a mean estimate of the parameter and a standard deviation of the estimates of the parameter. this validation method is time consuming. so to speed up the process, a 7-fold cross-validation is performed where the whole dataset is divided into random parts and iterations are performed to predict the accuracy. in every iteration parts are used for training and the other part is used for testing. in the following iteration the other parts are used for training and the remaining part is used for testing. authors' contributions fa was responsible for implementation and reviewing the research and manuscript. smb was responsible for coding, implementation and manuscript preparation. tj was responsible for implementation, coding and testing. mmh was responsible for biological information analysis and interpretation. ms proposed the idea for the research and generated the project specification. hk reviewed the research manuscript and inspired overall analysis. rw reviewed the research manuscript, revised previous versions, edited and formatted the final document.
the high sensitivity and specificity value of and for plant; and for drosophila; and for human; and for mouse and and for rat demonstrate that this technique is less prone to false positive results and exhibits better performance than many other tools. the increased availability of sequence data for various eukaryotic organisms in recent years has necessitated for better tools and techniques for the prediction and analysis of promoters in eukaryotic sequences. secondly improved performance can be achieved by enhancing the predictive ability of the machine learning algorithms used. the high sensitivity and specificity indicate that 4-mer frequencies in conjunction with supervised machine-learning methods can be beneficial in the identification of rna pol ii promoters comparative to other methods.
 eukaryotic promoter prediction using computational analysis techniques is one of the most difficult jobs in computational genomics that is essential for constructing and understanding genetic regulatory networks. the increased availability of sequence data for various eukaryotic organisms in recent years has necessitated for better tools and techniques for the prediction and analysis of promoters in eukaryotic sequences. many promoter prediction methods and tools have been developed to date but they have yet to provide acceptable predictive performance. one obvious criteria to improve on current methods is to devise a better system for selecting appropriate features of promoters that distinguish them from non-promoters. secondly improved performance can be achieved by enhancing the predictive ability of the machine learning algorithms used.  in this paper, a novel approach is presented in which 4-mer motifs in conjunction with a non-linear machine-learning algorithm utilising a support vector machine are used to distinguish between promoter and non-promoter dna sequences. by applying this approach to plant, drosophila, human, mouse and rat sequences, the classification model has showed 7-fold cross-validation percentage accuracies of %, %, %, % and % respectively. the high sensitivity and specificity value of and for plant; and for drosophila; and for human; and for mouse and and for rat demonstrate that this technique is less prone to false positive results and exhibits better performance than many other tools. moreover, this model successfully identifies location of promoter using tata weight matrix.  the high sensitivity and specificity indicate that 4-mer frequencies in conjunction with supervised machine-learning methods can be beneficial in the identification of rna pol ii promoters comparative to other methods. this approach can be extended to identify promoters in sequences for other eukaryotic genomes. 