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

An intuitionistic approach to scoring DNA sequences against transcription factor binding site motifs

Abstract

Background

Transcription factors (TFs) control transcription by binding to specific regions of DNA called transcription factor binding sites (TFBSs). The identification of TFBSs is a crucial problem in computational biology and includes the subtask of predicting the location of known TFBS motifs in a given DNA sequence. It has previously been shown that, when scoring matches to known TFBS motifs, interdependencies between positions within a motif should be taken into account. However, this remains a challenging task owing to the fact that sequences similar to those of known TFBSs can occur by chance with a relatively high frequency. Here we present a new method for matching sequences to TFBS motifs based on intuitionistic fuzzy sets (IFS) theory, an approach that has been shown to be particularly appropriate for tackling problems that embody a high degree of uncertainty.

Results

We propose SC intuit , a new scoring method for measuring sequence-motif affinity based on IFS theory. Unlike existing methods that consider dependencies between positions, SC intuit is designed to prevent overestimation of less conserved positions of TFBSs. For a given pair of bases, SC intuit is computed not only as a function of their combined probability of occurrence, but also taking into account the individual importance of each single base at its corresponding position. We used SC intuit to identify known TFBSs in DNA sequences. Our method provides excellent results when dealing with both synthetic and real data, outperforming the sensitivity and the specificity of two existing methods in all the experiments we performed.

Conclusions

The results show that SC intuit improves the prediction quality for TFs of the existing approaches without compromising sensitivity. In addition, we show how SC intuit can be successfully applied to real research problems. In this study the reliability of the IFS theory for motif discovery tasks is proven.

Background

Cells control the abundance of proteins by means of diverse mechanisms. One such mechanism is the regulation of transcription, which is a continuous process whereby many factors combine to ensure appropriate rates of protein synthesis. Understanding such complex processes is one of the main objectives in computational biology. In its early stages, transcription is controlled, among other mechanisms, by the binding of proteins called transcription factors (TFs) to specific regions of a given chromosome called transcription factor binding sites (TFBSs). These interactions between proteins and DNA usually take place upstream from the gene, close to the transcription start site (TSS), in the so-called promoter region of the gene.

One of the biggest issues in identifying TFBSs is that a single binding protein can bind to different DNA sequences. Related DNA sequences to which the same TF can bind are grouped together into a TFBS motif. The identification of TFBSs within a given set of DNA sequences is an active area of research. In this context there exist two main approaches: i) the de novo discovery of motifs, and ii) the detection of TFBSs using motifs that are already known.

De novo methods aim to find significant sub-sequence patterns within a set of TFBS sequences. Some of the most popular approaches are MEME [1], Gibbs sampling [2], AlignACE [3], and more recently PRIORITY and Trawler [4, 5]. For a review see [6].

Detection methods, on the other hand, focus on inferring new TFBSs from known binding motifs. Early detection methods assumed independence between positions within a putative TFBS sequence, e.g. in Patser [7] and ConSite [8]. However, it is now well established that this assumption is wrong [9–11], and some methods that consider position dependency for modeling and finding TFBSs using advanced HMM and Bayesian models have appeared [12, 13]. Likewise, two recent detection methods have been developed that take into account interdependencies between TFBS positions. Tomovic and Oakeley proposed a method that incorporates a measure of positional interdependence into the overall score [14]. More recently, Zare-Mirakabad et al. developed a method based on joint information content and mutual information [15]. In this method, positional dependencies are taken into account by considering all pairwise combinations of positions (see the Methods section for more information).

The fact that TFBS sequences are usually very short means that the same or very similar sequences tend to occur by chance at a relatively high frequency. Consequently one of the main goals in the prediction of TFBSs is to reduce the false positive rate without compromising sensitivity. Methods that take into account positional dependencies tend to be significantly more effective at meeting this challenge. However, there remains room for improvement. As we will show in the Results section, existing methods have some drawbacks, such as overlearning of the training data, arbitrary threshold selection for testing dependencies, etc. The purpose of the work presented here is to provide a new method for measuring sequence-motif affinity that improves on existing approaches.

Zadeh proposed fuzzy sets theory to mathematically model the imprecision inherent in certains concepts [16]. Briefly, fuzzy sets theory allows an object to partially belong to a set with a membership degree between 0 and 1. Classical set theory is a special case of its fuzzy counterpart in which membership and certainty degrees are restricted to either 0 or 1. Atanassov proposed intuitionistic fuzzy sets (IFS) theory as an extension of the fuzzy sets theory [17]. IFSs generalize the notion of a fuzzy set representing uncertainty with respect to both the degree of membership (μ) and non-membership (ν) of a set by allowing that the sum μ + ν ≤ 1.

Owing to the fact that IFSs are capable of modelling the uncertainty present in real-life situations, they have been widely applied during the past decades to a variety of problems (see the Methods section). In recent years, it has been seen that the inherent uncertainty and noise that characterize biological data cannot always be modeled sufficiently well using probabilistic approaches and that, as a consequence, alternative approaches to modelling this uncertainty may be required [18–21]. In addition to the usual problems of missing values and noisy data associated with biological data, there exist some additional hidden factors that affect binding affinities in the context of sequence-motif scoring, e.g. cooperative binding and chromatin structure [22]. Furthermore, the described motifs are subject to change as new experiments confirm new binding sites. In this work we make use of IFS theory to formally model the uncertainty associated with the problem of scoring DNA sequences against TFBS motifs.

Results

Case studies

First, we wanted to show the ability of our proposed method, SC intuit , to discriminate between the relative importance of poorly-conserved positions and well-conserved positions comparing it with the most representative scoring methods: i) SC indep , a probabilistic method that assumes positional independence; ii) SC dep , a scoring method proposed by Tomovic and Oakeley that take into account statistical interdependencies between TFBS positions [14]; and iii) SC mat , a scoring function proposed by Zare-Mirakabad et al. based on the dependency between all pairwise combinations of binding site positions [15].

