Diverse and Novel nifH and nifH-Like Gene Sequences in the Deep ...

APPLIED AND ENVIRONMENTAL MICROBIOLOGY, Apr. 2009, p. 2238–2245 0099-2240/09/$08.00⫹0 doi:10.1128/AEM.02556-08 Copyright © 2009, American Society for Microbiology. All Rights Reserved.

Vol. 75, No. 7

Diverse and Novel nifH and nifH-Like Gene Sequences in the Deep-Sea Methane Seep Sediments of the Okhotsk Sea䌤† Hongyue Dang,1* Xiwu Luan,2* Jingyi Zhao,1 and Jing Li1 Centre for Bioengineering and Biotechnology, China University of Petroleum (East China), Qingdao 266555, China,1 and Key Laboratory of Marine Geology and Environment, Institute of Oceanology, Chinese Academy of Sciences, Qingdao 266071, China2 Received 10 November 2008/Accepted 22 January 2009

Diverse nifH and nifH-like gene sequences were obtained from the deep-sea surface sediments of the methane hydrate-bearing Okhotsk Sea. Some sequences formed novel families of the NifH or NifH-like proteins, of currently unresolved bacterial or archaeal origin. Comparison with other marine environments indicates environmental specificity of some of the sequences, either unique to the methane seep sediments of the Okhotsk Sea or to the general deep-sea methane seep sedimentary environments. 16S rRNA genes were amplified with primers Arch21F and Arch958R (9). Methyl coenzyme M reductase A genes (mcrA) associated with methanogenic or anaerobic methane-oxidizing archaea were amplified with primers ME1 and ME2 (11). Primers AOM39_F and AOM40_R were also used to specifically amplify group b mcrA sequences (12). PCR products were cloned and sequenced by previously established protocols (6, 7), except that TaqI (MBI) was also used as a third restriction enzyme for nifH PCR product restriction fragment length polymorphism analysis. Diversity of nifH and nifH-like sequences. Diverse nifH sequences were obtained (81 restriction fragment length polymorphism sequence types with 41.8 to 98.5% sequence identity among each other) (Table 1). The protein sequences deduced shared 39.2 to 100.0% identity among each other, and 47.1 to 100.0% identity to the closest match GenBank NifH and NifHlike sequences. Using DOTUR at 5% sequence distance cutoff (28), 35 NifH operational taxonomic units (OTUs) were identified (see Fig. S1 in the supplemental material). Twelve OTUs were shared between the two sediment samples, while 14 OTUs were unique to station LV39-25H and 9 OTUs to station LV39-40H. Phylogenetic analysis indicated that almost all the NifH sequences might be obtained from currently uncultured or uncharacterized bacteria or archaea (Fig. 1). Seven NifH and NifH-like clusters were identified (Fig. 1), including three newly defined putative clusters (tentatively named clusters III-x, V, and VI) and all four previously established clusters (clusters I to IV) (23). Cluster V sequences were unique to the deep-sea methane seep sediments of the Okhotsk Sea, and cluster III-x sequences were unique to the general deep-sea methane seep sedimentary environments, including the Okhotsk Sea, Nankai Trough, and Eel River Basin so far studied (25). A single OTU (25H-0N-1) was affiliated within cluster I (Fig. 1), sharing 95.2% identity with a bacterial NifH sequence obtained from the dead biomass of Spartina alterniflora (18). Another single OTU (40H-0N-26) was affiliated within cluster II (Fig. 1), which comprised the NifH sequences from certain methanogenic archaea and the alternative nitrogenase reductases encoded by anfH and vnfH (23).

