Heterotrophic bacteria trigger transcriptome remodelling in the photosynthetic picoeukaryote Micromonas commoda

Marine biogeochemical cycles are built on interactions between surface ocean microbes, particularly those connecting phytoplankton primary producers to heterotrophic bacteria. Details of these associations are not well understood, especially in the case of direct influences of bacteria on phytoplankton physiology. Here we catalogue how the presence of three marine bacteria (Ruegeria pomeroyi DSS-3, Stenotrophomonas sp. SKA14 and Polaribacter dokdonensis MED152) individually and uniquely impact gene expression of the picoeukaryotic alga Micromonas commoda RCC 299. We find a dramatic transcriptomic remodelling by M. commoda after 8 h in co-culture, followed by an increase in cell numbers by 56 h compared with the axenic cultures. Some aspects of the algal transcriptomic response are conserved across all three bacterial co-cultures, including an unexpected reduction in relative expression of photosynthesis and carbon fixation pathways. Expression differences restricted to a single bacterium are also observed, with the Flavobacteriia P. dokdonensis uniquely eliciting changes in relative expression of algal genes involved in biotin biosynthesis and the acquisition and assimilation of nitrogen. This study reveals that M. commoda has rapid and extensive responses to heterotrophic bacteria in ways that are generalizable, as well as in a taxon specific manner, with implications for the diversity of phytoplankton-bacteria interactions ongoing in the surface ocean.

Functional annotation and importance of marine bacterial transporters of plankton exometabolites

Metabolite exchange within marine microbial communities transfers carbon and other major elements through global cycles and forms the basis of microbial interactions. Yet lack of gene annotations and concern about the quality of existing ones remain major impediments to revealing currencies of carbon flux. We employed an arrayed mutant library of the marine bacterium Ruegeria pomeroyi DSS-3 to experimentally annotate substrates of organic compound transporter systems, using mutant growth and compound drawdown analyses to link transporters to their cognate substrates. Mutant experiments verified substrates for thirteen R. pomeroyi transporters. Four were previously hypothesized based on gene expression data (taurine, glucose/xylose, isethionate, and cadaverine/putrescine/spermidine); five were previously hypothesized based on homology to experimentally annotated transporters in other bacteria (citrate, glycerol, N-acetylglucosamine, fumarate/malate/succinate, and dimethylsulfoniopropionate); and four had no previous annotations (thymidine, carnitine, cysteate, and 3-hydroxybutyrate). These bring the total number of experimentally-verified organic carbon influx transporters to 18 of 126 in the R. pomeroyi genome. In a longitudinal study of a coastal phytoplankton bloom, expression patterns of the experimentally annotated transporters linked them to different stages of the bloom, and also led to the hypothesis that citrate and 3-hydroxybutyrate were among the most highly available bacterial substrates. Improved functional annotation of the gatekeepers of organic carbon uptake is critical for deciphering carbon flux and fate in microbial ecosystems.

Dissolved organic metabolite extraction from high-salt media

We describe considerations and strategies for developing a nuclear magnetic resonance (NMR) sample preparation method to extract low molecular weight metabolites from high-salt spent media in a model coculture system of phytoplankton and marine bacteria. Phytoplankton perform half the carbon fixation and oxygen generation on Earth. A substantial fraction of fixed carbon becomes part of a metabolite pool of small molecules known as dissolved organic matter (DOM), which are taken up by marine bacteria proximate to phytoplankton. There is an urgent need to elucidate these metabolic exchanges due to widespread anthropogenic transformations on the chemical, phenotypic, and species composition of seawater. These changes are increasing water temperature and the amount of CO₂ absorbed by the ocean at energetic costs to marine microorganisms. Little is known about the metabolite-mediated, structured interactions occurring between phytoplankton and associated marine bacteria, in part because of challenges in studying high-salt solutions on various analytical platforms. NMR analysis is problematic due to the high-salt content of both natural seawater and culture media for marine microbes. High-salt concentration degrades the performance of the radio frequency coil, reduces the efficiency of some pulse sequences, limits signal-to-noise, and prolongs experimental time. The method described herein can reproducibly extract low molecular weight DOM from small-volume, high-salt cultures. It is a promising tool for elucidating metabolic flux between marine microorganisms and facilitates genetic screens of mutant microorganisms.

13C NMR metabolomics: J-resolved STOCSY meets INADEQUATE

Robust annotation of metabolites is a challenging task in metabolomics. Among available applications, ¹³C NMR experiment INADEQUATE determines direct ¹³C-¹³C connectivity unambiguously, offering indispensable information on molecular structure. Despite its great utility, it is not always practical to collect INADEQUATE data on every sample in a large metabolomics study because of its relatively long experiment time. Here, we propose an alternative approach that maintains the quality of information but saves experiment time. In this approach, individual samples in a study are first screened by ¹³C homonuclear J-resolved experiment (JRES). Next, JRES data are processed by statistical total correlation spectroscopy (STOCSY) to extract peaks that behave similarly among samples. Finally, INADEQUATE is collected on one internal pooled sample to select STOCSY peaks that originate from the same compound. We tested this concept using the ¹³C-labeled endometabolome of a model marine diatom strain incubated under various settings, intending to cover a range of metabolites produced under different external conditions. This scheme was able to extract known diatom metabolites proline, 2,3-dihydroxypropane-1-sulfonate (DHPS), β-1,3-glucan, choline, and glutamate. This pipeline also detected unknown compounds with structural information, which is valuable in metabolomics where a priori knowledge of metabolites is not always available. The ability of this scheme was seen even in sugar regions, which are usually challenging in ¹H NMR due to severe peak overlap. JRES and INADEQUATE were highly complementary; INADEQUATE provided directly-bonded ¹³C networks, whereas JRES linked INADEQUATE networks within the same compound but broken by nitrogen or sulfur atoms, highlighting the advantage of this integrated approach.

Bacterial transcriptional response to labile exometabolites from photosynthetic picoeukaryote Micromonas commoda

Dissolved primary production released into seawater by marine phytoplankton is a major source of carbon fueling heterotrophic bacterial production in the ocean. The composition of the organic compounds released by healthy phytoplankton is poorly known and difficult to assess with existing chemical methods. Here, expression of transporter and catabolic genes by three model marine bacteria (Ruegeria pomeroyi DSS-3, Stenotrophomonas sp. SKA14, and Polaribacter dokdonensis MED152) was used as a biological sensor of metabolites released from the picoeukaryote Micromonas commoda RCC299. Bacterial expression responses indicated that the three species together recognized 38 picoeukaryote metabolites. This was consistent with the Micromonas expression of genes for starch metabolism and synthesis of peptidoglycan-like intermediates. A comparison of the hypothesized Micromonas exometabolite pool with that of the diatom Thalassiosira pseudonana CCMP1335, analyzed previously with the same biological sensor method, indicated that both phytoplankton released organic acids, nucleosides, and amino acids, but differed in polysaccharide and organic nitrogen release. Future ocean conditions are expected to favor picoeukaryotic phytoplankton over larger-celled microphytoplankton. Results from this study suggest that such a shift could alter the substrate pool available to heterotrophic bacterioplankton.

