Annually reoccurring bacterial communities are predictable from ocean conditions

A novel method to accurately calculate statistical significance of local similarity analysis for high-throughput time series

In recent years, a large number of time series microbial community data has been produced in molecular biological studies, especially in metagenomics. Among the statistical methods for time series, local similarity analysis is used in a wide range of environments to capture potential local and time-shifted associations that cannot be distinguished by traditional correlation analysis. Initially, the permutation test is popularly applied to obtain the statistical significance of local similarity analysis. More recently, a theoretical method has also been developed to achieve this aim. However, all these methods require the assumption that the time series are independent and identically distributed. In this paper, we propose a new approach based on moving block bootstrap to approximate the statistical significance of local similarity scores for dependent time series. Simulations show that our method can control the type I error rate reasonably, while theoretical approximation and the permutation test perform less well. Finally, our method is applied to human and marine microbial community datasets, indicating that it can identify potential relationship among operational taxonomic units (OTUs) and significantly decrease the rate of false positives.

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.

Seasonal dynamics of the bacterioplankton community in a large, shallow, highly dynamic freshwater lake

The spatiotemporal shifts of the bacterioplankton community can mirror their transition of functional traits in an aquatic ecosystem. However, the spatiotemporal variation of the bacterioplankton community composition structure (BCCS) within a large, shallow, highly dynamic freshwater lake is still poorly understood. Here, we examined the seasonal and spatial variability of the BCCs within Poyang Lake by sequencing the 16S rRNA gene amplicon to explore how hydrological changes affect the BCCs. Principal coordinate analysis showed that the BCCs varied significantly among four sampling seasons, but not spatially. The seasonal changes of the BCCs were mainly attributed to the differences between autumn and spring–winter. Higher α diversity indices were observed in autumn. Redundancy analysis indicated that the BCCs co-variated with water level, pH, temperature, total phosphorus, ammoniacal nitrogen, electrical conductivity, total nitrogen, and turbidity. Among them, water level was the key determinant separating autumn BCCs from the BCCs in other seasons. A significantly lower relative abundance of Burkholderiales (betI and betVII) and a higher relative abundance of Actinomycetales (acI, acTH1, and acTH2) were found in autumn than in other seasons. Overall, our results suggest that water level changes associated with pH, temperature, and nutrient status shaped the seasonal patterns of the BCCs within Poyang Lake.

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.

A strong link between marine microbial community composition and function challenges the idea of functional redundancy

Marine microbes have tremendous diversity, but a fundamental question remains unanswered: why are there so many microbial species in the sea? The idea of functional redundancy for microbial communities has long been assumed, so that the high level of richness is often explained by the presence of different taxa that are able to conduct the exact same set of metabolic processes and that can readily replace each other. Here, we refute the hypothesis of functional redundancy for marine microbial communities by showing that a shift in the community composition altered the overall functional attributes of communities across different temporal and spatial scales. Our metagenomic monitoring of a coastal northwestern Mediterranean site also revealed that diverse microbial communities harbor a high diversity of potential proteins. Working with all information given by the metagenomes (all reads) rather than relying only on known genes (annotated orthologous genes) was essential for revealing the similarity between taxonomic and functional community compositions. Our finding does not exclude the possibility for a partial redundancy where organisms that share some specific function can coexist when they differ in other ecological requirements. It demonstrates, however, that marine microbial diversity reflects a tremendous diversity of microbial metabolism and highlights the genetic potential yet to be discovered in an ocean of microbes.

Optofluidic Single-Cell Genome Amplification of Sub-micron Bacteria in the Ocean Subsurface

Optofluidic single-cell genome amplification was used to obtain genome sequences from sub-micron cells collected from the euphotic and mesopelagic zones of the northwestern Sargasso Sea. Plankton cells were visually selected and manually sorted with an optical trap, yielding 20 partial genome sequences representing seven bacterial phyla. Two organisms, E01-9C-26 (Gammaproteobacteria), represented by four single cell genomes, and Opi.OSU.00C, an uncharacterized Verrucomicrobia, were the first of their types retrieved by single cell genome sequencing and were studied in detail. Metagenomic data showed that E01-9C-26 is found throughout the dark ocean, while Opi.OSU.00C was observed to bloom transiently in the nutrient-depleted euphotic zone of the late spring and early summer. The E01-9C-26 genomes had an estimated size of 4.76-5.05 Mbps, and contained "O" and "W"-type monooxygenase genes related to methane and ammonium monooxygenases that were previously reported from ocean metagenomes. Metabolic reconstruction indicated E01-9C-26 are likely versatile methylotrophs capable of scavenging C1 compounds, methylated compounds, reduced sulfur compounds, and a wide range of amines, including D-amino acids. The genome sequences identified E01-9C-26 as a source of "O" and "W"-type monooxygenase genes related to methane and ammonium monooxygenases that were previously reported from ocean metagenomes, but are of unknown function. In contrast, Opi.OSU.00C genomes encode genes for catabolizing carbohydrate compounds normally associated with eukaryotic phytoplankton. This exploration of optofluidics showed that it was effective for retrieving diverse single-cell bacterioplankton genomes and has potential advantages in microbiology applications that require working with small sample volumes or targeting cells by their morphology.

Prokaryotes in the WAIS Divide ice core reflect source and transport changes between Last Glacial Maximum and the early Holocene

We present the first long-term, highly resolved prokaryotic cell concentration record obtained from a polar ice core. This record, obtained from the West Antarctic Ice Sheet (WAIS) Divide (WD) ice core, spanned from the Last Glacial Maximum (LGM) to the early Holocene (EH) and showed distinct fluctuations in prokaryotic cell concentration coincident with major climatic states. The time series also revealed a ~1,500-year periodicity with greater amplitude during the Last Deglaciation (LDG). Higher prokaryotic cell concentration and lower variability occurred during the LGM and EH than during the LDG. A sevenfold decrease in prokaryotic cell concentration coincided with the LGM/LDG transition and the global 19 ka meltwater pulse. Statistical models revealed significant relationships between the prokaryotic cell record and tracers of both marine (sea-salt sodium [ssNa]) and burning emissions (black carbon [BC]). Collectively, these models, together with visual observations and methanosulfidic acid (MSA) measurements, indicated that the temporal variability in concentration of airborne prokaryotic cells reflected changes in marine/sea-ice regional environments of the WAIS. Our data revealed that variations in source and transport were the most likely processes producing the significant temporal variations in WD prokaryotic cell concentrations. This record provided strong evidence that airborne prokaryotic cell deposition differed during the LGM, LDG, and EH, and that these changes in cell densities could be explained by different environmental conditions during each of these climatic periods. Our observations provide the first ice-core time series evidence for a prokaryotic response to long-term climatic and environmental processes.