In Figure 1(A) we show the binding sequences of the motif Dof3 found in the JASPAR database. It can be observed how the first position is highly conserved while the fifth position is poorly conserved. We then considered two sequences (Figure 1(B)): i) a sequence with a mismatch in the conserved position; and ii) a sequence with a mismatch in the poorly-conserved position. As has been explained above, it would be desirable that the score obtained for the case of the mismatch in the conserved position is lower than the scoring for the other sequence, as it shares the similarities in the most conserved positions of the motif. In Table 1 we show the results for the three methods. It can be observed that our proposed scoring method discriminates between the two cases, while the other three approaches provide almost the same score for both sequences, missing the difference between the conservation level of the positions being compared.

Figure 1
figure 1

Interdependencies between TFBS positions. Motif Dof3. A) Shows the binding sequences found in JASPAR. B) Shows the logo representation of the motif (center), with the sequence having the mismatch in the conserved position (above), and the sequence having the mismatch in the poorly-conserved position (below).

Table 1 Scoring results

In the majority of cases, the sequences known to belong to a given TFBS motif have very similar nucleotide compositions and highly conserved positions. However, in the databases of known motifs there are a number of examples where individual sequences differ from the majority in highly-conserved positions. Such a binding sequence can be considered an outlier with respect to the motif, i.e. a binding site that is not closely related to the other binding sites in the motif. When scoring new sequences against a given TFBS motif, we should generally tolerate small, additional variations in the sequence with respect to non-outliers, but be far less tolerant of mutations to outlier sequences. Here we evaluate the extent to which each scoring method is able to discriminate between sequences belonging to these two categories. Take, as a preliminary example, the binding sequences for motif MZF1 in the JASPAR database, as shown in Figure 2(A). It can be observed how the highlighted outlier sequence GGAGGA does not contain the higly-conserved base G at the third position, while the highlighted sequence TGGGGA is clearly a non-outlier (see motif logo in Figure 2(B)). We selected the highlighted sequences and created two new sequences by mutating its sixth position giving GGAGGG (derived from an outlier) and TGGGGT (derived from a non-outlier). In order to observe the discrimination degrees of the different scoring methods, we scored each sequence against the motif by means of the different methods.

Figure 2
figure 2

Ratio of scoring values. Motif MZF1. A) Shows the binding sequences found in JASPAR. B) Shows the logo representation of the motif (center), with the mutated sequences for the non-outlier binding sequence (above), and the outlier binding sequence (below). C) Shows the ratio of the scoring results for each of the three methods when scoring the mutated sequence of the non-outlier binding sequence and the outlier binding sequence respectively.

In reality, it would be desirable that the scoring for the case of the mutated outlier sequence be lower than the scoring for the mutated non-oulier sequence. Results obtained by the SC mat , SC dep , and SC indep methods failed to capture the expected differences, giving the incorrect impression that binding is likely to occur. On the other hand, our proposed method obtained a more realistic distance between the sequences, providing a much lower score for the mutated outlier sequence (Figure 2(C)).

These insights are confirmed in the following sections where the experiments are extended to use large datasets, and the results are measured in terms of discovery rates.

Prediction of TFBSs

Synthetic sequences

In order to compare the performance of the different methods in predicting TFBSs, we used the non-redundant publicly available JASPAR motifs database for our experiments [23]. We selected all motifs for which binding sequences are available (not only matrix profiles), resulting in a dataset of 124 motifs. For each of these motifs, a random number between 2 and 6 binding sites were randomly selected and inserted in random sequences of a random length between 200 bp and 500 bp from a third-order Markov model background distribution obtained from the RSAT (Regulatory Sequence Analysis Tools) [24]. For each position of each sequence we computed the score for their corresponding motifs with an assumed known TFBS length (the length of the inserted motif).

Usually, methods have a high sensitivity (i.e. can detect true positives), so that the key difference between them is the number of false positives. Although our ultimate aim is not to rely on essentially arbitrary thresholds to assess performance, we began our analysis by following the recommendations of Tomovic and Oakeley in [14], selecting thresholds of 0.7 and 0.8 indicating a correct classification for a binding site. Table 2 shows the precision (TP/(TP + FP )) of the different methods. In the additional file 1: "Synthetic sequences experiment" we show thresholded results for the different methods. These indicate that our proposed scoring function performed best, giving the smallest number of false positives per TF whilst simultaneously giving a high number of true positives.

Table 2 Synthetic sequences precision.

In order not to rely on the selection of an arbitrary threshold for evaluating the results, we computed a precision-recall (PR) curve for each considered method. PR curves are commonly used in information retrieval for evaluating classification performance and give a more informative picture of a method's performance than ROC (Receiver Operating Characteristic) curves [25] when dealing with highly skewed datasets as is the case here [26]. Figure 3 shows the PR graphs. SC intuit produces a better PR graph than the remaining methods (see Table 3 for AUC values). In addition, ROC curves can be found in the the additional file 2: "ROC curves".

Figure 3
figure 3

Synthetic sequences results. PR curves for the three scoring methods. SC intuit provides a more consistent classification than the rest of the methods.

Table 3 AUC values for the synthetic and mutated sequence experiments.

Mutated sequences

To further evaluate our proposed method, we obtained a set of putative binding sites that are very similar to those that are already known. This is a common scenario in motif discovery, where the set of known sequences belonging to a given binding motif is incomplete. In order to simulate this situation, we proceeded in a similar way to our previous experiment; all the steps were the same except that we gave a single base mutation at a random position within the selected binding site for each motif. PR curves and AUC values were computed to compare the performance of the different methods (Figure 4 and Table 3). ROC curves for this experiment can be obtained from the additional file 2: "ROC curves". The ROC and precision-recall graphs shows how SC intuit gives consistently superior values, with a higher AUC value (Table 3). It can be observed that the improvement of the performance of our method compared to SC indep , SC mat and SC dep grew with respect to the synthetic sequences experiment discussed in the previous section.

Figure 4
figure 4

Mutated sequences results. PR curves for the three scoring methods. Again, SC intuit provides a better classification than the rest of the methods.

Real Data

