Skip to main content
  • Research article
  • Open access
  • Published:

Omni-PolyA: a method and tool for accurate recognition of Poly(A) signals in human genomic DNA

Abstract

Background

Polyadenylation is a critical stage of RNA processing during the formation of mature mRNA, and is present in most of the known eukaryote protein-coding transcripts and many long non-coding RNAs. The correct identification of poly(A) signals (PAS) not only helps to elucidate the 3′-end genomic boundaries of a transcribed DNA region and gene regulatory mechanisms but also gives insight into the multiple transcript isoforms resulting from alternative PAS. Although progress has been made in the in-silico prediction of genomic signals, the recognition of PAS in DNA genomic sequences remains a challenge.

Results

In this study, we analyzed human genomic DNA sequences for the 12 most common PAS variants. Our analysis has identified a set of features that helps in the recognition of true PAS, which may be involved in the regulation of the polyadenylation process. The proposed features, in combination with a recognition model, resulted in a novel method and tool, Omni-PolyA. Omni-PolyA combines several machine learning techniques such as different classifiers in a tree-like decision structure and genetic algorithms for deriving a robust classification model. We performed a comparison between results obtained by state-of-the-art methods, deep neural networks, and Omni-PolyA. Results show that Omni-PolyA significantly reduced the average classification error rate by 35.37% in the prediction of the 12 considered PAS variants relative to the state-of-the-art results.

Conclusions

The results of our study demonstrate that Omni-PolyA is currently the most accurate model for the prediction of PAS in human and can serve as a useful complement to other PAS recognition methods. Omni-PolyA is publicly available as an online tool accessible at www.cbrc.kaust.edu.sa/omnipolya/.

Background

Polyadenylation is an essential stage of RNA processing during the formation of mature mRNA and occurs in most of the known eukaryotic mRNA sequences [1] as well as in many long non-coding RNAs [2]. The polyadenylation process occurs during RNA processing and involves two stages: 1) cleavage of the primary transcript and 2) the polymerization of an adenosine tail at the downstream of the cleaved mRNA in the case of protein-coding transcripts [3]. The necessary proteins needed for an efficient and accurate cleavage/polyadenylation event include, among others, the cleavage and polyadenylation specific factor, cleavage stimulation factor, cleavage factors I and II, and poly(A) polymerase [4,5,6] (Fig. 1a). Cleavage and polyadenylation specific factor recognizes and binds to PAS upstream of the cleavage site [7]. Although isolated cleavage and polyadenylation specific factor binds to PAS, the strength of such binding considerably increases when acting along with cleavage stimulation factor [6]. Endonucleic cleavage at the 3′-end is performed by cleavage factors I and II [5]. Refer to reviews [6, 7] for a detailed description of these factors and their interactions.

Fig. 1
figure 1

Cleavage/polyadenylation event and general sequence consensus. a Diagram representing the proteins involved during the cleavage/polyadenylation event. b Schematic representation of the region surrounding the PAS [7]

Although PAS are necessary for the 3′-end polyadenylation, other sequence elements have shown to be necessary for a fully functional cleavage/polyadenylation event, especially the downstream sequence elements, characterized by a GU-rich sequence located up to 20 nucleotides (nt) downstream of the cleavage site [7]. The distance between the PAS and the GU-rich sequence primarily enhances and determines the 3′-end formation [8, 9]. Moreover, upstream sequence elements found upstream of the PAS, are U enriched and often affect the efficiency of the 3′-end formation [8, 10,11,12,13]. Finally, the short sequence in the 3′-end of the actual cleavage site may as well have an impact on the efficiency of the process [7, 14] (Fig. 1b). Contrary to the conserved PAS hexamers, downstream sequence elements and upstream sequence elements are highly variable in sequence composition and have not yet been adequately characterized [4]. This sequence variability of the regions flanking PAS causes a major problem in computational prediction of such signals in genomic DNA sequences. Moreover, there is evidence showing that over half of all protein-coding transcripts have alternative PAS [15], resulting in transcripts with variable 3′-end untranslated regions and multiple transcript isoforms [16,17,18,19]. Therefore, the correct identification of PAS not only helps in elucidating the 3′-end boundaries of a gene and regulatory mechanisms but also gives an insight into the multiple isoforms resulting from alternative PAS. Furthermore, an accurate predictive model of PAS would help in the identification of PAS for transcripts containing premature termination codons, which are degraded by cellular mechanisms [20].