Shallowly buried methane hydrates sustain a significant biomass and productivity in the deep-sea environment, potentially requiring a large supply of fixed nitrogen (25). The largest reservoir of nitrogen in the ocean is dissolved dinitrogen gas, which is abundant in deep seawater (23). Microbial nitrogen fixation may provide a source of reactive nitrogen in the carbon-rich and nitrogen-limited methane seep environment (15). Besides bacteria, archaea were also found to possess the ability of nitrogen fixation (5), potentially important in oligotrophic open seas (21) and in hydrothermal vent and other deep-sea extreme environments (22–24). Archaeal nitrogen fixation may be important in the deep-sea methane seep environment, although very little is known currently (25). To test this hypothesis, a molecular study of the putative nitrogen-fixing microbial community was carried out for the surface sediments of the methane hydrate-bearing Okhotsk Sea. Sampling and DNA extraction. Sediment core samples were collected from the LV39-25H and LV39-40H stations of the Okhotsk Sea during the CHAOS (hydro-Carbon Hydrate Accumulation in the Okhotsk Sea) international research expedition of 24 May to 18 June 2006 (20). Methane hydrates were discovered in shallow layers of both cores. Undisturbed surface sediments down to a 5-cm depth were sampled using sterile techniques and stored in liquid nitrogen during the cruise and at ⫺80°C after returning to the laboratory. Sediment DNA was extracted by a previously established procedure (8). PCR, cloning, and sequencing. Bacterial and archaeal nitrogenase reductase genes (nifH, including anfH and vnfH) were amplified with primers of Mehta et al. (23). Archaeal

* Corresponding author. Mailing address for Hongyue Dang: Centre for Bioengineering and Biotechnology, China University of Petroleum (East China), Qingdao 266555, China. Phone: 86-532-86981561. Fax: 86-532-86981318. E-mail: [email protected]. Mailing address for Xiwu Luan: Key Laboratory of Marine Geology and Environment, Institute of Oceanology, Chinese Academy of Sciences, Qingdao 266071, China. Phone and fax: 86-532-82898536. E-mail: xluan@ms .qdio.ac.cn. † Supplemental material for this article may be found at http://aem .asm.org/. 䌤 Published ahead of print on 30 January 2009. 2238

nifH AND nifH-LIKE GENES IN DEEP-SEA METHANE SEEP

VOL. 75, 2009

2239

TABLE 1. nifH, 16S rRNA, and mcrA gene clone libraries constructed for the deep-sea methane seep sediments of the Okhotsk Sea nifH Sampling station

LV39-25H LV39-40H

mcrA

16S rRNA

No. of clones

No. of RFLP phylotypes

No. of OTUsa

No. of clones


No. of OTUsa

No. of clones


No. of OTUsa

79 79

43 38

26 21

93 90

41 22

20 14

93 79

17 23

13 23

a OTUs were calculated using DOTUR (28) based on 5% sequence distance for the NifH protein sequences, 3% sequence distance for the 16S rRNA gene sequences, and uniqueness of the McrA protein sequences.

