David M. Needham

Short-term observations of marine bacterial and viral communities: patterns, connections and resilience

Observation of short-term temporal variation in bacterial and viral communities is important for understanding patterns of aquatic microbial diversity. We collected surface seawater once daily for 38 consecutive days with seven more samples interspersed over 40 more days at one location ∼2 km from Santa Catalina Island, California. Bacterial communities were analyzed by automated ribosomal intergenic spacer analysis (ARISA) and viral communities were analyzed by terminal restriction fragment length polymorphism (TRFLP) of the conserved T4-like myoviral gene encoding the major capsid protein (g23). Common bacterial and viral taxa were consistently dominant, and relatively few displayed dramatic increases/decreases or ‘boom/bust’ patterns that might be expected from dynamic predator-prey interactions. Association network analysis showed most significant covariations (associations) occurred among bacterial taxa or among viral taxa and there were several modular (highly-interconnected) associations (P⩽0.005). Associations observed between bacteria and viruses (P⩽0.005) occurred with a median time lag of 2 days. Regression of all pairwise Bray-Curtis similarities between samples indicated a rate of bacterial community change that slows from 2.1%–0.18% per day over a week to 2 months; the rate stays around 0.4% per day for viruses. Our interpretation is that, over the scale of days, individual bacterial and viral OTUs can be dynamic and patterned; resulting in statistical associations regarded as potential ecological interactions. However, over the scale of weeks, average bacterial community variation is slower, suggesting that there is strong community-level ecological resilience, that is, a tendency to converge towards a ‘mean’ microbial community set by longer-term controlling factors.

Temporal variability and coherence of euphotic zone bacterial communities over a decade in the Southern California Bight

Time-series are critical to understanding long-term natural variability in the oceans. Bacterial communities in the euphotic zone were investigated for over a decade at the San Pedro Ocean Time-series station (SPOT) off southern California. Community composition was assessed by Automated Ribosomal Intergenic Spacer Analysis (ARISA) and coupled with measurements of oceanographic parameters for the surface ocean (0–5 m) and deep chlorophyll maximum (DCM, average depth ∼30 m). SAR11 and cyanobacterial ecotypes comprised typically more than one-third of the measured community; diversity within both was temporally variable, although a few operational taxonomic units (OTUs) were consistently more abundant. Persistent OTUs, mostly Alphaproteobacteria (SAR11 clade), Actinobacteria and Flavobacteria, tended to be abundant, in contrast to many rarer yet intermittent and ephemeral OTUs. Association networks revealed potential niches for key OTUs from SAR11, cyanobacteria, SAR86 and other common clades on the basis of robust correlations. Resilience was evident by the average communities drifting only slightly as years passed. Average Bray-Curtis similarity between any pair of dates was ∼40%, with a slight decrease over the decade and obvious near-surface seasonality; communities 8–10 years apart were slightly more different than those 1–4 years apart with the highest rate of change at 0–5 m between communities <4 years apart. The surface exhibited more pronounced seasonality than the DCM. Inter-depth Bray-Curtis similarities repeatedly decreased as the water column stratified each summer. Environmental factors were better predictors of shifts in community composition than months or elapsed time alone; yet, the best predictor was community composition at the other depth (that is, 0–5 m versus DCM).

Pronounced daily succession of phytoplankton, archaea and bacteria following a spring bloom

Marine phytoplankton perform approximately half of global carbon fixation, with their blooms contributing disproportionately to carbon sequestration1, and most phytoplankton production is ultimately consumed by heterotrophic prokaryotes2. Therefore, phytoplankton and heterotrophic community dynamics are important in modelling carbon cycling and the impacts of global change3. In a typical bloom, diatoms dominate initially, transitioning over several weeks to smaller and motile phytoplankton4. Here, we show unexpected, rapid community variation from daily rRNA analysis of phytoplankton and prokaryotic community members following a bloom off southern California. Analysis of phytoplankton chloroplast 16S rRNA demonstrated ten different dominant phytoplankton over 18 days alone, including four taxa with animal toxin-producing strains. The dominant diatoms, flagellates and picophytoplankton varied dramatically in carbon export potential. Dominant prokaryotes also varied rapidly. Euryarchaea briefly became the most abundant organism, peaking over a few days to account for about 40% of prokaryotes. Phytoplankton and prokaryotic communities correlated better with each other than with environmental parameters. Extending beyond the traditional view of blooms being controlled primarily by physics and inorganic nutrients, these dynamics imply highly heterogeneous, continually changing conditions over time and/or space and suggest that interactions among microorganisms are critical in controlling plankton diversity, dynamics and fates.

