Search

Annually, half of all plant-derived carbon is added to soil where it is microbially respired to CO2. However, understanding of the microbiology of this process is limited because most culture-independent methods cannot link metabolic processes to the organisms present, and this link to causative agents is necessary to predict the results of perturbations on the system. We collected soil samples at two sub-root depths (10-20 cm and 30-40 cm) before and after a rainfall-driven nutrient perturbation event in a Northern California grassland that experiences a Mediterranean climate. From ten samples, we reconstructed 198 metagenome-assembled genomes that represent all major phylotypes. We also quantified 6,835 proteins and 175 metabolites and showed that after the rain event the concentrations of many sugars and amino acids approach zero at the base of the soil profile. Unexpectedly, the genomes of novel members of the Gemmatimonadetes and Candidate Phylum Rokubacteria phyla encode pathways for methylotrophy. We infer that these abundant organisms contribute substantially to carbon turnover in the soil, given that methylotrophy proteins were among the most abundant proteins in the proteome. Previously undescribed Bathyarchaeota and Thermoplasmatales archaea are abundant in deeper soil horizons and are inferred to contribute appreciably to aromatic amino acid degradation. Many of the other bacteria appear to breakdown other components of plant biomass, as evidenced by the prevalence of various sugar and amino acid transporters and corresponding hydrolyzing machinery in the proteome. Overall, our work provides organism-resolved insight into the spatial distribution of bacteria and archaea whose activities combine to degrade plant-derived organics, limiting the transport of methanol, amino acids and sugars into underlying weathered rock. The new insights into the soil carbon cycle during an intense period of carbon turnover, including biogeochemical roles to previ

Annually, half of all plant-derived carbon is added to soil where it is microbially respired to CO2. However, understanding of the microbiology of this process is limited because most culture-independent methods cannot link metabolic processes to the organisms present, and this link to causative agents is necessary to predict the results of perturbations on the system. We collected soil samples at two sub-root depths (10-20 cm and 30-40 cm) before and after a rainfall-driven nutrient perturbation event in a Northern California grassland that experiences a Mediterranean climate. From ten samples, we reconstructed 198 metagenome-assembled genomes that represent all major phylotypes. We also quantified 6,835 proteins and 175 metabolites and showed that after the rain event the concentrations of many sugars and amino acids approach zero at the base of the soil profile. Unexpectedly, the genomes of novel members of the Gemmatimonadetes and Candidate Phylum Rokubacteria phyla encode pathways for methylotrophy. We infer that these abundant organisms contribute substantially to carbon turnover in the soil, given that methylotrophy proteins were among the most abundant proteins in the proteome. Previously undescribed Bathyarchaeota and Thermoplasmatales archaea are abundant in deeper soil horizons and are inferred to contribute appreciably to aromatic amino acid degradation. Many of the other bacteria appear to breakdown other components of plant biomass, as evidenced by the prevalence of various sugar and amino acid transporters and corresponding hydrolyzing machinery in the proteome. Overall, our work provides organism-resolved insight into the spatial distribution of bacteria and archaea whose activities combine to degrade plant-derived organics, limiting the transport of methanol, amino acids and sugars into underlying weathered rock. The new insights into the soil carbon cycle during an intense period of carbon turnover, including biogeochemical roles to previ

In soil ecosystems, microorganisms produce diverse secondary metabolites such as antibiotics, antifungals and siderophores that mediate communication, competition and interactions with other organisms and the environment1,2. Most known antibiotics are derived from a few culturable microbial taxa 3 , and the biosynthetic potential of the vast majority of bacteria in soil has rarely been investigated 4 . Here we reconstruct hundreds of near-complete genomes from grassland soil metagenomes and identify microorganisms from previously understudied phyla that encode diverse polyketide and nonribosomal peptide biosynthetic gene clusters that are divergent from well-studied clusters. These biosynthetic loci are encoded by newly identified members of the Acidobacteria, Verrucomicobia and Gemmatimonadetes, and the candidate phylum Rokubacteria. Bacteria from these groups are highly abundant in soils5-7, but have not previously been genomically linked to secondary metabolite production with confidence. In particular, large numbers of biosynthetic genes were characterized in newly identified members of the Acidobacteria, which is the most abundant bacterial phylum across soil biomes 5 . We identify two acidobacterial genomes from divergent lineages, each of which encodes an unusually large repertoire of biosynthetic genes with up to fifteen large polyketide and nonribosomal peptide biosynthetic loci per genome. To track gene expression of genes encoding polyketide synthases and nonribosomal peptide synthetases in the soil ecosystem that we studied, we sampled 120 time points in a microcosm manipulation experiment and, using metatranscriptomics, found that gene clusters were differentially co-expressed in response to environmental perturbations. Transcriptional co-expression networks for specific organisms associated biosynthetic genes with two-component systems, transcriptional activation, putative antimicrobial resistance and iron regulation, linking metabolite biosynthesis to

