Skip to main content

Association of purine asymmetry, strand-biased gene distribution and PolC within Firmicutes and beyond: a new appraisal

Abstract

Background

The Firmicutes often possess three conspicuous genome features: marked Purine Asymmetry (PAS) across two strands of replication, Strand-biased Gene Distribution (SGD) and presence of two isoforms of DNA polymerase III alpha subunit, PolC and DnaE. Despite considerable research efforts, it is not clear whether the co-existence of PAS, PolC and/or SGD is an essential and exclusive characteristic of the Firmicutes. The nature of correlations, if any, between these three features within and beyond the lineages of Firmicutes has also remained elusive. The present study has been designed to address these issues.

Results

A large-scale analysis of diverse bacterial genomes indicates that PAS, PolC and SGD are neither essential nor exclusive features of the Firmicutes. PolC prevails in four bacterial phyla: Firmicutes, Fusobacteria, Tenericutes and Thermotogae, while PAS occurs only in subsets of Firmicutes, Fusobacteria and Tenericutes. There are five major compositional trends in Firmicutes: (I) an explicit PAS or G + A-dominance along the entire leading strand (II) only G-dominance in the leading strand, (III) alternate stretches of purine-rich and pyrimidine-rich sequences, (IV) G + T dominance along the leading strand, and (V) no identifiable patterns in base usage. Presence of strong SGD has been observed not only in genomes having PAS, but also in genomes with G-dominance along their leading strands – an observation that defies the notion of co-occurrence of PAS and SGD in Firmicutes. The PolC-containing non-Firmicutes organisms often have alternate stretches of R-dominant and Y-dominant sequences along their genomes and most of them show relatively weak, but significant SGD. Firmicutes having G + A-dominance or G-dominance along LeS usually show distinct base usage patterns in three codon sites of genes. Probable molecular mechanisms that might have incurred such usage patterns have been proposed.

Conclusion

Co-occurrence of PAS, strong SGD and PolC should not be regarded as a genome signature of the Firmicutes. Presence of PAS in a species may warrant PolC and strong SGD, but PolC and/or SGD not necessarily implies PAS.

Background

Three conspicuous genome features often co-occur in the Firmicutes. These are: (i) a pronounced Purine Asymmetry (PAS) with the dominance of purine bases (R = G/A) over pyrimidines (Y = C/T) along the entire leading strand of replication [1, 2], (ii) a strong Strand-specific bias in Gene Distribution (SGD), i.e., the presence of significantly larger population of genes, especially the essential and highly expressed ones, in the leading strand (LeS), as compared to that in the respective lagging strand (LaS) [35] and (iii) presence of two different isoforms of DNA polymerase III (PolIII) alpha subunit, PolC and DnaE, that are responsible for the synthesis of the LeS and LaS respectively [1, 3, 6]. Among these, the feature of SGD is not limited to the Firmicutes only. It exists in a large number of bacteria from diverse lineages, but the bias is the strongest in Firmicutes [3], reaching even 87% in some of its members such as Thermoanaerobacter tengcongensis [1, 7].

The other two genome features, PAS and PolC are believed to be the signature of the Firmicutes only [1]. Some stray cases of the existence of PolC in Fusobacteria and Themotogae were reported earlier [8], but these were taken as putative outcome of lateral gene transfer. Existence of PAS or G + A-dominance in LeS in any non-Firmicutes species is yet to be reported, though dominance of guanine along LeS is a common trait in bacteria [3, 9]. Earlier studies on Firmicutes attributed PAS to several factors [1, 1013]. A selection pressure exerted by PolC is believed to be the major contributor [11, 13]. Other plausible factors that might be responsible for PAS include an affinity in the genes to be co-oriented with the replicating fork [12], selective avoidance of stop codons and underrepresentation of costly amino acids [10]. A correlation between PAS and SGD might also exist [1]. It is worth mentioning at this point that a different type of strand-specific compositional bias - an enrichment of guanine and thymine (G + T) in the LeSs – has earlier been observed in many non-Firmicutes bacterial species [1416]. This trait, which is more frequent among the strictly host-associated endosymbionts or pathogens with reduced genomes [1720], has been attributed to the strand-biased deamination and 5-methylation of cytosine [9, 21].

All the studies on PAS, PolC and SGD reported so far, however, suffer from certain limitations. Some of these reports were based on limited number of genomes. For instance, the study proposing potential correlations between PAS, PolC and SGD [1] relied on a comparative analysis of only two model examples of Firmicutes and non-Firmicutes - Bacillus anthracis str. Ames 0581 and Francisella tularensis respectively. One may, however, argue whether the observations made in the study should be extrapolated to the entire bacterial kingdom or not. There were some large scale studies on strand-specific asymmetries in nucleotide composition and gene distribution in Firmicutes, which focused on the average biases in sequence composition at the whole genome levels [2, 12, 13]. However, none of these studies mentioned whether such global asymmetries also persist locally at smaller scales along the LeS or LaS of the respective genomes. There was also an effort towards the analysis of inter-strand variations in amino acid and codon usage in three DnaE-based groups of bacteria [2], but it focused only on the overall compositional features of those three groups. Additionally, the study did not pay attention to the preservence of the three features - PolC, PAS and SGD across the members within a group, especially when they thrive at diverse ecological conditions.

Studies on the Firmicutes, therefore, have left some pertinent questions unaddressed. Is PAS or G + A-dominance really an essential as well as exclusive feature of the Firmicutes? Do the usages of both guanine and adenine individually contribute to PAS across the whole genomes of Firmicutes species? Does the trait of PAS persist at local levels along all the LeS sequences of the Firmicutes? If yes, how does it influence the nucleotide usages in synonymous and non-synonymous codon sites of genes? Do PAS, PolC or SGD always co-occur in a bacterial genome? If not, how do they correlate with one another? In an attempt to address all these enigmatic issues, we have examined the status of PAS, SGD & PolC in diverse bacterial species (selected in a way to cover different genera of the phylum Firmicutes as well as other non-Firmicutes phyla of the bacterial world).

Our analysis reveals that co-existence of PAS, PolC and SGD is neither exclusive nor essential signature of the Firmicutes. These features co-exist only in a subset of the Firmicutes and also occur, either collectively or individually, in members of three other bacterial phyla - Fusobacteria, Tenericutes and Thermotogae. Almost all Firmicutes species contain PolC, but the usage of guanine and that of adenine do not always contribute individually to PAS across their whole genomes. A large number of Firmicutes members show the dominance of only guanine, but not of adenine, along their LeSs. Existence of some other trends like G + T dominance along LeS or presence of alternate segments of R and Y rich sequences along the genomes have also been observed. The study indicates that PAS might assure the presence of PolC and SGD, but the reverse is not true.

Results

PAS is neither an exclusive nor an essential feature of the Firmicutes

With a view to examine the status of PAS within and beyond the Firmicutes lineage, variations in local GC-skew and AT-skew values (averaged over 10 kb segments along the plus strands of the respective genomes) were studied in each of the organisms under study (Additional file 1: Table S1 and Additional file 2: Table S2). These skew trajectories may be classified into five distinct trends, as described in the Methods section. Some model examples of these five different trends in skew trajectories have been presented in Figures 1, 2, 3 and 4. In order to rule out any ambiguity while identifying such trends in skew trajectories, we have also examined the scatter plots of the local GC-skew and AT-skew values for each species under study. Some representative examples of such scatter plots are shown in Figure 5I-V.

Figure 1
figure 1

Instantaneous GC-skew (blue lines) and AT-skew (red lines) trajectories in model representatives of Trend I. (A) Bacillus anthracis str.Ames, (B) Listeria monocytogenes 07PF0776, (C) Staphylococcus aureus 04–02981, (D) Enterococcus faecalis V583, (E) Clostridium difficile CD196, (F) Thermoanaerobacter tengcongensis MB4, (G) Streptobacillus moniliformis DSM 12112, (H) Illyobacter polytropus DSM 2926.

Figure 2
figure 2

Instantaneous GC-skew (blue lines) and AT-skew (red lines) trajectories in model representatives of Trend II. (A) Streptococcus agalactiae NEM 316, (B) Acidaminococcus intestini RyC-MR95, (C) Geobacillus kaustophilus HTA426, (D) Veillonella parvula DSM 2008, (E) Thermodesulfobium narugense DSM 14796, (F) Clostridiales genomosp BVAB3 UPII9 5, (G) Acinetobacter sp. ADP1, (H) Candidatus Protochlamydia amoebophila UWE25.

Figure 3
figure 3

Instantaneous GC-skew (blue lines) and AT-skew (red lines) trajectories in model representatives of Trend III. (A) Ruminocococcus albus 7, (B) Fusobacterium nucleatum subsp. nucleatum ATCC 25586, (C) Leptotrichia buccalis C-1013-b, (D) Mycoplasma mycoides SC PG1, (E) Thermotoga maritima MSB8, (F) Aquifex aeolicus VF5.

Figure 4
figure 4

Instantaneous GC-skew (blue lines) and AT-skew (red lines) trajectories in model representatives of Trend IV & V. Trend IV - (A) Oenococcus oeni PSU1, (B) Sulfobacillus acidophilus DSM 10332, (C) Mycobacterium tuberculosis CDC 1511, (D) Bartonella henselae str. Houston-1. Trend V - (E) Mycoplasma synoviae 53, (F) Acidobacterium capsulatum ATCC 51196.