Five OTUs were affiliated within cluster III, which could be divided into two subclusters (IIIa and IIIb) (Fig. 1). Four OTUs were affiliated within subcluster IIIa and one within IIIb. Subcluster IIIa comprised sequences mainly from anaerobic bacteria, and IIIb comprised those mainly from Methanosarcina methanogens. OTU 40H-0N-5 shared 100.0% identity with its closest match GenBank NifH sequences (accession no. BAF96834, ACD50919, and ACD50920) obtained from the deep-sea methane seep sediments of the Nankai Trough and Eel River Basin (25). OTUs 40H-0N-5 and 40H-0N-6 also shared 86.0 to 86.8% identities to the NifH sequence of Desulfovibrio gigas. Thus, both OTUs were putatively obtained from deep-sea sulfate-reducing bacteria. OTU 40H0N-1 shared 89.1 to 89.9% identities with its closest match NifH sequences of Methanosarcina barkeri strain Fusaro, Methanosarcina acetivorans C2A, and Methanosarcina mazei Go1. This OTU was putatively obtained from a methanogen. Five NifH OTUs were affiliated within the putative cluster III-x, and they shared very high sequence identities (98.3 to 99.2%) with some sequences obtained from the methane seep sediments of the Nankai Trough and Eel River Basin (25) (Fig. 1). However, all of the sequences in this cluster had low identity (⬍79.7%) to the GenBank sequences obtained from other environments or microorganisms, indicating that these NifH sequences might be specific to the deep-sea methane seep sedimentary environments. This new NifH cluster is monophyletic and supported by bootstrapping statistics. Due to the lack of NifH sequences from culturable strains, no conclusive taxonomic information can be obtained for this NifH cluster. Within the DNA fragments analyzed, 10 highly conserved key amino acid residues previously identified to be potentially important in NifH structure and function were examined, including Lys15, Ser16, Arg100, Thr104, Asp125, Asp129, and four Cys residues (no. 38, 85, 97, and 132) (18, 27). All of these key amino acid residues were conserved in our cluster III-x NifH sequences, except for 25H-0N-12 and 25H-0N-9, each with a single-amino-acid residue substitution (see Fig. S2 in the supplemental material). The newly defined cluster III-x sequences are likely to encode functional nitrogenase reductases. Seven OTUs were affiliated within cluster IV. OTU 40H0N-24 shared 95.9 to 100.0% identities with some NifH sequences obtained from the methane seep sediments of the Nankai Trough and Eel River Basin (GenBank accession no. ACD50923 and BAF96826) (25). OTU 25H-0N-3 shared 97.3% identity with a Nankai Trough methane seep NifH sequence (GenBank accession no. BAF96825) and 87.1% identity with the NifH sequence of Methanococcoides burtonii (Fig. 1). Both OTUs were probably obtained from deep-sea archaea. However, it is difficult to determine the exact taxonomic

affiliations of most of our cluster IV NifH sequences, especially for the subcluster made exclusively by OTUs 25H0N-15, 40H-0N-33, and 40H-0N-37. Cluster IV NifH-like sequences are polyphyletic, including many atypical methanogen NifH-like sequences without detectable nitrogen fixation activity (23, 24, 26). This “cluster” is actually a sequence assemblage. Currently most of these sequences are not characterized (25, 26). Thirteen OTUs were affiliated within the putative cluster V, which comprised novel NifH-like sequences detected only from the methane seep sediments of the Okhotsk Sea. These sequences shared quite low identity (47.1 to 59.6%) with their closest match sequences, the hyperthermophilic methanogen FS406-22 NifH2 (GenBank accession no. ABK78685) (22), the Geobacter sp. strain FRC-32 NifH (GenBank accession no. EAT62753), and the Desulfitobacterium hafniense Y51 putative NifH (GenBank accession no. BAE85651). The detection of these sequences in both stations indicates that these sequences probably were not artifacts. Most of the 10 highly conserved key amino acid residues examined (18, 27) were present in these sequences, except for residues Cys85, Arg100, and Thr104. OTUs 25H-0N-39 and 25H-0N-40 had a replacement of Cys85 with Gly85. Substitutions at this position were also found in other environmental NifH sequences (18). All of the cluster V sequences had a double mutation: Arg100 was replaced by Lys100 and Thr104 by Glu104 (see Fig. S2 in the supplemental material). These two residues are located at the interaction interface between the Fe-protein and the MoFe-protein (27). A single replacement of Arg100 with Lys100 was found to render a loss of the NifH activity in Azotobacter vinelandii (30), while in the NifH-like light-independent protochlorophyllide reductase subunit BchL or ChlL, Arg100 was replaced by Tyr100 (4). Whether the double mutation would render any change to the enzyme activity or function is currently unknown. The cluster V NifH-like sequences also contained a 6-amino-acid insertion at position 78 (see Fig. S2 in the supplemental material). Previous studies also found insertions at this position (23). Molecular modeling with SwissModel (1) indicated that this insertion formed a loop on the exterior surface of the NifH homodimer structure, remote from the subunit-subunit interface (see Fig. S3 in the supplemental material) (4). For the NifH-like BchL, a recent study indicates that essentially no conservation of amino acids is observed on the surface of the BchL monomer outside of the BchL-BchL interface (3). BchL is not only a homolog but also a functional analog to NifH (26). The insertion at position 78 on the outside surface might not influence the NifH-like protein activity or function, if any ever exists for the cluster V sequences.

