Analysis of three genomes within the thermophilic bacterial species Caldanaerobacter subterraneus with a focus on carbon monoxide dehydrogenase evolution and hydrolase diversity.

BACKGROUND: The Caldanaerobacter subterraneus species includes thermophilic fermentative bacteria able to grow on carbohydrates substrates with acetate and L-alanine as the main products. In this study, comprehensive analysis of three genomes of C. subterraneus subspecies was carried in order to identify genes encoding key metabolic enzymes and to document the genomic basis for the evolution of these organisms. METHODS: Average nucleotide identity and in silico DNA relatedness were estimated for the studied C. subterraneus genomes. Genome synteny was evaluated using R2CAT software. Protein conservation was analyzed using mGenome Subtractor. Horizontal gene transfer was predicted through the GOHTAM pipeline (using tetranucleotide composition) and phylogenetic analyses (by maximum likelihood). Hydrolases were identified through the MEROPS and CAZy platforms. RESULTS: The three genomes of C. subterraneus showed high similarity, although there are substantial differences in their gene composition and organization. Each subspecies possesses a gene cluster encoding a carbon monoxide dehydrogenase (CODH) and an energy converting hydrogenase (ECH). The CODH gene is associated with an operon that resembles the Escherichia coli hydrogenase hyc/hyf operons, a novel genetic context distinct from that found in archetypical hydrogenogenic carboxydotrophs. Apart from the CODH-associated hydrogenase, these bacteria also contain other hydrogenases, encoded by ech and hyd genes. An Mbx ferredoxin:NADP oxidoreductase homolog similar to that originally described in the archaeon Pyrococcus furiosus was uniquely encoded in the C. subterraneus subsp. yonseiensis genome. Compositional analysis demonstrated that some genes of the CODH-ECH and mbx operons present distinct sequence patterns in relation to the majority of the other genes of each genome. Phylogenetic reconstructions of the genes from these operons and those from the ech operon are incongruent to the species tree. Notably, the cooS gene of C. subterraneus subsp. pacificus and its homologs in C. subterraneus subsp. tengcongensis and C. subterraneus subsp. yonseiensis form distinct clades. The strains have diverse hydrolytic enzymes and they appear to be proteolytic and glycolytic. Divergent glycosidases from 14 families, among them amylases, chitinases, alpha-glucosidases, beta-glucosidases, and cellulases, were identified. Each of the three genomes also contains around 100 proteases from 50 subfamilies, as well about ten different esterases. CONCLUSIONS: Genomic information suggests that multiple horizontal gene transfers conferred the adaptation of C. subterraneus subspecies to extreme niches throughout the carbon monoxide utilization and hydrogen production. The variety of hydrolases found in their genomes indicate the versatility of the species in obtaining energy and carbon from diverse substrates, therefore these organisms constitute a remarkable resource of enzymes with biotechnological potential.

Counterselection of prokaryotic ribosomal RNA during reverse transcription using non-random hexameric oligonucleotides.

Ribosomal RNA (rRNA) is the major component in total RNA extracts, interfering with the synthesis of cDNA corresponding to messenger RNA (mRNA). In this study, we present a novel strategy for selectively discriminating against rRNA and favoring mRNA from prokaryotes during synthesis of cDNA by reverse transcriptase. Our technique is based on the fact that rRNA sequences, in many species, are G+C rich relative to the genome at large, and highly conserved among prokaryotes. The sequence TTTT is therefore rarely found in rRNA sequences. However, TTTT priming sites are found at a much higher frequency in protein-encoding gene sequences. We designed specific hexamers (HD/DHTTTT) to prime reverse transcription reactions resulting in a selective synthesis of cDNA corresponding to mRNA from prokaryotic total RNA extractions.

[Characterization of membrane-bound Fe(III)-EDTA reductase activities of the thermophilic gram-positive dissimilatory iron-reducing bacterium Thermoterrabacterium ferrireducens].