In soil ecosystems, microorganisms produce diverse secondary metabolites such as antibiotics, antifungals and siderophores that mediate communication, competition and interactions with other organisms and the environment1,2. Most known antibiotics are derived from a few culturable microbial taxa 3 , and the biosynthetic potential of the vast majority of bacteria in soil has rarely been investigated 4 . Here we reconstruct hundreds of near-complete genomes from grassland soil metagenomes and identify microorganisms from previously understudied phyla that encode diverse polyketide and nonribosomal peptide biosynthetic gene clusters that are divergent from well-studied clusters. These biosynthetic loci are encoded by newly identified members of the Acidobacteria, Verrucomicobia and Gemmatimonadetes, and the candidate phylum Rokubacteria. Bacteria from these groups are highly abundant in soils5-7, but have not previously been genomically linked to secondary metabolite production with confidence. In particular, large numbers of biosynthetic genes were characterized in newly identified members of the Acidobacteria, which is the most abundant bacterial phylum across soil biomes 5 . We identify two acidobacterial genomes from divergent lineages, each of which encodes an unusually large repertoire of biosynthetic genes with up to fifteen large polyketide and nonribosomal peptide biosynthetic loci per genome. To track gene expression of genes encoding polyketide synthases and nonribosomal peptide synthetases in the soil ecosystem that we studied, we sampled 120 time points in a microcosm manipulation experiment and, using metatranscriptomics, found that gene clusters were differentially co-expressed in response to environmental perturbations. Transcriptional co-expression networks for specific organisms associated biosynthetic genes with two-component systems, transcriptional activation, putative antimicrobial resistance and iron regulation, linking metabolite biosynthesis to

Multicopper oxidases (MCOs) catalyze the oxidation of a diverse group of metal ions and organic substrates by successive single-electron transfers to O2 via four bound Cu ions. MnxG, which catalyzes MnO2 mineralization by oxidizing both Mn(II) and Mn(III), is unique among multicopper oxidases in that it carries out two energetically distinct electron transfers and is tightly bound to accessory proteins. There are two of these, MnxE and MnxF, both approximately 12kDa. Although their sequences are similar to those found in the genomes of several Mn-oxidizing Bacillus species, they are dissimilar to those of proteins with known function. Here, MnxE and MnxF are co-expressed independent of MnxG and are found to oligomerize into a higher order stoichiometry, likely a hexamer. They bind copper and heme, which have been characterized by electron paramagnetic resonance (EPR), X-ray absorption spectroscopy (XAS), and UV-visible (UV-vis) spectrophotometry. Cu is found in two distinct type 2 (T2) copper centers, one of which appears to be novel; heme is bound as a low-spin species, implying coordination by two axial ligands. MnxE and MnxF do not oxidize Mn in the absence of MnxG and are the first accessory proteins to be required by an MCO. This may indicate that Cu and heme play roles in electron transfer and/or Cu trafficking.

Multicopper oxidases (MCOs) catalyze the oxidation of a diverse group of metal ions and organic substrates by successive single-electron transfers to O2 via four bound Cu ions. MnxG, which catalyzes MnO2 mineralization by oxidizing both Mn(II) and Mn(III), is unique among multicopper oxidases in that it carries out two energetically distinct electron transfers and is tightly bound to accessory proteins. There are two of these, MnxE and MnxF, both approximately 12kDa. Although their sequences are similar to those found in the genomes of several Mn-oxidizing Bacillus species, they are dissimilar to those of proteins with known function. Here, MnxE and MnxF are co-expressed independent of MnxG and are found to oligomerize into a higher order stoichiometry, likely a hexamer. They bind copper and heme, which have been characterized by electron paramagnetic resonance (EPR), X-ray absorption spectroscopy (XAS), and UV-visible (UV-vis) spectrophotometry. Cu is found in two distinct type 2 (T2) copper centers, one of which appears to be novel; heme is bound as a low-spin species, implying coordination by two axial ligands. MnxE and MnxF do not oxidize Mn in the absence of MnxG and are the first accessory proteins to be required by an MCO. This may indicate that Cu and heme play roles in electron transfer and/or Cu trafficking.

