Alyse K. Hawley

Ecology and evolution of viruses infecting uncultivated SUP05 bacteria as revealed by single-cell- and meta- genomics

Viruses modulate microbial communities and alter ecosystem functions. However, due to cultivation bottlenecks, specific virus–host interaction dynamics remain cryptic. In this study, we examined 127 single-cell amplified genomes (SAGs) from uncultivated SUP05 bacteria isolated from a model marine oxygen minimum zone (OMZ) to identify 69 viral contigs representing five new genera within dsDNA Caudovirales and ssDNA Microviridae. Infection frequencies suggest that ∼1/3 of SUP05 bacteria is viral-infected, with higher infection frequency where oxygen-deficiency was most severe. Observed Microviridae clonality suggests recovery of bloom-terminating viruses, while systematic co-infection between dsDNA and ssDNA viruses posits previously unrecognized cooperation modes. Analyses of 186 microbial and viral metagenomes revealed that SUP05 viruses persisted for years, but remained endemic to the OMZ. Finally, identification of virus-encoded dissimilatory sulfite reductase suggests SUP05 viruses reprogram their hosts energy metabolism. Together, these results demonstrate closely coupled SUP05 virus–host co-evolutionary dynamics with the potential to modulate biogeochemical cycling in climate-critical and expanding OMZs. DOI: http://dx.doi.org/10.7554/eLife.03125.001

Metaproteomics reveals differential modes of metabolic coupling among ubiquitous oxygen minimum zone microbes

Significance Oxygen is an important organizing principle in marine ecosystems. As oxygen levels decline, energy is increasingly diverted away from higher trophic levels into microbial community metabolism causing changes in carbon and nutrient cycling. Here we use metagenomic and metaproteomic methods to chart in situ metabolic networks linking key microbial players driving carbon and nutrient cycling in a seasonally stratified fjord, Saanich Inlet, a model ecosystem for studying microbial responses to changing levels of water column oxygen deficiency. Based on this evidence, we develop a conceptual model that describes coupling of chemotrophic energy production with dark carbon fixation along defined redox gradients with implications for primary production and possibly carbon sedimentation in expanding marine oxygen minimum zones. Marine oxygen minimum zones (OMZs) are intrinsic water column features arising from respiratory oxygen demand during organic matter degradation in stratified waters. Currently OMZs are expanding due to global climate change with resulting feedback on marine ecosystem function. Here we use metaproteomics to chart spatial and temporal patterns of gene expression along defined redox gradients in a seasonally stratified fjord to better understand microbial community responses to OMZ expansion. The expression of metabolic pathway components for nitrification, anaerobic ammonium oxidation (anammox), denitrification, and inorganic carbon fixation were differentially expressed across the redoxcline and covaried with distribution patterns of ubiquitous OMZ microbes including Thaumarchaeota, Nitrospina, Nitrospira, Planctomycetes, and SUP05/ARCTIC96BD-19 Gammaproteobacteria. Nitrification and inorganic carbon fixation pathways affiliated with Thaumarchaeota dominated dysoxic waters, and denitrification, sulfur oxidation, and inorganic carbon fixation pathways affiliated with the SUP05 group of nitrate-reducing sulfur oxidizers dominated suboxic and anoxic waters. Nitrifier nitrite oxidation and anammox pathways affiliated with Nirospina, Nitrospira, and Planctomycetes, respectively, also exhibited redox partitioning between dysoxic and suboxic waters. The numerical abundance of SUP05 proteins mediating inorganic carbon fixation under anoxic conditions suggests that SUP05 will become increasingly important in global ocean carbon and nutrient cycling as OMZs expand.

Metabolic pathways for the whole community.

BackgroundA convergence of high-throughput sequencing and computational power is transforming biology into information science. Despite these technological advances, converting bits and bytes of sequence information into meaningful insights remains a challenging enterprise. Biological systems operate on multiple hierarchical levels from genomes to biomes. Holistic understanding of biological systems requires agile software tools that permit comparative analyses across multiple information levels (DNA, RNA, protein, and metabolites) to identify emergent properties, diagnose system states, or predict responses to environmental change.ResultsHere we adopt the MetaPathways annotation and analysis pipeline and Pathway Tools to construct environmental pathway/genome databases (ePGDBs) that describe microbial community metabolism using MetaCyc, a highly curated database of metabolic pathways and components covering all domains of life. We evaluate Pathway Tools’ performance on three datasets with different complexity and coding potential, including simulated metagenomes, a symbiotic system, and the Hawaii Ocean Time-series. We define accuracy and sensitivity relationships between read length, coverage and pathway recovery and evaluate the impact of taxonomic pruning on ePGDB construction and interpretation. Resulting ePGDBs provide interactive metabolic maps, predict emergent metabolic pathways associated with biosynthesis and energy production and differentiate between genomic potential and phenotypic expression across defined environmental gradients.ConclusionsThis multi-tiered analysis provides the user community with specific operating guidelines, performance metrics and prediction hazards for more reliable ePGDB construction and interpretation. Moreover, it demonstrates the power of Pathway Tools in predicting metabolic interactions in natural and engineered ecosystems.

