Metagenomic clustering links specific metabolic functions to globally relevant ecosystems.

Metagenomic sequencing has advanced our understanding of biogeochemical processes by providing an unprecedented view into the microbial composition of different ecosystems. While the amount of metagenomic data has grown rapidly, simple-to-use methods to analyze and compare across studies have lagged behind. Thus, tools expressing the metabolic traits of a community are needed to broaden the utility of existing data. Gene abundance profiles are a relatively low-dimensional embedding of a metagenome's functional potential and are, thus, tractable for comparison across many samples. Here, we compare the abundance of KEGG Ortholog Groups (KOs) from 6,539 metagenomes from the Joint Genome Institute's Integrated Microbial Genomes and Metagenomes (JGI IMG/M) database. We find that samples cluster into terrestrial, aquatic, and anaerobic ecosystems with marker KOs reflecting adaptations to these environments. For instance, functional clusters were differentiated by the metabolism of antibiotics, photosynthesis, methanogenesis, and surprisingly GC content. Using this functional gene approach, we reveal the broad-scale patterns shaping microbial communities and demonstrate the utility of ortholog abundance profiles for representing a rapidly expanding body of metagenomic data. IMPORTANCE: Metagenomics, or the sequencing of DNA from complex microbiomes, provides a view into the microbial composition of different environments. Metagenome databases were created to compile sequencing data across studies, but it remains challenging to compare and gain insight from these large data sets. Consequently, there is a need to develop accessible approaches to extract knowledge across metagenomes. The abundance of different orthologs (i.e., genes that perform a similar function across species) provides a simplified representation of a metagenome's metabolic potential that can easily be compared with others. In this study, we cluster the ortholog abundance profiles of thousands of metagenomes from diverse environments and uncover the traits that distinguish them. This work provides a simple to use framework for functional comparison and advances our understanding of how the environment shapes microbial communities.

Impact of aerobic granular sludge sizes and dissolved oxygen concentration on greenhouse gas N2O emission.

Aerobic granular sludge (AGS) at wastewater treatment plants (WWTPs) are known to produce nitrous oxide (N2O), a greenhouse gas which has a ∼300 times higher global warming potential than carbon dioxide. In this research, we studied N2O emissions from different sizes of AGS developed at a dissolved oxygen (DO) level of 2 mgO2/L while exposing them to disturbances at various DO concentrations ranging from 1 to 4 mgO2/L. Five different AGS size classes were studied: 212-600 µm, 600-1000 µm, 1000-1400 µm, 1400-2000 µm, and > 2000 µm. Metagenomic data showed N2O reductase genes (nosZ) were more abundant in the smaller AGS sizes which aligned with the observation of higher N2O reduction rates in small AGS under anaerobic conditions. However, when oxygen was present, the activity measurements of N2O emission showed an opposite trend compared to metagenomic data, smaller AGS (212 to 1000 µm) emitted significantly higher N2O (p < 0.05) than larger AGS (1000 µm to >2000 µm) at DO of 2, 3, and 4 mgO2/L. The N2O emission rate showed positive correlation with both oxygen levels and nitrification rate. This pattern indicates a connection between N2O emission and nitrification. In addition, the data suggested the penetration of oxygen into the anoxic zone of granules might have hindered nitrous oxide reduction, resulting in incomplete denitrification stopping at N2O and consequently contributing to an increase in N2O emissions. This work sets the stage to better understand the impacts of AGS size on N2O emissions in WWTPs under different disturbance of DO conditions, and thus ensure that wastewater treatment will comply with possible future regulations demanding lowering greenhouse gas emissions in an effort to combat climate change.

Coupling methanotrophic denitrification to anammox in a moving bed biofilm reactor for nitrogen removal under hypoxic conditions.