In this study, strain-resolved metagenomics was used to solve a mystery. A 6.4-Mbp complete closed genome was recovered from a soil metagenome and found to be astonishingly similar to that of Delftia acidovorans SPH-1, which was isolated in Germany a decade ago. It was suspected that this organism was not native to the soil sample because it lacked the diversity that is characteristic of other soil organisms; this suspicion was confirmed when PCR testing failed to detect the bacterium in the original soil samples. D. acidovorans was also identified in 16 previously published metagenomes from multiple environments, but detailed-scale single nucleotide polymorphism analysis grouped these into five distinct clades. All of the strains indicated as contaminants fell into one clade. Fragment length anomalies were identified in paired reads mapping to the contaminant clade genotypes only. This finding was used to establish that the DNA was present in specific size selection reagents used during sequencing. Ultimately, the source of the contaminant was identified as bacterial biofilms growing in tubing. On the basis of direct measurement of the rate of fixation of mutations across the period of time in which contamination was occurring, we estimated the time of separation of the contaminant strain from the genomically sequenced ancestral population within a factor of 2. This research serves as a case study of high-resolution microbial forensics and strain tracking accomplished through metagenomics-based comparative genomics. The specific case reported here is unusual in that the study was conducted in the background of a soil metagenome and the conclusions were confirmed by independent methods.IMPORTANCE It is often important to determine the source of a microbial strain. Examples include tracking a bacterium linked to a disease epidemic, contaminating the food supply, or used in bioterrorism. Strain identification and tracking are generally approached by using cultivation

In a natural geochemical cycle, manganese-oxide minerals (MnOx ) are principally formed through a microbial process, where a putative multicopper oxidase MnxG plays an essential role. Recent success in isolating the approximately 230 kDa, enzymatically active MnxEFG protein complex, has advanced our understanding of biogenic MnOx mineralization. Here, the kinetics of MnOx formation catalyzed by MnxEFG are examined using a quartz crystal microbalance (QCM), and the first electrochemical characterization of the MnxEFG complex is reported using Fourier transformed alternating current voltammetry. The voltammetric studies undertaken using near-neutral solutions (pH 7.8) establish the apparent reversible potentials for the Type 2 Cu sites in MnxEFG immobilized on a carboxy-terminated monolayer to be in the range 0.36-0.40 V versus a normal hydrogen electrode. Oxidative priming of the MnxEFG protein complex substantially enhances the enzymatic activity, as found by in situ electrochemical QCM analysis. The biogeochemical significance of this enzyme is clear, although the role of an oxidative priming of catalytic activity might be either an evolutionary advantage or an ancient relic of primordial existence.

In a natural geochemical cycle, manganese-oxide minerals (MnOx ) are principally formed through a microbial process, where a putative multicopper oxidase MnxG plays an essential role. Recent success in isolating the approximately 230 kDa, enzymatically active MnxEFG protein complex, has advanced our understanding of biogenic MnOx mineralization. Here, the kinetics of MnOx formation catalyzed by MnxEFG are examined using a quartz crystal microbalance (QCM), and the first electrochemical characterization of the MnxEFG complex is reported using Fourier transformed alternating current voltammetry. The voltammetric studies undertaken using near-neutral solutions (pH 7.8) establish the apparent reversible potentials for the Type 2 Cu sites in MnxEFG immobilized on a carboxy-terminated monolayer to be in the range 0.36-0.40 V versus a normal hydrogen electrode. Oxidative priming of the MnxEFG protein complex substantially enhances the enzymatic activity, as found by in situ electrochemical QCM analysis. The biogeochemical significance of this enzyme is clear, although the role of an oxidative priming of catalytic activity might be either an evolutionary advantage or an ancient relic of primordial existence.

In this study, strain-resolved metagenomics was used to solve a mystery. A 6.4-Mbp complete closed genome was recovered from a soil metagenome and found to be astonishingly similar to that of Delftia acidovorans SPH-1, which was isolated in Germany a decade ago. It was suspected that this organism was not native to the soil sample because it lacked the diversity that is characteristic of other soil organisms; this suspicion was confirmed when PCR testing failed to detect the bacterium in the original soil samples. D. acidovorans was also identified in 16 previously published metagenomes from multiple environments, but detailed-scale single nucleotide polymorphism analysis grouped these into five distinct clades. All of the strains indicated as contaminants fell into one clade. Fragment length anomalies were identified in paired reads mapping to the contaminant clade genotypes only. This finding was used to establish that the DNA was present in specific size selection reagents used during sequencing. Ultimately, the source of the contaminant was identified as bacterial biofilms growing in tubing. On the basis of direct measurement of the rate of fixation of mutations across the period of time in which contamination was occurring, we estimated the time of separation of the contaminant strain from the genomically sequenced ancestral population within a factor of 2. This research serves as a case study of high-resolution microbial forensics and strain tracking accomplished through metagenomics-based comparative genomics. The specific case reported here is unusual in that the study was conducted in the background of a soil metagenome and the conclusions were confirmed by independent methods.IMPORTANCE It is often important to determine the source of a microbial strain. Examples include tracking a bacterium linked to a disease epidemic, contaminating the food supply, or used in bioterrorism. Strain identification and tracking are generally approached by using cultiva