Verrucomicrobia are prevalent in north-temperate freshwater lakes and display class-level preferences between lake habitats

The bacterial phylum Verrucomicrobia was formally described two decades ago and originally believed to be a minor member of many ecosystems; however, it is now recognized as ubiquitous and abundant in both soil and aquatic systems. Nevertheless, knowledge of the drivers of its relative abundance and within-phylum habitat preferences remains sparse, especially in lake systems. Here, we documented the distribution of Verrucomicrobia in 12 inland lakes in Southeastern Michigan, a Laurentian Great Lake (Lake Michigan), and a freshwater estuary, which span a gradient in lake sizes, depths, residence times, and trophic states. A wide range of physical and geochemical parameters was covered by sampling seasonally from the surface and bottom of each lake, and by separating samples into particle-associated and free-living fractions. On average, Verrucomicrobia was the 4^th most abundant phylum (range 1.7–41.7%). Fraction, season, station, and depth explained up to 70% of the variance in Verrucomicrobia community composition and preference for these habitats was phylogenetically conserved at the class-level. When relative abundance was linearly modeled against environmental data, Verrucomicrobia and non-Verrucomicrobia bacterial community composition correlated to similar quantitative environmental parameters, although there were lake system-dependent differences and > 55% of the variance remained unexplained. A majority of the phylum exhibited preference for the particle-associated fraction and two classes (Opitutae and Verrucomicrobiae) were identified to be more abundant during the spring season. This study highlights the high relative abundance of Verrucomicrobia in north temperate lake systems and expands insights into drivers of within-phylum habitat preferences of the Verrucomicrobia.

Mixing alters the lytic activity of viruses in the dark ocean

In aquatic habitats, viral lysis of prokaryotic cells lowers the overall efficiency of the microbial loop, by which dissolved organic carbon is transfered to higher trophic levels. Mixing of water masses in the dark ocean occurs on a global scale and may have far reaching consequences for the different prokaryotic and virus communities found in these waters by altering the environmental conditions these communities experience. We hypothesize that mixing of deep ocean water masses enhances the lytic activity of viruses infecting prokaryotes. To address this hypothesis, major deep-sea water masses of the Atlantic Ocean such as North Atlantic Deep Water, Mediterranean Sea Overflow Water, Antarctic Intermediate Water, and Antarctic Bottom Water were sampled at five locations. Prokaryotic cells from these samples were collected by filtration and subsequently incubated in virus-reduced water from either the same (control) or a different water mass (transplantation treatment). Additionally, mixtures of prokaryotes obtained from two different water masses were incubated in a mixture of virus-reduced water from the same water masses (control) or in virus-reduced water from the source water masses separately (mixing treatments). Pronounced differences in productivity-related parameters (prokaryotic leucine incorporation, prokaryotic and viral abundance) between water masses caused strong changes in viral lysis of prokaryotes. Often, mixing of water masses increased viral lysis of prokaryotes, indicating that lysogenic viruses were induced into the lytic cycle. Mixing-induced changes in viral lysis had a strong effect on the community composition of prokaryotes and viruses. Our data show that mixing of deep-sea water masses alters levels of viral lysis of prokaryotes and in many cases weakens the efficiency of the microbial loop by enhancing the recycling of organic carbon in the deep ocean.

High resolution time series reveals cohesive but short-lived communities in coastal plankton

Because microbial plankton in the ocean comprise diverse bacteria, algae, and protists that are subject to environmental forcing on multiple spatial and temporal scales, a fundamental open question is to what extent these organisms form ecologically cohesive communities. Here we show that although all taxa undergo large, near daily fluctuations in abundance, microbial plankton are organized into clearly defined communities whose turnover is rapid and sharp. We analyze a time series of 93 consecutive days of coastal plankton using a technique that allows inference of communities as modular units of interacting taxa by determining positive and negative correlations at different temporal frequencies. This approach shows both coordinated population expansions that demarcate community boundaries and high frequency of positive and negative associations among populations within communities. Our analysis thus highlights that the environmental variability of the coastal ocean is mirrored in sharp transitions of defined but ephemeral communities of organisms.

Seasonal Shifts in Bacterial Community Responses to Phytoplankton-Derived Dissolved Organic Matter in the Western Antarctic Peninsula

Bacterial consumption of dissolved organic matter (DOM) drives much of the movement of carbon through the oceanic food web and the global carbon cycle. Understanding complex interactions between bacteria and marine DOM remains an important challenge. We tested the hypothesis that bacterial growth and community succession would respond differently to DOM additions due to seasonal changes in phytoplankton abundance in the environment. Four mesocosm experiments were conducted that spanned the spring transitional period (August-December 2013) in surface waters of the Western Antarctic Peninsula (WAP). Each mesocosm consisted of nearshore surface seawater (50 L) incubated in the laboratory for 10 days. The addition of DOM, in the form of cell-free exudates extracted from Thalassiosira weissflogii diatom cultures led to changes in bacterial abundance, production, and community composition. The timing of each mesocosm experiment (i.e., late winter vs. late spring) influenced the magnitude and direction of bacterial changes. For example, the same DOM treatment applied at different times during the season resulted in different levels of bacterial production and different bacterial community composition. There was a mid-season shift from Collwelliaceae to Polaribacter having the greatest relative abundance after incubation. This shift corresponded to a modest but significant increase in the initial relative abundance of Polaribacter in the nearshore seawater used to set up experiments. This finding supports a new hypothesis that starting community composition, through priority effects, influenced the trajectory of community succession in response to DOM addition. As strong inter-annual variability and long-term climate change may shift the timing of WAP phytoplankton blooms, and the corresponding production of DOM exudates, this study suggests a mechanism by which different seasonal successional patterns in bacterial communities could occur.

Recurrent patterns of microdiversity in a temperate coastal marine environment

Temperate coastal marine environments are replete with complex biotic and abiotic interactions that are amplified during spring and summer phytoplankton blooms. During these events, heterotrophic bacterioplankton respond to successional releases of dissolved organic matter as algal cells are lysed. Annual seasonal shifts in the community composition of free-living bacterioplankton follow broadly predictable patterns, but whether similar communities respond each year to bloom disturbance events remains unknown owing to a lack of data sets, employing high-frequency sampling over multiple years. We capture the fine-scale microdiversity of these events with weekly sampling using a high-resolution method to discriminate 16S ribosomal RNA gene amplicons that are >99% identical. Furthermore, we used 2 complete years of data to facilitate identification of recurrent sub-networks of co-varying microbes. We demonstrate that despite inter-annual variation in phytoplankton blooms and despite the dynamism of a coastal–oceanic transition zone, patterns of microdiversity are recurrent during both bloom and non-bloom conditions. Sub-networks of co-occurring microbes identified reveal that correlation structures between community members appear quite stable in a seasonally driven response to oligotrophic and eutrophic conditions.

