Analysis of intra-genomic GC content homogeneity within prokaryotes
© Bohlin et al. 2010
Received: 8 January 2010
Accepted: 6 August 2010
Published: 6 August 2010
Skip to main content
© Bohlin et al. 2010
Received: 8 January 2010
Accepted: 6 August 2010
Published: 6 August 2010
Bacterial genomes possess varying GC content (total guanines (Gs) and cytosines (Cs) per total of the four bases within the genome) but within a given genome, GC content can vary locally along the chromosome, with some regions significantly more or less GC rich than on average. We have examined how the GC content varies within microbial genomes to assess whether this property can be associated with certain biological functions related to the organism's environment and phylogeny. We utilize a new quantity GCVAR, the intra-genomic GC content variability with respect to the average GC content of the total genome. A low GCVAR indicates intra-genomic GC homogeneity and high GCVAR heterogeneity.
The regression analyses indicated that GCVAR was significantly associated with domain (i.e. archaea or bacteria), phylum, and oxygen requirement. GCVAR was significantly higher among anaerobes than both aerobic and facultative microbes. Although an association has previously been found between mean genomic GC content and oxygen requirement, our analysis suggests that no such association exits when phylogenetic bias is accounted for. A significant association between GCVAR and mean GC content was also found but appears to be non-linear and varies greatly among phyla.
Our findings show that GCVAR is linked with oxygen requirement, while mean genomic GC content is not. We therefore suggest that GCVAR should be used as a complement to mean GC content.
The knowledge of the chemical basis for nucleic acids goes back more than a hundred years, to the work of Miescher . By the early 1950's, it was known that the relative frequency of the four DNA bases ("base composition") was different for different organisms , and in general the number of A's was equal to the number of T's, and the number of G's was the same as the number of C's; this is known as 'Chargaff's first parity rule' . Further, for nearly all genomes studied, the parity rule appears to extend to each strand of the chromosome, when averaged over long distances , although in bacterial chromosomes, there is a clear bias of G's towards the replication leading strand, and for some genomes (many Firmicutes, for example) the A's are also biased towards the leading strand . For a circular chromosome with the replication origin and terminus on exactly opposite sides, this bias of G's towards the replication leading strand will average out to near zero, when one only looks at the DNA sequence in the GenBank file, and the sequence will appear to conform to Chargaff's second rule.
From the Genbank database at NCBI http://www.ncbi.nlm.nih.gov/genomes/lproks.cgi it can be seen that GC content in prokaryotes ranges from 16.6% in Carsonella ruddii strain Pv to 74.9% in Anaeromyxobacter dehalogenans Strain 2CP-C. Within a given genome, the GC content along the chromosome can vary, although since most bacterial genomes have a high coding density, usually the variation is less than that found in eukaryotes . The average genomic GC content is an important property in microbial genomes and has been associated with properties such as genome size , oxygen, and nitrogen exposure [8, 9] and specific habitats [10–13]. For instance, intracellular bacteria have, on average, smaller genomes and are mostly AT rich, while soil bacteria tend to have larger genomes and higher %GC . Higher AT content in intracellular bacteria may be attributed to a loss of repair genes; this loss will eventually lead to an increase in mutation rates from cytosine to thymine [15, 16]. Genes not expressed will eventually lead to reduced genome sizes [15, 16]. Higher GC content in soil bacteria may be due to the increased availability of nitrogen . However, increased nitrogen in the soil does not explain why GC rich bacteria often have larger genomes. The base composition in GC rich genomes might reflect stronger selective forces than AT rich genomes [17–19]. This may indicate that GC rich microbes live in more complex environments than intracellular bacteria . The reasons for stronger selective forces and GC richness is not known, but may be connected to the fact that considerably more energy is required to de-stack GC rich DNA sequences than AT rich DNA sequences .
Although GC content has been found to vary only slightly within prokaryotic genomes some regions differ more than others. A large region flanking the replication origin, for instance, is more GC rich than the average genomic GC content  whereas the region around replication terminus is more AT rich . Surface proteins and RNA genes often have GC content that differs from the average genomic GC content , and protein coding regions have been found to be, on average, approximately 5% more GC rich than non-coding regions . In addition to being more GC rich, coding regions have been found to be more homogeneous in terms of base composition than non-coding regions . The GC heterogeneity in coding regions has, however, been found to be associated with mean genomic AT content in non-coding regions [18, 23]. In other words, GC content variability tends to increase with higher mean genomic AT content in non-coding regions.
Horizontally transferred DNA may have a different fraction of GC than the host genome as a result of different evolutionary pressures [6, 24–26]. Since horizontally transferred DNA is often linked to pathogenesis in microbes , detection of such regions is of great importance. The GC content of foreign DNA will, however, become progressively more similar to the host genome in a process known as amelioration  making such regions more difficult to detect as time progress . The conformation of base compositional patterns from foreign DNA to host DNA may be related to the finding that a particular subunit of the DNA polymerase III, the Pol III α subunit, appears to be driving genomic GC content in prokaryotes .
There is a considerable amount of research and documentation related to mean genomic GC content in prokaryotes demonstrating that this property is the result of many factors interacting in a highly complex manner . On the other hand, analysis of genomic GC content variability within microbial chromosomes, has received much less attention. A more recent overview of methodology used to analyze GC content variation within genomes can be found in Bernaola-Galvan et. al., , and a study of how intra-genomic GC content variation affects codon usage is described by Daubin et. al. . In the present work, we introduce the GCVAR measure to examine GC content variability within prokaryotic genomes. The GCVAR metric gives a measure of how GC content varies within a given genome with respect to the mean genomic GC content. A low GCVAR thus points to little GC content variation, or GC content homogeneity, within the genome, while a high GCVAR designates varying GC content, or GC content heterogeneity.
To the best of our knowledge, no study has examined the interplay between environmental factors and GC content homogeneity in prokaryotes. In the present study the aim was therefore to examine whether GC content homogeneity in prokaryotes, measured here using the GCVAR measure, could be related to specific factors in the environment such as temperature and oxygen, as well as the broader properties implicated in phylogeny and GC content. To do this, regression analyses were performed using GCVAR as the response variable. The response variable was fitted to the following variables: oxygen requirement (a categorical variable defined as either aerobic, anaerobic or facultative), phylum, genomic GC content, genome size, growth temperature (a categorical variable used to define psychrophiles, mesophiles and thermophiles), pathogenicity (a dichotomous variable describing whether the microbe is pathogenic or not) and habitat (a categorical variable describing the environment where the microbe is found, i.e. aquatic, host-associated, multiple, specialized and terrestrial). The dataset consisted of 488 genomes (526 chromosomes) with similar strains and species removed from the analysis to reduce phylogenetic bias.
The coefficient estimates from the GCVAR regression model
Average size (mbp)
Table 1 shows that GCVAR is significantly influenced by phylum. We find that 10 phyla have GCVAR significantly above the average phylum, and 4 phyla have GCVAR significantly below average. The two phyla, Crenarchaeota and Euryarcheota, (both archaea) are among the groups with an above average GCVAR. The archaea domain, consisting predominantly of organisms living in extreme environments, had a significantly higher GCVAR than bacteria (p < 0.001). The highest GCVAR are found in the aquatic group Cyanobacteria, which is largely populated with species capable of photosynthesis . The lowest GCVAR are found in the phylum of the aquatic Planctomycetes, but this group is only based on one single genome, therefore no conclusions can be assumed at the phylum level.
To examine how GCVAR was affected by phylogenetic bias a regression model similar to the one described above was fitted, i.e. GCVAR was the response variable, with mean genomic GC content, oxygen requirement, habitat, optimal growth temperature, and genome size as predictors. In addition, an interaction term between GC content and phylum was added to account for more similar GC content within phyla (Equation (5) in the Methods section). Using this regression model we found that oxygen requirement was the only significant factor (p < 0.001). GCVAR was significantly higher in the genomes of anaerobic microbes (103 chromosomes) as compared to the genomes of aerobic microbes, meaning that the genomes of anaerobic microbes tend to have a more heterogeneous distribution of GC content than genomes of aerobic microbes (246 chromosomes). Facultative microbes were found to have GCVAR values in the region between aerobic and anaerobic microbes (see Figure 2).
The regression models described above indicates that aerobic microbes have genomes with more homogeneous GC content than those of organisms with facultative and anaerobic oxygen requirement (see Figure 2). It has been shown that GC rich genomes tend to be more homogeneous in terms of base composition than AT rich genomes [18, 19, 34]. Aerobic microbes have been associated with GC rich genomes . This result is supported by our linear regression model only when we ignore phylogenetic bias is (p < 0.001). However, adding phyla as a predictor (Equation (5) in the Methods section) fails to demonstrate such an association (p ~0.9).
There are many indications that mean genomic GC content is as much affected by the environment as by phyla [10–13]. It is also well known that chromosomally integrated foreign DNA may differ in base composition as compared to host DNA. The difference in base composition between foreign and host DNA is assumed to be the result of exposure to different selective pressures. It is thought that such genetic regions may be acquired from horizontal transfer or other means of DNA uptake [6, 24–26]. Since pathogenesis is often associated with horizontally transferred DNA, i.e. pathogenicity islands, , establishing a link between any genomic property and horizontal transfer is of considerably interest. However, no significant association (p ~ 0.25) was found between the dichotomous pathogenicity factor and GCVAR using the regression model that included all covariates discussed above (Equation (3) in the Methods section).
The introduction of atmospheric oxygen is presumed to have had profound effects on environment and life . Increase in atmospheric oxygen is believed to have influenced cellular compartmentalization and thus to have been instrumental in the evolution of eukaryotes . Prokaryotes were also affected by the introduction of oxygen  in that while some remained anaerobic others adapted to an aerobic metabolism .
The precise effect of increase in atmospheric oxygen on prokaryotic genomes is debated [37, 38]. A negative correlation has been found between proteomic oxygen content and genomic GC content . Although it has been suggested that genomic GC content is also affected by an aerobic lifestyle , the effects on prokaryote genome composition has remained unclear [37, 38]. Indeed, our own results presented above do not support any connection between genomic GC content and aerobiosis. Our results did, however, find a significant association between GCVAR and oxygen requirement. This greater GC content homogeneity found in aerobes implies that the genomes of these organisms have been subjected to stronger selective pressures than the genomes of anaerobes. This is supported by the recent report that metabolic networks of aerobic bacteria are more complex than those of anaerobic bacteria . From Figure 2 it can be seen that GCVAR appears to be progressively decreasing in facultative and aerobic prokaryotes, respectively.
In summary, we found that GCVAR was associated with oxygen requirement. It is possible that GCVAR is associated with GC content, but from Figure 3 it appears to be a highly non-linear relationship. Other factors such as genome size, habitat and growth temperature were not found significant in the GCVAR model. GCVAR was however found to be higher in archaea than bacteria. By adding an interaction term to model the closer similarity between the genomes in the same phylogenetic group, we found that oxygen requirement was not significantly associated with mean genomic GC content in microbes.
The different results obtained for the models describing GCVAR and mean genomic GC content imply that these properties are governed by different influences, or are interrelated in a non-linear manner. Thus, our findings suggest that GCVAR is linked with oxygen requirement, while mean genomic GC content is not.
All genomes and related information were gathered from the NCBI web site http://www.ncbi.nlm.nih.gov/genomes/lproks.cgi. The statistical package R was used for statistical analyses and graphical representations.
N is the maximum number of non-overlapping 100 bp sliding windows that can fit into the chromosome that is being analyzed. The log-transformation makes GCVAR s empirical distribution more Gaussian-like, for convenience in subsequent linear regression model fitting and statistical inference. Since the optimal sized sliding window varies from genome to genome , different window lengths were tested. The sliding window width of 100 bp was chosen to make the test as sensitive as possible. The other sliding window lengths tested contained 500, 1000, and 2000 bp. The 100 bp sliding window was found to be large enough to carry genome specific information without discarding weak genomic signals as noise. Since the aim of this study was to examine GC content difference within genomes, non-overlapping sliding windows were used to avoid bias and interactions from neighboring genetic regions.
where μ is the overall intercept, α y , y = 1,...,22, are the effects of phylum , δ o , o = 1,2,3, are the effects of oxygen requirement , κ t , t = 1,2,3, are the effects of growth temperature , λ a , a = 1,2, are the effects of pathogenicity and η h, h = 1,...4, are the effects of habitat. β and γ are the regression coefficients for the continuous variables GC and Size, respectively.
The α y in this model are defined as regression coefficients for each of the 22 phylum categories.
μ, α y , δ o are the same effects as those described for Equation (3).
The authors thank Colleen Ussery for critically reading the manuscript and the reviewers for their insightful suggestions.
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.