Annually reoccurring bacterial communities are predictable from ocean conditions

Estuarine gradients dictate spatiotemporal variations of microbiome networks in the Chesapeake Bay

BACKGROUND

Annually reoccurring microbial populations with strong spatial and temporal variations have been identified in estuarine environments, especially in those with long residence time such as the Chesapeake Bay (CB). However, it is unclear how microbial taxa cooccurr and how the inter-taxa networks respond to the strong environmental gradients in the estuaries.

RESULTS

Here, we constructed co-occurrence networks on prokaryotic microbial communities in the CB, which included seasonal samples from seven spatial stations along the salinity gradients for three consecutive years. Our results showed that spatiotemporal variations of planktonic microbiomes promoted differentiations of the characteristics and stability of prokaryotic microbial networks in the CB estuary. Prokaryotic microbial networks exhibited a clear seasonal pattern where microbes were more closely connected during warm season compared to the associations during cold season. In addition, microbial networks were more stable in the lower Bay (ocean side) than those in the upper Bay (freshwater side). Multivariate regression tree (MRT) analysis and piecewise structural equation modeling (SEM) indicated that temperature, salinity and total suspended substances along with nutrient availability, particulate carbon and Chl a, affected the distribution and co-occurrence of microbial groups, such as Actinobacteria, Bacteroidetes, Cyanobacteria, Planctomycetes, Proteobacteria, and Verrucomicrobia. Interestingly, compared to the abundant groups (such as SAR11, Saprospiraceae and Actinomarinaceae), the rare taxa including OM60 (NOR5) clade (Gammaproteobacteria), Micrococcales (Actinobacteria), and NS11-12 marine group (Bacteroidetes) contributed greatly to the stability of microbial co-occurrence in the Bay. Modularity and cluster structures of microbial networks varied spatiotemporally, which provided valuable insights into the 'small world' (a group of more interconnected species), network stability, and habitat partitioning/preferences.

CONCLUSION

Our results shed light on how estuarine gradients alter the spatiotemporal variations of prokaryotic microbial networks in the estuarine ecosystem, as well as their adaptability to environmental disturbances and co-occurrence network complexity and stability.

Bacterial community variations in the South China Sea driven by different chemical conditions

In this study, Illumina MiSeq sequencing of the 16 S rRNA gene was used to describe the bacterial communities in the South China Sea (SCS) during the southwest monsoon period. We targeted different regions in the SCS and showed that bacterial community was driven by the effects of the river, upwelling, and mesoscale eddy through changing the environmental factors (salinity, temperature, and nutrients). Distinct bacterial communities were observed among different chemical conditions, especially between the estuary and the open sea. The abundance of Burkholderiales, Frankiales, Flavobacteriales, and Rhodobacterales dominated the estuary and its adjacent waters. Bacteria in cyclonic eddy were dominated by Methylophilales and Pseudomonadales, whereas Prochlorococcus, SAR11 clade, and Oceanospirillales had relatively high abundance in the anticyclonic eddy. Overall, the abundance of specific phylotypes significantly varied among samples with different chemical conditions. Chemical conditions probably act as a driver that shapes and controls the diversity of bacteria in the SCS. This study suggests that the interaction between microbial and environmental conditions needs to be further considered to fully understand the diversity and function of marine microbes.

Prokaryotic responses to a warm temperature anomaly in northeast subarctic Pacific waters

Recent studies on marine heat waves describe water temperature anomalies causing changes in food web structure, bloom dynamics, biodiversity loss, and increased plant and animal mortality. However, little information is available on how water temperature anomalies impact prokaryotes (bacteria and archaea) inhabiting ocean waters. This is a nontrivial omission given their integral roles in driving major biogeochemical fluxes that influence ocean productivity and the climate system. Here we present a time-resolved study on the impact of a large-scale warm water surface anomaly in the northeast subarctic Pacific Ocean, colloquially known as the Blob, on prokaryotic community compositions. Multivariate statistical analyses identified significant depth- and season-dependent trends that were accentuated during the Blob. Moreover, network and indicator analyses identified shifts in specific prokaryotic assemblages from typically particle-associated before the Blob to taxa considered free-living and chemoautotrophic during the Blob, with potential implications for primary production and organic carbon conversion and export.

Seasonal and annual changes in the microbial communities of Ofunato Bay, Japan, based on metagenomics

Five years of datasets from 2015 to 2019 of whole genome shotgun sequencing for cells trapped on 0.2-µm filters of seawater collected monthly from Ofunato Bay, an enclosed bay in Japan, were analysed, which included the 2015 data that we had reported previously. Nucleotide sequences were determined for extracted DNA from three locations for both the upper (1 m) and deeper (8 or 10 m) depths. The biotic communities analysed at the domain level comprised bacteria, eukaryotes, archaea and viruses. The relative abundance of bacteria was over 60% in most months for the five years. The relative abundance of the SAR86 cluster was highest in the bacterial group, followed by Candidatus Pelagibacter and Planktomarina. The relative abundance of Ca. Pelagibacter showed no relationship with environmental factors, and those of SAR86 and Planktomarina showed positive correlations with salinity and dissolved oxygen, respectively. The bacterial community diversity showed seasonal changes, with high diversity around September and low diversity around January for all five years. Nonmetric multidimensional scaling analysis also revealed that the bacterial communities in the bay were grouped in a season-dependent manner and linked with environmental variables such as seawater temperature, salinity and dissolved oxygen.

Dietary Exposure of Pacific Oyster (Crassostrea gigas) Larvae to Compromised Microalgae Results in Impaired Fitness and Microbiome Shift

The Pacific oyster Crassostrea gigas is the world's most cultivated oyster and seed supply is heavily reliant on hatchery production where recurring mass mortality events are a major constraint. Outbreaks of bacterial infection via microalgal feed are frequently implicated in these mortalities. This study assessed the effects of feeding compromised microalgae to developing oyster larvae. Intentionally 'stressed' (high pH) or non-stressed microalgae were fed to 11 day-old oyster larvae at two feeding rations for 96 h, followed by a recovery period. Biological endpoints of larval performance were measured following the 96 h exposure and subsequent recovery. Bacterial communities associated with the microalgae feed, rearing seawater, and the oyster larvae, were characterized and correlated with effects on oyster fitness parameters. Feeding stressed algae to oyster larvae for 96 h increased the occurrence of deformities (>70% vs. 20% in control), reduced feeding and swimming ability, and slowed development. Following the recovery period, fewer larvae reached pediveliger stage (2.7% vs. 36% in control) and became spat (1.5% vs. 6.6% in control). The quantity of stressed algae supplied to oyster larvae also influenced overall larval performance, with high feeding rations generally causing greater impairment than low rations. Bacterial profiling using 16S rRNA showed that most bacterial families characterized in larval tissue were also present in larval rearing seawater and in the microalgae feed (98%). The rearing seawater showed the highest bacterial richness compared to the larval and the microalgal compartments, regardless of feeding regime. In larval tissue, bacterial richness was highest in stressed and high-feed treatments, and negatively correlated with larval fitness parameters. These results suggest significant dysbiosis induced by compromised feed and/or increased feed ration. Several bacterial genera (e.g., Halomonas, Marinomonas) were strongly associated with impaired larval performance while the presence of genera in larvae including Vibrio was closely associated with overfeeding. Our research demonstrated that metabarcoding can be effectively used to identify microbiota features associated with larval fitness.

A Combined Digital PCR and Next Generation DNA-Sequencing Based Approach for Tracking Nearshore Pollutant Dynamics Along the Southwest United States/Mexico Border

