Correlations between genomic GC levels and optimal growth temperatures are not 'robust'

Musto et al. [FEBS Lett. 573 (2004) 73] studied the correlations between GC levels and optimal growth temperatures in 20 prokaryotic families. They reported that positive correlations are generally observed, and many of these are significant. Here, we have shown that these correlations are not "robust," i.e., correlation coefficients and/or significance of correlations can be considerably influenced by exclusion of very few (even as small as one) species from each dataset. The sensitivity of correlations is assumed as a result of high levels of bias in the family datasets. We concluded that solely based on these data, one cannot establish that GC contents of prokaryotic genomes increase as a result of growth temperature increments.     

Correlations between genomic GC levels and optimal growth temperatures: some comments

Regarding the existence of any specific correlation between optimal growth temperature and genomic GC levels, Musto et al. [FEBS Lett. 573 (2004) 73] have recently performed analysis on 20 prokaryotic families and showed that in most of the families there exists a positive correlation between these two parameters. On the basis of these results they claimed that optimal growth temperature is one of the factors that influence genomic GC composition in prokaryotes. In a subsequent article, Marashi and Ghalanbor [Biochem. Biophys. Res. Commun. 325 (2004) 381] have demonstrated that the correlation values change substantially when very few points in some of the families were excluded from the data set of Musto et al. [FEBS Lett. 573 (2004) 73]. But Marashi and Ghalanbor have not provided any reason behind this. The points excluded by Marashi and Ghalanbor are actually the outliers in the data set, which strongly affect the correlation coefficients. But the presence of outliers in large data set hardly had any effect on the correlation values. Marashi and Ghalanbor have excluded points from only those families that have small sample sizes and observed a substantial change in correlation coefficient values. Therefore, we argue that any conclusion drawn for a small sample size having outliers is always questionable. Although Musto's approach is a novel one, but to make any generalization one needs to be careful about the flawlessness in the data set.     

Genomic GC level, optimal growth temperature, and genome size in prokaryotes

Two years ago, we showed that positive correlations between optimal growth temperature (T(opt)) and genome GC are observed in 15 out of the 20 families of prokaryotes we analyzed, thus indicating that "T(opt) is one of the factors that influence genomic GC in prokaryotes". Our results were disputed, but these criticisms were demonstrated to be mistaken and based on misconceptions. In a recent report, Wang et al. [H.C. Wang, E. Susko, A.J. Roger, On the correlation between genomic G+C content and optimal growth temperature in prokaryotes: data quality and confounding factors, Biochem. Biophys. Res. Commun. 342 (2006) 681-684] criticize our results by stating that "all previous simple correlation analyses of GC versus temperature have ignored the fact that genomic GC content is influenced by multiple factors including both intrinsic mutational bias and extrinsic environmental factors". This statement, besides being erroneous, is surprising because it applies in fact not to ours but to the authors' article. Here, we rebut the points raised by Wang et al. and review some issues that have been a matter of debate, regarding the influence of environmental factors upon GC content in prokaryotes. Furthermore, we demonstrate that the relationship that exists between genome size and GC level is valid for aerobic, facultative, and microaerophilic species, but not for anaerobic prokaryotes.     

On the correlation between genomic G+C content and optimal growth temperature in prokaryotes: data quality and confounding factors

The correlation between genomic G+C content and optimal growth temperature in prokaryotes has gained renewed interest after Musto et al. [H. Musto, H. Naya, A. Zavala, H. Romero, F. Alvarex-Valin, G. Bernardi, Correlations between genomic GC levels and optimal growth temperatures in prokaryotes, FEBS Lett. 573 (2004) 73-77], reported that positive correlations exist in 15 families studied. We have reanalyzed their data and found that when genome size and data quality were adjusted for, there was no significant evidence of relationship between optimal temperature and GC content for two of the families that had previously shown strongly significant correlations. Using updated temperature optima for Halobacteriaceae species we found the correlation is insignificant in this family. For the family Enterobacteriaceae when genome size and optimal temperature are included in a multiple linear regression, only genome size is significant as a predictor of GC content. We showed that more profound statistical methods than simple two factor correlation analysis should be used for analyzing complex intrinsic and extrinsic factors that affect genomic GC content. We further found that a positive correlation between temperature and genomic GC is only evident in free-living species of low optimal growth temperatures.     