Oxidative Stress Regulates a Pivotal Metabolic Switch in Dimethylsulfoniopropionate Degradation by the Marine Bacterium Ruegeria pomeroyi

Dimethylsulfoniopropionate (DMSP) is an abundant organic compound in marine surface water and source of dimethyl sulfide (DMS), the largest natural sulfur source to the upper atmosphere. Marine bacteria either mineralize DMSP through the demethylation pathway or transform it to DMS through the cleavage pathway. Factors that regulate which pathway is utilized are not fully understood. In chemostat experiments, the marine Roseobacter Ruegeria pomeroyi DSS-3 was exposed to oxidative stress either during growth with H2O2 or by mutation of the gene encoding catalase. Oxidative stress reduced expression of the genes in the demethylation pathway and increased expression of those encoding the cleavage pathway. These results are contrary to the sulfur demand hypothesis, which theorizes that DMSP metabolism is driven by sulfur requirements of bacterial cells. Instead, we find strong evidence consistent with oxidative stress control over the switch in DMSP metabolism from demethylation to DMS production in an ecologically relevant marine bacterium. IMPORTANCE Dimethylsulfoniopropionate (DMSP) is the most abundant low-molecular-weight organic compound in marine surface water and source of dimethyl sulfide (DMS), a climatically active gas that connects the marine and terrestrial sulfur cycles. Marine bacteria are the major DMSP consumers, either generating DMS or consuming DMSP as a source of reduced carbon and sulfur. However, the factors regulating the DMSP catabolism in bacteria are not well understood. Marine bacteria are also exposed to oxidative stress. RNA sequencing (RNA-seq) experiments showed that oxidative stress induced in the laboratory reduced expression of the genes encoding the consumption of DMSP via the demethylation pathway and increased the expression of genes encoding DMS production via the cleavage pathway in the marine bacterium Ruegeria pomeroyi. These results support a model where DMS production in the ocean is regulated in part by oxidative stress.

Divergent gene expression responses in two Baltic Sea heterotrophic model bacteria to dinoflagellate dissolved organic matter

Phytoplankton release massive amounts of dissolved organic matter (DOM) into the water column during recurring blooms in coastal waters and inland seas. The released DOM encompasses a complex mixture of both known and unknown compounds, and is a rich nutrient source for heterotrophic bacteria. The metabolic activity of bacteria during and after phytoplankton blooms can hence be expected to reflect the characteristics of the released DOM. We therefore investigated if bacterioplankton could be used as “living sensors” of phytoplankton DOM quantity and/or quality, by applying gene expression analyses to identify bacterial metabolisms induced by DOM. We used transcriptional analysis of two Baltic Sea bacterial isolates (Polaribacter sp. BAL334 [Flavobacteriia] and Brevundimonas sp. BAL450 [Alphaproteobacteria]) growing with DOM from axenic cultures of the dinoflagellate Prorocentrum minimum. We observed pronounced differences between the two bacteria both in growth and the expressed metabolic pathways in cultures exposed to dinoflagellate DOM compared with controls. Differences in metabolic responses between the two isolates were caused both by differences in gene repertoire between them (e.g. in the SEED categories for membrane transport, motility and photoheterotrophy) and the regulation of expression (e.g. fatty acid metabolism), emphasizing the importance of separating the responses of different taxa in analyses of community sequence data. Similarities between the bacteria included substantially increased expression of genes for Ton and Tol transport systems in both isolates, which are commonly associated with uptake of complex organic molecules. Polaribacter sp. BAL334 showed stronger metabolic responses to DOM harvested from exponential than stationary phase dinoflagellates (128 compared to 26 differentially expressed genes), whereas Brevundimonas sp. BAL450 responded more to the DOM from stationary than exponential phase dinoflagellates (33 compared to 6 differentially expressed genes). These findings suggest that shifts in bacterial metabolisms during different phases of phytoplankton blooms can be detected in individual bacterial species and can provide insights into their involvement in DOM transformations.

Diel investments in metabolite production and consumption in a model microbial system

Organic carbon transfer between surface ocean photosynthetic and heterotrophic microbes is a central but poorly understood process in the global carbon cycle. In a model community in which diatom extracellular release of organic molecules sustained growth of a co-cultured bacterium, we determined quantitative changes in the diatom endometabolome and the bacterial uptake transcriptome over two diel cycles. Of the nuclear magnetic resonance (NMR) peaks in the diatom endometabolites, 38% had diel patterns with noon or mid-afternoon maxima; the remaining either increased (36%) or decreased (26%) through time. Of the genes in the bacterial uptake transcriptome, 94% had a diel pattern with a noon maximum; the remaining decreased over time (6%). Eight diatom endometabolites identified with high confidence were matched to the bacterial genes mediating their utilization. Modeling of these coupled inventories with only diffusion-based phytoplankton extracellular release could not reproduce all the patterns. Addition of active release mechanisms for physiological balance and bacterial recognition significantly improved model performance. Estimates of phytoplankton extracellular release range from only a few percent to nearly half of annual net primary production. Improved understanding of the factors that influence metabolite release and consumption by surface ocean microbes will better constrain this globally significant carbon flux.

Microbial metabolites in the marine carbon cycle

One-quarter of photosynthesis-derived carbon on Earth rapidly cycles through a set of short-lived seawater metabolites that are generated from the activities of marine phytoplankton, bacteria, grazers and viruses. Here we discuss the sources of microbial metabolites in the surface ocean, their roles in ecology and biogeochemistry, and approaches that can be used to analyse them from chemistry, biology, modelling and data science. Although microbial-derived metabolites account for only a minor fraction of the total reservoir of marine dissolved organic carbon, their flux and fate underpins the central role of the ocean in sustaining life on Earth.

Growth-stage-related shifts in diatom endometabolome composition set the stage for bacterial heterotrophy

Phytoplankton-derived metabolites fuel a large fraction of heterotrophic bacterial production in the global ocean, yet methodological challenges have limited our understanding of the organic molecules transferred between these microbial groups. In an experimental bloom study consisting of three heterotrophic marine bacteria growing together with the diatom Thalassiosira pseudonana, we concurrently measured diatom endometabolites (i.e., potential exometabolite supply) by nuclear magnetic resonance (NMR) spectroscopy and bacterial gene expression (i.e., potential exometabolite uptake) by metatranscriptomic sequencing. Twenty-two diatom endometabolites were annotated, with nine increasing in internal concentration in the late stage of the bloom, eight decreasing, and five showing no variation through the bloom progression. Some metabolite changes could be linked to shifts in diatom gene expression, as well as to shifts in bacterial community composition and their expression of substrate uptake and catabolism genes. Yet an overall low match indicated that endometabolome concentration was not a good predictor of exometabolite availability, and that complex physiological and ecological interactions underlie metabolite exchange. Six diatom endometabolites accumulated to higher concentrations in the bacterial co-cultures compared to axenic cultures, suggesting a bacterial influence on rates of synthesis or release of glutamate, arginine, leucine, 2,3-dihydroxypropane-1-sulfonate, glucose, and glycerol-3-phosphate. Better understanding of phytoplankton metabolite production, release, and transfer to assembled bacterial communities is key to untangling this nearly invisible yet pivotal step in ocean carbon cycling.

