John M. Eppley

Integrated metatranscriptomic and metagenomic analyses of stratified microbial assemblages in the open ocean.

As part of an ongoing survey of microbial community gene expression in the ocean, we sequenced and compared ∼38 Mbp of community transcriptomes and ∼157 Mbp of community genomes from four bacterioplankton samples, along a defined depth profile at Station ALOHA in North Pacific subtropical gyre (NPSG). Taxonomic analysis suggested that the samples were dominated by three taxa: Prochlorales, Consistiales and Cenarchaeales, which comprised 36–69% and 29–63% of the annotated sequences in the four DNA and four cDNA libraries, respectively. The relative abundance of these taxonomic groups was sometimes very different in the DNA and cDNA libraries, suggesting differential relative transcriptional activities per cell. For example, the 125 m sample genomic library was dominated by Pelagibacter (∼36% of sequence reads), which contributed fewer sequences to the community transcriptome (∼11%). Functional characterization of highly expressed genes suggested taxon-specific contributions to specific biogeochemical processes. Examples included Roseobacter relatives involved in aerobic anoxygenic phototrophy at 75 m, and an unexpected contribution of low abundance Crenarchaea to ammonia oxidation at 125 m. Read recruitment using reference microbial genomes indicated depth-specific partitioning of coexisting microbial populations, highlighted by a transcriptionally active high-light-like Prochlorococcus population in the bottom of the photic zone. Additionally, nutrient-uptake genes dominated Pelagibacter transcripts, with apparent enrichment for certain transporter types (for example, the C4-dicarboxylate transport system) over others (for example, phosphate transporters). In total, the data support the utility of coupled DNA and cDNA analyses for describing taxonomic and functional attributes of microbial communities in their natural habitats.

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.

Pattern and synchrony of gene expression among sympatric marine microbial populations

Significance Microbial communities regulate the cycling of energy and matter in the environment, yet how they respond to environmental change is not well-known. We describe here a day in the life of wild planktonic microbial species using robotic sampling coupled with genome-wide gene expression analysis. Our results showed that closely related populations, as well as very different bacterial and archaeal species, displayed remarkably similar time-variable synchronous patterns of gene expression over 2 d. Our results suggest that specific environmental cues may elicit cross-species coordination of gene expression among diverse microbial groups, potentially enabling multispecies coupling of metabolic activity. Planktonic marine microbes live in dynamic habitats that demand rapid sensing and response to periodic as well as stochastic environmental change. The kinetics, regularity, and specificity of microbial responses in situ, however, are not well-described. We report here simultaneous multitaxon genome-wide transcriptome profiling in a naturally occurring picoplankton community. An in situ robotic sampler using a Lagrangian sampling strategy enabled continuous tracking and repeated sampling of coherent microbial populations over 2 d. Subsequent RNA sequencing analyses yielded genome-wide transcriptome profiles of eukaryotic (Ostreococcus) and bacterial (Synechococcus) photosynthetic picoplankton as well as proteorhodopsin-containing heterotrophs, including Pelagibacter, SAR86-cluster Gammaproteobacteria, and marine Euryarchaea. The photosynthetic picoplankton exhibited strong diel rhythms over thousands of gene transcripts that were remarkably consistent with diel cycling observed in laboratory pure cultures. In contrast, the heterotrophs did not cycle diurnally. Instead, heterotrophic picoplankton populations exhibited cross-species synchronous, tightly regulated, temporally variable patterns of gene expression for many genes, particularly those genes associated with growth and nutrient acquisition. This multitaxon, population-wide gene regulation seemed to reflect sporadic, short-term, reversible responses to high-frequency environmental variability. Although the timing of the environmental responses among different heterotrophic species seemed synchronous, the specific metabolic genes that were expressed varied from taxon to taxon. In aggregate, these results provide insights into the kinetics, diversity, and functional patterns of microbial community response to environmental change. Our results also suggest a means by which complex multispecies metabolic processes could be coordinated, facilitating the regulation of matter and energy processing in a dynamically changing environment.

Planktonic Euryarchaeota are a significant source of archaeal tetraether lipids in the ocean