Whole-cell suspensions of T. ferrireducens reduced Fe(III) citrate, Fe(III)-EDTA, and ferrihydrite with glycerol as an electron donor. After cell disruption, the highest activity was registered with Fe(III)-EDTA as the electron acceptor and NADH or NAD(P)H as electron donors. About 80% of the NAD(P)H-dependent Fe(III)-EDTA reductase activities were associated with the membrane fraction of the cells. Treatment of the membranes with lauryl maltoside led to complete solubilization of the NADH-dependent and 70% solubilization of the NADPH-dependent Fe(III)-EDTA reductase activities. After purification by ion-exchange chromatography, the NADH-dependent activity was concentrated 8-fold, and the NADPH-dependent activity was concentrated 11-fold, with a yield of about 10% for both activities. The Fe(III)-EDTA-reducing enzyme complex included c-type cytochromes and a protein with a molecular mass of ca. 115 k Da, consisting of two polypeptides. This is the first description of membrane-bound Fe(III)-reducing oxidoreductase activities from a gram-positive dissimilatory Fe(III)-reducing bacterium.

Development of two PCR-based techniques for detecting helical and coccoid forms of Helicobacter pylori.

The primary mode of transmission of Helicobacter pylori, a human pathogen carried by more than half the population worldwide, is still unresolved. Some epidemiological data suggest water as a possible transmission route. H. pylori in the environment transforms into a nonculturable, coccoid form, which frequently results in the failure to detect this bacterium in environmental samples by conventional culture techniques. To overcome limitations associated with culturing, molecular approaches based on DNA amplification by PCR have been developed and used for the detection of H. pylori in clinical and environmental samples. Our results showed the glmM gene as the most promising target for detection of H. pylori by PCR amplification. Under optimal amplification conditions, glmM-specific primers generated PCR-amplified products that were specific for H. pylori and some other Helicobacter species. Genome sequence analysis revealed the existence of a conserved region linked to a hypervariable region upstream of the 16S rRNA gene of H. pylori. Selective PCR primer sets targeting this sequence were evaluated for the specific detection of H. pylori. One primer set, Cluster2 and B1J99, were shown to be highly specific for H. pylori strains and did not produce any PCR products when other Helicobacter species and other bacterial species were analyzed. In tests with 32 strains of H. pylori, 6 strains of other Helicobacter species, 8 strains of Campylobacter jejuni, and 21 strains belonging to different genera, the primers for glmM were selective for the Helicobacter genus and the primers containing the region flanking the 16S rRNA gene were selective for H. pylori species only. The combination of two sensitive PCR-based methods, one targeting the glmM gene and the other targeting a hypervariable flanking region upstream of the 16S rRNA gene, are complementary to each other. Whereas the glmM-specific primers provide a rapid, sensitive presumptive assay for the presence of H. pylori and closely related Helicobacter spp., the primers for sequences flanking the 16S rRNA gene can confirm the presence of H. pylori and locate the potential source of this bacterium.

Multi-subunit assembly of the Pyrococcus furiosus small heat shock protein is essential for cellular protection at high temperature.

The hyperthermophilic archaeon, Pyrococcus furiosus, expresses a small, alpha-crystallin-like protein in response to exposure to extreme temperatures, above 103 degrees C. The P. furiosus small heat shock protein (Pfu-sHSP) forms large oligomeric complexes. Based on the available crystal structures of the Methanocaldococcus jannaschii and wheat sHSPs, the protruding carboxy terminal domain is probably involved in subunit interactions. We constructed Pfu-sHSP mutants to analyze chaperone function and to study multi-subunit assembly. The results confirmed that the carboxy terminus of Pfu-sHSP is involved in inter-dimer interactions, whereas the amino terminal deletion mutant still exhibited the wild-type assembly characteristics. The ability to form oligomeric complexes via the carboxy terminal domain was shown to be necessary for thermotolerance of Escherichia coli overexpressing Pfu-sHSP. The amino terminal domain was not required for inter-species thermotolerance.

Mechanism of pressure-induced thermostabilization of proteins: studies of glutamate dehydrogenases from the hyperthermophile Thermococcus litoralis.

In this study, we investigated the effect of pressure on protein structure and stability at high temperature. Thermoinactivation experiments at 5 and 500 atm were performed using the wild-type (WT) enzyme and two single mutants (D167T and T138E) of the glutamate dehydrogenase (GDH) from the hyperthermophile Thermococcus litoralis. All three GDHs were stabilized, although to different degrees, by the application of 500 atm. Interestingly, the degree of pressure stabilization correlated with GDH stability as well as the magnitude of electrostatic repulsion created by residues at positions 138 and 167. Thermoinactivation experiments also were performed in the presence of trehalose. Addition of the sugar stabilized all three GDHs; the degree of sugar-induced thermostabilization followed the same order as pressure stabilization. Previous studies suggested a mechanism whereby the enzyme adopts a more compact and rigid structure and volume fluctuations away from the native state are diminished under pressure. The present results on the three GDHs allowed us to further confirm and refine the proposed mechanism for pressure-induced thermostabilization. In particular, we propose that pressure stabilizes against thermoinactivation by shifting the equilibrium between conformational substates of the GDH hexamer, thus inhibiting irreversible aggregation.