The correlation between genomic G+C and optimal growth temperature of prokaryotes is robust: a reply to Marashi and Ghalanbor

We have recently shown that optimal growth temperature (T(opt)) is one of the factors that influence genomic GC in prokaryotes. Our results have been disputed by Marashi and Ghalanbor, who claim that the correlations we show are not "robust" because the elimination of some points (arbitrarily chosen) leads, in some families, to variations in the correlation coefficients and/or significance of correlations. Here, we test whether the correlation between T(opt) and genomic GC is robust by using two independent approaches: detection of possible outliers (using robust Mahalanobis distance) and usage of a non-parametric correlation coefficient that is not sensitive to the presence of outliers. The results presented here reinforce our previous proposal that T(opt) is correlated with genomic GC in prokaryotes.     

Uracil content of 16S rRNA of thermophilic and psychrophilic prokaryotes correlates inversely with their optimal growth temperatures

We report here the finding of a highly significant inverse correlation of the uracil content of 16S rRNA and the optimum growth temperature (T_opt) of cultured thermophilic and psychrophilic prokaryotes. This correlation was significantly different from the weaker correlations between the contents of other nucleotides and T_opt. Analysis of the 16S rRNA secondary structure regions revealed a fall in the A:U base-pair content in step with the increase in T_opt that was much steeper than that of mismatched base-pairs, which are thermodynamically less stable. These findings indicate that the 16S rRNA sequences of thermophiles and psychrophiles are under a strong thermo-adaptive pressure, and that structure–function constraints play a crucial role in determining their 16S rRNA nucleotide composition. The derived relationship between uracil content and T_opt was used to develop an algorithm to predict the T_opt values of uncultured prokaryotes lacking cultured close relatives and belonging to the phyla predominantly containing thermophiles. This algorithm may be useful in guiding the design of cultivation conditions for hitherto uncultured microbes.     

PGTdb: a database providing growth temperatures of prokaryotes

Included in Prokaryotic Growth Temperature database (PGTdb) are a total of 1334 temperature data from 1072 prokaryotic organisms, Bacteria and Archaea: PGTdb integrates microbial growth temperature data from literature survey with their nucleotide/protein sequence and protein structure data from related databases. A direct correlation is observed between the average growth temperature of an organism and the melting temperature of proteins from the organism. Therefore, this database is useful not only for microbiologists to obtain cultivation condition, but also for biochemists and structure biologists to study the correlation between protein sequences/structures and their thermostability. In addition, the taxonomy and ribosomal RNA sequence(s) of an organism are linked through NCBI Taxonomy and the Ribosomal RNA Operon Copy Number Database umdb, respectively. PGTdb is the only integrated database on the Internet to provide the growth temperature data of the prokaryotes and the combined information of their nucleotide/protein sequences, protein structures, taxonomy and phylogeny. AVAILABILITY: http://pgtdb.csie.ncu.edu.tw     

r/K selection of GC content in prokaryotes

The guanine/cytosine (GC) content of prokaryotic genomes is species-specific, taking values from 16% to 77%. This diversity of selection for GC content remains contentious. We analyse the correlations between GC content and a range of phenotypic and genotypic data in thousands of prokaryotes. GC content integrates well with these traits into r/K selection theory when phenotypic plasticity is considered. High GC-content prokaryotes are r-strategists with cheaper descendants thanks to a lower average amino acid metabolic cost, colonize unstable environments thanks to flagella and a bacillus form and are generalists in terms of resource opportunism and their defence mechanisms. Low GC content prokaryotes are K-strategists specialized for stable environments that maintain homeostasis via a high-cost outer cell membrane and endospore formation as a response to nutrient deprivation, and attain a higher nutrient-to-biomass yield. The lower proteome cost of high GC content prokaryotes is driven by the association between GC-rich codons and cheaper amino acids in the genetic code, while the correlation between GC content and genome size may be partly due to functional diversity driven by r/K selection. In all, molecular diversity in the GC content of prokaryotes may be a consequence of ecological r/K selection.     