We analyzed the performance of the proposed methods when dealing with real experimental data. In order to do so, we made use of the published ChIP-seq data on binding of TFs in embryonic stem cells from mouse by Chen et al. in [27], as provided in the supplementary material of [28]. We considered the three TFs (SMAD1, c-Myc, and STAT3) that have binding sequences available in the TRANSFAC database [29]. Thus, we obtained three sets of 200 bp sequence segments centered at TF binding locations, and we randomly selected 50 sequence segments from each set for our study (see additional file 3: "FASTA sequences"). We scanned each set of sequences using the 124 TFs from JASPAR for which binding sequences are available. The results demonstrate the superior performance of our new scoring method, as it gives the smallest number of false positives per nucleotide and per TF (Figure 5), and maintains an excellent true-positive rate (Table 4). Detailed results can be found in the additional file 4: "Motif statistics". It can be seen that our method presents consistently low false-positive rates with all three sets of sequences, whereas the performance of the other methods is much more variable.

Figure 5
figure 5

Real data results. Average false-positive ratio per TF for different thresholds for the proposed scoring methods.

Table 4 True positive rate for the real data experiment

Study of Single Nucleotide Polymorphisms in TNFR1 Gene for the Response against Aspergillus Fumigatus

Hematological patients are typically treated by chemotherapy and/or radiation. These treatments usually produce immunosuppression and severe neutropenia. This clinical situation can be exploited by opportunistic pathogens such as Aspergillus fumigatus to cause a deadly infection called Invasive Pulmonary Aspergillosis (IPA) [30, 31]. The importance of finding ways to combat this pathogen is evidenced by the fact that IPA occurs in roughly 10% to 40% of hematological patients, with overall mortality rates ranging from 50% to 90% [32, 33].

Tumor necrosis factor (TNF) activates T lymphocytes in response to fungal infections through TNF receptors. One of the most important TNF receptors is TNFR1, which plays a crucial role in immune regulation and host immune responses. Experimental studies with TNFR1 knockout mice indicate that TNFR1 is indispensable in host resistance against several infections [34]. Our hypothesis is that single nucleotide polymorphisms (SNPs) in the TNFR1 gene may influence the innate immune response against Aspergillus fumigatus.

The gene encoding TNFR1 contain numerous polymorphisms [35, 36]. By means of different experiments, we concluded that TNFR1-609(G/T) polymorphism is critical in the development of the response against Aspergillus because it might be regulating the cell-mediated Th1 immune response. Details on these experiments are out of the scope of this work and can be consulted in [37]. In this section, we use our proposed scoring method SC intuit to investigate whether the TNFR1-609(G/T) promoter polymorphism is involved in the disruption of the recognition of a potential binding site for a critical transcription factor that could influence TNFR1 transcription level.

TNFR1-609(G/T) Polymorphism Binding Affinity

For this experiment we used TFBS motifs found in TRANSFAC database [29], which has been widely used in research involving regulatory elements [38]. In order to find interesting dependencies between the TNFR1-609(G/T) SNP and TFs binding affinity we scored the human TRANSFAC TFBSs against the TNFR1-609(G/T) polymorphism by means of the SC intuit method.

TFs bind to short parts of the TNFR1 promoter region and, therefore, for each trial, we need to define a fragment of the promoter sequence containing the TNFR1-609(G/T) SNP that might be considered as the putative TFBS. To this end, we need to determine the length of the sub-sequences and the relative offset to the position of the TNFR1-609(G/T) SNP. For each of the 446 human TFs in TRANSFAC, we generated a set of putative binding sequences by using a window size of a fixed length equal to the number of position of the corresponding TF. Moving the window across the sequence in 5'-3' direction gave us the sub-sequences for the TNFR1-609(G/T) SNP that we considered to be putative TFBSs (see Figure 6 for an example). Next, we scored each pair of sub-sequences (one sub-sequence for the G allele, and one for the T allele) against the given TF applying the SC intuit method.

Figure 6
figure 6

Putative binding sequences for the G allele of the TNFR1 -609(G/T) . Best results were found for the highlighted SG-8sequence. A corresponding set was also obtained for the T allele.

We were interested in those sub-sequences that fulfil two properties: i) they have a high score in one allele (G or T) so they can be considered as candidates to be binding sites, and ii) the score is substantially lower when considering the remaining allele so the SNP may affect to the binding affinity. For our current research we chose a conservative cut-off of 0.7, and retained TFBSs with a score above this threshold for further analysis. The scores for the selected TFBSs with respect to their corresponding alleles (sequences with a G(T) instead of a T(G) at position -609) are shown in Table 5. Subsequently we will discuss these findings from a biological perspective and show that the most interesting insights arise in the context of the ICSBP TF, which represents the highest scoring of all the human TRANSFAC motifs.

Table 5 SC intuit scores for the two alleles.

Functional Effect of ICSBP/IRF-8 in the TNFR1-609(C/T) SNP

In the previous section, we obtained predictive results using our SC intuit scoring method and TRANSFAC database (Table 5). From them, we selected four candidates according to the two properties outlined in the previous section, i.e E2A, HNF4, ICSBP, and Pax-2. We did not find described relations between IPA response for any of E2A, HNF4, and Pax-2 TFs. Logos for these TFs are provided in Figure 7.

Figure 7
figure 7

Discarded TFs. No described relations with IPA response were found.

On the other hand, we found ICSBP (also known as IRF-8) to be directly related with the purpose of our study. ICSBP/IRF-8 shows a preference for binding the T allele (see Table 5). As a member of IRF family of transcription factors it is an important modulator of IFNγ signalling cascade and was identified in association on the promoter region of numerous macrophage essential genes such as IL12, IL1β, IL18, iNOS or ISG15 [39].

In addition, several genes regulated by ICSBP/IRF-8, such as MAP4K4, IL-17R, and SOCS7, are involved in different stages of the nuclear factor κ B (NFκ B) signaling pathway [39]. Therefore, we can hypothesize that ICSBP/IRF-8 transcription factor might be also regulating the NFκ B signaling pathway through the control of the first gene of this signalling cascade, the TNFR1 gene. In support of this hypothesis, Zhao et al. established that ICSBP/IRF-8 and TNFR1 are closely related genes [40]. They found ICSBP/IRF-8 to be associated with an enhanced ubiquination of TNFR associated factor 6 (TRAF6), a protein that mediate the signal transduction from members of the TNF receptor superfamily, and the activation of AP-1 and NFκ B transcription factors.

