A hot story from comparative genomics: reverse gyrase is the only hyperthermophile-specific protein

Growth temperature and genome size in bacteria are negatively correlated, suggesting genomic streamlining during thermal adaptation

Prokaryotic genomes are small and compact. Either this feature is caused by neutral evolution or by natural selection favoring small genomes—genome streamlining. Three separate prior lines of evidence argue against streamlining for most prokaryotes. We find that the same three lines of evidence argue for streamlining in the genomes of thermophile bacteria. Specifically, with increasing habitat temperature and decreasing genome size, the proportion of genomic DNA in intergenic regions decreases. Furthermore, with increasing habitat temperature, generation time decreases. Genome-wide selective constraints do not decrease as in the reduced genomes of host-associated species. Reduced habitat variability is not a likely explanation for the smaller genomes of thermophiles. Genome size may be an indirect target of selection due to its association with cell volume. We use metabolic modeling to demonstrate that known changes in cell structure and physiology at high temperature can provide a selective advantage to reduce cell volume at high temperatures.

Strain-level genomic variation in natural populations of Lebetimonas from an erupting deep-sea volcano

Chemolithoautotrophic Epsilonproteobacteria are ubiquitous in sulfidic, oxygen-poor habitats, including hydrothermal vents, marine oxygen minimum zones, marine sediments and sulfidic caves and have a significant role in cycling carbon, hydrogen, nitrogen and sulfur in these environments. The isolation of diverse strains of Epsilonproteobacteria and the sequencing of their genomes have revealed that this group has the metabolic potential to occupy a wide range of niches, particularly at dynamic deep-sea hydrothermal vents. We expand on this body of work by examining the population genomics of six strains of Lebetimonas, a vent-endemic, thermophilic, hydrogen-oxidizing Epsilonproteobacterium, from a single seamount in the Mariana Arc. Using Lebetimonas as a model for anaerobic, moderately thermophilic organisms in the warm, anoxic subseafloor environment, we show that genomic content is highly conserved and that recombination is limited between closely related strains. The Lebetimonas genomes are shaped by mobile genetic elements and gene loss as well as the acquisition of novel functional genes by horizontal gene transfer, which provide the potential for adaptation and microbial speciation in the deep sea. In addition, these Lebetimonas genomes contain two operons of nitrogenase genes with different evolutionary origins. Lebetimonas expressed nifH during growth with nitrogen gas as the sole nitrogen source, thus providing the first evidence of nitrogen fixation in any Epsilonproteobacteria from deep-sea hydrothermal vents. In this study, we provide a comparative overview of the genomic potential within the Nautiliaceae as well as among more distantly related hydrothermal vent Epsilonproteobacteria to broaden our understanding of microbial adaptation and diversity in the deep sea.

The common ancestor of archaea and eukarya was not an archaeon

It is often assumed that eukarya originated from archaea. This view has been recently supported by phylogenetic analyses in which eukarya are nested within archaea. Here, I argue that these analyses are not reliable, and I critically discuss archaeal ancestor scenarios, as well as fusion scenarios for the origin of eukaryotes. Based on recognized evolutionary trends toward reduction in archaea and toward complexity in eukarya, I suggest that their last common ancestor was more complex than modern archaea but simpler than modern eukaryotes (the bug in-between scenario). I propose that the ancestors of archaea (and bacteria) escaped protoeukaryotic predators by invading high temperature biotopes, triggering their reductive evolution toward the "prokaryotic" phenotype (the thermoreduction hypothesis). Intriguingly, whereas archaea and eukarya share many basic features at the molecular level, the archaeal mobilome resembles more the bacterial than the eukaryotic one. I suggest that selection of different parts of the ancestral virosphere at the onset of the three domains played a critical role in shaping their respective biology. Eukarya probably evolved toward complexity with the help of retroviruses and large DNA viruses, whereas similar selection pressure (thermoreduction) could explain why the archaeal and bacterial mobilomes somehow resemble each other.

The legacy of Carl Woese and Wolfram Zillig: from phylogeny to landmark discoveries

Two pioneers of twentieth century biology passed away during the past decade, Wolfram Zillig in April 2005 and Carl Woese in December 2012. Among several other accomplishments, Woese has been celebrated for the discovery of the domain Archaea and for establishing rRNA as the 'Rosetta Stone' of evolutionary and environmental microbiology. His work inspired many scientists in various fields of biology, and among them was Wolfram Zillig, who is credited with the discovery of several unique molecular features of archaea. In this Essay, we highlight the remarkable achievements of Woese and Zillig and consider how they have shaped the archaeal research landscape.

Genetic manipulation in Sulfolobus islandicus and functional analysis of DNA repair genes

Recently, a novel gene-deletion method was developed for the crenarchaeal model Sulfolobus islandicus, which is a suitable tool for addressing gene essentiality in depth. Using this technique, we have investigated functions of putative DNA repair genes by constructing deletion mutants and studying their phenotype. We found that this archaeon may not encode a eukarya-type of NER (nucleotide excision repair) pathway because depleting each of the eukaryal NER homologues XPD, XPB and XPF did not impair the DNA repair capacity in their mutants. However, among seven homologous recombination proteins, including RadA, Hel308/Hjm, Rad50, Mre11, HerA, NurA and Hjc, only the Hjc nuclease is dispensable for cell viability. Sulfolobus encodes redundant BER (base excision repair) enzymes such as two uracil DNA glycosylases and two putative apurinic/apyrimidinic lyases, but inactivation of one of the redundant enzymes already impaired cell growth, highlighting their important roles in archaeal DNA repair. Systematically characterizing these mutants and generating mutants lacking two or more DNA repair genes will yield further insights into the genetic mechanisms of DNA repair in this model organism.

How hyperthermophiles adapt to change their lives: DNA exchange in extreme conditions

Transfer of DNA has been shown to be involved in genome evolution. In particular with respect to the adaptation of bacterial species to high temperatures, DNA transfer between the domains of bacteria and archaea seems to have played a major role. In addition, DNA exchange between similar species likely plays a role in repair of DNA via homologous recombination, a process that is crucial under DNA damaging conditions such as high temperatures. Several mechanisms for the transfer of DNA have been described in prokaryotes, emphasizing its general importance. However, until recently, not much was known about this process in prokaryotes growing in highly thermophilic environments. This review describes the different mechanisms of DNA transfer in hyperthermophiles, and how this may contribute to the survival and adaptation of hyperthermophilic archaea and bacteria to extreme environments.

New mechanistic and functional insights into DNA topoisomerases