Bacterial responses to background organic pollutants in the northeast subarctic Pacific Ocean

Thousands of man-made synthetic chemicals are released to oceans and compose the anthropogenic dissolved organic carbon (ADOC). Little is known about the effects of this chronic pollution on marine microbiome activities. In this study, we measured the pollution level at three sites in the Northeast Subarctic Pacific Ocean (NESAP) and investigated how mixtures of three model families of ADOC at different environmentally relevant concentrations affected naturally occurring marine bacterioplankton communities' structure and metabolic functioning. The offshore northernmost site (North) had the lowest concentrations of hydrocarbons, as well as organophosphate ester plasticizers, contrasting with the two other continental shelf sites, the southern coastal site (South) being the most contaminated. At North, ADOC stimulated bacterial growth and promoted an increase in the contribution of some Gammaproteobacteria groups (e.g. Alteromonadales) to the 16 rRNA pool. These groups are described as fast responders after oil spills. In contrast, minor changes in South microbiome activities were observed. Gene expression profiles at Central showed the coexistence of ADOC degradation and stress-response strategies to cope with ADOC toxicities. These results show that marine microbial communities at three distinct domains in NESAP are influenced by background concentrations of ADOC, expanding previous assessments for polar and temperate waters.

Quantification of Amine- and Alcohol-Containing Metabolites in Saline Samples Using Pre-extraction Benzoyl Chloride Derivatization and Ultrahigh Performance Liquid Chromatography Tandem Mass Spectrometry (UHPLC MS/MS)

Dissolved metabolites serve as nutrition, energy, and chemical signals for microbial systems. However, the full scope and magnitude of these processes in marine systems are unknown, largely due to insufficient methods, including poor extraction of small, polar compounds using common solid-phase extraction resins. Here, we utilized pre-extraction derivatization and ultrahigh performance liquid chromatography electrospray ionization tandem mass spectrometry (UHPLC-ESI-MS/MS) to detect and quantify targeted dissolved metabolites in seawater and saline culture media. Metabolites were derivatized with benzoyl chloride by their primary and secondary amine and alcohol functionalities and quantified using stable isotope-labeled internal standards (SIL-ISs) produced from ¹³C₆-labeled benzoyl chloride. We optimized derivatization, extraction, and sample preparation for field and culture samples and evaluated matrix-derived biases. We have optimized this quantitative method for 73 common metabolites, of which 50 cannot be quantified without derivatization due to low extraction efficiencies. Of the 73 metabolites, 66 were identified in either culture media or seawater and 45 of those were quantified. This derivatization method is sensitive (detection limits = pM to nM), rapid (∼5 min per sample), and high throughput.

Transcriptional activity differentiates families of Marine Group II Euryarchaeota in the coastal ocean

Marine Group II Euryarchaeota (Candidatus Poseidoniales), abundant but yet-uncultivated members of marine microbial communities, are thought to be (photo)heterotrophs that metabolize dissolved organic matter (DOM), such as lipids and peptides. However, little is known about their transcriptional activity. We mapped reads from a metatranscriptomic time series collected at Sapelo Island (GA, USA) to metagenome-assembled genomes to determine the diversity of transcriptionally active Ca. Poseidoniales. Summer metatranscriptomes had the highest abundance of Ca. Poseidoniales transcripts, mostly from the O1 and O3 genera within Ca. Thalassarchaeaceae (MGIIb). In contrast, transcripts from fall and winter samples were predominantly from Ca. Poseidoniaceae (MGIIa). Genes encoding proteorhodopsin, membrane-bound pyrophosphatase, peptidase/proteases, and part of the ß-oxidation pathway were highly transcribed across abundant genera. Highly transcribed genes specific to Ca. Thalassarchaeaceae included xanthine/uracil permease and receptors for amino acid transporters. Enrichment of Ca. Thalassarchaeaceae transcript reads related to protein/peptide, nucleic acid, and amino acid transport and metabolism, as well as transcript depletion during dark incubations, provided further evidence of heterotrophic metabolism. Quantitative PCR analysis of South Atlantic Bight samples indicated consistently abundant Ca. Poseidoniales in nearshore and inshore waters. Together, our data suggest that Ca. Thalassarchaeaceae are important photoheterotrophs potentially linking DOM and nitrogen cycling in coastal waters.

Resource partitioning of phytoplankton metabolites that support bacterial heterotrophy

The communities of bacteria that assemble around marine microphytoplankton are predictably dominated by Rhodobacterales, Flavobacteriales, and families within the Gammaproteobacteria. Yet whether this consistent ecological pattern reflects the result of resource-based niche partitioning or resource competition requires better knowledge of the metabolites linking microbial autotrophs and heterotrophs in the surface ocean. We characterized molecules targeted for uptake by three heterotrophic bacteria individually co-cultured with a marine diatom using two strategies that vetted the exometabolite pool for biological relevance by means of bacterial activity assays: expression of diagnostic genes and net drawdown of exometabolites, the latter detected with mass spectrometry and nuclear magnetic resonance using novel sample preparation approaches. Of the more than 36 organic molecules with evidence of bacterial uptake, 53% contained nitrogen (including nucleosides and amino acids), 11% were organic sulfur compounds (including dihydroxypropanesulfonate and dimethysulfoniopropionate), and 28% were components of polysaccharides (including chrysolaminarin, chitin, and alginate). Overlap in phytoplankton-derived metabolite use by bacteria in the absence of competition was low, and only guanosine, proline, and N-acetyl-D-glucosamine were predicted to be used by all three. Exometabolite uptake pattern points to a key role for ecological resource partitioning in the assembly marine bacterial communities transforming recent photosynthate.

Ecological drivers of bacterial community assembly in synthetic phycospheres

In the nutrient-rich region surrounding marine phytoplankton cells, heterotrophic bacterioplankton transform a major fraction of recently fixed carbon through the uptake and catabolism of phytoplankton metabolites. We sought to understand the rules by which marine bacterial communities assemble in these nutrient-enhanced phycospheres, specifically addressing the role of host resources in driving community coalescence. Synthetic systems with varying combinations of known exometabolites of marine phytoplankton were inoculated with seawater bacterial assemblages, and communities were transferred daily to mimic the average duration of natural phycospheres. We found that bacterial community assembly was predictable from linear combinations of the taxa maintained on each individual metabolite in the mixture, weighted for the growth each supported. Deviations from this simple additive resource model were observed but also attributed to resource-based factors via enhanced bacterial growth when host metabolites were available concurrently. The ability of photosynthetic hosts to shape bacterial associates through excreted metabolites represents a mechanism by which microbiomes with beneficial effects on host growth could be recruited. In the surface ocean, resource-based assembly of host-associated communities may underpin the evolution and maintenance of microbial interactions and determine the fate of a substantial portion of Earth's primary production.