Integrating biogeochemistry with multiomic sequence information in a model oxygen minimum zone

Significance Modern molecular sequencing is beginning to provide great insight into microbial community structure and function at ecosystem scales. However, the quantitative integration of multiomic sequence information (DNA, mRNA, and protein) and geochemical processes has so far been elusive. Here, we develop a biogeochemical model that integrates geochemistry and multiomic sequence information to explain key metabolic processes in the oxygen-starved waters of Saanich Inlet, a model ecosystem for studying microbial community responses to oxygen minimum zone expansion. Our model largely explains DNA, mRNA, and protein distributions and sheds light on the metabolic networks coupling carbon, sulfur, and nitrogen transformations across a redox gradient. Our approach is extensible to other biogeochemical models incorporating feedbacks of global change on ecosystem functions. Microorganisms are the most abundant lifeform on Earth, mediating global fluxes of matter and energy. Over the past decade, high-throughput molecular techniques generating multiomic sequence information (DNA, mRNA, and protein) have transformed our perception of this microcosmos, conceptually linking microorganisms at the individual, population, and community levels to a wide range of ecosystem functions and services. Here, we develop a biogeochemical model that describes metabolic coupling along the redox gradient in Saanich Inlet—a seasonally anoxic fjord with biogeochemistry analogous to oxygen minimum zones (OMZs). The model reproduces measured biogeochemical process rates as well as DNA, mRNA, and protein concentration profiles across the redox gradient. Simulations make predictions about the role of ubiquitous OMZ microorganisms in mediating carbon, nitrogen, and sulfur cycling. For example, nitrite “leakage” during incomplete sulfide-driven denitrification by SUP05 Gammaproteobacteria is predicted to support inorganic carbon fixation and intense nitrogen loss via anaerobic ammonium oxidation. This coupling creates a metabolic niche for nitrous oxide reduction that completes denitrification by currently unidentified community members. These results quantitatively improve previous conceptual models describing microbial metabolic networks in OMZs. Beyond OMZ-specific predictions, model results indicate that geochemical fluxes are robust indicators of microbial community structure and reciprocally, that gene abundances and geochemical conditions largely determine gene expression patterns. The integration of real observational data, including geochemical profiles and process rate measurements as well as metagenomic, metatranscriptomic and metaproteomic sequence data, into a biogeochemical model, as shown here, enables holistic insight into the microbial metabolic network driving nutrient and energy flow at ecosystem scales.

Large Insert Environmental Genomic Library Production