DNA topoisomerases are nature's tools for resolving the unique problems of DNA entanglement that occur owing to unwinding and rewinding of the DNA helix during replication, transcription, recombination, repair, and chromatin remodeling. These enzymes perform topological transformations by providing a transient DNA break, formed by a covalent adduct with the enzyme, through which strand passage can occur. The active site tyrosine is responsible for initiating two transesterifications to cleave and then religate the DNA backbone. The cleavage reaction intermediate is exploited by cytotoxic agents, which have important applications as antibiotics and anticancer drugs. The reactions mediated by these enzymes can also be regulated by their binding partners; one example is a DNA helicase capable of modulating the directionality of strand passage, enabling important functions like reannealing denatured DNA and resolving recombination intermediates. In this review, we cover recent advances in mechanistic insights into topoisomerases and their various cellular functions.

The hyperthermophilic bacterium Aquifex aeolicus: from respiratory pathways to extremely resistant enzymes and biotechnological applications

Aquifex aeolicus isolated from a shallow submarine hydrothermal system belongs to the order Aquificales which constitute an important component of the microbial communities at elevated temperatures. This hyperthermophilic chemolithoautotrophic bacterium, which utilizes molecular hydrogen, molecular oxygen, and inorganic sulfur compounds to flourish, uses the reductive TCA cycle for CO(2) fixation. In this review, the intricate energy metabolism of A. aeolicus is described. As the chemistry of sulfur is complex and multiple sulfur species can be generated, A. aeolicus possesses a multitude of different enzymes related to the energy sulfur metabolism. It contains also membrane-embedded [NiFe] hydrogenases as well as oxidases enzymes involved in hydrogen and oxygen utilization. We have focused on some of these proteins that have been extensively studied and characterized as super-resistant enzymes with outstanding properties. We discuss the potential use of hydrogenases in an attractive H(2)/O(2) biofuel cell in replacement of chemical catalysts. Using complete genomic sequence and biochemical data, we present here a global view of the energy-generating mechanisms of A. aeolicus including sulfur compounds reduction and oxidation pathways as well as hydrogen and oxygen utilization.

Crystal structures of Thermotoga maritima reverse gyrase: inferences for the mechanism of positive DNA supercoiling

Reverse gyrase is an ATP-dependent topoisomerase that is unique to hyperthermophilic archaea and eubacteria. The only reverse gyrase structure determined to date has revealed the arrangement of the N-terminal helicase domain and the C-terminal topoisomerase domain that intimately cooperate to generate the unique function of positive DNA supercoiling. Although the structure has elicited hypotheses as to how supercoiling may be achieved, it lacks structural elements important for supercoiling and the molecular mechanism of positive supercoiling is still not clear. We present five structures of authentic Thermotoga maritima reverse gyrase that reveal a first view of two interacting zinc fingers that are crucial for positive DNA supercoiling. The so-called latch domain, which connects the helicase and the topoisomerase domains is required for their functional cooperation and presents a novel fold. Structural comparison defines mobile regions in parts of the helicase domain, including a helical insert and the latch that are likely important for DNA binding during catalysis. We show that the latch, the helical insert and the zinc fingers contribute to the binding of DNA to reverse gyrase and are uniquely placed within the reverse gyrase structure to bind and guide DNA during strand passage. A possible mechanism for positive supercoiling by reverse gyrases is presented.

Reverse gyrase transiently unwinds double-stranded DNA in an ATP-dependent reaction

Reverse gyrase is a unique DNA topoisomerase that catalyzes the introduction of positive supercoils into DNA in an ATP-dependent reaction. It consists of a helicase domain that functionally cooperates with a topoisomerase domain. Different models for the catalytic mechanism of reverse gyrase that predict a central role of the helicase domain have been put forward. The helicase domain acts as a nucleotide-dependent conformational switch that alternates between open and closed states with different affinities for single- and double-stranded DNA. It has been suggested that the helicase domain can unwind double-stranded regions, but helicase activity has not been demonstrated as yet. Here, we show that the isolated helicase domain and full-length reverse gyrase can transiently unwind double-stranded regions in an ATP-dependent reaction. The latch region of reverse gyrase, an insertion into the helicase domain, is required for DNA supercoiling. Strikingly, the helicase domain lacking the latch cannot unwind DNA, linking unwinding to DNA supercoiling. The unwinding activity may provide and stabilize the single-stranded regions required for strand passage by the topoisomerase domain, either de novo or by expanding already existing unpaired regions that may form at high temperatures.

Extreme Thermophiles: Moving beyond single-enzyme biocatalysis

Extremely thermophilic microorganisms have been sources of thermostable and thermoactive enzymes for over 30 years. However, information and insights gained from genome sequences, in conjunction with new tools for molecular genetics, have opened up exciting new possibilities for biotechnological opportunities based on extreme thermophiles that go beyond single-step biotransformations. Although the pace for discovering novel microorganisms has slowed over the past two decades, genome sequence data have provided clues to novel biomolecules and metabolic pathways, which can be mined for a range of new applications. Furthermore, recent advances in molecular genetics for extreme thermophiles have made metabolic engineering for high temperature applications a reality.

Synergic and opposing activities of thermophilic RecQ-like helicase and topoisomerase 3 proteins in Holliday junction processing and replication fork stabilization

RecQ family helicases and topoisomerase 3 enzymes form evolutionary conserved complexes that play essential functions in DNA replication, recombination, and repair, and in vitro, show coordinate activities on model recombination and replication intermediates. Malfunctioning of these complexes in humans is associated with genomic instability and cancer-prone syndromes. Although both RecQ-like and topoisomerase 3 enzymes are present in archaea, only a few of them have been studied, and no information about their functional interaction is available. We tested the combined activities of the RecQ-like helicase, Hel112, and the topoisomerase 3, SsTop3, from the thermophilic archaeon Sulfolobus solfataricus. Hel112 showed coordinate DNA unwinding and annealing activities, a feature shared by eukaryotic RecQ homologs, which resulted in processing of synthetic Holliday junctions and stabilization of model replication forks. SsTop3 catalyzed DNA relaxation and annealing. When assayed in combination, SsTop3 inhibited the Hel112 helicase activity on Holliday junctions and stimulated formation and stabilization of such structures. In contrast, Hel112 did not affect the SsTop3 DNA relaxation activity. RecQ-topoisomerase 3 complexes show structural similarity with the thermophile-specific enzyme reverse gyrase, which catalyzes positive supercoiling of DNA and was suggested to play a role in genome stability at high temperature. Despite such similarity and the high temperature of reaction, the SsTop3-Hel112 complex does not induce positive supercoiling and is thus likely to play different roles. We propose that the interplay between Hel112 and SsTop3 might regulate the equilibrium between recombination and anti-recombination activities at replication forks.