Correction: Single-cell genomics shedding light on marine Thaumarchaeota diversification

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Scientists' warning to humanity: microorganisms and climate change

In the Anthropocene, in which we now live, climate change is impacting most life on Earth. Microorganisms support the existence of all higher trophic life forms. To understand how humans and other life forms on Earth (including those we are yet to discover) can withstand anthropogenic climate change, it is vital to incorporate knowledge of the microbial 'unseen majority'. We must learn not just how microorganisms affect climate change (including production and consumption of greenhouse gases) but also how they will be affected by climate change and other human activities. This Consensus Statement documents the central role and global importance of microorganisms in climate change biology. It also puts humanity on notice that the impact of climate change will depend heavily on responses of microorganisms, which are essential for achieving an environmentally sustainable future.

Microdiversity and temporal dynamics of marine bacterial dimethylsulfoniopropionate genes

Dimethylsulfoniopropionate (DMSP) is an abundant organic sulfur metabolite produced by many phytoplankton species and degraded by bacteria via two distinct pathways with climate-relevant implications. We assessed the diversity and abundance of bacteria possessing these pathways in the context of phytoplankton community composition over a 3-week time period spanning September-October, 2014 in Monterey Bay, CA. The dmdA gene from the DMSP demethylation pathway dominated the DMSP gene pool and was harboured mostly by members of the alphaproteobacterial SAR11 clade and secondarily by the Roseobacter group, particularly during the second half of the study. Novel members of the DMSP-degrading community emerged from dmdA sequences recovered from metagenome assemblies and single-cell sequencing, including largely uncharacterized gammaproteobacteria and alphaproteobacteria taxa. In the DMSP cleavage pathway, the SAR11 gene dddK was the most abundant early in the study, but was supplanted by dddP over time. SAR11 members, especially those harbouring genes for both DMSP degradation pathways, had a strong positive relationship with the abundance of dinoflagellates, and DMSP-degrading gammaproteobacteria co-occurred with haptophytes. This in situ study of the drivers of DMSP fate in a coastal ecosystem demonstrates for the first time correlations between specific groups of bacterial DMSP degraders and phytoplankton taxa.

Transient O2 pulses direct Fe crystallinity and Fe(III)-reducer gene expression within a soil microbiome

BACKGROUND

Many environments contain redox transition zones, where transient oxygenation events can modulate anaerobic reactions that influence the cycling of iron (Fe) and carbon (C) on a global scale. In predominantly anoxic soils, this biogeochemical cycling depends on Fe mineral composition and the activity of mixed Fe(III)-reducer populations that may be altered by periodic pulses of molecular oxygen (O₂).

METHODS

We repeatedly exposed anoxic (4% H₂:96% N₂) suspensions of soil from the Luquillo Critical Zone Observatory to 1.05 × 10², 1.05 × 10³, and 1.05 × 10⁴ mmol O₂ kg^-1 soil h^-1 during pulsed oxygenation treatments. Metatranscriptomic analysis and ⁵⁷Fe Mössbauer spectroscopy were used to investigate changes in Fe(III)-reducer gene expression and Fe(III) crystallinity, respectively.

RESULTS

Slow oxygenation resulted in soil Fe-(oxyhydr)oxides of higher crystallinity (38.1 ± 1.1% of total Fe) compared to fast oxygenation (30.6 ± 1.5%, P < 0.001). Transcripts binning to the genomes of Fe(III)-reducers Anaeromyxobacter, Geobacter, and Pelosinus indicated significant differences in extracellular electron transport (e.g., multiheme cytochrome c, multicopper oxidase, and type-IV pilin gene expression), adhesion/contact (e.g., S-layer, adhesin, and flagellin gene expression), and selective microbial competition (e.g., bacteriocin gene expression) between the slow and fast oxygenation treatments during microbial Fe(III) reduction. These data also suggest that diverse Fe(III)-reducer functions, including cytochrome-dependent extracellular electron transport, are associated with type-III fibronectin domains. Additionally, the metatranscriptomic data indicate that Methanobacterium was significantly more active in the reduction of CO₂ to CH₄ and in the expression of class(III) signal peptide/type-IV pilin genes following repeated fast oxygenation compared to slow oxygenation.

CONCLUSIONS

This study demonstrates that specific Fe(III)-reduction mechanisms in mixed Fe(III)-reducer populations are uniquely sensitive to the rate of O₂ influx, likely mediated by shifts in soil Fe(III)-(oxyhydr)oxide crystallinity. Overall, we provide evidence that transient oxygenation events play an important role in directing anaerobic pathways within soil microbiomes, which is expected to alter Fe and C cycling in redox-dynamic environments.

Identifying labile DOM components in a coastal ocean through depleted bacterial transcripts and chemical signals

Understanding which compounds comprising the complex and dynamic marine dissolved organic matter (DOM) pool are important in supporting heterotrophic bacterial production remains a major challenge. We eliminated sources of labile phytoplankton products, advected terrestrial material and photodegradation products to coastal microbial communities by enclosing water samples in situ for 24 h in the dark. Bacterial genes for which expression decreased between the beginning and end of the incubation and chemical formulae that were depleted over this same time frame were used as indicators of bioavailable compounds, an approach that avoids augmenting or modifying the natural DOM pool. Transport- and metabolism-related genes whose relative expression decreased implicated osmolytes, carboxylic acids, fatty acids, sugars and organic sulfur compounds as candidate bioreactive molecules. FT-ICR MS analysis of depleted molecular formulae implicated functional groups ~ 30-40 Da in size cleaved from semi-polar components of DOM as bioreactive components. Both gene expression and FT-ICR MS analyses indicated higher lability of compounds with sulfur and nitrogen heteroatoms. Untargeted methodologies able to integrate biological and chemical perspectives can be effective strategies for characterizing the labile microbial metabolites participating in carbon flux.

Mosaic patterns of B-vitamin synthesis and utilization in a natural marine microbial community

Aquatic environments contain large communities of microorganisms whose synergistic interactions mediate the cycling of major and trace nutrients, including vitamins. B-vitamins are essential coenzymes that many organisms cannot synthesize. Thus, their exchange among de novo synthesizers and auxotrophs is expected to play an important role in the microbial consortia and explain some of the temporal and spatial changes observed in diversity. In this study, we analyzed metatranscriptomes of a natural marine microbial community, diel sampled quarterly over one year to try to identify the potential major B-vitamin synthesizers and consumers. Transcriptomic data showed that the best-represented taxa dominated the expression of synthesis genes for some B-vitamins but lacked transcripts for others. For instance, Rhodobacterales dominated the expression of vitamin-B₁₂ synthesis, but not of vitamin-B₇ , whose synthesis transcripts were mainly represented by Flavobacteria. In contrast, bacterial groups that constituted less than 4% of the community (e.g., Verrucomicrobia) accounted for most of the vitamin-B₁ synthesis transcripts. Furthermore, ambient vitamin-B₁ concentrations were higher in samples collected during the day, and were positively correlated with chlorophyll-a concentrations. Our analysis supports the hypothesis that the mosaic of metabolic interdependencies through B-vitamin synthesis and exchange are key processes that contribute to shaping microbial communities in nature.