The vast majority of microbes in nature currently remain inaccessible to traditional cultivation methods. Over the past decade, culture-independent environmental genomic (i.e. metagenomic) approaches have emerged, enabling researchers to bridge this cultivation gap by capturing the genetic content of indigenous microbial communities directly from the environment. To this end, genomic DNA libraries are constructed using standard albeit artful laboratory cloning techniques. Here we describe the construction of a large insert environmental genomic fosmid library with DNA derived from the vertical depth continuum of a seasonally hypoxic fjord. This protocol is directly linked to a series of connected protocols including coastal marine water sampling [1], large volume filtration of microbial biomass [2] and a DNA extraction and purification protocol [3]. At the outset, high quality genomic DNA is end-repaired with the creation of 5 -phosphorylated blunt ends. End-repaired DNA is subjected to pulsed-field gel electrophoresis (PFGE) for size selection and gel extraction is performed to recover DNA fragments between 30 and 60 thousand base pairs (Kb) in length. Size selected DNA is purified away from the PFGE gel matrix and ligated to the phosphatase-treated blunt-end fosmid CopyControl vector pCC1 (EPICENTRE http://www.epibio.com/item.asp?ID=385). Linear concatemers of pCC1 and insert DNA are subsequently headfull packaged into phage particles by lambda terminase, with subsequent infection of phage-resistant E. coli cells. Successfully transduced clones are recovered on LB agar plates under antibiotic selection and archived in 384-well plate format using an automated colony picking robot (Qpix2, GENETIX). The current protocol draws from various sources including the CopyControl Fosmid Library Production Kit from EPICENTRE and the published works of multiple research groups [4-7]. Each step is presented with best practice in mind. Whenever possible we highlight subtleties in execution to improve overall quality and efficiency of library production. The whole process of fosmid library production and automated colony picking takes at least 7-10 days as there are many incubation steps included. However, there are several stopping points possible which are mentioned within the protocol.

Seawater Sampling and Collection

This video documents methods for collecting coastal marine water samples and processing them for various downstream applications including biomass concentration. nucleic acid purification, cell abundance, nutrient and trace gas analyses. For todays demonstration samples were collected from the deck of the HMS John Strickland operating in Saanich Inlet. An A-frame derrick, with a multi-purpose winch and cable system, is used in combination with Niskin or Go-Flo water sampling bottles. A Conductivity, Temperature, and Depth (CTD) sensor is also be used to sample the underlying water mass. To minimize outgassing, trace gas samples are collected first. Then, nutrients, chemistry, and cell counts are determined. Finally, waters are collected for biomass filtration. The set-up and collection time for a single cast is approximately 1.5 hours at a maximum depth of 215 meters. Therefore, a total of 6 hours is generally needed to complete the four-part collection series described here.

Diverse marinimicrobia bacteria may mediate coupled biogeochemical cycles along eco-thermodynamic gradients

Microbial communities drive biogeochemical cycles through networks of metabolite exchange that are structured along energetic gradients. As energy yields become limiting, these networks favor co-metabolic interactions to maximize energy disequilibria. Here we apply single-cell genomics, metagenomics, and metatranscriptomics to study bacterial populations of the abundant “microbial dark matter” phylum Marinimicrobia along defined energy gradients. We show that evolutionary diversification of major Marinimicrobia clades appears to be closely related to energy yields, with increased co-metabolic interactions in more deeply branching clades. Several of these clades appear to participate in the biogeochemical cycling of sulfur and nitrogen, filling previously unassigned niches in the ocean. Notably, two Marinimicrobia clades, occupying different energetic niches, express nitrous oxide reductase, potentially acting as a global sink for the greenhouse gas nitrous oxide.Little is known about Marinimicrobia, a group of bacteria that are prevalent in the oceans. Here, the authors study global populations of Marinimicrobia using single-cell genomics, metagenomics and metatranscriptomics, showing potential co-metabolic interactions and participation in the sulfur and nitrogen cycles.

A compendium of geochemical information from the Saanich Inlet water column

Extensive and expanding oxygen minimum zones (OMZs) exist at variable depths in coastal and open ocean waters. As oxygen levels decline, nutrients and energy are increasingly diverted away from higher trophic levels into microbial community metabolism, resulting in fixed nitrogen loss and production of climate active trace gases including nitrous oxide and methane. While ocean deoxygenation has been reported on a global scale, our understanding of OMZ biology and geochemistry is limited by a lack of time-resolved data sets. Here, we present a historical dataset of oxygen concentrations spanning fifty years and nine years of monthly geochemical time series observations in Saanich Inlet, a seasonally anoxic fjord on the coast of Vancouver Island, British Columbia, Canada that undergoes recurring changes in water column oxygenation status. This compendium provides a unique geochemical framework for evaluating long-term trends in biogeochemical cycling in OMZ waters.

Methanotrophic community dynamics in a seasonally anoxic fjord: Saanich Inlet, British Columbia

Marine oxygen minimum zones (OMZs) play disproportionate roles in nutrient and climate active trace gas cycling including nitrous oxide and methane, in the ocean. OMZs are currently expanding due to climate change making it increasingly important to identify microbial controls on trace gas cycling at the individual, population and community levels. Here we present a two-year survey of the microbial community along seasonal redox gradients in Saanich Inlet focused on identifying microbial agents driving methane oxidation. Although methanotrophs were rare, we identified three uncultivated groups affiliated with particulate methane monooxygenase (pMMO) encoding phylogenetic groups (OPU), and methanotrophic symbionts as primary drivers of methane oxidation in Saanich Inlet. Distribution and activity patterns for these three groups were consistent with niche partitioning that became increasingly resolved during water column stratification. Moreover co-occurrence analysis combined with multi-level indicator species analysis revealed significant correlations between operational taxonomic units affiliated with Methylophaga, Methylophilales, SAR324, Verrucomicrobia, and Planctomycetes with OPUs and methanotrophic symbiont groups. Taken together these observations shed new light on the composition, dynamics, and potential interspecific interactions of microbes associated with CH4 cycling in the Saanich Inlet water column, provide a baseline for comparison between coastal and open ocean OMZs and support the potential role of OPUs, and methanotrophic symbiont groups as a widely distributed pelagic sink for CH4 along continental margins.

Molecular Tools for Investigating Microbial Community Structure and Function in Oxygen-Deficient Marine Waters

Water column oxygen (O2)-deficiency shapes food-web structure by progressively directing nutrients and energy away from higher trophic levels into microbial community metabolism resulting in fixed nitrogen loss and greenhouse gas production. Although respiratory O2 consumption during organic matter degradation is a natural outcome of a productive surface ocean, global-warming-induced stratification intensifies this process leading to oxygen minimum zone (OMZ) expansion. Here, we describe useful tools for detection and quantification of potential key microbial players and processes in OMZ community metabolism including quantitative polymerase chain reaction primers targeting Marine Group I Thaumarchaeota, SUP05, Arctic96BD-19, and SAR324 small-subunit ribosomal RNA genes and protein extraction methods from OMZ waters compatible with high-resolution mass spectrometry for profiling microbial community structure and functional dynamics.