2240

DANG ET AL.

APPL. ENVIRON. MICROBIOL.

VOL. 75, 2009


Three OTUs were affiliated within the putative cluster VI (Fig. 1). However, the low bootstrap support (14%) indicates the uncertainty about the phylogenetic position of this cluster. These OTUs shared low identity (⬍56.1%) to known NifH or NifH-like sequences. The occurrence of these sequences in both sampling stations indicated that they were probably not artifacts. The preservation of the most conserved key amino acid residues (18, 27), including Lys15, Ser16, Arg100, Asp125, Asp129, and four Cys residues (no. 38, 85, 97, and 132) (see Fig. S2 in the supplemental material), suggests that these sequences may encode functional nitrogenase reductases. Some studies indicate that the loss of nitrogen fixation capability in some of the related microorganisms, such as Fusobacterium nucleatum, may be caused by the lack of a complete nif operon (26), probably not by the nifH-encoded enzyme per se. However, the phylogenetic position and the functionality and activity of the “cluster VI” NifH-like proteins need to be determined. Environmental classification of the nifH-carrying microbial communities. UniFrac principal coordinate analysis and environment clustering analysis (19) of the nifH-carrying microbial assemblages in various marine environments indicated that environmental characteristics might have a strong influence on the composition and structure of the nifH-carrying microbial community (Fig. 2; and see Fig. S4 in the supplemental material). The nifH-carrying microbial assemblages of the deep-sea methane seep sediment environments from the Okhotsk Sea, Nankai Trough, and Eel River Basin (25) were distinct from those of the marine water and hydrothermal vent environments (21, 23, 24). The occurrence of the unique cluster V and the putative cluster VI NifH-like sequences in the Okhotsk Sea indicated the environmental specificity of these novel sequences and differentiated the nitrogen-fixing community of the Okhotsk Sea from those of the Nankai Trough and Eel River Basin methane seep environments (Fig. 2; and see Fig. S4 in the supplemental material). Archaeal 16S rRNA gene diversity. A total of 63 unique archaeal 16S rRNA gene sequences and 27 OTUs (based on 3% sequence distance cutoff) were obtained from the methane seep sediments of the Okhotsk Sea (Table 1; and see Fig. S1 in the supplemental material). These sequences fall into diverse archaeal clusters, including the uncultured marine group I; marine benthic groups B, C, and D; marine hydrothermal vent group; miscellaneous crenarchaeotic group; novel group of crenarchaea; terrestrial miscellaneous euryarchaeotic group; and ANME-2 (Fig. 3). The diversity of archaeal 16S rRNA gene sequences from the Okhotsk Sea is comparable to or higher than those of some of other deep-sea methane seeps or like environments (13, 14, 17).

2241

FIG. 2. Principal coordinate plot of UniFrac principal coordinate analysis of the nifH-containing microbial assemblages from various marine environments where nifH clone libraries were constructed with the PCR primers of Mehta et al. (23). 25H and 40H, Okhotsk Sea deep-sea methane seep sediments of stations LV39-25H and LV3940H (this study); NT, Nankai Trough deep-sea anoxic methane seep sediments (GenBank accession no. BAF96788, BAF96793, BAF96798, BAF96803, BAF96808, and BAF96824 to BAF96838); ERB, Eel River Basin deep-sea methane seep sediments (25); HV20, hydrothermal vent fluid of Juan de Fuca Ridge sampled in 2000; HV99, hydrothermal vent fluid of Juan de Fuca Ridge sampled in 1999; DVP, diffuse vent fluid near Puffer on Endeavor segment of the Juan de Fuca Ridge; DSW, deep seawater on Endeavor segment (23); NPD, Northeast Pacific deep seawater (24); H01, H04 and H05, surface seawater samples collected in stations H01, H04, and H05 of the eastern Mediterranean Sea (21).