Ocean currents, multiple fecal bacteria input sources, and jurisdictional boundaries can complicate pollution source tracking and associated mitigation and management efforts within the nearshore coastal environment. In this study, multiple microbial source tracking tools were employed to characterize the impact and reach of an ocean wastewater treatment facility discharge in Mexico northward along the coast and across the Southwest United States- Mexico Border. Water samples were evaluated for fecal indicator bacteria (FIB), Enterococcus by culture-based methods, and human-associated genetic marker (HF183) and Enterococcus by droplet digital polymerase chain reaction (ddPCR). In addition, 16S rRNA gene sequence analysis was performed and the SourceTracker algorithm was used to characterize the bacterial community of the wastewater treatment plume and its contribution to beach waters. Sampling dates were chosen based on ocean conditions associated with northern currents. Evidence of a gradient in human fecal pollution that extended north from the wastewater discharge across the United States/Mexico border from the point source was observed using human-associated genetic markers and microbial community analysis. The spatial extent of fecal contamination observed was largely dependent on swell and ocean conditions. These findings demonstrate the utility of a combination of molecular tools for understanding and tracking specific pollutant sources in dynamic coastal water environments.

Groundwater bacterial communities evolve over time in response to recharge

Time series analyses are a crucial tool for uncovering the patterns and processes shaping microbial communities and their functions, especially in aquatic ecosystems. Subsurface aquatic environments are perceived to be more stable than surface oceans and lakes, due to the lack of sunlight, the absence of photosysnthetically-driven primary production, low temperature variations, and oligotrophic conditions. However, periodic groundwater recharge should affect the structure and succession of groundwater microbiomes. To disentangle the long-term temporal changes in bacterial communities of shallow fractured bedrock groundwater, and identify the drivers of the observed patterns, we analysed bacterial 16S rRNA gene sequencing data for samples collected monthly from three groundwater wells over a six-year period (n = 230) along a hillslope recharge area. We showed that the bacterial communities in the groundwater of limestone-mudstone alternations were not stable over time and exhibited non-linear dissimilarity patterns which corresponded to periods of groundwater recharge. Further, we observed an increase in dissimilarity over time (generalized additive model P < 0.001) indicating that the successive recharge events result in communities that are increasingly more dissimilar to the initial reference time point. The sampling period was able to explain up to 29.5% of the variability in bacterial community composition and the impact of recharge events on the groundwater microbiome was linked to the strength of the recharge and local environmental selection. Many groundwater bacteria originated from the recharge-related sources (mean = 66.5%, SD = 15.1%) and specific bacterial taxa were identified as being either enriched or repressed during recharge events. Overall, similar to surface aquatic environments, the microbiomes in shallow fractured-rock groundwater vary through time, though we revealed groundwater recharges as unique driving factors for these patterns. The high temporal resolution employed here highlights the dynamics of bacterial communities in groundwater, which is an essential resource for the provision of clean drinking water; understanding the biological complexities of these systems is therefore crucial.

Seasonal dynamics of the microbial community in two full-scale wastewater treatment plants: Diversity, composition, phylogenetic group based assembly and co-occurrence pattern

The optimal operation and functional stability of a wastewater treatment plant (WWTP) strongly depend on the properties of its microbial community. However, a knowledge gap remains regarding the seasonal dynamics of microbial community properties, especially phylogenetic group based assembly and co-occurrence patterns. Accordingly, in this study, influent and activated sludge (AS) samples were weekly collected from 2 full-scale WWTPs for one year (89 influent and 103 AS samples in total) and examined by high-throughput Illumina-MiSeq sequencing. The results suggested that the microbial community diversity and composition in the influent fluctuated substantially with season, while those in the AS had a relatively more stable pattern throughout the year. The phylogenetic group based assembly mechanisms of AS community were identified by using "Infer Community Assembly Mechanisms by Phylogenetic-bin-based null model (iCAMP)". The results showed that drift accounted for the largest proportion (52.8%), while homogeneous selection (18.2%) was the most important deterministic process. Deterministic processes dominated in 47 microbial groups (bins), which were also found (~40%) in the AS core taxa dataset. Moreover, the results suggested that Nitrospira were more susceptible to stochastic processes in winter, which may provide a possible explanation for nitrification failure in winter. Network analysis results suggested that the network structure of the AS community could be more stable in summer and autumn. In addition, there were no identical keystone taxa found in different networks (constructed from different plants, sources, and seasons), which supported the context dependency theory. The results of this study deepened our understanding of the microbial ecology in AS systems and provided a foundation for further studies on the community regulation strategy of WWTPs.

The surface bacterial community of an Australian kelp shows cross-continental variation and relative stability within regions

Epiphytic microbial communities often have a close relationship with their eukaryotic host, assisting with defence, health, disease prevention and nutrient transfer. Shifts in the structure of microbial communities could therefore have negative effects on the individual host and indirectly impact the surrounding ecosystem, particularly for major habitat-forming hosts, such as kelps in temperate rocky shores. Thus, an understanding of the structure and dynamics of host-associated microbial communities is essential for monitoring and assessing ecosystem changes. Here, samples were taken from the ecologically important kelp, Ecklonia radiata, over a 17-month period, from six different sites in two distinct geographic regions (East and West coasts of Australia), separated by ∼3,300 kms, to understand variation in the kelp bacterial community and its potential environmental drivers. Differences were observed between kelp bacterial communities between the largely disconnected geographical regions. In contrast, within each region and over time the bacterial communities were considerably more stable, despite substantial seasonal changes in environmental conditions.

Bacterial communities in temperate and polar coastal sands are seasonally stable

Coastal sands are biocatalytic filters for dissolved and particulate organic matter of marine and terrestrial origin, thus, acting as centers of organic matter transformation. At high temporal resolution, we accessed the variability of benthic bacterial communities over two annual cycles at Helgoland (North Sea), and compared it with seasonality of communities in Isfjorden (Svalbard, 78°N) sediments, where primary production does not occur during winter. Benthic community structure remained stable in both, temperate and polar sediments on the level of cell counts and 16S rRNA-based taxonomy. Actinobacteriota of uncultured Actinomarinales and Microtrichales were a major group, with 8 ± 1% of total reads (Helgoland) and 31 ± 6% (Svalbard). Their high activity (frequency of dividing cells 28%) and in situ cell numbers of >10% of total microbes in Svalbard sediments, suggest Actinomarinales and Microtrichales as key heterotrophs for carbon mineralization. Even though Helgoland and Svalbard sampling sites showed no phytodetritus-driven changes of the benthic bacterial community structure, they harbored significantly different communities (p < 0.0001, r = 0.963). The temporal stability of benthic bacterial communities is in stark contrast to the dynamic succession typical of coastal waters, suggesting that pelagic and benthic bacterial communities respond to phytoplankton productivity very differently.

Urban wastewater bacterial communities assemble into seasonal steady states

BACKGROUND

Microorganisms in urban sanitary sewers exhibit community properties that suggest sewers are a novel ecosystem. Sewer microorganisms present both an opportunity as a control point for wastewater treatment and a risk to human health. If treatment processes are to be improved and health risks quantified, then it is necessary to understand microbial distributions and dynamics within this community. Here, we use 16S rRNA gene sequencing to characterize raw influent wastewater bacterial communities in a 5-year time series from two wastewater treatment plants in Milwaukee, WI; influent wastewater from 77 treatment plants across the USA; and wastewater in 12 Milwaukee residential sewers.

RESULTS

In Milwaukee, we find that in transit from residences to treatment plants, the human bacterial component of wastewater decreases in proportion and exhibits stochastic temporal variation. In contrast, the resident sewer community increases in abundance during transit and cycles seasonally according to changes in wastewater temperature. The result is a bacterial community that assembles into two distinct community states each year according to the extremes in wastewater temperature. Wastewater bacterial communities from other northern US cities follow temporal trends that mirror those in Milwaukee, but southern US cities have distinct community compositions and differ in their seasonal patterns.

CONCLUSIONS

Our findings provide evidence that environmental conditions associated with seasonal change and climatic differences related to geography predictably structure the bacterial communities residing in below-ground sewer pipes. Video abstract.

Role of Phylogenetic Structure in the Dynamics of Coastal Viral Assemblages