Analysing Microbial Community Composition through Amplicon Sequencing: From Sampling to Hypothesis Testing

Microbial ecology as a scientific field is fundamentally driven by technological advance. The past decade's revolution in DNA sequencing cost and throughput has made it possible for most research groups to map microbial community composition in environments of interest. However, the computational and statistical methodology required to analyse this kind of data is often not part of the biologist training. In this review, we give a historical perspective on the use of sequencing data in microbial ecology and restate the current need for this method; but also highlight the major caveats with standard practices for handling these data, from sample collection and library preparation to statistical analysis. Further, we outline the main new analytical tools that have been developed in the past few years to bypass these caveats, as well as highlight the major requirements of common statistical practices and the extent to which they are applicable to microbial data. Besides delving into the meaning of select alpha- and beta-diversity measures, we give special consideration to techniques for finding the main drivers of community dissimilarity and for interaction network construction. While every project design has specific needs, this review should serve as a starting point for considering what options are available.

Phylogenetically conserved resource partitioning in the coastal microbial loop

Resource availability influences marine microbial community structure, suggesting that population-specific resource partitioning defines discrete niches. Identifying how resources are partitioned among populations, thereby characterizing functional guilds within the communities, remains a challenge for microbial ecologists. We used proteomic stable isotope probing (SIP) and NanoSIMS analysis of phylogenetic microarrays (Chip-SIP) along with 16S rRNA gene amplicon and metagenomic sequencing to characterize the assimilation of six ¹³C-labeled common metabolic substrates and changes in the microbial community structure within surface water collected from Monterey Bay, CA. Both sequencing approaches indicated distinct substrate-specific community shifts. However, observed changes in relative abundance for individual populations did not correlate well with directly measured substrate assimilation. The complementary SIP techniques identified assimilation of all six substrates by diverse taxa, but also revealed differential assimilation of substrates into protein and ribonucleotide biomass between taxa. Substrate assimilation trends indicated significantly conserved resource partitioning among populations within the Flavobacteriia, Alphaproteobacteria and Gammaproteobacteria classes, suggesting that functional guilds within marine microbial communities are phylogenetically cohesive. However, populations within these classes exhibited heterogeneity in biosynthetic activity, which distinguished high-activity copiotrophs from low-activity oligotrophs. These results indicate distinct growth responses between populations that is not apparent by genome sequencing alone.

Characterising and predicting cyanobacterial blooms in an 8-year amplicon sequencing time course

Cyanobacterial blooms occur in lakes worldwide, producing toxins that pose a serious public health threat. Eutrophication caused by human activities and warmer temperatures both contribute to blooms, but it is still difficult to predict precisely when and where blooms will occur. One reason that prediction is so difficult is that blooms can be caused by different species or genera of cyanobacteria, which may interact with other bacteria and respond to a variety of environmental cues. Here we used a deep 16S amplicon sequencing approach to profile the bacterial community in eutrophic Lake Champlain over time, to characterise the composition and repeatability of cyanobacterial blooms, and to determine the potential for blooms to be predicted based on time course sequence data. Our analysis, based on 135 samples between 2006 and 2013, spans multiple bloom events. We found that bloom events significantly alter the bacterial community without reducing overall diversity, suggesting that a distinct microbial community-including non-cyanobacteria-prospers during the bloom. We also observed that the community changes cyclically over the course of a year, with a repeatable pattern from year to year. This suggests that, in principle, bloom events are predictable. We used probabilistic assemblages of OTUs to characterise the bloom-associated community, and to classify samples into bloom or non-bloom categories, achieving up to 92% classification accuracy (86% after excluding cyanobacterial sequences). Finally, using symbolic regression, we were able to predict the start date of a bloom with 78-92% accuracy (depending on the data used for model training), and found that sequence data was a better predictor than environmental variables.

Marine archaeal dynamics and interactions with the microbial community over 5 years from surface to seafloor

Marine archaea are critical contributors to global carbon and nitrogen redox cycles, but their temporal variability and microbial associations across the water column are poorly known. We evaluated seasonal variability of free living (0.2-1 μm size fraction) Thaumarchaea Marine Group I (MGI) and Euryarchaea Marine Group II (MGII) communities and their associations with the microbial community from surface to seafloor (890 m) over 5 years by 16S rRNA V4-V5 gene sequencing. MGI and MGII communities demonstrated distinct compositions at different depths, and seasonality at all depths. Microbial association networks at 150 m, 500 m and 890 m, revealed diverse assemblages of MGI (presumed ammonia oxidizers) and Nitrospina taxa (presumed dominant nitrite oxidizers, completing the nitrification process), suggesting distinct MGI-Nitrospina OTUs are responsible for nitrification at different depths and seasons, and depth- related and seasonal variability in nitrification could be affected by alternating MGI-Nitrospina assemblages. MGII taxa also showed distinct correlations to possibly heterotrophic bacteria, most commonly to members of Marine Group A, Chloroflexi, Marine Group B, and SAR86. Thus, both MGI and MGII likely have dynamic associations with bacteria based on similarities in activity or other interactions that select for distinct microbial assemblages over time. The importance of MGII taxa as members of the heterotrophic community previously reported for photic zone appears to apply throughout the water column.

Bacterial Community Composition and Dynamics Spanning Five Years in Freshwater Bog Lakes

Lakes are excellent systems for investigating bacterial community dynamics because they have clear boundaries and strong environmental gradients. The results of our research demonstrate that bacterial community composition varies by year, a finding which likely applies to other ecosystems and has implications for study design and interpretation. Understanding the drivers and controls of bacterial communities on long time scales would improve both our knowledge of fundamental properties of bacterial communities and our ability to predict community states. In this specific ecosystem, bog lakes play a disproportionately large role in global carbon cycling, and the information presented here may ultimately help refine carbon budgets for these lakes. Finally, all data and code in this study are publicly available. We hope that this will serve as a resource for anyone seeking to answer their own microbial ecology questions using a multiyear time series.