Induction of gene expression in bacteria at optimal growth temperatures

Traditional temperature-sensitive systems use either heat shock (40–42 °C) or cold shock (15–23 °C) to induce gene expression at temperatures that are not the optimal temperature for host cell growth (37 °C). This impacts the overall productivity and yield by disturbing cell growth and cellular metabolism. Here, we have developed a new system which controls gene expression in Escherichia coli at more permissive temperatures. The temperature-sensitive cI857-P _L system and the classic lacI-P _lacO system were connected in series to control the gene of interest. When the culture temperature was lowered, the thermolabile cI857 repressor was activated and blocked the expression of lacI from P _L. Subsequently, the decrease of LacI derepressed the expression of gene of interest from P _lacO . Using a green fluorescent protein marker, we demonstrated that (1) gene expression was tightly regulated at 42 °C and strongly induced by lowering temperature to 25–37 °C; (2) different levels of gene expression can be induced by varying culture temperature; and (3) gene expression after induction was sustained until the end of the log phase. We then applied this system in the biosynthesis of acetoin and demonstrated that high yield and production could be achieved using temperature induction. The ability to express proteins at optimal growth temperatures without chemical inducers is advantageous for large-scale and industrial fermentations.     

IslandPath: aiding detection of genomic islands in prokaryotes

Genomic islands (clusters of genes of potential horizontal origin in a prokaryotic genome) are frequently associated with a particular adaptation of a microbe that is of medical, agricultural or environmental importance, such as antibiotic resistance, pathogen virulence, or metal resistance. While many sequence features associated with such islands have been adopted separately in applications for analysis of genomic islands, including pathogenicity islands, there is no single application that integrates multiple features for island detection. IslandPath is a network service which incorporates multiple DNA signals and genome annotation features into a graphical display of a bacterial or archaeal genome, to aid the detection of genomic islands. AVAILABILITY: This application is available at http://www.pathogenomics.sfu.ca/islandpath and the source code is freely available, under GNU public licence, from the authors. SUPPLEMENTARY INFORMATION: An online help file, which includes analyses of the utility of IslandPath, can be found at http://www.pathogenomics.sfu.ca/islandpath/current/islandhelp.html     

Correlations between heparan sulfate metabolism and hepatoma growth

A rat hepatoma cell line (Gershenson et al., Science, 170:859-861, 1970) contains a dynamic steady-state pool of free heparan sulfate (HS) chains in the nucleus that increases in amount when growing cells reach confluence (Fedarko and Conrad, J. Cell Biol., 102:587-599, 1986). In logarithmically growing cells labeled with 35SO4(2-) steady-state levels of [35SO4]HS in the nucleus are altered by a variety of culture conditions. Rapidly dividing cells (doubling time = 18-22 h) growing under optimized conditions had steady-state levels of nuclear HS within the range of 40-50 pmol 35SO4 in nuclear HS/10(6) cells. The steady-state levels of nuclear HS were lowered by several changes in culture conditions, including 1) additions of 1 mM p-nitrophenyl-beta-D-xyloside, 0.25-0.5 mM (+)-catechin, 0.5 ng/ml transforming growth factor beta, 20 ng/ml phorbol-12-myristate-13-acetate, 1 mM dibutyryl cAMP, or 1 mM inositol-2-PO4; 2) decreased levels of D-glucose; or 3) deletions of serum, insulin, or inositol. In all cases lowering of the nuclear HS level was accompanied by an increase in the cell doubling times, suggesting a correlation in which nuclear HS levels must be optimized for maximal growth rates. When cells cultured under optimal growth conditions reached confluence, the level of nuclear HS increased threefold and the cells stopped dividing. The same culture conditions that lowered the steady-state levels of HS in the logarithmically growing cells prevented this rise in the nuclear HS as the cells reached confluence and resulted in loss of contact inhibition and overgrowth of the confluent cultures. These observations suggest a second correlation in which elevated nuclear HS levels are found when cell growth is inhibited at confluence; prevention of this rise results in continued growth. Consistent with this correlation between elevated nuclear HS and reduced growth rates, it was observed that addition of either 0.5 microgram/ml hydrocortisone or 0.05 microgram/ml retinoic acid to the culture medium of logarithmically growing cultures resulted in increases in steady-state levels of nuclear HS that were accompanied by increased cell doubling times. The two agents that increased the levels of nuclear HS in logarithmically growing cultures had little effect on levels of nuclear HS in confluent cells or on contact inhibition.(ABSTRACT TRUNCATED AT 400 WORDS)     