Ocean biogeochemistry modeled with emergent trait-based genomics

Marine ecosystem models have advanced to incorporate metabolic pathways discovered with genomic sequencing, but direct comparisons between models and "omics" data are lacking. We developed a model that directly simulates metagenomes and metatranscriptomes for comparison with observations. Model microbes were randomly assigned genes for specialized functions, and communities of 68 species were simulated in the Atlantic Ocean. Unfit organisms were replaced, and the model self-organized to develop community genomes and transcriptomes. Emergent communities from simulations that were initialized with different cohorts of randomly generated microbes all produced realistic vertical and horizontal ocean nutrient, genome, and transcriptome gradients. Thus, the library of gene functions available to the community, rather than the distribution of functions among specific organisms, drove community assembly and biogeochemical gradients in the model ocean.

Bacterial transcriptome remodeling during sequential co-culture with a marine dinoflagellate and diatom

In their role as primary producers, marine phytoplankton modulate heterotrophic bacterial activities through differences in the types and amounts of organic matter they release. This study investigates the transcriptional response of bacterium Ruegeria pomeroyi, a member of the Roseobacter clade known to affiliate with diverse phytoplankton groups in the ocean, during a shift in phytoplankton taxonomy. The bacterium was initially introduced into a culture of the dinoflagellate Alexandrium tamarense, and then experienced a change in phytoplankton community composition as the diatom Thalassiosira pseudonana gradually outcompeted the dinoflagellate. Samples were taken throughout the 30-day experiment to track shifts in bacterial gene expression informative of metabolic and ecological interactions. Transcriptome data indicate fundamental differences in the exometabolites released by the two phytoplankton. During growth with the dinoflagellate, gene expression patterns indicated that the main sources of carbon and energy for R. pomeroyi were dimethysulfoniopropionate (DMSP), taurine, methylated amines, and polyamines. During growth with the diatom, dihydroxypropanesulfonate (DHPS), xylose, ectoine, and glycolate instead appeared to fuel the bulk of bacterial metabolism. Expression patterns of genes for quorum sensing, gene transfer agent, and motility suggest that bacterial processes related to cell communication and signaling differed depending on which phytoplankton species dominated the co-culture. A remodeling of the R. pomeroyi transcriptome implicating more than a quarter of the genome occurred through the change in phytoplankton regime.

Bacterial Biogeography across the Amazon River-Ocean Continuum

Spatial and temporal patterns in microbial biodiversity across the Amazon river-ocean continuum were investigated along ∼675 km of the lower Amazon River mainstem, in the Tapajós River tributary, and in the plume and coastal ocean during low and high river discharge using amplicon sequencing of 16S rRNA genes in whole water and size-fractionated samples (0.2–2.0 μm and >2.0 μm). River communities varied among tributaries, but mainstem communities were spatially homogeneous and tracked seasonal changes in river discharge and co-varying factors. Co-occurrence network analysis identified strongly interconnected river assemblages during high (May) and low (December) discharge periods, and weakly interconnected transitional assemblages in September, suggesting that this system supports two seasonal microbial communities linked to river discharge. In contrast, plume communities showed little seasonal differences and instead varied spatially tracking salinity. However, salinity explained only a small fraction of community variability, and plume communities in blooms of diatom-diazotroph assemblages were strikingly different than those in other high salinity plume samples. This suggests that while salinity physically structures plumes through buoyancy and mixing, the composition of plume-specific communities is controlled by other factors including nutrients, phytoplankton community composition, and dissolved organic matter chemistry. Co-occurrence networks identified interconnected assemblages associated with the highly productive low salinity near-shore region, diatom-diazotroph blooms, and the plume edge region, and weakly interconnected assemblages in high salinity regions. This suggests that the plume supports a transitional community influenced by immigration of ocean bacteria from the plume edge, and by species sorting as these communities adapt to local environmental conditions. Few studies have explored patterns of microbial diversity in tropical rivers and coastal oceans. Comparison of Amazon continuum microbial communities to those from temperate and arctic systems suggest that river discharge and salinity are master variables structuring a range of environmental conditions that control bacterial communities across the river-ocean continuum.

Expression patterns of elemental cycling genes in the Amazon River Plume

Metatranscriptomics and metagenomics data sets benchmarked with internal standards were used to characterize the expression patterns for biogeochemically relevant bacterial and archaeal genes mediating carbon, nitrogen, phosphorus and sulfur uptake and metabolism through the salinity gradient of the Amazon River Plume. The genes were identified in 48 metatranscriptomic and metagenomic data sets summing to >500 million quality-controlled reads from six locations in the plume ecosystem. The ratio of transcripts per gene copy (a direct measure of expression made possible by internal standard additions) showed that the free-living bacteria and archaea exhibited only small changes in the expression levels of biogeochemically relevant genes through the salinity and nutrient zones of the plume. In contrast, the expression levels of genes in particle-associated cells varied over orders of magnitude among the stations, with the largest differences measured for genes mediating aspects of nitrogen cycling (nifH, amtB and amoA) and phosphorus acquisition (pstC, phoX and phoU). Taxa varied in their baseline gene expression levels and extent of regulation, and most of the spatial variation in the expression level could be attributed to changes in gene regulation after removing the effect of shifting taxonomic composition. We hypothesize that changes in microbial element cycling along the Amazon River Plume are largely driven by shifting activities of particle-associated cells, with most activities peaking in the mesohaline regions where N₂ fixation rates are elevated.

Spontaneous mutations of a model heterotrophic marine bacterium

Heterotrophic marine bacterioplankton populations display substantive genomic diversity that is commonly explained to be the result of selective forces imposed by resource limitation or interactions with phage and predators. Here we use a mutation-accumulation experiment followed by whole-genome sequencing of mutation lines to determine an unbiased rate and molecular spectrum of spontaneous mutations for a model heterotrophic marine bacterium in the globally important Roseobacter clade, Ruegeria pomeroyi DSS-3. We find evidence for mutational bias towards deletions over insertions, and this process alone could account for a sizable portion of genome size diversity among roseobacters and also implies that lateral gene transfer and/or selection must also play a role in maintaining roseobacters with large genome sizes. We also find evidence for a mutational bias in favor of changes from A/T to G/C nucleobases, which explains widespread occurrences of G/C-enriched Roseobacter genomes. Using the calculated mutation rate of 1.39 × 10^-10 per base per generation, we implement a 'mutation-rate clock' approach to date the evolution of roseobacters by assuming a constant mutation rate along their evolutionary history. This approach gives an estimated date of Roseobacter genome expansion in good agreement with an earlier fossil-based estimate of ~250 million years ago and is consistent with a hypothesis of a correlated evolutionary history between roseobacters and marine eukaryotic phytoplankton groups.

Small RNAs expressed during dimethylsulfoniopropionate degradation by a model marine bacterium