Bacteria play a key role in freshwater biogeochemical cycling, but long-term trends in freshwater bacterial community composition and dynamics are not yet well characterized. We used a multiyear time series of 16S rRNA gene amplicon sequencing data from eight bog lakes to census the freshwater bacterial community and observe annual and seasonal trends in abundance. The sites that we studied encompassed a range of water column mixing frequencies, which we hypothesized would be associated with trends in alpha and beta diversity. Each lake and layer contained a distinct bacterial community, with distinct levels of richness and indicator taxa that likely reflected the environmental conditions of each lake type sampled, including Actinobacteria in polymictic lakes (i.e., lakes with multiple mixing events per year), Methylophilales in dimictic lakes (lakes with two mixing events per year, usually in spring and fall), and “Candidatus Omnitrophica” in meromictic lakes (lakes with no recorded mixing events). The community present during each year at each site was also surprisingly unique. Despite unexpected interannual variability in community composition, we detected a core community of taxa found in all lakes and layers, including Actinobacteria tribe acI-B2 and Betaprotobacteria lineage PnecC. Although trends in abundance did not repeat annually, each freshwater lineage within the communities had a consistent lifestyle, defined by persistence, abundance, and variability. The results of our analysis emphasize the importance of long-term multisite observations, as analyzing only a single year of data or one lake would not have allowed us to describe the dynamics and composition of these freshwater bacterial communities to the extent presented here.

IMPORTANCE Lakes are excellent systems for investigating bacterial community dynamics because they have clear boundaries and strong environmental gradients. The results of our research demonstrate that bacterial community composition varies by year, a finding which likely applies to other ecosystems and has implications for study design and interpretation. Understanding the drivers and controls of bacterial communities on long time scales would improve both our knowledge of fundamental properties of bacterial communities and our ability to predict community states. In this specific ecosystem, bog lakes play a disproportionately large role in global carbon cycling, and the information presented here may ultimately help refine carbon budgets for these lakes. Finally, all data and code in this study are publicly available. We hope that this will serve as a resource for anyone seeking to answer their own microbial ecology questions using a multiyear time series.

Functional Gene Diversity and Metabolic Potential of the Microbial Community in an Estuary-Shelf Environment

Microbes play crucial roles in various biogeochemical processes in the ocean, including carbon (C), nitrogen (N), and phosphorus (P) cycling. Functional gene diversity and the structure of the microbial community determines its metabolic potential and therefore its ecological function in the marine ecosystem. However, little is known about the functional gene composition and metabolic potential of bacterioplankton in estuary areas. The East China Sea (ECS) is a dynamic marginal ecosystem in the western Pacific Ocean that is mainly affected by input from the Changjiang River and the Kuroshio Current. Here, using a high-throughput functional gene microarray (GeoChip), we analyzed the functional gene diversity, composition, structure, and metabolic potential of microbial assemblages in different ECS water masses. Four water masses determined by temperature and salinity relationship showed different patterns of functional gene diversity and composition. Generally, functional gene diversity [Shannon–Weaner’s H and reciprocal of Simpson’s 1/(1-D)] in the surface water masses was higher than that in the bottom water masses. The different presence and proportion of functional genes involved in C, N, and P cycling among the bacteria of the different water masses showed different metabolic preferences of the microbial populations in the ECS. Genes involved in starch metabolism (amyA and nplT) showed higher proportion in microbial communities of the surface water masses than of the bottom water masses. In contrast, a higher proportion of genes involved in chitin degradation was observed in microorganisms of the bottom water masses. Moreover, we found a higher proportion of nitrogen fixation (nifH), transformation of hydroxylamine to nitrite (hao) and ammonification (gdh) genes in the microbial communities of the bottom water masses compared with those of the surface water masses. The spatial variation of microbial functional genes was significantly correlated with salinity, temperature, and chlorophyll based on canonical correspondence analysis, suggesting a significant influence of hydrologic conditions on water microbial communities. Our data provide new insights into better understanding of the functional potential of microbial communities in the complex estuarine-coastal environmental gradient of the ECS.

Changes in Microbial (Bacteria and Archaea) Plankton Community Structure after Artificial Dispersal in Grazer-Free Microcosms

Microbes are considered to have a global distribution due to their high dispersal capabilities. However, our knowledge of the way geographically distant microbial communities assemble after dispersal in a new environment is limited. In this study, we examined whether communities would converge because similar taxa would be selected under the same environmental conditions, or would diverge because of initial community composition, after artificial dispersal. To this aim, a microcosm experiment was performed, in which the temporal changes in the composition and diversity of different prokaryoplankton assemblages from three distant geographic coastal areas (Banyuls-sur-Mer in northwest Mediterranean Sea, Pagasitikos Gulf in northeast Mediterranean and Woods Hole, MA, USA in the northwest Atlantic), were studied. Diversity was investigated using amplicon pyrosequencing of the V1–V3 hypervariable regions of the 16S rRNA. The three assemblages were grown separately in particle free and autoclaved Banyuls-sur-mer seawater at 18 °C in the dark. We found that the variability of prokaryoplankton community diversity (expressed as richness, evenness and dominance) as well as the composition were driven by patterns observed in Bacteria. Regarding community composition, similarities were found between treatments at family level. However, at the OTU level microbial communities from the three different original locations diverge rather than converge during incubation. It is suggested that slight differences in the composition of the initial prokaryoplankton communities, resulted in separate clusters the following days even when growth took place under identical abiotic conditions.

Bacterial Biogeography across the Amazon River-Ocean Continuum

Spatial and temporal patterns in microbial biodiversity across the Amazon river-ocean continuum were investigated along ∼675 km of the lower Amazon River mainstem, in the Tapajós River tributary, and in the plume and coastal ocean during low and high river discharge using amplicon sequencing of 16S rRNA genes in whole water and size-fractionated samples (0.2–2.0 μm and >2.0 μm). River communities varied among tributaries, but mainstem communities were spatially homogeneous and tracked seasonal changes in river discharge and co-varying factors. Co-occurrence network analysis identified strongly interconnected river assemblages during high (May) and low (December) discharge periods, and weakly interconnected transitional assemblages in September, suggesting that this system supports two seasonal microbial communities linked to river discharge. In contrast, plume communities showed little seasonal differences and instead varied spatially tracking salinity. However, salinity explained only a small fraction of community variability, and plume communities in blooms of diatom-diazotroph assemblages were strikingly different than those in other high salinity plume samples. This suggests that while salinity physically structures plumes through buoyancy and mixing, the composition of plume-specific communities is controlled by other factors including nutrients, phytoplankton community composition, and dissolved organic matter chemistry. Co-occurrence networks identified interconnected assemblages associated with the highly productive low salinity near-shore region, diatom-diazotroph blooms, and the plume edge region, and weakly interconnected assemblages in high salinity regions. This suggests that the plume supports a transitional community influenced by immigration of ocean bacteria from the plume edge, and by species sorting as these communities adapt to local environmental conditions. Few studies have explored patterns of microbial diversity in tropical rivers and coastal oceans. Comparison of Amazon continuum microbial communities to those from temperate and arctic systems suggest that river discharge and salinity are master variables structuring a range of environmental conditions that control bacterial communities across the river-ocean continuum.