Figure 5
figure 5

Scatter plots of Local GC-skew and AT-skew values in model representatives of organisms following different trends in purine usages. (I) Trend I - Bacillus anthracis str.Ames (L) and Streptobacillus moniliformis DSM 12112 (R); (II) Trend II - Geobacillus kaustophilus HTA 426 (L) and Acinetobacter sp. ADP1 (R); (III) Trend III - Ruminocococcus albus 7 (L) and Fusobacterium nucleatum subsp.nucleatum ATCC 25586 (R); (IV) Trend IV - Oenococcus oeni PSU1 (L) and Bacteroides fragilis 638R (R); (V) Trend V - Mycoplasma synoviae 53 (L) and Acidobacterium capsulatum ATCC 51196 (R).

Trend I – Explicit PAS with individual dominance of G and A along the entire LeS

Trend I refers to the cases, where both the purine bases (guanine and adenine) individually contribute to the purine-richness of the LeSs. Some typical examples of Trend I species are shown in Figure 1A–H, where local GC-skew and AT-skew values are, by and large, positive between the putative origin (oriC) and termination (ter) sites of replication along the plus strand, and negative in the other half of the genomes; with a sharp transition from the positive to negative values at ter (Figure 1). In most of the Trend I organisms, more than 70% of the 10 kb LeS segments have exceptionally high frequencies of both guanine and adenine as compared to cytosine and thymine respectively (Tables 1 and 2), while the number of LeS segments of other three possible combinations (b), (c) or (d) are significantly low in most cases. These observations indicate that the LeS sequences have explicit enrichment of both the purine bases (guanine and adenine) in all the organisms of Trend I. We shall henceforth refer to this trend as explicit PAS or simply PAS.

Table 1 Status of combinations (a) – (d), PAS, SGD and PolC in Firmicutes taken in this study
Table 2 Status of combinations (a) – (d), PAS, SGD and PolC in the non-Firmicutes organisms examined in this study

Presence of Trend I are found in more than 70% of the Firmicutes under study and it is predominant among the members of Bacillales, especially in those belonging to the genera of Bacillus, Listeria, Staphylococcus, Enterococcus and Thermoanaerobcter (Figure 1A-D, F). However, Bacillus is the only genus among Firmicutes, all members of which show predominance of both guanine and adenine along the LeS. Trend I has been observed in some members of Clostridia also (Figure 1E).

Interestingly enough, Trend I is not confined to the lineage of Firmicutes only. It has also been observed in some Fusobacteria and Tennericutes. Of the five Fusobacteria and twelve Tenericutes species studied (Additional file 2: Table S2), three Fusobacteria including S. moniliformis (Figure 1G), I. polytropus (Figure 1H) and five Tenericutes (Table 2) display explicit PAS, indicating that PAS is not an exclusive characteristic of Firmicutes only.

None of the non-Firmicutes, non-Fusobacteria and non-Tenericutes organisms under study exhibited unequivocal G + A-enrichment of LeS. It suggests that the presence of PAS might be confined only to the three bacterial phyla, Firmicutes, Fusobacteria and Tenericutes, which are thought to be closely related from the evolutionary point of view [8].

Trend II – Only G-dominance in LeS with no unequivocal trend in adenine usage

All Firmicutes genera except Bacillus include certain members, which show dominance of only guanine, but not of adenine along the LeS. This trend (Trend II) has also been observed in a number of non-Firmicutes species from diverse bacterial phyla. Some model examples of Trend II have been depicted in Figure 2, where the representatives of Firmicutes are shown in Figures 2A-F, and those of non-Firmicutes in Figure 2G-H. In all cases, the GC-skew trajectory exhibits a sharp transition in sign only once at the oriC/ter region, but AT-skew values undergo irregular oscillation around the null axis, showing no definite pattern. Cumulative GC and AT-skew trajectories and instantaneous RY skew values of the respective species are shown in Additional file 3: Figure S1. As expected, the cumulative GC-skew always increases between oriC and ter and decreases along the other half of the plus strand. But the nature of the cumulative AT-skew varies from species to species and in majority of the organisms following Trend II, hardly deviating from the null value (Additional file 3: Figure S1 CL, DL, EL, GL & HL). In all Firmicutes members of this category, the magnitude of GC-skew values is usually much higher than the respective AT-skew values. Hence the average local purine-content of LeS sequences remain higher than the respective pyrimidine content (Additional file 3: Figure S1), but the total contribution to such apparent purine-richness of LeS comes from the G-dominance only with little or no contribution from the adenine frequencies. However, in certain Trend II Firmicutes, the overall R- usage does not follow any definite strand-specific pattern (Additional file 3: Figure S1).

The differences between PAS (Trend I) and G-dominance (Trend II) can be clearly understood from Figure 5. In organisms having PAS (Trend I, Figure 5 IL, IR), the points from the segments between oriC and ter (blue points) usually lie in the first quadrant (barring a few exceptions). It re-confirms that both GC-skew and AT-skew values are in general positive. The points from the segments between ter and oriC (red) lie in the third quadrants indicating negative values for both the skews. On the contrary in Trend II organisms, the points corresponding to the LeS part of the plus strand are almost equally distributed in first and fourth quadrants (Figure 5, IIL, IIR, blue points), while those corresponding to the LaS parts (red points) are distributed among the second and third quadrants (red points). This indicates that the GC-skew values remain mostly positive along LeS and negative along LaS, but the AT-skew values fluctuates between positive and negative values along both the replicating strands. Fluctuations in AT-skew magnitudes along two replication strands of Trend II organisms are also apparent from Tables 1 and 2 - clearly indicating that in organisms following Trend II, frequencies of LeS segments with base usage combinations (a) and (b) both are significantly high and their values are often comparable to one another. However the presence of the other two combinations (c) and (d) are negligible, in general.Another distinct feature of Trend I is that the pairs of instantaneous GC-skew and AT-skew values exhibit significant positive correlations for both oriC-ter (blue points) and ter-oriC (red points) regions along the plus strand (Figure 5, IL, IR). In cases of Trend II (Figure 5, IIL, IIR), no significant positive correlations exist in general between the pairs of GC-skew and AT-skew values. Even if it exists, the magnitudes of the correlation coefficients are not as high as those observed in Trend I (Figure 5, IL, IR). All these observations clearly indicate that in organisms following Trend I, usages of guanine and adenine both contribute significantly to PAS. In Trend II organisms, only an apparent purine-richness often prevails along the LeS, where the sole contribution to purine enrichment comes from the G-dominance only, with the adenine usage hardly playing any role.

Within the Firmicutes phylum, Trend II prevails in the non-Bacilli classes like Clostridia or Negativicutes, along with certain Bacilli genera like Streptococcus, Geobacillus or Lactobacillus etc. A small number of exceptions from the order Bacillales also fall under this category.

Trend III - Presence of alternate stretches of R-dominant and Y-dominant sequences along both the replicating strands

There is one Firmicutes species, Ruminococous albus, which exhibits a conspicuous trend of purine usage (Trend III). In this species, instantaneous GC-skew and AT-skew trajectories toggle their signs frequently and simultaneously in a way such that the respective GC and AT-skew values remain, in most cases, of the same sign (Figure 3A). Though it shows an overrepresentation of R-dominant stretches (combination (a) ≈ 53%), the Y-dominant stretches also occurs with random frequency (combination (d) ≈ 22%) [Table 1]. This suggests that a major part of the genome of R. albus is comprised of alternate purine-rich and pyrimidine-rich segments. A similar trend is also observed in two Fusobacteria species, namely Fusobacterium nucleatum (Figure 3B) and Leptotrichia buccalis (Figure 3C). Majority of the Tenericutes members examined in the study, including certain Mycoplasma and Ureplasma species, also follow Trend III (Figure 3D).In organisms following Trend III, local GC-skew and AT-skew values bear strong positive correlations (Figure 5, IIIL, IIIR), as observed earlier in Trend I. However, there is a major difference between the scatter plots in two trends. In Trend I, points corresponding to LeS (blue) and LaS (red) parts of the plus strands are segregated in the first and third quadrants respectively. In Trend II organisms, on the contrary, points from both the LeS and LaS sequences are distributed uniformly in the first and third quadrants, implying that both guanine and adenine frequencies are oscillating simultaneously between positive and negative values along the replicating strands.

The presence of alternate genomic segments of R-rich and Y-rich sequences was reported earlier for thermophilic/hyperthermophilic bacteria [22]. A number of thermophiles in the current dataset, especially those belonging to the Aquificae and Thermotogae lineages show the presence of Trend III in their genomes (Table 2). Two typical examples of such thermophilic organisms Thermotoga maritima and Aquifex aeolicus are presented in Figure 3E and F. The amplitudes of purine-rich/pyrimidine-rich segments of the genomes are, in general, much smaller (Figure 3A-D), but the percentage occurrence of such segments are much higher in thermophiles, as compared to the Trend III Firmicutes, Fusobacteria or Tenericutes (Tables 1 and 2). It is worth mentioning at this point that all thermophiles/hyperthermophiles does not exhibit Trend III. A substantial part of them follow a distinct trend of G + T-enrichment along LeS (Trend IV) as described below.

