| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
,
Toshiaki Kudo,1
Sadaharu Ui,2 and
Moriya Ohkuma1*
Environmental Molecular Biology Laboratory, RIKEN, Wako, Saitama 351-0198,1 Division of Biotechnology, Interdisciplinary Graduate School of Medicine and Engineering, University of Yamanashi, Yamanashi 400-8510, Japan2
Received 12 July 2007/ Accepted 13 August 2007
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
Nevertheless, little is known about the detailed physiology and biochemistry of the gut symbiotic protists of termites; this is because most of these protists have formidably resisted cultivation in vitro. Culture-independent molecular studies, however, have been successfully applied to the gut dwellers of termites (for reviews, see reference 26); these studies were mainly molecular phylogenetic analyses. They revealed unexpected diversity among termite gut protists (12, 27-29). Genes involved in cellulose decomposition have also been identified previously (references 18 and 36 and references therein).
The gut symbiotic protists are unique to termites and related wood-feeding roaches and belong to the phylum Parabasalia or the order Oxymonadida (phylum Preaxostyla). The former includes most of the large gut protists that are responsible for the decomposition of highly polymerized cellulose. Parabasalids are typified by the presence of parabasal bodies and hydrogenosomes in the cell (4). Hydrogenosomes are H2-producing organelles that play a central role in anaerobic energy production (19) and are distributed in a variety of anaerobic eukaryotes; therefore, they are considered to have evolved multiple times independently from the mitochondria (8, 11).
All the known eukaryotic hydrogenases, whether they are located in the hydrogenosome or not, are classified as iron-only hydrogenases ([FeFe] hydrogenases; herein referred to simply as Fe hydrogenases), which are evolutionarily unrelated to the other major type of hydrogenases, [NiFe] hydrogenases (37). Fe hydrogenases are considered to play a role in H2 evolution by reducing protons; however, the eukaryotic enzymes are poorly understood because, thus far, only a few eukaryotic Fe hydrogenases, namely, those in Chlamydomonas and Trichomonas species, have been purified and characterized (13, 30), although a higher number of primary structures have been deduced from the genes of various anaerobic eukaryotes (1, 5, 14, 21, 38). As proposed previously (11), the gut-dwelling anaerobic eukaryotes of termites and other animals potentially harbor an enormous biodiversity of Fe hydrogenases due to their great phylogenetic diversity. Further research regarding the various anaerobic hydrogenosome-carrying eukaryotes and their hydrogenases is critical for understanding eukaryotic organelle evolution and the evolutionary origin of eukaryotic Fe hydrogenases (11).
In this study, we investigated the Fe hydrogenase genes and the characteristics of the encoded enzymes in parabasalian symbionts of the termite Coptotermes formosanus. The study aimed to understand the molecular mechanisms underlying the intensive H2 production in termites and to investigate the biodiversity of hydrogenases. The termite C. formosanus is a globally distributed and economically important pest of wood houses. It harbors only three species of gut protists, all of which belong to the phylogenetically monophyletic phlyum Parabasalia but are distantly related to the cultivated Parabasalia representatives Trichomonas vaginalis and Tritrichomonas foetus (28, 29). Among the three, the largest protist, Pseudotrichonympha grassii, is maximally responsible for the decomposition of highly polymerized cellulose (40). Recombinant enzyme forms of two distinct Fe hydrogenases of P. grassii were purified and characterized.
| MATERIALS AND METHODS |
|---|
|
|
|---|
The termites were fed an artificial diet that comprised cellulose and water in the ratio of 2:3 for 1 week. Ten termites were placed in 100-ml serum vials that contained 5 g of the artificial diet; the vials were then sealed with butyl rubber stoppers. The H2 concentration in the headspace of the vials was measured every 12 h for 72 h by gas chromatography with a Shimadzu 14C chromatograph and with a semiconductor detector and a Unibead C packed column (GL Science) operating at 60 and 80°C, respectively. The carrier was nitrogen gas flowing at 60 ml/min. Linear emission continued (R2 = 0.997) for at least 72 h.
Isolation of cDNA clones and 5'-end amplification. The sequences of the cDNAs encoding hydrogenase homologs were identified by BLAST searches of the nonredundant NCBI database as part of our ongoing expressed-sequence-tag (EST) survey. The cDNA library used for the EST analysis has been established from the mixed population of gut protists in C. formosanus (18). The truncated 5' ends of cDNAs were determined by 5' rapid amplification of cDNA end (RACE) reactions using a Gene Racer kit (Invitrogen) with sequence-specific primers (see Table S1 in the supplemental material). The 5' RACE products were cloned, and the DNA sequences of multiple clones were determined to confirm the identity of the overlapping cDNA region of the original clones. The entire coding sequences were amplified by reverse transcription-PCR (RT-PCR) with sequence-specific primers (see Table S1 in the supplemental material), and the PCR products were used for subsequent analyses. DNA sequencing was performed by primer walking using an automated DNA sequencer (ABI model 3700) with BigDye Terminator chemistry (Applied Biosystems).
Identification of protist origins. Cells of each species of symbiotic protists in the gut were manually isolated with the aid of a micromanipulator (Cell-Tram; Eppendorf) as described previously (27), and the cells were washed extensively. As described previously (29), a pool of five isolated cells of each protist species was used directly for specific amplification by RT-PCR with sequence-specific primers (see Table S1 in the supplemental material) to detect the species of origin of the sequences. The amplification product from each protist species was cloned, and the DNA sequences in multiple clones were determined to confirm the identity to the corresponding cDNA sequence. Whole-cell hybridizations using oligonucleotide probes specific for 09A82 and 06B09 sequences (see Table S1 in the supplemental material) were applied for the in situ identification of the organismal origins of the Fe hydrogenase genes. Fixation and hybridization were performed as described previously (18). Due to the small amount of mRNA for Fe hydrogenase genes, hybridization signals were detected by an enzyme-linked immunoassay. The cells that hybridized with probes 5'-end labeled with digoxigenin (DIG) were detected with alkaline phosphatase-conjugated anti-DIG antibodies and colorimetric substrates by using a DIG nucleic acid detection kit (Roche) and observed by light microscopy.
Phylogenetic analysis. The in silico-translated amino acid sequences corresponding to the genes for the Fe hydrogenase homologs were aligned with sequences of known Fe hydrogenases and those of their homologous proteins listed in the protein databases based on the Pfam seed alignment of Fe hydrogenases (http://www.sanger.ac.uk/Software/Pfam/), and the alignment was corrected manually. The program PHYML (version 2.4.4) (10) was used for inferring the maximum-likelihood (ML) tree using the Jones-Taylor-Thornton model, gamma distribution with four categories plus invariant sites, and the shape parameter and the fraction of invariable sites estimated from the data set. MrBayes (version 3.1.2) (15) was used for evaluating the posterior probabilities at the nodes of the ML tree. MrBayes was run for 1,000,000 generations with the Jones-Taylor-Thornton model, gamma-distributed rate variation, and a proportion of invariable sites; the burn-in value was set at 100,000 generations. The significance of the differences in tree topologies was examined by the approximately unbiased test implemented in CONSEL (34) with 1,000 bootstrap replicates. All the Fe hydrogenase sequences from Parabasalia were constrained to monophyly, and the topology of an ML tree inferred with MrBayes under the constraint was compared with that of an unconstrained ML tree.
Production of recombinant enzymes. Hydrogenosome-targeting signals in the amino acid sequences of the Fe hydrogenases of the termite gut symbionts were predicted as mitochondrial targeting signals by using the programs PSORT (http://psort.ims.u-tokyo.ac.jp/), TargetP (http://www.cbs.dtu.dk/services/SignalP/), and MitoProt (http://ihg.gsf.de/ihg/mitoprot.html). The coding regions of the Fe hydrogenase cDNAs of P. grassii without the predicted hydrogenosome-targeting signals were amplified by PrimeSTAR DNA polymerase (Takara) using sequence-specific primers (see Table S1 in the supplemental material), and the amplification products were subcloned into the pET47b vector (Novagen) by using the restriction sites (SmaI and NotI) introduced into these primers. Escherichia coli BL21(DE3) cells carrying the resulting plasmids were cultured at 20°C in Luria broth for the production of His6-tagged recombinant enzymes. The cells were grown first aerobically for 12 h and then anaerobically for 40 h. The grown cells were transferred into fresh medium and grown anaerobically, and recombinant enzyme production was induced by the addition of isopropyl-ß-D-thiogalactopyranoside (IPTG). IPTG was added at 48 and 60 h for the production of the 09A82 and 06B09 proteins, respectively. After the 12-h induction, the cells were harvested by centrifugation and disrupted by sonication to extract the cell lysate. The recombinant proteins were purified using a His trap nickel affinity purification column (GE Healthcare). All these manipulations were performed anaerobically in an anaerobic chamber (Bactron). Purity was assessed by measuring the absorption of the proteins at 280 and 400 nm, the latter of which was specific for Fe-S cluster proteins, as well as by performing sodium dodecyl sulfate-polyacrylamide gel electrophoresis.
Fractionation of protist cells and their cellular components. The termites were fed an artificial diet for 1 week, and 100 termites were used for fractionating cellular components of the protist species. All of the following manipulations were performed anaerobically. Termite guts were dissected using a fine-tipped forceps and gently crushed in the solution described previously (28). After the release of the gut contents and the removal of the gut tissues, the P. grassii cells were concentrated by low-speed centrifugation at 18 x g for 5 min in a Kubota 1-13 centrifuge and the pellet was resuspended and washed three times. The concentrated P. grassii cells were mildly broken by the addition of Nonidet P-40 (Nacalai Tesque) at a final concentration of 0.8%. A fraction of the large particles (nuclei, flagella, and axostyle fragments) was obtained by centrifugation at 500 x g for 4 min. From the supernatant produced by the 500 x g centrifugation, the hydrogenosome-enriched fraction was obtained by centrifugation at 2,000 x g for 10 min. The fraction enriched with the intracellular bacteria of P. grassii was obtained by centrifugation at 10,000 x g for 10 min. The presence of the bacteria was confirmed by 4',6-diamidino-2-phenylindole (DAPI) staining and fluorescence microscopy. Each fraction was homogenized using a sonic homogenizer (model UR-20P; Tomy Seiko) and used for the hydrogenase assay.
Hydrogenase assay. The hydrogenases were evaluated by measuring both H2 evolution and H2 uptake activities. The H2 evolution activity was measured by gas chromatography, as described above, using reduced methyl viologen (MV) as an electron donor. The H2 uptake activity was measured by spectroscopy analysis of the H2-dependent reduction of methylene blue (MB), and the value was used to calculate the molar H2 consumed. The reactions were performed in a tightly closed quartz cuvette filled with 500 µl of the reaction mixture (50 mM buffer, 2 mM dithiothreitol, 1 mM electron carrier, and dithionite). The dithionite concentrations used for the H2 evolution and H2 uptake assays were 15 and 1 mM, respectively. The reaction mixture was bubbled for 15 min with nitrogen for the H2 evolution assay or with H2 for the H2 uptake assay. Measurement was initiated by the addition of 10 µl of the enzyme solution. All these activities were determined in triplicate. The protein concentration was determined by the DC protein assay (Bio-Rad) using bovine serum albumin as a standard.
Optimum pH was examined using 50 mM concentrations of acetate buffer, phosphate buffer, Tris-HCl buffer, and glycine-NaOH buffer at pH ranges of 5.0 to 6.5, 6.0 to 7.5, 7.0 to 9.0, and 8.5 to 9.5, respectively. In the H2 uptake assay, the specificity of the electron carrier was examined with benzyl viologen, flavin adenine dinucleotide (FAD), flavin mononucleotide (FMN), NAD, NADH, NADP, NADPH, and ferredoxin from Clostridium pasteurianum (Sigma). Only benzyl viologen, FAD, FMN, and ferredoxin were examined in the H2 evolution assay. Enzyme stability was investigated by incubating the enzyme solution at different pHs or temperatures for 5 min prior to the assay. A temperature range from 5 to 60°C and a pH range from 5 to 10 were used to investigate enzyme stability.
The effects of high H2 partial pressure on the H2 evolution activity of the recombinant 09A82 enzyme were evaluated. The headspace of a quartz cuvette containing the reaction mixture described above was first filled and compressed appropriately with nitrogen, and then the H2 partial pressure was adjusted to 10 to 50% of the final total pressure of 200 kPa. The enzyme reaction was performed at 30°C, and the increase in the headspace H2 was measured by gas chromatography.
Nucleotide sequence accession numbers. The DNA sequences of four homologous Fe hydrogenase cDNAs have been deposited in the DDBJ database under the accession numbers AB331667 to AB331670.
| RESULTS |
|---|
|
|
|---|
Fe hydrogenase genes and their origin. In the EST analysis of the mixed population of gut symbiotic protists of C. formosanus, four cDNA clones (06B09, 09A82, 11B87, and 08B89) that encoded proteins homologous to known Fe hydrogenases were discovered. No sequence homologous to the [NiFe] hydrogenases was found in the ESTs examined. In 06B09, 09A82, and 11B87, the entire coding regions were identified by sequence-specific 5' RACE and subsequent RT-PCR, while the putative translational initiation codon was found in the original cDNA clone of 08B89. The open reading frames in 06B09, 09A82, 11B87, and 08B89 encoded 550, 476, 456, and 447 amino acid residues, respectively.
The flagellated protists in the gut of C. formosanus comprise exclusively the three parabasalian species P. grassii, Holomastigotoides mirabile, and Spirotrichonympha leidyi (28, 29). P. grassii belongs to the family Eucomonymphidae, while the latter two belong to the family Spirotrichonymphidae. The protist origins of the sequences 06B09 and 09A82 were examined by in situ hybridizations with sequence-specific probes (Fig. 1). Each probe gave positive signals only in the P. grassii cells, indicating that P. grassii was the origin of these sequences. The protist origin was also confirmed by sequence-specific RT-PCR using manually isolated cells of each protist species. The specific amplification of 06B09 and 09A82 from P. grassii and of 08B89 from H. mirabile was detected, but no amplification of 11B87 from any of the protist species was observed. The specific amplification products showed >99.6% identity to the respective cDNA sequences. Further in situ hybridization experiments are necessary to identify the protist origins of the 08B89 and 11B87 sequences.
|
A phylogenetic analysis of the conserved catalytic large subunit domain (Fig. 2) revealed that all the four sequences obtained from the gut symbionts belonged to the eukaryotic group of Fe hydrogenases that also contains some sequences from bacteria, such as Thermotoga maritima and Treponema denticola, but does not include the sequence from the Nyctotherus species ciliate. 09A82, 11B87, and 08B89 formed a robust monophyletic cluster with the sequences from T. vaginalis, indicating that they share a common ancestor as predicted by their classification in Parabasalia. However, the 06B09 sequence was distantly related to the cluster of parabasalid sequences and formed a lineage unique among the eukaryotic Fe hydrogenases. Monophyly of 06B09 with the other sequences of parabasalids was rejected by the approximately unbiased test (P < 0.05). The phylogenetic analyses of longer regions (the large and small subunits, the large subunit plus 4Fe-4S clusters, and all of these regions) always exhibited the monophyly of the 09A82, 08B89, and 11B87 enzymes and T. vaginalis Fe hydrogenases and the distant relationship of the 06B09 protein to this monophyletic cluster (data not shown).
|
|
|
Subcellular fractionation of hydrogenase activity. The localization and properties of native Fe hydrogenases were investigated by subcellular fractionation of the protist cells. The mixed population of hindgut microorganisms was enriched with P. grassii cells, which were used for the small-scale subcellular fractionation after the unassociated bacteria were washed out. Since P. grassii harbors a high density of intracellular bacteria (22, 23), the bacterial and organelle fractions were roughly fractionated by differential centrifugation (Table 2). The organelle-rich fraction (the 2,000 x g precipitate [ppt]) exhibited the highest H2 evolution activity at pH 8, while the highest H2 uptake activity of the bacterium-rich fraction was observed at both pH 8 and pH 6.
|
| DISCUSSION |
|---|
|
|
|---|
The levels of H2 evolution activity exhibited by the Fe hydrogenases of P. grassii, particularly that of the short-form 09A82 enzyme, were some of the highest among those demonstrated by known Fe hydrogenases. The activity of the P. grassii Fe hydrogenases, however, is probably underestimated, because the recombinant enzymes are considerably unstable. These enzymes lost their activity rapidly when exposed to ambient temperature even under strict anaerobic conditions. The loss of activity was irreversible, probably due to the decomposition of Fe-S clusters that was revealed by the decay of the absorption spectrum at around 400 nm. Electron paramagnetic resonance (EPR) spectroscopy analysis of the 09A82 enzyme indicated that only 44% of the protein formed an authentic H-cluster (the presence of an EPR signal with principal g factors of 2.032, 2.017, and 1.998 confirmed that the protein is a typical Fe hydrogenase).
It is generally believed that the production of heterologous Fe hydrogenase in E. coli is difficult since E. coli lacks auxiliary maturases for the metal clusters. In Chlamydomonas, such assembly factors have been shown to be necessary for active Fe hydrogenase production (32). In T. vaginalis, the homologous maturases are reported to localize in the hydrogenosomes (33), although their contribution to active Fe hydrogenase assembly has not been estimated. The P. grassii Fe hydrogenases probably assembled, albeit imperfectly, with the aid of unknown E. coli machinery or spontaneously without auxiliary maturases. Similarly, an Fe hydrogenase of Entamoeba histolytica is reported to be active when expressed in E. coli, although it shows only weak activity (21).
The presence of the two Fe hydrogenases in P. grassii and their distinct characteristics in terms of optimal pH and domain structure suggest that they play different physiological roles. The difference in the accessory domains in the N termini may affect the redox potentials of the Fe hydrogenases since these domains conduct electrons between the donors and the H2-activating center (20, 31). The in vivo electron donors may be ferredoxin reduced by pyruvate-ferredoxin oxidoreductase or involve NADH dehydrogenase as predicted previously for T. vaginalis (6). Further studies are necessary to address the physiological differences. The localization of both P. grassii enzymes in hydrogenosomes also needs to be ascertained.
A notable characteristic of the P. grassii Fe hydrogenases was that they preferentially catalyzed H2 evolution. Another interesting feature was that the evolution activity of the short-form 09A82 enzyme was resistant to high H2 partial pressure. In addition to their advantageous application for H2 production, these characteristics are suitable for the hydrogenosomal function of anaerobic energy production. Since the elimination of reducing equivalents is an important rate-limiting step for anaerobic fermentation, the efficient production of H2 stimulates cellulose fermentation by the protist cells. The reducing equivalents produced during the fermentation as forms of NADH and reduced ferredoxin are rapidly eliminated by reoxidizing these electron carriers through the intensive H2 production activity. Indeed, the hydrogenosome fraction of P. grassii cells exhibited strong H2 evolution activity at pH 8—the optimal pH for the short-form 09A82 enzyme. The simultaneous presence of strong H2 uptake activity in the endosymbiotic bacterial fraction (Table 2), however, appears to disturb the evolution activity. The H2 evolution activities of both P. grassii enzymes at pH 6 are lower than the evolution activity of the short-form enzyme at pH 8; similarly, H2 evolution activity in hydrogenosomes at pH 6 is probably lower than that at pH 8, and the uptake enzyme quickly utilizes the produced H2.
The potentially strong H2 evolution activity of the gut protists and the considerable H2 uptake activity of the endosymbiotic bacteria suggest that interspecies H2 transfer is involved in the symbiotic relationship between the host protist and the endosymbiont that has been shown to be a novel member of the "Bacteroidales" (22, 23). This novel aspect of the intracellular symbiotic relationship is beneficial to the host since, as shown for intracellular methanogens in anaerobic protists (9), the rapid removal of the reducing equivalent (H2) by the associated bacteria leads to enhanced fermentation by the host protists. The electron acceptor(s) for H2 utilization in the endosymbiotic bacteria remains to be identified in order to understand the true nature of this symbiotic relationship.
In addition to the H2 uptake by the endosymbiotic bacteria, the presence of hydrogenotrophic populations, such as methanogens and homoacetogens, in the gut microbial community of termites is well known (3). As apparent by the low level of activity in the unfractionated homogenate of P. grassii (Table 1), the H2 evolution ability of gut symbiotic protists is underestimated. Nevertheless, termites emit H2 at a considerable rate, as shown in this and previous studies (35), and this is a unique characteristic of natural ecosystems. Their potentially strong H2 evolution activity as well as their efficient cellulolytic ability indicates that termites and their gut symbiotic protists represent an attractive biological model of H2 production from cellulosic dead plant matter.
| ACKNOWLEDGMENTS |
|---|
J.-I.I. was supported by the Junior Research Associate Program from RIKEN. This work was partially supported by grants-in-aid for scientific research from JSPS (no. 16380065 and 19380055) and a grant from the Noda Institute for Scientific Research to M.O. and by grants for the Bioarchitect Research and the EcoMolecular Science Research programs from RIKEN.
| FOOTNOTES |
|---|
Published ahead of print on 31 August 2007. ![]()
Supplemental material for this article may be found at http://ec.asm.org/. ![]()
Present address: Kumamoto Industrial Research Institute and Kumamoto University, Kumamoto 862-0901, Japan. ![]()
| REFERENCES |
|---|
|
|
|---|
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| Appl. Environ. Microbiol. | Infect. Immun. | J. Bacteriol. |
|---|---|---|
| Mol. Cell Biol. | Microbiol. Mol. Biol. Rev. | ALL ASM JOURNALS |