Significance All three domains of life—Eukarya, Bacteria, and Archaea—are known to inhabit the marine water column from surface waters to great depth. Planktonic marine Archaea are comprised of two dominant groups—the Thaumarchaeaota and Euryarchaeota. The Thaumarchaeota contain characteristic biomarker lipids known as tetraethers that are thought to be diagnostic for this group, and are used as paleotemperature proxies, since these lipids are well preserved in marine sediments. In this study, we show that planktonic Euryarchaeota produce the same types of archaeal tetraether lipids as do Thaumarchaeaota. Our results have important implications for environmental surveys of marine Archaea, and the use of their lipids for interpretation of the sedimentary record. Archaea are ubiquitous in marine plankton, and fossil forms of archaeal tetraether membrane lipids in sedimentary rocks document their participation in marine biogeochemical cycles for >100 million years. Ribosomal RNA surveys have identified four major clades of planktonic archaea but, to date, tetraether lipids have been characterized in only one, the Marine Group I Thaumarchaeota. The membrane lipid composition of the other planktonic archaeal groups—all uncultured Euryarchaeota—is currently unknown. Using integrated nucleic acid and lipid analyses, we found that Marine Group II Euryarchaeota (MG-II) contributed significantly to the tetraether lipid pool in the North Pacific Subtropical Gyre at shallow to intermediate depths. Our data strongly suggested that MG-II also synthesize crenarchaeol, a tetraether lipid previously considered to be a unique biomarker for Thaumarchaeota. Metagenomic datasets spanning 5 y indicated that depth stratification of planktonic archaeal groups was a stable feature in the North Pacific Subtropical Gyre. The consistent prevalence of MG-II at depths where the bulk of exported organic matter originates, together with their ubiquitous distribution over diverse oceanic provinces, suggests that this clade is a significant source of tetraether lipids to marine sediments. Our results are relevant to archaeal lipid biomarker applications in the modern oceans and the interpretation of these compounds in the geologic record.

Microbial community phylogenetic and trait diversity declines with depth in a marine oxygen minimum zone

Oxygen minimum zones (OMZs) are natural physical features of the worlds oceans. They create steep physiochemical gradients in the water column, which most notably include a dramatic draw down in oxygen concentrations over small vertical distances (<100 m). Microbial communities within OMZs play central roles in ocean and global biogeochemical cycles, yet we still lack a fundamental understanding of how microbial biodiversity is distributed across OMZs. Here, we used metagenomic sequencing to investigate microbial diversity across a vertical gradient in the water column during three seasons in the Eastern Tropical South Pacific (ETSP) OMZ. Based on analysis of small subunit ribosomal RNA (SSU rRNA) gene fragments, we found that both taxonomic and phylogenetic diversity declined steeply along the transition from oxygen-rich surface water to the permanent OMZ. We observed similar declines in the diversity of protein-coding gene categories, suggesting a decrease in functional (trait) diversity with depth. Metrics of functional and trait dispersion indicated that microbial communities are phylogenetically and functionally more overdispersed in oxic waters, but clustered within the OMZ. These dispersion patterns suggest that community assembly drivers (e.g., competition, environmental filtering) vary strikingly across the oxygen gradient. To understand the generality of our findings, we compared OMZ results to two marine depth gradients in subtropical oligotrophic sites and found that the oligotrophic sites did not display similar patterns, likely reflecting unique features found in the OMZ. Finally, we discuss how our results may relate to niche theory, diversity-energy relationships and stress gradients.

Genetic Exchange Across a Species Boundary in the Archaeal Genus Ferroplasma

Speciation as the result of barriers to genetic exchange is the foundation for the general biological species concept. However, the relevance of genetic exchange for defining microbial species is uncertain. In fact, the extent to which microbial populations comprise discrete clusters of evolutionarily related organisms is generally unclear. Metagenomic data from an acidophilic microbial community enabled a genomewide, comprehensive investigation of variation in individuals from two coexisting natural archaeal populations. Individuals are clustered into species-like groups in which cohesion appears to be maintained by homologous recombination. We quantified the dependence of recombination frequency on sequence similarity genomewide and found a decline in recombination with increasing evolutionary distance. Both inter- and intralineage recombination frequencies have a log-linear dependence on sequence divergence. In the declining phase of interspecies genetic exchange, recombination events cluster near the origin of replication and are localized by tRNAs and short regions of unusually high sequence similarity. The breakdown of genetic exchange with increasing sequence divergence could contribute to, or explain, the establishment and preservation of the observed population clusters in a manner consistent with the biological species concept.

Community transcriptomics reveals universal patterns of protein sequence conservation in natural microbial communities