Marine microbes, including viruses, are an essential part of the marine ecosystem, forming the base of the food web and driving biogeochemical cycles. Within this system, the composition of viral assemblages changes markedly with time, and some of these changes are repeatable through time; however, the extent to which these dynamics are reflected within versus among evolutionarily related groups of viruses is largely unexplored. To examine these dynamics, changes in the composition of two groups of ecologically important viruses and communities of their potential hosts were sampled every 2 weeks for 13 months at a coastal site in British Columbia, Canada. We sequenced two marker genes for viruses-the gene encoding the major capsid protein of T4-like phages and their relatives (gp23) and the RNA-dependent RNA polymerase (RdRp) gene of marnavirus-like RNA viruses-as well as marker genes for their bacterial and eukaryotic host communities, the genes encoding 16S rRNA and 18S rRNA. There were strong lagged correlations between viral diversity and community similarity of putative hosts, implying that the viruses influenced the composition of the host communities. The results showed that for both viral assemblages, the dominant clusters of phylogenetically related viruses shifted over time, and this was correlated with environmental changes. Viral clusters contained many ephemeral taxa and few persistent taxa, but within a viral assemblage, the ephemeral and persistent taxa were closely related, implying ecological dynamics within these clusters. Furthermore, these dynamics occurred in both the RNA and DNA viral assemblages surveyed, implying that this structure is common in natural viral assemblages.IMPORTANCE Viruses are major agents of microbial mortality in marine systems, yet little is known about changes in the composition of viral assemblages in relation to those of the microbial communities that they infect. Here, we sampled coastal seawater every 2 weeks for 1 year and used high-throughput sequencing of marker genes to follow changes in the composition of two groups of ecologically important viruses, as well as the communities of bacteria and protists that serve as their respective hosts. Different subsets of genetically related viruses dominated at different times. These results demonstrate that although the genetic composition of viral assemblages is highly dynamic temporally, for the most part the shuffling of genotypes occurs within a few clusters of phylogenetically related viruses. Thus, it appears that even in temperate coastal waters with large seasonal changes, the highly dynamic shuffling of viral genotypes occurs largely within a few subsets of related individuals.

Temperate and tropical coastal waters share relatively similar microbial biofilm communities while free-living or particle-attached communities are distinct

Free-living (FL) marine microbial communities differ from those attached to particles (PA). Likewise, biofilms (B) colonizing artificial surfaces, including plastics or ship hulls, hardly resemble their planktonic surroundings. However, few studies have examined the effect of the environment on these lifestyles and on the source of organisms colonizing marine surfaces. Using 16S rRNA gene metabarcoding, we identified specificities of marine prokaryotic community lifestyles (FL, PA or B) sampled in three coastal polluted locations with dissimilar environmental conditions: the North-Western Mediterranean Sea and the Atlantic and Indian Oceans. Biofilms developed over polyvinyl chloride (PVC) were found to be significantly different from FL or PA collected during the immersions. Alpha-diversity increased from FL to PA and to B, illustrating the integrative aspect of the latter, with little proportion of operational taxonomic units shared with the first two. Beta-diversity clustered first the lifestyles and then the sites. FL and PA were more affected by water quality, especially by trace metal contamination, whereas B were as sensitive to trace metals as to nutrients. Although biofilms should be supplied by the planktonic (ultra) rare biosphere, source tracking could only detect small contributions of FL or PA taxa to B communities.

Seasonal marine microorganisms change neighbours under contrasting environmental conditions

Marine picoplankton contribute to global carbon sequestration and nutrient recycling. These processes are directly related to the composition of communities, which in turn depends on microbial interactions and environmental forcing. Under regular seasonal cycles, marine communities show strong predictable patterns of annual re-occurrences, but little is known about the effect of environmental perturbation on their organization. The aim of our study was to investigate the co-occurrence patterns of planktonic picoeukaryote, bacteria and archaea under contrasting environmental conditions. The study was designed to have high sampling frequency that could match both the biological rhythm of marine microbes and the short time scale of extreme weather events. Our results show that microbial networks changed from year to year depending on conditions. In addition, individual taxa became less interconnected and changed neighbours, which revealed an unfaithful relationship between marine microorganisms. This unexpected pattern suggests possible switches between organisms that have similar specific functions, or hints at the presence of organisms that share similar environmental niches without interacting. Despite the observed annual changes, the time series showed re-occurring communities that appear to recover from perturbations. Changing co-occurrence patterns between marine microorganisms may allow the long-term stability of ecosystems exposed to contrasting meteorological events.

Diversity of Bacterioplankton and Bacteriobenthos from the Veracruz Reef System, Southwestern Gulf of Mexico

Bacterial diversity was explored among field samples and cultured isolates from coral reefs within the Veracruz Reef System. Bacterioplankton and bacteriobenthos were characterized by pyrosequencing 16S rRNA genes. Identified sequences belonged to the kingdom Bacteria and classified into 33 phyla. Proteobacteria (likely SAR11 clade) dominated in collective field samples, whereas Firmicutes were the most abundant taxa among cultured isolates. Bioinformatic sorting of sequences to family level revealed 223 bacterial families. Pseudomonadaceae, Exiguobacteraceae and Bacillaceae were dominant among cultured isolates. Vibrionaceae, Alteromonadaceae, and Flavobacteriaceae dominated in reef-associated sediments, whereas Rickettsiaceae and Synechoccaceae were more highly represented in the water column. Bacterial communities from sediments were more diverse than from the water column. This study reveals cryptic bacterial diversity among microenvironmental components of marine microbial reef communities subject to differential influence of anthropogenic stressors. Such investigations are critical for constructing scenarios of environmentally induced shifts in bacterial biodiversity and species composition.

Inter- and Intra-Annual Bacterioplankton Community Patterns in a Deepwater Sub-Arctic Region: Persistent High Background Abundance of Putative Oil Degraders

Hydrocarbon-degrading bacteria naturally degrade and remove petroleum pollutants, yet baselines do not currently exist for these critical microorganisms in many regions where the oil and gas industry is active. Furthermore, understanding how a baseline community changes across the seasons and its potential to respond to an oil spill event are prerequisites for predicting their response to elevated hydrocarbon exposures. In this study, 16S rRNA gene-based profiling was used to assess the spatiotemporal variability of baseline bacterioplankton community composition in the Faroe-Shetland Channel (FSC), a deepwater sub-Arctic region where the oil and gas industry has been active for the last 40 years. Over a period of 2 years, we captured the diversity of the bacterioplankton community within distinct water masses (defined by their temperature and salinity) that have a distinct geographic origin (Atlantic or Nordic), depth, and direction of flow. We demonstrate that bacterioplankton communities were significantly different across water samples of contrasting origin and depth. Taxa of known hydrocarbon-degrading bacteria were observed at higher-than-anticipated abundances in water masses originating in the Nordic Seas, suggesting these organisms are sustained by an unconfirmed source of oil input in that region. In the event of an oil spill, our results suggest that the response of these organisms is severely hindered by the low temperatures and nutrient levels that are typical for the FSC.IMPORTANCE Oil spills at sea are one of the most disastrous anthropogenic pollution events, with the Deepwater Horizon spill providing a testament to how profoundly the health of marine ecosystems and the livelihood of its coastal inhabitants can be severely impacted by spilled oil. The fate of oil in the environment is largely dictated by the presence and activities of natural communities of oil-degrading bacteria. While a significant effort was made to monitor and track the microbial response and degradation of the oil in the water column in the wake of the Deepwater Horizon spill, the lack of baseline data on the microbiology of the Gulf of Mexico confounded scientists' abilities to provide an accurate assessment of how the system responded relative to prespill conditions. This data gap highlights the need for long-term microbial ocean observatories in regions at high risk of oil spills. Here, we provide the first microbiological baseline established for a subarctic region experiencing high oil and gas industry activity, the northeast Atlantic, but with no apparent oil seepage or spillage. We also explore the presence, relative abundances, and seasonal dynamics of indigenous hydrocarbon-degrading communities. These data will advance the development of models to predict the behavior of such organisms in the event of a major oil spill in this region and potentially impact bioremediation strategies by enhancing the activities of these organisms in breaking down the oil.