On the other hand, several studies demonstrated that ICSBP/IRF-8 promotes the differentiation and activation of dendritic cells and macrophages cells [41, 42], and that, at the same time, TNFR1 mRNA level is increased during this biological process [43].

Taken into account these observations, we hypothesize that the presence of TNFR1-609(G/T) promoter polymorphisms can modify the binding affinity to ICSBP/IRF-8 (see Figure 8) and, therefore, it could be used to predict susceptibility to infection and to facilitate risk stratification of hematological patients. However, the question of whether the TNFR1 polymorphisms have biological relevance regulating mRNA TNFR1 levels through ICSBP/IRF-8 transcription factor remains unanswered. Functional analysis should be performed to demonstrate the role of TNFR1-609(G/T) polymorphism mediating the binding of ICSBP/IRF-8 to TNFR1 promoter.

Figure 8
figure 8

ICSBP against TNFR1 -609(G/T) polymorphism. It can be seen how the fourth position in the ICSBP motif (middle) matches the T allele (bottom) and mismatches the G allele (top).

Discussion

We have introduced a new IFS-based approach for scoring DNA sequences against DNA motifs called SC intuit . In this work we review three scoring schemes. These approaches have several drawbacks. SC indep is based on an incorrect assumption that the nucleotides of a given TFBS are independent. In that context, SC dep extended the score in order to account for positional dependencies. The problems associated with unnormalized scores at each position have been pointed out [15]. In addition, the results vary depending on the choice of the method and parameters for testing the dependencies. The main drawback with SC mat is that it has a tendency to overlearn the training data and consequently its performance decreases when applied to real problems. There is therefore a need for a scoring method that accounts for positional dependencies without compromising either the consistency or the accuracy of the results.

As explained above, SC intuit is based on the IFS theory, which has been successfully applied to problems that suffers from noisy and imprecise data. IFS theory represents uncertainty with respect to both the degree of membership and non-membership. The uncertainty associated with the tasks of scoring DNA sequences against motifs makes intuitionistic concepts particularly suitable for handling this kind of data. Taking advantage of such properties, we define the membership and non-membership degrees of a given pair bases at a given position not only as a function of their combined probability of occurrence, but also taking into account the importance of each individual base at its corresponding position.

One of the biggest issues for this kind of scoring methods is giving high scores for the known binding sequences of the motifs without overfitting. Our proposed approach adequately solves the problem of computing the score of a given sequence against a given motif by considering the binding sequences that comprise the motif not only individually but also as part of such set of sequences. Simple experiments shows how other methods fail in capturing realistic differences, while SC intuit provides good results (Figure 1, 2). Our method assigned high scores for known binding sites, disfavouring mutations in the conserved positions of the binding site.

These insights are confirmed from experiments for predicting TFBSs in large datasets. We compared the performance of the proposed scoring methods on recognizing motifs in sets of random sequences from a third-order Markov model background distribution in two circumstances: i) when inserting known binding sequences, and ii) when inserting mutated binding sequences. In both situations we found that our proposed method gave the smallest number of false positives per TF whilst simultaneously giving a high number of true positives (Figures 3, 4). More importantly, our method outperforms the other approaches when dealing with real experimental data derived from Chip-seq assays. In this case, again, the number of false positive is significantly reduced (Figure 5). Finally, we validated our method studying the effect of known SNPs of TNFR1 gene in the binding affinity of TRANSFAC TFs for the response against Aspergillus fumigatus. We found the highest scoring for ICSBP TF among all the human TRANSFAC motifs. Although functional analysis should be performed, according to several previous studies, we hypothesize that the presence of TNFR1-609(G/T)polymorphisms could be used to predict susceptibility to infection of hematological patients.

In general, the obtained results on the different experiments demonstrated that the proposed intuitionistic approach provide a better and more accurate model for the detection of motifs and for the relationships between positions of the TFBSs.

Conclusions

In the present study, we have introduced SC intuit , a new scoring method for measuring sequence-motif affinity, based on IFS theory. Our main objective was to improve the prediction quality for TFs of the existing approaches, reducing the false positive rate without compromising sensitivity. We show that SC intuit outperforms other approaches in motif recognition tasks, and prove how it can be successfully applied to real research problems. We have used our approach as a scanning method for the prediction of TFBSs, but it also can be incorporated with methods for de novo discovery of motifs. As intuitionistic theory is specially suitable for problems that deal with imprecise concepts, we are currently working on a fuzzy approach that applies the proposed scoring in an ab initio method to find motifs in large sets of related DNA sequences.

Methods

Alternative approaches

In recent years, several scoring methods for the prediction of TFBSs have been proposed. In this section we give a brief overview of those methods that take account of positional dependencies, as they have been shown to outperform methods that assume independence. Let us first introduce the notation. Let B = {A, C, G, T} be the set of the four DNA nucleotides. Let D be a set of ordered DNA sequences on B of length n. Let us suppose that we have a motif M = S1,..., S t , where S i is a DNA sequence on D consisting of t aligned binding sites of length n. The problem is then reduced to assigning a score to the pair formed by a given putative TFBS, S ∈ D, and a given motif, M.

In what follows we will follow the notation proposed by Wasserman and Sandelin in [44], where F (b, i), for b ∈ B and 1 ≤ i ≤ n shows the occurrences of nucleotide b in position i, and P ( b , i ) = F ( b , i ) t + a ( b ) , for b∈ B and 1 ≤ i ≤ n is the corrected probability of base b at position i, where a(b) is a smoothing parameter (a(b) = 0.001). a(b) = 0.01 is usually reported but our experiments show that smaller values provide more accurate results.

Statistical dependencies

Tomovic and Oakeley extended the previous method that assumed positional independence [14]. The authors also followed the notation of [44] and defined Wb, ias a position weighted matrix (PWM) of base b in position i computed as:

W b , i = l o g 2 P ( b , i ) P ( b ) ,
(1)

where P(b) is the background probability of base b. In the case where independence is assumed, the score for a given DNA sequence S can be computed by summing all the values of Wb, ifor every base in S:

S C i n d e p ( s ) = ∑ i = 1 n W s i , i .
(2)

The first step for extending this score involves testing the dependencies between each pair of positions i and j. The authors introduced three different methods: i) χ2 test; ii) G statistics; and iii) Bayesian hypothesis testing. The authors used these three methods to calculate the dependencies between pairs of positions in the motifs available in the public database JASPAR [23]. The reader should note that the accurate computation of positional dependencies is still an open problem since different results are obtained depending on the method and parameters used in their computation (see Supplementary Material 2-4 in [14]). Further details about obtaining the position dependencies and multiple test corrections can be found in [14].

In order to compute the new score, the corrected probability for the bases b1b2 ... b m in the dependent positions i1i2 ... i m is defined by:

P ( b 1 , … , b m , i 1 , … , i m ) = F ( b 1 , … , b m , i 1 , … , i m ) t + + a ( b 1 , … , b m ) ,
(3)

where a(b1,..., b m ) = a(b1) ... a(b m ) is a smoothing parameter.

It is straightforward then to obtain values that correspond to the PWM values:

W b 1 , … , b m , i 1 , … , i m = log 2 ( P ( b 1 , … , b m , i 1 , … , i m ) P ( b 1 ) … P ( b m ) ) .
(4)

Finally, their proposed scoring function, which incorporates positional dependencies, can be computed as:

S C d e p ( S ) = ∑ i = 1 k 1 W S i , i + ∑ i = 1 k 2 W S j i , S j i + 1 , j i , j i + 1 + ⋯ + + ∑ i = 1 k m W S j i , … , S j i + m − 1 , j i , … , j i + m − 1 ,
(5)

where, k1 is the number of independent positions, k2 is the number of dependent positions of order 2 (nucleotides at positions j i and ji+1) and k m the number of dependent positions of order m (nucleotides at positions j i , ji+1,..., ji+m-1).

For both the SC indep and SC dep it is advisable to perform the following normalization:

N S C = S C − min ( S C ) max ( S C ) − min ( S C )
(6)

Matrix based

Zare-Mirakabad et al. proposed a new scoring function based on the dependency between all pairwise combinations of binding site positions [15]. Their method is based on the mutual information matrix, defined as: (see equation (7))

M i j = ∑ b i , b j P ( b i , b j , i , j ) log 2 ( P ( b i , b j , i , j ) P ( b i , i ) P ( b j , j ) ) ,
(7)

and on the joint information content (JIC), defined as:

J I C = ∑ i = 1 n − 1 ∑ j = i + 1 n ∑ b 1 ∈ B ∑ b 2 ∈ B P ( b 1 , b 2 , i , j ) log ( P ( b 1 , b 2 , i , j ) P ( b 1 ) P ( b 2 ) ) .
(8)

In order to compute their score, the authors defined a PWM, WPW , containing 16 rows and (n · (n - 1)/2) columns for all the pairwise combinations of the positions:

W b 1 , b 2 , i , j P W = log ( P ( b 1 , b 2 , i , j ) P ( b 1 ) P ( b 2 ) ) + log ( P ( b 1 , b 2 , i , j ) P ( b 1 , i ) P ( b 2 , j ) ) ,
(9)

where b1, b2 ∈ 2 B and 1 ≤ i, j ≤ n and i ≠ j. For more on this method see [15].

Finally, for a given DNA sequence S ∈ D of length n the score SC mat is computed as:

S C m a t = ∑ i = 1 n − 1 ∑ j = i + 1 n W S i , S j , i , j P W .
(10)

In order to obtain a normalized value for the score, equation (6) should be applied.

Intuitionistic fuzzy sets

Intuitionistic fuzzy sets (IFS) theory was proposed by Atanassov [17]. It has been applied in such diverse fields as decision making [45], logic programming [46] medical diagnosis [47, 48], pattern recognition [49], etc. IFS theory is an extension of the fuzzy sets theory previously proposed by Zadeh [16] that allows the degrees of membership and non-membership to be independently uncertain, which makes the representation more flexible at capturing the current state of our understanding given inconclusive data [50, 51]. Next, we introduce some basic IFS concepts.

Let X be the universe of discourse. An intuitionistic fuzzy set A in X is an object having the form:

A = { ( x , μ A ( x ) , ν A ( x ) ) : x ∈ X } ,
(11)

where μ A , ν A : X → 0[1] denote membership function and non-membership function of A, satisfying 0 ≤ μ A + ν A ≤ 1 for every x ∈ X. Therefore, the degree of uncertainty of x to A is π A (x) = 1 - μ A - ν A . For more on this topic please refer to [17, 50, 51].

Intuitionistic representation of motifs

For our approach, a given motif M is represented as the set of IFSs of all the pairwise combinations of its positions: I M = { I i , j M } , where 1 ≤ i, j ≤ n and i ≠ j. Each of the i, j combinations for the motif positions is then an IFS of 16 elements defined as:

I i , j M = { b , μ I i , j M ( b ) , v I i , j M ( b ) : b ∈ B × B } ,
(12)

where B × B is the universe of discourse, i.e. the set of all 16 possible combinations of bases for two given positions i and j (AA, AC,..., TT).

Membership degree computation

μ I l , j M

represents the degree of membership of the pairs for the basis b1, b2 ∈ B in a given pair of positions i, j in a motif M. It can be automatically computed as:

μ I i , j M ( b 1 , b 2 ) = P ( b 1 , b 2 , i , j ) + ( 1 − P ( b 1 , b 2 , i , j ) ) P ( b 1 , i ) + P ( b 2 , j ) 2 ,
(13)

where the above notation holds. As can be seen, the membership degree is a function of the probability of the pair of bases being compared and their individual conservation. Obviously, 0 ≤ μ I i , j M ( b 1 , b 2 ) ≤ 1 and the degree increases as do the corrected probabilities of bases b1 and b2 in positions i and j, as well as the individual corrected probabilities P (b1, i) and P (b2, j).

Non-membership degree computation

ν I i , j M

represents the non-membership degree of the pairs for the basis b1, b2 ∈ B in a given pair of positions i and j in a motif M. It can be automatically computed as:

ν I i , j M ( b 1 , b 2 ) = ( I C i b 1 + I C j b 2 2 ) ( 1 − μ I i , j M ( b 1 , b 2 ) ) ,
(14)

where I C p b = 2 + P ( b , p ) log 2 ( P ( b , p ) ) 2 is the normalized information content of base b in position p and ν I i , j M ( b 1 , b 2 ) is in the range 0 ≤ ν I i , j M ( b 1 , b 2 ) ≤ 1 . The reader should note that we make use of the IC of the two positions i, j to assign the allowed degree of uncertainty for such pair of positions. When the IC is high, the degree of uncertainty will be low and viceversa. Likewise, as the information content of the two basis increases, the sum of the membership degrees, μ I i , j M ( b 1 , b 2 ) + ν I i , j M ( b 1 , b 2 ) , gets closer to 1. Therefore, the non-membership degree in positions i, j is a function of the corresponding membership degree and the uncertainty level computed for b1 and b2. It is easy to prove that μ I i , j M ( b 1 , b 2 ) + ν I i , j M ( b 1 , b 2 ) ≤ 1

Scoring

In order to define our proposed score, we first introduce the simplest case of scoring a length-2 DNA subsequence D = b1, b2 in the positions i and j of a motif M:

S C i n t u i t i , j ( b 1 , b 2 ) = μ I i , j M ( b 1 , b 2 ) ( max ( ν I i , j M ) − ν I i , j M ( b 1 , b 2 ) ) ,
(15)

where max ( ν I i , j M ) is the maximum degree of non-membership in M found in the pair of positions i and j considering all the possible combination of basis b1, b2 ∈ B2, and μ I i , j M ( b 1 , b 2 ) and ν I i , j M ( b 1 , b 2 ) are the membership degree and non-membership degree of the pairs for the basis (b1, b2) ∈ B in the pair of positions i, j of M, computed as stated in sections and respectively.

As with the previously defined scores, a normalization step needs to be performed in order to obtain comparable results. The source code can obtained from http://genome.ugr.es/intuit.

N S C i n t u i t i , j ( b 1 , b 2 ) = S C i n t u i t i , j ( b 1 , b 2 ) − min ( S C i n t u i t i , j ) max ( S C i n t u i t i , j ) − min ( S C i n t u i t i , j ) ,
(16)

where are the min ( S C i n t u i t i , j ) and max ( S C i n t u i t i , j ) are the min/max possible scores in the positions (i, j) of the motif.

Finally, for a given DNA sequence S ∈ D of length n the score SC intuit is computed as:

S C i n t u i t = ∑ i = 1 n − 1 ∑ j = i + 1 n N S C i n t u i t i , j ( S i , S j ) .
(17)