Simultaneous removal of ammonium and nitrate was achieved in a methane-fed moving bed biofilm reactor (MBBR). In the reactor, methanotrophic microorganisms oxidized methane under hypoxic conditions likely to methanol, hence providing an electron donor to denitrifiers to reduce nitrate to nitrite that then allowed anaerobic ammonium oxidizing bacteria (Anammox) to remove excess ammonium as N2. The ammonium and nitrate removal rates reached 72.09 ± 5.81 mgNH4+-N/L/d and 62.61 ± 4.17 mgNO3--N/L/d when the MBBR was operated in continuous mode. Nitrate removal by the methane-fed mixed consortia was confirmed in a batch test revealing a CH4/NO3- molar removal ratio of 1.15. The functional populations were unveiled by FISH analysis and 16S rRNA gene sequencing, which showed that the biofilm was dominated by Anammox bacteria (Candidatus Kuenenia) and diverse taxa associated with the capacity for denitrification: aerobic methanotrophs (Methylobacter, Methylomonas, and unclassified Methylococcaceae), methylotrophic denitrifiers (Opitutaceae and Methylophilaceae), and other heterotrophic denitrifiers (Ignavibacteriaceae, Anaerolineaceae, Comamonadaceae, Rhodocyclaceae and Thauera). Neither DAMO archaea nor DAMO bacteria were found in the sequencing analysis, indicating that more unknown community members possess the metabolic capacity of methanotrophic denitrification.

Metagenomic Insights Into Competition Between Denitrification and Dissimilatory Nitrate Reduction to Ammonia Within One-Stage and Two-Stage Partial-Nitritation Anammox Bioreactor Configurations.

Anaerobic ammonia oxidizing bacteria (Anammox) are implemented in high-efficiency wastewater treatment systems operated in two general configurations; one-stage systems combine aerobic ammonia oxidizing bacteria (AOB) and Anammox within a single aerated reactor, whereas two-stage configurations separate these processes into discrete tanks. Within both configurations heterotrophic populations that perform denitrification or dissimilatory nitrate reduction to ammonia (DNRA) compete for carbon and nitrate or nitrite and can impact reactor performance because DNRA retains nitrogen in the system. Therefore, it is important to understand how selective pressures imposed by one-stage and two-stage reactor configurations impact the microbial community structure and associated nitrogen transforming functions. We performed 16S rRNA gene and metagenomic sequencing on different biomass fractions (granules, flocs, and suspended biomass) sampled from two facilities treating sludge dewatering centrate: a one-stage treatment facility (Chambers Creek, Tacoma, WA) and a two-stage system (Rotterdam, Netherlands). Similar microbial populations were identified across the different samples, but relative abundances differed between reactor configurations and biomass sources. Analysis of metagenome assembled genomes (MAGs) indicated different lifestyles for abundant heterotrophic populations. Acidobacteria, Bacteroidetes, and Chloroflexi MAGs had varying capacity for DNRA and denitrification. Acidobacteria MAGs possessed high numbers of glycosyl hydrolases and glycosyl transferases indicating a role in biomass degradation. Ignavibacteria and Phycosphaerae MAGs contributed to the greater relative abundance of DNRA associated nrf genes in the two-stage granules and contained genomic features suggesting a preference for an anoxic or microoxic niche. In the one-stage granules a MAG assigned to Burkholderiales accounted for much of the abundant denitrification genes and had genomic features, including the potential for autotrophic denitrification using reduced sulfur, that indicate an ability to adapt its physiology to varying redox conditions. Overall, the competition for carbon substrates between denitrifying and DNRA performing heterotrophs may be impacted by configuration specific selective pressures. In one-stage systems oxygen availability in the bulk liquid and the oxygen gradient within granules would provide a greater niche space for heterotrophic populations capable of utilizing both oxygen and nitrate or nitrite as terminal electron acceptors, compared to two-stage systems where a homogeneous anoxic environment would favor heterotrophic populations primarily adapted to anaerobic metabolism.

Application of pyritic sludge with an anaerobic granule consortium for nitrate removal in low carbon systems.

Granules recovered from a highly reduced anaerobic digester were capable of active nitrogen removal in the absence of exogenous electron donors, averaging 0.25 mg mgNO3--N /gVSS/d over 546 days of operation. Electron mass balance indicated that about half the influent nitrate was converted to ammonia via DNRA and another half denitrified. This capacity was associated with an onion-like structure of multiple layers enriched in reduced iron and sulfur, and a complex microbial community shown by metagenomic sequencing to consist of multiple physiological groups and associated activities, including methanogenesis, denitrification, dissimilatory nitrate reduction to ammonia (DNRA), iron oxidation and reduction, and sulfur reduction and oxidation. Nitrate reduction was supported by both entrained organic material and reduced iron and sulfur species, corresponding to 2.13 mg COD/gVSS/d. Batch incubations showed that approximately 15% of denitrified nitrate was coupled to the oxidation of sulfur derived from both sulfate respiration and granular material enriched in iron-sulfide. Inhibition of sulfate reduction resulted in redirection of electron flow to methanogenesis and, in combination with other batch tests, showed that these granules supported a complex microbial community in which cryptic redox cycles linked carbon, sulfur, and iron oxidation with nitrate, sulfate, iron, and carbon dioxide reduction. This system shows promise for treatment of nitrate contaminated ground water without addition of an external organic carbon source as well as wastewater treatment in combination with (granular) sludge elimination leading in a net reduction of solid treatment costs.