Extensive Microbial Processing of Polysaccharides in the South Pacific Gyre via Selfish Uptake and Extracellular Hydrolysis

Primary productivity occurs throughout the deep euphotic zone of the oligotrophic South Pacific Gyre (SPG), fueled largely by the regeneration of nutrients and thus recycling of organic matter. We investigated the heterotrophic capabilities of the SPG's bacterial communities by examining their ability to process polysaccharides, an important component of marine organic matter. We focused on the initial step of organic matter degradation by measuring the activities of extracellular enzymes that hydrolyze six different polysaccharides to smaller sizes. This process can occur by two distinct mechanisms: "selfish uptake," in which initial hydrolysis is coupled to transport of large polysaccharide fragments into the periplasmic space of bacteria, with little to no loss of hydrolysis products to the external environment, and "external hydrolysis," in which low molecular weight (LMW) hydrolysis products are produced in the external environment. Given the oligotrophic nature of the SPG, we did not expect high enzymatic activity; however, we found that all six polysaccharides were hydrolyzed externally and taken up selfishly in the central SPG, observations that may be linked to a comparatively high abundance of diatoms at the depth and location sampled (75 m). At the edge of the gyre and close to the center of the gyre, four of six polysaccharides were externally hydrolyzed, and a lower fraction of the bacterial community showed selfish uptake. One polysaccharide (fucoidan) was selfishly taken up without measurable external hydrolysis at two stations. Additional incubations of central gyre water from depths of 1,250 and 2,800 m with laminarin (an abundant polysaccharide in the ocean) led to extreme growth of opportunistic bacteria (Alteromonas), as tracked by cell counts and next generation sequencing of the bacterial communities. These Alteromonas appear to concurrently selfishly take up laminarin and release LMW hydrolysis products. Overall, extracellular enzyme activities in the SPG were similar to activities in non-oligotrophic regions, and a considerable fraction of the community was capable of selfish uptake at all three stations. A diverse set of bacteria responded to and are potentially important for the recycling of organic matter in the SPG.

The Evolutionary Success of the Marine Bacterium SAR11 Analyzed through a Metagenomic Perspective

As the most abundant bacteria in oceans, the Pelagibacterales order (here SAR11) plays an important role in the global carbon cycle, but the study of the evolutionary forces driving its evolution has lagged considerably due to the inherent difficulty of obtaining pure cultures. Multiple evolutionary models have been proposed to explain the diversification of distinct lineages within a population; however, the identification of many of these patterns in natural populations remains mostly enigmatic. We have used a metagenomic approach to explore microdiversity patterns in their natural habitats. Comparison with a collection of bacterial and archaeal groups from the same environments shows that SAR11 populations have a different evolutionary regime, where multiple genotypes coexist within the same population and remain stable over time. Widespread homologous recombination could be one of the main driving factors of this homogenization.

The SAR11 clade of Alphaproteobacteria is the most abundant group of planktonic cells in the near-surface epipelagic waters of the ocean, but the mechanisms underlying its exceptional success have not been fully elucidated. Here, we applied a metagenomic approach to explore microdiversity patterns by measuring the accumulation of synonymous and nonsynonymous mutations as well as homologous recombination in populations of SAR11 from different aquatic habitats (marine epipelagic, bathypelagic, and surface freshwater). The patterns of mutation accumulation and recombination were compared to those of other groups of representative marine microbes with multiple ecological strategies that share the same marine habitat, namely, Cyanobacteria (Prochlorococcus and Synechococcus), Archaea (“Candidatus Nitrosopelagicus” and Marine Group II Thalassoarchaea), and some heterotrophic marine bacteria (Alteromonas and Erythrobacter). SAR11 populations showed widespread recombination among distantly related members, preventing divergence leading to a genetically stable population. Moreover, their high intrapopulation sequence diversity with an enrichment in synonymous replacements supports the idea of a very ancient divergence and the coexistence of multiple different clones. However, other microbes analyzed seem to follow different evolutionary dynamics where processes of diversification driven by geographic and ecological instability produce a higher number of nonsynonymous replacements and lower intrapopulation sequence diversity. Together, these data shed light on some of the evolutionary and ecological processes that lead to the large genomic diversity in SAR11. Furthermore, this approach can be applied to other similar microbes that are difficult to culture in the laboratory, but abundant in nature, to investigate the underlying dynamics of their genomic evolution.

IMPORTANCE As the most abundant bacteria in oceans, the Pelagibacterales order (here SAR11) plays an important role in the global carbon cycle, but the study of the evolutionary forces driving its evolution has lagged considerably due to the inherent difficulty of obtaining pure cultures. Multiple evolutionary models have been proposed to explain the diversification of distinct lineages within a population; however, the identification of many of these patterns in natural populations remains mostly enigmatic. We have used a metagenomic approach to explore microdiversity patterns in their natural habitats. Comparison with a collection of bacterial and archaeal groups from the same environments shows that SAR11 populations have a different evolutionary regime, where multiple genotypes coexist within the same population and remain stable over time. Widespread homologous recombination could be one of the main driving factors of this homogenization.

Seasonality in Spatial Turnover of Bacterioplankton Along an Ecological Gradient in the East China Sea: Biogeographic Patterns, Processes and Drivers

Seasonal succession in bacterioplankton is a common process in marine waters. However, seasonality in their spatial turnover is largely unknown. Here, we investigated spatial turnover of surface bacterioplankton along a nearshore-to-offshore gradient in the East China Sea across four seasons. Although seasonality overwhelmed spatial variability of bacterioplankton composition, we found significant spatial turnover of bacterioplankton along the gradient as well as overall seasonal consistency in biogeographic patterns (including distance-decay relationship and covariation of community composition with distance to shore) with subtle changes. Bacterioplankton assembly was consistently dominated by deterministic mechanisms across seasons, with changes in specific processes. We found overall seasonal consistency in abiotic factors (mainly salinity and nitrogen and phosphorus nutrients) shaping bacterioplankton composition, while phytoplankton showed a similar influence as abiotic factors only in spring. Although key taxa responsible for bacterioplankton spatial turnover showed certain season-specificity, seasonal switching between closely related taxa occurred within most dominant families. Moreover, many close relatives showed different responding patterns to the environmental gradients in different seasons, suggesting their differences in both seasonally climatic and spatially environmental preferences. Our results provide insights into seasonal consistency and variability in spatial turnover of bacterioplankton in terms of biogeographic patterns, ecological processes, and external and internal drivers.

Predetermined clockwork microbial worlds: Current understanding of aquatic microbial diel response from model systems to complex environments

In the photic zone of aquatic ecosystems, microorganisms with different metabolisms and their viruses form complex interactions and food webs. Within these interactions, phototrophic microorganisms such as eukaryotic microalgae and cyanobacteria interact directly with sunlight, and thereby generate circadian rhythms. Diel cycling originally generated in microbial phototrophs is directly transmitted toward heterotrophic microorganisms utilizing the photosynthetic products as they are excreted or exuded. Such diel cycling seems to be indirectly propagated toward heterotrophs as a result of complex biotic interactions. For example, cell death of phototrophic microorganisms induced by viral lysis and protistan grazing provides additional resources of dissolved organic matter to the microbial community, and so generates diel cycling in other heterotrophs with different nutrient dependencies. Likewise, differences in the diel transmitting pathway via complex interactions among heterotrophs, and between heterotrophs and their viruses, may also generate higher variation and time lag diel rhythms in different heterotrophic taxa. Thus, sunlight and photosynthesis not only contribute energy and carbon supply, but also directly or indirectly control diel cycling of the microbial community through complex interactions in the photic zone of aquatic ecosystems.

Persistent El Niño driven shifts in marine cyanobacteria populations