BackgroundCombined metagenomic and metatranscriptomic datasets make it possible to study the molecular evolution of diverse microbial species recovered from their native habitats. The link between gene expression level and sequence conservation was examined using shotgun pyrosequencing of microbial community DNA and RNA from diverse marine environments, and from forest soil.ResultsAcross all samples, expressed genes with transcripts in the RNA sample were significantly more conserved than non-expressed gene sets relative to best matches in reference databases. This discrepancy, observed for many diverse individual genomes and across entire communities, coincided with a shift in amino acid usage between these gene fractions. Expressed genes trended toward GC-enriched amino acids, consistent with a hypothesis of higher levels of functional constraint in this gene pool. Highly expressed genes were significantly more likely to fall within an orthologous gene set shared between closely related taxa (core genes). However, non-core genes, when expressed above the level of detection, were, on average, significantly more highly expressed than core genes based on transcript abundance normalized to gene abundance. Finally, expressed genes showed broad similarities in function across samples, being relatively enriched in genes of energy metabolism and underrepresented by genes of cell growth.ConclusionsThese patterns support the hypothesis, predicated on studies of model organisms, that gene expression level is a primary correlate of evolutionary rate across diverse microbial taxa from natural environments. Despite their complexity, meta-omic datasets can reveal broad evolutionary patterns across taxonomically, functionally, and environmentally diverse communities.

Microbial community structure and function on sinking particles in the North Pacific Subtropical Gyre.

Sinking particles mediate the transport of carbon and energy to the deep-sea, yet the specific microbes associated with sedimenting particles in the oceans interior remain largely uncharacterized. In this study, we used particle interceptor traps (PITs) to assess the nature of particle-associated microbial communities collected at a variety of depths in the North Pacific Subtropical Gyre. Comparative metagenomics was used to assess differences in microbial taxa and functional gene repertoires in PITs containing a preservative (poisoned traps) compared to preservative-free traps where growth was allowed to continue in situ (live traps). Live trap microbial communities shared taxonomic and functional similarities with bacteria previously reported to be enriched in dissolved organic matter (DOM) microcosms (e.g., Alteromonas and Methylophaga), in addition to other particle and eukaryote-associated bacteria (e.g., Flavobacteriales and Pseudoalteromonas). Poisoned trap microbial assemblages were enriched in Vibrio and Campylobacterales likely associated with eukaryotic surfaces and intestinal tracts as symbionts, pathogens, or saprophytes. The functional gene content of microbial assemblages in poisoned traps included a variety of genes involved in virulence, anaerobic metabolism, attachment to chitinaceaous surfaces, and chitin degradation. The presence of chitinaceaous surfaces was also accompanied by the co-existence of bacteria which encoded the capacity to attach to, transport and metabolize chitin and its derivatives. Distinctly different microbial assemblages predominated in live traps, which were largely represented by copiotrophs and eukaryote-associated bacterial communities. Predominant sediment trap-assocaited eukaryotic phyla included Dinoflagellata, Metazoa (mostly copepods), Protalveolata, Retaria, and Stramenopiles. These data indicate the central role of eukaryotic taxa in structuring sinking particle microbial assemblages, as well as the rapid responses of indigenous microbial species in the degradation of marine particulate organic matter (POM) in situ in the oceans interior.

Time‐series analyses of Monterey Bay coastal microbial picoplankton using a ‘genome proxy’ microarray

To investigate the temporal, spatial and phylogenetic resolution of marine microbial community structure and variability, we designed and expanded a genome proxy array (an oligonucleotide microarray targeting marine microbial genome fragments and genomes), evaluated it against metagenomic sequencing, and applied it to time-series samples from the Monterey Bay. The expanded array targeted 268 microbial genotypes across much of the known diversity of cultured and uncultured marine microbes. The target abundances measured by the array were highly correlated to pyrosequence-based abundances (linear regression R(2) = 0.85-0.91, P < 0.0001). Fifty-seven samples from ∼4 years in Monterey Bay were examined with the array, spanning the photic zone (0 m), the base of the surface mixed layer (30 m) and the subphotic zone (200 m). A significant portion of the expanded genome proxy arrays targets showed signal (95 out of 268 targets present in ≥ 1 sample). The multi-year community survey showed the consistent presence of a core group of common and abundant targeted taxa at each depth in Monterey Bay, higher variability among shallow than deep samples, and episodic occurrences of more transient marine genotypes. The abundance of the most dominant genotypes peaked after strong episodic upwelling events. The genome-proxy arrays ability to track populations of closely related genotypes indicated population shifts within several abundant target taxa, with specific populations in some cases clustering by depth or oceanographic season. Although 51 cultivated organisms were targeted (representing 19% of the array) the majority of targets detected and of total target signal (85% and ∼92% respectively) were from uncultivated genotypes, often those derived from Monterey Bay. The array provided a relatively cost-effective approach (∼

Strainer: software for analysis of population variation in community genomic datasets

15 per array) for surveying the natural history of uncultivated lineages.