Trend IV - G + T dominance along the leading strands

In majority of the bacteria from non-Firmicutes, non-Fusobacteria, non-Tenericutes, non-Aquificae and non-Thermotogae lineages, a strand specific bias exists not in favour of G + A, but in favour of G + T usage along the entire LeS (Trend IV). Organisms following Trend IV include Proteobacteria, Actinobacteria, Bacteroides, Chloroflexi, Planctomycetes, Spirochetes etc. (Table 2). Two model examples of Trend IV genomes are shown in Figure 4(C, D), where the signs of GC-skew and AT-skew trajectories are of opposite signs. Both the skew trajectories change their signs simultaneously at oriC/ter regions, so that their LeSs have, in general, an over representation of guanine and thymine, as reported earlier [1720]. In free living organisms, the magnitudes of the instantaneous GC-skew and AT-skew values are often quite low (Figure 4A). However, in obligatory intracellular microbes undergoing genome reduction, both GC-skew and AT-skew values are, in general, of significantly higher magnitudes confirming the general notion of their parasitic adaptation [1720].

Though quite common among other bacteria, Trend IV is rarely seen within the Firmicutes. Among 102 Firmicutes in the dataset, only two organisms seem to follow Trend IV. These include Oenococcus oeni - a Lactobacillales species and Sulfobacillus acidophilus - a Clostridiales member (Table 1). Some typical examples of the scatter plot of local GC-skew and AT-skew values in organisms following Trend IV are shown in Figure 5 IVL and IVR. As expected, most of points from the LeS portion of the plus strand lie in the fourth quadrants (since GC-skews are positive and AT-skews are negative), but those from the LaS regions mostly appear in the second quadrants (as GC-skews are negative and AT-skews are positive, in most cases).

Trend V - No identifiable pattern in base usage

Lastly, there are a small number of bacterial genomes displaying random oscillation around the abscissa in both GC-skew and AT-skew trajectories. In these cases, no general trend can be detected either in the signs of GC-skew/AT-skew values or in the distribution of 10 kb segments among four combinations (a)–(d) (Table 2). Certain Tenericutes, Acidobacteria, Actinobacteria, Cyanobacteria etc. show ambiguous behavior in their GC-skew and AT-skew values (Figure 4E and F, Table 2). As expected, points in the scatter plots of GC and AT-skew values (Figure 5, VL, and VR) are also randomly distributed in all four quadrants, having no definite pattern or correlations.

PAS, SGD and PolC might not bear any definite correlation in Firmicutes or other bacteria

As indicated in the present analysis, PAS exists in a substantial fraction of the Firmicutes but it is not a signature trait of this phylum. On the other hand, there are certain Fusobacteria and Tenericutes that clearly show the presence of PAS. In view of a recent hypothesis in favor of a correlation between PAS and SGD, it will be intriguing to examine the correspondence between PAS, PolC and SGD in Firmicutes, Fusobacteria and other organisms under study. To this end, we have checked the status of SGD as well as of PolC across all bacterial species of our dataset. Outcomes of the study are provided in Tables 1 and 2. As can be seen from these files, all organisms having PAS (Trend I) show very strong SGD. If we consider the number of 10 kb segments with G > C and A > T as a measure of the strength of PAS in a Trend I organism (Tables 1 and 2), then the scattered plot of PAS and SGD shows a strong positive correlation between themselves, the correlation coefficient being 0.59 (the scattered plot not shown).

PolC is found to be present in almost all Firmicutes members as well as in all Fusobacteria, Tenericutes and Thermotogae members under study. There are only two exceptions – Desulfotomaculum acetoxidans and Thermaerobacter marianensis both belonging to the class Clostridia under the Firmicutes phylum. D. acetoxidans and T marianensis both possess marked SGD but no PAS. BLASTP search for PolC homolog could not detect the presence of PolC in these two organisms.

All PolC-containing Firmicutes, Fusobacteria and Tenericutes have shown statistically significant SGD, irrespective of the trends in their nucleotide usages (Table 1). PolC are also present in Thermotogae members, but they do not possess PAS. In most cases, they have alternate R and Y-dominant stretches along their genome sequences (Trend III, Table 2). Our analysis shows that five out of seven Thermotogae species do not display any significant SGD. On the contrary, a large fraction of non-PolC organisms following Trend IV (i.e., G + T-dominance along LeS) have shown significant SGD – an observation that comply with earlier reports [16, 21]. These observations re-confirm that the presence of PolC is neither a necessary nor a sufficient condition for SGD in bacteria.

The strength of SGD varies appreciably in organisms with different trends in nucleotide usages along their LeS/LaS, as can be seen from their SGD distribution profiles (Figure 6A) as well as from the individual SGD values (Table 1). Interestingly enough, the major peaks of the SGD distribution profiles of PolC-Trend I and PolC-Trend II organisms fall in the same range (~0.8) (Figure 6A), while the SGD profiles of the PolC-containing Trend III organisms, non-PolC Trend II organisms and non-PolC Trend IV organisms - all display peaks in the range of 0.55-0.6. In both Trend I and PolC-Trend II categories, SGD is greater than 0.7 for majority of the organisms in the dataset (Tables 1 and 2, Figure 6A). The only difference between two profiles is that in case of PolC-Trend II, there are a few genomes having SGD distribution profiles < 0.65, which could not be found in case of Trend I (Figure 6A). This observation indicates that organisms with only G-dominance may have relatively low SGD in some cases, but organisms showing explicit G + A-dominance are always characterized by a strong bias in gene orientation along replication direction. The strong resemblance between the distribution profile of Trend I organisms, characterized by PAS (and PolC) and that of PolC-Trend II organisms having G-dominance suggests that PAS asserts SGD, but SGD does not warrant PAS. For instance, the PolC-containing Thermoanaerobacterales species Ammonifex degensii KC4 or Selenomonadales species Veillonella parvula DSM 2008 do not show explicit PAS, but have extremely high SGD (>80% genes in LeS) (Table 1).The number of organisms in PolC-Trend III group is too low (one Firmicutes and twelve non-Firmicutes members) to provide any statistically significant pattern. Nevertheless, it is intriguing to find that the major peak of its SGD distribution profile comes in the same range as that of the non-PolC-Trend IV population. These distribution profiles give a hint that the average SGD of PolC-Trend III (and also of non-PolC-Trend II/non-PolC-Trend IV) organisms might not be as high as in cases of Trend I or PolC-Trend II (Figure 6A). In order to gain a conclusive picture on SGD profiles of PolC-Trend III genomes, one must wait for availability of complete genome sequence information for more number of species belonging to this category. Distribution profiles have not been plotted for PolC-Trend IV or Trend V organisms, since the current dataset contains only three organisms in Trend IV and four organisms in the Trend V categories.

Figure 6
figure 6

SGD and Genomic GC-content distribution profiles in organisms showing different trends. (A) Distribution of SGD in Trend I (red solid line), PolC-Trend II (blue solid line), non-PolC-Trend II (blue dotted line), PolC-Trend III (green solid line) and non-PolC Trend IV (violet dotted line) organisms; (B) Distribution of Genomic GC-content for Trend I and Trend II organisms.

Distinct trends in base usage in three codon sites of Les and LaS genes and intergenic regions in Trend I and Trend II Firmicutes

On the basis of strand-specific sequence composition, Firmicutes members may broadly be classified in two major categories: 1) the ones with G + A-dominance or PAS in LeS (Trend I) and 2) those having only G-dominance in LeS with no definite strand-specific bias in adenine usage (Trend II). There are some exceptions like R. albus or O. onei showing other conspicuous patterns in base usage (Trend III or Trend IV), but they are very few in number. Analysis of the distribution patterns of average genomic GC content of Trend I and Trend II organisms showed that the average GC-contents of Trend I organisms are usually significantly less than 50%, while the GC-content of Trend II genomes vary in much broader range (35 – 80%) (Figure 6B). It is not clear whether the relatively lower GC-content of the Trend I genomes could anyway be associated with PAS. This observation inspired us to further probe into the base usage patterns in three different codon sites of genes in two replicating strands of the Trend I and Trend II Firmicutes members of the current dataset. Figures 7, 8 and 9 represent three typical examples of the outcomes of this study. Figure 7 represents the trends in base usage in three individual codon sites and intergenic regions as well as in overall coding regions for all annotated genes in LeS (left panels) and LaS (right panels) of S. aureus. The organism is a typical representative of Trend I Firmicutes. Figures 8 and 9 depict the base usage patterns in S. agalactiae and G. kaustrophilus - two model representatives of Trend II Firmicutes with low and relatively high genomic GC-contents (35.6% and 52% respectively). Among the 102 Firmicutes species examined, only three species exhibited Trend IV. It is difficult to say whether the patterns observed in these three organisms typically represent the general trends in base usages within the PolC-containing Trend IV species of similar genomic G + C-content. Nevertheless, the base usage patterns in O. onei are shown in Figure 10 as a representative of these three species. The base usage in E. coli and B. henselae genes are depicted in Additional file 4: Figure S2 and Additional file 5: Figure S3 respectively, as the representatives of non-PolC organisms. E. coli represents Trend II non-PolC species, while B. henselae exemplifies Trend IV non-PolC organisms. There are usually no distinct strand-specific divergences in nucleotide usages in genes of Trend III or Trend V organisms (data not shown).