Regulation and mechanism of action of the small heat shock protein from the hyperthermophilic archaeon Pyrococcus furiosus.

The small heat shock protein (sHSP) from the hyperthermophile Pyrococcus furiosus was specifically induced at the level of transcription by heat shock at 105 degrees C. The gene encoding this protein was cloned and overexpressed in Escherichia coli. The recombinant sHSP prevented the majority of E. coli proteins from aggregating in vitro for up to 40 min at 105 degrees C. The sHSP also prevented bovine glutamate dehydrogenase from aggregating at 56 degrees C. Survivability of E. coli overexpressing the sHSP was enhanced approximately sixfold during exposure to 50 degrees C for 2 h compared with the control culture, which did not express the sHSP. Apparently, the sHSP confers a survival advantage on mesophilic bacteria by preventing protein aggregation at supraoptimal temperatures.

Carboxydobrachium pacificum gen. nov., sp. nov., a new anaerobic, thermophilic, CO-utilizing marine bacterium from Okinawa Trough.

A new anaerobic, thermophilic, CO-utilizing marine bacterium, strain JMT, was isolated from a submarine hot vent in Okinawa Trough. Cells of strain JMT were non-motile thin straight rods, sometimes branching, with a cell wall of the Gram-positive type, surrounded with an S-layer. Chains of three to five cells were often observed. The isolate grew chemolithotrophically on CO, producing equimolar quantities of H2 and CO2 (according to the equation CO+H2O-->CO2+H2) and organotrophically on peptone, yeast extract, starch, cellobiose, glucose, galactose, fructose and pyruvate, producing H2, acetate and CO2. Growth was observed from 50 to 80 degrees C with an optimum at 70 degrees C. The optimum pH was 6.8-7.1. The optimum concentration of sea salts in the medium was 20.5-25.5 g l(-1). The generation time under optimal conditions was 7.1 h. The DNA G+C content was 33 mol %. Growth of isolate JMT was not inhibited by penicillin, but ampicillin, streptomycin, kanamycin and neomycin completely inhibited growth. The results of 16S rDNA sequence analysis revealed that strain JMT belongs to the Thermoanaerobacter phylogenetic group within the Bacillus-Clostridium subphylum of Gram-positive bacteria but represents a separate branch of this group. On the basis of morphological and physiological features and phylogenetic data, this isolate should be assigned to a new genus, for which the name Carboxydobrachium is proposed. The type species is Carboxydobrachium pacificum; the type strain is JMT (= DSM 12653T).

Evidence of recent lateral gene transfer among hyperthermophilic archaea.

A total of 153 nucleotide differences were found over a contiguous 16 kb region between two hyperthermophilic Archaea, Pyrococcus furiosus and Thermococcus litoralis. The 16 kb region in P. furiosus is flanked by insertion sequence (IS) elements with inverted and direct repeats. Both IS elements contain a single open reading frame (ORF) encoding a putative protein of 233 amino acids identified as a transposase. This 16 kb region has the features of a typical bacterial composite transposon and represents a possible mechanism for lateral gene transfer between Archaea or possibly between Archaea and Bacteria. A total of 23 homologous IS elements was found in the genome sequence of P. furiosus, whereas no full-length IS elements were identified in the genomes of Pyrococcus abyssi and Pyrococcus horikoshii. Only one IS element was found in T. litoralis. In P. furiosus and T. litoralis, the 16 kb region contains an ABC transport system for maltose and trehalose that was characterized biochemically for T. litoralis. Regulation of expression studies showed that the malE gene, located on the transposon, and the encoded trehalose/maltose-binding protein (TMBP) are induced in the presence of maltose and trehalose in both P. furiosus and T. litoralis. The implications of transposition as a mechanism for lateral gene transfer among Archaea are discussed.

Genetic analysis of Carboxydothermus hydrogenoformans carbon monoxide dehydrogenase genes cooF and cooS.