During the past few decades, several models for the in-silico prediction of PAS in genomic DNA and mRNA sequences have been proposed. These models make use of the sequence surrounding the PAS to differentiate true PAS from pseudo-PAS, i.e., hexamers (AATAA, ATTAAA, etc.) that are identical to true PAS hexamers but have no link to the 3′-end polyadenylation. Yada et al. [21] performed a statistical analysis of human genomic sequences surrounding the true PAS covering the region from −80 to +48 where the AATAA hexamer corresponds to positions 0 to 5. From their analysis, the authors observed that C and T/C nucleotides are often found upstream and downstream from the PAS, respectively, and concluded that CAATAAA(T/C) might be regarded as a consensus sequence for PAS. Later, Kondrakhin et al. [22] developed a generalized consensus matrix from a set of 63 vertebrate pre-mRNAs. Elements of the matrix represented the absolute frequencies of nucleotide triplets at each site and were applied to each nucleotide sequence to discriminate PAS from pseudo-PAS. However, when tested on sequences from the adenovirus Ad2 genome, their method produced a high level of false positive predictions. For instance, when their parameters were set to recognize 8 out of 9 true PAS, the model predicted over 1000 pseudo-PAS as true PAS in the Ad2 genome that is of 35,937 nt in length [23]. Subsequently, Salamov and Solovyev [4] developed a model based on a linear discriminant function, from 8 variables defined from a window of 300 nt surrounding the PAS (−100,+200). These variables include, among others, scores from position weight matrices, hexanucleotide composition upstream and downstream, and positional triplet composition. Although the authors achieved better results compared to Yada et al. [21] and Kondrakhin et al. [22], the number of false positives remained relatively high (specificity of ~50%). Tabaska and Zhang [24] developed the polyadq tool consisting of two quadratic discriminant functions (one for each AATAAA and ATTAAA variants) derived from three variables. Their results on two new datasets outperformed the existing methods, especially in the reduction of false positives. In 2003, Legendre and Gautheret [25] developed the ERPIN method based on a probabilistic hidden Markov model. ERPIN used position weight matrices computed for each di-nucleotide in a window of 600 nt surrounding the PAS (−300, +300), and achieved a prediction specificity of 85% for a sensitivity of 56%, resulting in a specificity improvement of 9.7% relative to the polyadq method. Bajic et al. [26] developed the Dragon PolyA tool based on artificial neural networks and self-organized maps for predicting the two most common PAS variants in human (AATAAA and ATTAAA). Their tool improved both sensitivity and specificity by ~5% and 5% on AATAAA variant, respectively, and 11.3% and 7.9% on ATTAAA variant, respectively, relative to those obtained by polyadq. In addition, support vector machine (SVM) approaches have been proposed. Liu et al. [3] derived an SVM model from k-gram and artificially translated amino acid patterns from DNA sequences. Their method includes an entropy-based feature selection process to select the most discriminative features. Their results improved specificity for the three out of four considered datasets compared to ERPIN and polyadq. In 2006, Cheng et al. [27] proposed a polya_svm tool based on an SVM model. Polya_svm is derived from the 15 cis-regulatory elements previously found by Hu et al. [28], and achieved an improvement of sensitivity by 33.8% relative to polyadq while preserving the same specificity. In a subsequent study, Xu et al. [29] used SVM-based models for the prediction of PAS in the chromosomal data, i.e., human chromosome 21, and achieved an accuracy of 83%, sensitivity of 90%, specificity of 76%, and a precision of 80%. Akhtar et al. [30] developed the POLYAR tool based on a linear discriminant analysis model. The tool analyzes 600 nt sequences surrounding the PAS and extracts sequence characteristics using position weight matrices, pentamers composition downstream and upstream of PAS, and the distance between cis-elements, among others. Moreover, authors divided PAS signals into three categories: 1) PAS-strong, containing the two most common variants AATAAA or ATTAAA, 2) PAS-weak, containing any of the other ten remaining variants and, 3) PAS-less, referring to PAS not having any of the 12 most common variants. For PAS-strong, POLYAR made an improvement of sensitivity by a relative 23% compared to polya_svm, at the expense of reducing specificity by 5.6% relative to polyadq. Both POLYAR and polya_svm obtained similar but considerably lower specificity/sensitivity for PAS-weak and PAS-less sequences, demonstrating the need to characterize signals surrounding other PAS variants. In this direction, Kalkatawi et al. [31, 32] developed the Dragon PolyA Spotter tool (DPS) for PAS prediction in human genomic sequences for each of the 12 most common PAS variants separately. Their method used artificial neural networks and random forest models derived from a set of thermodynamic, compositional and statistical features. DPS method considerably outperformed other results obtained by polyadq, POLYAR and poly_svm tools on the most common PAS variant (AATAAA hexamer). Later, Xie et al. [33] used a hidden Markov model (HMM) to extract latent spectral features from DNA sequences, which were subsequently used as input for a linear SVM model (we refer to this model as HMM_SVM hereafter). The authors considered the same genomic DNA sequences for the 12 PAS variants as used in Kalkatawi et al. [31] and reduced the weighted average error rate by 25% relative to the results generated by the DPS tool. Although considerable progress has been made in the PAS prediction, the predictions still produce an unacceptable level of false positives. Moreover, new biological features surrounding PAS may be defined for the development of more efficient PAS recognition models and may provide a better understanding of the polyadenylation machinery. Furthermore, the application of other machine learning methods may result in more accurate and more robust prediction models.

In this study, we proposed and analyzed a new set of features surrounding the PAS in human genomic DNA sequences and in combination with a prediction model we developed a novel prediction method and tool, Omni-PolyA that predicts PAS in human genomic DNA sequences. Omni-PolyA combines different classification models in a tree-like structure. We implemented a general-purpose optimization technique, namely a parallel genetic algorithm [34] to optimize the Omni-PolyA model structure and its parameters for deriving more accurate results. We compared the performance of Omni-PolyA against results obtained by the methods proposed in Kalkatawi et al. [31] (as implemented in the DPS tool) and Xie et al. [33] (as implemented in the HMM_SVM tool), thereby demonstrating the utility of the proposed feature set and the Omni-PolyA model. Our comparison analysis shows that Omni-PolyA reduces the weighted average error rate by 35.37% relative to the state-of-the-art results [31, 33] for the 12 considered PAS variants.

Results

The key contributions of our study are the analysis of the genomic sequences for deriving a new set of features capturing elements important for the identification of PAS and the development of a novel method and tool, Omni-PolyA, designed for the prediction of PAS.

Feature mining from genomic DNA sequences

Tian et al. [15] suggested that polyadenylation events may be determined from the combination of the PAS, DNA elements surrounding PAS, and the binding factors. Clearly, an accurate tool for PAS prediction from genomic DNA sequences would be of great help for real applications, i.e., for finding computationally alternative PAS or as a component of gene finding tools. Therefore, several studies have focused on the identification of cis-elements and significant sequence patterns surrounding the PAS [3, 6, 9, 19, 24,25,26,27,28, 31,32,33]. Notably, Kalkatawi et al. [31] and Xie et al. [33] proposed the most discriminant features and models to date for the prediction of PAS from human genomic sequences. However, selecting the optimal combination of features from a big set and the type of classification model to efficiently utilize them is not trivial. Consequently, we analyzed the genomic DNA sequences flanking the PAS (200 nt long sequences compiled by Kalkatawi et al. [31], see Methods) to obtain new biologically significant features. For better understanding the differences and variations between PAS, we followed the PAS categorization proposed by Akhtar et al. [30]. As such, PAS are divided into two categories: 1) PAS-strong sequences containing AATAAA or ATTAAA variants, and 2) PAS-weak sequences with any of the ten other considered PAS variants. Fig. 2a and b show the different sequence composition between PAS-strong and PAS-weak signal surroundings, respectively. In general, PAS-weak variants show a notably higher enrichment of adenine downstream of PAS, while the upstream region presents greater variability of nucleotide composition compared to PAS-strong. Interestingly, the enrichment of nucleotides A/T downstream of pseudo-PAS (Fig. 2c) may suggest similarities to the functional PAS. Consistent with previous studies [15, 19] this may indicate that there may be functional PAS within coding DNA sequences. However, this hypothesis would have to be tested in a laboratory. It is important to note that Fig. 2 shows the averaged DNA distribution for all sequences. Consequently, not all sequences necessarily contain the same characteristics as observed in the figure. Nonetheless, these observations reveal the main differences between PAS-weak and PAS-strong sequences, which can be used to generate a new set of discriminant features. Based on this analysis, we propose a feature set represented by 218 numeric values (referred to as Omni-PolyA feature set). The Omni-PolyA feature set includes mono-nucleotide and di-nucleotides frequencies in particular regions of the genomic DNA sequences, i.e., downstream, upstream, and in-frame codons with respect to the PAS hexamer, among others. Novel and more specific features determined, make use of the entropy and positional information gain of the nucleotide content to determine the most relevant sequence positions. Moreover, we calculate a sequence score derived from 2-mer weight matrices, which capture the di-nucleotide characteristics for the entire sequence surrounding the PAS. Additional file 1: Table S1 shows the list of the 218 compiled features in Omni-polyA feature set.

