Scott M. Gifford

Genome characteristics of a generalist marine bacterial lineage

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 H2 oxidation, proteorhodopsin-based phototrophy, xanthorhodpsin-based phototrophy, and CO2 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.

Quantitative analysis of a deeply sequenced marine microbial metatranscriptome

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 × 1012 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 × 106 transcripts per liter) for 36 genes (for example, phosphonate metabolism gene phnH, dissimilatory nitrate reductase subunit napA) to >2.7 × 109 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.

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.

Expression patterns reveal niche diversification in a marine microbial assemblage.

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.

Multispecies diel transcriptional oscillations in open ocean heterotrophic bacterial assemblages

Up and down go the cyanobacteria Plankton move together in strikingly coordinated daily patterns, sinking at night to avoid being eaten and rising to the surface in daylight to photosynthesize. Otteson et al. found similar activity patterns in even the smallest of planktonic organisms, such as photosynthetic bacteria (see the Perspective by Armbrust). Because its hard to take regular samples in the open ocean, the authors built a robotic sampler and set it adrift for several days in the mid-Pacific. The captured bacteria showed immediate responses to changes in light, temperature, and salinity in ways that could affect the oceans carbon and nitrogen cycles. Science, this issue p. 207; see also p. 134 Multispecies’ daily waves of gene transcription are observed in open ocean microplankton. [Also see Perspective by Armbrust] Oscillating diurnal rhythms of gene transcription, metabolic activity, and behavior are found in all three domains of life. However, diel cycles in naturally occurring heterotrophic bacteria and archaea have rarely been observed. Here, we report time-resolved whole-genome transcriptome profiles of multiple, naturally occurring oceanic bacterial populations sampled in situ over 3 days. As anticipated, the cyanobacterial transcriptome exhibited pronounced diel periodicity. Unexpectedly, several different heterotrophic bacterioplankton groups also displayed diel cycling in many of their gene transcripts. Furthermore, diel oscillations in different heterotrophic bacterial groups suggested population-specific timing of peak transcript expression in a variety of metabolic gene suites. These staggered multispecies waves of diel gene transcription may influence both the tempo and the mode of matter and energy transformation in the sea.

Metatranscriptomic analysis of ammonia-oxidizing organisms in an estuarine bacterioplankton assemblage

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_1354–Nmar_1357, proposed to encode components of an alternative, nitroxyl-based ammonia oxidation pathway; however, reads from Nmar_1259 and Nmar_1667, annotated as encoding a multicopper oxidase with homology to nirK, were abundant. Reads assigned to two homologous ORFs (Nmar_1201 and Nmar_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.

Analyzing Gene Expression from Marine Microbial Communities using Environmental Transcriptomics

Analogous to metagenomics, environmental transcriptomics (metatranscriptomics) retrieves and sequences environmental mRNAs from a microbial assemblage without prior knowledge of what genes the community might be expressing. Thus it provides the most unbiased perspective on community gene expression in situ. Environmental transcriptomics protocols are technically difficult since prokaryotic mRNAs generally lack the poly(A) tails that make isolation of eukaryotic messages relatively straightforward 1 and because of the relatively short half lives of mRNAs 2. In addition, mRNAs are much less abundant than rRNAs in total RNA extracts, thus an rRNA background often overwhelms mRNA signals. However, techniques for overcoming some of these difficulties have recently been developed. A procedure for analyzing environmental transcriptomes by creating clone libraries using random primers to reverse-transcribe and amplify environmental mRNAs was recently described was successful in two different natural environments, but results were biased by selection of the random primers used to initiate cDNA synthesis 3. Advances in linear amplification of mRNA obviate the need for random primers in the amplification step and make it possible to use less starting material decreasing the collection and processing time of samples and thereby minimizing RNA degradation 4. In vitro transcription methods for amplifying mRNA involve polyadenylating the mRNA and incorporating a T7 promoter onto the 3 end of the transcript. Amplified RNA (aRNA) can then be converted to double stranded cDNA using random hexamers and directly sequenced by pyrosequencing 5. A first use of this method at Station ALOHA demonstrated its utility for characterizing microbial community gene expression 6.

Seasonal variation in the metratranscriptomes 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.

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.

High molecular weight dissolved organic matter enrichment selects for methylotrophs in dilution to extinction cultures

The role of bacterioplankton in the cycling of marine dissolved organic matter (DOM) is central to the carbon and energy balance in the ocean, yet there are few model organisms available to investigate the genes, metabolic pathways, and biochemical mechanisms involved in the degradation of this globally important carbon pool. To obtain microbial isolates capable of degrading semi-labile DOM for growth, we conducted dilution to extinction cultivation experiments using seawater enriched with high molecular weight (HMW) DOM. In total, 93 isolates were obtained. Amendments using HMW DOM to increase the dissolved organic carbon concentration 4x (280 μM) or 10x (700 μM) the ocean surface water concentrations yielded positive growth in 4–6% of replicate dilutions, whereas <1% scored positive for growth in non-DOM-amended controls. The majority (71%) of isolates displayed a distinct increase in cell yields when grown in increasing concentrations of HMW DOM. Whole-genome sequencing was used to screen the culture collection for purity and to determine the phylogenetic identity of the isolates. Eleven percent of the isolates belonged to the gammaproteobacteria including Alteromonadales (the SAR92 clade) and Vibrio. Surprisingly, 85% of isolates belonged to the methylotrophic OM43 clade of betaproteobacteria, bacteria thought to metabolically specialize in degrading C1 compounds. Growth of these isolates on methanol confirmed their methylotrophic phenotype. Our results indicate that dilution to extinction cultivation enriched with natural sources of organic substrates has a potential to reveal the previously unsuspected relationships between naturally occurring organic nutrients and the microorganisms that consume them.