Re-annotation of two hyperthermophilic archaea Pyrococcus abyssi GE5 and Pyrococcus furiosus DSM 3638

Pyrococcus abyssi GE5 (P. aby) and Pyrococcus furiosus DSM 3638 (P. fur) are two model hyperthermophilic archaea. However, their annotations in public databases are unsatisfactory. In this article, the two genomes were re-annotated according to the following steps. (i) All "hypothetical genes" in the original annotation were re-identified based on the Z-curve method, and some of them were recognized as non-coding open reading frames (ORFs). Evidence showed that the recognized non-coding ORFs were highly unlikely to encode proteins. (ii) The translation initiation sites (TISs) of all the annotated genes were re-located, and more than 10% of the TISs were shifted to 5'-upstream or 3'-downstream regions. (iii) The functions of the refined "hypothetical genes" were predicted using sequence alignment tools, more than 200 originally annotated "hypothetical genes" in either of the two hyperthermophiles were assigned functions. A large number of these functions have reference support or experimentally characterized homologues. All the refined information will serve as a valuable resource for research on P. aby and P. fur, which may be helpful in the exploration of thermal adaptation mechanisms. The complete re-annotation files of P. aby and P. fur are available at http://211.69.128.148/download/ .

A new hybrid fractal algorithm for predicting thermophilic nucleotide sequences

Knowledge of thermophilic mechanisms about some organisms whose optimum growth temperature (OGT) ranges from 50 to 80 degree plays a major role in helping design stable proteins. How to predict a DNA sequence to be thermophilic is a long but not fairly resolved problem. Chaos game representation (CGR) can investigate the patterns hiding in DNA sequences, and can visually reveal previously unknown structure. Fractal dimensions are good tools to measure sizes of complex, highly irregular geometric objects. In this paper, we convert every DNA sequence into a high dimensional vector by CGR algorithm and fractal dimension, and then predict the DNA sequence thermostability by these fractal features and support vector machine (SVM). We have conducted experiments on three groups: 17-dimensional vector, 65-dimensional vector, and 257-dimensional vector. Each group is evaluated by the 10-fold cross-validation test. For the results, the group of 257-dimensional vector gets the best results: the average accuracy is 0.9456 and average MCC is 0.8878. The results are also compared with the previous work with single CGR features. The comparison shows the high effectiveness of the new hybrid fractal algorithm.

Sulfolobus mutants, generated via PCR products, which lack putative enzymes of UV photoproduct repair

In order to determine the biological relevance of two S. acidocaldarius proteins to the repair of UV photoproducts, the corresponding genes (Saci_1227 and Saci_1096) were disrupted, and the phenotypes of the resulting mutants were examined by various genetic assays. The disruption used integration by homologous recombination of a functional but heterologous pyrE gene, promoted by short sequences attached to both ends via PCR. The phenotypic analyses of the disruptants confirmed that ORF Saci_1227 encodes a DNA photolyase which functions in vivo, but they could not implicate ORF Saci_1096 in repair of UV- or other externally induced DNA damage despite its similarity to genes encoding UV damage endonucleases. The success of the gene-disruption strategy, which used 5' extensions of PCR primers to target cassette integration, suggests potential advantages for routine construction of Sulfolobus strains.

Systematic computational prediction of protein interaction networks

Determining the network of physical protein associations is an important first step in developing mechanistic evidence for elucidating biological pathways. Despite rapid advances in the field of high throughput experiments to determine protein interactions, the majority of associations remain unknown. Here we describe computational methods for significantly expanding protein association networks. We describe methods for integrating multiple independent sources of evidence to obtain higher quality predictions and we compare the major publicly available resources available for experimentalists to use.

TopR2, the second reverse gyrase of Sulfolobus solfataricus, exhibits unusual properties

Whereas reverse gyrase is considered as a strong marker of thermophily, the function of this peculiar type IA topoisomerase still remains to be elucidated. The archaeon Sulfolobus solfataricus encodes two reverse gyrases, TopR1 and TopR2. This duplication seems to be important because most of Crenarcheota exhibit two copies of reverse gyrase. However, to date, while TopR1 has been well characterized, no characterization of TopR2 has been reported. In this study, we describe for the first time the activity of S. solfataricus TopR2 that appears as a new reverse gyrase. Indeed, in spite of the sequence similarities between TopR1 and TopR2, we evidence unexpected great differences between the two enzymes. While TopR1 exhibits ATP-independent relaxation activity, TopR2 does not, and its activity is strictly dependent on the presence of ATP. Whereas TopR1 is a distributive topoisomerase, TopR2 exhibits an amazing high intrinsic processivity compared to all the topoisomerases studied so far. TopR2 is able to introduce a very high number of positive superturns in DNA, while TopR1 generates weakly positively supercoiled DNA. Finally, TopR2 behaves differently from TopR1 when incubated at different assay temperatures. All the results presented in this study indicate that TopR1 and TopR2 have, in vitro, different activities suggesting different functions in vivo.

Nucleotide-driven conformational changes in the reverse gyrase helicase-like domain couple the nucleotide cycle to DNA processing

Reverse gyrase introduces positive supercoils into DNA in an ATP-dependent process. It has a modular structure comprising a helicase-like and a topoisomerase domain. The helicase-like domain consists of two RecA-like subdomains and thus structurally resembles members of the helicase superfamily 2. It is a nucleotide-dependent switch that alters between an ATP state with a slight preference for dsDNA, and an ADP state with a high preference for ssDNA. Inter-domain communication between the helicase-like and the topoisomerase domain is central for their functional cooperation in reverse gyrase. The latch, an insertion into the helicase-like domain, has been suggested as an important element in coordinating their activities. Here, we have dissected the nucleotide cycle of the reverse gyrase helicase-like domain in the absence and presence of different DNA substrates. With this comprehensive thermodynamic characterization of the nucleotide cycle of the helicase-like domain, in combination with single molecule FRET data on the conformation of the helicase-like domain at all stages of the catalytic cycle, a picture emerges as to how the helicase-like domain may guide ATP-dependent positive supercoiling by reverse gyrase.

Functional evaluation of four putative DNA-binding regions in Thermoanaerobacter tengcongensis reverse gyrase

