The major sinks for methane are biological oxidation at or near the sites of production (∼700 Tg · year
−1), uptake of methane from the atmosphere by aerobic soils (20 to 60 Tg · year
−1), and photochemical oxidation in the atmosphere (∼450 Tg · year
−1) (
27). Soil uptake of atmospheric methane is significant since the magnitude of the soil sink is equivalent to the observed annual increase in the methane concentration and it is more susceptible to disturbance by human activities (
16,
21,
24,
34). A change in the soil sink can have a significant effect on the atmospheric mixing ratios of methane.
Biological methane oxidation consists of both aerobic and anaerobic processes. The global methane sink is dominated by aerobic methane-oxidizing bacteria (MOB). The biochemical basis of methane oxidation in all known MOB is similar (
1,
9,
13,
20). All MOB possess a membrane-bound monooxygenase whose substrate range includes both methane and ammonia (note that some MOB contain an additional, biochemically distinct enzyme designated the soluble methane monooxygenase [sMMO]) (
1). The membrane-bound monooxygenases are thought to be evolutionarily related (
15). The MOB exhibit limited physiological, structural, and phyletic diversity compared to other functionally defined groups of bacteria (
13,
25). Of particular significance are differences in the fate of carbon, the kinetic properties of the monooxygenase, and the evolutionary separation of the four major phyletic groups.
On the basis of cell physiology, the MOB can be divided into the methane-assimilating bacteria (MAB) (methanotrophs) and bacteria which cooxidize methane (autotrophic ammonia-oxidizing bacteria [AAOB]). The former organisms use methane as a sole source of carbon and energy and are characterized by the presence of a complete pathway for methane oxidation, the ability to assimilate cell carbon as formaldehyde, and apparent
Km values for methane in the micromolar range (
1,
13). The AAOB use ammonia oxidation as an energy source for autotrophic growth; they are characterized by a complete pathway for oxidation of ammonia to nitrite and assimilation of cell carbon by the Benson-Calvin cycle. In most cases their apparent
Km values for methane are in the millimolar range and methane is cooxidized with no apparent benefit to the cells (
1).
Soil methane uptake has been demonstrated to be biological. Methane uptake activity shares many features with the known MOB activity but also exhibits traits which do not occur during methane oxidation by extant organisms. The differences include a >100-fold-greater affinity for methane but an apparently poor capacity for growth on this substrate (
2,
9,
20,
21,
31,
35). Perhaps the most significant difference is a much lower threshold concentration for sustained methane uptake. Several explanations have been proposed to account for this, including (i) mixotrophic growth of methanotrophs, (ii) cooxidation of methane by ammonia oxidizers, (iii) induction of a high-affinity enzyme system in response to starvation, and (iv) activity of novel methanotrophic bacteria (
2,
7,
9,
20,
31,
40). More recently, workers have shown that methanotrophs in pure cultures can exhibit sustained uptake of atmospheric methane at normal atmospheric concentrations if the cultures are supplemented with methanol (
3,
18). These workers proposed that the presence of methanol in soil may provide a physiological basis for methane uptake by conventional methanotrophs in soil.
RESULTS
Methane oxidation.
In all of the soils examined, the greatest atmospheric methane oxidation activity was found in the mineral soil below the organic horizon (depth, <4 cm). The potential activities for oxidation of atmospheric methane were 229, 752, and 21 pmol · g−1 · h−1 for the soils from the Rold, Maine, and Pantanal forests, respectively. The apparentKm for methane oxidation was determined for the Rold forest soil. Fresh soil exhibited an apparentKm of 10.7 ppm of CH4, which corresponded to approximately 14.7 nM CH4 in the soil water.
PLFA analysis.
Soil samples were incubated with
14CH
4 (<50 ppm) in order to radiolabel the microorganisms that metabolize methane at near-atmospheric concentrations. Phospholipids extracted from soil samples incubated with
14CH
4 were assayed for the presence of
14C-PLFAs. Figure
1A shows the relative levels of
14C-PLFAs in radiolabelled samples obtained from the Rold and Maine forest soils. Significant amounts of radioactivity were detected in 3 of the 15 fractions (fractions 5, 7, and 11). These PLFA fractions represent fatty acids with equivalent chain lengths of 15.7 to 16.1 (fraction 5), 16.5 to 16.8 (fraction 7), and 17.7 to 17.9 (fraction 11) (see references
29 and 30 for details).
The inferred PLFA signature for the soil methanotrophic community under these labelling conditions was different from the PLFA signature reported for known MOB, including autotrophic ammonia oxidizers and type I or type II methanotrophic bacteria (
12,
30,
41). The distinction from methanotrophs was confirmed by the results of a
14C-PLFA analysis of a representative type I methanotroph (
Methylomicrobium album BG8) and a type II methanotroph (
Methylosinus trichosporium OB3b) labelled under comparable conditions (Fig.
1B). The
14C-PLFA(s) present in fraction 7 (Fig.
1A) had a retention time that was significantly different from the retention time of any fatty acid labelled in the control methanotrophs (the values were >5 min different). The fatty acids in fraction 7 coeluted with the reference standard fatty acids i17:0, a17:0, and 17:1w8.
pmoA analysis.
The presence of uncharacterized MOB in soil was confirmed by using a molecular biology approach. The Rold library contained 14 clones of methane monooxygenase (MMO) type sequences, and 4 of these clones (designated RA9, RA14, RA21, and RA26) could not be placed in any known group of MMO or AMO sequences in our database. The remaining clones were identified asMethylococcus sp. (four clones), Nitrosospira sp. (three clones), Nitrosococcus sp. (one clone), orNitrosomonas sp. (two clones). Clones RA9, RA14, and RA26 were very similar, exhibiting 98% identity.
pmoA libraries were also constructed by using three soils obtained from the Rold, Maine, and Pantanal forests. Forty clones from each library were screened by PCR with RA14- and RA21-specific PCR primers. A total of 15 clones were identified and sequenced. Three clones each from the Rold, Maine, and Pantanal libraries were similar to RA14 (Rold 1, Rold 3, and Rold 5; Maine 6, Maine 8, and Maine 10; and Pantanal 11, Pantanal 13, and Pantanal 14), and a tenth clone which was amplified with the RA14 primers, Rold 4, was similar to thepmoA sequence of Methylocystis sp. strain M. However, the five clones identified with the RA21 primers were all similar to the Nitrosomonas amoA sequences.
Alignments of available PmoA and AmoA sequence information (Fig.
2) revealed 52 universally conserved amino acid residues. Members of the RA14 group possessed 50 of these 52 universally conserved residues, and clone RA21 had 45 matches with these 52 conserved amino acid residues. While the new sequences clearly belong to the same protein family, they are significantly different at sites which are otherwise highly conserved in physiologically and phyletically diverse organisms.
Sequence analysis.
Dendrograms constructed from the inferred peptide sequences of PmoA and AmoA showed that there were three strongly supported monophyletic groups: (i) the representatives of the α subclass of the class
Proteobacteria(α-
Proteobacteria) (type II methanotrophs) and RA14, (ii) the γ-
Proteobacteria representatives (type I methanotrophs and
Nitrosococcus oceanus), and (iii) the β-
Proteobacteria ammonia oxidizers (Fig.
3). The relationships among the cultivated strains were broadly similar to the relationships based on an analysis of several other genes, including the 16S rRNA gene (
26). We did not observe a significant relationship between RA21 and any of the other PmoA or AmoA sequences.
Phylogenetic analysis of the new Rold, Maine, and Pantanal PmoA sequences (Fig.
4) revealed a number of clones from each soil which grouped with the RA14 sequence, but we detected no clones which grouped with the RA21 sequence.
DISCUSSION
The MOB are a physiologically (functionally) defined group of bacteria. Four subgroups of MOB have been recognized; there are two groups of MAB (type I and II methanotrophs) and two groups of AAOB (β- and γ-subdivision ammonia oxidizers). Soil methane uptake activity experiments typically reveal biochemical characteristics with broad similarity to characteristics of the known members of these four phyletic groups (Table
1). The debate concerning the nature of the organisms responsible for soil methane uptake has focused on three distinct issues. Do all soil methane oxidizers possess the same biochemistry for methane oxidation? What is the physiological role of methane in soil methane oxidizers? Do soil methane oxidizers belong to any of the four currently known phyletic groups of MOB?
To determine how diverse the soil methane oxidizers are, we analyzed the distribution of radiolabelled PLFAs in soil samples after incubation with
14CH
4. This analysis was based on the assumption that atmospheric methane is assimilated into biomass by the same population of organisms that oxidize it. This has been shown previously to be the most likely scenario for methane oxidizers in a forest soil (
28). In the Rold, Maine, and Pantanal forest soils, 36, 39, and 35%, respectively, of the atmospheric methane consumed was assimilated into microbial biomass (
30). This relatively high level of carbon conversion efficiency is consistent with the hypothesis that soil methane uptake is dominated by MAB. The relatively simple
14C-PLFA profile obtained from soils supplied with low methane concentrations (Fig.
1) suggests that the diversity of the organisms that incorporate the label is low (i.e., all members of the population are likely to be related at the generic level). In broad terms the PLFA profile of the soil methanotroph community was comparable to the PLFA profiles of members of the α-
Proteobacteria, including the type II methanotrophs (Fig.
1). However, of particular significance is the peak associated with fraction 7, which represents 8 to 9% of the total. The labelled PLFAs in this fraction coeluted with reference standards i17:0, a17:0, and 17:1w8, which are scarce (<1%) or absent in previously studied methanotrophs and autotrophic ammonia oxidizers (
5,
12,
30,
41). We have recently obtained similar
14C-PLFA profiles for high-affinity methane oxidizers in other soils, including soils from Brazil, the United States, Denmark, and Greenland (
30).
Some bacteria are known to produce culture-dependent PLFA profiles, and it could be argued that the differences between the profiles obtained for the type II methanotrophs and the soil methane oxidizers (principally the fraction 7 fatty acid) reflect metabolic stresses on cells in the soil. However, soils labelled in the presence of saturating concentrations of methane (10,000 ppm) still contained significant amounts of these PLFAs, which could not be detected in pure cultures of methanotrophs labelled at low methane concentrations (
30).
Soil samples that actively oxidized atmospheric methane contained phyletically distinct pMMO-like enzymes. On the basis of the relationships between PmoA sequences, clone RA14 may be considered a deep-branching member of the type II methanotroph group. Clone RA21 is not significantly related to any known pMMO-like sequence (Fig.
3 and
4). While RA14 is strongly related to type II methanotroph PmoA sequences, the evolutionary distance of this relationship is relatively large compared to the distance between known type II methanotroph PmoA sequences. On this basis, while it is probable that the RA14 organism is also a member of the α-
Proteobacteria, it is unlikely to be sufficiently closely related to type II methanotrophs to form a monophyletic group. Other RA14 type sequences were detected from Maine, Pantanal, and additional Rold soil samples that had
14C-PLFA profiles similar to that of the Rold forest soil (Fig.
4). These sequences formed a cluster with the RA14 sequence, indicating that these types of sequences are found in a range of soils that oxidize atmospheric methane.
The topology of the monooxygenase trees appears to primarily reflect the phylogenetic relationships of the organisms and not the physiological roles of the enzymes (viz., the
N. oceanusAmoA clusters with PmoA sequences). However, a comparison of the monooxygenase tree with trees derived from 16S rRNA (
26,
42) or 23S rRNA (
10) did indicate that the monooxygenase tree may reflect some specialization for using either methane or ammonia as a substrate for the enzyme (
26). Phylogenies inferred from rRNAs indicate that the β- and γ-
Proteobacteria are sister groups, whereas the monooxygenase tree (Fig.
3A) revealed that there is a comparatively close relationship between the α- and γ-
Proteobacteria and that the β-
Proteobacteria appear to be a highly divergent group. This could reflect selection for ammonia-specialized enzymes in the β-
Proteobacteria and/or convergence to methane-specialized monooxygenases in the α- and γ-
Proteobacteria. This implies that some sites in the protein sequence alignments are subject to positive selection, which results in bias in the monooxygenase tree. This hypothesis is consistent with current data on apparent
Km and
Vmax values for methane for members of the α- and γ-
Proteobacteria(Table
1 and Fig.
3A) and the apparent ability of
N. oceanusto use methane as an alternative carbon source (
1,
19,
39). The relationship of RA14 to the α- and γ-
Proteobacteriasuggests that this clone represents a methane-specialized form. In the absence of pure cultures, confirmation of this will require demonstration that a radiotracer is incorporated into the RA14
pmoA gene or gene product. We have been unable to recover sufficient labelled genetic material from soil to demonstrate this to date.
Identification of amino acid residues which indicate that an organism is adapted to methane as the preferred substrate (i.e., signature residues for methane-specialized enzymes) may provide circumstantial evidence which supports the positive selection hypothesis. An amino acid was considered a putative signature residue if it satisfied the following criteria: it had to be universally conserved in all known members of the α- and γ-
Proteobacteria, and at the same position in the alignment there had to be universal conservation of a different residue in all known members of the β-
Proteobacteria. Sites fulfilling these criteria were found to be concentrated in two regions of the peptides. Interestingly, these two regions lie outside potential membrane-spanning domains (Fig.
2) and correlate closely with domains of the AMO which have been predicted to form the active site of the enzyme (
36). The RA14 cluster sequences contained 74% (16 of 21) of the putative MMO signature residues and only 10% (2 of 21) of the putative AMO signatures. In contrast, RA21 exhibited no clear bias toward either substrate; it contained 33% (7 of 21) of the putative MMO signatures and 24% (5 of 21) of the putative AMO signatures. The effect of this putative positive selection on tree topology was investigated by constructing trees for which all such residues were omitted from the analysis. The resulting tree (Fig.
3B) still strongly supported the hypothesis that RA14 is related to the α-
Proteobacteria(98% of bootstrap replicates), and RA21 branched at the predicted position for the root of the tree.
The existing biomarkers used to study MOB ecology include signature lipids for both type I and type II methanotrophs (
12,
13), as well as several phylogenetic group-specific oligonucleotide probes that target the 16S rRNA (
6,
14,
23,
37,
38). The evolutionary distances of the RA14 and RA21 PmoA sequences suggest that previously described phylogenetic group-specific biomarkers are not likely to detect these new groups of putative MOB.
In conclusion, using radiotracers and functional gene probes, we obtained substantial evidence that there are phyletically distinct populations of MOB in soils. The results of a comparative analysis of two distinct biomarkers (PLFA and PmoA) suggested that the high-affinity methane oxidizers are a novel group of α-Proteobacteria methanotrophs. Studies on factors that affect the activity, population dynamics, and contribution to global methane flux of atmospheric methane oxidizers should be facilitated by use of these markers.