Comparative genomics of Thermus thermophilus and Deinococcus radiodurans: divergent routes of adaptation to thermophily and radiation resistance
© Omelchenko et al; licensee BioMed Central Ltd. 2005
Received: 23 August 2005
Accepted: 20 October 2005
Published: 20 October 2005
Thermus thermophilus and Deinococcus radiodurans belong to a distinct bacterial clade but have remarkably different phenotypes. T. thermophilus is a thermophile, which is relatively sensitive to ionizing radiation and desiccation, whereas D. radiodurans is a mesophile, which is highly radiation- and desiccation-resistant. Here we present an in-depth comparison of the genomes of these two related but differently adapted bacteria.
By reconstructing the evolution of Thermus and Deinococcus after the divergence from their common ancestor, we demonstrate a high level of post-divergence gene flux in both lineages. Various aspects of the adaptation to high temperature in Thermus can be attributed to horizontal gene transfer from archaea and thermophilic bacteria; many of the horizontally transferred genes are located on the single megaplasmid of Thermus. In addition, the Thermus lineage has lost a set of genes that are still present in Deinococcus and many other mesophilic bacteria but are not common among thermophiles. By contrast, Deinococcus seems to have acquired numerous genes related to stress response systems from various bacteria. A comparison of the distribution of orthologous genes among the four partitions of the Deinococcus genome and the two partitions of the Thermus genome reveals homology between the Thermus megaplasmid (pTT27) and Deinococcus megaplasmid (DR177).
After the radiation from their common ancestor, the Thermus and Deinococcus lineages have taken divergent paths toward their distinct lifestyles. In addition to extensive gene loss, Thermus seems to have acquired numerous genes from thermophiles, which likely was the decisive contribution to its thermophilic adaptation. By contrast, Deinococcus lost few genes but seems to have acquired many bacterial genes that apparently enhanced its ability to survive different kinds of environmental stresses. Notwithstanding the accumulation of horizontally transferred genes, we also show that the single megaplasmid of Thermus and the DR177 megaplasmid of Deinococcus are homologous and probably were inherited from the common ancestor of these bacteria.
Deinococcus spp. and Thermus spp. are believed to belong to a distinct branch of bacteria called the Deinococcus-Thermus group . The common origin of these bacteria is supported by the fact that they consistently form a strongly supported clade in phylogenetic trees of ribosomal RNAs and several conserved proteins including ribosomal proteins, RNA polymerase subunits, RecA, and others [1–7]. The Thermus genus currently consists of 14 thermophilic species, whereas the Deinococcus genus includes at least eleven, mostly mesophilic species known for their extreme resistance to γ-irradiation and other agents causing DNA damage, particularly, desiccation [7, 8]. Interestingly, two new species, D. frigens sp. nov., D. saxicola, have been isolated recently from Antarctic rock and soil samples . D. geothermalis and D. murrayi, are considered to be thermophilic (topt~45–50°C) .
T. thermophilus (TT) (strain HB27, ATCC BAA-163)and D. radiodurans (DR) (strain R1, ATCC BAA-816) have similar general physiology, both being catalase-positive, red-pigmented, non-sporulating, aerobic chemoorganoheterotrophs . However, the two organisms are dramatically different in terms of stress resistance: DR is one of the most resistant to radiation and desiccation among the characterized organisms and, generally, can survive diverse types of oxidative stress, whereas TT is a thermophile that thrives under thermal stress conditions but is relatively sensitive to radiation and other forms of oxidative stress. Recently, the genomes of T. thermophilus HB27 and HB8 [12, 13] became available for comparison with the previously sequenced genome of D. raduodurans R1 .
Despite extensive research, genetic systems underlying both thermophilic adaptations and radiation resistance remain poorly understood. Attempts have been made to detect "thermophilic determinants" in the proteomes of thermophiles using a comparative-genomic approach. This resulted in the delineation of a set of proteins that might be associated with the thermophilic phenotype, although most of these are significantly enriched in thermophiles but not unique to these organisms [15–17]. Other distinctions between thermophilic and mesophilic proteins might be due to differences in their structural properties, such as different amino acid compositions, loop lengths, number of salt bridges, strength of hydrophobic interaction, number of disulfide bonds, and other features [18–25].
Various hypotheses also have been proposed to explain radiation resistance, some postulating the existence of specialized genetic systems, particularly those for DNA repair and stress response [7, 26]. Recently, however, alternative possibilities have been advanced. For example, the post-irradiation adjustment of metabolism of D. radiodurans might prevent production of reactive oxygen species (ROS) by decreasing the number of reactions involving oxygen [27, 28], and high intracellular manganese concentrations of Deinococcus spp. might help scavenge ROS generated during irradiation and post-irradiation recovery [28, 29]. However, these explanations of radiation resistance have received little direct support from comparative-genomic analyses [7, 27, 28].
Several evolutionary processes could potentially contribute to the genome differentiation of TT and DR subsequent to the divergence from the common ancestor: (i) differential gene loss and gain, (ii) acquisition of genes via horizontal gene transfer (HGT) which may be followed by loss of the ancestral orthologous gene (xenologous gene displacement (XGD)), (iii) lineage-specific expansion of paralogous gene families by duplication and/or acquisition of paralogs via HGT; (IV) modification of amino acid composition that could affect protein stability. Here, we experimentally characterize radiation and desiccation resistance of TT in comparison to DR, and, assess the contribution of different evolutionary processes to distinct adaptations of TT and DR, using a variety of comparative-genomic approaches and phylogenetic analysis. We identify the unique feature of the gene repertoires of TT and DR that might contribute to these phenotypic differences, which could be the subject of further experimental work. In addition, we describe the results of a detailed analysis of the proteins predicted to be involved in DNA repair and stress response functions, which are particularly relevant for adaptive evolution of resistance phenotypes.
Results and discussion
Experimental characterization of resistance to gamma-radiation and desiccation, and determination of intracellular Mn/Fe ratio for T. thermophilus
We recently reported a trend of intracellular Mn/Fe concentration ratios in bacterial IR resistance, where very high and very low Mn/Fe ratios correlated with very high and very low resistances, respectively . We have also shown that growing D. radiodurans in conditions which limited Mn(II) accumulation, significantly lowered the cells' IR resistance . These observations led to the hypothesis that the ratio of Mn to Fe in a cell might determine the relative abundance of different ROS induced during exposure to and recovery from IR [28, 29]. At high concentrations, Mn(II) can act as true catalyst of the dismutation of superoxide (O2•-), with Mn cycling between the divalent and trivalent states; Mn redox-cycling scavenges both O2•- and hydrogen peroxide . In this context, we determined the intracellular Mn/Fe ratio of TT using an inductively coupled plasma-mass spectrometry method (ICP-MS), as previously described . TT cells contained 0.211 (± 0.0254) nmol Mn/mg protein and 4.54 (± 0.778) nmol Fe/mg protein. The intracellular Mn/Fe ratio of TT (D10, ~0.8 kGy) is 0.047 compared to 0.0072 for E. coli (D10, ~0.7 kGy), 0.24 for D. radiodurans (D10, ~16 kGy), and <0.0001 for Pseudomonas putida (D10, ~0.25 kGy) . Thus, TT appears to be somewhat more sensitive to acute IR and desiccation than predicted by its Mn/Fe ratio, suggesting the possibility of more complex relationships between these two variables. Notably, scavenging of O2•- by Mn(II) is highly dependent on the availability of H2O2, which in TT would be expected to become limiting during recovery at 65°C because of thermal decomposition ; inefficient Mn redox-cycling can lead to Mn(III) accumulation, which is cytotoxic. Both DR and TT encode an ABC- type Mn transporter and a transcriptional regulator that probably regulates Mn homeostasis. Additionally, DR has a NRAMP family Mn transporter, for which there is no ortholog in TT.
Reconstruction of the gene-content tree and the gene repertoire of the common ancestor of Deinococcus and Thermus
Concordant and discordant phyletic patterns between DR, TT, Aquifex aeolicus (AA) and Thermotoga maritima (TM)
The majority of the genes shared by TT and DR encode house-keeping proteins and are widespread in bacteria. Among those, 14 COGs are unusual in that they are not found in any other bacteria but instead are shared by TT and DR with archaea or eukaryotes. This set includes 8 subunits of the Archaeal/vacuolar-type H+-ATPase and six other COGs that consist of characteristic archaeal genes (phosphoglycerate mutase, COG3635; 2-phosphoglycerate kinase, COG2074; SAR1-like GTPase, COG1100; GTP:adenosylcobinamide-phosphate guanylyltransferase, COG2266; Predicted membrane protein, COG3374 ; Predicted DNA modification methylase, COG1041).
As noticed previously, some DR genes that belong to families well-represented in both bacteria and archaea showed clear archaeal affinity . To assess how many genes of apparent "thermophilic" descent (either archaeal or bacterial) might already have been present already in the common DR-TT ancestor, we performed phylogenetic analysis of genes that were assigned to the DR-TT ancestor and had 3 of the 5 best hits to the genes from thermophiles (both archaeal and bacterial; see Additional file 1: "Phylogenetic analysis of the genes of the reconstructed DR-TT common ancestor"). We found that at least 122 genes (~10% of the predicted gene repertoire of the DR-TT common ancestor) of the originally selected 205 genes showed an affinity to thermophilic species, i.e., either a branch of two orthologous genes from DR and TT, or a DR or TT gene (in cases when the respective ortholog apparently was lost in the other lineage) clustered with thermophiles (see Additional file 1, table 1S; and Additional file 6). Due to the fact that many tree topologies are highly perturbed by multiple HGT events and may be inaccurate due to differences in evolutionary rates between lineages, this is only a rough estimate of the number of "thermophilic" genes in the DR-TT ancestor. In particular, it cannot be ruled out that some of the ancestral "thermophilic" genes, which are currently present in TT but not DR (10 such genes out of 122 tested genes), have been acquired by TT via XGD (see Additional file 1, table 1S). Taken together, these observations suggest the possibility of ecological contacts between the DR-TT ancestor and hyperthermophilic archaea and/or bacteria, leading to substantial acquisition of "thermophilic" genes via HGT.
Gene gain and loss in Thermus
Our reconstruction of the evolutionary events that occurred after the divergence of the TT and DR lineages from the common ancestor delineated the sets of genes that likely have been gained and lost by each lineage (Figure 4). We first consider in greater detail the pattern of gene loss and gain in TT. The absence of certain metabolic genes in TT creates gaps in its metabolic pathways, some of which are essential. However, TT is capable of synthesizing all amino acids, nucleotides and a majority of cofactors, suggesting that the gaps are filled by analogous or at least non-orthologous enzymes. Several such cases have been described. For example, TT and DR encode unrelated thymidylate synthases, DR2630 (COG0207) and TTC0731 (COG1351). The classical, folate-dependent thymidylate synthase (COG0207) present in DR is probably ancestral in bacteria and apparently was displaced via HGT in the Thermus lineage by the flavin-dependent thymidylate synthase, typical of archaea and bacterial thermophiles. In other cases, the substituted analogous enzymes or pathways remain uncharacterized. For instance, the displacement of the folate-dependent thymidylate synthase with the flavin-dependent type in the TT lineage and in other bacteria and archaea correlates with the apparent loss of dihydrofolate reductase (folA, COG0262), which catalyzes the last step of the tetrahydrofolate biosynthesis pathway. Since tetrahydrofolate is an essential cofactor, a displacement appears most likely. Recently, it has been shown that the halobacterial FolP-FolC fusion protein complemented a Haloferax volcanii folA mutant . Thus, it appears likely that FolC and FolP in TT and other organisms complement the activity of FolA although the existence of an unrelated, as yet uncharacterized dihydrofolate reductase cannot be ruled out.
In addition, TT does not encode two enzymes for pyridoxal phosphate biosynthesis (pdxK, COG0259 and pdxH, COG2240), while DR has a complete set of enzymes of this pathway. Our reconstructions suggest that pdxK was likely lost in the TT lineage, whereas DR probably independently acquired pdxH. However, the pathway is likely to be functional in both organisms. Since similar gaps in the pyridoxal phosphate biosynthesis are seen in a variety of prokaryotes , it appears that, for at least some steps of this pathway, there exists a set of distinct enzymes which remain unidentified.
Some systems apparently were completely lost in the TT lineage. These include the urease complex, the ramnose metabolism pathway, acyl CoA:acetate/3-ketoacid CoA transferase, fructose transport and utilization, and glycerol metabolism. Notably, most of these systems are also absent in thermophilic bacteria and many thermophilic archaea.
The Sox-like sulfur oxidation system is among the group of genes that were apparently acquired in the TT lineage. The TT Sox operon is partly similar to the one identified in AA (see Additional file 1, Figure 2S), and might have been horizontally transferred between the AA and TT lineages with subsequent local rearrangements. The presence of the Sox operon in TT suggests that this bacterium can use reduced sulfur compounds as a source of energy and sulfur. Another system likely acquired by TT is lactose utilization (at least three proteins: LacZ, COG3250; GalK, COG0153; GalT COG1085; GalA, COG3345). This system is also present in TM, which indicates that sugars can be utilized as carbon sources by thermophilic bacteria.
Gene gain and loss in Deinococcus
The DR lineage has apparently acquired many more genes than it has lost (Figure 4). The majority of the genes lost by DR encode enzymes of energy metabolism and biosynthesis of cofactors. One example is the loss of the three subunits of pyruvate:ferredoxin oxidoreductase, which is one of the several known systems for pyruvate oxidation, a key reaction of central metabolism. Another example of gene loss in DR involves the three subunits of NAD/NADP transhydrogenase, which is responsible for energy-dependent reduction of NADP+ . In addition, the DR lineage lost four enzymes of NAD biosynthesis and six enzymes of cobalamine biosynthesis, and consistently, DR is dependent on an exogenous source of NAD for growth [28, 47].
A conspicuous number of genes apparently acquired by the DR lineage encode systems of protein degradation and amino acid catabolism (e.g., urease, DRA0311-DRA0319, and a predicted urea transporter, DRA0320-DRA0324; histidine degradation system, DRA0147-DRA0150; monoamine oxidase, DRA0274; lysine 2,3-aminomutase, DRA0027; kynureninase and tryptophan-2,3-dioxygenase, DRA0338-DRA0339; peptidase E, DR1070, and carboxypeptidase C, DR0964; and D-aminopeptidase, DR1843 (see Additional file 4). A similar trend is observed for the expansion of several protein families in DR, such as secreted subtilisin-like proteases (see below). Additionally, DR acquired two three-subunit complexes of aerobic-type carbon monoxide dehydrogenase (DRA0231-DRA0233 and DRA0235-DRA0237); oxidation of CO by this enzyme might be used as an energy source as shown for some bacteria . Acquisition and expansion of these metabolic systems, together with the loss of certain biosynthetic capabilities, supports the possibility that metabolic restructuring could impact oxidative stress resistance in DR by decreasing the need for high-energy-dependent cellular activities. Energy production through the respiratory chain is one major source of free radicals in the cell .
DR has many more genes for proteins involved in inorganic ion transport and metabolism than TT. In particular, DR has acquired the multisubunit Na+/H+ antiporter (7 genes), K+-transporting ATPase (3 genes), and the FeoA/FeoB Fe transport system. This abundance of ion transport systems might be indirectly linked to oxidative stress resistance through regulation of membrane ion gradients and Mn/Fe homeostasis (see above).
DR is more dependent than TT on peptide-derived growth substrates  and has a more complex stress response circuitry. Consistent with this, the signal transduction systems of DR, as predicted by genome analysis, are considerably more elaborate than those of TT. In particular, DR has at least 33 COGs (15 probably acquired after the divergence from the common ancestor with TT) related to signal-transduction functions that are not represented in TT as compared to 12 (5 acquired) such COGs in TT. Furthermore, although most of the signal-transduction domains are shared by DR and TT, the domain architectures of the respective multidomain proteins are completely different ( and KSM, unpublished observations).
Among the genes apparently acquired by DR, two encode multidomain proteins containing distinct periplasmic ligand-binding sensor domains (DRA0202, COG5278 and DR1174, COG3614). Another protein, DRA0204, contains the CHASE3 domain  and is located in a predicted operon with superoxide dismutase (DRA0202), indicating a function in oxidative stress response. The protein DR0724 contains the SARP domain which is involved in apoptosis-related signaling pathways in eukaryotes  but its function in bacteria is unknown. The roles of the other signal transduction proteins of DR are even less clear, with the notable exception of a phytochrome-like protein (DRA0050) that apparently was acquired by DR from a bacterial source and has been implicated in UV resistance .
Additionally, DR has many genes (25 COGs) encoding systems for microbial defense; TT has only 14 COGs in this category, 13 of which are shared with DR. At least 7 genes for restriction-modification system subunits were specifically acquired by the DR lineage, along with several antibiotic-resistance enzymes. This difference might be linked to the reduced metabolic capabilities of DR, which is dependent on nutrient-rich conditions for growth and, perhaps, encounters more microbial species than TT.
The previous analysis of the DR genome revealed 15 genes that appear to have been horizontally transferred from unexpected sources, such as eukaryotes and viruses ; only two of these 15 genes are present in TT, the desiccation-related protein of the ferritin family and the Uma2-like family proteins (see discussion of these proteins below). Two desiccation-related proteins have been shown to be involved in desiccation but not radiation resistance ; Ro ribonucleoprotein is apparently involved in UV resistance , and topoisomerase IB, while active, has no known role in DR . So far, none of these genes has been linked experimentally with radioresistance in DR.
Identification of xenologous gene displacement by phylogenetic analysis
It is well-established that HGT has made major contributions to the gene repertoires of most thermophilic bacteria as supported by the presence of numerous genes with unexpectedly high similarity to and/or phylogenetic affinity with genes typical of hyperthermophilic archaea [56–60]. These cases include even those proteins that have orthologs in mesophilic bacteria but, as shown by phylogenetic analysis, have clear affinity to archaea or thermophilic bacteria from distant bacterial lineages, which is indicative of XGD . To investigate the impact of HGT from thermophiles leading to XGD on the evolution of the gene repertoire of TT, we determined the taxonomic affiliations of the proteins from the common gene core of TT and DR. We used the taxonomic distribution of best hits in BLAST searches for preliminary identification of HGT candidates, followed by a detailed phylogenetic analysis of selected genes.
Amino acid composition bias is known to affect not only sequence similarity searches but phylogeny reconstruction as well . We tested this effect on our data set by comparing the sequence-based maximum likelihood trees to the neighbor-joining trees reconstructed from the amino acid frequencies of corresponding proteins (see Additional file 1, "Influence of amino acid composition on phylogenetic reconstructions"). We found that, in the majority of cases (>80%), the topology of the sequence-based tree was not congruent with that of the amino acid composition tree; thus, the effect of the amino acid composition on the breakdown shown in Figure 5B is unlikely to be substantial. Taken together, these results suggest that XGD involving genes from thermophiles made a measurable contribution to the evolution of the core gene set of TT after the divergence from the common DR-TT ancestor; no such contribution was detected in the case of DR. While this interpretation seems most plausible given the ecological proximity of TT and other thermophiles, it cannot be ruled out that the observed patterns (Figure 5B) are partially explained by XGD with genes from mesophilic bacteria in the DR lineage. More generally, these results emphasize that the taxonomic distribution of best database hits can be taken only as a rough and preliminary indicator of HGT.
Examples of horizontally transferred genes in TT and DR
Selected horizontally transferred genes/operons
T. thermophilus Protein ID
D. radiodurans Protein ID
Xenologous gene displacement
Ribosomal protein L15 (COG0200)
Proteobacterial affinity of DR; Affinity to Thermatoga/Aquifex/other Gram-positive bacteria in TT
Ribosomal protein L30/L7E (COG1841)
Thermatogales/Aquifecales version in TT, affinity to other bacteria in DR
Thiamine biosynthesis protein ThiC (COG0422)
Proteobacterial affinity in DR; Cyanobacteria/Actinobacterial affinity in TT
Threonine synthase (COG0498)
Proteobacterial version in DR; Cyanobacteria/Actinobacterial affinity in TT
Trk system potassium uptake protein trkG, trkA (COG0168, COG0569)
DR1667, DR1668; DR1666
Archaeal version in TT; Gram-positive version in DR
Ribonucleoside-diphosphate reductase alpha chain and beta chains (COG0208, COG0209)
Proteobacterial affinity for TT; affinity with Gram-positive bacteria in DR
Paralogous genes acquisition
5-formyltetrahydrofolate cyclo-ligase (COG0212)
Gram-positive version in DR and TT; additional pseudoparalog TT1803 appears to be of archaeal origin
Non-orthologous gene displacement/acquisition
Represented by two non-homologous genes in DR (protobacterial version) and TT (orthologs in Treponema, Wolbachia)
Purine-nucleoside phosphorylase (purine degradation enzyme)
TTC1070 (COG0813) TTC0194 (COG0005)
TT has two non-homologous enzymes; TTC0194 probably was acquired from a thermopilic source
Expanded families of paralogs
Paralogous gene families expanded in DR
Number of representatives in DR
Number of representatives in TT
Widespread families expanded in DR
MutT-like phosphohydrolases (Nudix)
COG1768, COG1408, COG1692, COG0639
1; 1; 1; 9
1; 0; 1; 2
Lipase-like alpha/beta hydrolase
PadR-like transcriptional regulators (possibly involved in chemical stress response)
MOSC sulfur-carrier domains
FlaR like kinases
LigT phosphatases (may participate in RNA repair or methabolism)
TerZ family (could confer resistance to a variety of DNA-damaging agents)
PR1 family (stress response)
DinB family (DNA damage and stress inducible proteins)
COG2318; no COG
Unique DR families
GRXGG repeats containing protein
DR0082, DR2593, DR1748
Alpha/beta proteins, tryptophan-rich
Proteins with GXTXXXG and CXPXXXC motifs (DR0871 has duplication of the domain)
DR0871, DR1920, DR2360
Secreted alpha/beta proteins with a single conserved domain
DR1251, DR1319, DR1545
Predominantly alpha-helical proteins
DR0481, DR1195, DR1301
Predominantly alpha-helical proteins
Predicted metabolic regulator containing V4R domain
Predicted sirohydrochlorin cobaltochelatase
Conserved histidine rich protein (now also found in Caulobacter and Mesorhizobium)
Paralogous gene families expanded in TT
Number of representatives in DR
Number of representatives in TT
Uncharacterized protein conserved in cyanobacteria, Uma2 homolog
ABC-type sugar transport system, periplasmic component
ABC-type sugar transport system, permease component
ABC-type sugar transport systems, ATPase components
ABC-type sugar transport systems, permease components
ABC-type Fe3+ transport system, permease component
Fe2+/Zn2+ uptake regulation proteins
Minimal nucleotidyl transferases
COG1487, COG3744, COG4113, COG4374 COG1848
Antitoxin of toxin-antitoxin stability system
Nucleotide-binding proteins of the UspA family
TRAP-type mannitol/chloroaromatic compound transport system, periplasmic component
TRAP-type mannitol/chloroaromatic compound transport system, small permease component
TRAP-type mannitol/chloroaromatic compound transport system, large permease component
Arabinose efflux permease
Predicted phosphoesterases, related to the Icc protein
HEPN, Nucleotide-binding domain
Tfp pilus assembly protein FimT
Tfp pilus assembly protein PilE
Notably, several protein families that are expanded in TT but are absent in DR belong to the set of potential thermophilic determinants (HEPN nucleotide-binding domain; predicted phosphoesterases related to the Icc protein)  or are expanded in thermophylic archaea (PIN-like nuclease domain, minimal nucleotidyl transferases, UspA-like nucleotide-binding domain) , Table 4). In particular, TT has three paralogs of the archaea-specific tungsten-containing aldehyde ferredoxin oxidoreductase (TTC0012, TTC1834, TTP0122, TTP0212), which is the first occurrenceof this enzyme in thermophilic bacteria. However, these enzymes are present in several mesophilic bacteria, and have various substrate specificities and might be involved in sugar, amino acids or sulfur metabolism [72, 73].
Comparison of DNA repair and stress response systems
Comparison of general repair pathways in DR and TT
DNA repair pathways
DR – direct damage reversal
BER – base excision repair
NER – nucleotide excision repair
mMM – methylation-dependent mismatch repair
MM – mismatch repair
MMY – mutY-dependent repair
VSP – very short path mismatch repair
RER – recombinational repair
MP – multiple pathways
Total number (Genome fraction)
However, a conclusion that there are no important differences between the repair systems of TT and DR might be premature. Recently, several additional proteins of DR have been implicated in DNA or RNA repair, either in direct experiments or on the basis of up-regulation following irradiation, complemented with protein sequence analysis. These putative repair enzymes include DRB0094, an RNA ligase  that is strongly up-regulated in response to irradiation  and might be involved in an uncharacterized RNA repair process; a predicted double-strand break repair complex specific for recovery after irradiation, which consists of DRB0100, a DNA ligase, DRB0098, a protein containing an HD family phosphatase and polynucleotide kinase domains ; DRB0099, a predicted phosphatase of the H2Macro superfamily ( and KSM, unpublished observations); a double-stranded DNA-binding protein PprA (DRA0346), which stimulates the DNA end-joining reaction in vitro ; a predicted DNA single-strand annealing protein DdrA (DR0423) ; a regulator of radiation response IrrE (DR0167); a metal-dependent protease fused to a helix-turn-helix domain [79, 80]; and the uncharacterized protein DR0070 that has been shown to be essential for full resistance to acute irradiation . Among these poorly characterized (predicted) repair proteins of DR, only DdrA has an ortholog in the TT genome. In general, these putative repair genes are sparsely represented in bacteria, and it appears most unlikely that they were present in the common ancestor of TT and DR; most likely, these genes were acquired by the DR lineage via HGT after the divergence from the common ancestor with TT, and might have contributed to the evolution of the resistance phenotype. However, functional relevance of these genes to radiation resistance remains to be confirmed because most of the corresponding knockout mutants showed only relatively small to moderate decreases in radiation resistance [27, 78].
Among the unique (predicted) repair enzymes of TT, the most conspicuous ones are the components of the putative thermophile-specific repair system, which are predominantly encoded on the TT megaplasmid (Figure 7A). The functional features of the proteins encoded in this system (COG1203, a DNA helicase, often fused to a predicted HD-superfamily hydrolase; COG1468, a RecB-family exonuclease; COG1353, predicted polymerase) suggest that they are involved in an as yet uncharacterized DNA repair pathway. It has been hypothesized that this novel gene complex might be functionally analogous to the bacterial-eukaryotic system of translesion, mutagenic repair whose central components are DNA polymerases of the UmuC-DinB-Rad30-Rev1 superfamily, which typically are missing in thermophiles .
Comparison of proteins comprising various (predicted) systems involved in stress response reveals a greater number and diversity of such proteins in DR, which has 26 COGs with relevant functions that are not represented in TT compared to 3 such COGs in TT (see Additional file 1, Table 4S). Altogether, there are 147 proteins in DR in this category and 86 in TT, suggesting that some of them are additionally expanded in DR (see the section on "Expanded families").
Enzymatic systems of defense against oxidative stress predicted in TT and DR also show important differences. TT has one Mn-dependent superoxide dismutase (SodA)  and one Mn-dependent catalase (with no ortholog in DR), whereas DR encodes three superoxide dismutases (one of which is the ortholog of the SodA gene of TT) and three predicted catalases with no TT orthologs . Additionally, DR has a cytochrome C peroxidase (DRA0301) and a predicted iron-dependent peroxidase (DRA0145), enzymes that are likely to provide protection against toxic peroxides . Orthologs of these enzymes are rare among bacteria, suggesting that the Deinococcus lineage acquired them via HGT after the divergence of TT and DR from the common ancestor. Reduction of oxidized methionine residues in proteins is crucial for survival of cells under oxidative stress [82, 83]. Consistent with this idea, two peptide methionine sulfoxide reductases (PMSRs), MsrA (DR1849) and MsrB (DR1378), are encoded in the DR genome [84, 85], whereas none are present in TT. Interestingly, both PMSRs are also missing in Aquifex, Thermotoga and most thermophilic archaea, suggesting at least two possibilities: either this type of oxidative damage is rare at high temperatures or the known PMSRs are replaced by uncharacterized analogous enzymes due to the inefficiency of the former at high temperatures.
Oxidative stress defense mechanisms also might include control of Mn and Fe partitioning in the cell . Proteins of the Dps/ferritin family are required for the storage of iron in a non-reactive state, which prevents iron-catalyzed formation of hydroxyl radicals, thus protecting the cell from iron toxicity (Fenton-type chemistry) . Two Dps-related proteins are encoded in the DR genome (DR2263, DRB0092, COG1528), and it has been shown that one of them (DR2263) protects DNA from both hydroxyl radical cleavage and from DNase I-mediated cleavage . Some proteins homologous to DPS can non-specifically bind DNA and therefore are viewed as DNA-specific protectors . Like most thermophiles, TT has no proteins of this family but encodes a ferritin from another family (TTC1122).
Since desiccation also causes oxidative stress, proteins involved in desiccation resistance belong to the general cellular defense category . Desiccation-related proteins from at least two distinct families have been detected in DR . These desiccation resistance protein families (Lea 76 family, DR0105 and DR1172; and Lea14 family, DR1372) are not represented in TT. However, three TT proteins (TTP0170, TTP0166, TTP0169) are homologs of another desiccation resistance protein that was originally characterized in a plant, Craterostigma plantagineum ; DR also has two proteins of this family, DRB0118 and DRA0258. These proteins are distantly related to COG1633 and belong to the ferritin family of iron storage proteins (KSM, unpublished). Highly conserved homologs of these proteins are also present in thermophilic bacteria and archaea. Two desiccation-related proteins (DR1172 and DRB0118) appear to be essential for desiccation resistance but not for radiation resistance in DR .
Comparison of the genome partitions of TT and DR
Homology between the DR and TT megaplasmids
A. Number of orthologs of TT proteins, encoded in DR genome partitions
P(χ 2 )
7 × 10 -15
B. Number of orthologs of DR proteins encoded in TT genome partitions
P(χ 2 )
89 (87. 9)
3 × 10 -24
Notably, among the putative thermophilic determinants of TT, ~50% are encoded on the megaplasmid (18 out of 36). Of these, 11 belong to the putative mobile thermophile-specific DNA repair system [16, 45] (Figure 7A). Additionally, the megaplasmid carries at least four other genes associated with this system, which have not yet been assigned to COGs (Figure 7A). Moreover, the TT megaplasmid also carries a pseudogene for reverse gyrase, the most conspicuous signature protein of hyperthermophiles [15, 17].
Recently, the genome of another strain of TT (HB8) has been completely sequenced and became available in public databases . A preliminary comparison of the two TT strains (HB27 and HB8) revealed considerable differences in the gene orders and contents of the megaplasmids (see Additional file 1, Figure 3S). Interestingly, these differences in gene content are derived mostly from genes that appear to be associated with the thermophylic lifestyle. In particular, strain HB8 encodes an intact reverse gyrase. Thus, it appears most likely that the gene for reverse gyrase was acquired from a hyperthermophilic source by the TT lineage and was present in the common ancestor of HB27 and HB8 but decayed in the former. Conversely, HB27 encodes a unique, three-domain fusion protein (DnaQ endonuclease, DinG helicase and RecQ helicase), whereas HB8 lacks the DinG and RecQ orthologs (Figure 7B). Furthermore, there are unexpected differences in the organization of predicted thermophile-specific repair systems between the two strains of TT. Specifically, HB27 contains a "gram-positive version" (TTP0132-TTP0136), whereas HB8 has a "proteobacterial version" of these genes (TTHB186-TTHB194) (Figure 7A).
Furthermore, a nearly complete draft genome sequence of Deinococcus geothermalis (DG) has recently become publicly available . Since DG is closely related to DR but is moderately thermophilic, we searched for "thermophilic" genes in DG genome. Using the "thermophilic" protein sequences ( and see above) of TT and DR as queries, we identified orthologs of 5 of the 6 DR proteins from this set (four of these are also present in both TT strains) and orthologs of 7 of the remaining 23 "thermophilic" proteins of TT (all components of the predicted thermophilic DNA repair system). These 7 proteins had the respective TT proteins as the best hits, and their monophyly was supported by phylogenetic analysis (data not shown), suggesting that these genes were already present in the genome of the DR-TT common ancestor.
These observations give rise to two hypotheses: (i) the TT megaplasmid is essential for the survival of the organism at high temperatures. Consistent with this idea, we detected an expansion of a two-component toxin-antitoxin system, which consists of a PIN-like nuclease (toxin) and a MazE family transcriptional regulator (antitoxin) . Such a system is known to be responsible for the segregational stability of antibiotic resistance plasmids and other plasmids via selective elimination of cells that have failed to acquire a plasmid copy , and/or exclusion of competing plasmids ; and (ii) the TT megaplasmid is a dynamic genome compartment and a veritable sink for horizontally transferred genes, some of which might affect the thermophilic phenotype of this bacterium. This is compatible with the considerable differences in gene content observed between the two TT strains.
Specific roles of plasmid-borne genes in recovery from DNA damage have been proposed previously, including class Ib ribonucleotide reductase, periplasmic alkaline phosphatase, and extracellular nuclease, and subsequent analyses revealed at least five other genes implicated in this process (two operons, DRB0098-DRB0100, DRB0094-DRB0084, see above in "Comparison of DNA repair and stress response systems") [27, 76]. There is also a toxin-antitoxin system operon in the DR177 megaplasmid (DRB0012a and another gene located immediately upstream of DRB0012a and currently absent from the genome annotation). Additionally, there are five other toxin-antitoxin systems encoded on DR412, the smaller chromosome of DR, which might be related to maintenance of DR177 megaplasmid in DR. The DR412 chromosome appears to have some special features as well. Numerous genes that apparently have been acquired by DR via HGT after the divergence from the common ancestor with TT and are implicated in various processes of amino acid and nucleotide degradation, map to this genome partition. Thus, the megaplasmids of TT and DR (and the smaller chromosome of DR, DR412) appear to have participated in extensive HGT, which might have been important for the evolution of thermophily and radioresistance, although the repertoires of the respective acquired genes are completely different.
TT and DR share a large core of genes and form a clade in the gene-content tree, which supports the idea that these bacteria form a distinct clade, as indicated previously by phylogenetic analysis of rRNA and various proteins, and that the evolution of their gene complements was, roughly, clock-like. However, major differences between the gene repertoires of TT and DR were observed, indicating that both genomes lost numerous ancestral genes and acquired distinct sets of new genes primarily via HGT. In addition, numerous lineage-specific expansions of paralogous gene families were identified, particularly, in DR.
Some of the differences in the gene repertoires of TT and DR can be linked to the distinctive adaptive strategies of these bacteria. For example, TT appears to have acquired many genes from (hyper)thermophilic bacteria and archaea, whereas DR apparently acquired various genes involved in oxidative stress response and other "house-cleaning" functions from diverse bacterial sources.
The gene content of the TT megaplasmid (pTT27) and the DR megaplasmid (DR177) are sufficiently similar to conclude that they evolved from a common ancestor. To our knowledge, this is the first evidence of persistence of a megaplasmid beyond the genus level. However, the TT megaplasmid also carries many genes whose functions are implicated in the thermophylic phenotype; in particular, components of the predicted thermophile-specific repair system. These megaplasmids are likely to be essential for the survival of both TT and DR, with their maintenance controlled by toxin-antitoxin systems. Furthermore, the substantial differences between the gene repertoires of the megaplasmids of TT strains HB27 and HB8 indicate that this genome partition has been highly dynamic, with high rates of gene loss and HGT events occurring during evolution.
The evolutionary reconstruction based on the parsimony principle is generally compatible with the idea that the common ancestor of TT and DR was a mesophilic bacterium, whereas the thermophylic phenotype of TT evolved gradually via HGT of genes from thermophiles. Conversely, the radiation-desiccation resistance phenotype of DR might have gradually evolved via HGT of genes from other mesophiles, particularly, with highly developed oxidative stress response systems. However, it should be noted that the TT-DR common gene core includes dozens of genes of apparent archaeal origin or, at least, genes with thermophilic affiliation. Moreover, the DG genome encodes a few additional "thermophilic determinants" that are missing in DR but are unlikely to have been transferred from a thermophilic source independently of TT, as shown by comparative-genomic and phylogenetic analyses described here. Thus, acquisition of a considerable number of archaeal genes might have occurred along the evolutionary branch leading to the common ancestor of TT and DR. Accordingly, at this stage, we cannot rule out the possibility that this ancestor was a moderate thermophile rather than a mesophile. Further sequencing of bacterial genomes of the Thermus-Deinococcus clade should allow more definitive comparative-genomic analysis to elucidate the nature of the common ancestor of these bacteria.
Irradiation and desiccation
Irradiations. Three TT colony-isolates were inoculated individually in liquid TGY (10 g/L Bactotryptone, 1 g/L glucose, 5 g/L yeast extract) and incubated at 70°C. Cells were harvested at OD600 ~0.9, which corresponds to 107 – 108 colony forming units (CFU)/ml; 1 TT cell/CFU. TT cells grown in TGY were examined for their total Mn and Fe content by ICP-MS (see main text). For radiation resistance assays, cells were irradiated without change of broth on ice with 60Co at 6.8 kGy/hour (60Co Gammacell irradiation unit [J. L. Shepard and Associates, Model 109]). At the indicated doses, cultures (3 biological replicates) were appropriately diluted and plated on solid medium (8 g/L Bactotryptone, 4 g/L yeast extract, 3 g/L NaCl, pH 7.3, 2.8% Bactoagar), and CFU counts were determined after 2 days' incubation at 65°C.
Desiccation. Five separate colony-isolates were pre-grown in TGY as for irradiation trials. Cell samples 106–107 cells (25 μl) were transferred to microtiter plates, which were transferred to desiccation chambers containing anhydrous calcium sulfate (drierite) and incubated at room temperature or 65°C. At the indicated times, cells were re-suspended in TGY, and CFU-survival frequencies were determined by dilution-plating on solid medium (65°C).
The sets of predicted proteins of TT and DR were searched against each other for symmetrical and non-symmetrical hits using PSI-BLAST with expectation (E) value threshold of 10-5. Taxonomic affiliation was determined by best hits in non-redundant database of protein sequences at the National Center of Biotechnology Information (NIH, Bethesda) using BLASTP program  with default expectation value (0.01). Assignments to COGs were performed using the COGNITOR program  and CDD-search against COG-based profiles . Contradictory assignments were resolved manually. Lineage-specific expansions (LSE) were identified as described previously [66, 97]. The common genomic core was determined as follows: among genes that were not assigned to any COG, orthology relationships between TT and DR were determined via symmetrical best hits. Genes belonging to COGs, shared between TT and DR and having only one ortholog from each of the two genomes were directly assigned to the core. For multiple-paralog COGs, symmetrical best hits between GOG members were used to refine the relationships between TT and DR proteins. Members of the corresponding lineage-specific expansions were added to SymBeT pairs to form many-to-many core clusters. LBCA gene set was determined using an empirical parsimony procedure based on COG phyletic patterns (See Additional file 1: "Reconstruction of the gene set of the Last Bacterial Common Ancestor") which assigned a COG to LBCA if it was present in several diverse bacterial clades. All COGs that were present in LBCA and in TT and/or DR were assigned to the DR-TT ancestor.
Multiple alignments for phylogenetic analysis were constructed using the MUSCLE program ; columns containing gaps in >30% of the sequences were discarded. Maximum likelihood trees were constructed using the ProtML program of the MOLPHY package by optimizing the least-square trees with local rearrangements . Trees based on amino acid content were constructed from the matrix of Euclidean distances between frequency vectors using the NEIGHBOR program of the PHYLIP package . Support for particular arrangements of species (relationships between DR, TT, thermophiles and mesophiles) was calculated using the bipartition analysis of bootstrap samples from original sequences (see Additional file 1, "Influence of amino acid composition on phylogenetic reconstructions"). The gene-content tree based on COG patterns was constructed using the NEIGHBOR program of the PHYLIP package ; the number of COGs shared between two genomes was normalized by the smaller genome size .
This research was partially supported by the Office of Science (Biological and Environmental Research), U. S. Department of Energy (DOE) grant number DE-FG02-04ER63918 awarded to MJD. We are grateful to James K. Fredrickson and Heather M. Kostandarithes at Pacific Northwest National Laboratory for Fe and Mn ICP-MS analyses of T. thermophilus.
- Weisburg WG, Giovannoni SJ, Woese CR: The Deinococcus-Thermus phylum and the effect of rRNA composition on phylogenetic tree construction. Syst Appl Microbiol. 1989, 11: 128-134.View ArticlePubMedGoogle Scholar
- Eisen JA: The RecA protein as a model molecule for molecular systematic studies of bacteria: comparison of trees of RecAs and 16S rRNAs from the same species. J Mol Evol. 1995, 41 (6): 1105-1123. 10.1007/BF00173192.PubMed CentralView ArticlePubMedGoogle Scholar
- Rainey FA, Nobre MF, Schumann P, Stackebrandt E, da Costa MS: Phylogenetic diversity of the deinococci as determined by 16S ribosomal DNA sequence comparison. Int J Syst Bacteriol. 1997, 47 (2): 510-514.View ArticlePubMedGoogle Scholar
- Gupta RS, Bustard K, Falah M, Singh D: Sequencing of heat shock protein 70 (DnaK) homologs from Deinococcus proteolyticus and Thermomicrobium roseum and their integration in a protein-based phylogeny of prokaryotes. J Bacteriol. 1997, 179 (2): 345-357.PubMed CentralPubMedGoogle Scholar
- Gupta RS: Protein phylogenies and signature sequences: A reappraisal of evolutionary relationships among archaebacteria, eubacteria, and eukaryotes. Microbiol Mol Biol Rev. 1998, 62 (4): 1435-1491.PubMed CentralPubMedGoogle Scholar
- Huang YP, Ito J: DNA polymerase C of the thermophilic bacterium Thermus aquaticus: classification and phylogenetic analysis of the family C DNA polymerases. J Mol Evol. 1999, 48 (6): 756-769.View ArticlePubMedGoogle Scholar
- Makarova KS, Aravind L, Wolf YI, Tatusov RL, Minton KW, Koonin EV, Daly MJ: Genome of the extremely radiation-resistant bacterium Deinococcus radiodurans viewed from the perspective of comparative genomics. Microbiol Mol Biol Rev. 2001, 65 (1): 44-79. 10.1128/MMBR.65.1.44-79.2001.PubMed CentralView ArticlePubMedGoogle Scholar
- Battista JR, Earl AM, Park MJ: Why is Deinococcus radiodurans so resistant to ionizing radiation?. Trends Microbiol. 1999, 7 (9): 362-365. 10.1016/S0966-842X(99)01566-8.View ArticlePubMedGoogle Scholar
- Hirsch P, Gallikowski CA, Siebert J, Peissl K, Kroppenstedt R, Schumann P, Stackebrandt E, Anderson R: Deinococcus frigens sp. nov., Deinococcus saxicola sp. nov., and Deinococcus marmoris sp. nov., low temperature and draught-tolerating, UV-resistant bacteria from continental Antarctica. Syst Appl Microbiol. 2004, 27 (6): 636-645. 10.1078/0723202042370008.View ArticlePubMedGoogle Scholar
- Ferreira AC, Nobre MF, Rainey FA, Silva MT, Wait R, Burghardt J, Chung AP, da Costa MS: Deinococcus geothermalis sp. nov. and Deinococcus murrayi sp. nov., two extremely radiation-resistant and slightly thermophilic species from hot springs. Int J Syst Bacteriol. 1997, 47 (4): 939-947.View ArticlePubMedGoogle Scholar
- Holt JG, Krieg NR, Sneath PH, Staley JT, Williams ST: Bergy's Manual of Determinative Bacteriology. 1997, London , Williams@Wilkins, 9Google Scholar
- Henne A, Bruggemann H, Raasch C, Wiezer A, Hartsch T, Liesegang H, Johann A, Lienard T, Gohl O, Martinez-Arias R, Jacobi C, Starkuviene V, Schlenczeck S, Dencker S, Huber R, Klenk HP, Kramer W, Merkl R, Gottschalk G, Fritz HJ: The genome sequence of the extreme thermophile Thermus thermophilus. Nat Biotechnol. 2004, 22 (5): 547-553. 10.1038/nbt956.View ArticlePubMedGoogle Scholar
- NCBI, Thermus thermophilus HB8, complete genome. [http://www.ncbi.nlm.nih.gov/genomes/framik.cgi?db=Genome&gi=530]
- White O, Eisen JA, Heidelberg JF, Hickey EK, Peterson JD, Dodson RJ, Haft DH, Gwinn ML, Nelson WC, Richardson DL, Moffat KS, Qin H, Jiang L, Pamphile W, Crosby M, Shen M, Vamathevan JJ, Lam P, McDonald L, Utterback T, Zalewski C, Makarova KS, Aravind L, Daly MJ, Fraser CM: Genome sequence of the radioresistant bacterium Deinococcus radiodurans R1. Science. 1999, 286 (5444): 1571-1577. 10.1126/science.286.5444.1571.PubMed CentralView ArticlePubMedGoogle Scholar
- Forterre P: A hot story from comparative genomics: reverse gyrase is the only hyperthermophile-specific protein. Trends Genet. 2002, 18 (5): 236-237. 10.1016/S0168-9525(02)02650-1.View ArticlePubMedGoogle Scholar
- Makarova KS, Aravind L, Grishin NV, Rogozin IB, Koonin EV: A DNA repair system specific for thermophilic Archaea and bacteria predicted by genomic context analysis. Nucleic Acids Res. 2002, 30 (2): 482-496. 10.1093/nar/30.2.482.PubMed CentralView ArticlePubMedGoogle Scholar
- Makarova KS, Wolf YI, Koonin EV: Potential genomic determinants of hyperthermophily. Trends Genet. 2003, 19 (4): 172-176. 10.1016/S0168-9525(03)00047-7.View ArticlePubMedGoogle Scholar
- McDonald JH, Grasso AM, Rejto LK: Patterns of temperature adaptation in proteins from Methanococcus and Bacillus. Mol Biol Evol. 1999, 16 (12): 1785-1790.View ArticlePubMedGoogle Scholar
- Haney PJ, Badger JH, Buldak GL, Reich CI, Woese CR, Olsen GJ: Thermal adaptation analyzed by comparison of protein sequences from mesophilic and extremely thermophilic Methanococcus species. Proc Natl Acad Sci U S A. 1999, 96 (7): 3578-3583. 10.1073/pnas.96.7.3578.PubMed CentralView ArticlePubMedGoogle Scholar
- Thompson MJ, Eisenberg D: Transproteomic evidence of a loop-deletion mechanism for enhancing protein thermostability. J Mol Biol. 1999, 290 (2): 595-604. 10.1006/jmbi.1999.2889.View ArticlePubMedGoogle Scholar
- Sterner R, Liebl W: Thermophilic adaptation of proteins. Crit Rev Biochem Mol Biol. 2001, 36 (1): 39-106. 10.1080/20014091074174.View ArticlePubMedGoogle Scholar
- Chakravarty S, Varadarajan R: Elucidation of factors responsible for enhanced thermal stability of proteins: a structural genomics based study. Biochemistry. 2002, 41 (25): 8152-8161. 10.1021/bi025523t.View ArticlePubMedGoogle Scholar
- Vieille C, Zeikus GJ: Hyperthermophilic enzymes: sources, uses, and molecular mechanisms for thermostability. Microbiol Mol Biol Rev. 2001, 65 (1): 1-43. 10.1128/MMBR.65.1.1-43.2001.PubMed CentralView ArticlePubMedGoogle Scholar
- Mallick P, Boutz DR, Eisenberg D, Yeates TO: Genomic evidence that the intracellular proteins of archaeal microbes contain disulfide bonds. Proc Natl Acad Sci U S A. 2002, 99 (15): 9679-9684. 10.1073/pnas.142310499.PubMed CentralView ArticlePubMedGoogle Scholar
- Hickey DA, Singer GA: Genomic and proteomic adaptations to growth at high temperature. Genome Biol. 2004, 5 (10): 117-10.1186/gb-2004-5-10-117.PubMed CentralView ArticlePubMedGoogle Scholar
- Tanaka M, Earl AM, Howell HA, Park MJ, Eisen JA, Peterson SN, Battista JR: Analysis of Deinococcus radiodurans's transcriptional response to ionizing radiation and desiccation reveals novel proteins that contribute to extreme radioresistance. Genetics. 2004, 168 (1): 21-33. 10.1534/genetics.104.029249.PubMed CentralView ArticlePubMedGoogle Scholar
- Liu Y, Zhou J, Omelchenko MV, Beliaev AS, Venkateswaran A, Stair J, Wu L, Thompson DK, Xu D, Rogozin IB, Gaidamakova EK, Zhai M, Makarova KS, Koonin EV, Daly MJ: Transcriptome dynamics of Deinococcus radiodurans recovering from ionizing radiation. Proc Natl Acad Sci U S A. 2003, 100 (7): 4191-4196. 10.1073/pnas.0630387100.PubMed CentralView ArticlePubMedGoogle Scholar
- Ghosal D, Omelchenko MV, Gaidamakova EK, Matrosova VY, Vasilenko A, Venkateswaran A, Zhai M, Kostandarithes HM, Brim H, Makarova K, Wackett LP, Fredrickson JK, Daly MJ: How radiation kills cells: Survival of Deinococcus radiodurans and Shewanella oneidensis under oxidative stress. FEMS Microbiol Lett. 2005,Google Scholar
- Daly MJ, Gaidamakova EK, Matrosova VY, Vasilenko A, Zhai M, Venkateswaran A, Hess M, Omelchenko MV, Kostandarithes HM, Makarova KS, Wackett LP, Fredrickson JK, Ghosal D: Accumulation of Mn(II) in Deinococcus radiodurans facilitates gamma-radiation resistance. Science. 2004, 306 (5698): 1025-1028. 10.1126/science.1103185.View ArticlePubMedGoogle Scholar
- Battista JR: Against all odds: the survival strategies of Deinococcus radiodurans. Annu Rev Microbiol. 1997, 51: 203-224. 10.1146/annurev.micro.51.1.203.View ArticlePubMedGoogle Scholar
- Archibald FS, Fridovich I: . Arch Biochem Biophys. 1982, 214: 452-10.1016/0003-9861(82)90049-2.View ArticlePubMedGoogle Scholar
- Lee Y, Lee C, Yoon J: High temperature dependence of 2,4-dichlorophenoxyacetic acid degradation by Fe3+/H(2)O(2) system. Chemosphere. 2003, 51 (9): 963-971. 10.1016/S0045-6535(03)00043-2.View ArticlePubMedGoogle Scholar
- Snel B, Bork P, Huynen MA: Genome phylogeny based on gene content. Nat Genet. 1999, 21 (1): 108-110. 10.1038/5052.View ArticlePubMedGoogle Scholar
- Wolf YI, Rogozin IB, Grishin NV, Koonin EV: Genome trees and the tree of life. Trends Genet. 2002, 18 (9): 472-479. 10.1016/S0168-9525(02)02744-0.View ArticlePubMedGoogle Scholar
- Wolf YI, Rogozin IB, Grishin NV, Tatusov RL, Koonin EV: Genome trees constructed using five different approaches suggest new major bacterial clades. BMC Evol Biol. 2001, 1 (1): 8-10.1186/1471-2148-1-8.PubMed CentralView ArticlePubMedGoogle Scholar
- Korbel JO, Snel B, Huynen MA, Bork P: SHOT: a web server for the construction of genome phylogenies. Trends Genet. 2002, 18 (3): 158-162. 10.1016/S0168-9525(01)02597-5.View ArticlePubMedGoogle Scholar
- Tatusov RL, Fedorova ND, Jackson JD, Jacobs AR, Kiryutin B, Koonin EV, Krylov DM, Mazumder R, Mekhedov SL, Nikolskaya AN, Rao BS, Smirnov S, Sverdlov AV, Vasudevan S, Wolf YI, Yin JJ, Natale DA: The COG database: an updated version includes eukaryotes. BMC Bioinformatics. 2003, 4 (1): 41-10.1186/1471-2105-4-41.PubMed CentralView ArticlePubMedGoogle Scholar
- Phillips RW, Wiegel J, Berry CJ, Fliermans C, Peacock AD, White DC, Shimkets LJ: Kineococcus radiotolerans sp. nov., a radiation-resistant, gram-positive bacterium. Int J Syst Evol Microbiol. 2002, 52 (Pt 3): 933-938. 10.1099/ijs.0.02029-0.PubMedGoogle Scholar
- Ferreira AC, Nobre MF, Moore E, Rainey FA, Battista JR, da Costa MS: Characterization and radiation resistance of new isolates of Rubrobacter radiotolerans and Rubrobacter xylanophilus. Extremophiles. 1999, 3 (4): 235-238. 10.1007/s007920050121.View ArticlePubMedGoogle Scholar
- Billi D, Friedmann EI, Hofer KG, Caiola MG, Ocampo-Friedmann R: Ionizing-radiation resistance in the desiccation-tolerant cyanobacterium Chroococcidiopsis. Appl Environ Microbiol. 2000, 66 (4): 1489-1492. 10.1128/AEM.66.4.1489-1492.2000.PubMed CentralView ArticlePubMedGoogle Scholar
- Mirkin BG, Fenner TI, Galperin MY, Koonin EV: Algorithms for computing parsimonious evolutionary scenarios for genome evolution, the last universal common ancestor and dominance of horizontal gene transfer in the evolution of prokaryotes. BMC Evol Biol. 2003, 3 (1): 2-10.1186/1471-2148-3-2.PubMed CentralView ArticlePubMedGoogle Scholar
- Olendzenski L, Liu L, Zhaxybayeva O, Murphey R, Shin DG, Gogarten JP: Horizontal transfer of archaeal genes into the deinococcaceae: detection by molecular and computer-based approaches. J Mol Evol. 2000, 51 (6): 587-599.View ArticlePubMedGoogle Scholar
- Levin I, Giladi M, Altman-Price N, Ortenberg R, Mevarech M: An alternative pathway for reduced folate biosynthesis in bacteria and halophilic archaea. Mol Microbiol. 2004, 54 (5): 1307-1318. 10.1111/j.1365-2958.2004.04339.x.View ArticlePubMedGoogle Scholar
- Mittenhuber G: Phylogenetic analyses and comparative genomics of vitamin B6 (pyridoxine) and pyridoxal phosphate biosynthesis pathways. J Mol Microbiol Biotechnol. 2001, 3 (1): 1-20.PubMedGoogle Scholar
- Jansen R, Embden JD, Gaastra W, Schouls LM: Identification of genes that are associated with DNA repeats in prokaryotes. Mol Microbiol. 2002, 43 (6): 1565-1575. 10.1046/j.1365-2958.2002.02839.x.View ArticlePubMedGoogle Scholar
- Sauer U, Canonaco F, Heri S, Perrenoud A, Fischer E: The soluble and membrane-bound transhydrogenases UdhA and PntAB have divergent functions in NADPH metabolism of Escherichia coli. J Biol Chem. 2004, 279 (8): 6613-6619. 10.1074/jbc.M311657200.View ArticlePubMedGoogle Scholar
- Venkateswaran A, McFarlan SC, Ghosal D, Minton KW, Vasilenko A, Makarova K, Wackett LP, Daly MJ: Physiologic determinants of radiation resistance in Deinococcus radiodurans. Appl Environ Microbiol. 2000, 66 (6): 2620-2626. 10.1128/AEM.66.6.2620-2626.2000.PubMed CentralView ArticlePubMedGoogle Scholar
- Ferry JG: CO dehydrogenase. Annu Rev Microbiol. 1995, 49: 305-333. 10.1146/annurev.mi.49.100195.001513.View ArticlePubMedGoogle Scholar
- Imlay JA: Pathways of oxidative damage. Annu Rev Microbiol. 2003, 57: 395-418. 10.1146/annurev.micro.57.030502.090938.View ArticlePubMedGoogle Scholar
- Zhulin IB, Nikolskaya AN, Galperin MY: Common extracellular sensory domains in transmembrane receptors for diverse signal transduction pathways in bacteria and archaea. J Bacteriol. 2003, 185 (1): 285-294. 10.1128/JB.185.1.285-294.2003.PubMed CentralView ArticlePubMedGoogle Scholar
- Melkonyan HS, Chang WC, Shapiro JP, Mahadevappa M, Fitzpatrick PA, Kiefer MC, Tomei LD, Umansky SR: SARPs: a family of secreted apoptosis-related proteins. Proc Natl Acad Sci U S A. 1997, 94 (25): 13636-13641. 10.1073/pnas.94.25.13636.PubMed CentralView ArticlePubMedGoogle Scholar
- Davis SJ, Vener AV, Vierstra RD: Bacteriophytochromes: phytochrome-like photoreceptors from nonphotosynthetic eubacteria. Science. 1999, 286 (5449): 2517-2520. 10.1126/science.286.5449.2517.View ArticlePubMedGoogle Scholar
- Battista JR, Park MJ, McLemore AE: Inactivation of two homologues of proteins presumed to be involved in the desiccation tolerance of plants sensitizes Deinococcus radiodurans R1 to desiccation. Cryobiology. 2001, 43 (2): 133-139. 10.1006/cryo.2001.2357.View ArticlePubMedGoogle Scholar
- Chen X, Quinn AM, Wolin SL: Ro ribonucleoproteins contribute to the resistance of Deinococcus radiodurans to ultraviolet irradiation. Genes Dev. 2000, 14 (7): 777-782.PubMed CentralPubMedGoogle Scholar
- Krogh BO, Shuman S: A poxvirus-like type IB topoisomerase family in bacteria. Proc Natl Acad Sci U S A. 2002, 99 (4): 1853-1858. 10.1073/pnas.032613199.PubMed CentralView ArticlePubMedGoogle Scholar
- Koonin EV, Makarova KS: Archaeal genomics: how much have we learned in six years and what's next?. Genome Biol. 2003, in press:Google Scholar
- Koonin EV, Makarova KS, Aravind L.: Horizontal gene transfer in prokaryotes - quantification and classification. AnnuRevMicrobiol. 2001Google Scholar
- Nelson KE, Clayton RA, Gill SR, Gwinn ML, Dodson RJ, Haft DH, Hickey EK, Peterson JD, Nelson WC, Ketchum KA, McDonald L, Utterback TR, Malek JA, Linher KD, Garrett MM, Stewart AM, Cotton MD, Pratt MS, Phillips CA, Richardson D, Heidelberg J, Sutton GG, Fleischmann RD, Eisen JA, Fraser CM: Evidence for lateral gene transfer between Archaea and bacteria from genome sequence of Thermotoga maritima. Nature. 1999, 399 (6734): 323-329. 10.1038/20601.View ArticlePubMedGoogle Scholar
- Aravind L, Tatusov RL, Wolf YI, Walker DR, Koonin EV: Evidence for massive gene exchange between archaeal and bacterial hyperthermophiles. Trends Genet. 1998, 14 (11): 442-444. 10.1016/S0168-9525(98)01553-4.View ArticlePubMedGoogle Scholar
- Worning P, Jensen LJ, Nelson KE, Brunak S, Ussery DW: Structural analysis of DNA sequence: evidence for lateral gene transfer in Thermotoga maritima. Nucleic Acids Res. 2000, 28 (3): 706-709. 10.1093/nar/28.3.706.PubMed CentralView ArticlePubMedGoogle Scholar
- Koski LB, Golding GB: The closest BLAST hit is often not the nearest neighbor. J Mol Evol. 2001, 52 (6): 540-542.View ArticlePubMedGoogle Scholar
- Cambillau C, Claverie JM: Structural and genomic correlates of hyperthermostability. J Biol Chem. 2000, 275 (42): 32383-32386. 10.1074/jbc.C000497200.View ArticlePubMedGoogle Scholar
- McDonald JH: Patterns of temperature adaptation in proteins from the bacteria Deinococcus radiodurans and Thermus thermophilus. Mol Biol Evol. 2001, 18 (5): 741-749.View ArticlePubMedGoogle Scholar
- Felsenstein J: Inferring phylogenies. 2004, Sunderland, MA , Sinauer Associates, Inc.Google Scholar
- Omelchenko MV, Makarova KS, Wolf YI, Rogozin IB, Koonin EV: Evolution of mosaic operons by horizontal gene transfer and gene displacement in situ. Genome Biol. 2003, 4 (9): R55-10.1186/gb-2003-4-9-r55.PubMed CentralView ArticlePubMedGoogle Scholar
- Jordan IK, Makarova KS, Wolf YI, Koonin EV: Gene conversions in genes encoding outer-membrane proteins in H. pylori and C. pneumoniae. Trends Genet. 2001, 17 (1): 7-10. 10.1016/S0168-9525(00)02151-X.View ArticlePubMedGoogle Scholar
- Gury J, Barthelmebs L, Tran NP, Divies C, Cavin JF: Cloning, deletion, and characterization of PadR, the transcriptional repressor of the phenolic acid decarboxylase-encoding padA gene of Lactobacillus plantarum. Appl Environ Microbiol. 2004, 70 (4): 2146-2153. 10.1128/AEM.70.4.2146-2153.2004.PubMed CentralView ArticlePubMedGoogle Scholar
- Anantharaman V, Aravind L: MOSC domains: ancient, predicted sulfur-carrier domains, present in diverse metal-sulfur cluster biosynthesis proteins including Molybdenum cofactor sulfurases. FEMS Microbiol Lett. 2002, 207 (1): 55-61. 10.1016/S0378-1097(01)00515-8.PubMedGoogle Scholar
- Anantharaman V, Aravind L: New connections in the prokaryotic toxin-antitoxin network: relationship with the eukaryotic nonsense-mediated RNA decay system. Genome Biol. 2003, 4 (12): R81-10.1186/gb-2003-4-12-r81.PubMed CentralView ArticlePubMedGoogle Scholar
- Oshima T: [Comparative studies on biochemical properties of an extreme thermophile, Thermus thermophilus HB 8 (author's transl)]. Seikagaku. 1974, 46 (10): 887-907.PubMedGoogle Scholar
- Makarova KS, Aravind L, Galperin MY, Grishin NV, Tatusov RL, Wolf YI, Koonin EV: Comparative genomics of the Archaea (Euryarchaeota): evolution of conserved protein families, the stable core, and the variable shell. Genome Res. 1999, 9 (7): 608-628.PubMedGoogle Scholar
- Roy R, Menon AL, Adams MW: Aldehyde oxidoreductases from Pyrococcus furiosus. Methods Enzymol. 2001, 331: 132-144.View ArticlePubMedGoogle Scholar
- Roy R, Adams MW: Characterization of a fourth tungsten-containing enzyme from the hyperthermophilic archaeon Pyrococcus furiosus. J Bacteriol. 2002, 184 (24): 6952-6956. 10.1128/JB.184.24.6952-6956.2002.PubMed CentralView ArticlePubMedGoogle Scholar
- Xu W, Shen J, Dunn CA, Desai S, Bessman MJ: The Nudix hydrolases of Deinococcus radiodurans. Mol Microbiol. 2001, 39 (2): 286-290. 10.1046/j.1365-2958.2001.02267.x.View ArticlePubMedGoogle Scholar
- Narumi I, Satoh K, Kikuchi M, Funayama T, Yanagisawa T, Kobayashi Y, Watanabe H, Yamamoto K: The LexA protein from Deinococcus radiodurans is not involved in RecA induction following gamma irradiation. J Bacteriol. 2001, 183 (23): 6951-6956. 10.1128/JB.183.23.6951-6956.2001.PubMed CentralView ArticlePubMedGoogle Scholar
- Martins A, Shuman S: An RNA ligase from Deinococcus radiodurans. J Biol Chem. 2004, 279 (49): 50654-50661. 10.1074/jbc.M407657200.View ArticlePubMedGoogle Scholar
- Narumi I, Satoh K, Cui S, Funayama T, Kitayama S, Watanabe H: PprA: a novel protein from Deinococcus radiodurans that stimulates DNA ligation. Mol Microbiol. 2004, 54 (1): 278-285. 10.1111/j.1365-2958.2004.04272.x.View ArticlePubMedGoogle Scholar
- Harris DR, Tanaka M, Saveliev SV, Jolivet E, Earl AM, Cox MM, Battista JR: Preserving genome integrity: the DdrA protein of Deinococcus radiodurans R1. PLoS Biol. 2004, 2 (10): e304-10.1371/journal.pbio.0020304.PubMed CentralView ArticlePubMedGoogle Scholar
- Earl AM, Mohundro MM, Mian IS, Battista JR: The IrrE protein of Deinococcus radiodurans R1 is a novel regulator of recA expression. J Bacteriol. 2002, 184 (22): 6216-6224. 10.1128/JB.184.22.6216-6224.2002.PubMed CentralView ArticlePubMedGoogle Scholar
- Hua Y, Narumi I, Gao G, Tian B, Satoh K, Kitayama S, Shen B: PprI: a general switch responsible for extreme radioresistance of Deinococcus radiodurans. Biochem Biophys Res Commun. 2003, 306 (2): 354-360. 10.1016/S0006-291X(03)00965-3.View ArticlePubMedGoogle Scholar
- Ludwig ML, Metzger AL, Pattridge KA, Stallings WC: Manganese superoxide dismutase from Thermus thermophilus. A structural model refined at 1.8 A resolution. J Mol Biol. 1991, 219 (2): 335-358. 10.1016/0022-2836(91)90569-R.View ArticlePubMedGoogle Scholar
- Moskovitz J: Methionine sulfoxide reductases: ubiquitous enzymes involved in antioxidant defense, protein regulation, and prevention of aging-associated diseases. Biochim Biophys Acta. 2005, 1703 (2): 213-219.View ArticlePubMedGoogle Scholar
- Morgan PE, Dean RT, Davies MJ: Protective mechanisms against peptide and protein peroxides generated by singlet oxygen. Free Radic Biol Med. 2004, 36 (4): 484-496. 10.1016/j.freeradbiomed.2003.11.021.View ArticlePubMedGoogle Scholar
- Moskovitz J, Weissbach H, Brot N: Cloning the expression of a mammalian gene involved in the reduction of methionine sulfoxide residues in proteins. Proc Natl Acad Sci U S A. 1996, 93 (5): 2095-2099. 10.1073/pnas.93.5.2095.PubMed CentralView ArticlePubMedGoogle Scholar
- Grimaud R, Ezraty B, Mitchell JK, Lafitte D, Briand C, Derrick PJ, Barras F: Repair of oxidized proteins. Identification of a new methionine sulfoxide reductase. J Biol Chem. 2001, 276 (52): 48915-48920. 10.1074/jbc.M105509200.View ArticlePubMedGoogle Scholar
- van Vliet AH, Ketley JM, Park SF, Penn CW: The role of iron in Campylobacter gene regulation, metabolism and oxidative stress defense. FEMS Microbiol Rev. 2002, 26 (2): 173-186. 10.1016/S0168-6445(02)00095-5.View ArticlePubMedGoogle Scholar
- Grove A, Wilkinson SP: Differential DNA Binding and Protection by Dimeric and Dodecameric forms of the Ferritin Homolog Dps from Deinococcus radiodurans. J Mol Biol. 2005, 347 (3): 495-508. 10.1016/j.jmb.2005.01.055.View ArticlePubMedGoogle Scholar
- Grant RA, Filman DJ, Finkel SE, Kolter R, Hogle JM: The crystal structure of Dps, a ferritin homolog that binds and protects DNA. Nat Struct Biol. 1998, 5 (4): 294-303. 10.1038/nsb0498-294.View ArticlePubMedGoogle Scholar
- Mattimore V, Battista JR: Radioresistance of Deinococcus radiodurans: functions necessary to survive ionizing radiation are also necessary to survive prolonged desiccation. J Bacteriol. 1996, 178 (3): 633-637.PubMed CentralPubMedGoogle Scholar
- Piatkowski D, Schneider K, Salamini F, Bartels D: Characterization of five abscisic acid-responsive cDNA clones isolated from the desiccation-tolerant plant Craterostigma plantagineum and their relationship to other water-stress genes. Plant Physiol. 1990, 94 (1682-88):Google Scholar
- JGI, Deinococcus geothermalis draft genome. [http://genome.jgi-psf.org/draft_microbes/deige/deige.home.html]. [http://genome.jgi-psf.org/draft_microbes/deige/deige.home.html]
- Grady R, Hayes F: Axe-Txe, a broad-spectrum proteic toxin-antitoxin system specified by a multidrug-resistant, clinical isolate of Enterococcus faecium. Mol Microbiol. 2003, 47 (5): 1419-1432. 10.1046/j.1365-2958.2003.03387.x.View ArticlePubMedGoogle Scholar
- Cooper TF, Heinemann JA: Postsegregational killing does not increase plasmid stability but acts to mediate the exclusion of competing plasmids. Proc Natl Acad Sci U S A. 2000, 97 (23): 12643-12648. 10.1073/pnas.220077897.PubMed CentralView ArticlePubMedGoogle Scholar
- Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997, 25 (17): 3389-3402. 10.1093/nar/25.17.3389.PubMed CentralView ArticlePubMedGoogle Scholar
- Tatusov RL, Natale DA, Garkavtsev IV, Tatusova TA, Shankavaram UT, Rao BS, Kiryutin B, Galperin MY, Fedorova ND, Koonin EV: The COG database: new developments in phylogenetic classification of proteins from complete genomes. Nucleic Acids Res. 2001, 29 (1): 22-28. 10.1093/nar/29.1.22.PubMed CentralView ArticlePubMedGoogle Scholar
- Marchler-Bauer A, Anderson JB, Cherukuri PF, DeWeese-Scott C, Geer LY, Gwadz M, He S, Hurwitz DI, Jackson JD, Ke Z, Lanczycki CJ, Liebert CA, Liu C, Lu F, Marchler GH, Mullokandov M, Shoemaker BA, Simonyan V, Song JS, Thiessen PA, Yamashita RA, Yin JJ, Zhang D, Bryant SH: CDD: a Conserved Domain Database for protein classification. Nucleic Acids Res. 2005, 33 Database Issue: D192-6.Google Scholar
- Lespinet O, Wolf YI, Koonin EV, Aravind L: The role of lineage-specific gene family expansion in the evolution of eukaryotes. Genome Res. 2002, 12 (7): 1048-1059. 10.1101/gr.174302.PubMed CentralView ArticlePubMedGoogle Scholar
- Edgar RC: MUSCLE: a multiple sequence alignment method with reduced time and space complexity. BMC Bioinformatics. 2004, 5 (1): 113-10.1186/1471-2105-5-113.PubMed CentralView ArticlePubMedGoogle Scholar
- Hasegawa M, Kishino H, Saitou N: On the maximum likelihood method in molecular phylogenetics. J Mol Evol. 1991, 32 (5): 443-445. 10.1007/BF02101285.View ArticlePubMedGoogle Scholar
- Felsenstein J: Inferring phylogenies from protein sequences by parsimony, distance, and likelihood methods. Methods Enzymol. 1996, 266: 418-427.View ArticlePubMedGoogle Scholar
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.