Reverse gyrase (RG) is an ATP-dependent type I DNA topoisomerase that introduces positive supercoils into DNA in thermophiles. Four regions of RG, i.e., the N-terminal zinc-finger motif, the β-hairpin in subdomain H1, the "latch", and the C-terminal zinc-finger motif, were predicted to be involved in DNA binding previously. In this paper, the functions of these regions in the enzymatic activity were evaluated by mutational analysis of the Thermoanaerobacter tengcongensis reverse gyrase (TtRG). We demonstrated that TtRG exhibited positive-supercoiling activity only at high temperature (>50°C) and low salt concentration (~30 mM NaCl), and three of these four regions (except for the "latch") were involved in DNA binding. Notably, mutations in the "latch" and β-hairpin regions of TtRG strongly impaired the ATPase activity, while mutations in the two zinc-finger motifs dramatically affected its thermal stability besides significant impairment of the DNA-binding ability. Accordingly, all of these four regions were found to be indispensable for the positive-supercoiling activity of TtRG. Taken together, we revealed that these putative DNA-contact regions affect the enzymatic activity of RG in different ways, and provided new insights into the structure and function of RG.

Insights into the evolution of Archaea and eukaryotic protein modifier systems revealed by the genome of a novel archaeal group

The domain Archaea has historically been divided into two phyla, the Crenarchaeota and Euryarchaeota. Although regarded as members of the Crenarchaeota based on small subunit rRNA phylogeny, environmental genomics and efforts for cultivation have recently revealed two novel phyla/divisions in the Archaea; the 'Thaumarchaeota' and 'Korarchaeota'. Here, we show the genome sequence of Candidatus 'Caldiarchaeum subterraneum' that represents an uncultivated crenarchaeotic group. A composite genome was reconstructed from a metagenomic library previously prepared from a microbial mat at a geothermal water stream of a sub-surface gold mine. The genome was found to be clearly distinct from those of the known phyla/divisions, Crenarchaeota (hyperthermophiles), Euryarchaeota, Thaumarchaeota and Korarchaeota. The unique traits suggest that this crenarchaeotic group can be considered as a novel archaeal phylum/division. Moreover, C. subterraneum harbors an ubiquitin-like protein modifier system consisting of Ub, E1, E2 and small Zn RING finger family protein with structural motifs specific to eukaryotic system proteins, a system clearly distinct from the prokaryote-type system recently identified in Haloferax and Mycobacterium. The presence of such a eukaryote-type system is unprecedented in prokaryotes, and indicates that a prototype of the eukaryotic protein modifier system is present in the Archaea.

Gene-centric association analysis for the correlation between the guanine-cytosine content levels and temperature range conditions of prokaryotic species

Background

The environment has been playing an instrumental role in shaping and maintaining the morphological, physiological and biochemical diversities of prokaryotes. It has been debatable whether the whole-genome Guanine-Cytosine (GC) content levels of prokaryotic organisms are correlated with their optimal growth temperatures. Since the GC content is variable within a genome, we here focus on the correlation between the genic GC content levels and the temperature range conditions of prokaryotic organisms.

Results

The GC content levels in the coding regions of four genes were consistently identified as correlated with the temperature range condition when the association analysis was applied to (i) the 722 mesophilic and 93 thermophilic/hyperthermophilic organisms regardless of their phylogeny, oxygen requirement, salinity, or habitat conditions, and (ii) partial lists of organisms when organisms with certain phylogeny, oxygen requirement, salinity or habitat conditions were excluded. These four genes are K01251 (adenosylhomocysteinase), K03724 (DNA repair and recombination proteins), K07588 (LAO/AO transport system kinase), and K09122 (hypothetical protein).

To further validate the identified correlation relationships, we examined to what extent the temperature range condition of an organism can be predicted based on the GC content levels in the coding regions of the selected genes. The 84.52% accuracy for the complete genomes, the 84.09% accuracy for the in-progress genomes, and 82.70% accuracy for the metagenomes, especially when being compared to the 50% accuracy rendered by random guessing, suggested that the temperature range condition of a prokaryotic organism can generally be predicted based on the GC content levels of the selected genomic regions.

Conclusions

The results rendered by various statistical tests and prediction tests indicated that the GC content levels of the coding/non-coding regions of certain genes are highly likely to be correlated with the temperature range conditions of prokaryotic organisms. Therefore, it is promising to carry out “reverse ecology” and to complete the ecological characterizations of prokaryotic organisms, i.e., to infer their temperature range conditions based on the GC content levels of certain genomic regions.

Separate and combined biochemical activities of the subunits of a naturally split reverse gyrase

Reverse gyrase reanneals denatured DNA and induces positive supercoils in DNA, an activity that is critical for life at very high temperatures. Positive supercoiling occurs by a poorly understood mechanism involving the coordination of a topoisomerase domain and a helicase-like domain. In the parasitic archaeon Nanoarchaeum equitans, these domains occur as separate subunits. We express the subunits, and characterize them both in isolation and as a heterodimer. Each subunit tightly associates and interacts with the other. The topoisomerase subunit enhances the catalytic specificity of the DNA-dependent ATPase activity of the helicase-like subunit, and the helicase-like subunit inhibits the relaxation activity of the topoisomerase subunit while promoting positive supercoiling. DNA binding preference for both single- and double-stranded DNA is partitioned between the subunits. Based on a sensitive topological shift assay, the binding preference of helicase-like subunit for underwound DNA is modulated by its binding with ATP cofactor. These results provide new insight into the mechanism of positive supercoil induction by reverse gyrase.

Deep phylogeny--how a tree can help characterize early life on Earth

The Darwinian concept of biological evolution assumes that life on Earth shares a common ancestor. The diversification of this common ancestor through speciation events and vertical transmission of genetic material implies that the classification of life can be illustrated in a tree-like manner, commonly referred to as the Tree of Life. This article describes features of the Tree of Life, such as how the tree has been both pruned and become bushier throughout the past century as our knowledge of biology has expanded. We present current views that the classification of life may be best illustrated as a ring or even a coral with tree-like characteristics. This article also discusses how the organization of the Tree of Life offers clues about ancient life on Earth. In particular, we focus on the environmental conditions and temperature history of Precambrian life and show how chemical, biological, and geological data can converge to better understand this history."You know, a tree is a tree. How many more do you need to look at?"--Ronald Reagan (Governor of California), quoted in the Sacramento Bee, opposing expansion of Redwood National Park, March 3, 1966.

Discontinuity and limited linkage in the homologous recombination system of a hyperthermophilic archaeon