The fate of the sulfur moiety of dimethylsulfoniopropionate (DMSP) depends on the 'bacterial switch', a regulatory point between two metabolic pathways with different biogeochemical endpoints. Studies have focused on transcriptional patterns of known genes to determine physiological and environmental factors affecting this switch, but post-transcriptional regulation has been under-studied. Here we use a model bacterium containing both pathways to look for transcription of non-coding regulatory small RNAs (sRNAs) during DMSP metabolism. RNA-seq analysis of Ruegeria pomeroyi DSS-3 grown with DMSP, metabolic intermediates of DMSP degradation (MMPA or acetate), or methionine revealed 182 putative sRNAs, with 46 showing differential expression during growth on DMSP. A knockout mutant constructed for an upregulated sRNA had a phenotype that differed in its use of the two degradation pathways. Because transcription patterns of many differentially expressed sRNAs were not correlated with the transcription of their putative target gene, their effects on DMSP degradation would not be observable in the transcriptome. Overall, our results indicate that sRNAs are crucial but largely cryptic actors in regulating DMSP metabolism in this model marine bacterium and potentially other bacterial groups involved in the surface ocean sulfur cycle.

Low genome content diversity of marine planktonic Thaumarchaeota

Members of Thaumarchaeota are responsible for much of the ammonia oxidation occurring in the ocean. Recent studies showed that marine Thaumarchaeota have versatile metabolic capabilities, but sequencing additional genomes has not significantly increased the gene content ascribed to this group. We used the assembly-free dN pipeline software in combination with phylogenetic analyses to interrogate shotgun metagenomic data sets to gain a better understanding of the genomic diversity of Thaumarchaeota populations. The program confidently assigned ∼3,000 paired-end reads to Thaumarchaeota, independent of homologies to any known Thaumarchaeota genome sequence. Only 2% of these reads potentially harbor new genes that were absent from the genome of 'Candidatus Nitrosopumilus maritimus' str. SCM1, even though this strain was isolated from a marine aquarium rather than directly from the ocean. One of these novel genes encode proteins associated with the CRISPR/Cas system, Cas1, suggesting that phage defense through CRISPR may be also present in planktonic Thaumarchaeota lineages. Our results suggest that marine Thaumarchaeota populations have very low diversity in genome content, which is corroborated using computer simulation analyses of two bacterial lineages with known genome content diversity.

Experimental Identification of Small Non-Coding RNAs in the Model Marine Bacterium Ruegeria pomeroyi DSS-3

In oligotrophic ocean waters where bacteria are often subjected to chronic nutrient limitation, community transcriptome sequencing has pointed to the presence of highly abundant small RNAs (sRNAs). The role of sRNAs in regulating response to nutrient stress was investigated in a model heterotrophic marine bacterium Ruegeria pomeroyi grown in continuous culture under carbon (C) and nitrogen (N) limitation. RNAseq analysis identified 99 putative sRNAs. Sixty-nine were cis-encoded and located antisense to a presumed target gene. Thirty were trans-encoded and initial target prediction was performed computationally. The most prevalent functional roles of genes anti-sense to the cis-sRNAs were transport, cell-cell interactions, signal transduction, and transcriptional regulation. Most sRNAs were transcribed equally under both C and N limitation, and may be involved in a general stress response. However, 14 were regulated differentially between the C and N treatments and may respond to specific nutrient limitations. A network analysis of the predicted target genes of the R. pomeroyi cis-sRNAs indicated that they average fewer connections than typical protein-encoding genes, and appear to be more important in peripheral or niche-defining functions encoded in the pan genome.

Metagenomic and metatranscriptomic inventories of the lower Amazon River, May 2011

BACKGROUND

The Amazon River runs nearly 6500 km across the South American continent before emptying into the western tropical North Atlantic Ocean. In terms of both volume and watershed area, it is the world's largest riverine system, affecting elemental cycling on a global scale.

RESULTS

A quantitative inventory of genes and transcripts benchmarked with internal standards was obtained at five stations in the lower Amazon River during May 2011. At each station, metagenomes and metatranscriptomes were obtained in duplicate for two microbial size fractions (free-living, 0.2 to 2.0 μm; particle-associated, 2.0 to 297 μm) using 150 × 150 paired-end Illumina sequencing. Forty eight sample datasets were obtained, averaging 15 × 10(6) potential protein-encoding reads each (730 × 10(6) total). Prokaryotic metagenomes and metatranscriptomes were dominated by members of the phyla Actinobacteria, Planctomycetes, Betaproteobacteria, Verrucomicrobia, Nitrospirae, and Acidobacteria. The actinobacterium SCGC AAA027-L06 reference genome recruited the greatest number of reads overall, with this single bin contributing an average of 50 billion genes and 500 million transcripts per liter of river water. Several dominant taxa were unevenly distributed between the free-living and particle-associated size fractions, such as a particle-associated bias for reads binning to planctomycete Schlesneria paludicola and a free-living bias for actinobacterium SCGC AAA027-L06. Gene expression ratios (transcripts to gene copy ratio) increased downstream from Óbidos to Macapá and Belém, indicating higher per cell activity of Amazon River bacteria and archaea as river water approached the ocean.

CONCLUSION

This inventory of riverine microbial genes and transcripts, benchmarked with internal standards for full quantitation, provides an unparalleled window into microbial taxa and functions in the globally important Amazon River ecosystem.

How do divergent ecological strategies emerge among marine bacterioplankton lineages?

Heterotrophic bacteria in pelagic marine environments are frequently categorized into two canonical ecological groups: patch-associated and free-living. This framework provides a conceptual basis for understanding bacterial utilization of oceanic organic matter. Some patch-associated bacteria are ecologically linked with eukaryotic phytoplankton, and this observation fits with predicted coincidence of their genome expansion with marine phytoplankton diversification. By contrast, free-living bacteria in today's oceans typically live singly with streamlined metabolic and regulatory functions that allow them to grow in nutrient-poor seawater. Recent analyses of marine Alphaproteobacteria suggest that some free-living bacterioplankton lineages evolved from patch-associated ancestors up to several hundred million years ago. While evolutionary analyses agree with the hypothesis that natural selection has maintained these distinct ecological strategies and genomic traits in present-day populations, they do not rule out a major role for genetic drift in driving ancient ecological switches. These two evolutionary forces may have acted on ocean bacteria at different geological time scales and under different geochemical constraints, with possible implications for future adaptations to a changing ocean. New evolutionary models and genomic data are leading to a more comprehensive understanding of marine bacterioplankton evolutionary history.

Phenotypic plasticity in heterotrophic marine microbial communities in continuous cultures

Phenotypic plasticity (PP) is the development of alternate phenotypes of a given taxon as an adaptation to environmental conditions. Methodological limitations have restricted the quantification of PP to the measurement of a few traits in single organisms. We used metatranscriptomic libraries to overcome these challenges and estimate PP using the expressed genes of multiple heterotrophic organisms as a proxy for traits in a microbial community. The metatranscriptomes captured the expression response of natural marine bacterial communities grown on differing carbon resource regimes in continuous cultures. We found that taxa with different magnitudes of PP coexisted in the same cultures, and that members of the order Rhodobacterales had the highest levels of PP. In agreement with previous studies, our results suggest that continuous culturing may have specifically selected for taxa featuring a rather high range of PP. On average, PP and abundance changes within a taxon contributed equally to the organism's change in functional gene abundance, implying that both PP and abundance mediated observed differences in community function. However, not all functional changes due to PP were directly reflected in the bulk community functional response: gene expression changes in individual taxa due to PP were partly masked by counterbalanced expression of the same gene in other taxa. This observation demonstrates that PP had a stabilizing effect on a community's functional response to environmental change.

