Skip to main content

G4PromFinder: an algorithm for predicting transcription promoters in GC-rich bacterial genomes based on AT-rich elements and G-quadruplex motifs

Abstract

Background

Over the last few decades, computational genomics has tremendously contributed to decipher biology from genome sequences and related data. Considerable effort has been devoted to the prediction of transcription promoter and terminator sites that represent the essential “punctuation marks” for DNA transcription. Computational prediction of promoters in prokaryotes is a problem whose solution is far from being determined in computational genomics. The majority of published bacterial promoter prediction tools are based on a consensus-sequences search and they were designed specifically for vegetative σ70 promoters and, therefore, not suitable for promoter prediction in bacteria encoding a lot of σ factors, like actinomycetes.

Results

In this study we investigated the possibility to identify putative promoters in prokaryotes based on evolutionarily conserved motifs, and focused our attention on GC-rich bacteria in which promoter prediction with conventional, consensus-based algorithms is often not-exhaustive. Here, we introduce G4PromFinder, a novel algorithm that predicts putative promoters based on AT-rich elements and G-quadruplex DNA motifs. We tested its performances by using available genomic and transcriptomic data of the model microorganisms Streptomyces coelicolor A3(2) and Pseudomonas aeruginosa PA14. We compared our results with those obtained by three currently available promoter predicting algorithms: the σ70consensus-based PePPER, the σ factors consensus-based bTSSfinder, and PromPredict which is based on double-helix DNA stability. Our results demonstrated that G4PromFinder is more suitable than the three reference tools for both the genomes. In fact our algorithm achieved the higher accuracy (F1-scores 0.61 and 0.53 in the two genomes) as compared to the next best tool that is PromPredict (F1-scores 0.46 and 0.48). Consensus-based algorithms produced lower performances with the analyzed GC-rich genomes.

Conclusions

Our analysis shows that G4PromFinder is a powerful tool for promoter search in GC-rich bacteria, especially for bacteria coding for a lot of σ factors, such as the model microorganism S. coelicolor A3(2). Moreover consensus-based tools and, in general, tools that are based on specific features of bacterial σ factors seem to be less performing for promoter prediction in these types of bacterial genomes.

Background

In all living organisms the flow of genetic information starts with gene transcription, an essential process that is tightly regulated at each step (initiation, elongation, termination). In bacteria and archaea a single RNA polymerase (RNAP) carries out this process, whereas in eukaryotes multiple different RNAP are responsible for transcription of different classes of genes [1]. Despite a lot of differences in transcription machinery among the three domains of life (including RNAP subunit composition with five subunits in bacterial (α2ββ′ω) and more than 12 subunits in archaeal and eukaryotic RNAP), evolutionary conserved features such as similar overall shape in RNAP, highly conserved active centers and similar contact to the nucleic acid chains have been recognized [2, 3]. Structural and functional similarities also extend to several accessory factors modulating the different steps of the transcription cycle.

In bacteria, specific transcription initiation requires sigma (σ) factors that, when bound to RNAP, recognize and melt promoters [4]. Based on sequence, structural and functional similarity, the bacterial σ-factors can be grouped into two families, the σ70- and σ54-family (this latter existing in most but not all bacteria), with little if any sequence identity between them [5, 6]. In contrast to bacterial RNAP, archaeal RNAP and eukaryotic RNAP II utilizes two key basal factors, transcription factor B (TFB for archaeal RNAP, TFIIB for eukaryotic RNAP II) and TATA-binding protein (TBP) rather than σ factors for transcription initiation. TFB/TFIIB and TBP bind to DNA and subsequently recruit RNAP and additional factors to form a core initiation complex [7, 8]. However, recently, structural comparison of initiating RNAP complexes and structure-based amino acid sequence alignment have provided evidence of structural and functional analogies, and evolutionary relatedness between bacterial σ70-family factors and archaeal/eukaryotic TFB/TFIIB suggesting a simple model for promoter evolution and genesis of transcription systems [9, 10]. The model is based on apparent conservation of helix-turn-helix (HTH) motifs in archaeal/eukaryotic TFB/TFIIB and bacterial σ70-family factors. These HTH motifs are involved in recognition of the structural promoter elements: a GC-rich “anchor sequence” (corresponding to bacterial − 35 element and archaeal/eukaryotic BREup) and a downstream located “AT-rich element” (corresponding to bacterial − 10 element [TATAAT, Pribnow box] and TATAAAAG boxes). Contact to double strand anchor DNA maintains the position of the most C-terminal HTH domain, while more N-terminal HTH domains facilitate bubble opening and initiation [9, 10].

Recently, G-quadruplex motifs, tertiary structures formed by nucleic acid sequences that are rich in guanine via non-Watson-Crick base pairing, have received a great deal of attention because of their putative role in promoter function [11]. In these dynamic structures, four guanine residues can associate through Hoogsteen hydrogen binding to form a square planar structure called a guanine tetrad, and two or more guanine tetrads can stack on top of each other to form a G-quadrulpex [12,13,14,15]. Interestingly, more than 40% of human gene promoters contain one or more G-quadruplex motifs [16]. In fungi G-quadruplex DNA motifs are significantly associated with promoter regions and to a lesser extent with open reading frames (ORFs) [17], and these DNA motifs are more conserved than expected from a random distribution among related fungi suggesting in vivo functions that are under evolutionary constraint [18]. Conserved G-quadruplex DNA motifs have been also reported in promoters of orthologous gene across phylogenetically distant prokaryotes [19], and, very recently, a conserved putative G-quadruplex-Hairpin-Duplex switch has been described [20].