Carboxydothermus hydrogenoformans is an extremely thermophilic, Gram-positive bacterium growing on carbon monoxide (CO) as single carbon and energy source and producing only H(2) and CO(2). Carbon monoxide dehydrogenase is a key enzyme for CO metabolism. The carbon monoxide dehydrogenase genes cooF and cooS from C. hydrogenoformans were cloned and sequenced. These genes showed the highest similarity to the cooF genes from the archaeon Archaeoglobus fulgidus and the cooS gene from the bacterium Rhodospirillum rubrum, respectively. The cooS gene was identified immediately downstream of cooF, however, the cooF and cooS genes from C. hydrogenoformans have substantially different codon usage, and the cooF gene Arg codon usage pattern, dominated by AGA and AGG, resembles the archaeal pattern. The data therefore suggest lateral transfer of these genes, possibly from different donor species.

Enzymes of hydrogen metabolism in Pyrococcus furiosus.

The genome of Pyrococcus furiosus contains the putative mbhABCDEFGHIJKLMN operon for a 14-subunit transmembrane complex associated with a Ni-Fe hydrogenase. Ten ORFs (mbhA-I and mbhM) encode hydrophobic, membrane-spanning subunits. Four ORFs (mbhJKL and mbhN) encode putative soluble proteins. Two of these correspond to the canonical small and large subunit of Ni-Fe hydrogenase, however, the small subunit can coordinate only a single iron-sulfur cluster, corresponding to the proximal [4Fe-4S] cubane. The structural genes for the small and the large subunits, mbhJ and mbhL, are separated in the genome by a third ORF, mbhK, encoding a protein of unknown function without Fe/S binding. The fourth ORF, mbhN, encodes a 2[4Fe-4S] protein. With P. furiosus soluble [4Fe-4S] ferredoxin as the electron donor the membranes produce H2, and this activity is retained in an extracted core complex of the mbh operon when solubilized and partially purified under mild conditions. The properties of this membrane-bound hydrogenase are unique. It is rather resistant to inhibition by carbon monoxide. It also exhibits an extremely high ratio of H2 evolution to H2 uptake activity compared with other hydrogenases. The activity is sensitive to inhibition by dicyclohexylcarbodiimide, an inhibitor of NADH dehydrogenase (complex I). EPR of the reduced core complex is characteristic for interacting iron-sulfur clusters with Em approximately -0.33 V. The genome contains a second putative operon, mbxABCDFGHH'MJKLN, for a multisubunit transmembrane complex with strong homology to the mbh operon, however, with a highly unusual putative binding motif for the Ni-Fe-cluster in the large hydrogenase subunit. Kinetic studies of membrane-bound hydrogenase, soluble hydrogenase and sulfide dehydrogenase activities allow the formulation of a comprehensive working hypothesis of H2 metabolism in P. furiosus in terms of three pools of reducing equivalents (ferredoxin, NADPH, H2) connected by devices for transduction, transfer, recovery and safety-valving of energy.

Novel structure and redox chemistry of the prosthetic groups of the iron-sulfur flavoprotein sulfide dehydrogenase from Pyrococcus furiosus; evidence for a [2Fe-2S] cluster with Asp(Cys)3 ligands.

The consecutive structural genes for the iron-sulfur flavoenzyme sulfide dehydrogenase, sudB and sudA, have been identified in the genome of Pyrococcus furiosus. The translated sequences encode a heterodimeric protein with an alpha-subunit, SudA, of 52598 Da and a beta-subunit, SudB, of 30686 Da. The alpha-subunit carries a FAD, a putative nucleotide binding site for NADPH, and a [2Fe-2S]2+,+ prosthetic group. The latter exhibit EPR g-values, 2.035, 1.908, 1.786, and reduction potential, Em,8 = +80 mV, reminiscent of Rieske-type clusters; however, comparative sequence analysis indicates that this cluster is coordinated by a novel motif of one Asp and three Cys ligands. The motif is not only found in the genome of hyperthermophilic archaea and hyperthermophilic bacteria, but also in that of mesophilic Treponema pallidum. The beta-subunit of sulfide dehydrogenase contains another FAD, another putative binding site for NADPH, a [3Fe-4S]+,0 cluster, and a [4Fe-4S]2+,+ cluster. The 3Fe cluster has an unusually high reduction potential, Em,8 = +230 mV. The reduced 4Fe cluster exhibits a complex EPR signal, presumably resulting from magnetic interaction of its S = 1/2 spin with the S=2 spin of the reduced 3Fe cluster. The 4Fe cluster can be reduced with deazaflavin/EDTA/light but not with sodium dithionite; however, it is readily reduced with NADPH. SudA is highly homologous to KOD1-GO-GAT (or KOD1-GltA), a single-gene encoded protein in Pyrococcus kodakaraensis, which has been putatively identified as hyperthermophilic glutamate synthase. However, P. furiosus sulfide dehydrogenase does not have glutamate synthase activity. SudB is highly homologous to HydG, the gamma-subunit of P. furiosus NiFe hydrogenase. The latter enzyme also has sulfide dehydrogenase activity. The P. furiosus genome contains a second set of consecutive genes, sudY and sudX, with very high homology to the sudB and sudA genes, and possibly encoding a sulfide dehydrogenase isoenzyme. Each subunit of sulfide dehydrogenase is a primary structural paradigm for a different class of iron-sulfur flavoproteins.