References

  1. Bailey T, Elkan C: Fitting a mixture model by expectation maximization to discover motifs in biopolymers. Proc Int Conf Intell Syst Mol Biol 1994, 2: 28–36. Citeseer Citeseer

    CAS  PubMed  Google Scholar 

  2. Lawrence C, Altschul S, Boguski M, Liu J, Neuwald A, Wootton J: Detecting subtle sequence signals: a Gibbs sampling strategy for multiple alignment. Chem Rev 1993, 93: 741. 10.1021/cr00018a006

    Article  Google Scholar 

  3. Hughes J, Estep P, Tavazoie S, Church G: Computational identification of cis-regulatory elements associated with groups of functionally related genes in Saccharomyces cerevisiae. Journal of molecular biology 2000, 296(5):1205–1214. 10.1006/jmbi.2000.3519

    Article  CAS  PubMed  Google Scholar 

  4. Gordan R, Narlikar L, Hartemink A: Finding regulatory DNA motifs using alignment-free evolutionary conservation information. Nucleic Acids Research 2010, 38(6):e90. 10.1093/nar/gkp1166

    Article  PubMed  PubMed Central  Google Scholar 

  5. Haudry Y, Ramialison M, Paten B, Wittbrodt J, Ettwiller L: Using Trawler standalone to discover overrepresented motifs in DNA and RNA sequences derived from various experiments including chromatin immunoprecipitation. Nature Protocols 2010, 5(2):323–334. 10.1038/nprot.2009.158

    Article  CAS  PubMed  Google Scholar 

  6. Das M, Dai H: A survey of DNA motif finding algorithms. BMC Bioinformatics 2007, 8(Suppl 7):S21. 10.1186/1471-2105-8-S7-S21

    Article  PubMed  PubMed Central  Google Scholar 

  7. Hertz G, Hartzell G III, Stormo G: Identification of consensus patterns in unaligned DNA sequences known to be functionally related. Bioinformatics 1990, 6(2):81. 10.1093/bioinformatics/6.2.81

    Article  CAS  Google Scholar 

  8. Sandelin A, Wasserman W, Lenhard B: ConSite: web-based prediction of regulatory elements using cross-species comparison. Nucleic acids research 2004, (32 Web Server):W249. 10.1093/nar/gkh372

    Google Scholar 

  9. Benos P, Lapedes A, Stormo G: Probabilistic code for DNA recognition by proteins of the EGR family. Journal of molecular biology 2002, 323(4):701–727. 10.1016/S0022-2836(02)00917-8

    Article  CAS  PubMed  Google Scholar 

  10. Bulyk M, Johnson P, Church G: Nucleotides of transcription factor binding sites exert interdependent effects on the binding affinities of transcription factors. Nucleic acids research 2002, 30(5):1255. 10.1093/nar/30.5.1255

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Eisen M: All motifs are not created equal: structural properties of transcription factor-DNA interactions and the inference of sequences specificity. Genome Biology 2005, 6(5):P7. 10.1186/gb-2005-6-5-p7

    Article  Google Scholar 

  12. Won K, Sandelin A, Marstrand T, Krogh A: Modeling promoter grammars with evolving hidden Markov models. Bioinformatics 2008, 24(15):1669. 10.1093/bioinformatics/btn254

    Article  CAS  PubMed  Google Scholar 

  13. Kim N, Tharakaraman K, Mariño-Ramírez L, Spouge J: Finding sequence motifs with Bayesian models incorporating positional information: an application to transcription factor binding sites. BMC bioinformatics 2008, 9: 262. 10.1186/1471-2105-9-262

    Article  PubMed  PubMed Central  Google Scholar 

  14. Tomovic A, Oakeley E: Position dependencies in transcription factor binding sites. Bioinformatics 2007, 23(8):933. 10.1093/bioinformatics/btm055

    Article  CAS  PubMed  Google Scholar 

  15. Zare-Mirakabad F, Ahrabian H, Sadeghi M, Nowzari-Dalini A, Goliaei B: New scoring schema for finding motifs in DNA Sequences. BMC Bioinformatics 2009, 10: 93. 10.1186/1471-2105-10-93

    Article  PubMed  PubMed Central  Google Scholar 

  16. Zadeh L: Fuzzy sets*. Information and control 1965, 8(3):338–353. 10.1016/S0019-9958(65)90241-X

    Article  Google Scholar 

  17. Atanassov K: Intuitionistic fuzzy sets. Fuzzy Sets and Systems 1986, 20: 87–96. 10.1016/S0165-0114(86)80034-3

    Article  Google Scholar 

  18. Garcia F, Lopez F, Cano C, Blanco A: FISim: A new similarity measure between transcription factor binding sites based on the fuzzy integral. BMC Bioinformatics 2009, 10: 224. 10.1186/1471-2105-10-224

    Article  PubMed  PubMed Central  Google Scholar 

  19. Lopez F, Blanco A, Garcia F, Cano C, Marin A: Fuzzy association rules for biological data analysis: a case study on yeast. BMC Bioinformatics 2008, 9: 107. 10.1186/1471-2105-9-107

    Article  PubMed  PubMed Central  Google Scholar 

  20. Liang L, Mandal V, Lu Y, Kumar D: MCM-test: a fuzzy-set-theory-based approach to differential analysis of gene pathways. BMC Bioinformatics 2008, 9(6):S16. 10.1186/1471-2105-9-S6-S16

    Article  PubMed  PubMed Central  Google Scholar 

  21. Chen J, Bardes E, Aronow B, Jegga A: ToppGene Suite for gene list enrichment analysis and candidate gene prioritization. Nucleic Acids Research 2009, (37 Web Server):W305. 10.1093/nar/gkp427

  22. Lam F, Steger D, O'Shea E: Chromatin decouples promoter threshold from dynamic range. Nature 2008, 453(7192):246–250. 10.1038/nature06867

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Sandelin A, Alkema W, Engstrom P, Wasserman W, Lenhard B: JASPAR: an open-access database for eukaryotic transcription factor binding profiles. Nucleic Acids Research 2004, (32 Database):D91. 10.1093/nar/gkh012

    Google Scholar 

  24. Van Helden J: Regulatory sequence analysis tools. Nucleic acids research 2003, 31(13):3593. 10.1093/nar/gkg567

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Hanley J, McNeil B: The meaning and use of the area under a receiver operating characteristic (ROC) curve. Radiology 1982, 143: 29–36.

    Article  CAS  PubMed  Google Scholar 

  26. Fawcett T: An introduction to ROC analysis. Pattern recognition letters 2006, 27(8):861–874. 10.1016/j.patrec.2005.10.010

    Article  Google Scholar 

  27. Chen X, Xu H, Yuan P, Fang F, Huss M, Vega V, Wong E, Orlov Y, Zhang W, Jiang J, et al.: Integration of external signaling pathways with the core transcriptional network in embryonic stem cells. Cell 2008, 133(6):1106–1117. 10.1016/j.cell.2008.04.043

    Article  CAS  PubMed  Google Scholar 

  28. Sharov A, Ko M: Exhaustive Search for Over-represented DNA Sequence Motifs with CisFinder. DNA Research 2009, 16:(5):261–73. 10.1093/dnares/dsp014

    Article  PubMed  PubMed Central  Google Scholar 

  29. Matys V, Kel-Margoulis O, Fricke E, Liebich I, Land S, Barre-Dirrie A, Reuter I, Chekmenev D, Krull M, Hornischer K, et al.: TRANSFAC (R) and its module TRANSCompel (R): transcriptional gene regulation in eukaryotes. Nucleic acids research 2006, (34 Database):D108. 10.1093/nar/gkj143

  30. Denning D: Invasive aspergillosis. Clinical infectious diseases 1998, 781–803. 10.1086/513943

    Google Scholar 

  31. Offner F, Cordonnier C, Ljungman P, Prentice H, Engelhard D, Bacquer D, Pauw B: Impact of previous aspergillosis on the outcome of bone marrow transplantation. Clinical infectious diseases 1998, 26(5):1098–1103. 10.1086/520274

    Article  CAS  PubMed  Google Scholar 

  32. Chamilos G, Luna M, Lewis R, Bodey G, Chemaly R, Tarrand J, Safdar A, Raad I, Kontoyiannis D: Invasive fungal infections in patients with hematologic malignancies in a tertiary care cancer center: an autopsy study over a 15-year period (1989–2003). Haematologica 2006, 91(7):986.

    PubMed  Google Scholar 

  33. Diop G, Spadoni J, Do H, Hirtzig T, Coulonges C, Labib T, Issing W, Rappaport J, Therwath A, Lathrop M, et al.: Genomic approach of AIDS pathogenesis: exhaustive genotyping of the TNFR1 gene in a French AIDS cohort. Biomedicine & Pharmacotherapy 2005, 59(8):474–480.

    Article  CAS  Google Scholar 

  34. Hehlgans T, Pfeffer K: The intriguing biology of the tumour necrosis factor/tumour necrosis factor receptor superfamily: players, rules and the games. Immunology 2005, 115: 1. 10.1111/j.1365-2567.2005.02143.x

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Bochud P, Chien J, Marr K, Leisenring W, Upton A, Janer M, Rodrigues S, Li S, Hansen J, Zhao L, et al.: Toll-like receptor 4 polymorphisms and aspergillosis in stem-cell transplantation. New England Journal of Medicine 2008, 359(17):1766. 10.1056/NEJMoa0802629

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. Baker E, Chen L, Smith C, Callen D, Goodwin R, Sutherland G: Chromosomal location of the human tumor necrosis factor receptor genes. Cytogenetics and cell genetics 1991, 57:(2–3):117.

    Article  Google Scholar 

  37. Sainz J, Salas I, Lopez E, Olmedo C, Comino A, Garcia F, Blanco A, Oyonarte S, Bueno P, Jurado M: TNFR1 mRNA Expression Level and TNFR1 Gene Polymorphisms Are Predictive Markers for Susceptibility to Develop Invasive Pulmonary Aspergillosis. International Journal of Immunopathology and Pharmacology 2009, 22(3):557–565.

    Google Scholar 

  38. Wingender E: TheTRANSFAC project as an example of framework technology that supports the analysis of genomic regulation. Brief Bioinform 2008, 9: 326–332. 10.1093/bib/bbn016

    Article  CAS  PubMed  Google Scholar 

  39. Dror N, Alter-Koltunoff M, Azriel A, Amariglio N, Jacob-Hirsch J, Zeligson S, Morgenstern A, Tamura T, Hauser H, Rechavi G, et al.: Identification of IRF-8 and IRF-1 target genes in activated macrophages. Molecular immunology 2007, 44(4):338–346. 10.1016/j.molimm.2006.02.026

    Article  CAS  PubMed  Google Scholar 

  40. Zhao J, Kong H, Li H, Huang B, Yang M, Zhu C, Bogunovic M, Zheng F, Mayer L, Ozato K, et al.: IRF-8/interferon (IFN) consensus sequence-binding protein is involved in Toll-like receptor (TLR) signaling and contributes to the cross-talk between TLR and IFN-γ signaling pathways. Journal of Biological Chemistry 2006, 281(15):10073. 10.1074/jbc.M507788200

    Article  CAS  PubMed  Google Scholar 

  41. Tamura T, Ozato K: Review: ICSBP/IRF-8: its regulatory roles in the development of myeloid cells. Journal of Interferon & Cytokine Research 2002, 22: 145–152.

    Article  CAS  Google Scholar 

  42. Tamura T, Nagamura-Inoue T, Shmeltzer Z, Kuwata T, Ozato K: ICSBP directs bipotential myeloid progenitor cells to differentiate into mature macrophages. Immunity 2000, 13(2):155–165. 10.1016/S1074-7613(00)00016-9

    Article  CAS  PubMed  Google Scholar 

  43. Schling P, Rudolph C, Heimerl S, Fruth S, Schmitz G: Expression of tumor necrosis factor alpha and its receptors during cellular differentiation. Cytokine 2006, 33(5):239–245. 10.1016/j.cyto.2006.02.007

    Article  CAS  PubMed  Google Scholar 

  44. Wasserman W, Sandelin A: Applied bioinformatics for the identification of regulatory elements. Nature Reviews Genetics 2004, 5(4):276–287. 10.1038/nrg1315

    Article  CAS  PubMed  Google Scholar 

  45. Szmidt E, Kacprzyk J: Intuitionistic fuzzy sets in group decision making. Notes on IFS 1996, 2: 11–14.

    Google Scholar 

  46. Atanassov K, Gargov G: Intuitionistic fuzzy logic. Compt Rend Acad Bulg Sci 1990, 43: 9–12.

    Google Scholar 

  47. De S, Biswas R, Roy A: An application of intuitionistic fuzzy sets in medical diagnosis. Fuzzy Sets and Systems 2001, 117(2):209–213. 10.1016/S0165-0114(98)00235-8

    Article  Google Scholar 

  48. Khatibi V, Montazer G: Intuitionistic fuzzy set vs. fuzzy set application in medical pattern recognition. Artificial Intelligence In Medicine 2009, 47: 43–52. 10.1016/j.artmed.2009.03.002

    Article  PubMed  Google Scholar 

  49. Hung W, Yang M: Similarity measures of intuitionistic fuzzy sets based on Hausdorff distance. Pattern Recognition Letters 2004, 25(14):1603–1611. 10.1016/j.patrec.2004.06.006

    Article  Google Scholar 

  50. Atanassov K: New operations defined over the intuitionistic fuzzy sets. Fuzzy Sets and Systems 1994, 61(2):142. 10.1016/0165-0114(94)90229-1

    Article  Google Scholar 

  51. Atanassov K: Intuitionistic Fuzzy Sets: Theory and Applications. Physica-Verlag, Heidelberg, New York 1999.

    Google Scholar 