In the California Current Ecosystem, El Niño acts as a natural phenomenon that is partially representative of climate change impacts on marine bacteria at timescales relevant to microbial communities. Between 2014–2016, the North Pacific warm anomaly (a.k.a., the “blob”) and an El Niño event resulted in prolonged ocean warming in the Southern California Bight (SCB). To determine whether this “marine heatwave” resulted in shifts in microbial populations, we sequenced the rpoC1 gene from the biogeochemically important picocyanobacteria Prochlorococcus and Synechococcus at 434 time points from 2009–2018 in the MICRO time series at Newport Beach, CA. Across the time series, we observed an increase in the abundance of Prochlorococcus relative to Synechococcus as well as elevated frequencies of ecotypes commonly associated with low-nutrient and high-temperature conditions. The relationships between environmental and ecotype trends appeared to operate on differing temporal scales. In contrast to ecotype trends, most microdiverse populations were static and possibly reflect local habitat conditions. The only exceptions were microdiversity from Prochlorococcous HLI and Synechococcus Clade II that shifted in response to the 2015 El Niño event. Overall, Prochlorococcus and Synechococcus populations did not return to their pre-heatwave composition by the end of this study. This research demonstrates that extended warming in the SCB can result in persistent changes in key microbial populations.

Environmental stability impacts the differential sensitivity of marine microbiomes to increases in temperature and acidity

Ambient conditions shape microbiome responses to both short- and long-duration environment changes through processes including physiological acclimation, compositional shifts, and evolution. Thus, we predict that microbial communities inhabiting locations with larger diel, episodic, and annual variability in temperature and pH should be less sensitive to shifts in these climate-change factors. To test this hypothesis, we compared responses of surface ocean microbes from more variable (nearshore) and more constant (offshore) sites to short-term factorial warming (+3 °C) and/or acidification (pH -0.3). In all cases, warming alone significantly altered microbial community composition, while acidification had a minor influence. Compared with nearshore microbes, warmed offshore microbiomes exhibited larger changes in community composition, phylotype abundances, respiration rates, and metatranscriptomes, suggesting increased sensitivity of microbes from the less-variable environment. Moreover, while warming increased respiration rates, offshore metatranscriptomes yielded evidence of thermal stress responses in protein synthesis, heat shock proteins, and regulation. Future oceans with warmer waters may enhance overall metabolic and biogeochemical rates, but they will host altered microbial communities, especially in relatively thermally stable regions of the oceans.

Seasonal dynamics of prokaryotes and their associations with diatoms in the Southern Ocean as revealed by an autonomous sampler

The Southern Ocean remains one of the least explored marine environments. The investigation of temporal microbial dynamics has thus far been hampered by the limited access to this remote ocean. We present here high-resolution seasonal observations of the prokaryotic community composition during phytoplankton blooms induced by natural iron fertilization. A total of 18 seawater samples were collected by a moored remote autonomous sampler over 4 months at 5-11 day intervals in offshore surface waters (central Kerguelen Plateau). Illumina sequencing of the 16S rRNA gene revealed that among the most abundant amplicon sequence variants, SAR92 and Aurantivirga were the first bloom responders, Pseudomonadaceae, Nitrincolaceae and Polaribacter had successive peaks during the spring bloom decline, and Amylibacter increased in relative abundance later in the season. SAR11 and SUP05 were abundant prior to and after the blooms. Using network analysis, we identified two groups of diatoms representative of the spring and summer bloom that had opposite correlation patterns with prokaryotic taxa. Our study provides the first seasonal picture of microbial community dynamics in the open Southern Ocean and thereby offers biological insights to the cycling of carbon and iron, and to an important puzzling issue that is the modest nitrate decrease associated to iron fertilization.

Seasonal Variation of Bacterial Diversity Along the Marine Particulate Matter Continuum

Seasonal dynamics of ocean prokaryotic communities in the free-living fraction have been widely described, but less is known about the seasonality of prokaryotes inhabiting marine particles. We describe the seasonality of bacterial communities in the particulate matter continuum by sampling monthly over two years in a temperate oligotrophic coastal ecosystem and using a serial filtration (including six size-fractions spanning from 0.2 to 200 μm). We observed that bacterial communities in the particulate matter continuum had annual changes following harmonic seasonal oscillations, where alpha, beta, and gamma diversity increased during the warm period and decreased during the cold period. Communities in each size-fraction changed gradually over time, being the communities in larger size-fractions the ones with stronger annual changes. Annual community changes were driven mainly by day length and sea surface temperature, and each size-fraction was additionally affected by other variables (e.g., smaller size-fractions by dissolved PO₄ and larger size-fractions by turbidity). While some taxonomic groups mantained their preference for a given size fraction during most of the year, others varied their distribution into different size fractions over time, as e.g., SAR11, which increased its presence in particles during the cold period. Our results indicate that the size-fractionation scheme provides novel seasonal patterns that are not possible to unveil by analyzing only free-living bacteria, and that help to better understand the temporal dynamics of prokaryotes.

Depth-Differentiation and Seasonality of Planktonic Microbial Assemblages in the Monterey Bay Upwelling System

Coastal upwelling regions are hotspots of biological productivity, supporting diverse communities of microbial life and metabolisms. Monterey Bay (MB), a coastal ocean embayment in central California, experiences seasonal upwelling of cold, nutrient-rich waters that sustain episodes of high phytoplankton production in surface waters. While productivity in surface waters is intimately linked to metabolisms of diverse communities of Archaea and Bacteria, a comprehensive understanding of the microbial community in MB is missing thus far, particularly in relation to the distinct hydrographic seasons characteristic of the MB system. Here we present the results of a 2-year microbial time-series survey in MB, investigating community composition and structure across spatiotemporal gradients. In deciphering these patterns, we used unique sequence variants (SVs) of the 16S rRNA gene (V4–V5 region), complemented with metagenomes and metatranscriptomes representing multiple depth profiles. We found clear depth-differentiation and recurring seasonal abundance patterns within planktonic communities, particularly when analyzed at finer taxonomic levels. Compositional changes were more pronounced in the upper 0–40 m of the water column, whereas deeper depths were characterized by temporally stable populations. In accordance with the dynamic nutrient profiles, the system appears to change from a Bacteroidetes- and Rhodobacterales-dominated upwelling period to an oceanic season dominated by oligotrophic groups such as SAR11 and picocyanobacteria. The cascade of environmental changes brought about by upwelling and relaxation events thus impacts microbial community structure in the bay, with important implications for the temporal variability of nutrient and energy fluxes within the MB ecosystem. Our observations emphasize the need for continued monitoring of planktonic microbial communities in order to predict and manage the behavior of this sensitive marine sanctuary ecosystem, over projected intensification of upwelling in the region.

Tara Oceans: towards global ocean ecosystems biology

A planetary-scale understanding of the ocean ecosystem, particularly in light of climate change, is crucial. Here, we review the work of Tara Oceans, an international, multidisciplinary project to assess the complexity of ocean life across comprehensive taxonomic and spatial scales. Using a modified sailing boat, the team sampled plankton at 210 globally distributed sites at depths down to 1,000 m. We describe publicly available resources of molecular, morphological and environmental data, and discuss how an ecosystems biology approach has expanded our understanding of plankton diversity and ecology in the ocean as a planetary, interconnected ecosystem. These efforts illustrate how global-scale concepts and data can help to integrate biological complexity into models and serve as a baseline for assessing ecosystem changes and the future habitability of our planet in the Anthropocene epoch.

Timescales of variation in diversity and production of bacterioplankton assemblages in the Lower Mississippi River