The tree of life is one of the most important organizing principles in biology(1). Gene surveys suggest the existence of an enormous number of branches(2), but even an approximation of the full scale of the tree has remained elusive. Recent depictions of the tree of life have focused either on the nature of deep evolutionary relationships(3-5) or on the known, well-classified diversity of life with an emphasis on eukaryotes(6). These approaches overlook the dramatic change in our understanding of life's diversity resulting from genomic sampling of previously unexamined environments. New methods to generate genome sequences illuminate the identity of organisms and their metabolic capacities, placing them in community and ecosystem contexts(7,8). Here, we use new genomic data from over 1,000 uncultivated and little known organisms, together with published sequences, to infer a dramatically expanded version of the tree of life, with Bacteria, Archaea and Eukarya included. The depiction is both a global overview and a snapshot of the diversity within each major lineage. The results reveal the dominance of bacterial diversification and underline the importance of organisms lacking isolated representatives, with substantial evolution concentrated in a major radiation of such organisms. This tree highlights major lineages currently underrepresented in biogeochemical models and identifies radiations that are probably important for future evolutionary analyses.

The tree of life is one of the most important organizing principles in biology(1). Gene surveys suggest the existence of an enormous number of branches(2), but even an approximation of the full scale of the tree has remained elusive. Recent depictions of the tree of life have focused either on the nature of deep evolutionary relationships(3-5) or on the known, well-classified diversity of life with an emphasis on eukaryotes(6). These approaches overlook the dramatic change in our understanding of life's diversity resulting from genomic sampling of previously unexamined environments. New methods to generate genome sequences illuminate the identity of organisms and their metabolic capacities, placing them in community and ecosystem contexts(7,8). Here, we use new genomic data from over 1,000 uncultivated and little known organisms, together with published sequences, to infer a dramatically expanded version of the tree of life, with Bacteria, Archaea and Eukarya included. The depiction is both a global overview and a snapshot of the diversity within each major lineage. The results reveal the dominance of bacterial diversification and underline the importance of organisms lacking isolated representatives, with substantial evolution concentrated in a major radiation of such organisms. This tree highlights major lineages currently underrepresented in biogeochemical models and identifies radiations that are probably important for future evolutionary analyses.

The dynamics of manganese solid formation (as MnOx) by the multicopper oxidase (MCO)-containing Mnx protein complex were examined by electron paramagnetic resonance (EPR) spectroscopy. Continuous-wave (CW) EPR spectra of samples of Mnx, prepared in atmosphere and then reacted with Mn(II) for times ranging from 7 to 600 s, indicate rapid oxidation of the substrate manganese (with two-phase pseudo-first-order kinetics modeled using rate coefficients of: k(1obs) = 0.205 ± 0.001 s(-1) and k(2obs) = 0.019 ± 0.001 s(-1)). This process occurs on approximately the same time scale as in vitro solid MnOx formation when there is a large excess of Mn(II). We also found CW and pulse EPR spectroscopic evidence for at least three classes of Mn(II)-containing species in the reaction mixtures: (i) aqueous Mn(II), (ii) a specifically bound mononuclear Mn(II) ion coordinated to the Mnx complex by one nitrogenous ligand, and (iii) a weakly exchange-coupled dimeric Mn(II) species. These findings provide new insights into the molecular mechanism of manganese mineralization.

The dynamics of manganese solid formation (as MnOx) by the multicopper oxidase (MCO)-containing Mnx protein complex were examined by electron paramagnetic resonance (EPR) spectroscopy. Continuous-wave (CW) EPR spectra of samples of Mnx, prepared in atmosphere and then reacted with Mn(II) for times ranging from 7 to 600 s, indicate rapid oxidation of the substrate manganese (with two-phase pseudo-first-order kinetics modeled using rate coefficients of: k(1obs) = 0.205 ± 0.001 s(-1) and k(2obs) = 0.019 ± 0.001 s(-1)). This process occurs on approximately the same time scale as in vitro solid MnOx formation when there is a large excess of Mn(II). We also found CW and pulse EPR spectroscopic evidence for at least three classes of Mn(II)-containing species in the reaction mixtures: (i) aqueous Mn(II), (ii) a specifically bound mononuclear Mn(II) ion coordinated to the Mnx complex by one nitrogenous ligand, and (iii) a weakly exchange-coupled dimeric Mn(II) species. These findings provide new insights into the molecular mechanism of manganese mineralization.