Seasonal and interannual variability of the marine bacterioplankton community throughout the water column over ten years.

Microbial activities that affect global oceanographic and atmospheric processes happen throughout the water column, yet the long-term ecological dynamics of microbes have been studied largely in the euphotic zone and adjacent seasonally mixed depths. We investigated temporal patterns in the community structure of free-living bacteria, by sampling approximately monthly from 5 m, the deep chlorophyll maximum (∼15–40 m), 150, 500 and 890 m, in San Pedro Channel (maximum depth 900 m, hypoxic below ∼500 m), off the coast of Southern California. Community structure and biodiversity (inverse Simpson index) showed seasonal patterns near the surface and bottom of the water column, but not at intermediate depths. Inverse Simpson’s index was highest in the winter in surface waters and in the spring at 890 m, and varied interannually at all depths. Biodiversity appeared to be driven partially by exchange of microbes between depths and was highest when communities were changing slowly over time. Meanwhile, communities from the surface through 500 m varied interannually. After accounting for seasonality, several environmental parameters co-varied with community structure at the surface and 890 m, but not at the intermediate depths. Abundant and seasonally variable groups included, at 890 m, Nitrospina, Flavobacteria and Marine Group A. Seasonality at 890 m is likely driven by variability in sinking particles, which originate in surface waters, pass transiently through the middle water column and accumulate on the seafloor where they alter the chemical environment. Seasonal subeuphotic groups are likely those whose ecology is strongly influenced by these particles. This surface-to-bottom, decade-long, study identifies seasonality and interannual variability not only of overall community structure, but also of numerous taxonomic groups and near-species level operational taxonomic units.

Ecological dynamics and co-occurrence among marine phytoplankton, bacteria and myoviruses shows microdiversity matters

Numerous ecological processes, such as bacteriophage infection and phytoplankton–bacterial interactions, often occur via strain-specific mechanisms. Therefore, studying the causes of microbial dynamics should benefit from highly resolving taxonomic characterizations. We sampled daily to weekly over 5 months following a phytoplankton bloom off Southern California and examined the extent of microdiversity, that is, significant variation within 99% sequence similarity clusters, operational taxonomic units (OTUs), of bacteria, archaea, phytoplankton chloroplasts (all via 16S or intergenic spacer (ITS) sequences) and T4-like-myoviruses (via g23 major capsid protein gene sequence). The extent of microdiversity varied between genes (ITS most, g23 least) and only temporally common taxa were highly microdiverse. Overall, 60% of taxa exhibited microdiversity; 59% of these had subtypes that changed significantly as a proportion of the parent taxon, indicating ecologically distinct taxa. Pairwise correlations between prokaryotes and myoviruses or phytoplankton (for example, highly microdiverse Chrysochromulina sp.) improved when using single-base variants. Correlations between myoviruses and SAR11 increased in number (172 vs 9, Spearman>0.65) and became stronger (0.61 vs 0.58, t-test: P<0.001) when using SAR11 ITS single-base variants vs OTUs. Whole-community correlation between SAR11 and myoviruses was much improved when using ITS single-base variants vs OTUs, with Mantel rho=0.49 vs 0.27; these results are consistent with strain-specific interactions. Mantel correlations suggested >1 μm (attached/large) prokaryotes are a major myovirus source. Consideration of microdiversity improved observation of apparent host and virus networks, and provided insights into the ecological and evolutionary factors influencing the success of lineages, with important implications to ecosystem resilience and microbial function.

Cross-depth analysis of marine bacterial networks suggests downward propagation of temporal changes