Rivers are characterized by rapid and continuous one-way directional fluxes of flowing, aqueous habitat, chemicals, suspended particles, and resident plankton. Therefore, at any particular location in such systems there is the potential for continuous, and possibly abrupt, changes in diversity and metabolic activities of suspended biota. As microorganisms are the principal catalysts of organic matter degradation and nutrient cycling in rivers, examination of their assemblage dynamics is fundamental to understanding system-level biogeochemical patterns and processes. However, there is little known of the dynamics of microbial assemblage composition or production of large rivers along a time interval gradient. We quantified variation in alpha and beta diversity and production of particle-associated and free-living bacterioplankton assemblages collected at a single site on the Lower Mississippi River (LMR), the final segment of the largest river system in North America. Samples were collected at timescales ranging from days to weeks to months up to a year. For both alpha and beta diversity, there were similar patterns of temporal variation in particle-associated and free-living assemblages. Alpha diversity, while always higher on particles, varied as much at a daily as at a monthly timescale. Beta diversity, in contrast, gradually increased with time interval of sampling, peaking between samples collected 180 days apart, before gradually declining between samples collected up to one year apart. The primary environmental driver of the temporal pattern in beta diversity was temperature, followed by dissolved nitrogen and chlorophyll a concentrations. Particle-associated bacterial production corresponded strongly to temperature, while free-living production was much lower and constant over time. We conclude that particle-associated and free-living bacterioplankton assemblages of the LMR vary in richness, composition, and production at distinct timescales in response to differing sets of environmental factors. This is the first temporal longitudinal study of microbial assemblage structure and dynamics in the LMR.

Deforestation impacts network co-occurrence patterns of microbial communities in Amazon soils

Co-occurrence networks allow for the identification of potential associations among species, which may be important for understanding community assembly and ecosystem functions. We employed this strategy to examine prokaryotic co-occurrence patterns in the Amazon soils and the response of these patterns to land use change to pasture, with the hypothesis that altered microbial composition due to deforestation will mirror the co-occurrence patterns across prokaryotic taxa. In this study, we calculated Spearman correlations between operational taxonomic units (OTUs) as determined by 16S rRNA gene sequencing, and only robust correlations were considered for network construction (-0.80 ≥ P ≥ 0.80, adjusted P < 0.01). The constructed network represents distinct forest and pasture components, with altered compositional and topological features. A comparative analysis between two representative modules of these contrasting ecosystems revealed novel information regarding changes to metabolic pathways related to nitrogen cycling. Our results showed that soil physicochemical properties such as temperature, C/N and H++Al3+ had a significant impact on prokaryotic communities, with alterations to network topologies. Taken together, changes in co-occurrence patterns and physicochemical properties may contribute to ecosystem processes including nitrification and denitrification, two important biogeochemical processes occurring in tropical forest systems.

Prokaryotic Population Dynamics and Viral Predation in a Marine Succession Experiment Using Metagenomics

We performed an incubation experiment of seawater confined in plastic bottles with samples collected at three depths (15, 60, and 90 m) after retrieval from a single offshore location in the Mediterranean Sea, from a late summer stratified water column. Two samples representative of each depth were collected and stored in opaque bottles after two periods of 7 h. We took advantage of the "bottle effect" to investigate changes in the natural microbial communities (abundant and rare). We recovered 94 metagenome-assembled genomes (MAGs) and 1089 metagenomic viral contigs and examined their abundance using metagenomic recruitment. We detected a significant fast growth of copiotrophic bacteria such as Alteromonas or Erythrobacter throughout the entire water column with different dynamics that we assign to "clonal," "polyclonal," or "multispecies" depending on the recruitment pattern. Results also showed a marked ecotype succession in the phototropic picocyanobacteria that were able to grow at all the depths in the absence of light, highlighting the importance of their mixotrophic potential. In addition, "wall-chain-reaction" hypothesis based on the study of phage-host dynamics showed the higher impact of viral predation on archaea in deeper waters, evidencing their prominent role during incubations. Our results provide a step forward in understanding the mechanisms underlying dynamic patterns and ecology of the marine microbiome and the importance of processing the samples immediately after collection to avoid changes in the community structure.

Community response of microbial primary producers to salinity is primarily driven by nutrients in lakes

Higher microbial diversity was frequently observed in saline than fresh waters, but the underlying mechanisms remains unknown, particularly in microbial primary producers (MPP). MPP abundance and activity are notably constrained by high salinity, but facilitated by high nutrients. It remains to be ascertained whether and how nutrients regulate the salinity constraints on MPP abundance and community structure. Here we investigated the impact of nutrients on salinity constraints on MPP abundance and diversity in undisturbed lakes with a wide salinity range on the Tibetan Plateau. MPP community was explored using quantitative PCR, terminal restriction fragment length polymorphism and sequencing of cloning libraries targeting form IC cbbL gene. The MPP community structure was sorted by salinity into freshwater (salinity<1‰), saline (1‰ < salinity<29‰) and hypersaline (salinity>29‰) lakes. Furthermore, while MPP abundance, diversity and richness were significantly constrained with increasing salinity, these constraints were mitigated by enhancing total organic carbon (TOC) and total nitrogen (TN) contents in freshwater and saline lakes. In contrast, the MPP diversity increased significantly with the salinity in hypersaline lakes, due to the mitigation of enhancing TOC and TN contents and salt-tolerant MPP taxa. The mitigating effect of nutrients was more pronounced in saline than in freshwater and hypersaline lakes. The MPP compositions varied along salinity, with Betaproteobacteria dominating both the freshwater and saline lakes and Gammaproteobacteria dominating the hypersaline lakes. We concluded that high nutrients could mitigate the salinity constraining effects on MPP abundance, community richness and diversity. Our findings offer a novel insight into the salinity effects on primary producers and highlight the interactive effects of salinity and nutrients on MPP in lakes. These findings can be used as a baseline to illuminate the effects of increased anthropogenic activities altering nutrient dynamics on the global hydrological cycle and the subsequent responses thereof by MPP communities.

Gene Expression Changes and Community Turnover Differentially Shape the Global Ocean Metatranscriptome

Ocean microbial communities strongly influence the biogeochemistry, food webs, and climate of our planet. Despite recent advances in understanding their taxonomic and genomic compositions, little is known about how their transcriptomes vary globally. Here, we present a dataset of 187 metatranscriptomes and 370 metagenomes from 126 globally distributed sampling stations and establish a resource of 47 million genes to study community-level transcriptomes across depth layers from pole-to-pole. We examine gene expression changes and community turnover as the underlying mechanisms shaping community transcriptomes along these axes of environmental variation and show how their individual contributions differ for multiple biogeochemically relevant processes. Furthermore, we find the relative contribution of gene expression changes to be significantly lower in polar than in non-polar waters and hypothesize that in polar regions, alterations in community activity in response to ocean warming will be driven more strongly by changes in organismal composition than by gene regulatory mechanisms. VIDEO ABSTRACT.

Recurrent seasonal changes in bacterial growth efficiency, metabolism and community composition in coastal waters

The microbial response to environmental changes in coastal waters of the eastern Cantabrian Sea was explored for four years by analysing a broad set of environmental variables along with bacterial community metabolism and composition. A recurrent seasonal cycle emerged, consisting of two stable periods, characterized by low bacterial metabolic activity (winter) from October to March, and high bacterial metabolic activity (summer) from May to August. These two contrasting periods were linked by short transition periods in April (T_A ) and September (T_S ). The phylogenetic groups Alphaproteobacteria and Bacteroidetes were dominant during winter and summer respectively, and their recurrent alternation was mainly driven by the bloom of eukaryotic phytoplankton before T_A and the bloom of prokaryotic phytoplankton before T_S . Bacterial growth efficiency remained high and stable during the winter and summer periods but dropped during the two short transition periods. Our results suggest that bacterial growth efficiency should be considered a very resilient property that reflects different stages in the adaptation of the bacterial community composition to the environmental changes occurring throughout the seasonal cycle in this coastal ecosystem.

Transient exposure to novel high temperatures reshapes coastal phytoplankton communities

Average sea surface temperatures are expected to rise 4° this century, and marine phytoplankton and bacterial community composition, biogeochemical rates, and trophic interactions are all expected to change in a future warmer ocean. Thermal experiments typically use constant temperatures; however, weather and hydrography cause marine temperatures to fluctuate on diel cycles and over multiple days. We incubated natural communities of phytoplankton collected from California coastal waters during spring, summer, and fall under present-day and future mean temperatures, using thermal treatments that were either constant or fluctuated on a 48 h cycle. As assayed by marker-gene sequencing, the emergent microbial communities were consistent within each season, except when culture temperatures exceeded the highest temperature recorded in a 10-year local thermal dataset. When temperature treatments exceeded the 10-year maximum the phytoplankton community shifted, becoming dominated by diatom amplicon sequence variants (ASVs) not seen at lower temperatures. When mean temperatures were above the 10-year maximum, constant and fluctuating regimes each selected for different ASVs. These findings suggest coastal microbial communities are largely adapted to the current range of temperatures they experience. They also suggest a general hypothesis whereby multiyear upper temperature limits may represent thresholds, beyond which large community restructurings may occur. Now inevitable future temperature increases that exceed these environmental thresholds, even temporarily, may fundamentally reshape marine microbial communities and therefore the biogeochemical cycles that they mediate.