Fig. 2
figure 2

DNA sequence composition. Nucleotide distribution from positions −100 to 100 with respect to a poly(a) hexamer. a The nucleotide distribution for the genomic DNA sequences in the PAS-strong category. b The nucleotide distribution for the genomic DNA in the PAS-weak category. c The nucleotide distribution for sequences containing pseudo-PAS. PAS hexamers were removed from the sequence analysis

Omni-PolyA method

In contrast to typical univariate decision trees (where decision splits are usually made from a single feature threshold), omnivariate decision trees allow the tree to use a more complex classification model at each non-terminal node [35]. Each non-terminal node in the omnivariate decision tree may use a different classification model depending on the data reaching the node. Therefore, this procedure offers the advantages of multiple classification models to learn from different subsets of data. While the method is generic, in our implementation we considered the following classification models for non-terminal nodes: 1) C4.5 univariate decision tree, 2) artificial neural network, 3) random forest and 4) multinomial logistic regression model (Fig. 3a). However, the key challenges rely on determining the depth of the tree and the classification models used at each non-terminal node [36,37,38]. For this, we implemented a genetic algorithm for optimizing the following processes: 1) tree pruning, where non-terminal nodes with little or no contribution to the classification performance are deleted, and 2) classification model selection and parameters tuning for each non-terminal node. Figure 3b shows an illustration of the resulting tree after the genetic algorithm optimization.

Fig. 3
figure 3

Omnivariate decision tree. a Illustration of the omnivariate decision tree structure where each non-terminal node may be an artificial neural network, random forest, multinomial logistic regression or a C4.5 model. b An example of the simplified tree resulting from the genetic algorithm optimization

Comparison of performances in PAS prediction

Results reported in this section were obtained from a 5-fold cross-validation in agreement with benchmark results published for DPS and HMM_SVM state-of-the-art tools for poly(A) signal prediction (see Methods). Table 1 shows the prediction error rates achieved by deep neural networks (DNN), Omni-PolyA and, DPS [31, 32]. It is important to note that DPS, DNN, and Omni-PolyA models were derived by using the same 274 features proposed by Kalkatawi et al. [31] and we refer to this set as DPS feature set. Results in Table 1 show that the DPS feature set may be used to derive more accurate classification models than those reported by Kalkatawi et al. [31, 32]. For example, DNNs consist of a multi-layered architecture that transforms the data representation at one level to a higher and more abstract level [39]. Such an approach has shown to outperform different models applied to solve some challenging tasks in the fields of bioinformatics and cheminformatics [39,40,41]. Although DNNs typically require a considerably large training set in order to tune a network properly, DNN reduced the weighted average error rate even for the less common PAS variants and reduced the weighted average of the error rate by 11.32% compared to DPS that used the same feature set. These results demonstrate that different models may be better suited for each of the PAS variants. Furthermore, the Omni-PolyA model derived by using the DPS feature set, reduced the weighted average of the error rate by 26.85% and 33.60%, compared to DPS and HMM_SVM, respectively (Additional file 2: Table S2 shows the results obtained by DPS, DNN and Omni-PolyA using the DPS feature and HMM_SVM). We did not consider the latent spectral features proposed by Xie et al. [33] used to derive HMM_SVM due to the high number of features (~10,000) that would result in an unacceptable compute time for non-linear models used as decision nodes in our Omni-PolyA setup based on genetic algorithms. Therefore, we analyzed the genomic sequences and proposed the Omni-PolyA feature set (see the Feature mining from genomic DNA sequences subsection). To assess the discrimination capabilities of the Omni-PolyA feature set, we used these features to derive an Omni-PolyA model for each PAS variant. Results in Table 1 show that the classification error rate, reduced for 8 of the 12 variants with an average improvement of 7.74%, relative to Omni-PolyA derived by using the existing DPS feature set (fifth and sixth columns, Table 1). Overall, this represents a relative error rate improvement of 32.51% and 38.74% over DPS and HMM_SVM state-of-the-art tools, respectively (Additional file 3: Table S3 shows the false positive and false negative rates). Nevertheless, Omni-PolyA achieved an inferior performance on some of the variants from the PAS-weak category compared to HMM_SVM. This may be due to the limited amount of training data for variants with only few hundred sequences. In general, overfitting the training and validation data, outliers, and noise, are among the most relevant problems when deriving classification models from small datasets. For example, in the less frequent PAS variants, the Omni-PolyA algorithm would train a model using a very limited number of sequences, which may derive an unstable and inaccurate classifier. This data limitation is known to considerably affect the performance of the classification models [42, 43]. The creation of synthetic samples, data processing (i.e., noise removal, feature selection, etc.), and data pooling (from similar sub-problems) are among the possible solutions to address the small data problems in deriving more stable classifiers. Notably, the straightforward possibility for deriving robust models for the less frequent PAS variants is to pool the data from different PAS variants with similar nucleotide distributions. Additional file 4: Figure S1 shows the nucleotide distribution for the true PAS sequences for the PAS-weak variants in which we observe a consistent nucleotide enrichment pattern in some sequence regions (i.e., enrichment of T nucleotide in the downstream region 25–50 and the overall A and T enrichment in the 200 nt flanking the PAS). Therefore, we expanded the training data by pooling the PAS-weak variants. For instance, to derive a model for AATAGA variant, the data from the nine remaining PAS-weak variants and the respective training portion from AATAGA variant are used for model tuning (see Methods). Table 1 (seventh column) shows that by pooling the PAS-weak variants, the error rate reduced for eight out of the 10 PAS-weak variants. Considering the best performing Omni-PolyA model for each variant, the classification error rate reduced by 35.37% and 41.34% relative to DPS and HMM_SVM, respectively (Table 2). Notably, the largest error reduction was obtained for the two most common PAS variants (PAS-strong), representing an improvement of 36.95% and 49.50% relative to results achieved by the DPS and HMM_SVM tools.

Table 1 Error rate comparison between DPS, DNN and Omni-PolyA derived by using different feature sets from benchmark poly(A) dataset
Table 2 Error rate comparison between best performing Omni-PolyA model and state-of-the-art results in benchmark poly(A) dataset from [31]

The results discussed in Table 1 and Table 2 are obtained by using the data collected by Kalkatawi et al. [31], which has been used as a benchmark in subsequent studies. However, to account for the current annotation for GRCh37, we used the GENCODE Poly(A) feature annotation (release 19) [44] to extract the true PAS sequences (see Methods). Consistent with the previously discussed results, Table 3 shows that Omni-PolyA considerably outperformed the state-of-the-art methods, reducing the error rate by 6.86% and 20.63% relative to the results achieved by DPS and HMM_SVM tools, respectively. Notably, the largest reduction in the error rate was observed in the PAS-weak variants, in which data were pooled to increase the volume of the training data. For PAS-weak, Omni-PolyA reduced the weighted average of the error rate by 11.50% and 32.64% relative to DPS and HMM_SVM tools, respectively. We want to highlight that when the error rates of a predictive system are less than 20%, it appears to be significantly harder to reduce the error further.