Phytoplankton exudates and lysates support distinct microbial consortia with specialized metabolic and ecophysiological traits.

Blooms of marine phytoplankton fix complex pools of dissolved organic matter (DOM) that are thought to be partitioned among hundreds of heterotrophic microbes at the base of the food web. While the relationship between microbial consumers and phytoplankton DOM is a key component of marine carbon cycling, microbial loop metabolism is largely understood from model organisms and substrates. Here, we took an untargeted approach to measure and analyze partitioning of four distinct phytoplankton-derived DOM pools among heterotrophic populations in a natural microbial community using a combination of ecogenomics, stable isotope probing (SIP), and proteomics. Each 13C-labeled exudate or lysate from a diatom or a picocyanobacterium was preferentially assimilated by different heterotrophic taxa with specialized metabolic and physiological adaptations. Bacteroidetes populations, with their unique high-molecular-weight transporters, were superior competitors for DOM derived from diatom cell lysis, rapidly increasing growth rates and ribosomal protein expression to produce new relatively high C:N biomass. Proteobacteria responses varied, with relatively low levels of assimilation by Gammaproteobacteria populations, while copiotrophic Alphaproteobacteria such as the Roseobacter clade, with their diverse array of ABC- and TRAP-type transporters to scavenge monomers and nitrogen-rich metabolites, accounted for nearly all cyanobacteria exudate assimilation and produced new relatively low C:N biomass. Carbon assimilation rates calculated from SIP data show that exudate and lysate from two common marine phytoplankton are being used by taxonomically distinct sets of heterotrophic populations with unique metabolic adaptations, providing a deeper mechanistic understanding of consumer succession and carbon use during marine bloom events.

An investigation into the optimal granular sludge size for simultaneous nitrogen and phosphate removal.

An aerobic granular sludge (AGS) pilot plant fed with a mixture of acetate amended centrate and secondary effluent was used to investigate the optimal granule size range for simultaneous nitrification and denitrification (SND) and ortho-phosphate removal. The anaerobic phase was mixed to understand how AGS will perform if integrated with a continuous flow activated sludge system that cannot feed the influent through the settled sludge bed. Five different granule size fractions were taken from the pilot (operated at DO setpoint of 2mgO2/L) and each size was subjected to activity tests in a well-controlled lab-scale AGS reactor at four dissolved oxygen (DO) concentrations of 1, 2, 3, and 4 mgO2/L. The size fractions were: 212 - 600 µm, 600 - 1000 µm, 1000 - 1400 µm, 1400 - 2000 µm, and >2000 µm. The smallest size range (212 - 600 µm) had the highest nitrification and phosphate removal rates at DO setpoints from 1 - 3 mgO2/L, which was attributed to the higher aerobic volume fraction in small granules and hence a higher abundance of phosphorus accumulating organisms (PAO) and ammonia oxidizing bacteria (AOB). In comparison, large granules (>1000 µm) had 1.4 - 4.7 times lower ammonia oxidation rates than the smallest size range, which aligned with their lower AOB abundance relative to granule biomass. The granules with the highest anoxic volume fraction had the highest abundance of nitrite reductase genes (nir gene) but did not show the highest specific nitrogen removal rate. Instead, smaller granules (212 - 600 and 600 - 1000 µm), which had a lower nir gene abundance, had the highest specific nitrogen removal rates (1.2 - 3.1 times higher than larger granules) across all DO values except at 4 mgO2/L. At a DO setpoint of 4 mgO2/L, nitrite production by ammonia oxidation (ammonia monooxygenase) exceeded nitrite reduction by nitrite reductase in granules smaller than 1000 µm, in addition, some denitrifying heterotrophs switched to oxygen utilization in deeper layers hence suppressing denitrification activity. At the DO range of 2 - 4 mg/L, granular size had a greater effect on nutrient removal than DO. Therefore, for AGS developed at an average DO setpoint of 2 mgO2/L, selecting for size fractions in the range of 212 - 1000 µm and avoiding DO values higher than 3 mgO2/L can achieve both a higher nitrogen removal capacity and energy savings. This study is the first to investigate the influence of different DO values on SND and biological phosphorus removal performance of different aerobic granular sludge sizes.