Figure 7
figure 7

Trends in individual base usages in Staphylococcus aureus 04–02981 for genes encoded by both LeS and LaS. Subscripts 1, 2, 3 indicate the percentage of occurrences of the respective base at 1st ( A , B ), 2nd ( C , D ) and 3rd ( E , F ) codon sites, intergenic indicate the percentage of intergenic regions ( G , H ) and the subscript T stands for the total percentage ( I , J ) of occurrence of the base in individual genes of the organism.

Figure 8
figure 8

Trends in individual base usages in Streptococcus agalactiae NEM316 for genes encoded by both LeS and LaS. Subscripts are same as in Figure 7.

Figure 9
figure 9

Trends in individual base usages in Geobacillus kaustophilus HTA426 for genes encoded by both LeS and LaS. Subscripts are same as in Figure 7.

Figure 10
figure 10

Trends in individual base usages in Oenococcus oeni PSU 1 for genes encoded by both LeS and LaS. Subscripts are same as in Figure 7.

As revealed in Figures 7, 8, 9 and 10 and Additional file 4: Figure S2 and Additional file 5: Figure S3, there are some common features in base usages in organisms in general irrespective of their compositional trends. For instance, in most of the cases, G1 > C1 and A1 > T1, while G2 < C2 and A2 ≥ T2 in both LeS and LaS genes - an observation that conform with the existing notion of the universal three-base periodical pattern (G-non-G-N) of mRNA sequences [23]. Inter-group differences in base preferences are more apparent in the third codon sites of both LeS and LaS genes. There are some general patterns observed in 3rd codon sites of genes in PolC-containing organisms following Trend I - Trend III, as given below,

In Trend I species:

A3 ~ T3 > G3 > C3 (LeS genes), T3 ≥ A3 > C3 > G3 (LaS genes)

where, N3 indicates the average frequency of the nucleotide N in the 3rd codon sites of genes in the respective strands of the species under study.

In A + T-rich Trend II species:

T3 > A3 > G3 > C3 (LeS genes), T3 ≥ A3 > C3 > G3 (LaS genes)

In G + C-rich Trend II species:

G3 ≥ C3 > T3 > A3 (LeS genes), C3 ≥ G3 > T3 ~ A3 (LaS genes)

In O. onei, which represents the group of Trend IV organisms, especially of the A + T-rich ones:

T3 > A3 > G3 > C3 (LeS genes), T3 ≥ A3 > C3 ≥ G3 (LaS genes)

As shown in Additional file 4: Figure S2 and Additional file 5: Figure S3, trends in 3rd codon sites base usages in non-PolC organisms (both Trend II and Trend IV) are, by and large, similar to those observed in the PolC-containing Trend II organisms of similar G + C-bias, though the actual frequencies of different bases vary from one species to another.

In intergenic regions, usages of A and T are usually higher than those of G and C in most of the organisms (except in some highly G + C-rich organisms, where usages of A or T are comparable to usage of G or C). It was expected because of the presence of A + T-rich promoter sequences (TATA box etc.) in intergenic regions. Nevertheless, some specific biases in the base usages in the intergenic regions could be observed. For instance, in Trend I organisms, Aintergenic ~ Tintergenic along the LeS, but Tintergenic ≥ Aintergenic in LaS. This pattern is similar to that observed in the 3rd codon sites of the respective species. Furthermore, in most of the species, Gintergenic > Cintergenic along LeS, but Cintergenic > Gintergenic along LaS - a pattern observed in the 3rd codon sites the genic regions of the bacteria, in general, irrespective of their trends in base usages (Figures 7, 8, 9 and 10).

The overall base frequencies follow the trends, as given below.

In Trend I species,

AT > TT > GT > CT (LeS genes), A T ≥ TT > CT ≥ GT (LaS genes)

In A + T-rich Trend II species:

AT ~ TT > GT > CT (LeS genes), AT ~ TT > C T > GT (LaS genes)

In G + C-rich Trend II species,

GT > CT > AT ~ TT (LeS genes), CT ≥ GT > TT ~ AT (LaS genes)

In O. onei (Trend IV),

TT ≥ AT > GT > CT (LeS genes) AT > TT > GT ~ CT (LaS genes)These trends are in complete agreement with the GC-skew and AT-skew trajectories shown in Figures 1 and 2. Needless to say, a finite number of genes in each organism under study stand out as exceptions.

At a first glance, it may appear that base usage patterns in non-synonymous sites are quite similar across the two replicating strands of a particular species. However, a careful examination reveals some subtle differences. For instance, G1 in LeS genes is, in general, significantly higher than that in LaS genes of the same organism. On the contrary, C1 is, significantly lower in LeS genes as compared to that in LaS genes (in many cases, but not in all) (data not shown). Appreciable cross-strand differences in nucleotide selection have also been observed in the second codon sites of genes in a substantial number of PolC-containing organisms of the dataset (data not shown). The most prominent cross-strand difference in base usage is the preference for G over C by LeS genes and for C over G by LaS genes at their third codon sites (C3 ~ G3 in LaS genes in some cases, especially in GC-rich organisms).

Discussion

The present study examines the status of PAS, SGD & PolC in Firmicutes and other bacterial species from diverse lineages. Co-existence of PAS, SGD and PolC in Firmicutes has earlier been reported by various investigators and several molecular mechanisms have been put forward as plausible explanations of this co-existence [1, 6, 10, 12]. Among these, the most accepted hypothesis is that the R-richness on the LeS and R-poorness on the LaS might be a type of sequence signature of the heterodimeric DNA polymerase III alpha subunit in Firmicutes [24]. It was also proposed that the presence of PolC might have exerted a selection pressure in favour of R-enrichment in LeS in order to prevent nonspecific RNA–RNA interactions and formation of excessive double-stranded RNA [22]. This, in turn, has led to the emergence of a strong SGD through preferential localization of R-rich genes in LeS during random genetic exchange across two strands [25]. On contrary to these existing notions, the present analysis clearly demonstrates that PAS or G + A-dominance in LeS is neither an essential feature of the Firmicutes, nor a sequence signature of PolC and/or SGD. It exists only in a subset of the Firmicutes, especially in those belong to the order Bacillales. There are an appreciable number of non-Bacillales Firmicutes (e.g., Streptococcus, Geobacillus or Lactobacillus), which contain PolC and have strong SGD.They do not show any definite strand-specific bias in their adenine usage patterns. In most of these Firmicutes, the cumulative R-content is significantly higher in the LeS than that in the LaS, but the sole contribution to R-asymmetry comes from the guanine bias, with little or no role of the adenine content. There is also a Firmicutes species R. albus that despite having PolC does not show strand-specific purine asymmetry. It rather contains alternate stretches of R-rich and Y-rich segments. Certain Firmicutes also exhibit G + T-dominance in their LeS sequences. It may therefore be said that PAS is not an essential feature of Firmicutes.

PAS is not an exclusive characteristic of the Firmicutes either. It has been observed in some Fusobacteria and Tenericutes species also. Among five Fusobacteria under study, three organisms namely S. moniliformis, I. polytropus and S. termitidis, exhibit strong PAS and strong SGD. The other two Fusobacteria members have alternate stretches of R-rich and R-poor regions along both the strands of replication, though all five members of the phylum possess PolC. Similarly, among twelve PolC-containing Tenericutes members of the dataset (Table 2), five species display strong PAS as well as highly significant SGD.

Observations made in the present study also suggest that the existence of PAS or G + A-richness of LeS is usually associated with PolC and a strong SGD, but the reverse may not be true. There are four bacterial phyla, namely Firmicutes, Fusobacteria, Tenericutes and Thermotogae, members of which contain PolC. Among these, PAS or G + A-richness of LeS prevails only in a certain fraction of Firmicutes and in three Fusobacteria, all of which carry PolC and almost all of which show strong SGD. However, there are a number of non-PAS Firmicutes, especially the ones exhibiting Trend II which also display equally strong SGD. It is therefore suggested that presence of a strong SGD does not necessarily imply PAS.

It was proposed earlier that PolC might play a role in maintenance of SGD in Firmicutes. The present study concords with this notion in the sense that majority of the PolC-containing genomes have significant SGD. However, the presence of PolC alone might not lead to a strong SGD (>70%). Most of the Trend III Firmicutes, Fusobacteria and Tenricutes members examined so far have shown relatively weak SGD (<70%). Interestingly enough, three Firmicutes species B. tusiae, O. oeni, S, acidophilus, having strong G + T dominance along their LeSs, exhibit the presence of strong SGD. It is, therefore, tempting to postulate that it might not be PolC alone, but a coupling between PolC and the G-dominance in LeS that has led to a strong SGD in the Firmicutes/Fusobacteria. Again, there are some exceptions. Two Clostridial species, T. marianensis and D. acetoxidans have SGD, but not PAS and PolC. It is intriguing to note that all Thermotogae members possess PolC and follow Trend III, but do not have any significant SGD. This observation indicates that the suggested correlation between PolC and SGD did not hold well in Thermotogae.

