Multispecies diel transcriptional oscillations in open ocean heterotrophic bacterial assemblages

See allHide authors and affiliations

Science  11 Jul 2014:
Vol. 345, Issue 6193, pp. 207-212
DOI: 10.1126/science.1252476

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 it's 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 ocean's carbon and nitrogen cycles.

Science, this issue p. 207; see also p. 134


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.

The coordination of biological activities into daily periodic cycles is a common feature of eukaryotes and is widespread among plants, fungi, and animals, including man (1). Among single celled noneukaryotic microbes, diel cycles have been well documented in cyanobacterial isolates (24), one halophilic archaeon (5), and bacterial symbionts of fish and squid (6, 7). Some evidence for diel cycling in microbial plankton has also been suggested on the basis of bulk community amino acid incorporation, viral production, or metabolite consumption (810). However, the existence of regular diel oscillations in free-living heterotrophic bacterial species has rarely been assessed.

Microbial community RNA sequencing techniques now allow simultaneous determination of whole-genome transcriptome profiles among multiple cooccurring species (11, 12), enabling high-frequency, time-resolved analyses of microbial community dynamics (12, 13). To better understand temporal transcriptional dynamics in oligotrophic bacterioplankton communities, we conducted a high-resolution multiday time series of bacterioplankton sampled from the North Pacific Subtropical Gyre (14).

To facilitate repeated sampling of the same planktonic microbial populations through time, automated Lagrangian sampling of bacterioplankton was performed every 2 hours over 3 days by using a free-drifting robotic Environmental Sample Processor (ESP) (13, 15) (fig. S1). After instrument recovery, planktonic microbial RNA was extracted, purified, converted to cDNA, and sequenced to assess whole-genome transcriptome dynamics of predominant planktonic microbial populations (tables S1 and S2). The recovered cDNAs were dominated by transcripts from Prochlorococcus and several proteorhodopsin-containing or photoheterotrophic bacteria, including members of the Pelagibacter (SAR11), Roseobacter, SAR116, SAR86, and SAR324 clades (fig. S2).

Phylogenetic analysis of gene transcripts in the most abundant taxa revealed the presence of some microdiversity (figs. S3 to S8). However, the most abundant transcripts sampled at any given time point were dominated by only a few genotypes within each population that persisted throughout the sampling period. An exception was Roseobacter, with transcripts for two different genes (groEL and dnaK) indicating the presence of genotypes that started at a very low abundance and increased in representation over the course of the time series. This variability could be due to an injection of a new population as water masses mixed during the latter portion of the time series or to an alteration in the relative transcriptional activities of two ecotypes that are responding to changes in the surrounding environment.

Transcriptional activity in Prochlorococcus was highly dependent on the time of day. Harmonic regression analyses indicated that nearly half (1491) of all Procholorococcus population transcripts were significantly periodic (Table 1, table S3, and Fig. 1). The expression patterns observed were similar to those of monocultures growing in controlled laboratory settings (4), but there were also notable differences (Fig. 1). For example, photosystem I gene expression exhibited a double peak in the wild Prochlorococcus transcriptome around noon (Fig. 1). In contrast, under laboratory conditions most photosystem I genes, for example, psaL and psaF, were found to peak just before noon, and psaA and psaB peaked shortly after noon (4).

Table 1 Harmonic regression results.

Sequence reads are the total number assigned to each taxon bin. Transcripts are the total number of unique ortholog clusters (see table S1) with at least one mapped sequence. Periodic is the total number of sequences identified as showing 24-hour periodicity using harmonic regression. The next row shows the proportion of variance explained by 24-hour periodicity in constrained principal components analysis (PCA). The fifth row indicates the Procrustes correlation between the first two principal components from Prochlorococcus and other taxa (unconstrained principal components analysis as shown in Fig. 3); P values based on 999 permutations. Bottom row is the correlation between pairwise transcriptome similarities at each time point from Prochlorococcus and heterotrophic taxa based on Mantel test on Euclidean distance matrices; P value based on 999 permutations.

View this table:
Fig. 1 Laboratory versus field comparisons of periodic expression patterns in Prochlorococcus populations.