The evolutionary relatedness and/or functional analogies in transcription initiation mechanisms between all three domains of life prompted us to explore the possibility of recognizing promoter elements in prokaryotic genomes based on conserved structured motifs. In particular, we focused our attention on GC-rich bacterial genomes where promoter prediction with conventional, consensus-based algorithms is often difficult and certainly not exhaustive [21]. Promoter prediction is especially problematic in actinomycetes, a group of mycelial organisms with complex transcriptional patterns because their large genomes may encode more than 60 sigma factors [22], although consensus sequences have been proposed for computer assisted promoter identification and classification in Streptomyces spp. [23]. In this study we have developed an algorithm to identify putative promoters based on AT-rich elements and G-quadruplex DNA motifs in the GC-rich “anchor sequence”, and tested its performances by using available genomic and transcriptomic data of the model microorganisms Streptomyces coelicolor A3(2) [22] and Pseudomonas aeruginosa PA14 [24]. Results were compared with those obtained by some currently available tools for bacterial promoter prediction. Currently available tools for prokaryote promoter prediction include BPROM [25], NNPP2 [26], PePPER [27], PromPredict [28] and bTSSfinder [29]. BPROM, NNPP2 and PePPER are tools for prediction of prokaryote promoter elements based on a consensus-sequences search and they were designed specifically for vegetative σ70 promoters. bTSSfinder is the most recent consensus-based promoter prediction algorithm in prokaryotes. It extends the concept of consensus prediction to five classes of σ factors in E. coli70, σ38, σ32, σ28 and σ24) and to five classes of σ factors in CyanobacteriaA, σC, σH, σG and σF). This tool performed successfully in E. coli genome achieving very high accuracy values (F1-score = 0.93) [29]. PromPredict instead identifies promoter regions on the basis of DNA double helix stability, therefore using a different strategy than consensus-based algorithms. In fact, PromPredict algorithm is based on the general observation that promoter regions are less stable than flanking regions [21, 30]. For this reason, PromPredict is a more general tool than consensus-based tools and could be more suitable in GC-rich bacteria featuring diverse σ factors. For comparison, we focused our attention on PromPredict [31], on the most recent consensus-based tools PePPER [32] and on bTSSfinder [33]. We excluded from the comparison BPROM and NNPP2 because they work similarly to the most recent PePPER. All these tools are designed for a genome-wide prediction. PePPER was optimized for E. coli, PromPredict for both E. coli and B. subtilis, while bTSSfinder for E. coli and Cyanobacteria.

Implementation

Programming language and data sets

G4PromFinder algorithm was implemented in Python (v.3.5) [34], and works as a genome-wide promoter predictor taking as input bacterial genome-sequences. In particular, we used available genomic sequences (from National Center for Biotechnology Information) of the model microorganisms S. coelicolor A3(2) (accession code NC_003888.3) and P. aeruginosa PA14 (accession code NC_008463.1) (see below) for promoter predictions and their genomic annotation together with transcriptomic data [22, 24] for the prediction quality evaluation. The method to identify putative promoter elements is described below.

Method to identify putative promoters

A two-step procedure was used to detect putative promoters (Fig. 1). The first step consisted in the identification of the putative promoter “AT-rich element”. To this purpose, the algorithm slides a window of 25 bp over the query sequence, 1 bp at a time, until the AT% content of the window reaches the threshold value of 40%. Afterwards, by scanning a window of 75 bp (starting from the position where the threshold value of the AT% content was reached), the 25 bp long region with maximal AT content (herein referred to as AT-rich element) is selected. The second step was the identification of putative G-quadruplex motifs extended up to 50 bp upstream from the 5′-end of the selected AT-rich element. Motif GxNyGxNyGxNyGx with 2 ≤ x ≤ 4, 1 ≤ y ≤ 10 and maximum length of 30 bp is commonly used to predict the presence of G-quadruplexes [35]. G-quadruplexes could have an influence on gene expression also when localized on the reverse strand relative to transcription direction [36]. For this reason we searched for putative G-quadruplex motifs on either sense or antisense strand (motif CxNyCxNyCxNyCx with 2 ≤ x ≤ 4, 1 ≤ y ≤ 10 and maximum length of 30 bp was used to predict putative G-quadruplex on antisense strand). We considered as a single prediction all the predictions that were within 35 bp from each other, because most signals involved in determining TSS are located in the short region between the − 35 and the − 10 boxes. It is also relevant to point out that the length of the regions evaluated in the first step (25 bp) and in the second step (50 bp) are arbitrary; they were determined experimentally to optimize our search, and were compatible with the overall geometry of bacterial RNAP-promoter complex and with the proposed model for the genesis of transcription systems [9, 10]. Finally, two additional features of our algorithm are: i) the possibility of predicting multiple putative promoters in a single query region and ii) the possibility of searching for promoters in both strands.

Fig. 1
figure 1

Method used for the prediction and the validation of putative promoters

TSS global map datasets

