Brandon M. Satinsky

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.

Microspatial gene expression patterns in the Amazon River Plume

Significance The microbial community of the Amazon River Plume determines the fate of the world’s largest input of terrestrial carbon and nutrients to the ocean. By benchmarking with internal standards during sample collection, we determined that each liter of plume seawater contains 1 trillion genes and 50 billion transcripts from thousands of bacterial, archaeal, and eukaryotic taxa. Gene regulation by taxa inhabiting distinct microenvironments provides insights into micron-scale patterns of transformations in the marine carbon, nitrogen, phosphorus, and sulfur cycles in this globally important ecosystem. 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 × 1013 genes and 4.7 × 1011 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.

Use of internal standards for quantitative metatranscriptome and metagenome analysis

Next generation sequencing-enabled metatranscriptomic and metagenomic datasets are providing unprecedented insights into the functional diversity of microbial communities, allowing detection of the genes present in a community as well as differentiation of those being actively transcribed. An emerging challenge of meta-omics approaches is how to quantitatively compare metagenomes and metatranscriptomes collected across spatial and temporal scales, or among treatments in experimental manipulations. Here, we describe the use of internal DNA and mRNA standards in meta-omics methodologies, and highlight how data collected in an absolute framework (per L or per cell) provides increased comparative power and insight into underlying causes of differences between samples.

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.

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.

Ocean biogeochemistry modeled with emergent trait-based genomics

Functional ocean biogeography Marine ecosystems are well represented in metagenomic and transcriptomic data. These data are not routinely used to test ecosystem models that explore ocean biogeography or biogeochemistry. Coles et al. built a model in which genes for a range of functions were assigned to different suites of simulated microbes (see the Perspective by Rynearson). Communities emerged from the model with realistic biogeographical and biogeochemical profiles when compared to microbial data collected from the Amazon River plume. However, functional composition trumped the details of taxonomy, and different, coevolving community compositions emerged that provided similar biogeochemical outcomes. Science, this issue p. 1149; see also p. 1129 Results of genomic-based microbe community modeling resemble real-world observations from the Amazon River plume. 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.

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.

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.

Patterns of Bacterial and Archaeal Gene Expression through the Lower Amazon River

Analysis of metatranscriptomic and metagenomic datasets from the lower reaches of the Amazon River between Obidos and the river mouth revealed microbial transcript and gene pools dominated by Actinobacteria, Thaumarchaeota, Bacteroidetes, Acidobacteria, Betaproteobacteria, and Planctomycetes. Three mainstem stations spanning a 625 km reach had similar gene expression patterns (transcripts gene copy-1) across a diverse suite of element cycling genes, but two tributary-influenced stations at the mouth of the Tapajos River and near the Tocantins River at Belem had distinct transcriptome composition and expression ratios, particularly for genes encoding light-related energy capture (higher) and iron acquisition and ammonia oxidation (lower). Environmental parameters that were useful predictors of gene expression ratios included concentrations of lignin phenols, suspended sediments, nitrate, phosphate, and particulate organic carbon and nitrogen. Similar to the gene expression data, these chemical properties reflected highly homogeneous mainstem stations punctuated by distinct tributary-influenced stations at Tapajos and Belem. Although heterotrophic processes were expected to dominate in the lower Amazon, transcripts from photosynthetic bacteria were abundant in tributary-influenced regions, and transcripts from Thaumarcheota taxa genetically capable of chemosynthetic ammonia oxidation accounted for up to 21% of the transcriptome at others. Based on regressions of transcript numbers against gene numbers, expression ratios of Thaumarchaeota populations were largely unchanged within the mainstem, suggesting a relatively minor role for gene regulation. These quantitative gene and transcript inventories detail a diverse array of energy acquisition strategies and metabolic capabilities for bacteria and archaea populations of the world’s largest river system.