Search

Unlike biologically available nitrogen and phosphorus, which are often at limiting concentrations in surface seawater, sulfur in the form of sulfate is plentiful and not considered to constrain marine microbial activity. Nonetheless, in a model system in which a marine bacterium obtains all of its carbon from co-cultured phytoplankton, bacterial gene expression suggests that at least seven dissolved organic sulfur (DOS) metabolites support bacterial heterotrophy. These labile exometabolites of marine dinoflagellates and diatoms include taurine, N-acetyltaurine, isethionate, choline-O-sulfate, cysteate, 2,3-dihydroxypropane-1-sulfonate (DHPS), and dimethylsulfoniopropionate (DMSP). Leveraging from the compounds identified in this model system, we assessed the role of sulfur metabolites in the ocean carbon cycle by mining the Tara Oceans dataset for diagnostic genes. In the 1.4 million bacterial genome equivalents surveyed, estimates of the frequency of genomes harboring the capability for DOS metabolite utilization ranged broadly, from only 1 out of every 190 genomes (for the C2 sulfonate isethionate) to 1 out of every 5 (for the sulfonium compound DMSP). Bacteria able to participate in DOS transformations are dominated by Alphaproteobacteria in the surface ocean, but by SAR324, Acidimicrobiia, and Gammaproteobacteria at mesopelagic depths, where the capability for utilization occurs in higher frequency than in surface bacteria for more than half the sulfur metabolites. The discovery of an abundant and diverse suite of marine bacteria with the genetic capacity for DOS transformation argues for an important role for sulfur metabolites in the pelagic ocean carbon cycle.

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

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.

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.

Resolving the ecological niches of coexisting marine microbial taxa is challenging due to the high species richness of microbial communities and the apparent functional redundancy in bacterial genomes and metagenomes. Here, we generated over 11 million Illumina reads of protein-encoding transcripts collected from well-mixed southeastern US coastal waters to characterize gene expression patterns distinguishing the ecological roles of hundreds of microbial taxa sharing the same environment. The taxa with highest in situ growth rates (based on relative abundance of ribosomal protein transcripts) were typically not the greatest contributors to community transcription, suggesting strong top-down ecological control, and their diverse transcriptomes indicated roles as metabolic generalists. The taxa with low in situ growth rates typically had low diversity transcriptomes dominated by specialized metabolisms. By identifying protein-encoding genes with atypically high expression for their level of conservation, unique functional roles of community members emerged related to substrate use (such as complex carbohydrates, fatty acids, methanesulfonate, taurine, tartrate, ectoine), alternative energy-conservation strategies (proteorhodopsin, AAnP, V-type pyrophosphatases, sulfur oxidation, hydrogen oxidation) and mechanisms for negotiating a heterogeneous environment (flagellar motility, gliding motility, adhesion strategies). On average, the heterotrophic bacterioplankton dedicated 7% of their transcriptomes to obtaining energy by non-heterotrophic means. This deep sequencing of a coastal bacterioplankton transcriptome provides the most highly resolved view of bacterioplankton niche dimensions yet available, uncovering a spectrum of unrecognized ecological strategies.

Quantitative PCR (qPCR) analysis revealed elevated relative abundance (1.8% of prokaryotes) of marine group 1 Crenarchaeota (MG1C) in two samples of southeastern US coastal bacterioplankton, collected in August 2008, compared with samples collected from the same site at different times (mean 0.026%). We analyzed the MG1C sequences in metatranscriptomes from these samples to gain an insight into the metabolism of MG1C population growing in the environment, and for comparison with ammonia-oxidizing bacteria (AOB) in the same samples. Assemblies revealed low diversity within sequences assigned to most individual MG1C open reading frames (ORFs) and high homology with 'Candidatus Nitrosopumilus maritimus' strain SCM1 genome sequences. Reads assigned to ORFs for ammonia uptake and oxidation accounted for 37% of all MG1C transcripts. We did not recover any reads for Nmar&#95;1354-Nmar&#95;1357, proposed to encode components of an alternative, nitroxyl-based ammonia oxidation pathway; however, reads from Nmar&#95;1259 and Nmar&#95;1667, annotated as encoding a multicopper oxidase with homology to nirK, were abundant. Reads assigned to two homologous ORFs (Nmar&#95;1201 and Nmar&#95;1547), annotated as hypothetical proteins were also abundant, suggesting that their unknown function is important to MG1C. Superoxide dismutase and peroxiredoxin-like transcripts were more abundant in the MG1C transcript pool than in the complete metatranscriptome, suggesting an enhanced response to oxidative stress by the MG1C population. qPCR indicated low AOB abundance (0.0010% of prokaryotes), and we found no transcripts related to ammonia oxidation and only one RuBisCO transcript among the transcripts assigned to AOB, suggesting they were not responding to the same environmental cues as the MG1C population.