To evaluate the reliability of our promoter predictions we used TSS global maps obtained by dRNAseq experiments. For S. coelicolor A3(2) 3570 TSSs were identified [22], and were categorized by their positions relative to known coding sequences (CDSs) giving 2771 primary TSSs (P) associated with currently annotated genes, which corresponds to 35.0% of the total genes in the S. coelicolor genome. In addition to P, 333 secondary TSSs (S) were identified revealing a total of 297 transcription units initiated by more than one TSS. 256 TSSs mapped in the antisense strand (A) of 241 genes, while 79 internal TSSs (I) were detected within 73 genes. Finally, 131 TSSs were mapped to IRs with no previously associated genes (N). For P. aeruginosa PA14, 2117 TSSs were predicted spanning 3325 protein coding genes (55% of all protein coding genes) [24]. In this last study, TSSs were not categorized.

Generation of the positive and negative sets of sequences

Using the publicly available genomes of S. coelicolor A3(2) and P. aeruginosa PA14 (accessions above) and the above-indicated TSS annotations, we created, for each of the two genomes, a promoter set (positive set) consisting of 251 bp long sequences covering the region from − 200 bp to + 51 bp with respect to each experimentally annotated TSS. Hence, for S. coelicolor A3(2) and P. aeruginosa PA14 genomes, the positive set consisted of 3570 and 2117 sequences, respectively.

As a negative set of sequences we considered all the IRs < = 251 bp and > = 50 bp in length in which promoters are not expected. Precisely, we considered all the previous IRs that separated two convergently oriented CDSs. In order to compare regions of same length, we decided to extend IRs < 251 bp equally from both their extremities until the length of 251 bp was reached. Finally, we assessed the total absence of annotated TSSs in the previous regions. For S. coelicolor A3(2) genome, the negative set consisted of 548 sequences, while for P. aeruginosa PA14 genome it consisted of 338 sequences.

Evaluation of the performances of the promoter predictor

To estimate the performances of G4PromFinder, we used the following statistical measures:

  • Recall (sensitivity or the true positive rate) = TP/(TP + FN)

  • Precision (the positive predictive value) = TP/(TP + FP)

  • Specificity (the true negative rate) = TN/(TN + FP)

  • Accuracy (the fraction of samples correctly classified) = (TP + TN)/(TP + TN + FP + FN)

  • F1-score (the harmonic mean of Precision and Accuracy) = 2*Precision*Recall/(Precision+Recall)

  • where TP = True positives, FP = False positives, FN = False negatives and TN = True negatives.

In accordance with the validation strategies adopted in previous studies [29], we considered a predicted promoter as a true positive (TP) if it started within 50 bp from an experimentally derived TSS (upstream or downstream). It is important to point out that at most one TP was considered for each sequence of the positive set. We considered as false positives (FP) all the samples of the negative set in which the algorithm predicted at least a promoter, as true negatives (TN) the sequences of the negative set in which the algorithm did not predict promoters and, finally, as false negatives (FN) the sequences of the positive set in which TPs were absent.

Results

Genome statistics, IRs and promoter prediction with consensus sequence-based algorithm

Data sources in this study were: i.) the annotated, whole-genome sequences of S. coelicolor A3(2) [37] and P. aeruginosa PA14 [38] for promoter prediction; ii) the TSS list obtained by dRNAseq experiments [22, 24] for promoter prediction validation. The complete genome of S. coelicolor A3(2) has a GC-content of 72.1% and consists of putative 7825 genes (mean GC-content 72.1%), with the median IR length of 118 bp, the first quartile IR length of 13 bp and the third of 162 bp (Fig. 2). Regarding P. aeruginosa PA14, its genome, with a GC-content of 66.3%, consists of 5973 annotated genes (mean GC-content 66.2%) with the median IR length of 118 bp, the first quartile IR length of 61 bp and the third of 211 bp (Fig. 2).

Fig. 2
figure 2

Boxplot of IRs length for S. coelicolor (a) and P. aeruginosa (b)

Preliminarily, in all IRs we searched for the presence of the “-35 consensus sequence” (TTGAC for S. coelicolor A3(2) and TTGNC for P. aeruginosa PA14) and the “-10 consensus sequence” (TANNNT for S. coelicolor A3(2) and TANAAT for P. aeruginosa PA14) for σ70-family factors [23, 39] separated by a sequence ranging in length from 16 to 22 bp. We detected both these sequences only in 1.6% and 0.4% of IRs of S. coelicolor A3(2) and P. aeruginosa PA14, respectively. This result further shows how difficult is to predict promoters by using consensus-based methods in GC-rich bacterial genomes.

Promoter prediction by AT-rich element and G-quadruplex motif-based algorithm and evaluation

Statistics of putative promoters that were predicted by G4PromFinder algorithm in the positive set are summarized in Table 1. Preliminarly we considered all the putative promoters predicted by G4PromFinder, without further constraints. G4PromFinder predicted, respectively, for S. coelicolor A3(2) and P. aeruginosa PA14, putative promoters in almost all the examined regions, precisely in 91.2% and 91.5% of such regions (Table 1). Overall, the algorithm predicted, respectively, 3751 and 2305 putative promoters in the two genomes. Therefore multiple putative promoters were associated with some of the samples. Precisely, in 13.8 and 17.4% of examined regions in S. coelicolor A3(2) and P. aeruginosa PA14, respectively, more than one predicted promoter could be found.

Table 1 Statistics of predicted promoters by G4PromFinder algorithm