Archaeal mcrA gene diversity. A total of 33 unique mcrA gene sequences and 28 McrA OTUs were obtained based on sequence uniqueness from the methane seep sediments of the Okhotsk Sea (Table 1; and see Fig. S1 in the supplemental material). Diverse tentative ANME mcrA gene sequences, including groups a, c, d, and e (12), were detected via phylogenetic analysis (Fig. 4). No group b or f mcrA sequence could be obtained with the primers used (11, 12, 17). The mcrA analysis results (Fig. 4) detected much higher ANME diversity of both sampling stations LV39-25 and LV39-40H than the archaeal 16S rRNA gene analysis results (Fig. 3). For the sampling

FIG. 1. Phylogenetic analysis of the NifH and NifH-like sequences obtained from the methane seep sediments of the Okhotsk Sea. Sequences were aligned using the Clustal_X program (29), and the phylogenetic tree was constructed with the amino acid residues corresponding to positions 17 to 130 of the Klebsiella pneumoniae NifH sequence (GenBank DNA accession no. J01740) using the software package Phylip via the distance and neighbor-joining methods (10). The tree branch distances represent the amino acid substitution rate, and the scale bar represents the expected number of changes per homologous position. Bootstrap values (no less than 50%) of 100 resamplings are shown near nodes. The Rhodobacter capsulatus chlorophyllide reductase iron protein subunit BchX sequence was used as an outgroup. The NifH and NifH-like sequences obtained in this study are shown in bold. *, the NifH sequences Bca2c_nif1g, Bca2c_nif3h, and Bca2c_nif9h recovered from the isolated ANME-2 consortia of the Eel River Basin deep-sea methane seep sediments (25) are highlighted in bold. ERB, Eel River Basin deep-sea methane seep sediments; NT, Nankai Trough deep-sea anoxic methane seep sediments.

FIG. 3. Phylogenetic analysis of the archaeal 16S rRNA gene sequences obtained from the deep-sea methane seep sediments of the Okhotsk Sea. Sequences were aligned using the Clustal_X program, and the phylogenetic tree was constructed with 964 aligned nucleotide positions using the software package Phylip via the distance and neighbor-joining methods. The tree branch distances represent nucleotide substitution rate, and the scale bar represents the expected number of changes per homologous position. Bootstrap values no less than 70% of 100 resamplings are shown with solid circles on the corresponding nodes, and those less than 70% but greater than or equal to 50% are shown with open circles on the corresponding nodes. The bacterial 16S rRNA gene sequence from Thermotoga maritima MSB8 was used as an outgroup. The archaeal 16S rRNA gene sequences obtained in this study are shown in boldface. DSAG, deep-sea archaeal group; MBG, marine benthic group; MCG, miscellaneous crenarchaeotic group; MG I, marine group I; MHVG, marine hydrothermal vent group; NGC, novel group of crenarchaea; TMEG, terrestrial miscellaneous euryarchaeotic group. 2242

VOL. 75, 2009


2243

FIG. 4. Phylogenetic analysis of the archaeal McrA sequences obtained from the methane seep sediments of the Okhotsk Sea. Sequences were aligned using the Clustal_X program, and the phylogenetic tree was constructed with 244 aligned amino acid residues using the software package Phylip via the distance and neighbor-joining methods. The tree branch distances represent the amino acid substitution rate, and the scale bar represents the expected number of changes per homologous position. Bootstrap values no less than 70% of 100 resamplings are shown with solid circles on the corresponding nodes and those less than 70% but greater than or equal to 50% are shown with open circles on the corresponding nodes. The McrA sequence from Methanopyrus kandleri AV19 was used as an outgroup. McrA sequences obtained in this study are shown in boldface.

2244

DANG ET AL.