Cloning and sequence analysis of the mercury resistance operon of Streptomyces sp. Strain CHR28 reveals a novel putative second regulatory gene.

A DNA library of pRJ28, a large linear plasmid encoding mercury resistance, was constructed, and the mercury resistance genes were cloned. The 5,921-bp sequence was analyzed and showed a high degree of similarity to the Streptomyces lividans 1326 mercury resistance operon. Genes merR, merT, merP, and orfIV were found in a similar order and in a single transcription unit. merA and merB were found to be transcribed in the opposite direction to genes merR, merT, merP, and orfIV, as in S. lividans 1326. A novel putative regulatory gene, orfX, was found 22 bp downstream of merA. orfX encodes a 137-amino acid protein with a potential helix-turn-helix motif in the N-terminal domain, characteristic of the MerR family of transcriptional regulators. Transcriptional studies showed that orfX is cotranscribed with merA and merB. It is hypothesized that orfX plays a role in the regulation of the mercury resistance operon, probably by binding at the MerR operator site.

Polyphasic taxonomy of a novel Halobacillus, Halobacillus thailandensis sp. nov. isolated from fish sauce.

In order to speed up fish sauce production, a more complete understanding of the microorganisms associated with the fermentation was needed. This study was undertaken to meet that need. A bacterium was isolated from a fish sauce production line containing 25% NaCl. It is a Gram-positive, rod-shaped bacillus with pointed ends, occurring as single cells, pairs, or short chains. Endospores are produced on a low nutrient medium and, in old cultures, the cells round up, even when undergoing division. The cell wall is relatively amorphous and similar to that of Gram-positive bacteria in structure and composition. Cells grown in a medium containing 10-20% salt possess thicker cell walls than those grown in a medium with 3% salt. Based on 16S rRNA sequence and DNA/DNA hybridization data, we conclude that the bacterium is a species of Halobacillus. This bacterium shares 99.2% and 97.2% 16S rRNA similarity with Halobacillus litoralis and Halobacillus halophilus respectively and DNA/DNA homology was lower than 70%, considered indicative of species similarity. Three highly expressed extra-cellular proteolytic enzymes with M(r) of approximately 100 kDa, 42 kDa and 17 kDa, respectively, were detected by gelatin-polyacrylamide gel electrophoresis. Activity of the 100 kDa and 17 kDa proteases was inhibited by phenylmethanesulphonyl fluoride (PMSF), without being affected by L-trans epoxysuccinyl-leucylamide 4-guanidino-butane (E-64), pepstatin, EDTA, or 1, 10-phenanthroline, leading to the conclusion that these enzymes are serine proteases. The 42-kDa protease was inhibited by EDTA and 1,10-phenanthroline, but not by PMSF, thus, being classified a metalloprotease. The strain has been successfully employed to improve fermentation in industrial production of fish sauce in Thailand.

Structure determination of the glutamate dehydrogenase from the hyperthermophile Thermococcus litoralis and its comparison with that from Pyrococcus furiosus.