Nearly a decade-long repeatable seasonal diversity patterns of bacterioplankton communities in the eutrophic Lake Donghu (Wuhan, China)

Uncovering which environmental factors govern community diversity patterns and how ecological processes drive community turnover are key questions related to understand the community assembly. However, the ecological mechanisms regulating long-term variations of bacterioplankton communities in lake ecosystems remain poorly understood. Here we present nearly a decade-long study of bacterioplankton communities from the eutrophic Lake Donghu (Wuhan, China) using 16S rRNA gene amplicon sequencing with MiSeq platform. We found strong repeatable seasonal diversity patterns in terms of both common (detected in more than 50% samples) and dominant (relative abundance >1%) bacterial taxa turnover. Moreover, community composition tracked the seasonal temperature gradient, indicating that temperature is a key environmental factor controlling observed diversity patterns. Total phosphorus also contributed significantly to the seasonal shifts in bacterioplankton composition. However, any spatial pattern of bacterioplankton communities across the main lake areas within season was overwhelmed by their temporal variabilities. Phylogenetic analysis further indicated that 75%-82% of community turnover was governed by homogeneous selection due to consistent environmental conditions within seasons, suggesting that the microbial communities in Lake Donghu are mainly controlled by niche-based processes. Therefore, dominant niches available within seasons might be occupied by similar combinations of bacterial taxa with modest dispersal rates throughout different lake areas.

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.

Impact of hypoxia on gene expression patterns by the human pathogen, Vibrio vulnificus, and bacterial community composition in a North Carolina estuary

Estuarine environments are continuously being shaped by both natural and anthropogenic sources which directly/indirectly influence the organisms that inhabit these important niches on both individual and community levels. Human infections caused by pathogenic Vibrio species are continuing to rise, and factors associated with global climate change have been suggested to be impacting their abundance and geographical range. Along with temperature, hypoxia has also increased dramatically in the last 40 years, which has led to persistent dead zones worldwide in areas where these infections are increasing. Thus, utilizing membrane diffusion chambers, we investigated the impact of in situ hypoxia on the gene expression of one such bacterium, Vibrio vulnificus, which is an inhabitant of these vulnerable areas worldwide. By coupling these data with multiple abiotic factors, we were able to demonstrate that genes involved in numerous functions, including those involved in virulence, environmental persistence, and stressosome production, were negatively correlated with dissolved oxygen. Furthermore, comparing 16S ribosomal RNA, we found similar overall community compositions during both hypoxia and normoxia. However, unweighted beta diversity analyses revealed that although certain classes of bacteria dominate in both low- and high-oxygen environments, there is the potential for quantitative shifts in lower abundant species, which may be important for effective risk assessment in areas that are becoming increasingly more hypoxic. This study emphasizes the importance of investigating hypoxia as a trigger for gene expression changes by marine Vibrio species and highlights the need for more in depth community analyses during estuarine hypoxia.

Microbe biogeography tracks water masses in a dynamic oceanic frontal system

Dispersal limitation, not just environmental selection, plays an important role in microbial biogeography. The distance-decay relationship is thought to be weak in habitats where dispersal is high, such as in the pelagic environment, where ocean currents facilitate microbial dispersal. Most studies of microbial community composition to date have observed little geographical heterogeneity on a regional scale (100 km). We present a study of microbial communities across a dynamic frontal zone in the southwest Indian Ocean and investigate the spatial structure of the microbes with respect to the different water masses separated by these fronts. We collected 153 samples of free-living microorganisms from five seamounts located along a gradient from subtropical to subantarctic waters and across three depth layers: (i) the sub-surface chlorophyll maximum (approx. 40 m), (ii) the bottom of the euphotic zone (approx. 200 m), and (iii) the benthic boundary layer (300-2000 m). Diversity and abundance of microbial operational taxonomic units (OTUs) were assessed by amplification and sequencing of the 16S rRNA gene on an Illumina MiSeq platform. Multivariate analyses showed that microbial communities were structured more strongly by depth than by latitude, with similar phyla occurring within each depth stratum across seamounts. The deep layer was homogeneous across the entire survey area, corresponding to the spread of Antarctic intermediate water. However, within both the sub-surface layer and the intermediate depth stratum there was evidence for OTU turnover across fronts. The microbiome of these layers appears to be divided into three distinct biological regimes corresponding to the subantarctic surface water, the convergence zone and subtropical. We show that microbial biogeography across depth and latitudinal gradients is linked to the water masses the microbes persist in, resulting in regional patterns of microbial biogeography that correspond to the regional scale physical oceanography.

Annual community patterns are driven by seasonal switching between closely related marine bacteria

Marine microbes exhibit seasonal cycles in community composition, yet the key drivers of these patterns and microbial population fidelity to specific environmental conditions remain to be determined. To begin addressing these questions, we characterized microbial dynamics weekly for 3 years at a temperate, coastal site with dramatic environmental seasonality. This high-resolution time series reveals that changes in microbial community composition are not continuous; over the duration of the time series, the community instead resolves into distinct summer and winter profiles with rapid spring and fall transitions between these states. Here, we show that these community shifts involve switching between closely related strains that exhibit either summer or winter preferences. Moreover, taxa repeat this process annually in both this and another temperate coastal time series, suggesting that this phenomenon may be widespread in marine ecosystems. To address potential biogeochemical impacts of these community changes, PICRUSt-based metagenomes predict seasonality in transporters, photosynthetic proteins, peptidases and carbohydrate metabolic pathways in spite of closely related summer- and winter-associated taxa. Thus, even small temperature shifts, such as those predicted by climate change models, could affect both the structure and function of marine ecosystems.

Metapopulation theory identifies biogeographical patterns among core and satellite marine bacteria scaling from tens to thousands of kilometers