Genetic transformation of Sulfolobus acidocaldarius by a multiply marked pyrE gene provided a high-resolution assay of homologous recombination in a hyperthermophilic archaeon. Analysis of 100 Pyr(+) transformants revealed that this recombination system could transfer each of 23 nonselected base pair substitutions to the recipient chromosome along with the selected marker. In 30% of the recombinants, donor markers were transferred as multiple blocks. In at least 40% of the recombinants, donor markers separated by 5 or 6 bp segregated from each other, whereas similar markers separated by 2 bp did not segregate. Among intermarker intervals, the frequency of recombination tract endpoints varied 40-fold, but in contrast to other recombination systems, it did not correlate with the length of the interval. The average length of donor tracts (161 bp) and the frequent generation of multiple tracts seemed generally consistent with the genetic properties observed previously in S. acidocaldarius conjugation. The efficiency with which short intervals of diverged pyrE sequence were incorporated into the genome raises questions about the threat of ectopic recombination in Sulfolobus spp. mediated by this apparently efficient yet permissive system.

Helical chirality: a link between local interactions and global topology in DNA

DNA supercoiling plays a major role in many cellular functions. The global DNA conformation is however intimately linked to local DNA-DNA interactions influencing both the physical properties and the biological functions of the supercoiled molecule. Juxtaposition of DNA double helices in ubiquitous crossover arrangements participates in multiple functions such as recombination, gene regulation and DNA packaging. However, little is currently known about how the structure and stability of direct DNA-DNA interactions influence the topological state of DNA. Here, a crystallographic analysis shows that due to the intrinsic helical chirality of DNA, crossovers of opposite handedness exhibit markedly different geometries. While right-handed crossovers are self-fitted by sequence-specific groove-backbone interaction and bridging Mg(2+) sites, left-handed crossovers are juxtaposed by groove-groove interaction. Our previous calculations have shown that the different geometries result in differential stabilisation in solution, in the presence of divalent cations. The present study reveals that the various topological states of the cell are associated with different inter-segmental interactions. While the unstable left-handed crossovers are exclusively formed in negatively supercoiled DNA, stable right-handed crossovers constitute the local signature of an unusual topological state in the cell, such as the positively supercoiled or relaxed DNA. These findings not only provide a simple mechanism for locally sensing the DNA topology but also lead to the prediction that, due to their different tertiary intra-molecular interactions, supercoiled molecules of opposite signs must display markedly different physical properties. Sticky inter-segmental interactions in positively supercoiled or relaxed DNA are expected to greatly slow down the slithering dynamics of DNA. We therefore suggest that the intrinsic helical chirality of DNA may have oriented the early evolutionary choices for DNA topology.

Helicase-appended topoisomerases: new insight into the mechanism of directional strand transfer

DNA strand passage through an enzyme-mediated gate is a key step in the catalytic cycle of topoisomerases to produce topological transformations in DNA. In most of the reactions catalyzed by topoisomerases, strand passage is not directional; thus, the enzyme simply provides a transient DNA gate through which DNA transport is allowed and thereby resolves the topological entanglement. When studied in isolation, the type IA topoisomerase family appears to conform to this rule. Interestingly, type IA enzymes can carry out directional strand transport as well. We examined here the biochemical mechanism for directional strand passage of two type IA topoisomerases: reverse gyrase and a protein complex of topoisomerase III alpha and Bloom helicase. These enzymes are able to generate vectorial strand transport independent of the supercoiling energy stored in the DNA molecule. Reverse gyrase is able to anneal single strands, thereby increasing linkage number of a DNA molecule. However, topoisomerase III alpha and Bloom helicase can dissolve DNA conjoined with a double Holliday junction, thus reducing DNA linkage. We propose here that the helicase or helicase-like component plays a determinant role in the directionality of strand transport. There is thus a common biochemical ground for the directional strand passage for the type IA topoisomerases.

Inhibition of translesion DNA polymerase by archaeal reverse gyrase

Reverse gyrase is a unique DNA topoisomerase endowed with ATP-dependent positive supercoiling activity. It is typical of microorganisms living at high temperature and might play a role in maintenance of genome stability and repair. We have identified the translesion DNA polymerase SsoPolY/Dpo4 as one partner of reverse gyrase in the hyperthermophilic archaeon Sulfolobus solfataricus. We show here that in cell extracts, PolY and reverse gyrase co-immunoprecipitate with each other and with the single strand binding protein, SSB. The interaction is confirmed in vitro by far-western and Surface Plasmon Resonance. In functional assays, reverse gyrase inhibits PolY, but not the S. solfataricus B-family DNA polymerase PolB1. Mutational analysis shows that inhibition of PolY activity depends on both ATPase and topoisomerase activities of reverse gyrase, suggesting that the intact positive supercoiling activity is required for PolY inhibition. In vivo, reverse gyrase and PolY are degraded after induction of DNA damage. Inhibition by reverse gyrase and degradation might act as a double mechanism to control PolY and prevent its potentially mutagenic activity when undesired. Inhibition of a translesion polymerase by topoisomerase-induced modification of DNA structure may represent a previously unconsidered mechanism of regulation of these two-faced enzymes.

Detection of biochemical pathways by probabilistic matching of phyletic vectors

A phyletic vector, also known as a phyletic (or phylogenetic) pattern, is a binary representation of the presences and absences of orthologous genes in different genomes. Joint occurrence of two or more genes in many genomes results in closely similar binary vectors representing these genes, and this similarity between gene vectors may be used as a measure of functional association between genes. Better understanding of quantitative properties of gene co-occurrences is needed for systematic studies of gene function and evolution. We used the probabilistic iterative algorithm Psi-square to find groups of similar phyletic vectors. An extended Psi-square algorithm, in which pseudocounts are implemented, shows better sensitivity in identifying proteins with known functional links than our earlier hierarchical clustering approach. At the same time, the specificity of inferring functional associations between genes in prokaryotic genomes is strongly dependent on the pathway: phyletic vectors of the genes involved in energy metabolism and in de novo biosynthesis of the essential precursors tend to be lumped together, whereas cellular modules involved in secretion, motility, assembly of cell surfaces, biosynthesis of some coenzymes, and utilization of secondary carbon sources tend to be identified with much greater specificity. It appears that the network of gene coinheritance in prokaryotes contains a giant connected component that encompasses most biosynthetic subsystems, along with a series of more independent modules involved in cell interaction with the environment.

Integration of phenotypic metadata and protein similarity in Archaea using a spectral bipartitioning approach