We evaluated G4PromFinder performances on promoter prediction using a positive sequence set including all regions surrounding by a dRNAseq verified TSS and a negative sequence set composed by short IRs located between two convergently oriented CDSs, for both S. coelicolor A3(2) and P. aeruginosa PA14 genomes (see Implementation section for details). To fairly compare the positive and negative sets, which originally do not have the same size, we decided to randomly select 548 and 338 regions of the positive sets, respectively in S. coelicolor A3(2) and P. aeruginosa PA14 genomes, and we repeated the testing 10 times on different series of randomly selected sequences to obtain the mean values reported in Table 2 (column 1 and 2). We observed good performances in both bacterial genomes. In fact the F1-scores obtained in S. coelicolor A3(2) and P. aeruginosa PA14 were 0.61 and 0.53, respectively (Table 2). Recall-values obtained were high, about 70% (precisely 70.1 and 69.0%), while precision-values were lower (54.3 and 43.1%). Interestingly, in S. coelicolor A3(2) about 40% of validated promoters contained the “-10 consensus sequence” (TANNNT) that was previously proposed for σ70-family factors in Streptomycetes [23] (Table 3). In contrast, only a low percentage of validated promoters (6.1%) contained the “-35 consensus sequence” (TTGAC) for σ70-family factors [22] (Table 3). In P. aeruginosa PA14, the “-10 consensus sequence” (TANAAT) and the “-35 consensus sequence” (TTGNC) for σ70-family factors were contained in 7.4% and 28.2% of validated promoters, respectively (Table 3). Moreover, the mean AT content of the validated promoter AT-rich elements obtained was rather higher than the threshold value of 40% (see Implementation section), 48.5% and 53.3% in S. coelicolor A3(2) and P. aeruginosa PA14, respectively. These values are also higher than the mean AT content of the total validated promoters (Table 3).

Table 2 - Testing results of G4PromFindera
Table 3 – Some features of the validated promoters

A negative control: Specificity control of the promoter element G-quadruplex

In these GC-rich bacterial genomes, the G-quadruplex motif (see Implementation section) occurs very frequently. In fact we found about 120,000 and 70,000 instances of the G-quadruplex motif in S. coelicolor A3(2) and P. aeruginosa PA14 genomes. For this reason, we decided to carry out a negative control, in order to assess the not-random presence of a G-Quadruplex motif in a promoter. This control consisted in the searching for a random sequence motif with the same frequency as the G-quadruplex motif in the genome, in the identification of putative promoters as AT-rich elements that were preceded by this random sequence motif, and in their subsequent validation by using the same procedure adopted for G4PromFinder prediction. We used as random motif tetranucleotides sequences with a GC content similar to that of the entire genomes, preceded and followed by 13 random bp (in order to have a motif of length similar to G-quadruplex motif). In S. coelicolor A3(2) we carried out two controls, each with two pairs of tetranucleotide sequences (GCAG and GCTG; CACG and TCGC) that together have the same frequency to the G-quadruplex motif, while in P. aeruginosa PA14 we carried out one control with a pair of tetranucleotide sequences (GACG and ACGC). The random approach achieved lower accuracy (in S. coelicolor A3(2), F1-score for the first pair of tetranucleotides 0.56, F1-score for the second pair of tetranucleotides 0.51; in P. aeruginosa PA14 F1-score 0.45) compared to G4PromFinder (F1-score 0.61 and 0.53 in S. coelicolor A3(2) and P. aeruginosa PA14, Table 2). From these tests resulted that the fraction of the positive results obtained by G4PromFinder that surely were not a consequence of random chance is 0.12 in S. coelicolor A3(2) (1–0.535/0.61; 0.535 mean value of F1-scores of the two negative controls) and 0.15 in P. aeruginosa PA14 (1–0.45/0.53). The presence of AT-rich element in our negative control is the most probable reason for the relatively high performances in promoter predictions obtained by it. Indeed AT-rich element has by itself a well-known role in promoter regions definition. In any case the random approach achieved lower accuracy compared to G4PromFinder, and we can conclude that the G-quadruplex element has a higher specificity in the association with the AT-rich element compared to a random sequence with the same frequency of G-quadruplex and with a GC-content similar to that of the whole genome.

Comparison with PePPER, PromPredict and bTSSfinder tools.

We compared our results with those obtained by PePPER [27], PromPredict [28] and bTSSfinder [29] tools. PePPER predicts prokaryote promoters based on a consensus-sequences search. Precisely, PePPER software looks for the “-35 consensus sequence” and “-10 consensus sequence” for σ70-family factors of Escherichia coli allowing a certain degree of variability for the bases belonging to the consensus sequences. PePPER takes the annotated bacterial genome sequence as input and provides as output the positions of putative TSS, “-35 consensus sequence” and “-10 consensus sequence”, with a score assigned to them that indicates the probability that the extracted region actually corresponds to a promoter. In contrast, PromPredict predicts prokaryote promoters based on differences in DNA double helix stability in promoter and non-promoter regions, taking as input bacterial genome sequences and providing as output promoter coordinates, with a reliability level assigned to them. bTSSfinder, instead, predicts putative promoters for five classes of σ factors in CyanobacteriaA, σC, σH, σG and σF) and for five classes of σ factors in E. coli70, σ38, σ32, σ28 and σ24) taking as input bacterial genome sequences and providing as output TSS coordinates [29]. Preliminarily we run the three comparison tools on the whole genome of S. coelicolor A3(2) and P. aeruginosa PA14 considering all the identified promoter regions independently from their score. Table S1 (in Additional file 1) shows the global numbers of prediction obtained by each tool in comparison to G4PromFinder whole genome predictions. Then we intersected the three tool genome wide predictions with the positive and negative region sets already used for the evaluation of G4PromFinder performances (see Implementation). We considered as positive intersections only the predicted promoters falling for their entire length within those regions (i.e. the predicted promoters by PePPER and PromPredict whose coordinates falling within those regions and all the predicted TSSs by bTSSfinder falling within those regions). We defined true promoters those whose position difference between the experimentally derived TSS and the predicted promoter regions was included in the range ± 50 bp, as we have done for G4PromFinder.

