Digital Microbe: a genome-informed data integration framework for team science on emerging model organisms.

The remarkable pace of genomic data generation is rapidly transforming our understanding of life at the micron scale. Yet this data stream also creates challenges for team science. A single microbe can have multiple versions of genome architecture, functional gene annotations, and gene identifiers; additionally, the lack of mechanisms for collating and preserving advances in this knowledge raises barriers to community coalescence around shared datasets. "Digital Microbes" are frameworks for interoperable and reproducible collaborative science through open source, community-curated data packages built on a (pan)genomic foundation. Housed within an integrative software environment, Digital Microbes ensure real-time alignment of research efforts for collaborative teams and facilitate novel scientific insights as new layers of data are added. Here we describe two Digital Microbes: 1) the heterotrophic marine bacterium Ruegeria pomeroyi DSS-3 with > 100 transcriptomic datasets from lab and field studies, and 2) the pangenome of the cosmopolitan marine heterotroph Alteromonas containing 339 genomes. Examples demonstrate how an integrated framework collating public (pan)genome-informed data can generate novel and reproducible findings.

A mutant fitness assay identifies bacterial interactions in a model ocean hot spot.

Bacteria that assemble in phycospheres surrounding living phytoplankton cells metabolize a substantial proportion of ocean primary productivity. Yet the type and extent of interactions occurring among species that colonize these micron-scale "hot spot" environments are challenging to study. We identified genes that mediate bacterial interactions in phycosphere communities by culturing a transposon mutant library of copiotrophic bacterium Ruegeria pomeroyi DSS-3 with the diatom Thalassiosira pseudonana CCMP1335 as the sole source of organic matter in the presence or absence of other heterotrophic bacterial species. The function of genes having significant effects on R. pomeroyi fitness indicated explicit cell-cell interactions initiated in the multibacterial phycospheres. We found that R. pomeroyi simultaneously competed for shared substrates while increasing reliance on substrates that did not support the other species' growth. Fitness outcomes also indicated that the bacterium competed for nitrogen in the forms of ammonium and amino acids; obtained purines, pyrimidines, and cofactors via crossfeeding; both initiated and defended antagonistic interactions; and sensed an environment with altered oxygen and superoxide levels. The large genomes characteristic of copiotrophic marine bacteria are hypothesized to enable responses to dynamic ecological challenges occurring at the scale of microns. Here, we discover >200 nonessential genes implicated in the management of fitness costs and benefits of membership in a globally significant bacterial community.

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.

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.

Genome Sequences and Metagenome-Assembled Genome Sequences of Microbial Communities Enriched on Phytoplankton Exometabolites.

We report 11 bacterial draft genome sequences and 38 metagenome-assembled genomes (MAGs) from marine phytoplankton exometabolite enrichments. The genomes and MAGs represent marine bacteria adapted to the metabolite environment of phycospheres, organic matter-rich regions surrounding phytoplankton cells, and are useful for exploring functional and taxonomic attributes of phytoplankton-associated bacterial communities.

Microbial metagenomes and metatranscriptomes during a coastal phytoplankton bloom.

Metagenomic and metatranscriptomic time-series data covering a 52-day period in the fall of 2016 provide an inventory of bacterial and archaeal community genes, transcripts, and taxonomy during an intense dinoflagellate bloom in Monterey Bay, CA, USA. The dataset comprises 84 metagenomes (0.8 terabases), 82 metatranscriptomes (1.1 terabases), and 88 16S rRNA amplicon libraries from samples collected on 41 dates. The dataset also includes 88 18S rRNA amplicon libraries, characterizing the taxonomy of the eukaryotic community during the bloom. Accompanying the sequence data are chemical and biological measurements associated with each sample. These datasets will facilitate studies of the structure and function of marine bacterial communities during episodic phytoplankton blooms.

Recognition cascade and metabolite transfer in a marine bacteria-phytoplankton model system.