In order to simplify and meaningfully categorize large sets of protein sequence data, it is commonplace to cluster proteins based on the similarity of those sequences. However, it quickly becomes clear that the sequence flexibility allowed a given protein varies significantly among different protein families. The degree to which sequences are conserved not only differs for each protein family, but also is affected by the phylogenetic divergence of the source organisms. Clustering techniques that use similarity thresholds for protein families do not always allow for these variations and thus cannot be confidently used for applications such as automated annotation and phylogenetic profiling. In this work, we applied a spectral bipartitioning technique to all proteins from 53 archaeal genomes. Comparisons between different taxonomic levels allowed us to study the effects of phylogenetic distances on cluster structure. Likewise, by associating functional annotations and phenotypic metadata with each protein, we could compare our protein similarity clusters with both protein function and associated phenotype. Our clusters can be analyzed graphically and interactively online.

Reverse gyrase and genome stability in hyperthermophilic organisms

Reverse gyrase is a DNA topoisomerase that is peculiar in many aspects: it has the unique ability to introduce positive supercoils into DNA molecules; it comprises a type IA topoisomerase fused to a helicase-like domain; although it is a type IA topoisomerase, its reaction is ATP-dependent; and it is the only hyperthermophile-specific protein. All these features have made reverse gyrase the subject of biochemical, structural and functional studies, although they have not shed complete light on the evolution, mechanism and function of this distinctive enzyme. In the present article, we review the latest progress on structure-function relationships of reverse gyrase, and discuss old and recent data linking reverse gyrase to DNA stability, protection and repair in hyperthermophilic organisms.

Phylogenomics of DNA topoisomerases: their origin and putative roles in the emergence of modern organisms

Topoisomerases are essential enzymes that solve topological problems arising from the double-helical structure of DNA. As a consequence, one should have naively expected to find homologous topoisomerases in all cellular organisms, dating back to their last common ancestor. However, as observed for other enzymes working with DNA, this is not the case. Phylogenomics analyses indicate that different sets of topoisomerases were present in the most recent common ancestors of each of the three cellular domains of life (some of them being common to two or three domains), whereas other topoisomerases families or subfamilies were acquired in a particular domain, or even a particular lineage, by horizontal gene transfers. Interestingly, two groups of viruses encode topoisomerases that are only distantly related to their cellular counterparts. To explain these observations, we suggest that topoisomerases originated in an ancestral virosphere, and that various subfamilies were later on transferred independently to different ancient cellular lineages. We also proposed that topoisomerases have played a critical role in the origin of modern genomes and in the emergence of the three cellular domains.

Complete genome sequence of the anaerobic, protein-degrading hyperthermophilic crenarchaeon Desulfurococcus kamchatkensis

Desulfurococcus kamchatkensis is an anaerobic organotrophic hyperthermophilic crenarchaeon isolated from a terrestrial hot spring. Its genome consists of a single circular chromosome of 1,365,223 bp with no extrachromosomal elements. A total of 1,474 protein-encoding genes were annotated, among which 205 are exclusive for D. kamchatkensis. The search for a replication origin site revealed a single region coinciding with a global extreme of the nucleotide composition disparity curve and containing a set of crenarchaeon-type origin recognition boxes. Unlike in most archaea, two genes encoding homologs of the eukaryotic initiator proteins Orc1 and Cdc6 are located distantly from this site. A number of mobile elements are present in the genome, including seven transposons representing IS607 and IS200/IS605 families and multiple copies of miniature inverted repeat transposable elements. Two large clusters of regularly interspaced repeats are present; none of the spacer sequences matches known archaeal extrachromosomal elements, except one spacer matches the sequence of a resident gene of D. kamchatkensis. Many of the predicted metabolic enzymes are associated with the fermentation of peptides and sugars, including more than 30 peptidases with diverse specificities, a number of polysaccharide degradation enzymes, and many transporters. Consistently, the genome encodes both enzymes of the modified Embden-Meyerhof pathway of glucose oxidation and a set of enzymes needed for gluconeogenesis. The genome structure and content reflect the organism's nutritionally diverse, competitive natural environment, which is periodically invaded by viruses and other mobile elements.

Novel multiprotein complexes identified in the hyperthermophilic archaeon Pyrococcus furiosus by non-denaturing fractionation of the native proteome

Virtually all cellular processes are carried out by dynamic molecular assemblies or multiprotein complexes, the compositions of which are largely undefined. They cannot be predicted solely from bioinformatics analyses nor are there well defined techniques currently available to unequivocally identify protein complexes (PCs). To address this issue, we attempted to directly determine the identity of PCs from native microbial biomass using Pyrococcus furiosus, a hyperthermophilic archaeon that grows optimally at 100 degrees C, as the model organism. Novel PCs were identified by large scale fractionation of the native proteome using non-denaturing, sequential column chromatography under anaerobic, reducing conditions. A total of 967 distinct P. furiosus proteins were identified by mass spectrometry (nano LC-ESI-MS/MS), representing approximately 80% of the cytoplasmic proteins. Based on the co-fractionation of proteins that are encoded by adjacent genes on the chromosome, 106 potential heteromeric PCs containing 243 proteins were identified, only 20 of which were known or expected. In addition to those of unknown function, novel and uncharacterized PCs were identified that are proposed to be involved in the metabolism of amino acids (10), carbohydrates (four), lipids (two), vitamins and metals (three), and DNA and RNA (nine). A further 30 potential PCs were classified as tentative, and the remaining potential PCs (13) were classified as weakly interacting. Some major advantages of native biomass fractionation for PC identification are that it provides a road map for the (partial) purification of native forms of novel and uncharacterized PCs, and the results can be utilized for the recombinant production of low abundance PCs to provide enough material for detailed structural and biochemical analyses.

Genomics of bacteria and archaea: the emerging dynamic view of the prokaryotic world

The first bacterial genome was sequenced in 1995, and the first archaeal genome in 1996. Soon after these breakthroughs, an exponential rate of genome sequencing was established, with a doubling time of approximately 20 months for bacteria and approximately 34 months for archaea. Comparative analysis of the hundreds of sequenced bacterial and dozens of archaeal genomes leads to several generalizations on the principles of genome organization and evolution. A crucial finding that enables functional characterization of the sequenced genomes and evolutionary reconstruction is that the majority of archaeal and bacterial genes have conserved orthologs in other, often, distant organisms. However, comparative genomics also shows that horizontal gene transfer (HGT) is a dominant force of prokaryotic evolution, along with the loss of genetic material resulting in genome contraction. A crucial component of the prokaryotic world is the mobilome, the enormous collection of viruses, plasmids and other selfish elements, which are in constant exchange with more stable chromosomes and serve as HGT vehicles. Thus, the prokaryotic genome space is a tightly connected, although compartmentalized, network, a novel notion that undermines the ‘Tree of Life’ model of evolution and requires a new conceptual framework and tools for the study of prokaryotic evolution.