Table 3 Error rate comparison between Omni-PolyA and state-of-the-art methods in GENCODE poly(A) data

Discussion

Survey of PAS-associated genomic features

To capture the sequence variations (Fig. 2a-c) and essential elements involved in the polyadenylation event, we proposed Omni-PolyA feature set containing 218 numeric features (see Additional file 1: Table S1). This set includes features such as positional information gain, scored derived from 2-mer weight matrices, numerical DNA structural profiles, among others. Here, we present a brief survey of the most discriminative features and their biological interpretation.

Positional information gain

It is clear that PAS-weak/strong and pseudo-PAS sequences, on average, differ in sequence composition, most notably in the region (−20, 40) surrounding PAS (Fig. 2a-c). However, we asked if other regions and independent positions may also be relevant for the polyadenylation event. Therefore, to accurately detect the most discriminative positions in a systematic manner, we calculated the information gain independently for each position of the DNA sequence surrounding PAS (see Methods). Consequently, positional information gain detects the positions within the sequence with the highest contribution for differentiating PAS from pseudo-PAS. In agreement with the PAS consensus shown in Fig. 1, positional information gain identified regions in the proximity of PAS for most of the variants (Fig. 4), where proteins involved in the mRNA cleavage/polyadenylation are expected to bind. Interestingly, the AATACA variant revealed the importance of the upstream region, showing that the downstream segment does not contain significant differences with respect to the pseudo-PAS sequences. Finally, the least common variants, e.g., AATAGA, CATAAA, and AGTAAA, show that the relevant discriminatory positions are spread over the 200 nt sequence and show no apparent significant segment. However, a non-linear classification model may be able to use the relevant positions for classification of PAS for each of the variants. As such, we considered the positional information gain to 1) calculate an overall sequence score, and 2) compute the nucleotide frequency for the top ranked positions based on the information gain (see Methods).

Fig. 4
figure 4

Positional information gain. Blue bars represent the 50 most discriminative positions for each variant

Numerical structural profile and PAS

DNA sequences may be converted into a numerical representation to characterize different physical and chemical interactions, i.e., double helix organization, predisposition to interact with other proteins, etc. [45]. These structural profiles have been used in the literature for characterizing genomic signals, i.e., promoter regions [46, 47] and PAS [31, 32]. In total, we used 16 different DNA numeric conversions (conversion tables were obtained from [46]) to define the 200 nt sequences flanking PAS. Notably, A-philicity [48] and protein DNA twist [45] profiles (Fig. 5a-b) show a clear distinguishing pattern for both PAS-strong and PAS-weak compared to pseudo-PAS. A-philicity profile represents the propensity of the DNA to form an A-DNA double helix. Recently, DiMaio et al. [49] suggested that protein binding causes the DNA to adopt an A-form. Similarly, a protein DNA twist profile indicates that both PAS-strong/weak sequences are likely to be deformed by proteins only when within the proximity of PAS (peak around PAS in Fig. 5). Although these structural profiles show a similar pattern, they are, in fact, capturing different information of the sequences and, both profiles accurately detect the region where cleavage and polyadenylation specific factor, cleavage stimulation factor, cleavage factors and poly(A) polymerase proteins are expected to bind. Interestingly, these profiles suggest that the upstream segment in PAS-weak variants is more irregular as opposed to PAS-strong variants. Other numerical conversions were also considered to describe the sequence surrounding PAS. Namely, propeller twist [50], bendability [51], duplex stability free energy [52], DNA bending stiffness [53], stability energy of Z-DNA [54], DNA denaturation [55, 56], nucleosome position preference [57], and base stacking energy [58]. However, B-DNA twist [59], and duplex stability disrupt energy [60] were the least contributing numerical conversions (Additional file 5: Figure S2 and Additional file 6: Figure S3 show the 16 considered structural profiles for PAS-strong and PAS-weak, respectively).

Fig. 5
figure 5

Sequence structural profiles. a A-philicity and protein DNA twist numerical representation of the genomic sequences around the PAS-strong variants (sequences of 200 nt in length). Numerical representations show the actual average values over all sequences for each position. b Similarly to (a), this shows A-philicity and protein DNA twist numerical representations of the genomic sequences for the PAS-weak variants. c The short segments used to calculate a granular structural profile

The next key factor is to determine how to use the information from different structural profiles for capturing relevant information. One option is to consider each numerical position in the sequence as independent input for a classification model (representing 199 features for a di-nucleotide structural profile around the 200 nt sequence). As such, we asked if A-philicity numerical representation alone could be used to identify PAS correctly. However, Omni-PolyA achieved ~24.5% error rate for the prediction of PAS-strong variants (compared to 14% when the DPS feature set is used, see Table 1). Although A-philicity numerical representation can moderately discriminate PAS from pseudo-PAS, a combination of several numerical profiles may grant better discrimination capabilities. However, considering the 199 numeric values for each numerical representation would result in a considerably large feature set (~3100 from the 16 numerical representations). The high number of features would lead to complex models trained to use many irrelevant features (assuming that not all positions contribute to the correct classification of PAS). For this, Kalkatawi et al. [31, 32] combined all di/tri-nucleotide representations from each structural profile into a single score, resulting in one feature per structural profile. Arguably, representing all numeric representations in a sequence by a single score may incur a loss of information. For instance, the discriminant information in the region surrounding the PAS (−20, 40) may be diluted by information from less relevant regions. Therefore, we divided the sequence into sub-sequences of 25 nt and calculated the average of each of these (Fig. 5c), resulting in 8 features for each numerical representation. This procedure not only reduces the number of features but also captures different downstream/upstream elements while minimizing the noise of the independent positions.

Conclusions