Metapopulation theory developed in terrestrial ecology provides applicable frameworks for interpreting the role of local and regional processes in shaping species distribution patterns. Yet, empirical testing of metapopulation models on microbial communities is essentially lacking. We determined regional bacterioplankton dynamics from monthly transect sampling in the Baltic Sea Proper using 16S rRNA gene sequencing. A strong positive trend was found between local relative abundance and occupancy of populations. Notably, the occupancy-frequency distributions were significantly bimodal with a satellite mode of rare endemic populations and a core mode of abundant cosmopolitan populations (e.g. Synechococcus, SAR11 and SAR86 clade members). Temporal changes in population distributions supported several theoretical frameworks. Still, bimodality was found among bacterioplankton communities across the entire Baltic Sea, and was also frequent in globally distributed datasets. Datasets spanning waters with widely different physicochemical characteristics or environmental gradients typically lacked significant bimodal patterns. When such datasets were divided into subsets with coherent environmental conditions, bimodal patterns emerged, highlighting the importance of positive feedbacks between local abundance and occupancy within specific biomes. Thus, metapopulation theory applied to microbial biogeography can provide novel insights into the mechanisms governing shifts in biodiversity resulting from natural or anthropogenically induced changes in the environment.

Regional synchrony in full-scale activated sludge bioreactors due to deterministic microbial community assembly

Seasonal community structure and regionally synchronous population dynamics have been observed in natural microbial ecosystems, but have not been well documented in wastewater treatment bioreactors. Few studies of community dynamics in full-scale activated sludge systems facing similar meteorological conditions have been done to compare the importance of deterministic and neutral community assembly mechanisms. We subjected weekly activated sludge samples from six regional full-scale bioreactors at four wastewater treatment plants obtained over 1 year to Illumina sequencing of 16S ribosomal RNA genes, resulting in a library of over 17 million sequences. All samples derived from reactors treating primarily municipal wastewater. Despite variation in operational characteristics and location, communities displayed temporal synchrony at the individual operational taxonomic unit (OTU), broad phylogenetic affiliation and community-wide scale. Bioreactor communities were dominated by 134 abundant and highly regionally synchronized OTU populations that accounted for over 50% of the total reads. Non-core OTUs displayed abundance-dependent population synchrony. Alpha diversity varied by reactor, but showed a highly reproducible and synchronous seasonal fluctuation. Community similarity was dominated by seasonal changes, but individual reactors maintained minor stable differences after 1 year. Finally, the impacts of mass migration driven by direct biomass transfers between reactors was investigated, but had no significant effect on community similarity or diversity in the sink community. Our results show that population dynamics in activated sludge bioreactors are consistent with niche-driven assembly guided by seasonal temperature fluctuations.

Potential Environmental Factors Affecting Oil-Degrading Bacterial Populations in Deep and Surface Waters of the Northern Gulf of Mexico

Understanding bacterial community dynamics as a result of an oil spill is important for predicting the fate of oil released to the environment and developing bioremediation strategies in the Gulf of Mexico. In this study, we aimed to elucidate the roles of temperature, water chemistry (nutrients), and initial bacterial community in selecting oil degraders through a series of incubation experiments. Surface (2 m) and bottom (1537 m) waters, collected near the Deepwater Horizon site, were amended with 200 ppm light Louisiana sweet crude oil and bacterial inoculums from surface or bottom water, and incubated at 4 or 24°C for 50 days. Bacterial community and residual oil were analyzed by pyrosequencing and gas chromatography-mass spectrometry (GC-MS), respectively. The results showed that temperature played a key role in selecting oil-degrading bacteria. Incubation at 4°C favored the development of Cycloclasticus, Pseudoalteromonas, Sulfitobacter, and Reinekea, while 24°C incubations enhanced Oleibacter, Thalassobius, Phaeobacter, and Roseobacter. Water chemistry and the initial community also had potential roles in the development of hydrocarbon-degrading bacterial communities. Pseudoalteromonas, Oleibacter, and Winogradskyella developed well in the nutrient-enriched bottom water, while Reinekea and Thalassobius were favored by low-nutrient surface water. We revealed that the combination of 4°C, crude oil and bottom inoculum was a key factor for the growth of Cycloclasticus, while the combination of surface inoculum and bottom water chemistry was important for the growth of Pseudoalteromonas. Moreover, regardless of the source of inoculum, bottom water at 24°C was a favorable condition for Oleibacter. Redundancy analysis further showed that temperature and initial community explained 57 and 19% of the variation observed, while oil and water chemistry contributed 14 and 10%, respectively. Overall, this study revealed the relative roles of temperature, water chemistry, and initial bacterial community in selecting oil degraders and regulating their evolution in the northern Gulf of Mexico.

Bacterial abundance and diversity in pond water supplied with different feeds

The abundance and diversity of bacteria in two types of ponds were investigated by quantitative PCR and Illumina MiSeq sequencing. The results revealed that the abundance of bacterial 16S rRNA genes in D ponds (with grass carp fed sudan grass) was significantly lower than that in E ponds (with grass carp fed commercial feed). The microbial communities were dominated by Proteobacteria, Cyanobacteria, Bacteroidetes, and Actinobacteria in both E and D ponds, while the abundance of some genera was significantly different between the two types of ponds. Specifically, some potential pathogens such as Acinetobacter and Aeromonas were found to be significantly decreased, while some probiotics such as Comamonadaceae unclassified and Bacillales unclassified were significantly increased in D ponds. In addition, water quality of D ponds was better than that of E ponds. Temperature, dissolved oxygen and nutrients had significant influence on bacterial communities. The differences in bacterial community compositions between the two types of ponds could be partially explained by the different water conditions.

Coupling Bacterioplankton Populations and Environment to Community Function in Coastal Temperate Waters

Bacterioplankton play a key role in marine waters facilitating processes important for carbon cycling. However, the influence of specific bacterial populations and environmental conditions on bacterioplankton community performance remains unclear. The aim of the present study was to identify drivers of bacterioplankton community functions, taking into account the variability in community composition and environmental conditions over seasons, in two contrasting coastal systems. A Least Absolute Shrinkage and Selection Operator (LASSO) analysis of the biological and chemical data obtained from surface waters over a full year indicated that specific bacterial populations were linked to measured functions. Namely, Synechococcus (Cyanobacteria) was strongly correlated with protease activity. Both function and community composition showed seasonal variation. However, the pattern of substrate utilization capacity could not be directly linked to the community dynamics. The overall importance of dissolved organic matter (DOM) parameters in the LASSO models indicate that bacterioplankton respond to the present substrate landscape, with a particular importance of nitrogenous DOM. The identification of common drivers of bacterioplankton community functions in two different systems indicates that the drivers may be of broader relevance in coastal temperate waters.

Riverine Bacterial Communities Reveal Environmental Disturbance Signatures within the Betaproteobacteria and Verrucomicrobia