The reverse gyrase helicase-like domain is a nucleotide-dependent switch that is attenuated by the topoisomerase domain

Reverse gyrase is a topoisomerase that introduces positive supercoils into DNA in an ATP-dependent manner. It is unique to hyperthermophilic archaea and eubacteria, and has been proposed to protect their DNA from damage at high temperatures. Cooperation between its N-terminal helicase-like and the C-terminal topoisomerase domain is required for positive supercoiling, but the precise role of the helicase-like domain is currently unknown. Here, the characterization of the isolated helicase-like domain from Thermotoga maritima reverse gyrase is presented. We show that the helicase-like domain contains all determinants for nucleotide binding and ATP hydrolysis. Its intrinsic ATP hydrolysis is significantly stimulated by ssDNA, dsDNA and plasmid DNA. During the nucleotide cycle, the helicase-like domain switches between high- and low-affinity states for dsDNA, while its affinity for ssDNA in the ATP and ADP states is similar. In the context of reverse gyrase, the differences in DNA affinities of the nucleotide states are smaller, and the DNA-stimulated ATPase activity is strongly reduced. This inhibitory effect of the topoisomerase domain decelerates the progression of reverse gyrase through the nucleotide cycle, possibly providing optimal coordination of ATP hydrolysis with the complex reaction of DNA supercoiling.

Transcriptional analysis of the two reverse gyrase encoding genes of Sulfolobus solfataricus P2 in relation to the growth phases and temperature conditions

Sulfolobus solfataricus, a hyperthermophilic crenarchaeon, contains two genes encoding reverse gyrases, topR1 and topR2. The steady-state level of their transcripts were quantified during the growth phases for cells maintained either at 72, or 80 degrees C, and after temperature changes from one to the other temperature. The transcripts of both genes are weakly expressed, but the highest level is observed in actively dividing cells, and is almost undetectable in cells in decline phase. During the temperature shift experiments, there is no significant topR2 variation. By contrast, there is a maximum 2.4-fold increase in topR1 transcripts within 30 min after the downshift. After 1 h, the transcript level reaches the level characteristic of cells adapted to the new temperature. After an upward shift, the topR1 expression pattern is inversely regulated with a transient decrease with the same time course. The topR1 expression profile is completely different from that of topR2 after temperature shift experiments; this suggests a different regulation process for the two reverse gyrase genes. The fine tuning of the topR1 transcript expression within a short interval of time after a temperature shift illustrates a rapid adaptation response to temperature change.

The rate and character of spontaneous mutation in Thermus thermophilus

Selection of spontaneous, loss-of-function mutations at two chromosomal loci (pyrF and pyrE) enabled the first molecular-level analysis of replication fidelity in the extremely thermophilic bacterium Thermus thermophilus. Two different methods yielded similar mutation rates, and mutational spectra determined by sequencing of independent mutants revealed a variety of replication errors distributed throughout the target genes. The genomic mutation rate estimated from these targets, 0.00097 +/- 0.00052 per replication, was lower than corresponding estimates from mesophilic microorganisms, primarily because of a low rate of base substitution. However, both the rate and spectrum of spontaneous mutations in T. thermophilus resembled those of the thermoacidophilic archaeon Sulfolobus acidocaldarius, despite important molecular differences between these two thermophiles and their genomes.

Dissection of reverse gyrase activities: insight into the evolution of a thermostable molecular machine

Reverse gyrase is a peculiar DNA topoisomerase, specific of thermophilic microorganisms, which induces positive supercoiling into DNA molecules in an ATP-dependent reaction. It is a modular enzyme and comprises an N-terminal helicase-like module fused to a C-terminal topoisomerase IA-like domain. The exact molecular mechanism of this unique reaction is not understood, and a fundamental mechanistic question is how its distinct steps are coordinated. We studied the cross-talk between the components of this molecular motor and probed communication between the DNA-binding sites and the different activities (DNA relaxation, ATP hydrolysis and positive supercoiling). We show that the isolated ATPase and topoisomerase domains of reverse gyrase form specific physical interactions, retain their own DNA binding and enzymatic activities, and when combined cooperate to achieve the unique ATP-dependent positive supercoiling activity. Our results indicate a mutual effect of both domains on all individual steps of the reaction. The C-terminal domain shows ATP-independent topoisomerase activity, which is repressed by the N-terminal domain in the full-length enzyme; experiments with the isolated domains showed that the C-terminal domain has stimulatory influence on the ATPase activity of the N-terminal domain. In addition, the two domains showed a striking reciprocal thermostabilization effect.

Mutational analysis of the helicase-like domain of Thermotoga maritima reverse gyrase

Reverse gyrase is a unique type IA topoisomerase that is able to introduce positive supercoils into DNA in an ATP-dependent process. ATP is bound to the helicase-like domain of the enzyme that contains most of the conserved motifs found in helicases of the SF1 and SF2 superfamilies. In this paper, we have investigated the role of the conserved helicase motifs I, II, V, VI, and Q by generating mutants of the Thermotoga maritima reverse gyrase. We show that mutations in motifs I, II, V, and VI completely eliminate the supercoiling activity of reverse gyrase and that a mutation in the Q motif significantly reduces this activity. Further analysis revealed that for most mutants, the DNA binding and cleavage properties are not significantly changed compared with the wild type enzyme, whereas their ATPase activity is impaired. These results clearly show that the helicase motifs are tightly involved in the coupling of ATP hydrolysis to the topoisomerase activity. The zinc finger motif located at the N-terminal end of reverse gyrases was also mutated. Our results indicate that this motif plays an important role in DNA binding.

How does gene expression level contribute to thermophilic adaptation of prokaryotes? An exploration based on predictors

By analyzing the predicted gene expression levels of 33 prokaryotes with living temperature span from <10 degrees C to >100 degrees C, a universal positive correlation was found between the percentage of predicted highly expressed genes and the organisms' optimal growth temperature. A physical interpretation of the correlation revealed that highly expressed genes are statistically more thermostable than lowly expressed genes. These findings show the possibility of the significant contribution of gene expression level to the prokaryotic thermal adaptation and provide evidence for the translational selection pressure on the thermostability of natural proteins during evolution.