(A) Scatter plot shows time of peak abundance for 973 transcripts identified as significantly periodic in both studies. Histograms show the total number of genes peaking in 1-hour intervals in this study (top) and the laboratory experiment (side). Black bars represent genes identified as significantly periodic in both studies, gray bars represent genes expressed in both studies but significantly periodic in only one, and white bars represent significantly periodic transcripts that were not detected in the other data set. For this comparison, we used published significance cutoffs from the laboratory study (4) but for consistency generated new peak times by using our harmonic regression approach and the published normalized mean expression levels for each time point. In general, the peak times generated using our approach closely matched published values for that data set. (B and C) Plots showing relative expression (normalized to mean expression level) over time for our metatranscriptome (top trace) and in microarray data (bottom trace) for selected transcripts. For comparison, experimental midnights (24 and 48 hours) from the microarray study are aligned with the 12:00 midnight samples from 9 and 10 September, respectively. All adenosine triphosphate (ATP) synthase subunits (B) and selected subunits from photosystem I (C) are shown.

The largest discrepancy between Prochlorococcus laboratory studies and our field observations was that a considerable number of Prochlorococcus transcripts in our field populations peaked around midday (Fig. 1). Some of these genes did exhibit periodicity in cultures but peaked at a different time of day than in field populations. A larger fraction of the field midday-peaking transcripts were either not periodically expressed or were not present in the culture experiments. In addition, 62% of the 10 a.m.–to–4 p.m. peaking transcripts in our field study lacked Kyoto Encyclopedia of Genes and Genomes (KEGG) orthology annotations, as opposed to those peaking in the evening or late at night.

A number of factors may be responsible for differences in transcript dynamics between in laboratory cultures versus field Prochlorococcus populations. Maximal light levels at our study site at 23-m depth were frequently twofold higher (450 μmol Q m−2 s−1) than those used in laboratory microarray experiments (232 μmol Q m−2 s−1) (4). Fundamental genetic differences between our field populations and the Prochlorococcus strain used in laboratory culture experiments likely also contribute to the differences we observed. Other variables, including nutrient composition and organismal interactions, may also be factors in the observed differences. Although we could not identify obvious trends in the type or function of field population transcripts showing peak expression during the midday period, they did include a wide range of enzymatic functions that are more consistent with nutrient-responsive metabolic changes rather than a simple high-light stress response.

An abundant Roseobacter population also showed strong diel oscillations in its transcriptome profile, most notably in expressed genes involved in bacteriochlorophyll-associated aerobic anoxygenic photosynthesis. Overall, a large fraction of Roseobacter transcripts were periodically expressed (Table 1). Of these, the majority peaked during daylight hours, with only a few gene transcripts peaking at night (Fig. 2). Although this pattern contrasts with that observed in Prochlorococcus, where most diurnally regulated transcripts peaked at dawn or dusk, it was consistent with transcriptional regulation recently reported in Dinoroseobacter shibae (16).

Fig. 2 Timing of periodically expressed transcripts.

For each population, a histogram showing the number of periodically expressed diel transcripts with peak expression within 1-hour intervals throughout the day is shown (left). (Right) Time of peak expression of all transcripts assigned to selected KEGG pathways is plotted (gray circles). Red circles denote transcripts identified as significantly periodic (24-hour period). The transporters category includes both the ATP binding cassette transporters KEGG pathway and the transporters BRITE hierarchy; the photosynthesis category includes both the photosynthesis KEGG pathway and the BRITE photosynthesis proteins categorizations. Carbon fixation refers to genes assigned to the carbon fixation in photosynthetic organisms KEGG pathway. The photosynthesis and carbon fixation categories are present in heterotrophic organisms because of cross-assignment of ATP synthase genes and pentose phosphate cycle genes. Black and yellow bars depict the daily photoperiod (based on sunrise and sunset times). SAGs, single amplified genome-similar transcripts; TCA, tricarboxylic acid.

Thirty-five of the 40 significantly periodic Roseobacter transcripts that peaked between 11 p.m. and 7 a.m. encoded genes belonging to a large photosynthetic “superoperon” (fig. S9). Nightly expression of these genes, followed by immediate repression upon light onset, is consistent with the D. shibae study (16) and may be preparing cells for efficient solar energy harvest in the early morning hours. Functions that peaked during the daytime hours included ribosomal proteins, respiratory transcripts, genes involved in amino acid metabolism, and transporters (Fig. 2).

