Joshua A. Steele

Annually reoccurring bacterial communities are predictable from ocean conditions.

Factors influencing patterns in the distribution and abundance of plant and animal taxa modulate ecosystem function and ecosystem response to environmental change, which is often taken to infer low functional redundancy among such species, but such relationships are poorly known for microbial communities. Using high-resolution molecular fingerprinting, we demonstrate the existence of extraordinarily repeatable temporal patterns in the community composition of 171 operational taxonomic units of marine bacterioplankton over 4.5 years at our Microbial Observatory site, 20 km off the southern California coast. These patterns in distribution and abundance of microbial taxa were highly predictable and significantly influenced by a broad range of both abiotic and biotic factors. These findings provide statistically robust demonstration of temporal patterning in marine bacterial distribution and abundance, which suggests that the distribution and abundance of bacterial taxa may modulate ecosystem function and response and that a significant subset of the bacteria exhibit low levels of functional redundancy as documented for many plant and animal communities.

A latitudinal diversity gradient in planktonic marine bacteria

For two centuries, biologists have documented a gradient of animal and plant biodiversity from the tropics to the poles but have been unable to agree whether it is controlled primarily by productivity, temperature, or historical factors. Recent reports that find latitudinal diversity gradients to be reduced or absent in some unicellular organisms and attribute this to their high abundance and dispersal capabilities would suggest that bacteria, the smallest and most abundant organisms, should exhibit no latitudinal pattern of diversity. We used amplified ribosomal intergenic spacer analysis (ARISA) whole-assemblage genetic fingerprinting to quantify species richness in 103 near-surface samples of marine bacterial plankton, taken from tropical to polar in both hemispheres. We found a significant latitudinal gradient in richness. The data can help to evaluate hypotheses about the cause of the gradient. The correlations of richness with latitude and temperature were similarly strong, whereas correlations with parameters relating to productivity (chlorophyll, annual primary productivity, bacterial abundance) and other variables (salinity and distance to shore) were much weaker. Despite the high abundance and potentially high dispersal of bacteria, they exhibit geographic patterns of species diversity that are similar to those seen in other organisms. The latitudinal gradient in marine bacteria supports the hypothesis that the kinetics of metabolism, setting the pace for life, has strong influence on diversity.

Marine bacterial, archaeal and protistan association networks reveal ecological linkages

Microbes have central roles in ocean food webs and global biogeochemical processes, yet specific ecological relationships among these taxa are largely unknown. This is in part due to the dilute, microscopic nature of the planktonic microbial community, which prevents direct observation of their interactions. Here, we use a holistic (that is, microbial system-wide) approach to investigate time-dependent variations among taxa from all three domains of life in a marine microbial community. We investigated the community composition of bacteria, archaea and protists through cultivation-independent methods, along with total bacterial and viral abundance, and physico-chemical observations. Samples and observations were collected monthly over 3 years at a well-described ocean time-series site of southern California. To find associations among these organisms, we calculated time-dependent rank correlations (that is, local similarity correlations) among relative abundances of bacteria, archaea, protists, total abundance of bacteria and viruses and physico-chemical parameters. We used a network generated from these statistical correlations to visualize and identify time-dependent associations among ecologically important taxa, for example, the SAR11 cluster, stramenopiles, alveolates, cyanobacteria and ammonia-oxidizing archaea. Negative correlations, perhaps suggesting competition or predation, were also common. The analysis revealed a progression of microbial communities through time, and also a group of unknown eukaryotes that were highly correlated with dinoflagellates, indicating possible symbioses or parasitism. Possible ‘keystone’ species were evident. The network has statistical features similar to previously described ecological networks, and in network parlance has non-random, small world properties (that is, highly interconnected nodes). This approach provides new insights into the natural history of microbes.

Virus and prokaryote enumeration from planktonic aquatic environments by epifluorescence microscopy with SYBR Green I

Viruses are the most abundant biological entities in aquatic environments, typically exceeding the abundance of bacteria by an order of magnitude. The reliable enumeration of virus-like particles in marine microbiological investigations is a key measurement parameter. Although the size of typical marine viruses (20–200 nm) is too small to permit the resolution of details by light microscopy, such viruses can be visualized by epifluorescence microscopy if stained brightly. This can be achieved using the sensitive DNA dye SYBR Green I (Molecular Probes–Invitrogen). The method relies on simple vacuum filtration to capture viruses on a 0.02-μm aluminum oxide filter, and subsequent staining and mounting to prepare slides. Virus-like particles are brightly stained and easily observed for enumeration, and prokaryotic cells can easily be counted on the same slides. The protocol provides an inexpensive, rapid (30 min) and reliable technique for obtaining counts of viruses and prokaryotes simultaneously.