The results for S. coelicolor A3(2) and P. aeruginosa PA14 are shown in Table 4. We evaluated the performances of the four methods by using precision, recall and F1-score. This comparison clearly indicates that G4PromFinder has significantly higher prediction accuracy, both in S. coelicolor A3(2) and in P. aeruginosa PA14. In fact, as presented in Table 4, G4PromFinder produced the best performance (F1-score 0.61 in S. coelicolor A3(2), F1-score 0.53 in P. aeruginosa PA14), followed by PromPredict (F1-score 0.46 in S. coelicolor A3(2), F1-score 0.48 in P.aeruginosa PA14), bTSSfinder for E. coli σ factors (F1-score 0.38 in S. coelicolor A3(2), F1-score 0.36 in P.aeruginosa PA14), and finally bTSSfinder for Cyanobacteria σ factors (F1-score 0.28 in S. coelicolor A3(2), F1-score 0.28 in P.aeruginosa PA14) and PePPER (F1-score 0.32 in S. coelicolor A3(2), F1-score 0.42 in P.aeruginosa PA14). Therefore consensus-based algorithms (PePPER and bTSSfinder) produced lower performances with the analyzed GC-rich genomes. Actually PePPER algorithm provided the highest precision values, but it also produced the lowest recall values and, for this reason its F1-scores were very low.

Table 4 Comparison between G4PromFinder, PePPER, PromPredict and bTSSfinder testing resultsa

Moreover, we carried out another comparison between G4PromFinder and the available promoter prediction programs. To perform this analysis, we considered all the predicted promoters by the four programs that were within the regions of the positive set, considering now as false positives the predictions whose distance from the annotated TSS was more than 50 bp. The results of the comparison are presented in Table 5, and again they clearly show that G4PromFinder has the highest prediction accuracy in these bacterial genomes.

Table 5 Comparison between G4PromFinder and available promoter prediction programs assessed on all the samples of the positive sets

In Fig. 3 we, instead, show the distributions of distances occurring between the 3′-end points of validated promoters of G4PromFinder and the TSSs used for validation. In all examined cases, we noticed a peak of distribution around the “-10” value.

Fig. 3
figure 3

Distribution of validated promoters in S. coelicolor A3(2) (a) and P. aeruginosa PA14 (b) as a function of their distance from the TSSs obtained by dRNAseq experiments and used for validation.Predicted promoters are grouped based on distances between the AT-rich element 3′-end points and the annotated TSS. A: predicted promoters in S. coelicolor A3(2); B: predicted promoters in P. aeruginosa PA14

Discussion and conclusions

In this study we investigated the possibility of predicting prokaryotic promoters by detecting evolutionarily conserved motifs. We focused on possible G-quadruplex structures upstream of AT-rich elements. The rationale started from the evidence that in human, yeast and bacterial genomes G-quadruplexes are overrepresented in promoter-proximal regions [18, 19, 40, 41]. In this study we showed that an AT-rich element preceded by a G-quadruplex motif is within ±50 bp from an experimentally identified TSS in 75.6 and 73.4% of total cases, in S. coelicolor A3(2) and P. aeruginosa PA14 genomes, respectively (Table 5). These high percentages support the idea that G-quadruplex is a prototypical motif involved in general promoter function/regulation.

G-quadruplex are highly dynamic structures whose thermal stability is affected by a number of features including the number of G-quartets present in the structure, the length and the composition of the loops formed by non-guanine bases [42]. Many G-quadruplex DNA structures, once folded, are more thermodynamically stable than double-strand DNA in vitro, and their unfolding kinetics are much slower than those of DNA or RNA hairpin structures [43]. As G-quadruplexes are likely to inhibit DNA and RNA metabolism, their formation must be regulated, and recently, a number of proteins that specifically regulate G-quadruplex folding and unfolding have been identified [41].

There is evidence that G-quadruplex formation in promoter “anchor” (− 35 sequence) elements could impair transcription initiation by RNA polymerase, or if present in the antisense strand of bacterial σ70 promoter between “anchor” and “AT-rich” (− 10 sequence) element could impair the initiation-elongation transition (the so-called promoter clearance) [11, 36]. On one hand, recognition of double strand “anchor” sequence in promoters may be strongly influenced by G-quadruplex that could create a physical barrier that hinders RNAP binding or complicates promoter recognition by σ factors. On the other hand, RNAP binding might also facilitate G-quadruplex formation on antisense strand after promoter melting, which ultimately might hamper the initiation-elongation transition [36]. Regulation of G-quadruplex folding and unfolding by G-quadruplex-binding proteins might represent a general mechanism to modulate promoter activity.