Interactions among microbes and stratification across depths are both believed to be important drivers of microbial communities, though little is known about how microbial associations differ between and across depths. We have monitored the free-living microbial community at the San Pedro Ocean Time-series station, monthly, for a decade, at five different depths: 5 m, the deep chlorophyll maximum layer, 150 m, 500 m and 890 m (just above the sea floor). Here, we introduce microbial association networks that combine data from multiple ocean depths to investigate both within- and between-depth relationships, sometimes time-lagged, among microbes and environmental parameters. The euphotic zone, deep chlorophyll maximum and 890 m depth each contain two negatively correlated ‘modules’ (groups of many inter-correlated bacteria and environmental conditions) suggesting regular transitions between two contrasting environmental states. Two-thirds of pairwise correlations of bacterial taxa between depths lagged such that changes in the abundance of deeper organisms followed changes in shallower organisms. Taken in conjunction with previous observations of seasonality at 890 m, these trends suggest that planktonic microbial communities throughout the water column are linked to environmental conditions and/or microbial communities in overlying waters. Poorly understood groups including Marine Group A, Nitrospina and AEGEAN-169 clades contained taxa that showed diverse association patterns, suggesting these groups contain multiple ecological species, each shaped by different factors, which we have started to delineate. These observations build upon previous work at this location, lending further credence to the hypothesis that sinking particles and vertically migrating animals transport materials that significantly shape the time-varying patterns of microbial community composition.

Genome and epigenome of a novel marine Thaumarchaeota strain suggest viral infection, phosphorothioation DNA modification and multiple restriction systems

Summary Marine Thaumarchaeota are abundant ammonia‐oxidizers but have few representative laboratory‐cultured strains. We report the cultivation of Candidatus Nitrosomarinus catalina SPOT01, a novel strain that is less warm‐temperature tolerant than other cultivated Thaumarchaeota. Using metagenomic recruitment, strain SPOT01 comprises a major portion of Thaumarchaeota (4–54%) in temperate Pacific waters. Its complete 1.36 Mbp genome possesses several distinguishing features: putative phosphorothioation (PT) DNA modification genes; a region containing probable viral genes; and putative urea utilization genes. The PT modification genes and an adjacent putative restriction enzyme (RE) operon likely form a restriction modification (RM) system for defence from foreign DNA. PacBio sequencing showed >98% methylation at two motifs, and inferred PT guanine modification of 19% of possible TGCA sites. Metagenomic recruitment also reveals the putative virus region and PT modification and RE genes are present in 18–26%, 9–14% and <1.5% of natural populations at 150 m with ≥85% identity to strain SPOT01. The presence of multiple probable RM systems in a highly streamlined genome suggests a surprising importance for defence from foreign DNA for dilute populations that infrequently encounter viruses or other cells. This new strain provides new insights into the ecology, including viral interactions, of this important group of marine microbes.

Dynamics and interactions of highly resolved marine plankton via automated high-frequency sampling

Short timescale observations are valuable for understanding microbial ecological processes. We assessed dynamics in relative abundance and potential activities by sequencing the small sub-unit ribosomal RNA gene (rRNA gene) and rRNA molecules (rRNA) of Bacteria, Archaea, and Eukaryota once to twice daily between March 2014 and May 2014 from the surface ocean off Catalina Island, California. Typically Ostreococcus, Braarudosphaera, Teleaulax, and Synechococcus dominated phytoplankton sequences (including chloroplasts) while SAR11, Sulfitobacter, and Fluviicola dominated non-phytoplankton Bacteria and Archaea. We observed short-lived increases of diatoms, mostly Pseudo-nitzschia and Chaetoceros, with quickly responding Bacteria and Archaea including Flavobacteriaceae (Polaribacter & Formosa), Roseovarius, and Euryarchaeota (MGII), notably the exact amplicon sequence variants we observed responding similarly to another diatom bloom nearby, 3 years prior. We observed correlations representing known interactions among abundant phytoplankton rRNA sequences, demonstrating the biogeochemical and ecological relevance of such interactions: (1) The kleptochloroplastidic ciliate Mesodinium 18S rRNA gene sequences and a single Teleaulax taxon (via 16S rRNA gene sequences) were correlated (Spearman r = 0.83) yet uncorrelated to a Teleaulax 18S rRNA gene OTU, or any other taxon (consistent with a kleptochloroplastidic or karyokleptic relationship) and (2) the photosynthetic prymnesiophyte Braarudosphaera bigelowii and two strains of diazotrophic cyanobacterium UCYN-A were correlated and each taxon was also correlated to other taxa, including B. bigelowii to a verrucomicrobium and a dictyochophyte phytoplankter (all r > 0.8). We also report strong correlations (r > 0.7) between various ciliates, bacteria, and phytoplankton, suggesting interactions via currently unknown mechanisms. These data reiterate the utility of high-frequency time series to show rapid microbial reactions to stimuli, and provide new information about in situ dynamics of previously recognized and hypothesized interactions.