High bacterial diversity in nearshore and oceanic biofilms and their influence on larval settlement by Hydroides elegans (Polychaeta)

Settlement of many benthic marine invertebrates is stimulated by bacterial biofilms, although it is not known if patterns of settlement reflect microbial communities that are specific to discrete habitats. Here, we characterized the taxonomic and functional gene diversity (16S rRNA gene amplicon and metagenomic sequencing analyses), as well as the specific bacterial abundances, in biofilms from diverse nearby and distant locations, both inshore and offshore, and tested them for their ability to induce settlement of the biofouling tubeworm Hydroides elegans, an inhabitant of bays and harbours around the world. We found that compositions of the bacterial biofilms were site specific, with the greatest differences between inshore and offshore sites. Further, biofilms were highly diverse in their taxonomic and functional compositions across inshore sites, while relatively low diversity was found at offshore sites. Hydroides elegans settled on all biofilms tested, with settlement strongly correlated with bacterial abundance. Bacterial density in biofilms was positively correlated with biofilm age. Our results suggest that the localized distribution of H. elegans is not determined by 'selection' to locations by specific bacteria, but it is more likely linked to the prevailing local ecology and oceanographic features that affect the development of dense biofilms and the occurrence of larvae.

Multi-year dynamics of fine-scale marine cyanobacterial populations are more strongly explained by phage interactions than abiotic, bottom-up factors

Currently defined ecotypes in marine cyanobacteria Prochlorococcus and Synechococcus likely contain subpopulations that themselves are ecologically distinct. We developed and applied high-throughput sequencing for the 16S-23S rRNA internally transcribed spacer (ITS) to examine ecotype and fine-scale genotypic community dynamics for monthly surface water samples spanning 5 years at the San Pedro Ocean Time-series site. Ecotype-level structure displayed regular seasonal patterns including succession, consistent with strong forcing by seasonally varying abiotic parameters (e.g. temperature, nutrients, light). We identified tens to thousands of amplicon sequence variants (ASVs) within ecotypes, many of which exhibited distinct patterns over time, suggesting ecologically distinct populations within ecotypes. Community structure within some ecotypes exhibited regular, seasonal patterns, but not for others, indicating other more irregular processes such as phage interactions are important. Network analysis including T4-like phage genotypic data revealed distinct viral variants correlated with different groups of cyanobacterial ASVs including time-lagged predator-prey relationships. Variation partitioning analysis indicated that phage community structure more strongly explains cyanobacterial community structure at the ASV level than the abiotic environmental factors. These results support a hierarchical model whereby abiotic environmental factors more strongly shape niche partitioning at the broader ecotype level while phage interactions are more important in shaping community structure of fine-scale variants within ecotypes.

Diel population and functional synchrony of microbial communities on coral reefs

On coral reefs, microorganisms are essential for recycling nutrients to primary producers through the remineralization of benthic-derived organic matter. Diel investigations of reef processes are required to holistically understand the functional roles of microbial players in these ecosystems. Here we report a metagenomic analysis characterizing microbial communities in the water column overlying 16 remote forereef sites over a diel cycle. Our results show that microbial community composition is more dissimilar between day and night samples collected from the same site than between day or night samples collected across geographically distant reefs. Diel community differentiation is largely driven by the flux of Psychrobacter sp., which is two-orders of magnitude more abundant during the day. Nighttime communities are enriched with species of Roseobacter, Halomonas, and Alteromonas encoding a greater variety of pathways for carbohydrate catabolism, further illustrating temporal patterns of energetic provisioning between different marine microbes. Dynamic diel fluctuations of microbial populations could also support the efficient trophic transfer of energy posited in coral reef food webs.

Composition and Dynamics of the Activated Sludge Microbiome during Seasonal Nitrification Failure

Wastewater treatment plants in temperate climate zones frequently undergo seasonal nitrification failure in the winter month yet maintain removal efficiency for other contaminants. We tested the hypothesis that nitrification failure can be correlated to shifts in the nitrifying microbial community. We monitored three parallel, full-scale sequencing batch reactors over the course of a year with respect to reactor performance, microbial community composition via 16S rRNA gene amplicon sequencing, and functional gene abundance using qPCR. All reactors demonstrated similar changes to their core microbiome, and only subtle variations among seasonal and transient taxa. We observed a decrease in species richness during the winter, with a slow recovery of the activated sludge community during spring. Despite the change in nitrification performance, ammonia monooxygenase gene abundances remained constant throughout the year, as did the relative sequence abundance of Nitrosomonadacae. This suggests that nitrification failure at colder temperatures might result from different reaction kinetics of nitrifying taxa, or that other organisms with strong seasonal shifts in population abundance, e.g. an uncultured lineage of Saprospiraceae, affect plant performance in the winter. This research is a comprehensive analysis of the seasonal microbial community dynamics in triplicate full-scale sequencing batch reactors and ultimately strengthens our basic understanding of the microbial ecology of activated sludge communities by revealing seasonal succession patterns of individual taxa that correlate with nutrient removal efficiency.

Constructing the Microbial Association Network from Large-Scale Time Series Data Using Granger Causality

The increasing availability of large-scale time series data allows the inference of microbial community dynamics by association network analysis. However, correlation-based association network analyses are noninformative of causal, mediating and time-dependent relationships between microbial community functional factors. To address this insufficiency, we introduced the Granger causality model to the analysis of a recent marine microbial time series dataset. We systematically constructed a directed acyclic network, representing both internal and external causal relationships among the microbial and environmental factors. We further optimized the network by removing false causal associations using the conditional Granger causality. The final network was visualized as a Granger graph, which was analyzed to identify causal relationships driven by key functional operators in the environment, such as Gammaproteobacteria, which was Granger caused by total organic nitrogen and primary production (p < 0.05 and Q < 0.05).

Changes in Bacterioplankton Communities Resulting From Direct and Indirect Interactions With Trace Metal Gradients in an Urbanized Marine Coastal Area

Unraveling the relative importance of both environmental conditions and ecological processes regulating bacterioplankton communities is a central goal in microbial ecology. Marine coastal environments are among the most urbanized areas and as a consequence experience environmental pressures. The highly anthropized Toulon Bay (France) was considered as a model system to investigate shifts in bacterioplankton communities along natural and anthropogenic physicochemical gradients during a 1-month survey. In depth geochemical characterization mainly revealed strong and progressive Cd, Zn, Cu, and Pb contamination gradients between the entrance of the Bay and the north-western anthropized area. On the other hand, low-amplitude natural gradients were observed for other environmental variables. Using 16S rRNA gene sequencing, we observed strong spatial patterns in bacterioplankton taxonomic and predicted function structure along the chemical contamination gradient. Variation partitioning analysis demonstrated that multiple metallic contamination explained the largest part of the spatial biological variations observed, but DOC and salinity were also significant contributors. Network analysis revealed that biotic interactions were far more numerous than direct interactions between microbial groups and environmental variables. This suggests indirect effects of the environment, and especially trace metals, on the community through a few taxonomic groups. These spatial patterns were also partially found for predicted bacterioplankton functions, thus indicating a limited functional redundancy. All these results highlight both potential direct influences of trace metals contamination on coastal bacterioplankton and indirect forcing through biotic interactions and cascading.

Probabilistic Modeling of Microbial Metabolic Networks for Integrating Partial Quantitative Knowledge Within the Nitrogen Cycle