Proteorhodopsin-containing photoheterotrophs including Pelagibacter (SAR11), SAR116, and SAR86 also showed evidence of diel periodicity in many of their gene transcripts (Fig. 2 and Table 1). All opsin-containing bacteria analyzed (SAR11, SAR116, SAR86, and SAR324) exhibited statistically significant diel oscillations in their proteorhodopsin gene transcripts (table S3 and fig. S10). Peak expression of the opsin transcripts occurred near dawn in all these populations (fig. S10), potentially optimizing solar energy capture by the light-driven, proton-pumping rhodopsins.

Principal components analysis distinguished time series samples for each heterotroph by time of day (Fig. 3) and showed significant correlation with the light-driven behavior of Prochlorococcus (Table 1). Overall, these data are consistent with genome-wide transcriptional changes across the day-night cycle for each population. In addition, coclustering of transcripts using GeneARMA (GA) (14, 17) revealed suites of gene transcripts that exhibited similar expression patterns among different taxa (Fig. 4, figs. S11 to S14, and table S4). For example, a group of transcripts that fit highly similar GA expression models across multiple species included Pro GA5, Pro GA7, Pro GA9, Pro GA23, SAR11 GA6, SAR11 GA18, SAR116 GA2, and Roseobacter GA8 (Fig. 4, fig. S14, and table S4). These multispecies, day-peaking transcripts (fig. S14 and table S4) included gene products associated with respiration (Procholorococcus, SAR11, SAR116, and Roseobacter), nitrogen metabolism (Procholorococcus, SAR11, and SAR116), glycine metabolism (Procholorococcus, SAR11, and Roseobacter), carbon monoxide metabolism (SAR116 and Roseobacter), and DNA synthesis (Procholorococcus and Roseobacter). The coclustering of gene transcripts revealed a complex pattern of expression through the day and across the time series and provides evidence for parallel trends in gene expression across multiple species (Fig. 4, fig. S11 to S14, and table S4).

Fig. 3 Principal component (PC) analyses of population transcriptional profiles.

Transcript abundances were normalized to total transcripts assigned to each population at each time point and arcsin-transformed to approximate normality (19). Symbol color denotes time of day, and shape denotes day of collection. Gray lines connect samples to centroids for selected sample groupings that separate points well. Roseobacter SAG represents samples that were were collected between 7 a.m. and 9 p.m. (versus 9 p.m to 7 a.m.); Pelagibacter, SAR116, and SAR86 clusters represent samples that were collected between 9 a.m. and 6 p.m. (and vice versa); SAR324 cluster represents samples that were collected before or after 9 September 4 p.m. All factor correlations shown were highly significant (P = 0.001). Alternative time of day categories were also highly significant for Roseobacter SAGs, Pelagibacter, SAR116, and SAR86. SAR116 (r2 = 0.10, P = 0.037) and SAR86 (r2 = 0.12, P = 0.019) also correlated weakly with the grouping shown for SAR324. All analyses carried out by using functions in the vegan software package (20).

Fig. 4 Timing of expression of functional gene clusters among different taxa, clustered by the similarities of their temporal expression patterns.

Heat map shows cluster models for all GA clusters, colored by mean-centered relative expression (red, high; blue, low). Black and yellow bars show the daily photoperiod. Each box represents a single sampling event; for sample times, see table S1. Dendrograms show cluster model similarity (Pearson correlations, average linkage clustering, scale bar at upper right represents a correlation of 0.5). The total number of genes (column A), significantly periodic genes (B), and genes associated with photosynthesis (C), ribosome (D), oxidative phosphorylation (E), amino acid metabolism (F), and transport (G) (defined as for Fig. 2), are listed for each cluster. See table S4 for the identities of genes found within any specific GeneARMA (GA) cluster.