The trophic linkage between marine bacteria and phytoplankton in the surface ocean is a key step in the global carbon cycle, with almost half of marine primary production transformed by heterotrophic bacterioplankton within hours to weeks of fixation. Early studies conceptualized this link as the passive addition and removal of organic compounds from a shared seawater reservoir. Here, we analysed transcript and intracellular metabolite patterns in a two-member model system and found that the presence of a heterotrophic bacterium induced a potential recognition cascade in a marine phytoplankton species that parallels better-understood vascular plant response systems. Bacterium Ruegeria pomeroyi DSS-3 triggered differential expression of >80 genes in diatom Thalassiosira pseudonana CCMP1335 that are homologs to those used by plants to recognize external stimuli, including proteins putatively involved in leucine-rich repeat recognition activity, second messenger production and protein kinase cascades. Co-cultured diatoms also downregulated lipid biosynthesis genes and upregulated chitin metabolism genes. From differential expression of bacterial transporter systems, we hypothesize that nine diatom metabolites supported the majority of bacterial growth, among them sulfonates, sugar derivatives and organic nitrogen compounds. Similar recognition responses and metabolic linkages as observed in this model system may influence carbon transformations by ocean plankton.

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 N2 fixation rates are elevated.

Patterns of Transcript Abundance of Eukaryotic Biogeochemically-Relevant Genes in the Amazon River Plume.

The Amazon River has the largest discharge of all rivers on Earth, and its complex plume system fuels a wide array of biogeochemical processes, across a large area of the western tropical North Atlantic. The plume thus stimulates microbial processes affecting carbon sequestration and nutrient cycles at a global scale. Chromosomal gene expression patterns of the 2.0 to 156 µm size-fraction eukaryotic microbial community were investigated in the Amazon River Plume, generating a robust dataset (more than 100 million mRNA sequences) that depicts the metabolic capabilities and interactions among the eukaryotic microbes. Combining classical oceanographic field measurements with metatranscriptomics yielded characterization of the hydrographic conditions simultaneous with a quantification of transcriptional activity and identity of the community. We highlight the patterns of eukaryotic gene expression for 31 biogeochemically significant gene targets hypothesized to be valuable within forecasting models. An advantage to this targeted approach is that the database of reference sequences used to identify the target genes was selectively constructed and highly curated optimizing taxonomic coverage, throughput, and the accuracy of annotations. A coastal diatom bloom highly expressed nitrate transporters and carbonic anhydrase presumably to support high growth rates and enhance uptake of low levels of dissolved nitrate and CO2. Diatom-diazotroph association (DDA: diatoms with nitrogen fixing symbionts) blooms were common when surface salinity was mesohaline and dissolved nitrate concentrations were below detection, and hence did not show evidence of nitrate utilization, suggesting they relied on ammonium transporters to aquire recently fixed nitrogen. These DDA blooms in the outer plume had rapid turnover of the photosystem D1 protein presumably caused by photodegradation under increased light penetration in clearer waters, and increased expression of silicon transporters as silicon became limiting. Expression of these genes, including carbonic anhydrase and transporters for nitrate and phosphate, were found to reflect the physiological status and biogeochemistry of river plume environments. These relatively stable patterns of eukaryotic transcript abundance occurred over modest spatiotemporal scales, with similarity observed in sample duplicates collected up to 2.45 km in space and 120 minutes in time. These results confirm the use of metatranscriptomics as a valuable tool to understand and predict microbial community function.

Erratum to: An updated genome annotation for the model marine bacterium Ruegeria pomeroyi DSS-3.

[This corrects the article DOI: 10.1186/1944-3277-9-11.].

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.

Cryptic carbon and sulfur cycling between surface ocean plankton.

About half the carbon fixed by phytoplankton in the ocean is taken up and metabolized by marine bacteria, a transfer that is mediated through the seawater dissolved organic carbon (DOC) pool. The chemical complexity of marine DOC, along with a poor understanding of which compounds form the basis of trophic interactions between bacteria and phytoplankton, have impeded efforts to identify key currencies of this carbon cycle link. Here, we used transcriptional patterns in a bacterial-diatom model system based on vitamin B12 auxotrophy as a sensitive assay for metabolite exchange between marine plankton. The most highly up-regulated genes (up to 374-fold) by a marine Roseobacter clade bacterium when cocultured with the diatom Thalassiosira pseudonana were those encoding the transport and catabolism of 2,3-dihydroxypropane-1-sulfonate (DHPS). This compound has no currently recognized role in the marine microbial food web. As the genes for DHPS catabolism have limited distribution among bacterial taxa, T. pseudonana may use this sulfonate for targeted feeding of beneficial associates. Indeed, DHPS was both a major component of the T. pseudonana cytosol and an abundant microbial metabolite in a diatom bloom in the eastern North Pacific Ocean. Moreover, transcript analysis of the North Pacific samples provided evidence of DHPS catabolism by Roseobacter populations. Other such biogeochemically important metabolites may be common in the ocean but difficult to discriminate against the complex chemical background of seawater. Bacterial transformation of this diatom-derived sulfonate represents a previously unidentified and likely sizeable link in both the marine carbon and sulfur cycles.