An Updated genome annotation for the model marine bacterium Ruegeria pomeroyi DSS-3

When the genome of Ruegeria pomeroyi DSS-3 was published in 2004, it represented the first sequence from a heterotrophic marine bacterium. Over the last ten years, the strain has become a valuable model for understanding the cycling of sulfur and carbon in the ocean. To ensure that this genome remains useful, we have updated 69 genes to incorporate functional annotations based on new experimental data, and improved the identification of 120 protein-coding regions based on proteomic and transcriptomic data. We review the progress made in understanding the biology of R. pomeroyi DSS-3 and list the changes made to the genome.

The Marine Microbial Eukaryote Transcriptome Sequencing Project (MMETSP): illuminating the functional diversity of eukaryotic life in the oceans through transcriptome sequencing

Current sampling of genomic sequence data from eukaryotes is relatively poor, biased, and inadequate to address important questions about their biology, evolution, and ecology; this Community Page describes a resource of 700 transcriptomes from marine microbial eukaryotes to help understand their role in the world's oceans.

Linking activity and function to ecosystem dynamics in a coastal bacterioplankton community

For bacterial communities containing hundreds to thousands of distinct populations, connecting functional processes and environmental dynamics at high taxonomic resolution has remained challenging. Here we use the expression of ribosomal proteins (%RP) as a proxy for in situ activity of 200 taxa within 20 metatranscriptomic samples in a coastal ocean time series encompassing both seasonal variability and diel dynamics. %RP patterns grouped the taxa into seven activity clusters with distinct profiles in functional gene expression and correlations with environmental gradients. Clusters 1-3 had their highest potential activity in the winter and fall, and included some of the most active taxa, while Clusters 4-7 had their highest potential activity in the spring and summer. Cluster 1 taxa were characterized by gene expression for motility and complex carbohydrate degradation (dominated by Gammaproteobacteria and Bacteroidetes), and Cluster 2 taxa by transcription of genes for amino acid and aromatic compound metabolism and aerobic anoxygenic phototrophy (Roseobacter). Other activity clusters were enriched in transcripts for proteorhodopsin and methylotrophy (Cluster 4; SAR11 and methylotrophs), photosynthesis and attachment (Clusters 5 and 7; Synechococcus, picoeukaryotes, Verucomicrobia, and Planctomycetes), and sulfur oxidation (Cluster 7; Gammaproteobacteria). The seasonal patterns in activity were overlain, and sometimes obscured, by large differences in %RP over shorter day-night timescales. Seventy-eight taxa, many of them heterotrophs, had a higher %RP activity index during the day than night, indicating a strong diel activity rhythm at this coastal site. Emerging from these taxonomically- and time-resolved estimates of in situ microbial activity are predictions of specific ecological groupings of microbial taxa in a dynamic coastal environment.

Comparing effective population sizes of dominant marine alphaproteobacteria lineages

A fundamental question in marine microbial ecology is how microbes adapt to ocean environments. Although numerically dominant populations are typically considered more successful, higher census population sizes (Nc) do not equate directly to a greater capability for adaptation. Instead, effective population size (Ne) determines the fate of deleterious and favourable mutations, and thus is a key parameter for determining the adaptive potential of a population. In the case of the SAR11 and Roseobacter lineages, two abundant heterotrophic bacteria in ocean surface waters with contrasting life history strategies, culture-independent population surveys suggest that SAR11s have greater Nc than Roseobacters. To determine relative Ne, we compared the ratio of nonsynonymous to synonymous substitution rates (ω) of recently diverged lineages of these taxa. Values of ω associated with several of the Roseobacter subclades were lower than for SAR11 subclades, suggesting greater Ne in these cases. Most Roseobacter lineages also had smaller ω values compared with an atypical basal Roseobacter lineage with a large Nc. This finding provides insight into variability in Ne across two important marine bacterial lineages, and provides an evolutionary context for considering how heterotrophic marine bacteria may differ in their ability to adapt to changing ocean habitats.

Evolutionary analysis of a streamlined lineage of surface ocean Roseobacters

The vast majority of surface ocean bacteria are uncultivated. Compared with their cultured relatives, they frequently exhibit a streamlined genome, reduced G+C content and distinct gene repertoire. These genomic traits are relevant to environmental adaptation, and have generally been thought to become fixed in marine bacterial populations through selection. Using single-cell genomics, we sequenced four uncultivated cells affiliated with the ecologically relevant Roseobacter clade and used a composition-heterogeneous Bayesian phylogenomic model to resolve these single-cell genomes into a new clade. This lineage has no representatives in culture, yet accounts for ∼35% of Roseobacters in some surface ocean waters. Analyses of multiple genomic traits, including genome size, G+C content and percentage of noncoding DNA, suggest that these single cells are representative of oceanic Roseobacters but divergent from isolates. Population genetic analyses showed that substitution of physicochemically dissimilar amino acids and replacement of G+C-rich to G+C-poor codons are accelerated in the uncultivated clade, processes that are explained equally well by genetic drift as by the more frequently invoked explanation of natural selection. The relative importance of drift vs selection in this clade, and perhaps in other marine bacterial clades with streamlined G+C-poor genomes, remains unresolved until more evidence is accumulated.

Quantitative microbial metatranscriptomics

The direct retrieval and sequencing of environmental RNA is emerging as a powerful technique to elucidate the in situ activities of microbial communities. Here we provide a metatranscriptomic protocol describing environmental sample collection, rRNA depletion, mRNA amplification, cDNA synthesis, and bioinformatic analysis. In addition, the preparation of internal RNA standards and their addition to the sample are discussed, providing a method by which transcript numbers can be expressed as absolute abundances in the environment and more readily compared to other biogeochemical and ecological measurements.

The Amazon continuum dataset: quantitative metagenomic and metatranscriptomic inventories of the Amazon River plume, June 2010

BACKGROUND

The Amazon River is by far the world's largest in terms of volume and area, generating a fluvial export that accounts for about a fifth of riverine input into the world's oceans. Marine microbial communities of the Western Tropical North Atlantic Ocean are strongly affected by the terrestrial materials carried by the Amazon plume, including dissolved (DOC) and particulate organic carbon (POC) and inorganic nutrients, with impacts on primary productivity and carbon sequestration.

RESULTS