Together, the transcriptional profiles of Roseobacter, SAR11, SAR116, and SAR86 indicate diel cycling of metabolic gene transcripts and suggest a multispecies wavelike progression of up-regulated gene suites across the day-night cycle (Fig. 4). Most conspicuously, a regular diel succession of translational, transcriptional, and respiratory gene transcripts was followed by peaks in transporter transcripts that possibly reflect a metabolic recovery phase (Fig. 2). Many of these metabolic pathway transcripts peaked earlier in the day in Roseobacter field populations relative to other bacterial heterotrophs (Figs. 2 and 4 and table S4).

The overall transcriptional profile of SAR324 did not show as many transcript diel oscillations as other heterotrophic taxa. Instead, principal components analysis clustered SAR324 transcripts according to the specific day that they were collected (Fig. 3). In particular, the SAR324 group showed a strong separation between the first portion of the time series and the second in principal components analysis (Fig. 3). This split appears to be associated with the increases in temperature and salinity observed across the time series (fig. S1).

The diurnal patterns reported here for open-ocean heterotrophic bacterioplankton were different from those observed in a previous study of phylogenetically related coastal bacterioplankton using similar methods (12). For example, coastal versus open-ocean SAR11 populations revealed differential expression levels among several orthologous transcript categories (fig. S15). Additionally, whereas the open-ocean SAR11 populations reported here exhibited statistically significant diel oscillations for many gene transcripts (Fig. 2), the coastal SAR11 populations did not.

Currently available data are insufficient to provide definitive mechanistic explanations for the diel behaviors we observed in different heterotrophic bacterioplankton species. It is possible that photoreceptors in these bacteria are involved in regulating light-dark cycles of transcriptional activity. Marine Roseobacter species have previously been shown to regulate their global transcriptional behavior in response to light (16), and laboratory cultures of Pelagibacter also exhibit light-responsive metabolic behaviors (18). However, differences between the behaviors of SAR11 coastal versus open ocean field populations (fig. S15), as well as comparisons of several taxa in our field study versus laboratory experiments on related cultivated isolates (Fig. 1), suggest that other factors may be at play in regulating diel behavior among these different bacterioplankton populations.

Previous studies have proposed that tight metabolic coupling between primary producers and consumers in microbial plankton might elicit conspicuous diel cycling in heterotrophic bacterial activities (8). The diel cycling we observed among different bacterioplankton species is consistent with this hypothesis, with multiple coexisting heterotroph populations exhibiting diurnal oscillations resembling those of their photoautotrophic neighbors. We postulate that the tightly coupled multispecies temporal expression patterns observed may elicit corresponding waves of species-specific metabolic responses at regular time intervals, potentially coordinating diverse biogeochemical activities in these complex microbial communities. Such temporal coordination of biogeochemical activities among multiple species may be important regulators of both the tempo and the mode of microbial matter and energy transformation in the sea.

Supplementary Materials

Materials and Methods

Figs. S1 to S14

Tables S1 to S4

References (2135)

References and Notes

  1. Materials and methods are available as supplementary materials on Science Online.
  2. Acknowledgments: We thank the officers and crew of the Kilo Moana, J. Robidart, S. Wilson, and the ESP engineering and science team (J. Ryan, J. Birch, C. Preston, G. Massion, S. Jensen, and B. Roman) for all of the able assistance. This work was supported by grants from the Gordon and Betty Moore Foundation GBMF nos. 492.01 and 3777 (E.F.D.) and NSF grant EF0424599 (E.F.D.). Development of the ESP was supported by NSF grant OCE-0314222 (to C.A.S.), NASA Astrobiology grants NNG06GB34G and NNX09AB78G (to C.A.S.), the Gordon and Betty Moore Foundation nos. 731 and 2728 (C.A.S.), and the David and Lucile Packard Foundation. This work is a contribution of C-MORE. Sequences reported in this paper have been deposited in the GenBank database (accession no. SRP041215). The Monterey Bay Aquarium Research Institute holds rights to C. A. Scholin et al., U.S. Patent 6187530 (2001). C.A.S. is disqualified from receiving any royalties that might arise from licensing agreements. The ESP is available commercially from Spyglass Technologies and MacLane Research Laboratories; C.A.S. has no financial interest in either company and is not compensated in any way for giving advice on ESP technology transfer.
View Abstract

Stay Connected to Science

Navigate This Article