Isothermal Titration Calorimetry for the Screening of Aflatoxin B1 Surface-Enhanced Raman Scattering Sensor Affinity Agents.

In this work, isothermal titration calorimetry (ITC) is employed as an affinity agent screening method for the surface-enhanced Raman scattering (SERS) detection of aflatoxin B1 (AFB1). AFB1, a potent carcinogen produced by a fungus that infects crops, is an important target due to the monitoring required based on its FDA regulation. Polymer affinity agents, like those studied here, have the potential to enable separation and detection of relevant small molecules such as pesticides, drugs, and biological toxins, like AFB1, especially when paired with a vibrational spectroscopy technique such as SERS. Herein, seven homopolymers were synthesized to be evaluated as AFB1 affinity agents based on hypothetical hydrogen bonding interactions. Nitrogen-inclusive poly( N-(2-aminoethyl) methacrylamide) (pAEMA) polymers and their oxygen analogs, poly(2-hydroxyethyl methacrylate) (pHEMA) were evaluated. ITC was demonstrated as an effective method for rapid screening among the polymer affinity agents. Chain lengths between seven and 39 repeat units were synthesized to study length-based variance in affinity agent performance. An ITC method was optimized and used for the rapid screening of polymer affinity agents. The results were compared to those generated by SERS. Good agreement between the ITC results and follow-up SERS sensing experiments showcased ITC's screening potential for analytical applications such as separation and detection.

Microbial Community Structure-Function Relationships in Yaquina Bay Estuary Reveal Spatially Distinct Carbon and Nitrogen Cycling Capacities.

Linking microbial community structure to ecological processes requires understanding of the functional roles among individual populations and the factors that influence their distributions. These structure-function relationships are particularly difficult to disentangle in estuaries, due to highly variable physico-chemical conditions. Yet, examining microbe-mediated turnover of resources in these "bioreactor" ecosystems is critical for understanding estuarine ecology. In this study, a combined metagenomics and metaproteomics approach was used to show that the unequal distribution of microbial populations across the Yaquina Bay estuary led to a habitat-specific taxonomic and functional structure and a clear spatial distribution in microbe-mediated capacities for cycling of carbon and nitrogen. For example, size-fractionation revealed that communities inhabiting suspended particulate material encoded more diverse types of metabolisms (e.g., fermentation and denitrification) than those with a planktonic lifestyle, suggesting that the metabolic reactions can differ between size fractions of the same parcel of an estuarine water column. Similarly, communities inhabiting oligotrophic conditions in the lower estuary were enriched in genes involved in central carbon metabolism (e.g., TCA cycle), while communities in the upper estuary were enriched in genes typical of copiotrophic populations (e.g., cell growth, cell division). Integrating gene and protein data revealed that abundant populations of Flavobacteriales and Rhodobacterales encoded similar genomic functions, yet differed significantly in protein expression, dedicating a large proportion of their respective proteomes to rapid growth and division versus metabolic versatility and resource acquisition. This suggested potentially distinct life-strategies between these two co-occurring lineages and was concomitant with differing patterns of positive evolutionary selection on their encoded genes. Microbial communities and their functions across Yaquina Bay appear to be structured by population-level habitat preferences, resulting in spatially distinct elemental cycling, while within each community, forces such as competitive exclusion and evolutionary selection influence species life-strategies and may help maintain microbial diversity.

Phylogenetically conserved resource partitioning in the coastal microbial loop.

Resource availability influences marine microbial community structure, suggesting that population-specific resource partitioning defines discrete niches. Identifying how resources are partitioned among populations, thereby characterizing functional guilds within the communities, remains a challenge for microbial ecologists. We used proteomic stable isotope probing (SIP) and NanoSIMS analysis of phylogenetic microarrays (Chip-SIP) along with 16S rRNA gene amplicon and metagenomic sequencing to characterize the assimilation of six 13C-labeled common metabolic substrates and changes in the microbial community structure within surface water collected from Monterey Bay, CA. Both sequencing approaches indicated distinct substrate-specific community shifts. However, observed changes in relative abundance for individual populations did not correlate well with directly measured substrate assimilation. The complementary SIP techniques identified assimilation of all six substrates by diverse taxa, but also revealed differential assimilation of substrates into protein and ribonucleotide biomass between taxa. Substrate assimilation trends indicated significantly conserved resource partitioning among populations within the Flavobacteriia, Alphaproteobacteria and Gammaproteobacteria classes, suggesting that functional guilds within marine microbial communities are phylogenetically cohesive. However, populations within these classes exhibited heterogeneity in biosynthetic activity, which distinguished high-activity copiotrophs from low-activity oligotrophs. These results indicate distinct growth responses between populations that is not apparent by genome sequencing alone.