Transcriptome analysis of infection of the archaeon Sulfolobus solfataricus with Sulfolobus turreted icosahedral virus

Microarray analysis of infection by Sulfolobus turreted icosahedral virus (STIV) revealed insights into the timing and extent of virus transcription, as well as differential regulation of host genes. Using a microarray containing genes from both the host and the virus, the infection cycle of STIV was studied. Following infection of Sulfolobus solfataricus strain 2-2-12 with STIV, transcription of virus genes was first detected at 8 h postinfection (p.i.), with a peak at 24 h p.i. Lysis of cells was first detected at 32 h p.i. There was little temporal control of the transcription of virus genes, although the three open reading frames on the noncoding strand were transcribed later in the infection process. During the infection, 177 host genes were determined to be differentially expressed, with 124 genes up-regulated and 53 genes down-regulated. The up-regulated genes were dominated by genes associated with DNA replication and repair and those of unknown function, while the down-regulated genes, mostly detected at 32 h p.i., were associated with energy production and metabolism. Examination of infected cells by transmission electron microscopy revealed alterations in cell ultrastructure consistent with the microarray analysis. The observed patterns of transcription suggest that up-regulated genes are likely used by the virus to reprogram the cell for virus replication, while the down-regulated genes reflect the imminent lysis of the cells.

Clusters of orthologous genes for 41 archaeal genomes and implications for evolutionary genomics of archaea

Background

An evolutionary classification of genes from sequenced genomes that distinguishes between orthologs and paralogs is indispensable for genome annotation and evolutionary reconstruction. Shortly after multiple genome sequences of bacteria, archaea, and unicellular eukaryotes became available, an attempt on such a classification was implemented in Clusters of Orthologous Groups of proteins (COGs). Rapid accumulation of genome sequences creates opportunities for refining COGs but also represents a challenge because of error amplification. One of the practical strategies involves construction of refined COGs for phylogenetically compact subsets of genomes.

Results

New Archaeal Clusters of Orthologous Genes (arCOGs) were constructed for 41 archaeal genomes (13 Crenarchaeota, 27 Euryarchaeota and one Nanoarchaeon) using an improved procedure that employs a similarity tree between smaller, group-specific clusters, semi-automatically partitions orthology domains in multidomain proteins, and uses profile searches for identification of remote orthologs. The annotation of arCOGs is a consensus between three assignments based on the COGs, the CDD database, and the annotations of homologs in the NR database. The 7538 arCOGs, on average, cover ~88% of the genes in a genome compared to a ~76% coverage in COGs. The finer granularity of ortholog identification in the arCOGs is apparent from the fact that 4538 arCOGs correspond to 2362 COGs; ~40% of the arCOGs are new. The archaeal gene core (protein-coding genes found in all 41 genome) consists of 166 arCOGs. The arCOGs were used to reconstruct gene loss and gene gain events during archaeal evolution and gene sets of ancestral forms. The Last Archaeal Common Ancestor (LACA) is conservatively estimated to possess 996 genes compared to 1245 and 1335 genes for the last common ancestors of Crenarchaeota and Euryarchaeota, respectively. It is inferred that LACA was a chemoautotrophic hyperthermophile that, in addition to the core archaeal functions, encoded more idiosyncratic systems, e.g., the CASS systems of antivirus defense and some toxin-antitoxin systems.

Conclusion

The arCOGs provide a convenient, flexible framework for functional annotation of archaeal genomes, comparative genomics and evolutionary reconstructions. Genomic reconstructions suggest that the last common ancestor of archaea might have been (nearly) as advanced as the modern archaeal hyperthermophiles. ArCOGs and related information are available at: .

Reviewers

This article was reviewed by Peer Bork, Patrick Forterre, and Purificacion Lopez-Garcia.

Widespread distribution of archaeal reverse gyrase in thermophilic bacteria suggests a complex history of vertical inheritance and lateral gene transfers

Reverse gyrase, an enzyme of uncertain funtion, is present in all hyperthermophilic archaea and bacteria. Previous phylogenetic studies have suggested that the gene for reverse gyrase has an archaeal origin and was transferred laterally (LGT) to the ancestors of the two bacterial hyperthermophilic phyla, Thermotogales and Aquificales. Here, we performed an in-depth analysis of the evolutionary history of reverse gyrase in light of genomic progress. We found genes coding for reverse gyrase in the genomes of several thermophilic bacteria that belong to phyla other than Aquificales and Thermotogales. Several of these bacteria are not, strictly speaking, hyperthermophiles because their reported optimal growth temperatures are below 80 degrees C. Furthermore, we detected a reverse gyrase gene in the sequence of the large plasmid of Thermus thermophilus strain HB8, suggesting a possible mechanism of transfer to the T. thermophilus strain HB8 involving plasmids and transposases. The archaeal part of the reverse gyrase tree is congruent with recent phylogenies of the archaeal domain based on ribosomal proteins or RNA polymerase subunits. Although poorly resolved, the complete reverse gyrase phylogeny suggests an ancient acquisition of the gene by bacteria via one or two LGT events, followed by its secondary distribution by LGT within bacteria. Finally, several genes of archaeal origin located in proximity to the reverse gyrase gene in bacterial genomes have bacterial homologues mostly in thermophiles or hyperthermophiles, raising the possibility that they were co-transferred with the reverse gyrase gene. Our new analysis of the reverse gyrase history strengthens the hypothesis that the acquisition of reverse gyrase may have been a crucial evolutionary step in the adaptation of bacteria to high-temperature environments. However, it also questions the role of this enzyme in thermophilic bacteria and the selective advantage its presence could provide.

Robust determinants of thermostability highlighted by a codon frequency index capable of discriminating thermophilic from mesophilic genomes

Can genome analysis tell us about the lifestyle of an organism? We ask this question considering a thorough cross comparison of thermophilic and mesophilic genomes, since presently the number of available genomes is enough to ensure statistical significance of the results. We analyze, by means of principal component analysis (PCA), the codon composition of a database comprising 116 genomes, selected so as to include one species for each genus and show that a cross genomic approach can allow the extraction of common determinants of thermostability at the genome level. The results of our analysis indicate that all the known features of thermostability can be found in the 64 component loadings of the second principal axis of PCA. By this, we develop an index of thermostability whose discriminative power between mesophiles and thermophiles scores with 98% accuracy at the genome level and with 95% accuracy at the protein sequence level. We also prove that these results are not due to phylogenetic differences between archaea and bacteria.