station LV39-25H, a methanogen mcrA sequence (25H-0A-8) affiliated within Methanosarcinales was also obtained (Fig. 4). However, for both mcrA clone libraries, sequences putatively affiliated within ANME were much more predominant than those affiliated within methanogens. In both stations LV39-25H and LV39-40H, methane hydrates were directly obtained from shallow sediments, usually less than 2 m below the seafloor. The existence of intense methane plumes above the seafloor indicated a sufficient in situ CH4 supply for the sedimentary microbial ecosystem and a potential oxygen-limited environment of the surface sediments (20). ANME might be a major functional component of the in situ microbial community. Previous studies indicate that cultivated archaeal nitrogen fixers are all affiliated within methanogenic Euryarchaeota, including Methanosarcinales, Methanococcales, Methanomicrobiales, and Methanobacteriales (16). However, recent studies also indicate that some of the deep-sea novel nifH and nifH-like sequences might be obtained from ANME (25). In line with this, most of our nifH sequences had quite high identity with sequences obtained from anaerobic methanotrophic consortia of other deep-sea methane seep sedimentary environments, such as the Nankai Trough and Eel River Basin (25). The recent GenBank publications of several nearly complete nif operons putatively affiliated with anaerobic methane-oxidizing archaea from the Nankai Trough methane seep sediments (GenBank DNA accession no. AB362194, AB362195, and AB362197; corresponding NifH protein accession no. BAF96793, BAF96798, and BAF96808) indicate that some methanotrophic archaea may possess the genetic potential for nitrogen fixation (Ken Takai, personal communication). This was further demonstrated by several NifH sequences, namely Bca2c_nif1g, Bca2c_nif3h, and Bca2c_nif9h, recently obtained from the isolated ANME-2 consortia of the Eel River Basin methane seep sediments (25). All of these NifH sequences show a high degree of identity to the cluster III-x sequences obtained from our study (Fig. 1). In Pernthaler and colleagues’ study (25), in situ nitrogen fixation by methane seep sediment methanotrophic consortia was demonstrated via a 15N2 incubation study, which also showed that ANME might be the major contributors to the diazotrophic activity of the consortia. The ANME microorganisms still escape from being cultivated for physiological and genetic research. Cultivation-independent molecular approaches have significantly advanced our understanding of their distribution and ecophysiology in natural and engineering environments. In bacterial methanotrophs, the genetic potential and nitrogen fixation activity have already been identified (2). It is reasonable to speculate that the genetic potential and diazotrophic activity might also exist very commonly in anaerobic archaeal methanotrophs, especially in deep-sea methane-rich sedimentary environments. Recent studies also indicate the possibility that some of the deep-sea novel nifH sequences might be obtained from other archaeal or proteobacterial lineages (24, 25). Further investigations should be carried out to test these hypotheses. In conclusion, our work obtained diverse and novel nifH and nifH-like gene sequences in the deep-sea methane seep sediments of the Okhotsk Sea. Almost all of the sequences were obtained from uncultured or uncharacterized bacteria or archaea. Comparison with other marine environments indicates certain environmental specificity of some of the sequences. Methane seep sediments may harbor diverse and unique mi-

APPL. ENVIRON. MICROBIOL.