Riverine bacterial communities play an essential role in the biogeochemical coupling of terrestrial and marine environments, transforming elements and organic matter in their journey from land to sea. However, precisely due to the fact that rivers receive significant terrestrial input, the distinction between resident freshwater taxa vs. land-derived microbes can often become ambiguous. Furthermore, ecosystem perturbations could introduce allochthonous microbial groups and reshape riverine bacterial communities. Using full- and partial-length 16S ribosomal RNA gene sequences, we analyzed the composition of bacterial communities in the Tar River of North Carolina from November 2010 to November 2011, during which a natural perturbation occurred: the inundation of the lower reaches of an otherwise drought-stricken river associated with Hurricane Irene, which passed over eastern North Carolina in late August 2011. This event provided the opportunity to examine the microbiological, hydrological, and geochemical impacts of a disturbance, defined here as the large freshwater influx into the Tar River, superimposed on seasonal changes or other ecosystem variability independent of the hurricane. Our findings demonstrate that downstream communities are more taxonomically diverse and temporally variable than their upstream counterparts. More importantly, pre- vs. post-disturbance taxonomic comparison of the freshwater-dominant Betaproteobacteria class and the phylum Verrucomicrobia reveal a disturbance signature of previously undetected taxa of diverse origins. We use known traits of closely-related taxa to interpret the ecological function of disturbance-associated bacteria, and hypothesize that carbon cycling was enhanced post-disturbance in the Tar River, likely due to the flux of organic carbon into the system associated with the large freshwater pulse. Our analyses demonstrate the importance of geochemical and hydrological alterations in structuring bacterial communities, and illustrate the response of temperate riverine bacteria on fine taxonomic scales to a disturbance.

Bacterial diversity across a highly stratified ecosystem: A salt-wedge Mediterranean estuary

Highly stratified Mediterranean estuaries are unique environments where the tidal range is low and the tidal currents are almost negligible. The main characteristics of these environments are strong salinity gradients and other environmental parameters. In this study, 454 pyrosequencing of the 16S rRNA gene in combination with catalyzed reporter deposition-fluorescence in situ hybridization (CARD-FISH) was used to estimate the bacterial diversity across the Krka estuary in February and July 2013. The comparison of the data derived from these two techniques resulted in a significant but weak positive correlation (R=0.28) indicating a substantial difference in the bacterial community structure, depending on the applied method. The phytoplankton bloom observed in February was identified as one of the main factors shaping the bacterial community structure between the two environmentally contrasting sampling months. Roseobacter, Bacteroidetes and Gammaproteobacteria differed substantially between February and July. Typical freshwater bacterial classes (Actinobacteria and Betaproteobacteria) showed strong vertical distribution patterns depending on the salinity gradient. Cyanobacteria decreased in abundance in February due to competition with phytoplankton, while the SAR11 clade increased its abundance in July as a result of a better adaptation toward more oligotrophic conditions. The results provided the first detailed insight into the bacterial diversity in a highly stratified Mediterranean karstic estuary.

Longitudinal Analysis of Microbiota in Microalga Nannochloropsis salina Cultures

Large-scale open microalgae cultivation has tremendous potential to make a significant contribution to replacing petroleum-based fuels with biofuels. Open algal cultures are unavoidably inhabited with a diversity of microbes that live on, influence, and shape the fate of these ecosystems. However, there is little understanding of the resilience and stability of the microbial communities in engineered semicontinuous algal systems. To evaluate the dynamics and resilience of the microbial communities in microalgae biofuel cultures, we conducted a longitudinal study on open systems to compare the temporal profiles of the microbiota from two multigenerational algal cohorts, which include one seeded with the microbiota from an in-house culture and the other exogenously seeded with a natural-occurring consortia of bacterial species harvested from the Pacific Ocean. From these month-long, semicontinuous open microalga Nannochloropsis salina cultures, we sequenced a time-series of 46 samples, yielding 8804 operational taxonomic units derived from 9,160,076 high-quality partial 16S rRNA sequences. We provide quantitative evidence that clearly illustrates the development of microbial community is associated with microbiota ancestry. In addition, N. salina growth phases were linked with distinct changes in microbial phylotypes. Alteromonadeles dominated the community in the N. salina exponential phase whereas Alphaproteobacteria and Flavobacteriia were more prevalent in the stationary phase. We also demonstrate that the N. salina-associated microbial community in open cultures is diverse, resilient, and dynamic in response to environmental perturbations. This knowledge has general implications for developing and testing design principles of cultivated algal systems.

Biogeographic Variation in Host Range Phenotypes and Taxonomic Composition of Marine Cyanophage Isolates

Despite the important role of phages in marine systems, little is understood about how their diversity is distributed in space. Biogeographic patterns of marine phages may be difficult to detect due to their vast genetic diversity, which may not be accurately represented by conserved marker genes. To investigate the spatial biogeographic structure of marine phages, we isolated over 400 cyanophages on Synechococcus host strain WH7803 at three coastal locations in the United States (Rhode Island, Washington, and southern California). Approximately 90% of the cyanophage isolates were myoviruses, while the other 10% were podoviruses. The diversity of isolates was further characterized in two ways: (i) taxonomically, using conserved marker genes and (ii) phenotypically, by testing isolates for their ability to infect a suite of hosts, or their "host range." Because host range is a highly variable trait even among closely related isolates, we hypothesized that host range phenotypes of cyanophage isolates would vary more strongly among locations than would taxonomic composition. Instead, we found evidence for strong biogeographic variation both in taxonomic composition and host range phenotypes, with little taxonomic overlap among the three coastal regions. For both taxonomic composition and host range phenotypes, cyanophage communities from California and Rhode Island were the most dissimilar, while Washington communities exhibited similarity to each of the other two locations. These results suggest that selection imposed by spatial variation in host dynamics influence the biogeographic distribution of cyanophages.

Wind and sunlight shape microbial diversity in surface waters of the North Pacific Subtropical Gyre

Few microbial time-series studies have been conducted in open ocean habitats having low seasonal variability such as the North Pacific Subtropical Gyre (NPSG), where surface waters experience comparatively mild seasonal variation. To better describe microbial seasonal variability in this habitat, we analyzed rRNA amplicon and shotgun metagenomic data over two years at the Hawaii Ocean Time-series Station ALOHA. We postulated that this relatively stable habitat might reveal different environmental factors that influence planktonic microbial community diversity than those previously observed in more seasonally dynamic habitats. Unexpectedly, the data showed that microbial diversity at 25 m was positively correlated with average wind speed 3 to 10 days prior to sampling. In addition, microbial community composition at 25 m exhibited significant correlations with solar irradiance. Many bacterial groups whose relative abundances varied with solar radiation corresponded to taxa known to exhibit strong seasonality in other oceanic regions. Network co-correlation analysis of 25 m communities showed seasonal transitions in composition, and distinct successional cohorts of co-occurring phylogenetic groups. Similar network analyses of metagenomic data also indicated distinct seasonality in genes originating from cyanophage, and several bacterial clades including SAR116 and SAR324. At 500 m, microbial community diversity and composition did not vary significantly with any measured environmental parameters. The minimal seasonal variability in the NPSG facilitated detection of more subtle environmental influences, such as episodic wind variation, on surface water microbial diversity. Community composition in NPSG surface waters varied in response to solar irradiance, but less dramatically than reported in other ocean provinces.