Noticeably, less than half of validated promoters that were identified by the algorithm in S. coelicolor A3(2) genome contained the “-10 consensus sequence” (TANNNT) (Table 3) that was previously proposed for σ70-family factors in Streptomycetes [23]. In contrast, a very small percentage of putative promoters in S. coelicolor contained the proposed “-35 consensus sequence” (TTGAC) in these bacteria [23] (Table 3). This finding, which may be explained by the occurrence of a huge number of σ-factors in Streptomycetes, confirms how difficult is to identify promoters in Streptomycetes with conventional, consensus-based algorithms.

The evaluation of G4PromFinder performances on S. coelicolor A3(2) and P. aeruginosa PA14 show high recall values 70.1% and 69.0% respectively (Table 2), but also low specificity values. This is particularly striking for P. aeruginosa PA14 genome whose specificity results only 8.9%, (Table 2) if compared to that obtained in S. coelicolor A3(2) genome 40.8% (Table 2). We would have expected a lower specificity in S. coelicolor A3(2), where the G-quadruplex motif has a higher density (see Results section). This result instead could suggest that G4PromFinder specificity could be linked to the genome GC richness, because the GC-content of S. coelicolor A3(2) (72.1%) is higher than that of P. aeruginosa PA14 (66.3%). Also at the genome-wide level (Additional file 1: Table S1) we can see number of predictions higher than the number of annotated TSSs but, compared to the other tools, G4PromFinder shows numbers of predictions of the same order of magnitude. The only exception is PePPER which seems the most restricitive. The high number of predictions is probably due to multiple causes, such as, for example, the lack of complete and precise TSS maps and the existence of unknown repression mechanisms for which some computationally predicted promoters are not used in vivo.

The comparison of G4PromFinder predictions with those obtained by PePPER [27], PromPredict [28] and bTSSfinder [29] tools, highlighted also its reliability. Indeed our analysis showed that G4PromFinder produces the best performances in both the genomes, obtaining as F1-score 0.61 and 0.53 in S. coelicolor A3(2) and P. aeruginosa PA14, compared to the next best tool PromPredict (F1-score 0.46 and 0.48 in S. coelicolor A3(2) and P. aeruginosa PA14). The σ factors consensus-based tool bTSSfinder and especially the consensus-based PePPER, which was designed specifically for vegetative σ70 promoters, achieved the lowest accuracy (Table 4). This is a further confirmation that a promoter prediction with conventional, consensus-based algorithms is often difficult in this type of bacteria, especially in bacteria coding for several σ factors, like actinomycetes. In fact PePPER, despite achieved the highest precision values, identified very few promoters in the examined genomes; as consequence its recall values were very low (0.20 and 0.31, Table 4). Even bTSSfinder, that in E. coli achieved the highest accuracy (F1-score 0.93 [29]) compared to the available tools, failed in these genomes. The same results were obtained also when we tested G4PromFinder and the three comparison tools only on the sequences of the positive set (Table 5). G4PromFinder produced the best results and PromPredict was the best of the three available tools used for comparison. Compared to the results reported in Table 4, the only difference was a general increase in F1-scores.

On the basis of these findings, we believed that G4PromFinder is a very powerful tool in GC rich bacterial genomes when compared to currently available tools, which are instead suitable in predicting promoter regions in other genomes, especially E. coli for which they were optimized.

Abbreviations

A:

Antisense TSS

CDS:

Coding sequence

dRNAseq:

Differential RNAseq

FN:

False negatives

FP:

False positives

HTH:

Helix-turn-helix

I:

Internal TSS

IR:

Intergenic region

N:

TSS with no previously associated genes

ORF:

Open reading frame

P:

Primary TSS

RNAP:

RNA polymerase

S:

Secondary TSS

TBP:

TATA-binding protein

TFB:

Transcription factor B for archaeal RNAP

TFIIB:

Transcription factor B for eukaryotic RNAP II)

TP:

True positives

TSS:

Transcription start site

σ factor:

Sigma factor