The potential of metatranscriptomic sequencing to provide insights into the environmental factors that regulate microbial activities depends on how fully the sequence libraries capture community expression (that is, sample-sequencing depth and coverage depth), and the sensitivity with which expression differences between communities can be detected (that is, statistical power for hypothesis testing). In this study, we use an internal standard approach to make absolute (per liter) estimates of transcript numbers, a significant advantage over proportional estimates that can be biased by expression changes in unrelated genes. Coastal waters of the southeastern United States contain 1 × 10(12) bacterioplankton mRNA molecules per liter of seawater (~200 mRNA molecules per bacterial cell). Even for the large bacterioplankton libraries obtained in this study (~500,000 possible protein-encoding sequences in each of two libraries after discarding rRNAs and small RNAs from >1 million 454 FLX pyrosequencing reads), sample-sequencing depth was only 0.00001%. Expression levels of 82 genes diagnostic for transformations in the marine nitrogen, phosphorus and sulfur cycles ranged from below detection (<1 × 10(6) transcripts per liter) for 36 genes (for example, phosphonate metabolism gene phnH, dissimilatory nitrate reductase subunit napA) to >2.7 × 10(9) transcripts per liter (ammonia transporter amt and ammonia monooxygenase subunit amoC). Half of the categories for which expression was detected, however, had too few copy numbers for robust statistical resolution, as would be required for comparative (experimental or time-series) expression studies. By representing whole community gene abundance and expression in absolute units (per volume or mass of environment), 'omics' data can be better leveraged to improve understanding of microbially mediated processes in the ocean.

Dimethylsulfoniopropionate (DMSP) is an important source of reduced sulfur and carbon for marine microbial communities, as well as the precursor of the climate-active gas dimethylsulfide (DMS). In this study, we used metatranscriptomic sequencing to analyze gene expression profiles of a bacterial assemblage from surface waters at the Bermuda Atlantic Time-series Study (BATS) station with and without a short-term enrichment of DMSP (25 nM for 30 min). An average of 303 143 reads were obtained per treatment using 454 pyrosequencing technology, of which 51% were potential protein-encoding sequences. Transcripts from Gammaproteobacteria and Bacteroidetes increased in relative abundance on DMSP addition, yet there was little change in the contribution of two bacterioplankton groups whose cultured members harbor known DMSP degradation genes, Roseobacter and SAR11. The DMSP addition led to an enrichment of transcripts supporting heterotrophic activity, and a depletion of those encoding light-related energy generation. Genes for the degradation of C3 compounds were significantly overrepresented after DMSP addition, likely reflecting the metabolism of the C3 component of DMSP. Mapping these transcripts to known biochemical pathways indicated that both acetyl-CoA and succinyl-CoA may be common entry points of this moiety into the tricarboxylic acid cycle. In a short time frame (30 min) in the extremely oligotrophic Sargasso Sea, different gene expression patterns suggest the use of DMSP by a diversity of marine bacterioplankton as both carbon and sulfur sources.

Members of the marine Roseobacter lineage have been characterized as ecological generalists, suggesting that there will be challenges in assigning well-delineated ecological roles and biogeochemical functions to the taxon. To address this issue, genome sequences of 32 Roseobacter isolates were analyzed for patterns in genome characteristics, gene inventory, and individual gene/pathway distribution using three predictive frameworks: phylogenetic relatedness, lifestyle strategy and environmental origin of the isolate. For the first framework, a phylogeny containing five deeply branching clades was obtained from a concatenation of 70 conserved single-copy genes. Somewhat surprisingly, phylogenetic tree topology was not the best model for organizing genome characteristics or distribution patterns of individual genes/pathways, although it provided some predictive power. The lifestyle framework, established by grouping isolates according to evidence for heterotrophy, photoheterotrophy or autotrophy, explained more of the gene repertoire in this lineage. The environment framework had a weak predictive power for the overall genome content of each strain, but explained the distribution of several individual genes/pathways, including those related to phosphorus acquisition, chemotaxis and aromatic compound degradation. Unassembled sequences in the Global Ocean Sampling metagenomic data independently verified this global-scale geographical signal in some Roseobacter genes. The primary findings emerging from this comparative genome analysis are that members of the lineage cannot be easily collapsed into just a few ecologically differentiated clusters (that is, there are almost as many clusters as isolates); the strongest framework for predicting genome content is trophic strategy, but no single framework gives robust predictions; and previously unknown homologs to genes for H(2) oxidation, proteorhodopsin-based phototrophy, xanthorhodpsin-based phototrophy, and CO(2) fixation by Form IC ribulose-1,5-bisphosphate carboxylase/oxygenase (RuBisCO) expand the possible mechanisms for energy and carbon acquisition in this remarkably versatile bacterial lineage.