A comparison of the trends in base usages within different codon sites in PolC–containing Firmicutes (Figures 7, 8, 9 and 10) with those in non-PolC bacteria like E. coli (Additional file 4: Figure S2) or B. henselae (Additional file 5: Figure S3) reveals that the non-synonymous sites of genes follow certain general trends in most of these species; whereas the actual nucleotide frequencies vary from species to species depending on their average genomic GC-bias. However, a conspicuous trend that differentiates Trend I Firmicutes, Fusobacteria and Tenericutes from all other organisms; is similar or even higher usage of A3 as compared to that of T3 in LeS genes. It is in contrast to the earlier observations on preferences of pyrimidines over purines in third codon sites [26]. However, in all other organisms under study, usage of T3 is higher than that of A3 in LeS. These observations point to the existence of a unique selection pressure in Trend I Firmicutes in favour of adenine over thymine individually in all three codon sites, especially in the third ones. This unique feature of Trend I organisms seems to have a major contribution to the PAS.

Molecular processes that may incur strand-specific compositional biases in bacterial genomes include DNA replication, transcription coupled repair (TCR) [1, 3, 2729] and the process of deamination and 5-methylation of cytosine [9, 21]. When a gene is located on the leading strand of a PolC-containing species, the mutational bias at the replication level and the bias at the transcription level both tend to increase its G + A-content; but the process of cytosine methylation generates a LeS-wide bias towards increasing G + T-content. On the contrary, genes on the LaS experience a mutational bias towards increasing C + A-content during the replicational process, a bias in favour of increasing G + A-content during TCR as well as a bias towards increasing C + T-content owing to the cytosine methylation. The resultant base composition of the LeS/LaS genes would depend on the relative intensities of these biases in the respective species. If all three processes remain significantly active in a genome, their collective effect is expected to create an unequivocal dominance of G over C in LeS genes of the organisms, as observed in Figures 7, 8, 9 and 10. If the mutational biases during replication and/or transcription dominate over the deamination/methylation bias, the frequencies of A would be higher than T. Thus it is tempting to propose that this might be the cases in Trend I organisms (Figure 7). On the other hand, if the G + T-bias owing to cytosine deamination be strong enough to nullify or even outshine the G + A-bias of replication/transcription processes, the LeS genes might exhibit Trend II or even Trend IV traits. Similar arguments may also be put forward to explain the compositional skews of LaS genes in Figures 7, 8, 9 and 10. Reports on the presence of a high level of α/β-type small, acid-soluble spore proteins (SASPs) in Bacillus subtilis [30] and in many other members of the orders Bacillales and Clostridiales [31, 32] suppressing cytosine deamination to uracil in native DNA are in good agreement with our proposition. Future investigations on the status and activities of the α/β-type SASPs in Trend II and Trend IV, which is out of the scope of the present analysis, may help in further validation of this notion.

In the entire dataset, there are only two Firmicutes members, which are devoid of two conspicuous features of the phyla, i.e., PAS and PolC. Considering the fact that bacterial genomes are highly dynamic in nature and they are continuously undergoing the processes of gene loss and gene gain, one could presume that the gene encoding PolC had been lost from these two Firmicutes members. Hence they did not experience any selection pressure in favour of PAS. Presence of SGD in these two organisms re-affirms that the existence of PAS or PolC is not an essential pre-requisite of SGD.

Among the non-Firmicutes, existence of PolC was reported earlier in F. nucleatum and T. maritima as potential cases of horizontal gene transfer [8, 33]. The present analysis indicates that PolC is present not only in these two species, but it is also shared with all other Fusobacteria and Thermotogae members examined in this study. In fact, among all non-Firmicutes in the current dataset, presence of PolC could so far be detected in three lineages – Fusobacteria, Mollicutes or Tenericutes and Thermotogae. Surprisingly enough, most the members of these three lineages exhibit strong explicit PAS (both G- and A-dominance in LeS) or have alternate R- and Y-dominance along their genomes (with a few exceptions that exhibit Trend V). It would not therefore be irrational to presume that the presence of PolC and the emergence of R-rich/Y-rich genome segments in some of these organisms might have some common link. It may be mentioned in this context that some of the earlier evolutionary studies pointed towards a plausible close evolutionary relationship among Firmicutes, Fusobacteria and Mollicutes. The ribosomal molecular phylogeny and core genome contents of Fusobacteria members indicated that this lineage might have branched out at the base of Firmicutes.

Mollicutes were previously thought to be a class within Firmicutes, but later on the basis of their unique phenotypic properties such as the lack of rigid cell walls and other evidences, they have been placed under a new phylum called Tenericutes [34]. However, the phylogenetic analysis based on phosphoglycerate kinase (Pgk) amino acid sequences indicated a monophyletic origin of the Mollicutes within Firmicutes [35]. The same study also had placed Fusobacteria (and even Thermotogae) within the Firmicutes – an observation that completely conforms to the findings made in the present study. One cannot, therefore, rule out the possibility that the feature of PAS was not horizontally acquired by the Fusobacteria or Mollicutes, but inherited normally from their Firmicutes like ancestors. Some members of Fusobacteria like S. moniliformis, I. polytropus, are still bearing the ancestral signature of PAS in their LeS sequences. However, their fellow members and the Mollicutes species might have undergone a series of genome reshuffling, recombination and local strand reversal processes in course of their evolution. As a consequence, their original ancestral genome architecture with R-rich LeS and R-poor LaS might have gradually been turned into the present-day genome structures having a mosaic of alternate R-rich and R-poor segments along both the strands. These processes of genome reshuffling or recombination might have also altered the gene orientation along two replicating strands. It would have been intriguing to study the correlations, if any between the processes of genome reshuffling and the evolution of gene orientation. However, it is beyond the scope of the present analysis.

The organisms showing Trend III or Trend V often exhibit zig-zag patterns in their GC-skew and other skew curves and it sometimes becomes difficult to identify the ter regions of their chromosomes unambiguously. One may argue that in such cases, a random pattern in base usages along two strands (Trend V) may arise due to an error in assignment of the ter region and hence among the LeS and LaS sequences. With a view to check whether it is mere shift in the ter region or mixing up of ancestral LeS and LaS sequences owing to genomic recombination that may alter the basic trend in base usage along LeS and LaS sequences, we have examined the GC-skew and AT-skew patterns (Additional file 6: Figure S4) in eight Yerisina pestis strains, which are known for having undergone drastic changes in the relative positions and directions of discrete genome segments following extensive genomic rearrangements [36]. In all strains except Y. pestis Pestoides F, putative oriC have been found near the start point of the reported plus strand sequences and the putative ter point, despite having finite displacement along plus strand, appeared to be located close to the mid-point of the plus strand. In Y. pestis D182038 and Y. pestis biovar Microtus 91001 yielding zig-zag cumulative GC-skew curves with multiple extrema, putative ter points were determined from the extremum point closest to the point representing the putative oriC plus half of the chromosome length (as described in the Methods section). Y. pestis Pestoides F is the only strain, where the putative oriC and ter regions (as detected from the unique extremum point of cumulative GC-skew) both have shifted in an uneven manner and as a consequence, the distances between oriC and ter points along two strands become significantly different (Additional file 6: Figure S4, HR). All the predicted locations of oriC and ter regions conform well to the findings made earlier by Liang et al. (Figure three of [36]). Interestingly enough, seven out of eight strains unambiguously exhibit Trend IV (Additional file 6: Figure S4, left panel, Table S3) and these include even Y. pestis Pestoides F having asymmetric locations of oriC and ter along the plus strand and Y. pestis D182038 showing a zig-zag skew curve. The only exceptional case that displayed Trend V (Additional file 6: Table S3) is Y. pestis biovar Microtus 91001 – the strain exhibiting maximum number of genomic rearrangement – translocation and/or inversion of 54 out of 61 genome plates with respect to the Y. pestis CO92 genome, as reported in Figure 3 of Liang et al. [36]. This observation clearly indicated that it is neither an asymmetric location of oriC and ter regions, nor any ambiguity in the prediction of the ter point, but the specific types of genomic rearrangements leading to a substantial mixing up of LeS and LaS sequences that may result in a change in the trends in local base usages in bacterial genomes.

As already mentioned, the situation might have been quite different in case of Thermotogae. The exact position of Thermotogae within the tree of life is also not clear yet. Different markers have yielded varying results, which place Thermotogae and other hyperthermophiles like Aquificae either close to the root of the tree of life [37] or a little “up” from the root close to Fusobacteria [38] or to Bacillus and Mycoplasma species [39]. A significant degree of horizontal acquisition of genes by Thermotogae from other species, especially from archaea, has made the situation even more confusing. As already mentioned, the Pgk-based phylogeny, which was otherwise congruent with 16S rRNA data placed Thermotogales closer to Firmicutes than to any other phylum. In the light of all these studies, it may be said that there could be multiple events leading to the current architectures of Thermotogae genomes. PolC might have horizontally (or even vertically) acquired by an ancestral species prior to the branching of the lineage of Thermotogae and the current architecture of R-rich and R-poor segments of Thermotogae might be the relics of their ancestral PAS like sequence signature of the PolC. Alternately, considering the fact that Thermotogae are hyperthermophile in nature and that they are believed to be close enough to Aquificae, it is more likely that the presence of purine-rich and pyrimidine-rich stretches in Thermotogae rather reflects their molecular adaptation to high temperature.

Conclusions