Detection and characterization of horizontal transfers in prokaryotes using genomic signature

Horizontal DNA transfer is an important factor of evolution and participates in biological diversity. Unfortunately, the location and length of horizontal transfers (HTs) are known for very few species. The usage of short oligonucleotides in a sequence (the so-called genomic signature) has been shown to be species-specific even in DNA fragments as short as 1 kb. The genomic signature is therefore proposed as a tool to detect HTs. Since DNA transfers originate from species with a signature different from those of the recipient species, the analysis of local variations of signature along recipient genome may allow for detecting exogenous DNA. The strategy consists in (i) scanning the genome with a sliding window, and calculating the corresponding local signature (ii) evaluating its deviation from the signature of the whole genome and (iii) looking for similar signatures in a database of genomic signatures. A total of 22 prokaryote genomes are analyzed in this way. It has been observed that atypical regions make up ~6% of each genome on the average. Most of the claimed HTs as well as new ones are detected. The origin of putative DNA transfers is looked for among ~12000 species. Donor species are proposed and sometimes strongly suggested, considering similarity of signatures. Among the species studied, Bacillus subtilis, Haemophilus Influenzae and Escherichia coli are investigated by many authors and give the opportunity to perform a thorough comparison of most of the bioinformatics methods used to detect HTs.     

Protein kinases and protein phosphatases in prokaryotes: a genomic perspective

For many years, the regulation of protein structure and function by phosphorylation and dephosphorylation was considered a relatively recent invention that arose independently in each phylogenetic domain. Over time, however, incidents of apparent domain trespass involving the presence of 'eukaryotic' protein kinases or protein phosphatases in prokaryotic organisms were reported with increasing frequency. Today, genomics has provided the means to examine the phylogenetic distribution of 'eukaryotic' protein kinases and protein phosphatases in a comprehensive and systematic manner. The results of these genome searches challenge previous conceptions concerning the origins and evolution of this versatile regulatory mechanism.     

Correlations between polyamine ratios and growth patterns in seedling roots

The levels of putrescine, cadaverine, spermidine and spermine were determined in seedling roots of pea, tomato, millet and corn, as well as in corn coleoptiles and pea internodes. In all roots, putrescine content increased as elongation progressed, and the putrescine/spermine ratio closely paralleled the sigmoid growth curve up until the time of lateral root initiation. Spermidine and spermine were most abundant near the apices and declined progressively with increasing age of the cells. In the zone of differentiation of root hairs in pea roots, putrescine rose progressively with increasing age, while cadaverine declined. In both pea internodes and corn coleoptiles, the putrescine/spermidine ratio rises with increasing age and elongation. Thus, a block in the conversion of the diamine putrescine to the triamine spermidine may be an important step in the change from cell division to cell elongation.Supported by a Fellowship of the Peoples' Republic of China.Aided by grant 5-RO1-AGO2742 from the National Institutes of Health to A.W.G.     

Genomic GC content prediction in prokaryotes from a sample of genes

GC level is a key feature in prokaryotic genomes. Widely employed in evolutionary studies, new insights appear however limited because of the relatively low number of characterized genomes. Since public databases mainly comprise several hundreds of prokaryotes with a low number of sequences per genome, a reliable prediction method based on available sequences may be useful for studies that need a trustworthy estimation of whole genomic GC. As the analysis of completely sequenced genomes shows a great variability in distributional shapes, it is of interest to compare different estimators. Our analysis shows that the mean of GC values of a random sample of genes is a reasonable estimator, based on simplicity of the calculation and overall performance. However, usually sequences come from a process that cannot be considered as random sampling. When we analyzed two introduced sources of bias (gene length and protein functional categories) we were able to detect an additional bias in the estimation for some cases, although the precision was not affected. We conclude that the mean genic GC level of a sample of 10 genes is a reliable estimator of genomic GC content, showing comparable accuracy with many widely employed experimental methods.     