We inventoried genes and transcripts at six stations in the Amazon River plume during June 2010. At each station, internal standard-spiked metagenomes, non-selective metatranscriptomes, and poly(A)-selective metatranscriptomes were obtained in duplicate for two discrete size fractions (0.2 to 2.0 μm and 2.0 to 156 μm) using 150 × 150 paired-end Illumina sequencing. Following quality control, the dataset contained 360 million reads of approximately 200 bp average size from Bacteria, Archaea, Eukarya, and viruses. Bacterial metagenomes and metatranscriptomes were dominated by Synechococcus, Prochlorococcus, SAR11, SAR116, and SAR86, with high contributions from SAR324 and Verrucomicrobia at some stations. Diatoms, green picophytoplankton, dinoflagellates, haptophytes, and copepods dominated the eukaryotic genes and transcripts. Gene expression ratios differed by station, size fraction, and microbial group, with transcription levels varying over three orders of magnitude across taxa and environments.

CONCLUSIONS

This first comprehensive inventory of microbial genes and transcripts, benchmarked with internal standards for full quantitation, is generating novel insights into biogeochemical processes of the Amazon plume and improving prediction of climate change impacts on the marine biosphere.

Single-cell genomics shedding light on marine Thaumarchaeota diversification

Previous studies based on analysis of amoA, 16S ribosomal RNA or accA gene sequences have established that marine Thaumarchaeota fall into two phylogenetically distinct groups corresponding to shallow- and deep-water clades, but it is not clear how water depth interacts with other environmental factors, including light, temperature and location, to affect this pattern of diversification. Earlier studies focused on single-gene distributions were not able to link phylogenetic structure to other aspects of functional adaptation. Here, we analyzed the genome content of 46 uncultivated single Thaumarchaeota cells sampled from epi- and mesopelagic waters of subtropical, temperate and polar oceans. Phylogenomic analysis showed that populations diverged by depth, as expected, and that mesopelagic populations from different locations were well mixed. Functional analysis showed that some traits, including putative DNA photolyase and catalase genes that may be related to adaptive mechanisms to reduce light-induced damage, were found exclusively in members of the epipelagic clade. Our analysis of partial genomes has thus confirmed the depth differentiation of Thaumarchaeota populations observed previously, consistent with the distribution of putative mechanisms to reduce light-induced damage in shallow- and deep-water populations.

Seasonal variation in the metatranscriptomes of a Thaumarchaeota population from SE USA coastal waters

We used a combination of metatranscriptomic analyses and quantitative PCR (qPCR) to study seasonal changes in Thaumarchaeota populations from a salt marsh-dominated estuary. Surface waters (0.5 m depth) were sampled quarterly at Marsh Landing, Sapelo Island, GA, USA over a 3-year period. We found a mid-summer peak in Thaumarchaeota abundance measured by qPCR of either 16S rRNA or amoA genes in each of the 3 years. Thaumarchaeota were 100-1000-fold more abundant during the peak than at other times of the year, whereas the abundance of ammonia- and nitrite-oxidizing Bacteria varied <10-fold over the same period. Analysis of the microdiversity of several highly transcribed genes in 20 metatranscriptomes from a 1-year subset of these samples showed that the transcriptionally active population consisted of 2 or 3 dominant phylotypes that differed between successive summers. This shift appeared to have begun during the preceding winter and spring. Transcripts from the same genes dominated the Thaumarchaeota mRNA pool throughout the year, with genes encoding proteins believed to be involved in nitrogen uptake and oxidation, and two hypothetical proteins being the most abundant transcripts in all libraries. Analysis of individual genes over the seasonal cycle suggested that transcription was tied more closely to variation in growth rates than to seasonal changes in environmental conditions. Day-night differences in the relative abundance of transcripts for ribosomal proteins suggested diurnal variation in Thaumarchaeota growth.

Assembly-free metagenomic analysis reveals new metabolic capabilities in surface ocean bacterioplankton

Uncovering the metabolic capabilities of microbes is key to understanding global energy flux and nutrient transformations. Since the vast majority of environmental microorganisms are uncultured, metagenomics has become an important tool to genotype the microbial community. This study uses a recently developed computational method to confidently assign metagenomic reads to microbial clades without the requirement of metagenome assembly by comparing the evolutionary pattern of nucleotide sequences at non-synonymous sites between metagenomic and orthologous reference genes. We found evidence for new, ecologically relevant metabolic pathways in several lineages of surface ocean bacterioplankton using the Global Ocean Survey (GOS) metagenomic data, including assimilatory sulfate reduction and alkaline phosphatase capabilities in the alphaproteobacterial SAR11 clade, and proteorhodopsin-like genes in the cyanobacterial genus Prochlorococcus. These findings raise new hypotheses about microbial roles in energy flux and organic matter transformation in the ocean.

Transcriptional response of bathypelagic marine bacterioplankton to the Deepwater Horizon oil spill

The Deepwater Horizon blowout released a massive amount of oil and gas into the deep ocean between April and July 2010, stimulating microbial blooms of petroleum-degrading bacteria. To understand the metabolic response of marine microorganisms, we sequenced ≈ 66 million community transcripts that revealed the identity of metabolically active microbes and their roles in petroleum consumption. Reads were assigned to reference genes from ≈ 2700 bacterial and archaeal taxa, but most assignments (39%) were to just six genomes representing predominantly methane- and petroleum-degrading Gammaproteobacteria. Specific pathways for the degradation of alkanes, aromatic compounds and methane emerged from the metatranscriptomes, with some transcripts assigned to methane monooxygenases representing highly divergent homologs that may degrade either methane or short alkanes. The microbial community in the plume was less taxonomically and functionally diverse than the unexposed community below the plume; this was due primarily to decreased species evenness resulting from Gammaproteobacteria blooms. Surprisingly, a number of taxa (related to SAR11, Nitrosopumilus and Bacteroides, among others) contributed equal numbers of transcripts per liter in both the unexposed and plume samples, suggesting that some groups were unaffected by the petroleum inputs and blooms of degrader taxa, and may be important for re-establishing the pre-spill microbial community structure.

Prokaryotic Super Program Advisory Committee DOE Joint Genome Institute, Walnut Creek, CA, March 27, 2013

The Prokaryotic Super Program Advisory Committee met on March 27, 2013 for their annual review the Prokaryotic Super Program at the DOE Joint Genome Institute. As is the case with any site visit or program review, the objective is to evaluate progress in meeting organizational objectives, provide feedback to from the user-community and to assist the JGI in formulating plans for the coming year. The advisors want to commend the JGI for its central role in developing new technologies and capabilities, and for catalyzing the formation of new collaborative user communities. Highlights of the post-meeting exchanges among the advisors focused on the importance of programmatic initiatives including:

• GEBA, which serves as a phylogenetic “base-map” on which our knowledge of functional diversity can be layered.

• FEBA, which promises to provide new insights into the physiological capabilities of prokaryotes under highly standardized conditions.

• Single-cell genomics technology, which is seen to significantly enhance our ability to interpret genomic and metagenomic data and broaden the scope of the GEBA program to encompass at least a part of the microbial “dark-matter”.

• IMG, which is seen to play a central role in JGI programs and is viewed as a strategically important asset in the JGI portfolio.

On this latter point, the committee encourages the formation of a strategic relationship between IMG and the Kbase to ensure that the intelligence, deep knowledge and experience captured in the former is not lost. The committee strongly urges the DOE to continue its support for maintaining this critical resource.