Extended local similarity analysis (eLSA) of microbial community and other time series data with replicates

BackgroundThe increasing availability of time series microbial community data from metagenomics and other molecular biological studies has enabled the analysis of large-scale microbial co-occurrence and association networks. Among the many analytical techniques available, the Local Similarity Analysis (LSA) method is unique in that it captures local and potentially time-delayed co-occurrence and association patterns in time series data that cannot otherwise be identified by ordinary correlation analysis. However LSA, as originally developed, does not consider time series data with replicates, which hinders the full exploitation of available information. With replicates, it is possible to understand the variability of local similarity (LS) score and to obtain its confidence interval.ResultsWe extended our LSA technique to time series data with replicates and termed it extended LSA, or eLSA. Simulations showed the capability of eLSA to capture subinterval and time-delayed associations. We implemented the eLSA technique into an easy-to-use analytic software package. The software pipeline integrates data normalization, statistical correlation calculation, statistical significance evaluation, and association network construction steps. We applied the eLSA technique to microbial community and gene expression datasets, where unique time-dependent associations were identified.ConclusionsThe extended LSA analysis technique was demonstrated to reveal statistically significant local and potentially time-delayed association patterns in replicated time series data beyond that of ordinary correlation analysis. These statistically significant associations can provide insights to the real dynamics of biological systems. The newly designed eLSA software efficiently streamlines the analysis and is freely available from the eLSA homepage, which can be accessed at http://meta.usc.edu/softs/lsa.

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).

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.

Co-occurrence patterns for abundant marine archaeal and bacterial lineages in the deep chlorophyll maximum of coastal California

Microorganisms remineralize and respire half of marine primary production, yet the niches occupied by specific microbial groups, and how these different groups may interact, are poorly understood. In this study, we identify co-occurrence patterns for marine Archaea and specific bacterial groups in the chlorophyll maximum of the Southern California Bight. Quantitative PCR time series of marine group 1 (MG1) Crenarchaeota 16S rRNA genes varied substantially over time but were well-correlated (r2=0.94, P<0.001) with ammonia monooxygenase subunit A (amoA) genes, and were more weakly related to 16S rRNA genes for all Archaea (r2=0.39), indicating that other archaeal groups (for example, Euryarchaeota) were numerically important. These data sets were compared with variability in bacterial community composition based on automated ribosomal intergenic spacer analysis (ARISA). We found that archaeal amoA gene copies and a SAR11 (or Pelagibacter) group Ib operational taxonomic unit (OTU) displayed strong co-variation through time (r2=0.55, P<0.05), and archaeal amoA and MG1 16S rRNA genes also co-occurred with two SAR86 and two Bacteroidetes OTUs. The relative abundance of these groups increased and decreased in synchrony over the course of the time series, and peaked during periods of seasonal transition. By using a combination of quantitative and relative abundance estimates, our findings show that abundant microbial OTUs—including the marine Crenarchaeota, SAR11, SAR86 and the Bacteroidetes—co-occur non-randomly; they consequently have important implications for our understanding of microbial community ecology in the sea.

Beta diversity of marine bacteria depends on temporal scale.

Factors controlling the spatial distribution of bacterial diversity have been intensely studied, whereas less is known about temporal changes. To address this, we tested whether the mechanisms that underlie bacterial temporal beta-diversity vary across different scales in three marine microbial communities. While seasonal turnover was detected, at least 73% of the community variation occurred at intra-seasonal temporal scales, suggesting that episodic events are important in structuring marine microbial communities. In addition, turnover at different temporal scales appeared to be driven by different factors. Intra-seasonal turnover was significantly correlated to environmental variables such as phosphate and silicate concentrations, while seasonal and interannual turnover were related to nitrate concentration and temporal distance. We observed a strong link between the magnitude of environmental variation and bacterial beta-diversity in different communities. Analogous to spatial biogeography, we found different rates of community changes across temporal scales.