crobial genetic potentials, important in nitrogen fixation and new production of the deep-sea ecosystem. Nucleotide sequence accession numbers. The nifH gene sequences reported in this study have been deposited in GenBank under accession no. EU713922 to EU713990 and EU713992 to EU714003, the archaeal 16S rRNA gene sequences under accession no. EU713859 to EU713921, and the mcrA gene sequences under accession no. FJ403593 to FJ403625. We thank the reviewers for valuable comments and suggestions. We also thank Linbao Zhang, Ying Zhang, Jian Sun, and Jin Sun for assistance with the project. This work was supported by China National Natural Science Foundation grants 40576069 and 40776032, National Basic Research Program of China (973 Program) grant 2007CB411702, Hi-Tech Research and Development Program of China grant 2007AA091903, and China Ocean Mineral Resources R&D Association grants DYXM-115-022-6 and DYXM-115-02-2-20. REFERENCES 1. Arnold, K., L. Bordoli, J. Kopp, and T. Schwede. 2006. The SWISS-MODEL workspace: a web-based environment for protein structure homology modelling. Bioinformatics 22:195–201. 2. Auman, A. J., C. C. Speake, and M. E. Lidstrom. 2001. nifH sequences and nitrogen fixation in type I and type II methanotrophs. Appl. Environ. Microbiol. 67:4009–4016. 3. Bro ¨cker, M. J., S. Virus, S. Ganskow, P. Heathcote, D. W. Heinz, W.-D. Schubert, D. Jahn, and J. Moser. 2008. ATP-driven reduction by darkoperative protochlorophyllide oxidoreductase from Chlorobium tepidum mechanistically resembles nitrogenase catalysis. J. Biol. Chem. 283:10559– 10567. 4. Burke, D. H., J. E. Hearst, and A. Sidow. 1993. Early evolution of photosynthesis: clues from nitrogenase and chlorophyll iron proteins. Proc. Natl. Acad. Sci. USA 90:7134–7138. 5. Cabello, P., M. D. Rolda ´n, and C. Moreno-Vivia ´n. 2004. Nitrate reduction and the nitrogen cycle in archaea. Microbiology 150:3527–3546. 6. Dang, H., and C. R. Lovell. 2000. Bacterial primary colonization and early succession on surfaces in marine waters as determined by amplified rRNA gene restriction analysis and sequence analysis of 16S rRNA genes. Appl. Environ. Microbiol. 66:467–475. 7. Dang, H., T. Li, M. Chen, and G. Huang. 2008. Cross-ocean distribution of Rhodobacterales bacteria as primary surface colonizers in temperate coastal marine waters. Appl. Environ. Microbiol. 74:52–60. 8. Dang, H. Y., X. X. Zhang, J. Sun, T. G. Li, Z. N. Zhang, and G. P. Yang. 2008. Diversity and spatial distribution of sediment ammonia-oxidizing crenarchaeota in response to estuarine and environmental gradients in the Changjiang Estuary and East China Sea. Microbiology 154:2084–2095. 9. DeLong, E. F. 1992. Archaea in coastal marine environments. Proc. Natl. Acad. Sci. USA 89:5685–5689. 10. Felsenstein, J. 1989. PHYLIP—Phylogeny Inference Package (version 3.2). Cladistics 5:164–166. 11. Hales, B. A., C. Edwards, D. A. Ritchie, G. Hall, R. W. Pickup, and J. R. Saunders. 1996. Isolation and identification of methanogen-specific DNA from blanket bog peat by PCR amplification and sequence analysis. Appl. Environ. Microbiol. 62:668–675. 12. Hallam, S. J., P. R. Girguis, C. M. Preston, P. M. Richardson, and E. F. DeLong. 2003. Identification of methyl coenzyme M reductase A (mcrA) genes associated with methane-oxidizing archaea. Appl. Environ. Microbiol. 69:5483–5491. 13. Heijs, S. K., R. R. Haese, P. W. van der Wielen, L. J. Forney, and J. D. van Elsas. 2007. Use of 16S rRNA gene based clone libraries to assess microbial communities potentially involved in anaerobic methane oxidation in a Mediterranean cold seep. Microb. Ecol. 53:384–398. 14. Inagaki, F., T. Nunoura, S. Nakagawa, A. Teske, M. Lever, A. Lauer, M. Suzuki, K. Takai, M. Delwiche, F. S. Colwell, K. H. Nealson, K. Horikoshi, S. D’Hondt, and B. B. Jørgensen. 2006. Biogeographical distribution and diversity of microbes in methane hydrate-bearing deep marine sediments on the Pacific Ocean Margin. Proc. Natl. Acad. Sci. USA 103:2815–2820. 15. Joye, S. B., A. Boetius, B. N. Orcutt, J. P. Montoya, H. N. Schulz, M. J. Erickson, and S. K. Lugo. 2004. The anaerobic oxidation of methane and sulfate reduction in sediments from Gulf of Mexico cold seeps. Chem. Geol. 205:219–238. 16. Leigh, J. A. 2000. Nitrogen fixation in methanogens: the archaeal perspective. Curr. Issues Mol. Biol. 2:125–131. 17. Lo ¨sekann, T., K. Knittel, T. Nadalig, B. Fuchs, H. Niemann, A. Boetius, and R. Amann. 2007. Diversity and abundance of aerobic and anaerobic methane