PAS, strong SGD and PolC should not be regarded as the signatures of the phylum of Firmicutes, as these features co-exist only in a subset of its members. Moreover, the features may occur, either collectively or individually in members of Fusobacteria, Tenericutes and Thermotogae as well. The study indicates that PAS might warrant the presence of PolC and strong SGD, but the presence of PolC or that of SGD not necessarily implies PAS. In other words, PAS might be a probable, but not an ordained outcome of PolC and strong SGD.

Methods

Sequence retrieval

All predicted protein coding sequences and the complete genome sequences of 102 Firmicutes members were retrieved from the NCBI GenBank. The organisms were chosen in a way to include representatives from all major subphyla and/or classes of the phylum of Firmicutes (Additional file 1: Table S1). Care had also been taken to keep the selection of organisms as varied as possible in terms of their characteristics lifestyle, habitat and genomic G + C-content. However, due to non uniform distribution of organisms of known genome sequences across different families of the Firmicutes, members from some family got overrepresented. Similarly 90 representative organisms of varying G + C-content and niche specificity from all other non-Firmicutes taxa (Additional file 2: Table S2) were also downloaded. All basic information of those organisms were collected from NCBI [40] and BacMap [41] databases.

For each organism under study, presumed duplicates, transposons and the annotated ORFs having less than 300 base pairs have been excluded from the dataset in order to reduce the stochastic errors,

Segregation of two strands of replication (LeS and LaS) and evaluation of SGD in organisms under study

In order to segregate the LeS and LaS genes, one needs to determine the replication origin (oriC) or termination (ter) of the respective genome. It is well known that in bacteria, the base composition of each chromosomal strand changes at the origin and terminus of replication [1315, 4245], which is reflected in the change in sign in the cumulative GC-skew [(G-C)/G + C)] and other skew plots at oriC and ter [4649]. With a view to determine oriC, the cumulative GC-skew analysis was performed with the help of an in-house developed program, using a sliding window of 10 Kb along the entire genome sequence of each species under examination. The oriC predicted from the extrema of the cumulative GC-skew were validated by checking the neighbouring gene organization along with the presence of DnaA boxes in their vicinity [46, 50], and also by comparing the same with the oriC sites of the respective genomes, as annotated in the DoriC database [51]. In most of the cases, the GenBank reference start point of the genome sequence turned out as the putative oriC, though there were a few exceptions.

The putative ter was then calculated as the location of the predicted oriC plus half of the length of the respective chromosome, as done previously by Mao et al. [52]. In majority of the organisms under study, the cumulative GC-skew changed the sign in the neighbourhood of the predicted ter, validating thereby the location of the ter region.

In some exceptional cases, especially in organisms following Trend III or Trend V, the cumulative GC-skew showed zig-zag trajectories with multiple extrema. The chromosomes of these organisms might have undergone large-scale genomic recombination, rearrangements and/or inversions, leading to a mixing of leading and lagging strands of replication and the zig-zag patterns of the cumulative GC-skew might be attributed to such genome rearrangement events. In such cases, the extremum point closest to the point representing the putative oriC plus half of the chromosome length was taken as the putative ter point. It may be argued that the oriC and ter sites in these organisms might undergo a shift from their original positions (i.e., prior to genetic rearrangements) and hence, the predicted oriC plus half of the chromosome length may not always represent the actual ter sites. However, shifting of ter sites would not change the general trends in base usage in such cases. A shift in oriC and/or ter would merely toggle the signs of local GC-skew and AT-skew. Since in Trend III organisms, most of the 10 kb windows have either both the skews positive or both negative and there would be no change in overall trend, if the skews toggle their signs simultaneously. On the other hand, the group of Trend V organisms includes all atypical cases of base combinations with no definite pattern and it is very unlikely that a shift in the oriC/ter sites would change an irregular pattern into a regular or well-defined one. This point has further been elaborated in the Discussion section, along with an example of Yerisina pestis strains, which have reportedly undergone substantial genetic rearrangements.

Based on the predicted oriC and ter sites, the two strands of replication were segregated by joining the oriC to ter region of one half of the plus strand with the ter to oriC region of the minus strand and vice-versa. The numbers of coding regions in two strands of replication were calculated for each genome and the strand with higher frequency of coding regions were taken as the LeS, following the usual convention [3, 52].

In order to ascertain SGD, a 2 × 2 chi-square contingency test was done with number of genes encoded by LeS and LaS, using STATISTICA (version 6.0, published by Statsoft Inc., Tulsa, Oklahoma, USA). Average G + C-content of each genome has also been calculated.

Determination of instantaneous GC-skew, AT-skew and RY-skew for the sequenced genomes used in the study

The total purine-pyrimidine skew values [(R-Y)/(R + Y)] and instantaneous AT-skew values [(A-T)/(A + T)] were also calculated for a sliding window of 10 kb, using an in-house program and subsequent plots have been made.Instantaneous GC-skew (blue color) and AT-skew (red color) values were plotted together against the respective windows along the genome sequence of each organism, in order to find out the distinct trends in purine/pyrimidine distributions. Some representatives of these plots are shown in Figures 1, 2, 3 and 4.The scatter plots of the instantaneous GC-skew and AT-skew values were also drawn in an attempt to affirm the nature of the trends in strand-specific purine and pyrimidine usages in LeS (blue color) and LaS (red color) of each genome, some representatives of which were shown in Figure 5.

Classification of genomes according to the trends in base usage along the respective LeS and LaS sequences

With a view to classify the genomes under study according to the trends in base usage along their two strands of replication, the individual base frequencies were calculated for each sliding window of 10 kb along the LeS sequences. There could be four different combination of base usage in these LeS sequence segments as given below.

  1. (a)

    frequency of G > frequency of C AND frequency of A > frequency of T.

  2. (b)

    frequency of G > frequency of C AND frequency of A ≤ frequency of T.

  3. (c)

    frequency of G ≤ frequency of C AND frequency of A > frequency of T.

  4. (d)

    frequency of G ≤ frequency of C AND frequency of A ≤ frequency of T.

If there had been no strand-specific bias in base usage, the distribution of 10 kb LeS segments among these four possible combinations should have been uniform (around 25%), whatever be their average genomic GC-composition. But all genomes examined in the study showed distinct biases in distribution patterns of LeS segments among four groups. On the basis of observed biases in distribution of 10 kb LeS segments among above four groups, the organisms were classified into five distinct categories, as shown in Tables 1 and 2. The criteria for such classification are given below. Considering up to 5% deviations from the expected frequency of occurrence as normal stochastic variations, ‘random’ refers to frequencies in the normal range, i.e., (25 ± 5%), while ‘high’ and ‘low’ refer to frequencies >30% and <20% respectively.

Trend I: (a) high, (b) random or low, (c) & (d) low → Enrichment of both G and A along LeS.

Trend II: (a) & (b) high, (c) & (d) low → Only G-enrichment along LeS.

Trend III: (a) high, (d) high or random, (b) & (c) low → Presence of both R-dominant & Y-dominant stretches along LeS.

Trend IV: (b) high, (a) random or low, (c) & (d) low → G + T-richness of LeS.

Trend V: all other possible cases such as (a)–(d) all random or (a) high, (b) & (d) random, or (b) high, (c) random etc. → No definite strand-specific bias.

Since these categorization criteria are based on the relative usages of G versus C and A versus T, they hold good for all types of genomes, irrespective of their average G + C-content.

Determination of PolC orthologues in bacteria by BLASTP search

The annotation of PolC in all genomes under study was checked individually from their respective protein tables. There were three possibilities. In most of the PolC-containing species, the genes encoding PolC were unambiguously annotated and hence, could be taken as an evidence of presence of PolC in these organisms. In a few cases, products of some specific genes were marked as “putative DNA polymerase III alpha subunit” or “DNA polymerase III PolC-type”. In these cases, a BLASTP search was carried out with these particular gene sequences against a database of genomes belonging to the genus of the respective organism. Lastly, in cases where no PolC/ PolC-type/DNA Polymerase III alpha subunit gene or gene product could be found, we have taken the annotated PolC sequence(s) from other organisms (from closely related ones, wherever available) and a BLASTP search is carried against the whole genome sequence of the target organism. In both the cases, database hits with e value 0 to 10-e20, if any, were retained and considered as evidences of existence of PolC in the respective organisms.

Determination of base usages at three codon positions and total sequences of individual genes and intergenic regions in leading and lagging strand of replication