Download references

Acknowledgements

We would like to thank the anonymous reviewers for their helpful comments which significantly improved the quality of the paper. This work has been carried out as part of projects P08-TIC-4299 of J. A., Sevilla and TIN2009-13489 of DGICT, Madrid.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Fernando Garcia-Alcalde.

Additional information

Authors' contributions

FG designed the study, designed and implemented the intuitionistic scoring, performed the experiments, helped with the analysis of the results and drafted the paper. AB assisted with the design of the study and helped to draft the paper. AS provided help with the analysis of the results and assisted in drafting the paper. All authors read and approved the final manuscript.

Electronic supplementary material

12859_2010_4134_MOESM1_ESM.XLS

Additional file 1: Synthetic sequences experiment. This files contains thresholded results for the different methods for the synthetic sequences experiment. (XLS 34 KB)

12859_2010_4134_MOESM2_ESM.PDF

Additional file 2: ROC curves. This file contains the ROC curves associated to the synthetic and mutated sequences experiments. (PDF 60 KB)

12859_2010_4134_MOESM3_ESM.ZIP

Additional file 3: FASTA sequences. This file contains the 50 sequence segments in FASTA format for each one of the motifs SMAD1, Myc, and STAT3. (ZIP 11 KB)

12859_2010_4134_MOESM4_ESM.XLS

Additional file 4: Motif statistics. This file contains some statistics for the experiment discussed in the section Real data. (XLS 9 KB)

Authors’ original submitted files for images

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Garcia-Alcalde, F., Blanco, A. & Shepherd, A.J. An intuitionistic approach to scoring DNA sequences against transcription factor binding site motifs. BMC Bioinformatics 11, 551 (2010). https://doi.org/10.1186/1471-2105-11-551

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1471-2105-11-551

Keywords