Large-scale comparative genomic analyses of cytoplasmic membrane transport systems in prokaryotes

The recent advancements in genome sequencing make it possible for the comparative analyses of essential cellular processes like transport in organisms across the three domains of life. Membrane transporters play crucial roles in fundamental cellular processes and functions in prokaryotic systems. Between 3 and 16% of open reading frames in prokaryotic genomes were predicted to encode membrane transport proteins, emphasizing the importance of transporters in their lifestyles. Hierarchical clustering of phylogenetic profiles of transporter families, which are derived from the presence or absence of a certain transporter family, showed distinct clustering patterns for obligate intracellular organisms, plant/soil-associated microbes and autotrophs. Obligate intracellular organisms possess the fewest types and number of transporters presumably due to their relatively stable living environment, while plant/soil-associated organisms generally encode the largest variety and number of transporters. A group of autotrophs are clustered together largely due to their absence of transporters for carbohydrate and organic nutrients and the presence of transporters for inorganic nutrients. Inside of each group, organisms are further clustered by their phylogenetic properties. These findings strongly suggest the correlation of transporter profiles to both evolutionary history and the overall physiology and lifestyles of the organisms.     

Evolution of the genomic systems of prokaryotes and its momentous consequences

The earliest self-reproducing cell on Earth, our common ancestor, was probably as small as present-day bacteria. It gave rise to a very large and durable clone whose descendants must have been the only living occupants of the oceans for about one thousand million years. They reached astronomical numbers of separate, disjunct cells, and synthesized many new genes. Their small volume could not accommodate ever larger genomes and useful new genes replaced resident, less successful sequences, thus increasing diversity and the number of strains with highly specialized, distinct, bioenergetic potentialities. Also, selective pressure favored strains able to participate successfully in division of labor and in the sharing of diverse abilities in mixed communities, counterbalancing the limited capacities of individual genomes. Lateral gene transfer mechanisms appeared and were progressively improved, furthering the development of diversity. The prokaryotes' constructive evolution resulted in the formation of a worldwide web of genetic information, and a global bacterial superbiosystem (superorganism). By contrast, eukaryotic evolution of organisms has been typically Darwinian. Diversification of eukaryotic organisms was, however, considerably enriched and accelerated by symbioses with prokaryotes. The more broadly diversified bioenergetic potential of prokaryotes considerably increased the diversity of eukaryotes. Without their participation, our biosphere would have remained much less diverse and less dynamic. Environmental homeostasis has been maintained all along by guided bacterial evolution.     

Structure-dependent relationships between growth temperature of prokaryotes and the amino acid frequency in their proteins

We studied the amino acid frequency and substitution patterns between homologues of prokaryotic species adapted to temperatures in the range 0–102°C, and found a significant temperature-dependent difference in frequency for many of the amino acids. This was particularly clear when we analysed the surface and core residues separately. The difference between the surface and the core is getting more pronounced in proteins adapted to warmer environments, with a more hydrophobic core, and more charged and long-chained amino acids on the surface of the proteins. We also see that mesophiles have a more similar amino acid composition to psychrophiles than to thermophiles, and that archea appears to have a slightly different pattern of substitutions than bacteria. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Links between DNA replication and recombination in prokaryotes

Recombination plays a crucial role in underpinning genome duplication, ensuring that replication blocks are removed or bypassed, and that the replication machinery is subsequently reloaded back onto the DNA. Recent studies have identified a surprising variety of ways in which damaged replication forks are repaired and have shown that the mechanism used depends on the nature of the original blocking lesion. Indeed, an emerging theme is that a single recombination enzyme or complex can perform highly varied tasks, depending on the context of the recombination reaction.