Exhaustive base composition analysis was carried out to find out the individual base frequencies in three codon positions of each protein-coding regions of each genome under study, using the program CODONW 1.4.2 (written by John Peden and available at (http://sourceforge.net/projects/codonw/)). The individual purine (G + A) and pyrimidine (C + T) contents and the base frequencies for the total sequence of individual genes (GT, AT, CT, TT) were also calculated. The base usage patterns in intergenic regions (of length ≥ 100 bases) in LeS and LaS sequences of the genomes have also been determined. Since the intergenic regions flanked by the convergently or divergently transcribed genes cannot be unambiguously assigned to any specific strand of replication, only the non-coding sequences existing between two co-oriented genes (i.e., the flanking genes are either both transcribed from the leading strand or both from the lagging strand of replication) have been considered. Each of these base frequencies were then plotted against the respective orders of genes along LeS and LaS of the respective organisms (Figures 7, 8, 9 and 10, Additional file 4: Figure S2 and Additional file 5: Figure S3).Distribution curves of SGD and the histograms of the genomic G + C-contents (Figure 6A and B) were also plotted for different groups of organisms showing distinct trends in purine usage.

Availability of supporting data

The data sets supporting the results of this article are included within the article (and its additional files).

Abbreviations

Les:

Leading strand

LaS:

Lagging strand

PAS:

Purine asymmetry

SGD:

Strand-specific bias in gene distribution

TCR:

Transcription-coupled repair.

References

  1. Hu J, Zhao X, Yu J: Replication-associated purine asymmetry may contribute to strand-biased gene distribution. Genomics. 2007, 90 (2): 186-194. 10.1016/j.ygeno.2007.04.002.

    Article  CAS  PubMed  Google Scholar 

  2. Qu H, Wu H, Zhang T, Zhang Z, Hu S, Yu J: Nucleotide compositional asymmetry between the leading and lagging strands of eubacterial genomes. Res Microbiol. 2010, 161 (10): 838-846. 10.1016/j.resmic.2010.09.015.

    Article  CAS  PubMed  Google Scholar 

  3. Rocha E: Is there a role for replication fork asymmetry in the distribution of genes in bacterial genomes?. Trends Microbiol. 2002, 10 (9): 393-395. 10.1016/S0966-842X(02)02420-4.

    Article  CAS  PubMed  Google Scholar 

  4. Rocha EP, Danchin A: Gene essentiality determines chromosome organisation in bacteria. Nucleic Acids Res. 2003, 31 (22): 6570-6577. 10.1093/nar/gkg859.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  5. Rocha EP, Danchin A: Essentiality, not expressiveness, drives gene-strand bias in bacteria. Nat Genet. 2003, 34 (4): 377-378. 10.1038/ng1209.

    Article  CAS  PubMed  Google Scholar 

  6. Dervyn E, Suski C, Daniel R, Bruand C, Chapuis J, Errington J, Janniere L, Ehrlich SD: Two essential DNA polymerases at the bacterial replication fork. Science. 2001, 294 (5547): 1716-1719. 10.1126/science.1066351.

    Article  CAS  PubMed  Google Scholar 

  7. Bao Q, Tian Y, Li W, Xu Z, Xuan Z, Hu S, Dong W, Yang J, Chen Y, Xue Y, Xu Y, Lai X, Huang L, Dong X, Ma Y, Ling L, Tan H, Chen R, Wang J, Yu J, Yang H: A complete sequence of the T. tengcongensis genome. Genome Res. 2002, 12 (5): 689-700. 10.1101/gr.219302.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  8. Mira A, Pushker R, Legault BA, Moreira D, Rodriguez-Valera F: Evolutionary relationships of Fusobacterium nucleatum based on phylogenetic analysis and comparative genomics. BMC Evol Biol. 2004, 4: 50-10.1186/1471-2148-4-50.

    Article  PubMed Central  PubMed  Google Scholar 

  9. Lobry JR, Sueoka N: Asymmetric directional mutation pressures in bacteria. Genome Biol. 2002, 3 (10): 1-14.

    Article  Google Scholar 

  10. Charneski CA, Honti F, Bryant JM, Hurst LD, Feil EJ: Atypical AT skew in Firmicute genomes results from selection and not from mutation. PLoS Genet. 2011, 7 (9): e1002283-10.1371/journal.pgen.1002283.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  11. Paul S, Million-Weaver S, Chattopadhyay S, Sokurenko E, Merrikh H: Accelerated gene evolution through replication-transcription conflicts. Nature. 2013, 495 (7442): 512-515. 10.1038/nature11989.

    Article  CAS  PubMed  Google Scholar 

  12. Wu H, Qu H, Wan N, Zhang Z, Hu S, Yu J: Strand-biased gene distribution in bacteria is related to both horizontal gene transfer and strand-biased nucleotide composition. Genomics Proteomics Bioinformatics. 2012, 10 (4): 186-196. 10.1016/j.gpb.2012.08.001.

    Article  CAS  PubMed  Google Scholar 

  13. Zhao X, Zhang Z, Yan J, Yu J: GC content variability of eubacteria is governed by the pol III alpha subunit. Biochem Biophys Res Commun. 2007, 356 (1): 20-25. 10.1016/j.bbrc.2007.02.109.

    Article  CAS  PubMed  Google Scholar 

  14. Francino MP, Ochman H: Strand asymmetries in DNA evolution. Trends Genet. 1997, 13 (6): 240-245. 10.1016/S0168-9525(97)01118-9.

    Article  CAS  PubMed  Google Scholar 

  15. Freeman JM, Plasterer TN, Smith TF, Mohr SC: Patterns of genome organization in bacteria. Science. 1998, 279 (5358): 1827-1827. 10.1126/science.279.5358.1827a.

    Article  Google Scholar 

  16. McLean MJ, Wolfe KH, Devine KM: Base composition skews, replication orientation, and gene orientation in 12 prokaryote genomes. J Mol Evol. 1998, 47 (6): 691-696. 10.1007/PL00006428.

    Article  CAS  PubMed  Google Scholar 

  17. McInerney JO: Replicational and transcriptional selection on codon usage in Borrelia burgdorferi. Proc Natl Acad Sci USA. 1998, 95 (18): 10698-10703. 10.1073/pnas.95.18.10698.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  18. Romero H, Zavala A, Musto H: Codon usage in Chlamydia trachomatis is the result of strand-specific mutational biases and a complex pattern of selective forces. Nucleic Acids Res. 2000, 28 (10): 2084-2090. 10.1093/nar/28.10.2084.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  19. Das S, Paul S, Chatterjee S, Dutta C: Codon and amino acid usage in two major human pathogens of genus Bartonella–optimization between replicational-transcriptional selection, translational control and cost minimization. DNA Res. 2005, 12 (2): 91-102. 10.1093/dnares/12.2.91.

    Article  CAS  PubMed  Google Scholar 

  20. Das S, Paul S, Dutta C: Evolutionary constraints on codon and amino acid usage in two strains of human pathogenic actinobacteria Tropheryma whipplei. J Mol Evol. 2006, 62 (5): 645-658. 10.1007/s00239-005-0164-6.

    Article  CAS  PubMed  Google Scholar 

  21. Frank A, Lobry J: Asymmetric substitution patterns: a review of possible underlying mutational or selective mechanisms. Gene. 1999, 238 (1): 65-77. 10.1016/S0378-1119(99)00297-8.

    Article  CAS  PubMed  Google Scholar 

  22. Lao PJ, Forsdyke DR: Thermophilic bacteria strictly obey Szybalski’s transcription direction rule and politely purine-load RNAs with both adenine and guanine. Genome Res. 2000, 10 (2): 228-236. 10.1101/gr.10.2.228.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  23. Trifonov E: Translation framing code and frame-monitoring mechanism as suggested by the analysis of mRNA and 16 S rRNA nucleotide sequences. J Mol Biol. 1987, 194 (4): 643-652. 10.1016/0022-2836(87)90241-5.

    Article  CAS  PubMed  Google Scholar 

  24. Engelen S, Vallenet D, Medigue C, Danchin A: Distinct co-evolution patterns of genes associated to DNA polymerase III DnaE and PolC. BMC Genomics. 2012, 13 (1): 69-10.1186/1471-2164-13-69.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  25. Bohlin J, Hardy S, Ussery D: Stretches of alternating pyrimidine/purines and purines are respectively linked with pathogenicity and growth temperature in prokaryotes. BMC Genomics. 2009, 10 (1): 346-10.1186/1471-2164-10-346.

    Article  PubMed Central  PubMed  Google Scholar 

  26. Rapoport AE, Trifonov EN: Excessive Clustering of Third Codon Position Pyrimidines in Prokaryotes. J Biomol Struct Dyn. 2008, 25 (6): 647-653. 10.1080/07391102.2008.10507211.

    Article  CAS  PubMed  Google Scholar 

  27. Svejstrup JQ: Mechanisms of transcription-coupled DNA repair. Nat Rev Mol Cell Biol. 2002, 3 (1): 21-29. 10.1038/nrm703.

    Article  CAS  PubMed  Google Scholar 

  28. Baran RH, Ko H, Jernigan RW: Methods for comparing sources of strand compositional asymmetry in microbial chromosomes. DNA Res. 2003, 10 (3): 85-95. 10.1093/dnares/10.3.85.

    Article  CAS  PubMed  Google Scholar 

  29. Necşulea A, Lobry JR: A new method for assessing the effect of replication on DNA base composition asymmetry. Mol Biol Evol. 2007, 24 (10): 2169-2179. 10.1093/molbev/msm148.

    Article  PubMed  Google Scholar 

  30. Sohail A, Hayes CS, Divvela P, Setlow P, Bhagwat AS: Protection of DNA by alpha/beta-type small, acid-soluble proteins from Bacillus subtilis spores against cytosine deamination. Biochemistry. 2002, 41 (38): 11325-11330. 10.1021/bi026332t.

    Article  CAS  PubMed  Google Scholar 

  31. Setlow P: I will survive: DNA protection in bacterial spores. Trends Microbiol. 2007, 15 (4): 172-180. 10.1016/j.tim.2007.02.004.

    Article  CAS  PubMed  Google Scholar 

  32. Paredes-Sabja D, Setlow P, Sarker MR: Germination of spores of Bacillales and Clostridiales species: mechanisms and proteins involved. Trends Microbiol. 2011, 19 (2): 85-94. 10.1016/j.tim.2010.10.004.

    Article  CAS  PubMed  Google Scholar 

  33. Nelson KE, Clayton RA, Gill SR, Gwinn ML, Dodson RJ, Haft DH, Hickey EK, Peterson JD, Nelson WC, Ketchum KA: Evidence for lateral gene transfer between Archaea and bacteria from genome sequence of Thermotoga maritima. Nature. 1999, 399 (6734): 323-329. 10.1038/20601.

    Article  CAS  PubMed  Google Scholar 

  34. Ludwig W, Schleifer K-H, Whitman WB: Revised road map to the phylum Firmicutes. Bergey’s Manual® of Systematic Bacteriology. 2009, Springer New York, 1-13.

    Google Scholar 

  35. Wolf M, Muller T, Dandekar T, Pollack JD: Phylogeny of Firmicutes with special reference to Mycoplasma (Mollicutes) as inferred from phosphoglycerate kinase amino acid sequence data. Int J Syst Evol Microbiol. 2004, 54 (3): 871-875. 10.1099/ijs.0.02868-0.

    Article  CAS  PubMed  Google Scholar 

  36. Liang Y, Hou X, Wang Y, Cui Z, Zhang Z, Zhu X, Xia L, Shen X, Cai H, Wang J: Genome rearrangements of completely sequenced strains of Yersinia pestis. J Clin Microbiol. 2010, 48 (5): 1619-1623. 10.1128/JCM.01473-09.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  37. Brown JR, Douady CJ, Italia MJ, Marshall WE, Stanhope MJ: Universal trees based on large combined protein sequence data sets. Nat Genet. 2001, 28 (3): 281-285. 10.1038/90129.

    Article  CAS  PubMed  Google Scholar 

  38. Brochier C, Philippe H: Phylogeny: a non-hyperthermophilic ancestor for bacteria. Nature. 2002, 417 (6886): 244-10.1038/417244a.

    Article  CAS  PubMed  Google Scholar 

  39. Daubin V, Gouy M, Perriere G: A phylogenomic approach to bacterial phylogeny: evidence of a core of genes sharing a common history. Genome Res. 2002, 12 (7): 1080-1090. 10.1101/gr.187002.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  40. Pruitt KD, Tatusova T, Maglott DR: NCBI reference sequences (RefSeq): a curated non-redundant sequence database of genomes, transcripts and proteins. Nucleic Acids Res. 2007, 35 (suppl 1): D61-D65.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  41. Stothard P, Van Domselaar G, Shrivastava S, Guo A, O’Neill B, Cruz J, Ellison M, Wishart DS: BacMap: an interactive picture atlas of annotated bacterial genomes. Nucleic Acids Res. 2005, 33 (suppl 1): D317-D320.

    CAS  PubMed Central  PubMed  Google Scholar 

  42. Arakawa K, Suzuki H, Tomita M: Quantitative analysis of replication-related mutation and selection pressures in bacterial chromosomes and plasmids using generalised GC-skew index. BMC Genomics. 2009, 10 (1): 640-10.1186/1471-2164-10-640.

    Article  PubMed Central  PubMed  Google Scholar 

  43. Grigoriev A: Analyzing genomes with cumulative skew diagrams. Nucleic Acids Res. 1998, 26 (10): 2286-2290. 10.1093/nar/26.10.2286.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  44. Lobry JR: Asymmetric substitution patterns in the two DNA strands of bacteria. Mol Biol Evol. 1996, 13 (5): 660-665. 10.1093/oxfordjournals.molbev.a025626.

    Article  CAS  PubMed  Google Scholar 

  45. Lobry J: A simple vectorial representation of DNA sequences for the detection of replication origins in bacteria. Biochimie. 1996, 78 (5): 323-326. 10.1016/0300-9084(96)84764-X.

    Article  CAS  PubMed  Google Scholar 

  46. Sernova NV, Gelfand MS: Identification of replication origins in prokaryotic genomes. Brief Bioinform. 2008, 9 (5): 376-391. 10.1093/bib/bbn031.

    Article  CAS  PubMed  Google Scholar 

  47. Song J, Ware A, Liu S-L: Wavelet to predict bacterial ori and ter: a tendency towards a physical balance. BMC Genomics. 2003, 4 (1): 17-10.1186/1471-2164-4-17.

    Article  PubMed Central  PubMed  Google Scholar 

  48. Zhang R, Zhang C-T: Identification of replication origins in archaeal genomes based on the Z-curve method. Archaea. 2005, 1 (5): 335-346. 10.1155/2005/509646.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  49. Grigoriev A: Graphical genome comparison: rearrangements and replication origin of Helicobacter pylori. Trends Genet. 2000, 16 (9): 376-378. 10.1016/S0168-9525(00)02087-4.

    Article  CAS  PubMed  Google Scholar 

  50. Mackiewicz P, Zakrzewska-Czerwinska J, Zawilak A, Dudek MR, Cebrat S: Where does bacterial replication start? Rules for predicting the oriC region. Nucleic Acids Res. 2004, 32 (13): 3781-3791. 10.1093/nar/gkh699.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  51. Gao F, Luo H, Zhang C-T: DoriC 5.0: an updated database of oriC regions in both bacterial and archaeal genomes. Nucleic Acids Res. 2013, 41 (D1): D90-D93. 10.1093/nar/gks990.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  52. Mao X, Zhang H, Yin Y, Xu Y: The percentage of bacterial genes on leading versus lagging strands is influenced by multiple balancing forces. Nucleic Acids Res. 2012, 40 (17): 8210-8218. 10.1093/nar/gks605.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

Download references

Acknowledgements

We are grateful to Dr. Sandip Paul, Department of Microbiology, and University of Washington, US for thoughtful and constructive suggestions and critical reading of the manuscript. We wish to thank Ms. Anindya Roy Chowdhury for some technical assistance during the progress of the study. We thankfully acknowledge the infrastructural support obtained from the DBT Bioinformatics Centre of this institute. This work was supported by the Council of Scientific and Industrial Research (CSIR), Govt. of India. (CSIR Network Project GENESIS, BSC0121). SKS and AG are supported by Senior Research Fellowships from CSIR and University Grants Commission (UGC), Govt. of India, respectively.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Chitra Dutta.

Additional information

Competing interests

The authors declare that they have no competing interests.

Author’s contributions

SKS and AG carried out all the computational analyses, prepared the numerical as well as graphical presentations of the data (Additional files and Figures) and made critical reading of the manuscript. CD conceived the project, designed the strategy, coordinated the study and wrote the manuscript. All authors read and approved the final manuscript.

Sanjoy Kumar Saha, Aranyak Goswami contributed equally to this work.

Electronic supplementary material

Additional file 1: Table S1: General features of Firmicutes used in this study.s (PDF 229 KB)

Additional file 2: Table S2: General features of non-Firmicutes used in this study. (PDF 89 KB)

12864_2013_6136_MOESM3_ESM.pdf

Additional file 3: Figure S1: (L) Cumulative GC-skew (blue lines) and AT-skew (red lines) and (R) purine/pyrimidine skews (black lines) in some model representatives of Trend II organisms. (A) Streptococcus agalactiae NEM316, (B) Acidaminococcus intestini RyC-MR95, (C) Geobacillus kaustophilus HTA426, (D) Veillonella parvula DSM 2008, (E) Thermodesulfobium narugense DSM 14796, (F) Clostridiales genomosp BVAB3 UPII9 5, (G) Acinetobacter sp.ADP1, (H) Candidatus Protochlamydia amoebophila UWE25. (PDF 528 KB)

12864_2013_6136_MOESM4_ESM.pdf

Additional file 4: Figure S2: Trends in individual base usages in Escherichia coli str. K-12 substr. MG1655 for genes encoded by both LeS and LaS. Subscripts are same as in Figure 7. (PDF 688 KB)

12864_2013_6136_MOESM5_ESM.pdf

Additional file 5: Figure S3: Trends in individual base usages in Bartonella henselae str.Houston-1 for genes encoded by both leading and lagging strands. Subscripts are same as in Figure 7. (PDF 1 MB)

12864_2013_6136_MOESM6_ESM.pdf

Additional file 6: Figure S4: (L) Instantaneous GC-skew (blue lines) and AT-skew (red lines) and (R) Cumulative GC-skew (blue lines) and AT-skew (red lines) in Yerisina pestis strains. (A) Yersinia pestis CO92, (B) Yersinia pestis D106004, (C) Yersinia pestis D106004 (D) Yersinia pestis Antiqua, (E) Yersinia pestis Nepal516, (F) Yersinia pestis KIM 10, (G) Yersinia pestis biovar Microtus 91001, (H) Yersinia pestis Pestoides F. Table S3. Status of combinations (a) – (d) in Y. pestis strains under study. (PDF 1 MB)

Authors’ original submitted files for images

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://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

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Saha, S.K., Goswami, A. & Dutta, C. Association of purine asymmetry, strand-biased gene distribution and PolC within Firmicutes and beyond: a new appraisal. BMC Genomics 15, 430 (2014). https://doi.org/10.1186/1471-2164-15-430

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1471-2164-15-430

Keywords