In summary, this study shows a comparison of various tools and models applied to the prediction of the 12 most common PAS variants in human genomic DNA sequences. Moreover, by analyzing the differences between PAS-strong, PAS-weak, and pseudo-PAS sequences, we have identified a set of relevant features that may be involved in the regulation of the polyadenylation machinery. In agreement with the consensus of the mammalian PAS (Fig. 1), positional information gain identified relevant regions in the proximity of most of the PAS variants (in the seven most common PAS variants). Conversely, positional information gain showed no clear segments in the rest of the less common variants, possibly indicating the weaker presence of cis-regulatory elements in such variants. Interestingly, the AATACA variant revealed the importance of the upstream region. These observations suggest that the polyadenylation mechanisms behind each of the PAS variants may be considerably different. With these points in mind, we proposed a new set of features along with Omni-PolyA, a novel model for PAS prediction implemented as an online tool. To derive a robust model for each of the PAS variants, Omni-PolyA consists of a set of different classification models organized in a tree-like structure. To evaluate the performance of our model, we derived an Omni-PolyA model by using the DPS feature set proposed in [31], showing that Omni-PolyA consistently outperformed reported results by DPS (by 26.85%, Table 1). Next, we showed the performance of the model using the novel Omni-PolyA feature set, which reduced the average error rate by 35.37 and 41.34% compared to DPS and HMM_SVM state-of-the-art tools, respectively. Notably, the prediction of PAS-strong variants showed the most significant improvement, reducing the error rate by 36.95% and 49.50% compared to DPS and HMM_SVM, respectively. Finally, we used the GENCODE annotation (release 19) to obtain the recent curated human poly(A) data for GRCh37. Results in Table 3 show that Omni-PolyA consistently reduced the weighted average error rate by 6.86% and 20.63% compared to DPS and HMM_SVM, with the largest error reduction for the PAS-weak variants (11.50% and 32.64% relative to DPS and HMM_SVM tools, respectively).

Methods

Datasets

We considered two different datasets to assess the performance of the Omni-PolyA method. A PAS sequence is considered to be a genomic DNA sequence of 206 nt in length (100 nt downstream and 100 nt upstream flanking a PAS hexamer). The first dataset, proposed by Kalkatawi et al. [31, 32], is considered as a benchmark to compare against the state-of-the-art methods. This dataset contains 14,470 PAS-like sequences (7370 sequences with true PAS motif and 7370 pseudo-PAS sequences). The pseudo-PAS sequences contain canonical PAS hexamers (from the 12 PAS variants we considered) but with no links to the polyadenylation process. For each PAS variant, the number of sequences with true PAS hexamer and pseudo-PAS hexamer is selected to be the same. The true PAS sequences were obtained by mapping human mRNA sequences to the human genome (hg19). The pseudo-PAS sequences were randomly selected from human chromosome 21 after excluding the true PAS sequences.