References

  1. Burton ZF. The old and new testaments of gene regulation: evolution of multi-subunit RNA polymerases and co-evolution of eukaryote complexity with the RNAP II CTD. Transcription. 2014;5:1–12.

    Google Scholar 

  2. Iyer LM, Koonin EV, Aravind L. Evolutionary connection between the catalytic subunits of DNA-dependent RNA polymerases and eukaryotic RNA-dependent RNA polymerases and the origin of RNA polymerases. BMC Struct Biol. 2003;3:1.

    Article  PubMed  PubMed Central  Google Scholar 

  3. Iyer LM, Aravind L. Insights from the architecture of the bacterial transcription apparatus. J Struct Biol. 2012;179:299–319.

    Article  CAS  PubMed  Google Scholar 

  4. Feklístov A, Sharon BD, Darst SA, Gross CA. Bacterial sigma factors: a historical, structural, and genomic perspective. Annu Rev Microbiol. 2014;68:357–76.

    Article  PubMed  Google Scholar 

  5. Wösten MM. Eubacterial sigma-factors. FEMS Microbiol Rev. 1998;22:127–50.

    Article  PubMed  Google Scholar 

  6. Paget MS, Helmann JD. The sigma70 family of sigma factors. Genome Biol. 2003;4:203.

    Article  PubMed  PubMed Central  Google Scholar 

  7. Werner F, Grohmann D. Evolution of multisubunit RNA polymerases in the three domains of life. Nat Rev Microbiol. 2011;9:85–98.

    Article  CAS  PubMed  Google Scholar 

  8. Mühlbacher W, Sainsbury S, Hemann M, Hantsche M, Neyer S, Herzog F, Cramer P. Conserved architecture of the core RNA polymerase II initiation complex. Nat Commun. 2014;5:4310.

    Article  PubMed  Google Scholar 

  9. Burton SP, Burton ZF. The σ enigma: bacterial σ factors, archaeal TFB and eukaryotic TFIIB are homologs. Transcription. 2014;5:e967599.

    Article  PubMed  PubMed Central  Google Scholar 

  10. Burton ZF, Opron K, Wei G, Geiger JHA. Model for genesis of transcription systems. Transcription. 2016;7:1–13.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Rhodes D, Lipps HJ. G-Quadruplexes and their regulatory roles in biology. Nucleic Acids Res. 2015;43:8627–37.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Gellert M, Lipsett MN, Davies DR. Helix formation by guanylicacid. Proc Natl Acad Sci U S A. 1962;48:2013–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. Simonsson T. G-quadruplex DNA structures–variations on a theme. Biol Chem. 2001;382:621–8.

    Article  CAS  PubMed  Google Scholar 

  14. Burge S, Parkinson GN, Hazel P, Todd AK, Neidle S. Quadruplex DNA: sequence, topology and structure. Nucleic Acids Res. 2006;34:5402–15.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Stegle O, Payet L, Mergny JL, MacKay DJ, Leon JH. Predicting and understanding the stability of G-quadruplexes. Bioinformatics. 2009;25:374–82.

    Article  Google Scholar 

  16. Huppert JL, Balasubramanian S. G-quadruplexes in promoters throughout the human genome. Nucleic Acids Res. 2007;35:406–13.

    Article  CAS  PubMed  Google Scholar 

  17. Hershman SG, Chen Q, Lee JY, Kozak ML, Yue P, et al. Genomic distribution and functional analyses of potential G-quadruplex-forming sequences in Saccharomyces cerevisiae. Nucleic Acids Res. 2008;36:144–56.

    Article  CAS  PubMed  Google Scholar 

  18. Capra JA, Paeschke K, Singh M, Zakian VA. G-quadruplex DNA sequences are evolutionarily conserved and associated with distinct genomic features in Saccharomyces cerevisiae. PLoS Comput Biol. 2010;6:e1000861.

    Article  PubMed  PubMed Central  Google Scholar 

  19. Rawal P, Kummarasetti VB, Ravindran J, Kumar N, Halder K, et al. Genome-wide prediction of G4 DNA as regulatory motifs: role in Escherichia coli global regulation. Genome Res. 2006;16:644–55.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Kaplan OI, Berber B, Hekim N, Doluca O. G- quadruplex prediction in E. coli genome reveals a conserved putative G-quadruplex-hairpin-duplex switch. Nucleic Acids Res. 2016;44:9083–95.

    CAS  PubMed  PubMed Central  Google Scholar 

  21. Rangannan V, Bansal M. High-quality annotation of promoter regions for 913 bacterial genomes. Bioinformatics. 2010;26:3043–50.

    Article  CAS  PubMed  Google Scholar 

  22. JeongY, Kim JN, Kim MW, Bucca G, Cho S, Yoon YJ, Kim BG, Roe JH, Kim SC, Smith CP and Cho BK. The dynamic transcriptional and translational landscape of the model antibiotic producer Streptomyces coelicolor A3(2). Nature Commun. 2016;7:11605.

  23. Bourn WR, Babb B. Computer assisted identification and classification of streptomycete promoters. Nucleic Acids Res. 1995;23:3696–703.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. Wurtzel O, Yoder-Himes DR, Han K, Dandekar AA, Edelheit S, Greenberg EP, Sorek R, Lor S. The single-nucleotide resolution Transcriptome of Pseudomonas Aeruginosa grown in body temperature. PLoS Pathog. 2012;8(9):e1002945. https://0-doi-org.brum.beds.ac.uk/10.1371/journal.ppat.1002945.

    Article  PubMed  PubMed Central  Google Scholar 

  25. Solovyev V, Salamov A. Automatic annotation of microbial genomes and metagenomic sequences. In: Li RW, editor. Metagenomics and its applications in agriculture, biomedicine and environmental studies. Hauppauge: Nova Science Publishers; 2011. p. 61–78.

    Google Scholar 

  26. Reese MG. Application of a time-delay neural network to promoter annotation in the Drosophila Melanogaster genome. Comput Chem. 2001;26:51–6.

    Article  CAS  PubMed  Google Scholar 

  27. De Jong A, Pietersma H, Cordes M, Kuipers OP, Kok J. PePPER: a webserver for prediction of prokaryote promoter elements and regulons. BMC Genomics. 2012;13:299.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Rangannan V, Relative BM. Stability of DNA as a generic criterion for promoter prediction: whole genome annotation of microbial genomes with varying nucleotide base composition. Mol BioSyst. 2009;5:1758–69.

    Article  CAS  PubMed  Google Scholar 

  29. Shahmuradov IA, Razali RM, Bougouffa S, Radovanovic A, Bajic VB. bTSSfinder: a novel tool for the prediction of promoters in cyanobacteria and Escherichia Coli. Bioinformatics. 2017;33(3):334–40.

    PubMed  Google Scholar 

  30. Rangannan V, Identification BM. Annotation of promoter regions in microbial genome sequences on the basis of DNA stability. J. Biosci. 2007;32:851–62.

    Article  CAS  Google Scholar 

  31. PromPredict (Web server for promoter identification in genomic DNA sequence). http://nucleix.mbu.iisc.ac.in/prompredict/prompredict.html. Accessed 25 May 2017.

  32. PePPER, prediction of prokaryote promoters. http://genome2d.molgenrug.nl/index.php/prokaryote-promoters. Accessed 20 Apr 2017.

  33. bTSSfinder. http://www.cbrc.kaust.edu.sa/btssfinder/.

  34. Python. https://www.python.org. Accessed 20 Nov 2016.

  35. Kikin O, D’Antonio L, Bagga PS. QGRS Mapper: a web-based server for predicting G-quadruplexes in nucleotide sequences. Nucleic Acids Res. 2006;34:W676–82.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. Holder IT, Hartig JS. A matter of location: influence of G-quadruplexes on Escherichia coli gene expression. Chem Biol. 2014;21:1511–21.

    Article  CAS  PubMed  Google Scholar 

  37. Bentley SD, Chater KF, Cerdeno-Tarraga AM, Challis GL, Thomson NR, James KD, Harris DE, Quail MA, Kieser H, Harper D, Bateman A, Brown S, Chandra G, Chen CW, Collins M, Cronin A, Fraser A, Goble A, Hidalgo J, Hornsby T, Howarth S, Huang CH, Kieser T, Larke L, Murphy L, Oliver K, O'Neil S, Rabbinowitsch E, Rajandream MA, Rutherford K, Rutter S, Seeger K, Saunders D, Sharp S, Squares R, Squares S, Taylor K, Warren T, Wietzorrek A, Woodward J, Barrell BG, Parkhill J, Hopwood DA. Complete genome sequence of the model actinomycete Streptomyces Coelicolor A3(2). Nature. 2002;417:141–7.

    Article  PubMed  Google Scholar 

  38. Lee DG, Urbach JM, Wu G, Liberati NT, Feinbaum RL, Miyata S, Diggins LT, He J, Saucier M, Déziel E, Friedman L, Li L, Grills G, Montgomery K, Kucherlapati R, Rahme LG, Ausubel FM. Genomic analysis reveals that Pseudomonas Aeruginosa virulence is combinatorial. Genome Biol. 2006;7:R90.

    Article  PubMed  PubMed Central  Google Scholar 

  39. Potvin E, Sanschagrin S, Levesque RC. Sigma factors in Pseudomonas aeruginosa. FEMS Microbial Rev. 2008:38–55.

  40. Verma A, Halder K, Halder R, Yadav VK, Rawal P, Thakur RK, Mohd F, Sharma A, Genome-wide CS. Computational and expression analyses reveal G-quadruplex DNA motifs as conserved cis-regulatory elements in human and related species. J Med Chem. 2008;51:5641–9.

    Article  CAS  PubMed  Google Scholar 

  41. Huppert JL, Balasubramanian S. Prevalence of quadruplexes in the human genome. Nucleic Acids Res. 2005;33:2908–16.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. König SLB, Evans AC, Huppert JL. Seven essential questions on G-quadruplex. Biomol Concepts. 2010;1:197–213.

    Article  PubMed  Google Scholar 

  43. Lane AN, Chaires JB, Gray RD, Trent JO. Stability and kinetics of G-quadruplex structures. Nucleic Acids Res. 2008;36:5482–515.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