Connectivity between surface and deep waters determines prokaryotic diversity in the North Atlantic Deep Water

To decipher the influence of depth stratification and surface provincialism on the dark ocean prokaryotic community composition, we sampled the major deep-water masses in the eastern North Atlantic covering three biogeographic provinces. Their diversity was evaluated using ordination and canonical analysis of 454 pyrotag sequences. Variance partitioning suggested that 16% of the variation in the bacterial community composition was based on depth stratification while 9% of the variation was due to geographic location. General linear mixed effect models showed that the community of the subsurface waters was connected to the dark ocean prokaryotic communities in different biogeographic provinces. Cluster analysis indicated that some prokaryotic taxa are specific to distinct regions in bathypelagic water masses. Taken together, our data suggest that the dark ocean prokaryotic community composition of the eastern North Atlantic is primed by the formation and the horizontal transport of water masses.

Insensitivity of Diverse and Temporally Variable Particle-Associated Microbial Communities to Bulk Seawater Environmental Parameters

There is a growing recognition of the roles of marine microenvironments as reservoirs of biodiversity and as sites of enhanced biological activity and in facilitating biological interactions. Here, we examine the bacterial community inhabiting free-living and particle-associated seawater microenvironments at the Pivers Island Coastal Observatory (PICO). 16S rRNA gene libraries from monthly samples (July 2013 to August 2014) were used to identify microbes in seawater in four size fractions: >63 μm (zooplankton and large particles), 63 to 5 μm (particles), 5 to 1 μm (small particles/dividing cells), and <1 μm (free-living prokaryotes). Analyses of microbial community composition highlight the importance of the microhabitat (e.g., particle-associated versus free-living lifestyle) as communities cluster by size fraction, and the microhabitat explains more of the community variability than measured environmental parameters, including pH, particle concentration, projected daily insolation, nutrients, and temperature. While temperature is statistically associated with community changes in the <1-μm and 5- to 1-μm fractions, none of the measured bulk seawater environmental variables are statistically significant in the larger-particle-associated fractions. These results, combined with high particle-associated community variability, especially in the largest size fraction (i.e., >63 μm), suggest that particle composition, including eukaryotes and their associated microbiomes, may be an important factor in selecting for specific particle-associated bacteria.

IMPORTANCE

By comparing levels of particle-associated and free-living bacterial diversity at a coastal location over the course of 14 months, we show that bacteria associated with particles are generally more diverse and appear to be less responsive to commonly measured environmental variables than free-living bacteria. These diverse and highly variable particle-associated communities are likely driven by differences in particle substrates both within the water column at a single time point and due to seasonal changes over the course of the year.

Effects of temperature and nutrients on changes in genetic diversity of bacterioplankton communities in a semi-closed bay, South Korea

Bacterioplankton communities in a semi-closed bay (Jangmok Bay, South Korea) were analysed using a 16S rDNA multiplex 454 pyrosequencing approach. Diversity and operational taxonomic units of bacterioplankton communities in the Jangmok Bay are highest in cold water seasons and lowest in warm water ones. During cold seasons, α-proteobacteria respond rapidly to pulses of the concentration of inorganic nutrients, while γ-proteobacteria during warm water seasons are the most active type of bacterioplankton resent in the prevailing conditions, which include high dissolved organic carbon, chemical oxygen demand and primary production. Cyanobacteria, a minor group constituting 4.58% of the total bacterioplankton, are more abundant at low temperature. Flavobacteria are more abundant in nutrient-rich conditions and the abundance of this group also demonstrated a delayed decline following summer phytoplankton blooms. The pronounced seasonal oscillations in phosphorus concentration and temperature exert strong selection pressure on bacterioplankton communities.

More, smaller bacteria in response to ocean's warming?

Heterotrophic bacteria play a major role in organic matter cycling in the ocean. Although the high abundances and relatively fast growth rates of coastal surface bacterioplankton make them suitable sentinels of global change, past analyses have largely overlooked this functional group. Here, time series analysis of a decade of monthly observations in temperate Atlantic coastal waters revealed strong seasonal patterns in the abundance, size and biomass of the ubiquitous flow-cytometric groups of low (LNA) and high nucleic acid (HNA) content bacteria. Over this relatively short period, we also found that bacterioplankton cells were significantly smaller, a trend that is consistent with the hypothesized temperature-driven decrease in body size. Although decadal cell shrinking was observed for both groups, it was only LNA cells that were strongly coherent, with ecological theories linking temperature, abundance and individual size on both the seasonal and interannual scale. We explain this finding because, relative to their HNA counterparts, marine LNA bacteria are less diverse, dominated by members of the SAR11 clade. Temperature manipulation experiments in 2012 confirmed a direct effect of warming on bacterial size. Concurrent with rising temperatures in spring, significant decadal trends of increasing standing stocks (3% per year) accompanied by decreasing mean cell size (−1% per year) suggest a major shift in community structure, with a larger contribution of LNA bacteria to total biomass. The increasing prevalence of these typically oligotrophic taxa may severely impact marine food webs and carbon fluxes by an overall decrease in the efficiency of the biological pump.

Photochemical alteration of dissolved organic matter and the subsequent effects on bacterial carbon cycling and diversity

The impact of solar radiation on dissolved organic matter (DOM) derived from 3 different sources (seawater, eelgrass leaves and river water) and the effect on the bacterial carbon cycling and diversity were investigated. Seawater with DOM from the sources was first either kept in the dark or exposed to sunlight (4 days), after which a bacterial inoculum was added and incubated for 4 additional days. Sunlight exposure reduced the coloured DOM and carbon signals, which was followed by a production of inorganic nutrients. Bacterial carbon cycling was higher in the dark compared with the light treatment in seawater and river samples, while higher levels were found in the sunlight-exposed eelgrass experiment. Sunlight pre-exposure stimulated the bacterial growth efficiency in the seawater experiments, while no impact was found in the other experiments. We suggest that these responses are connected to differences in substrate composition and the production of free radicals. The bacterial community that developed in the dark and sunlight pre-treated samples differed in the seawater and river experiments. Our findings suggest that impact of sunlight exposure on the bacterial carbon transfer and diversity depends on the DOM source and on the sunlight-induced production of inorganic nutrients.