Microspatial gene expression patterns in the Amazon River Plume.

We investigated expression of genes mediating elemental cycling at the microspatial scale in the ocean's largest river plume using, to our knowledge, the first fully quantitative inventory of genes and transcripts. The bacterial and archaeal communities associated with a phytoplankton bloom in Amazon River Plume waters at the outer continental shelf in June 2010 harbored ∼ 1.0 × 10(13) genes and 4.7 × 10(11) transcripts per liter that mapped to several thousand microbial genomes. Genomes from free-living cells were more abundant than those from particle-associated cells, and they generated more transcripts per liter for carbon fixation, heterotrophy, nitrogen and phosphorus uptake, and iron acquisition, although they had lower expression ratios (transcripts ⋅ gene(-1)) overall. Genomes from particle-associated cells contributed more transcripts for sulfur cycling, aromatic compound degradation, and the synthesis of biologically essential vitamins, with an overall twofold up-regulation of expression compared with free-living cells. Quantitatively, gene regulation differences were more important than genome abundance differences in explaining why microenvironment transcriptomes differed. Taxa contributing genomes to both free-living and particle-associated communities had up to 65% of their expressed genes regulated differently between the two, quantifying the extent of transcriptional plasticity in marine microbes in situ. In response to patchiness in carbon, nutrients, and light at the micrometer scale, Amazon Plume microbes regulated the expression of genes relevant to biogeochemical processes at the ecosystem scale.

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.

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.

Sizing up metatranscriptomics.

A typical marine bacterial cell in coastal seawater contains only ∼200 molecules of mRNA, each of which lasts only a few minutes before being degraded. Such a surprisingly small and dynamic cellular mRNA reservoir has important implications for understanding the bacterium's responses to environmental signals, as well as for our ability to measure those responses. In this perspective, we review the available data on transcript dynamics in environmental bacteria, and then consider the consequences of a small and transient mRNA inventory for functional metagenomic studies of microbial communities.

Transcriptional changes underlying elemental stoichiometry shifts in a marine heterotrophic bacterium.

Marine bacteria drive the biogeochemical processing of oceanic dissolved organic carbon (DOC), a 750-Tg C reservoir that is a critical component of the global C cycle. Catabolism of DOC is thought to be regulated by the biomass composition of heterotrophic bacteria, as cells maintain a C:N:P ratio of ∼50:10:1 during DOC processing. Yet a complicating factor in stoichiometry-based analyses is that bacteria can change the C:N:P ratio of their biomass in response to resource composition. We investigated the physiological mechanisms of resource-driven shifts in biomass stoichiometry in continuous cultures of the marine heterotrophic bacterium Ruegeria pomeroyi (a member of the Roseobacter clade) under four element limitation regimes (C, N, P, and S). Microarray analysis indicated that the bacterium scavenged for alternate sources of the scarce element when cells were C-, N-, or P-limited; reworked the ratios of biomolecules when C- and P- limited; and exerted tighter control over import/export and cytoplasmic pools when N-limited. Under S limitation, a scenario not existing naturally for surface ocean microbes, stress responses dominated transcriptional changes. Resource-driven changes in C:N ratios of up to 2.5-fold and in C:P ratios of up to sixfold were measured in R. pomeroyi biomass. These changes were best explained if the C and P content of the cells was flexible in the face of shifting resources but N content was not, achieved through the net balance of different transcriptional strategies. The cellular-level metabolic trade-offs that govern biomass stoichiometry in R. pomeroyi may have implications for global carbon cycling if extendable to other heterotrophic bacteria. Strong homeostatic responses to N limitation by marine bacteria would intensify competition with autotrophs. Modification of cellular inventories in C- and P-limited heterotrophs would vary the elemental ratio of particulate organic matter sequestered in the deep ocean.