The second dataset considered is based on the experimentally validated GENCODE annotation and is extracted from the human genome (hg19). We used the GENCODE PolyA feature annotation Release 19 (GRCh37.p13) [44], which contains polyA features manually annotated by the HAVANA group (http://www.sanger.ac.uk/science/groups/vertebrate-annotation). We used the information in the annotation file (GTF) to extract true PAS from the genome. In total, 18,786 sequences with true PAS were extracted for the 12 most frequent PAS variants in human. For each PAS variant, the same number of pseudo-PAS sequences was generated from human chromosome 21 after excluding all the true PAS sequences contained in that chromosome.

Data normalization and cross-validation splits

We used the 5-fold cross-validation technique to validate the performance of all considered models. In the k-fold cross-validation, the original data is partitioned into k (approximately) equal-sized subsets. For each of the cross-validation folds, one of the subsets is used for testing the model while the remaining k-1 subsets are used to derive the classification model. Moreover, we reserved 15% of the training set for each fold as a validation set to optimize Omni-PolyA, DNN, DPS and, HMM_SVM model parameters. Therefore, the test set is exclusively used to assess the model performance in the final testing phase. Finally, feature values for all PAS variants were normalized to have values within the range of (−1, 1) according to

$$ {norm}_i=\frac{x_i-\left(\mathrm{ma}{\mathrm{x}}_i+\mathrm{mi}{\mathrm{n}}_i\right)/2}{\left(\mathrm{ma}{\mathrm{x}}_i-\mathrm{mi}{\mathrm{n}}_i\right)/2} $$

where max and min refer to the maximum and minimum values for i-th feature and x is the feature value that will be normalized. To avoid biased predictions, it is important to note that max and min values are obtained exclusively from the training data and are used as part of the model for the normalization of validation and test data (Fig. 6).

Fig. 6
figure 6

Data normalization procedure. Schematic representation of data normalization for fold 1 in a 5-fold cross-validation

Model training and data configuration

DPS model

We used DPS feature set to derive a random forest model as specified by Kalkatawi et al. [31]. Model parameters (i.e., the number of trees and number of randomly selected features) were determined based on the validation set (see Additional file 7: Table S4 for model parameters for each PAS variant).

HMM_SVM model

We derived HMM_SVM models by using the code provided by Xie et al. [33]. We optimized the model parameters (number of observations to combine into mega-state and the number of singular vectors to keep) by using a grid search method as specified by authors (see Additional file 7: Table S4 for model parameters for each PAS variant).

DNN model

We used MATLAB and Neural Network Toolbox release 2016a to derive DNN models with two autoencoder layers and one softmax layer. The number of units in the autoencoder layers was experimentally found by optimizing the error rate based on the validation set (see Additional file 7: Table S4 for model parameters for each PAS variant). DNN results in Table 1 and Table 3 show the performance of the models derived by using the DPS feature set (274 numeric features) [31].

Omni-PolyA model

Omni-PolyA uses four different classification models, namely, artificial neural network, random forest, C4.5 decision tree, and MLR models from WEKA v3.6.12 [61] and it is available as a MATLAB toolbox and as an online tool accessible at www.cbrc.kaust.edu.sa/omnipolya/. Moreover, columns sixth and seventh from Table 1 show the results obtained by the Omni-PolyA model derived by using the Omni-PolyA feature set (218 numeric features, listed in S1 Table). Finally, we pooled all PAS variants from the PAS-weak category to expand the training data. Therefore, the Omni-PolyA model was trained by using data from 10 PAS-weak variants and tested on the separate test set for a given variant (Fig. 7). Additional file 8: Table S5 shows the Omni-PolyA model parameters (determined from the validation set) for each PAS variant.

Fig. 7
figure 7

Schematic representation of Omni-PolyA model derived from pooled PAS-weak variants. The illustration shows the first fold of a 5-fold cross-validation technique for variant AAGAAA

Measures for assessing the model performance

As a representative measure of model performance, we used the classification error rate defined as

$$ errorrate= 1\hbox{-} \frac{TP+ TN}{TP+ TN+ FN+ FP}, $$

where TP, TN, FP and FN stand for the number of true positive predictions, true negative predictions, false positive predictions, and false negative predictions, respectively.

Model-derived features

Here we show a brief description of features that are derived from a model using a portion of the training data. These features were inspired by those used by Magana-Mora et al. [62, 63] in models aimed to recognize translation initiation sites in plants genomic DNA.

PAS-score and pseudo PAS-score from 2-mer weight matrices

Using frequencies of the 16 di-nucleotide combinations (AA, AT, AC, AG, etc.), we derived two 2-mer weight matrices from the training data to represent the characteristics of both PAS and pseudo-PAS sequences separately. Both 2-mer weight matrices are then used to calculate a score (see below) indicating the likelihood of the sequence to be a functional PAS and a pseudo-PAS. Consequently, a PAS score and pseudo-PAS score is computed for each DNA sequence in the dataset and are calculated as follows: let S(aj) be a DNA sequence of length L, and P(pij) be a 2-mer weight matrix of L-1 columns and 16 rows (for each di-nucleotide combination). The PAS score and pseudo-PAS score are given by

$$ \left[ PAS| pseudoPAS\right] score={\sum}_{i=1}^{16}{\sum}_{j=1}^{L-1}\mathrm{lo}{\mathrm{g}}_2\left(\frac{p_{ij}\times {a}_j{a}_{j+1}}{Pb_i}\right) $$

and

$$ {p}_{ij}\times {a}_j{a}_{j+1}=\left\{\begin{array}{cc}\hfill {p}_{ij},\hfill & \hfill {a}_j{a}_{j+1}={r}_i\hfill \\ {}\hfill 1,\hfill & \hfill {a}_j{a}_{j+1}\ne {r}_i\hfill \end{array}\right., $$

where Pb i refers to the background probability from a uniform distribution.

Positional information gain score

We calculated the information gain independently for each position of the genomic sequence surrounding PAS. For this, we first computed the entropy of each position as follows: for a given position P in a training sequence we calculate the entropy for a nucleotide X as:

$$ E\left(P,X\right)=-\frac{c_1}{c_1+{c}_2}\mathrm{lo}{\mathrm{g}}_2\frac{c_1}{c_1+{c}_2}-\frac{c_2}{c_2+{c}_1}\log \frac{c_2}{c_2+{c}_1}, $$

where c 1 represents the number of occurrences of nucleotide X (A, C, T or G) at position P in PAS sequences and c 2 represents the number occurrences of the same nucleotide at position P in pseudo-PAS sequences. We also introduce another entropy measure at position P that adjusts for the proportion of PAS and pseudo-PAS samples in the training set in the following way

$$ E(P)=-\frac{c_1}{c_1+{c}_2}\mathrm{lo}{\mathrm{g}}_2\frac{c_1}{c_1+{c}_2}-\frac{c_2}{c_2+{c}_1}\log \frac{c_2}{c_2+{c}_1}, $$

where c 1 and c 2 indicate the number of PAS and pseudo PAS sequences in the training set, respectively. Finally, we calculated the information gain for a position P as defined in Russel and Norvig [64]:

$$ Gain(P)=E(P)-E\left(A,P\right)-E\left(C,P\right)-E\left(G,P\right)-E\left(T,P\right), $$

where A, C, G, and T refer to the four nucleotides. Finally, the sum of information gain for each position in the entire sequence (information gain score) is then used as one single feature. Therefore, samples with high and low information gain scores suggest PAS or pseudo-PAS sequences, respectively.

Nucleotide frequency of the most discriminative positions

We used the positional information gain (described above) for selecting the 40 most discriminant positions (20 from the upstream and 20 for the downstream regions relative to the PAS). We then counted the frequency of A, C, G, and T nucleotides in the 20 selected positions in the downstream and upstream, separately. Consequently, this results in eight numeric features denoting the frequency of A, C, G, and T in the most discriminant downstream and upstream positions.

Abbreviations

DNN:

Deep neural networks

DPS:

Dragon PolyA Spotter tool

FN:

False negatives

FP:

False positives

HMM_SVM:

Tool based on hidden Markov model and support vector machine for poly(A) prediction

Nt:

Nucleotides

PAS:

Polyadenylation signal

SVM:

Support vector machine

TN:

True negatives

TP:

True positives

References

  1. Proudfoot NJ. Poly(A) signals. Cell. 1991;64(4):671–4.

    Article  CAS  PubMed  Google Scholar 

  2. Kung JTY, Colognori D, Lee JT. Long Noncoding RNAs: Past, Present, and Future. Genetics. 2013;193(3):651–69.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  3. Liu H, Han H, Li J, Wong L. An in-silico method for prediction of polyadenylation signals in human sequences. Genome Inform. 2003;14:84–93.

    CAS  PubMed  Google Scholar 

  4. Salamov AA, Solovyev VV. Recognition of 3 -processing sites of human mRNA precursors. Bioinformatics. 1997;13(1):23–8.

    Article  CAS  Google Scholar 

  5. Wilusz J, Shenk T, Takagaki Y, Manley JL. A multicomponent complex is required for the AAUAAA-dependent cross-linking of a 64-kilodalton protein to polyadenylation substrates. Mol Cell Biol. 1990;10(3):1244–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Colgan DF, Manley JL. Mechanism and regulation of mRNA polyadenylation. Genes Dev. 1997;11:2755–66.

    Article  CAS  PubMed  Google Scholar 

  7. Proudfoot NJ. Ending the message: poly(A) signals then and now. Genes Dev. 2011;25(17):1770–82.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Gil A, Proudfoot NJ. Position-dependent sequence elements downstream of AAUAAA are required for efficient rabbit beta-globin mRNA 3′ end formation. Cell. 1987;49(3):399–406.

    Article  CAS  PubMed  Google Scholar 

  9. McLauchlan J, Gaffney D, Whitton JL, Clements B. The consensus sequence YGTGTTYY located downstream from the AATAAA signal is required for efficient formation of mRNA 3′ termini. Nucleic Acids Res. 1985;13(4):1347–68.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Carswell S, Alwine JC. Efficiency of utilization of the simian virus 40 late polyadenylation site: effects of upstream sequences. Mol Cell Biol. 1989;9(10):4248–58.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Valsamakis A, Zeichner S, Carswell S, Alwine JC. The human immunodeficiency virus type 1 polyadenylylation signal: a 3′ long terminal repeat element upstream of the AAUAAA necessary for efficient polyadenylylation. Proc Natl Acad Sci. 1991;88:2108–12.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. DeZazzo JD, Kilpatrick JE, Imperiale MJ. Involvement of long terminal repeat U3 sequences overlapping the transcription control region in human immunodeficiency virus type 1 mRNA 3′ end formation. Mol Cell Biol. 1991;11(3):1624–30.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. Wahle E. 3′-end cleavage and polyadelanytion of mRNA precursor. Biochim Biophys Acta. 1995;1261(2):183–94.

    Article  PubMed  Google Scholar 

  14. Chen F, MacDonald CC, Wilusz J. Cleavage site determinants in the mammalian polyadenylation signal. Nucleic Acids Res. 1995;23(14):2614–20.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Tian B, Hu J, Zhang H, Lutz CS. A large-scale analysis of mRNA polyadenylation of human and mouse genes. Nucleic Acids Res. 2005;33(1):201–12.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. Edwalds-Gilbert G, Veraldi KL, Milcarek C. Alternative poly(A) site selection in complex transcription units: means to an end? Nucleic Acids Res. 1997;25(13):2547–61.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Tian B, Pan Z, Lee JY. Widespread mRNA polyadenylation events in introns indicate dynamic interplay between polyadenylation and splicing. Genome Res. 2007;17(2):156–65.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Lutz CS. Alternative polyadenylation: a twist on mRNA 3′ end formation. ACS Chem Biol. 2008;17(3):609–17.

    Article  Google Scholar 

  19. Yan J, Marr TG. Computational analysis of 3′-ends of ESTs shows four classes of alternative polyadenylation in human, mouse, and rat. Genome Res. 2005;15(3):369–75.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Frischmeyer PA, Av H, O'Donnell K, Guerrerio AL, Parker R, Dietz HC. An mRNA surveillance mechanism that eliminates transcripts lacking termination codons. Science. 2002;295(5563):2258–61.

    Article  CAS  PubMed  Google Scholar 

  21. Yada T, Ishikawa M, Totoki Y, Okubo K. Statistical analysis of human DNA sequences in the vicinity of poly(A) signal. Institute for New Generation Computer Technology Technical Report TR-876. 1994.

    Google Scholar 

  22. Kondrakhin YV, Shamin VV, Kolchanov NA. Construction of a generalized consensus matrix for recognition of vertebrate pre-mRNA 3 - terminal processing sites. Comput Applic Biosci. 1994;10:597–603.

    CAS  Google Scholar 

  23. Roberts RJ, O'Neill KE, Yen CT. DNA sequence from the adenovirus 2 genome. J Biol Chem. 1984;259(22):13968–75.

    CAS  PubMed  Google Scholar 

  24. Tabaska JE, Zhang MQ. Detection of polyadenylation signals in human DNA sequences. Gene. 1999;231(1–2):77–86.

    Article  CAS  PubMed  Google Scholar 

  25. Legendre M, Gautheret D. Sequence determinants in human polyadenylation site selection. BMC Genomics. 2003;4:7.

    Article  PubMed  PubMed Central  Google Scholar 

  26. Bajic VB, Charn TH, Xu JX, Panda SK, Krishnan SPT. Prediction Models for DNA Transcription Termination Based on SOM Networks. Conf Proc IEEE Eng Med Biol Soc. 2005;5:4791–4.

    CAS  PubMed  Google Scholar 

  27. Cheng Y, Miura RM, Bin T. Prediction of mRNA polyadenylation sites by support vector machine. Bioinformatics. 2006;22(19):2320–5.

    Article  CAS  PubMed  Google Scholar 

  28. Hu J, Lutz CS, Wilusz J, Tian B. Bioinformatic identification of candidate cis-regulatory elements involved in human mRNA polyadenylation. RNA. 2005;11(10):1485–93.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Xu J-X, Ashok B, Panda SK, Bajic V. Modeling Transcription Termination of Selected Gene Groups Using Support Vector Machine, IEEE International Joint Conference on Neural Networks (IJCNN). Hong Kong, China: IEEE; 2008. p. 384–9.

    Google Scholar 

  30. Akhtar MN, Bukhari SA, Fazal Z, Qamar R, Shahmuradov IA. POLYAR, a new computer program for prediction of poly(A) sites in human sequences. BMC Genomics. 2010;11

  31. Kalkatawi M, Rangkuti F, Schramm M, Jankovic BR, Kamau A, Chowdary R, et al. Dragon PolyA Spotter: predictor of poly(A) motifs within human genomic DNA sequences. Bioinformatics. 2013;29(11):1484.

    Article  CAS  PubMed  Google Scholar 

  32. Kalkatawi M, Rangkuti F, Schramm M, Jankovic BR, Kamau A, Chowdary R, et al. Dragon PolyA Spotter: predictor of poly(A) motifs within human genomic DNA sequences. Bioinformatics. 2012;28(1):127–9.

    Article  CAS  PubMed  Google Scholar 

  33. Xie B, Jankovic BR, Bajic VB, Song L, Gao X. Poly(A) motif prediction using spectral latent features from human DNA sequences. Bioinformatics. 2013;29(13):i316–25.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Mitchell M. An Introduction to Genetic Algorithms. Cambridge, MA: MIT Press; 1996.

    Google Scholar 

  35. Yildiz OT, Alpaydin E. Omnivariate decision trees. IEEE transactions on Neural Neutworks. 2001;12(6):1539–46.

    Article  CAS  Google Scholar 

  36. Yildiz OT. Model selection in omnivariate decision trees using Structural Risk Minimization. Inf Sci. 2011;181(23):5214–26.

    Article  Google Scholar 

  37. Yildiz OT, Alpaydin E. In: Gama J, Camacho R, Brazdil P, Jorge A, Torgo L, editors. Model Selection in Omnivariate Decision Trees. ECML: Springer; 2005. p. 473–84.

    Google Scholar 

  38. Magana-Mora A, Bajic VB. OmniGA: Optimized omnivariate decision trees for generalizable classification models. Scientific Reports. 2017;7(1):3898.

  39. LeCun Y, Bengio Y, Hinton G. Deep learning. Nature. 2015;521:436-44.

  40. Ciodaro T, Deva D, Seixas JM, Damazio D. Online particle detection with neural networks based on topological calorimetry information. J Phys Conf Ser. 2012;368:1.

    Article  Google Scholar 

  41. Ma J, Sheridan RP, Liaw A, Dahl GE, Svetnik V. Deep neural nets as a method for quantitative structure-activity relationships. J Chem Inf Model. 2015;55(2):263–74.

    Article  CAS  PubMed  Google Scholar 

  42. Forman G, Cohen I. Learning from Little: Comparison of Classifiers Given Little Training. In: J F B, F E, F G, D P, editors. Knowledge Discovery in Databases: PKDD 2004: Lecture Notes in Computer Science, vol 3202. Springer, Berlin, Heidelberg; 2004.

  43. Basavanhally A, Viswanath S, Madabhushi A. Predicting Classifier Performance with Limited Training Data: Applications to Computer-Aided Diagnosis in Breast and Prostate Cancer. PLoS One. 2015;10:15.

    Article  Google Scholar 

  44. Harrow J, Frankish A, Gonzalez J, Tapanari E, Diekhans M, Kokocinski F, et al. GENCODE: the reference human genome annotation for The ENCODE Project. Genome Res. 2012;22(9):1760–74.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  45. Olson WK, Gorin AA, Lu X-J, Hock LM, Zhurkin VB. DNA sequence-dependent deformability deduced from protein-DNA crystal complexes. Proc Natl Acad Sci. 1998;95:1163–8.

    Article  Google Scholar 

  46. Abeel T, Saeys Y, Bonnet E, Rouzé P. Peer YVd. Generic eukaryotic core promoter prediction using structural features of DNA. Genome Res. 2008;18(2):310–23.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  47. Florquin K, Saeys Y, Degroeve S, Rouzé P. Peer YVd. Large-scale structural analysis of the core promoter in mammalian and plant genomes. Nucleic Acids Res. 2005;33(13):4255–64.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  48. Ivanov VI, Minchenkova LE. The A-form of DNA: in search of the biological role. Mol Biol. 1994;28(6):1258–71.

    CAS  Google Scholar 

  49. DiMaio F, Yu X, Rensen E, Krupovic M, Prangishvili D, Egelman EH. A virus that infects a hyperthermophile encapsidates A-form DNA. Science. 2015;348(6237):914–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  50. Elhassan MA, Calladine CR. Propeller-twisting of base-pairs and the conformational mobility of dinucleotide steps in DNA. J Mol Biol. 1996;259(1):95–103.

    Article  CAS  Google Scholar 

  51. Brukner I, Sánchez R, Suck D, Pongor S. Sequence-dependent bending propensity of DNA as revealed by DNase I: parameters for trinucleotides. EMBO J. 1995;14(8):1812–8.

    CAS  PubMed  PubMed Central  Google Scholar 

  52. Sugimoto N, Nakano S, Yoneyama M, Honda K. Improved thermodynamic parameters and helix initiation factor to predict stability of DNA duplexes. Nucleic Acids Res. 1996;24(22):4501–5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  53. Sivolob AV, Khrapunov SN. Translational positioning of nucleosomes on DNA: the role of sequence-dependent isotropic DNA bending stiffness. J Mol Biol. 1995;247(5):918–31.

    Article  CAS  PubMed  Google Scholar 

  54. Ho PS, Zhou GW, Clark LB. Polarized electronic spectra of Z-DNA single crystals. Biopolymers. 1990;30(1):151–63.

    Article  CAS  PubMed  Google Scholar 

  55. Blake RD, Delcourt SG. Thermal stability of DNA. Nucleic Acids Res. 1998;26(14):3323–32.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  56. Blake RD, Bizzaro JW, Blake JD, Day GR, Delcourt SG, Knowles J, et al. Statistical mechanical simulation of polymeric DNA melting with MELTSIM. Bioinformatics. 1990;15(5):370–5.

    Article  Google Scholar 

  57. Satchwell SC, Drew HR, Travers AA. Sequence periodicities in chicken nucleosome core DNA. J Mol Biol. 1986;191(4):659–75.

    Article  CAS  PubMed  Google Scholar 

  58. Ornstein LR, Rein R, Breen DL, Macelroy RD. An optimized potential function for the calculation of nucleic acid interaction energies I. Base stacking. Biopolymers. 1978;17(10):2341–60.

    Article  CAS  PubMed  Google Scholar 

  59. Gorin AA, Zhurkin VB, Olson WK. B-DNA twisting correlates with base-pair morphology. Mol Biol. 1995;247(1):34–48.

    Article  CAS  Google Scholar 

  60. Breslauer KJ, Frank R, Blöcker H, Marky LA. Predicting DNA duplex stability from the base sequence. Proc Natl Acad Sci U S A. 1986;83(11):3746–50.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  61. Witten IH, Frank E, Trigg L, Hall M, Holmes G, Jo S. Weka: Practical Machine Learning Tools and Techniques with Java Implementations. 1999.

    Google Scholar 

  62. Magana-Mora A, Ashoor H, Jankovic BR, Kamau A, Awara K, Chowdary R, et al. Dragon TIS Spotter: an Arabidopsis-derived predictor of translation initiation sites in plants. Bioinformatics. 2013;29(1):117–8.

    Article  CAS  PubMed  Google Scholar 

  63. Ashoor H, Magana-Mora A, Jankovic BR, Kamau A, Awara K, Chowdary R, et al. Recognition of Translation Initiation Sites in Arabidopsis Thaliana. In: Lecca P, Tulpan D, Rajaraman K, editors. Systemic Approaches in Bioinformatics and Computational Systems Biology: Recent Advances: IGI Global; 2011. p. 105–16.

    Google Scholar 

  64. Russel S, Norvig P. Artificial Intelligence: A Modern Approach. second edition ed: Prentice-Hall; 2003.

    Google Scholar 

Download references

Acknowledgements

This research made use of the resources of CBRC and IT Research Computing at King Abdullah University of Science and Technology (KAUST), Thuwal, Saudi Arabia.

Funding

This work was supported by King Abdullah University of Science and Technology (KAUST) through the baseline fund BAS/1/1606–01-01 of VBB.

Availability of data and materials

The datasets generated and/or analyzed during the current study are publicly available online at www.cbrc.kaust.edu.sa/omnipolya/download.html or https://doi.org/10.5281/zenodo.321282.

Omni-PolyA code can be obtained from and accessed as an online tool at www.cbrc.kaust.edu.sa/omnipolya/index.html.

Author information

Authors and Affiliations

Authors

Contributions

AMM designed the study, analyzed the data, implemented the software and wrote the manuscript. MK performed the data curation, validation and, contributed to the discussion. VBB conceptualized and designed the study, contributed to the discussion, and reviewed/edited the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Vladimir B. Bajic.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional files

Additional file 1: Table S1.

Omni-PolyA feature set. List of the 218 numerical features. (PDF 104 kb)

Additional file 2: Table S2.

Comparison of performances achieved by DPS, HMM_SVM, DNN, and Omni-PolyA. (PDF 110 kb)

Additional file 3: Table S3.

False positive and false negative rates comparison between DPS, DNN, and Omni-polyA derived by using different feature sets. (PDF 105 kb)

Additional file 4: Figure S1.

Nucleotide distribution for PAS variants in the PAS-weak category. These plots show the frequency of nucleotides for true PAS sequences in the 10 variants from the PAS-weak category. (PDF 1696 kb)

Additional file 5: Figure S2.

DNA structural profiles of the PAS-strong variants. These plots represent the 16 considered structural profiles. Each structural profile is the average over all sequences from the PAS-strong variants (AATAAA and ATTAAA). These plots show the actual average values (y axis) over all sequences for each position (x axis). (PDF 2541 kb)

Additional file 6: Figure S3.

DNA structural profiles of the PAS-weak variants. Each structural profile is the average over all sequences from the PAS-weak variants (10 PAS variants). These plots show the actual average values (y axis) over all sequences for each position (x axis). (PDF 2412 kb)

Additional file 7: Table S4.

DPS, HMM_SVM and DNN model parameters. Parameters were determined from the validation set. (PDF 87 kb)

Additional file 8: Table S5.

Omni-PolyA model parameters. Genetic algorithm parameters and feature set configuration determined from the validation set. (PDF 89 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Magana-Mora, A., Kalkatawi, M. & Bajic, V.B. Omni-PolyA: a method and tool for accurate recognition of Poly(A) signals in human genomic DNA. BMC Genomics 18, 620 (2017). https://doi.org/10.1186/s12864-017-4033-7

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12864-017-4033-7

Keywords