VOL. 75, 2009

18.

19. 20. 21. 22. 23.

24.


oxidizers at the Haakon Mosby Mud Volcano, Barents Sea. Appl. Environ. Microbiol. 73:3348–3362. Lovell, C. R., M. J. Friez, J. W. Longshore, and C. E. Bagwell. 2001. Recovery and phylogenetic analysis of nifH sequences from diazotrophic bacteria associated with dead aboveground biomass of Spartina alterniflora. Appl. Environ. Microbiol. 67:5308–5314. Lozupone, C. A., M. Hamady, S. T. Kelley, and R. Knight. 2007. Quantitative and qualitative ␤ diversity measures lead to different insights into factors that structure microbial communities. Appl. Environ. Microbiol. 73:1576–1585. Luan, X. W., Y. K. Jin, O. Anatoly, and B. J. Yue. 2008. Characteristics of shallow gas hydrate in Okhotsk Sea. Sci. China Ser. D Earth Sci. 51:415–421. Man-Aharonovich, D., N. Kress, E. B. Zeev, I. Berman-Frank, and O. Be´ja `. 2007. Molecular ecology of nifH genes and transcripts in the eastern Mediterranean Sea. Environ. Microbiol. 9:2354–2363. Mehta, M. P., and J. A. Baross. 2006. Nitrogen fixation at 92°C by a hydrothermal vent archaeon. Science 314:1783–1786. Mehta, M. P., D. A. Butterfield, and J. A. Baross. 2003. Phylogenetic diversity of nitrogenase (nifH) genes in deep-sea and hydrothermal vent environments of the Juan de Fuca Ridge. Appl. Environ. Microbiol. 69:960–970. Mehta, M. P., J. A. Huber, and J. A. Baross. 2005. Incidence of novel and

25.

26. 27. 28. 29.

30.

2245

potentially archaeal nitrogenase genes in the deep Northeast Pacific Ocean. Environ. Microbiol. 7:1525–1534. Pernthaler, A., A. E. Dekas, C. T. Brown, S. K. Goffredi, T. Embaye, and V. J. Orphan. 2008. Diverse syntrophic partnerships from deep-sea methane vents revealed by direct cell capture and metagenomics. Proc. Natl. Acad. Sci. USA 105:7052–7057. Raymond, J., J. L. Siefert, C. R. Staples, and R. E. Blankenship. 2004. The natural history of nitrogen fixation. Mol. Biol. Evol. 21:541–554. Schindelin, H., C. Kisker, J. Schlessman, J. Howard, and D. Rees. 1997. Structure of ADP ⫻ AIF4(⫺)-stabilized nitrogenase complex and its implications for signal transduction. Nature 387:370–376. Schloss, P. D., and J. Handelsman. 2005. Introducing DOTUR, a computer program for defining operational taxonomic units and estimating species richness. Appl. Environ. Microbiol. 71:1501–1506. Thompson, J. D., T. J. Gibson, F. Plewniak, F. Jeanmougin, and D. G. Higgins. 1997. The CLUSTAL_X windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucleic Acids Res. 25:4876–4882. Wolle, D., C. Kim, D. Dean, and J. B. Howard. 1992. Ionic interactions in the nitrogenase complex. Properties of Fe-protein containing substitutions for Arg-100. J. Biol. Chem. 267:3667–3673.