Glutamate dehydrogenase catalyses the oxidative deamination of glutamate to 2-oxoglutarate with concomitant reduction of NAD(P)(+), and has been shown to be widely distributed in nature across species ranging from psychrophiles to hyperthermophiles. Extensive characterisation of this enzyme isolated from hyperthermophilic organisms has led to its adoption as a model system for analysing the determinants of thermal stability. The crystal structure of the extremely thermostable glutamate dehydrogenase from Thermococcus litoralis has been determined at 2.5 A resolution, and has been compared to that from the hyperthermophile Pyrococcus furiosus. The two enzymes are 87 % identical in sequence, yet differ 16-fold in their half-lives at 104 degrees C. This is the first reported comparative analysis of the structures of a multisubunit enzyme from two closely related yet distinct hyperthermophilies. The less stable T. litoralis enzyme has a decreased number of ion pair interactions; modified patterns of hydrogen bonding resulting from isosteric sequence changes; substitutions that decrease packing efficiency; and substitutions which give rise to subtle but distinct shifts in both main-chain and side-chain elements of the structure. This analysis provides a rational basis to test ideas on the factors that confer thermal stability in proteins through a combination of mutagenesis, calorimetry, and structural studies.

Rapid extraction of plasmid pGT5 from the hyperthermophilic archaeon Pyrococcus abyssi.

Hyperthermophilic archaea, specifically Pyrococcus spp., are the target of current efforts in developing heterologous expression systems. However, the published plasmid purification and plasmid screening protocols are long and tedious. We describe a fast, simple protocol for plasmid purification from Pyrococcus spp. developed while extracting the plasmid pGT5 from Pyrococcus abyssi cells. The protocol is modified from the procedures for commercial plasmid minipreps and is completed in about 20 min. The DNA is easily digested by restriction enzymes and can be used in sequencing reactions without additional purification.

DNA repair systems in archaea: mementos from the last universal common ancestor?

DNA repair in the Archaea is relevant to the consideration of genome maintenance and replication fidelity in the last universal common ancestor (LUCA) from two perspectives. First, these prokaryotes embody a mix of bacterial and eukaryal molecular features. Second, DNA repair proteins would have been essential in LUCA to maintain genome integrity, regardless of the environmental temperature. Yet we know very little of the basic molecular mechanisms of DNA damage and repair in the Archaea in general. Many studies on DNA repair in archaea have been conducted with hyperthermophiles because of the additional stress imposed on their macromolecules by high temperatures. In addition, of the six complete archaeal genome sequences published so far, five are thermophilic archaea. We have recently shown that the hyperthermophile Pyrococcus furiosus has an extraordinarily high capacity for repair of radiation-induced double-strand breaks and we have identified and sequenced several genes involved in DNA repair in P. furiosus. At the sequence level, only a few genes share homology with known bacterial repair genes. For instance, our phylogenetic analysis indicates that archaeal recombinases occur in two paralogous gene families, one of which is very deeply branched, and both recombinases are more closely related to the eukaryotic RAD51 and Dmc1 gene families than to the Escherichia coli recA gene. We have also identified a gene encoding a repair endo/exonuclease in the genomes of several Archaea. The archaeal sequences are highly homologous to those of the eukaryotic Rad2 family and they cluster with genes of the FEN-1 subfamily, which are known to be involved in DNA replication and repair in eukaryotes. We argue that there is a commonality of mechanisms and protein sequences, shared between prokaryotes and eukaryotes for several modes of DNA repair, reflecting diversification from a minimal set of genes thought to represent the genome of the LUCA.

Divergence of the hyperthermophilic archaea Pyrococcus furiosus and P. horikoshii inferred from complete genomic sequences.

Divergence of the hyperthermophilic Archaea, Pyrococcus furiosus and Pyrococcus horikoshii, was assessed by analysis of complete genomic sequences of both species. The average nucleotide identity between the genomic sequences is 70-75% within ORFs. The P. furiosus genome (1.908 mbp) is 170 kbp larger than the P. horikoshii genome (1.738 mbp) and the latter displays significant deletions in coding regions, including the trp, his, aro, leu-ile-val, arg, pro, cys, thr, and mal operons. P. horikoshii is auxotrophic for tryptophan and histidine and is unable to utilize maltose, unlike P. furiosus. In addition, the genomes differ considerably in gene order, displaying displacements and inversions. Six allelic intein sites are common to both Pyrococcus genomes, and two intein insertions occur in each species and not the other. The bacteria-like methylated chemotaxis proteins form a functional group in P. horikoshii, but are absent in P. furiosus. Two paralogous families of ferredoxin oxidoreductases provide evidence of gene duplication preceding the divergence of the Pyrococcus species.