Vertical and Seasonal Patterns Control Bacterioplankton Communities at Two Horizontally Coherent Coastal Upwelling Sites off Galicia (NW Spain)

Analysis of seasonal patterns of marine bacterial community structure along horizontal and vertical spatial scales can help to predict long-term responses to climate change. Several recent studies have shown predictable seasonal reoccurrence of bacterial assemblages. However, only a few have assessed temporal variability over both horizontal and vertical spatial scales. Here, we simultaneously studied the bacterial community structure at two different locations and depths in shelf waters of a coastal upwelling system during an annual cycle. The most noticeable biogeographic patterns observed were seasonality, horizontal homogeneity, and spatial synchrony in bacterial diversity and community structure related with regional upwelling–downwelling dynamics. Water column mixing eventually disrupted bacterial community structure vertical heterogeneity. Our results are consistent with previous temporal studies of marine bacterioplankton in other temperate regions and also suggest a marked influence of regional factors on the bacterial communities inhabiting this coastal upwelling system. Bacterial-mediated carbon fluxes in this productive region appear to be mainly controlled by community structure dynamics in surface waters, and local environmental factors at the base of the euphotic zone.

Taxon Disappearance from Microbiome Analysis Reinforces the Value of Mock Communities as a Standard in Every Sequencing Run

Despite the routine use of standards and blanks in virtually all chemical or physical assays and most biological studies (a kind of “control”), microbiome analysis has traditionally lacked such standards. Here we show that unexpected problems of unknown origin can occur in such sequencing runs and yield completely incorrect results that would not necessarily be detected without the use of standards. Assuming that the microbiome sequencing analysis works properly every time risks serious errors that can be detected by the use of mock communities. ABSTRACT Mock communities have been used in microbiome method development to help estimate biases introduced in PCR amplification and sequencing and to optimize pipeline outputs. Nevertheless, the strong value of routine mock community analysis beyond initial method development is rarely, if ever, considered. Here we report that our routine use of mock communities as internal standards allowed us to discover highly aberrant and strong biases in the relative proportions of multiple taxa in a single Illumina HiSeqPE250 run. In this run, an important archaeal taxon virtually disappeared from all samples, and other mock community taxa showed >2-fold high or low abundance, whereas a rerun of those identical amplicons (from the same reaction tubes) on a different date yielded “normal” results. Although obvious from the strange mock community results, we could have easily missed the problem had we not used the mock communities because of natural variation of microbiomes at our site. The “normal” results were validated over four MiSeqPE300 runs and three HiSeqPE250 runs, and run-to-run variation was usually low. While validating these “normal” results, we also discovered that some mock microbial taxa had relatively modest, but consistent, differences between sequencing platforms. We strongly advise the use of mock communities in every sequencing run to distinguish potentially serious aberrations from natural variations. The mock communities should have more than just a few members and ideally at least partly represent the samples being analyzed to detect problems that show up only in some taxa and also to help validate clustering. IMPORTANCE Despite the routine use of standards and blanks in virtually all chemical or physical assays and most biological studies (a kind of “control”), microbiome analysis has traditionally lacked such standards. Here we show that unexpected problems of unknown origin can occur in such sequencing runs and yield completely incorrect results that would not necessarily be detected without the use of standards. Assuming that the microbiome sequencing analysis works properly every time risks serious errors that can be detected by the use of mock communities.