Not applicable.

Funding

This work was partially supported by the Italian Ministry for Education,

Universities and Research (Grant number PON01_02093).

Availability of data and materials

The datasets supporting the conclusions of this article are included within article (and its Additional files 1, 2, 3, and 4).

Author information

Authors and Affiliations

Authors

Contributions

MDS created and implemented the algorithm, as well as analyzing and interpreting the results. The study was designed, directed and coordinated by PA. Validation of promoter predictions was designed by EP and CP. The manuscript was drafted by MDS and PA. The article was critical revised by AT, EP, MF and CP. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Pietro Alifano.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare 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:

Statistics of predicted promoters by G4PromFinder, PromPredict, PePPER and bTSSfinder in the whole genomes of S. coelicolor A3(2) and P. aeruginosa PA14. (DOCX 13 kb)

Additional file 2:

Coordinates of the validated promoters identified by G4PromFinder in S. coelicolor A3(2). (XLSX 100 kb)

Additional file 3:

Coordinates of the validated promoters identified by G4PromFinder in P. aeruginosa PA14. (XLSX 63 kb)

Additional file 4:

G4PromFinder algorithm. (PY 8 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

Di Salvo, M., Pinatel, E., Talà, A. et al. G4PromFinder: an algorithm for predicting transcription promoters in GC-rich bacterial genomes based on AT-rich elements and G-quadruplex motifs. BMC Bioinformatics 19, 36 (2018). https://0-doi-org.brum.beds.ac.uk/10.1186/s12859-018-2049-x

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://0-doi-org.brum.beds.ac.uk/10.1186/s12859-018-2049-x

Keywords