Understanding the interactions between microbial communities and their environment sufficiently to predict diversity on the basis of physicochemical parameters is a fundamental pursuit of microbial ecology that still eludes us. However, modeling microbial communities is problematic, because (i) communities are complex, (ii) most descriptions are qualitative, and (iii) quantitative understanding of the way communities interact with their surroundings remains incomplete. One approach to overcoming such complications is the integration of partial qualitative and quantitative descriptions into more complex networks. Here we outline the development of a probabilistic framework, based on Event Transition Graph (ETG) theory, to predict microbial community structure across observed chemical data. Using reverse engineering, we derive probabilities from the ETG that accurately represent observations from experiments and predict putative constraints on communities within dynamic environments. These predictions can feedback into the future development of field experiments by emphasizing the most important functional reactions, and associated microbial strains, required to characterize microbial ecosystems.

Redox traits characterize the organization of global microbial communities

Biological communities are conventionally described as assemblages of species, whose ecological roles are known or predictable from their observable morphology. In microbial ecology, such a taxonomic approach is hindered by limited capacity to discriminate among different microbes, which bear highly dynamic genomes and establish complex associations. Approaches based on culture-independent functional genes profiling might overcome these problems, but a set of usable established genes in a general situation is still lacking. We show that genes related to reduction-oxidation (redox) processes separate microbial communities into their corresponding biomes. This redox-based characterization is linked to the microbial energetics of ecosystems and to most biogeochemical cycles and might be useful for assessing the impact of environmental degradation on the ecosystem services, underpinned by microorganisms.

The structure of biological communities is conventionally described as profiles of taxonomic units, whose ecological functions are assumed to be known or, at least, predictable. In environmental microbiology, however, the functions of a majority of microorganisms are unknown and expected to be highly dynamic and collectively redundant, obscuring the link between taxonomic structure and ecosystem functioning. Although genetic trait-based approaches at the community level might overcome this problem, no obvious choice of gene categories can be identified as appropriate descriptive units in a general ecological context. We used 247 microbial metagenomes from 18 biomes to determine which set of genes better characterizes the differences among biomes on the global scale. We show that profiles of oxidoreductase genes support the highest biome differentiation compared with profiles of other categories of enzymes, general protein-coding genes, transporter genes, and taxonomic gene markers. Based on oxidoreductases’ description of microbial communities, the role of energetics in differentiation and particular ecosystem function of different biomes become readily apparent. We also show that taxonomic diversity is decoupled from functional diversity, e.g., grasslands and rhizospheres were the most diverse biomes in oxidoreductases but not in taxonomy. Considering that microbes underpin biogeochemical processes and nutrient recycling through oxidoreductases, this functional diversity should be relevant for a better understanding of the stability and conservation of biomes. Consequently, this approach might help to quantify the impact of environmental stressors on microbial ecosystems in the context of the global-scale biome crisis that our planet currently faces.

Rhythmicity of coastal marine picoeukaryotes, bacteria and archaea despite irregular environmental perturbations

Seasonality in marine microorganisms has been classically observed in phytoplankton blooms, and more recently studied at the community level in prokaryotes, but rarely investigated at the scale of individual microbial taxa. Here we test if specific marine eukaryotic phytoplankton, bacterial and archaeal taxa display yearly rhythms at a coastal site impacted by irregular environmental perturbations. Our seven-year study in the Bay of Banyuls (North Western Mediterranean Sea) shows that despite some fluctuating environmental conditions, many microbial taxa displayed significant yearly rhythms. The robust rhythmicity was found in both autotrophs (picoeukaryotes and cyanobacteria) and heterotrophic prokaryotes. Sporadic meteorological events and irregular nutrient supplies did, however, trigger the appearance of less common non-rhythmic taxa. Among the environmental parameters that were measured, the main drivers of rhythmicity were temperature and day length. Seasonal autotrophs may thus be setting the pace for rhythmic heterotrophs. Similar environmental niches may be driving seasonality as well. The observed strong association between Micromonas and SAR11, which both need thiamine precursors for growth, could be a first indication that shared nutritional niches may explain some rhythmic patterns of co-occurrence.

High Frequency Multi-Year Variability in Baltic Sea Microbial Plankton Stocks and Activities

Marine bacterioplankton are essential in global nutrient cycling and organic matter turnover. Time-series analyses, often at monthly sampling frequencies, have established the paramount role of abiotic and biotic variables in structuring bacterioplankton communities and productivities. However, fine-scale seasonal microbial activities, and underlying biological principles, are not fully understood. We report results from four consecutive years of high-frequency time-series sampling in the Baltic Proper. Pronounced temporal dynamics in most investigated microbial variables were observed, including bacterial heterotrophic production, plankton biomass, extracellular enzyme activities, substrate uptake rate constants of glucose, pyruvate, acetate, amino acids, and leucine, as well as nutrient limitation bioassays. Spring blooms consisting of diatoms and dinoflagellates were followed by elevated bacterial heterotrophic production and abundances. During summer, bacterial productivity estimates increased even further, coinciding with an initial cyanobacterial bloom in early July. However, bacterial abundances only increased following a second cyanobacterial bloom, peaking in August. Uptake rate constants for the different measured carbon compounds varied seasonally and inter-annually and were highly correlated to bacterial productivity estimates, temperature, and cyanobacterial abundances. Further, we detected nutrient limitation in response to environmental conditions in a multitude of microbial variables, such as elevated productivities in nutrient bioassays, changes in enzymatic activities, or substrate preferences. Variations among biotic variables often occurred on time scales of days to a few weeks, yet often spanning several sampling occasions. Such dynamics might not have been captured by sampling at monthly intervals, as compared to more predictable transitions in abiotic variables such as temperature or nutrient concentrations. Our study indicates that high resolution analyses of microbial biomass and productivity parameters can help out in the development of biogeochemical and food web models disentangling the microbial black box.

Bacterioplankton Metacommunity Processes across Thermal Gradients: Weaker Species Sorting but Stronger Niche Segregation in Summer than in Winter in a Subtropical Bay

Thermal effluents from nuclear power plants greatly change the environmental and ecological conditions of the receiving marine water body, but knowledge about their impact on microbial ecology is limited. Here we used high-throughput sequencing of the 16S rRNA gene to examine marine bacterioplankton metacommunity assembly across thermal gradients in two representative seasons (i.e., winter and summer) in a subtropical bay located on the northern coast of the South China Sea. We found high heterogeneity in bacterioplankton community compositions (BCCs) across thermal gradients and between seasons. The spatially structured temperature gradient created by thermal effluents was the key determinant of BCCs, but its influence differed by season. Using a metacommunity approach, we found that in the thermal discharge area, i.e., where water is frequently exchanged with surrounding seawater and thermal effluent water, the BCC spatial patterns were shaped by species sorting rather than by mass effects from surrounding seawater or by dilution of thermal effluent water by surrounding seawater. However, this effect of species sorting was weaker in summer than in winter seawater. In both seasons, the bacterioplankton community structure was predominately determined by niche sharing; however, the relative importance of niche segregation was enhanced in summer seawater. Our findings suggest that for the seasonal differences in metacommunity processes, the BCCs of subtropical summer seawater were more sensitive to temperature and were more difficult to predict than those of winter seawater in the face of different intensities of thermal impacts.IMPORTANCE Understanding the mechanisms of bacterial community assembly across environmental gradients is one of the major goals of marine microbial ecology. Thermal effluents from two nuclear power plants have been present in the subtropical Daya Bay for more than 20 years and have generated a comparatively stable and long thermal gradient (a temperature increase from 0 to 10°C). The environmental patches across thermal gradients are heterogeneous and very strongly interconnected on a microbial scale; thus, this is a useful model for the study of the metacommunity processes (i.e., patch dynamics, species sorting, mass effects, and neutral processes) that underlie marine bacterioplankton assembly. The significance of our research is to reveal how environmental conditions and dispersal-related processes interact to influence bacterioplankton metacommunity processes and their seasonal differences across thermal gradients. Our results may advance the understanding of marine microbial ecology under future conditions of global warming.