Proteomic Stable Isotope Probing Reveals Taxonomically Distinct Patterns in Amino Acid Assimilation by Coastal Marine Bacterioplankton.

Heterotrophic marine bacterioplankton are a critical component of the carbon cycle, processing nearly a quarter of annual primary production, yet defining how substrate utilization preferences and resource partitioning structure microbial communities remains a challenge. In this study, proteomic stable isotope probing (proteomic SIP) was used to characterize population-specific assimilation of dissolved free amino acids (DFAAs), a major source of dissolved organic carbon for bacterial secondary production in aquatic environments. Microcosms of seawater collected from Newport, Oregon, and Monterey Bay, California, were incubated with 1 µM 13C-labeled amino acids for 15 and 32 h. The taxonomic compositions of microcosm metaproteomes were highly similar to those of the sampled natural communities, with Rhodobacteriales, SAR11, and Flavobacteriales representing the dominant taxa. Analysis of 13C incorporation into protein biomass allowed for quantification of the isotopic enrichment of identified proteins and subsequent determination of differential amino acid assimilation patterns between specific bacterioplankton populations. Proteins associated with Rhodobacterales tended to have a significantly high frequency of 13C-enriched peptides, opposite the trend for Flavobacteriales and SAR11 proteins. Rhodobacterales proteins associated with amino acid transport and metabolism had an increased frequency of 13C-enriched spectra at time point 2. Alteromonadales proteins also had a significantly high frequency of 13C-enriched peptides, particularly within ribosomal proteins, demonstrating their rapid growth during incubations. Overall, proteomic SIP facilitated quantitative comparisons of DFAA assimilation by specific taxa, both between sympatric populations and between protein functional groups within discrete populations, allowing an unprecedented examination of population level metabolic responses to resource acquisition in complex microbial communities. IMPORTANCE An estimated 50 gigatons of carbon is annually fixed within marine systems, of which heterotrophic microbial populations process nearly half. These communities vary in composition and activity across spatial and temporal scales, so understanding how these changes affect global processes requires the delineation of functional roles for individual members. In a step toward ascertaining these roles, we applied proteomic stable isotope probing to quantify the assimilation of organic carbon from DFAAs into microbial protein biomass, since the turnover of DFAAs accounts for a substantial fraction of marine microbial carbon metabolism that is directed into biomass production. We conducted experiments at two coastal North Pacific locations and found taxonomically distinct responses. This approach allowed us to compare amino acid assimilation by specific bacterioplankton populations and characterize their allocation of this substrate among cellular functions.

A novel sister clade to the enterobacteria microviruses (family Microviridae) identified in methane seep sediments.

Methane seep microbial communities perform a key ecosystem service by consuming the greenhouse gas methane prior to its release into the hydrosphere, minimizing the impact of marine methane sources on our climate. Although previous studies have examined the ecology and biochemistry of these communities, none has examined viral assemblages associated with these habitats. We employed virus particle purification, genome amplification, pyrosequencing and gene/genome reconstruction and annotation on two metagenomic libraries, one prepared for ssDNA and the other for all DNA, to identify the viral community in a methane seep. Similarity analysis of these libraries (raw and assembled) revealed a community dominated by phages, with a significant proportion of similarities to the Microviridae family of ssDNA phages. We define these viruses as the Eel River Basin Microviridae (ERBM). Assembly and comparison of 21 ERBM closed circular genomes identified five as members of a novel sister clade to the Microvirus genus of Enterobacteria phages. Comparisons among other metagenomes and these Microviridae major-capsid sequences indicated that this clade of phages is currently unique to the Eel River Basin sediments. Given this ERBM clade's relationship to the Microviridae genus Microvirus, we define this sister clade as the candidate genus Pequeñovirus.