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

Timescales and drivers of chlorophyll variability in a subtropical, long residence time estuary (Baffin Bay, Texas, USA)

Little is known about the role of short-term (hourly to daily) processes on estuarine phytoplankton dynamics, though these can be important timescales of physical-chemical variability in coastal systems. This study coupled high-frequency (daily; 2015-2016) and low frequency (monthly; 2014-2018) sampling to quantify timescales of chlorophyll variability and potential drivers of phytoplankton blooms in a subtropical, long residence time estuary (Baffin Bay, Texas, USA). The dominant timescale of chlorophyll variability in the system was sub-monthly, which accounted for ~ 37% of variability, followed by interannual at ~ 30% and seasonal at 21%. From bloom events identified in the high frequency dataset, wind was a predominant factor related to short-term (daily resolution) blooms, often positively correlated with chlorophyll concentrations. Results expand upon the currently limited knowledge of short-term chlorophyll variability in subtropical, long residence time estuaries. Additionally, findings also offer insight into design considerations for sampling programs in similar systems.

Microbial functional guilds respond cohesively to rapidly fluctuating environments

Microbial communities experience environmental fluctuations across timescales from rapid changes in moisture, temperature, or light levels to long-term seasonal or climactic variations. Understanding how microbial populations respond to these changes is critical for predicting the impact of perturbations, interventions, and climate change on communities. Since communities typically harbor tens to hundreds of distinct taxa, the response of microbial abundances to perturbations is potentially complex. However, while taxonomic diversity is high, in many communities taxa can be grouped into functional guilds of strains with similar metabolic traits. These guilds effectively reduce the complexity of the system by providing a physiologically motivated coarse-graining. Here, using a combination of simulations, theory, and experiments, we show that the response of guilds to nutrient fluctuations depends on the timescale of those fluctuations. Rapid changes in nutrient levels drive cohesive, positively correlated abundance dynamics within guilds. For slower timescales of environmental variation, members within a guild begin to compete due to similar resource preferences, driving negative correlations in abundances between members of the same guild. Our results provide a route to understanding the relationship between functional guilds and community response to changing environments, as well as an experimental approach to discovering functional guilds via designed nutrient perturbations to communities.

Seasonal recurrence and modular assembly of an Arctic pelagic marine microbiome

Deciphering how microbial communities are shaped by environmental variability is fundamental for understanding the structure and function of ocean ecosystems. While seasonal environmental gradients have been shown to structure the taxonomic dynamics of microbiomes over time, little is known about their impact on functional dynamics and the coupling between taxonomy and function. Here, we demonstrate annually recurrent, seasonal structuring of taxonomic and functional dynamics in a pelagic Arctic Ocean microbiome by combining autonomous samplers and in situ sensors with long-read metagenomics and SSU ribosomal metabarcoding. Specifically, we identified five temporal microbiome modules whose succession within each annual cycle represents a transition across different ecological states. For instance, Cand. Nitrosopumilus, Syndiniales, and the machinery to oxidise ammonia and reduce nitrite are signatures of early polar night, while late summer is characterised by Amylibacter and sulfur compound metabolism. Leveraging metatranscriptomes from Tara Oceans, we also demonstrate the consistency in functional dynamics across the wider Arctic Ocean during similar temporal periods. Furthermore, the structuring of genetic diversity within functions over time indicates that environmental selection pressure acts heterogeneously on microbiomes across seasons. By integrating taxonomic, functional and environmental information, our study provides fundamental insights into how microbiomes are structured under pronounced seasonal changes in understudied, yet rapidly changing polar marine ecosystems.

Gene content of seawater microbes is a strong predictor of water chemistry across the Great Barrier Reef

BACKGROUND

Seawater microbes (bacteria and archaea) play essential roles in coral reefs by facilitating nutrient cycling, energy transfer, and overall reef ecosystem functioning. However, environmental disturbances such as degraded water quality and marine heatwaves, can impact these vital functions as seawater microbial communities experience notable shifts in composition and function when exposed to stressors. This sensitivity highlights the potential of seawater microbes to be used as indicators of reef health. Microbial indicator analysis has centered around measuring the taxonomic composition of seawater microbial communities, but this can obscure heterogeneity of gene content between taxonomically similar microbes, and thus, microbial functional genes have been hypothesized to have more scope for predictive potential, though empirical validation for this hypothesis is still pending. Using a metagenomics study framework, we establish a functional baseline of seawater microbiomes across offshore Great Barrier Reef (GBR) sites to compare the diagnostic value between taxonomic and functional information in inferring continuous physico-chemical metrics in the surrounding reef.

RESULTS

Integrating gene-centric metagenomics analyses with 17 physico-chemical variables (temperature, salinity, and particulate and dissolved nutrients) across 48 reefs revealed that associations between microbial functions and environmental parameters were twice as stable compared to taxonomy-environment associations. Distinct seasonal variations in surface water chemistry were observed, with nutrient concentrations up to threefold higher during austral summer, explained by enhanced production of particulate organic matter (POM) by photoautotrophic picocyanobacteria, primarily Synechococcus. In contrast, nutrient levels were lower in winter, and POM production was also attributed to Prochlorococcus. Additionally, heterotrophic microbes (e.g., Rhodospirillaceae, Burkholderiaceae, Flavobacteriaceae, and Rhodobacteraceae) were enriched in reefs with elevated dissolved organic carbon (DOC) and phytoplankton-derived POM, encoding functional genes related to membrane transport, sugar utilization, and energy metabolism. These microbes likely contribute to the coral reef microbial loop by capturing and recycling nutrients derived from Synechococcus and Prochlorococcus, ultimately transferring nutrients from picocyanobacterial primary producers to higher trophic levels.

CONCLUSION

This study reveals that functional information in reef-associated seawater microbes more robustly associates with physico-chemical variables than taxonomic data, highlighting the importance of incorporating microbial function in reef monitoring initiatives. Our integrative approach to mine for stable seawater microbial biomarkers can be expanded to include additional continuous metrics of reef health (e.g., benthic cover of corals and macroalgae, fish counts/biomass) and may be applicable to other large-scale reef metagenomics datasets beyond the GBR. Video Abstract.

Inferring resource competition in microbial communities from time series

The competition for resources is a defining feature of microbial communities. In many contexts, from soils to host-associated communities, highly diverse microbes are organized into metabolic groups or guilds with similar resource preferences. The resource preferences of individual taxa that give rise to these guilds are critical for understanding fluxes of resources through the community and the structure of diversity in the system. However, inferring the metabolic capabilities of individual taxa, and their competition with other taxa, within a community is challenging and unresolved. Here we address this gap in knowledge by leveraging dynamic measurements of abundances in communities. We show that simple correlations are often misleading in predicting resource competition. We show that spectral methods such as the cross-power spectral density (CPSD) and coherence that account for time-delayed effects are superior metrics for inferring the structure of resource competition in communities. We first demonstrate this fact on synthetic data generated from consumer-resource models with time-dependent resource availability, where taxa are organized into groups or guilds with similar resource preferences. By applying spectral methods to oceanic plankton time-series data, we demonstrate that these methods detect interaction structures among species with similar genomic sequences. Our results indicate that analyzing temporal data across multiple timescales can reveal the underlying structure of resource competition within communities.

Dynamic microbiome diversity shaping the adaptation of sponge holobionts in coastal waters

The microbial communities associated with sponges contribute to the adaptation of hosts to environments, which are essential for the trophic transformation of benthic-marine coupling. However, little is known about the symbiotic microbial community interactions and adaptative strategies of high- and low-microbial abundance (HMA and LMA) sponges, which represent two typical ecological phenotypes. Here, we compared the 1-year dynamic patterns of microbiomes with the HMA sponge Spongia officinalis and two LMA sponge species Tedania sp. and Haliclona simulans widespread on the coast of China. Symbiotic bacterial communities with the characteristic HMA-LMA dichotomy presented higher diversity and stability in S. officinalis than in Tedania sp. and H. simulans, while archaeal communities showed consistent diversity across all sponges throughout the year. Dissolved oxygen, dissolved inorganic phosphorus, dissolved organic phosphorus, and especially temperature were the major factors affecting the seasonal changes in sponge microbial communities. S. officinalis-associated microbiome had higher diversity, stronger stability, and closer interaction, which adopted a relatively isolated strategy to cope with environmental changes, while Tedania sp. and H. simulans were more susceptible and shared more bacterial Amplicon Sequence Variants (ASVs) with surrounding waters, with an open way facing the uncertainty of the environment. Meta-analysis of the microbiome in composition, diversity, and ecological function from 13 marine sponges further supported that bacterial communities associated with HMA and LMA sponges have evolved two distinct environmental adaptation strategies. We propose that the different adaptive ways of sponges responding to the environment may be responsible for their successful evolution and their competence in global ocean change.

IMPORTANCE

During long-term evolution, sponge holobionts, among the oldest symbiotic relationships between microbes and metazoans, developed two distinct phenotypes with high- and low-microbial abundance (HMA and LMA). Despite sporadic studies indicating that the characteristic microbial assemblages present in HMA and LMA sponges, the adaptation strategies of symbionts responding to environments are still unclear. This deficiency limits our understanding of the selection of symbionts and the ecological functions during the evolutionary history and the adaptative assessment of HMA and LMA sponges in variable environments. Here, we explored symbiotic communities with two distinct phenotypes in a 1-year dynamic environment and combined with the meta-analysis of 13 sponges. The different strategies of symbionts in adapting to the environment were basically drawn: microbes with LMA were more acclimated to environmental changes, forming relatively loose-connected communities, while HMA developed relatively tight-connected and more similar communities beyond the divergence of species and geographical location.

Mechanism of bacterial predation via ixotrophy

Ixotrophy is a contact-dependent predatory strategy of filamentous bacteria in aquatic environments for which the molecular mechanism remains unknown. We show that predator-prey contact can be established by gliding motility or extracellular assemblages we call "grappling hooks." Cryo-electron microscopy identified the grappling hooks as heptamers of a type IX secretion system substrate. After close predator-prey contact is established, cryo-electron tomography and functional assays showed that puncturing by a type VI secretion system mediated killing. Single-cell analyses with stable isotope-labeled prey revealed that prey components are taken up by the attacker. Depending on nutrient availability, insertion sequence elements toggle the activity of ixotrophy. A marine metagenomic time series shows coupled dynamics of ixotrophic bacteria and prey. We found that the mechanism of ixotrophy involves multiple cellular machineries, is conserved, and may shape microbial populations in the environment.

Diversity in the utilization of different molecular classes of dissolved organic matter by heterotrophic marine bacteria

Heterotrophic marine bacteria utilize and recycle dissolved organic matter (DOM), impacting biogeochemical cycles. It is currently unclear to what extent distinct DOM components can be used by different heterotrophic clades. Here, we ask how a natural microbial community from the Eastern Mediterranean Sea (EMS) responds to different molecular classes of DOM (peptides, amino acids, amino sugars, disaccharides, monosaccharides, and organic acids) comprising much of the biomass of living organisms. Bulk bacterial activity increased after 24 h for all treatments relative to the control, while glucose and ATP uptake decreased or remained unchanged. Moreover, while the per-cell uptake rate of glucose and ATP decreased, that of Leucin significantly increased for amino acids, reflecting their importance as common metabolic currencies in the marine environment. Pseudoalteromonadaceae dominated the peptides treatment, while different Vibrionaceae strains became dominant in response to amino acids and amino sugars. Marinomonadaceae grew well on organic acids, and Alteromonadaseae on disaccharides. A comparison with a recent laboratory-based study reveals similar peptide preferences for Pseudoalteromonadaceae, while Alteromonadaceae, for example, grew well in the lab on many substrates but dominated in seawater samples only when disaccharides were added. We further demonstrate a potential correlation between the genetic capacity for degrading amino sugars and the dominance of specific clades in these treatments. These results highlight the diversity in DOM utilization among heterotrophic bacteria and complexities in the response of natural communities.

IMPORTANCE

A major goal of microbial ecology is to predict the dynamics of natural communities based on the identity of the organisms, their physiological traits, and their genomes. Our results show that several clades of heterotrophic bacteria each grow in response to one or more specific classes of organic matter. For some clades, but not others, growth in a complex community is similar to that of isolated strains in laboratory monoculture. Additionally, by measuring how the entire community responds to various classes of organic matter, we show that these results are ecologically relevant, and propose that some of these resources are utilized through common uptake pathways. Tracing the path between different resources to the specific microbes that utilize them, and identifying commonalities and differences between different natural communities and between them and lab cultures, is an important step toward understanding microbial community dynamics and predicting how communities will respond to perturbations.

Understanding phycosomal dynamics to improve industrial microalgae cultivation

Algal-bacterial interactions are ubiquitous in both natural and industrial systems, and the characterization of these interactions has been reinvigorated by potential applications in biosystem productivity. Different growth conditions can be used for operational functions, such as the use of low-quality water or high pH/alkalinity, and the altered operating conditions likely constrain microbial community structure and function in unique ways. However, research is necessary to better understand whether consortia can be designed to improve the productivity, processing, and sustainability of industrial-scale cultivations through different controls that can constrain microbial interactions for maximal light-driven outputs. The review highlights current knowledge and gaps for relevant operating conditions, as well as suggestions for near-term and longer-term improvements for large-scale cultivation and polyculture engineering.

Mechanisms of microbial coexistence in a patchy ecosystem: Differences in ecological niche overlap and species fitness between rhythmic and non-rhythmic species

Resource patchiness caused by external events breaks the continuity and homogeneity of resource distribution in the original ecosystem. For local organisms, this leads to drastic changes in the availability of resources, breaks down the co-existence of species, and reshuffles the local ecosystem. West Lake is a freshwater lake with resource patchiness caused by multiple exogenous disturbances that has strong environmental heterogeneity that prevents clear observation of seasonal changes in the microbial communities. Despite this, the emergence of rhythmic species in response to irregular changes in the environment has been helpful for observing microbial communities dynamics in patchy ecosystems. We investigated the ecological mechanisms of seasonal changes in microbial communities in West Lake by screening rhythmic species based on the ecological niche and modern coexistence theories. The results showed that rhythmic species were the dominant factors in microbial community changes and the effects of most environmental factors on the microbial community were indirectly realised through the rhythmic species. Random forest analyses showed that seasonal changes in the microbial community were similarly predicted by the rhythmic species. In addition, we incorporated species interactions and community phylogenetic patterns into stepwise multiple regression analyses, the results of which indicate that ecological niches and species fitness may drive the coexistence of these subcommunities. Thus, this study extends our understanding of seasonal changes in microbial communities and provides new ways for observing seasonal changes in microbial communities, especially in ecosystems with resource patches. Our study also show that combining community phylogenies with co-occurrence networks based on ecological niches and modern coexistence theory can further help us understand the ecological mechanisms of interspecies coexistence.

A global survey of host, aquatic, and soil microbiomes reveals shared abundance and genomic features between bacterial and fungal generalists

Environmental change, coupled with alteration in human lifestyles, is profoundly impacting the microbial communities critical to the health of the Earth and its inhabitants. To identify bacteria and fungi that are resistant and susceptible to habitat change, we analyze thousands of genera detected in 1,580 host, soil, and aquatic samples. This large-scale analysis identifies 48 bacterial and 4 fungal genera that are abundant across the three biomes, demonstrating fitness in diverse environmental conditions. Samples containing these generalists have significantly higher alpha diversity. These generalists play a significant role in shaping cross-kingdom community structure, boasting larger genomes with more secondary metabolism and antimicrobial resistance genes. Conversely, 30 bacterial and 19 fungal genera are only found in a single habitat, suggesting a limited ability to adapt to different and changing environments. These findings contribute to our understanding of microbial niche breadth and its consequences for global biodiversity loss.

Chaotic turnover of rare and abundant species in a strongly interacting model community

The composition of ecological communities varies not only between different locations but also in time. Understanding the fundamental processes that drive species toward rarity or abundance is crucial to assessing ecosystem resilience and adaptation to changing environmental conditions. In plankton communities in particular, large temporal fluctuations in species abundances have been associated with chaotic dynamics. On the other hand, microbial diversity is overwhelmingly sustained by a "rare biosphere" of species with very low abundances. We consider here the possibility that interactions within a species-rich community can relate both phenomena. We use a Lotka-Volterra model with weak immigration and strong, disordered, and mostly competitive interactions between hundreds of species to bridge single-species temporal fluctuations and abundance distribution patterns. We highlight a generic chaotic regime where a few species at a time achieve dominance but are continuously overturned by the invasion of formerly rare species. We derive a focal-species model that captures the intermittent boom-and-bust dynamics that every species undergoes. Although species cannot be treated as effectively uncorrelated in their abundances, the community's effect on a focal species can nonetheless be described by a time-correlated noise characterized by a few effective parameters that can be estimated from time series. The model predicts a nonunitary exponent of the power-law abundance decay, which varies weakly with ecological parameters, consistent with observation in marine protist communities. The chaotic turnover regime is thus poised to capture relevant ecological features of species-rich microbial communities.

Microfluidic approaches in microbial ecology

Microbial life is at the heart of many diverse environments and regulates most natural processes, from the functioning of animal organs to the cycling of global carbon. Yet, the study of microbial ecology is often limited by challenges in visualizing microbial processes and replicating the environmental conditions under which they unfold. Microfluidics operates at the characteristic scale at which microorganisms live and perform their functions, thus allowing for the observation and quantification of behaviors such as growth, motility, and responses to external cues, often with greater detail than classical techniques. By enabling a high degree of control in space and time of environmental conditions such as nutrient gradients, pH levels, and fluid flow patterns, microfluidics further provides the opportunity to study microbial processes in conditions that mimic the natural settings harboring microbial life. In this review, we describe how recent applications of microfluidic systems to microbial ecology have enriched our understanding of microbial life and microbial communities. We highlight discoveries enabled by microfluidic approaches ranging from single-cell behaviors to the functioning of multi-cellular communities, and we indicate potential future opportunities to use microfluidics to further advance our understanding of microbial processes and their implications.

Structured community transitions explain the switching capacity of microbial systems

Microbial systems appear to exhibit a relatively high switching capacity of moving back and forth among few dominant communities (taxon memberships). While this switching behavior has been mainly attributed to random environmental factors, it remains unclear the extent to which internal community dynamics affect the switching capacity of microbial systems. Here, we integrate ecological theory and empirical data to demonstrate that structured community transitions increase the dependency of future communities on the current taxon membership, enhancing the switching capacity of microbial systems. Following a structuralist approach, we propose that each community is feasible within a unique domain in environmental parameter space. Then, structured transitions between any two communities can happen with probability proportional to the size of their feasibility domains and inversely proportional to their distance in environmental parameter space-which can be treated as a special case of the gravity model. We detect two broad classes of systems with structured transitions: one class where switching capacity is high across a wide range of community sizes and another class where switching capacity is high only inside a narrow size range. We corroborate our theory using temporal data of gut and oral microbiota (belonging to class 1) as well as vaginal and ocean microbiota (belonging to class 2). These results reveal that the topology of feasibility domains in environmental parameter space is a relevant property to understand the changing behavior of microbial systems. This knowledge can be potentially used to understand the relevant community size at which internal dynamics can be operating in microbial systems.

Low-level resource partitioning supports coexistence among functionally redundant bacteria during successional dynamics

Members of microbial communities can substantially overlap in substrate use. However, what enables functionally redundant microorganisms to coassemble or even stably coexist remains poorly understood. Here, we show that during unstable successional dynamics on complex, natural organic matter, functionally redundant bacteria can coexist by partitioning low-concentration substrates even though they compete for one simple, dominant substrate. We allowed ocean microbial communities to self-assemble on leachates of the brown seaweed Fucus vesiculosus and then analyzed the competition among 10 taxonomically diverse isolates representing two distinct stages of the succession. All, but two isolates, exhibited an average of 90% ± 6% pairwise overlap in resource use, and functional redundancy of isolates from the same assembly stage was higher than that from between assembly stages, leading us to construct a simpler four-isolate community with two isolates from each of the early and late stages. We found that, although the short-term dynamics of the four-isolate communities in F. vesiculosus leachate was dependent on initial isolate ratios, in the long term, the four isolates stably coexist in F. vesiculosus leachate, albeit with some strains at low abundance. We therefore explored the potential for nonredundant substrate use by genomic content analysis and RNA expression patterns. This analysis revealed that the four isolates mainly differed in peripheral metabolic pathways, such as the ability to degrade pyrimidine, leucine, and tyrosine, as well as aromatic substrates. These results highlight the importance of fine-scale differences in metabolic strategies for supporting the frequently observed coexistence of large numbers of rare organisms in natural microbiomes.

Tailless and filamentous prophages are predominant in marine Vibrio

Although tailed bacteriophages (phages) of the class Caudoviricetes are thought to constitute the most abundant and ecologically relevant group of phages that can integrate their genome into the host chromosome, it is becoming increasingly clear that other prophages are widespread. Here, we show that prophages derived from filamentous and tailless phages with genome sizes below 16 kb make up the majority of prophages in marine bacteria of the genus Vibrio. To estimate prophage prevalence unaffected by database biases, we combined comparative genomics and chemical induction of 58 diverse Vibrio cyclitrophicus isolates, resulting in 107 well-curated prophages. Complemented with computationally predicted prophages, we obtained 1158 prophages from 931 naturally co-existing strains of the family Vibrionaceae. Prophages resembling tailless and filamentous phages predominated, accounting for 80% of all prophages in V. cyclitrophicus and 60% across the Vibrionaceae. In our experimental model, prophages of all three viral realms actively replicated upon induction indicating their ability to transfer to new hosts. Indeed, prophages were rapidly gained and lost, as suggested by variable prophage content between closely related V. cyclitrophicus. Prophages related to filamentous and tailless phages were integrated into only three genomic locations and restored the function of their integration site. Despite their small size, they contained highly diverse accessory genes that may contribute to host fitness, such as phage defense systems. We propose that, like their well-studied tailed equivalent, tailless and filamentous temperate phages are active and highly abundant drivers of host ecology and evolution in marine Vibrio, which have been largely overlooked.

Stochasticity-driven weekly fluctuations distinguished the temporal pattern of particle-associated microorganisms from its free-living counterparts in temperate coastal seawater

Microbial community dynamics directly determine their ecosystem functioning. Despite the well-known annual recurrence pattern, little is known how different lifestyles affect the temporal variation and how community assembly mechanisms change over different temporal scales. Here, through a high-resolution observation of size fractionated samples over 60 consecutive weeks, we investigate the distinction in weekly distribution pattern and assembly mechanism between free-living (FL) and particle-associated (PA) communities in highly dynamic coastal environments. A clear pattern of annual recurrence was observed, which was more pronounced in FL compared to PA, resulting in higher temporal specificity in the former samples. Both the two size fractions displayed significant temporal distance-decay patterns, yet the PA community showed a higher magnitude of community variation between adjacent weeks, likely caused by sudden, drastic and long-lived blooms of heterotrophic bacteria. Generally, determinism (environmental selection) had a greater effect on the community assembly than stochasticity (random birth, death, and dispersal events), with significant contributions from temperature and inorganic nutrients. However, a clear shift in the temporal assembly pattern was observed, transitioning from a prevalence of stochastic processes driving short-term (within a month) fluctuations to a dominance of deterministic processes over longer time intervals. Between adjacent weeks, stochasticity was more important in the community assembly of PA than FL. This study revealed that stochastic processes can lead to rapid, dramatic and irregular PA community fluctuations, indicating weak resistance and resilience to disturbances, which considering the role of PA microbes in carbon processing would significantly affect the coastal carbon cycle. Our results provided a new insight into the microbial community assembly mechanisms in the temporal dimension.

Vertical connectivity of microbiome and metabolome reveals depth-dependent variations across a deep cold-seep water column

Deciphering the vertical connectivity of oceanic microbiome and metabolome is crucial for understanding the carbon sequestration and achieving the carbon neutrality. However, we lack a systematic view of the interplay among particle transport, microbial community, and metabolic trait across depths. Through integrating the biogeochemical, microbial, and metabolic characteristics of a deep cold-seep water column (∼1989 m), we find the altered connectivity of microbial community and dissolved organic matter (DOM) across depths. Both the microbial communities (bacteria and protists) and DOM show a clear compositional connectivity from surface to the depth of 1000 m, highlighting the controls of sinking particle over microbial connectivity from the epipelagic to mesopelagic zone. However, due to the biological migration and ocean mixing, the fecal-associated bacteria and protistan consumers unexpectedly emerge and the degradation index of DOM substantially alters around 1000-1200 m. Collectively, we unveil the significance of multi-faceted particle dispersion, which supports the connectivity and variability of deep ocean microbial communities.

Turnover of strain-level diversity modulates functional traits in the honeybee gut microbiome between nurses and foragers

BACKGROUND

Strain-level diversity is widespread among bacterial species and can expand the functional potential of natural microbial communities. However, to what extent communities undergo consistent shifts in strain composition in response to environmental/host changes is less well understood.

RESULTS

Here, we used shotgun metagenomics to compare the gut microbiota of two behavioral states of the Western honeybee (Apis mellifera), namely nurse and forager bees. While their gut microbiota is composed of the same bacterial species, we detect consistent changes in strain-level composition between nurses and foragers. Single nucleotide variant profiles of predominant bacterial species cluster by behavioral state. Moreover, we identify strain-specific gene content related to nutrient utilization, vitamin biosynthesis, and cell-cell interactions specifically associated with the two behavioral states.

CONCLUSIONS

Our findings show that strain-level diversity in host-associated communities can undergo consistent changes in response to host behavioral changes modulating the functional potential of the community.

A Gulf Stream frontal eddy harbors a distinct microbiome compared to adjacent waters

Mesoscale oceanographic features, including eddies, have the potential to alter productivity and other biogeochemical rates in the ocean. Here, we examine the microbiome of a cyclonic, Gulf Stream frontal eddy, with a distinct origin and environmental parameters compared to surrounding waters, in order to better understand the processes dominating microbial community assembly in the dynamic coastal ocean. Our microbiome-based approach identified the eddy as distinct from the surround Gulf Stream waters. The eddy-associated microbial community occupied a larger area than identified by temperature and salinity alone, increasing the predicted extent of eddy-associated biogeochemical processes. While the eddy formed on the continental shelf, after two weeks both environmental parameters and microbiome composition of the eddy were most similar to the Gulf Stream, suggesting the effect of environmental filtering on community assembly or physical mixing with adjacent Gulf Stream waters. In spite of the potential for eddy-driven upwelling to introduce nutrients and stimulate primary production, eddy surface waters exhibit lower chlorophyll a along with a distinct and less even microbial community, compared to the Gulf Stream. At the population level, the eddy microbiome exhibited differences among the cyanobacteria (e.g. lower Trichodesmium and higher Prochlorococcus) and in the heterotrophic alpha Proteobacteria (e.g. lower relative abundances of specific SAR11 phylotypes) versus the Gulf Stream. However, better delineation of the relative roles of processes driving eddy community assembly will likely require following the eddy and surrounding waters since inception. Additionally, sampling throughout the water column could better clarify the contribution of these mesoscale features to primary production and carbon export in the oceans.

Circular and L50-like leaderless enterocins share a common ABC-transporter immunity gene

Microbes live within complex communities of interacting populations, either free-living in waters and soils or symbionts of animals and plants. Their interactions include the production of antimicrobial peptides (bacteriocins) to antagonize competitors, and these producers must carry their own immunity gene for self-protection. Whether other coexisting populations are sensitive or resistant to the bacteriocin producer will be key for the population dynamics within the microbial community. The immunity gene frequently consists of an ABC transporter to repel its own bacteriocin but rarely protects against a nonrelated bacteriocin. A case where this cross-resistance occurs mediated by a shared ABC transporter has been shown between enterocins MR10A/B and AS-48. The first is an L50-like leaderless enterocin, while AS-48 is a circular enterocin. In addition, L50-like enterocins such as MR10A/B have been found in E. faecalis and E. faecium, but AS-48 appears only in E. faecalis. Thus, using the ABC transporter of the enterocin MR10A/B gene cluster of Enterococcus faecalis MRR10-3 as a cross-resistance model, we aimed to unravel to what extent a particular ABC transporter can be shared across multiple bacteriocinogenic bacterial populations. To this end, we screened the MR10A/B-ABC transporters in available microbial genomes and analyzed their sequence homologies and distribution. Overall, our main findings are as follows: (i) the MR10A/B-ABC transporter is associated with multiple enterocin gene clusters; (ii) the different enterocins associated with this transporter have a saposin-like fold in common; (iii) the Mr10E component of the transporter is more conserved within its associated enterocin, while the Mr10FGH components are more conserved within the carrying species. This is the least known component of the transporter, but it has shown the greatest specificity to its corresponding enterocin. Bacteriocins are now being investigated as an alternative to antibiotics; hence, the wider or narrower distribution of the particular immunity gene should be taken into account for clinical applications to avoid the selection of resistant strains. Further research will be needed to investigate the mechanistic interactions between the Mr10E transporter component and the bacteriocin as well as the specific ecological and evolutionary mechanisms involved in the spread of the immunity transporter across multiple bacteriocins.

Epilithic Bacterial Assemblages on Subtidal Rocky Reefs: Variation Among Alternative Habitats at Ambient and Enhanced Nutrient Levels

Temperate rocky reefs often support mosaics of alternative habitats such as macroalgal forests, algal turfs and sea urchin barrens. Although the composition of epilithic microbial biofilms (EMBs) is recognized as a major determinant of macroalgal recruitment, their role in regulating the stability of alternative habitats on temperate rocky reefs remains unexplored. On shallow rocky reefs of the Island of Capraia (NW Mediterranean), we compared EMB structure among canopy stands formed by the fucoid Ericaria brachycarpa, algal turfs, and urchin barrens under ambient versus experimentally enhanced nutrient levels. The three habitats shared a core microbial community consisting of 21.6 and 25.3% of total ASVs under ambient and enhanced nutrient conditions, respectively. Although Gammaproteobacteria, Alphaproteobacteria and Flavobacteriia were the most abundant classes across habitats, multivariate analyses at the ASV level showed marked differences in EMB composition among habitats. Enhancing nutrient level had no significant effect on EMBs, although it increased their similarity between macroalgal canopy and turf habitats. At both ambient and enriched nutrient levels, ASVs mostly belonging to Proteobacteria and Bacteroidetes were more abundant in EMBs from macroalgal canopies than barrens. In contrast, ASVs belonging to the phylum of Proteobacteria and, in particular, to the families of Rhodobacteraceae and Flavobacteriaceae at ambient nutrient levels and of Rhodobacteraceae and Bacteriovoracaceae at enhanced nutrient levels were more abundant in turf than canopy habitats. Our results show that primary surfaces from alternative habitats that form mosaics on shallow rocky reefs in oligotrophic areas host distinct microbial communities that are, to some extent, resistant to moderate nutrient enhancement. Understanding the role of EMBs in generating reinforcing feedback under different nutrient loading regimes appears crucial to advance our understanding of the mechanisms underpinning the stability of habitats alternative to macroalgal forests as well as their role in regulating reverse shifts.

The Reemergence of Phycopathology: When Algal Biology Meets Ecology and Biosecurity

Viruses, bacteria, and eukaryotic symbionts interact with algae in a variety of ways to cause disease complexes, often shaping marine and freshwater ecosystems. The advent of phyconomy (a.k.a. seaweed agronomy) represents a need for a greater understanding of algal disease interactions, where underestimated cryptic diversity and lack of phycopathological basis are prospective constraints for algal domestication. Here, we highlight the limited yet increasing knowledge of algal pathogen biodiversity and the ecological interaction with their algal hosts. Finally, we discuss how ecology and cultivation experience contribute to and reinforce aquaculture practice, with the potential to reshape biosecurity policies of seaweed cultivation worldwide.

Multiyear analysis uncovers coordinated seasonality in stocks and composition of the planktonic food web in the Baltic Sea proper

The planktonic realm from bacteria to zooplankton provides the baseline for pelagic aquatic food webs. However, multiple trophic levels are seldomly included in time series studies, hampering a holistic understanding of the influence of seasonal dynamics and species interactions on food web structure and biogeochemical cycles. Here, we investigated plankton community composition, focusing on bacterio-, phyto- and large mesozooplankton, and how biotic and abiotic factors correlate at the Linnaeus Microbial Observatory (LMO) station in the Baltic Sea from 2011 to 2018. Plankton communities structures showed pronounced dynamic shifts with recurring patterns. Summarizing the parts of the planktonic microbial food web studied here to total carbon, a picture emerges with phytoplankton consistently contributing > 39% while bacterio- and large mesozooplankton contributed ~ 30% and ~ 7%, respectively, during summer. Cyanophyceae, Actinobacteria, Bacteroidetes, and Proteobacteria were important groups among the prokaryotes. Importantly, Dinophyceae, and not Bacillariophyceae, dominated the autotrophic spring bloom whereas Litostomatea (ciliates) and Appendicularia contributed significantly to the consumer entities together with the more traditionally observed mesozooplankton, Copepoda and Cladocera. Our findings of seasonality in both plankton composition and carbon stocks emphasize the importance of time series analyses of food web structure for characterizing the regulation of biogeochemical cycles and appropriately constraining ecosystem models.

Photobacterium predominate the microbial communities of muscle of European plaice (Pleuronectes platessa) caught in the Norwegian sea independent of skin and gills microbiota, fishing season, and storage conditions

The aim of this study was to investigate seasonal variations (September, December and April) in the initial microbial communities of skin and gills' external mucosal tissues (EMT) and muscle of European plaice (Pleuronectes platessa). Moreover, a potential relationship between EMT and fresh muscle microbiota was examined. The microbial community succession in plaice muscle as a function of fishing season and storage conditions was also investigated. The selected seasons for the storage experiment were September and April. Investigated storage conditions were; fillets packaged in either vacuum or modified atmosphere (70 % CO2, 20 % N2, 10 % O2) and chilled/refrigerated conditions (4 °C). Whole fish stored on ice (0 °C) was selected as a commercial standard. Seasonal variations were detected in the initial microbial communities of EMT and plaice muscle. The highest microbial diversity was found in EMT and muscle of April-caught plaice, followed by December and September catch indicating the important role of environmental factors in shaping the initial EMT and muscle microbial communities. The EMT microbial communities were more diverse than fresh muscle samples. The low number of shared taxa between EMT and initial muscle microbial communities indicates that only a minor part of the muscle microbiota came from the EMT. Psychrobacter and Photobacterium were the predominant genera in the EMT microbial communities in all seasons. Photobacterium dominated the initial muscle microbial communities with a gradual seasonal reduction of its abundance from September to April. Storage time and storage conditions shaped a less diverse and distinct community compared to the fresh muscle. However, no clear separation was seen between the communities at the middle and end of storage time. Regardless of EMT microbiota, fishing season and storage conditions, Photobacterium dominated the microbial communities of stored muscle samples. The Photobacterium prevalence as the primary specific spoilage organism (SSO) could be attributed to its high relative abundance in the initial microbiota of muscle and its CO2-tolerance. The findings of this study indicate the important contribution of Photobacterium to the microbial spoilage of plaice. Thus, the development of innovative preservation techniques addressing the rapid growth of Photobacterium could contribute to the production of high-quality and shelf-stable convenient retail plaice products.

Wavelet clustering analysis as a tool for characterizing community structure in the human microbiome

Human microbiome research is helped by the characterization of microbial networks, as these may reveal key microbes that can be targeted for beneficial health effects. Prevailing methods of microbial network characterization are based on measures of association, often applied to limited sampling points in time. Here, we demonstrate the potential of wavelet clustering, a technique that clusters time series based on similarities in their spectral characteristics. We illustrate this technique with synthetic time series and apply wavelet clustering to densely sampled human gut microbiome time series. We compare our results with hierarchical clustering based on temporal correlations in abundance, within and across individuals, and show that the cluster trees obtained by using either method are significantly different in terms of elements clustered together, branching structure and total branch length. By capitalizing on the dynamic nature of the human microbiome, wavelet clustering reveals community structures that remain obscured in correlation-based methods.

Aging by Near-Extinctions in Many-Variable Interacting Populations

Models of many-species ecosystems, such as the Lotka-Volterra and replicator equations, suggest that these systems generically exhibit near-extinction processes, where population sizes go very close to zero for some time before rebounding, accompanied by a slowdown of the dynamics (aging). Here, we investigate the connection between near-extinction and aging by introducing an exactly solvable many-variable model, where the time derivative of each population size vanishes at both zero and some finite maximal size. We show that aging emerges generically when random interactions are taken between populations. Population sizes remain exponentially close (in time) to the absorbing values for extended periods of time, with rapid transitions between these two values. The mechanism for aging is different from the one at play in usual glassy systems: At long times, the system evolves in the vicinity of unstable fixed points rather than marginal ones.

Monitoring trophic status using in situ data and Sentinel-2 MSI algorithm: lesson from Lake Malombe, Malawi

With excessive nutrient enrichment exacerbated by anthropogenic drivers, many standing water bodies are changing from oligotrophic to mesotrophic, eutrophic, and finally hypertrophic-negatively affecting ecosystem functioning, biodiversity, and human populations. Efforts have been devoted to developing novel algorithms for estimating chlorophyll-a (chl-a), cyno-blooms, and floating vegetation. However, to this date, little research has focused on freshwater lakes in the data-scarce Sub-Saharan African countries such as Malawi. We, therefore, estimated the trophic status of Lake Malombe in Malawi-a lake likely to be affected by eutrophication and algal bloom-emerging threats to freshwater ecosystem functioning globally-especially with the onset of climatic and anthropogenic drivers. We integrated in situ data with high-resolution Sentinel-2 Multispectral Imagery Analysis (MSI). We independently assessed the remote sensing technique using in situ data and tested the model at multiple stages. The scatter plot showed that most points were in the 95% confidence interval. The validation results between the measured in situ chl-a concentrations and the Sentinel-2 MSI-based chl-a retrieval had a root mean square error (RMSE) of 2.88 µg/L. The chl-a concentrations retrieved from MSI images were consistent with in situ data, indicating that the normalized difference chlorophyll index (NDCI) algorithm estimated chl-a concentrations in Lake Malombe with acceptable accuracy. Dissolved oxygen (DO), sulfate (SO₄^2-), nitrite [Formula: see text], soluble reactive phosphorous [Formula: see text]), total dissolved solids (TDS), and chl-a, except for temperatures from the hot-dry-season, cold-dry-windy-season, and rainy-season, were significantly different (P < 0.05). The Sentinel-2 MSI imagery analysis also depicted similar results, with high chl-a concentration reported in March (rainy season) and October (hot-dry season) and the lowest from May to August (cold-dry-windy season). On the contrary, the ANOVA results for water quality parameters from all five points had P > 0.05. The correlation matrix showed coefficients of (0.798 < r < 0.930, n = 30, P < 0.005), suggesting that Lake Malombe is homogenous. Our results demonstrate that integrating remote sensing based on MSI imagery and in situ data to estimate chl-a can provide an effective tool for monitoring eutrophication in small, medium, and large standing waterbodies-crucial information required to respond to global ecological and climatic dynamics.

Planktonic ciliate community driven by environmental variables and cyanobacterial blooms: A 9-year study in two subtropical reservoirs

It is well-established that environmental variability and cyanobacterial blooms have major effects on the assembly and functioning of bacterial communities in both marine and freshwater habitats. It remains unclear, however, how the ciliate community responds to such changes over the long-term, particularly in subtropical lake and reservoir ecosystems. We analysed 9-year planktonic ciliate data series from the surface water of two subtropical reservoirs to elucidate the role of cyanobacterial bloom and environmental variabilities on the ciliate temporal dynamics. We identified five distinct periods of cyanobacterial succession in both reservoirs. Using multiple time-scale analyses, we found that the interannual variability of ciliate communities was more strongly related to cyanobacterial blooms than to other environmental variables or to seasonality. Moreover, the percentage of species turnover across cyanobacterial bloom and non-bloom periods increased significantly with time over the 9-year period. Phylogenetic analyses further indicated that 84 %-86 % of ciliate community turnover was governed by stochastic dispersal limitation or undominated processes, suggesting that the ciliate communities in subtropical reservoirs were mainly controlled by neutral processes. However, short-term blooms increased the selection pressure and drove 30 %-53 % of the ciliate community turnover. We found that the ciliate community composition was influenced by environmental conditions with nutrients, cyanobacterial biomass and microzooplankton having direct and/or indirect significant effects on the ciliate taxonomic or functional community dynamics. Our results provide new insights into the long-term temporal dynamics of planktonic ciliate communities under cyanobacterial bloom disturbance.

Microbial Community Dynamics during a Harmful Chrysochromulina leadbeateri Bloom in Northern Norway

A harmful algal bloom occurred in late spring 2019 across multiple, interconnected fjords and bays in northern Norway. The event was caused by the haptophyte Chrysochromulina leadbeateri and led to severe fish mortality at several salmon aquaculture facilities. This study reports on the spatial and temporal succession dynamics of the holistic marine microbiome associated with this bloom by relating all detectable 18S and 16S rRNA gene amplicon sequence variants to the relative abundance of the C. leadbeateri focal taxon. A k-medoid clustering enabled inferences on how the causative focal taxon cobloomed with diverse groups of bacteria and microeukaryotes. These coblooming patterns showed high temporal variability and were distinct between two geographically separated time series stations during the regional harmful algal bloom. The distinct blooming patterns observed with respect to each station were poorly connected to environmental conditions, suggesting that other factors, such as biological interactions, may be at least as important in shaping the dynamics of this type of harmful algal bloom. A deeper understanding of microbiome succession patterns during these rare but destructive events will help guide future efforts to forecast deviations from the natural bloom cycles of the northern Norwegian coastal marine ecosystems that are home to intensive aquaculture activities. IMPORTANCE The 2019 Chrysochromulina leadbeateri bloom in northern Norway had a major impact on the local economy and society through its devastating effect on the aquaculture industry. However, many fail to remember that C. leadbeateri is, in fact, a common member of the seasonal marine microbiome and the same spring phytoplankton blooms that support the marine ecosystem. It is challenging to draw any conclusions about exact causation behind the harmful bloom of 2019, especially since the natural bloom cycles of C. leadbeateri are not well understood. This study begins to fill major knowledge gaps that may lead to future forecasting abilities, by providing a molecular-based investigation of the destructive 2019 bloom that presents new insights into a seasonal marine microbial ecosystem during one of these sporadically reoccurring events.

Fine-Scale Structuring of Planktonic Vibrio spp. in the Chinese Marginal Seas

Vibrio is ubiquitous in marine environments with high metabolism flexibility and genome plasticity. Studies have investigated the ecological distribution of Vibrio spp. in several narrow zones, but a broad scale pattern of distribution and community assembly is still lacking. Here, we elucidated the distribution of Vibrio spp. in seawater along the Chinese marginal seas with a high spatial range. Comparison of Vibrio abundance between 3- and 0.2-μm-pore-size membranes showed distinction in preferential lifestyle. Vibrio spp. in the Yellow Sea (YS) was low in abundance and adopted a particle-associated lifestyle, whereas that in the East China Sea (ECS) and South China Sea (SCS) was more abundant and was likely in a temporary free-living state as a strategy to cope with nutrient limitation. Vibrio community compositions were also separated by sampling area, with different dominant groups in YS (Vibrio chagasii and Vibrio harveyi), ECS and SCS (Vibrio japonicus and V. chagasii). The community niche breadth was significantly wider in ECS and SCS than that of YS. Among species, V. chagasii and V. harveyi had the largest niche breadths likely reflecting strong competitive positions. Stochastic processes played important roles in shaping the geographical pattern of the vibrionic community. Environmental selection (e.g., temperature, salinity, and dissolved oxygen) had a much greater impact on the community in surface than in bottom water. The large proportions of unexplained variations (78.9%) imply complex mechanisms in their community assembly. Our study provides insights into the spatial distribution patterns and underlying assembly mechanisms of Vibrio at a broad spatial scale. IMPORTANCE Vibrio spp. may exert large impacts on biogeochemical cycling in coastal habitats, and their ecological importance has drawn increasing attention. Here, we investigated the spatial distribution pattern and community assembly of Vibrio populations along the Chinese marginal seas, spanning a wide spatial scale. Our results showed that the abundances of the Vibrio population increased with decreasing latitude and their preferential lifestyle differed among adjacent coastal areas. The compositions of Vibrio spp. were also separated by geographical location, which was mainly attributable to stochastic processes. Overall, this work contributes to the understanding of the ecological distribution patterns and the community assembly mechanisms of marine vibrios at a high spatial range. The large proportion of unexplained variations indicates the existence of complex mechanisms in the assembly of vibrionic community which should be considered comprehensively in future.

Shared and contrasting associations in the dynamic nano- and picoplankton communities of two close but contrasting sites from the Bay of Biscay

Pico- and nanoplankton are key players in the marine ecosystems due to their implication in the biogeochemical cycles, nutrient recycling and the pelagic food webs. However, the specific dynamics and niches of most bacterial, archaeal and eukaryotic plankton remain unknown, as well as the interactions between them. Better characterization of these is critical for understanding and predicting ecosystem functioning under anthropogenic pressures. We used environmental DNA metabarcoding across a 6-year time series to explore the structure and seasonality of pico- and nanoplankton communities in two sites of the Bay of Biscay, one coastal and one offshore, and construct association networks to reveal potential keystone and connector taxa. Temporal trends in alpha diversity were similar between the two sites, and concurrent communities more similar than within the same site at different times. However, we found differences between the network topologies of the two sites, with both shared and site-specific keystones and connectors. For example, Micromonas, with lower abundance in the offshore site is a keystone here, indicating a stronger effect of associations such as resource competition. This study provides an example of how time series and association network analysis can reveal how similar communities may function differently despite being geographically close.

The Ecological Differentiation of Particle-Attached and Free-Living Bacterial Communities in a Seasonal Flooding Lake-the Poyang Lake

Particle-attached (PA) and free-living (FL) bacterial communities play essential roles in the biogeochemical cycling of essential nutrients in aquatic environments. However, little is known about the factors that drive the differentiation of bacterial lifestyles, especially in flooding lake systems. Here we assessed the compositional and functional similarities between the FL and PA bacterial fractions in a typical flooding lake-the Poyang Lake (PYL) of China. The results revealed that PA communities had significantly different compositions and functions from FL communities in every hydrological period, and the diversity of both PA and FL communities was affected mainly by the water regime rather than bacterial lifestyles. PA communities were more diverse and enriched with Proteobacteria and Bacteroidetes, while FL communities had more Actinobacteria. There was a higher abundance of photosynthetic and nitrogen-cycling bacterial groups in PA communities, but a higher abundance of members involved in hydrocarbon degradation, aromatic hydrocarbon degradation, and methylotrophy in FL communities. Water properties (e.g., temperature, pH, total phosphorus) significantly regulated the lifestyle variations of PA and FL bacteria in PYL. Collectively, our results have demonstrated a clear ecological differentiation of PA and FL bacterial communities in flooding lakes, suggesting that the connectivity between FL and PA bacterial fractions is water property-related rather than water regime-related.

Estuarine microbial networks and relationships vary between environmentally distinct communities

Microbial interactions have profound impacts on biodiversity, biogeochemistry, and ecosystem functioning, and yet, they remain poorly understood in the ocean and with respect to changing environmental conditions. We applied hierarchical clustering of an annual 16S and 18S amplicon dataset in the Skidaway River Estuary, which revealed two similar clusters for prokaryotes (Bacteria and Archaea) and protists: Cluster 1 (March-May and November-February) and Cluster 2 (June-October). We constructed co-occurrence networks from each cluster to explore how microbial networks and relationships vary between environmentally distinct periods in the estuary. Cluster 1 communities were exposed to significantly lower temperature, sunlight, NO3, and SiO4; only NH4 was higher at this time. Several network properties (e.g., edge number, degree, and centrality) were elevated for networks constructed with Cluster 1 vs. 2 samples. There was also evidence that microbial nodes in Cluster 1 were more connected (e.g., higher edge density and lower path length) compared to Cluster 2, though opposite trends were observed when networks considered Prokaryote-Protist edges only. The number of Prokaryote-Prokaryote and Prokaryote-Protist edges increased by >100% in the Cluster 1 network, mainly involving Flavobacteriales, Rhodobacterales, Peridiniales, and Cryptomonadales associated with each other and other microbial groups (e.g., SAR11, Bacillariophyta, and Strombidiida). Several Protist-Protist associations, including Bacillariophyta correlated with Syndiniales (Dino-Groups I and II) and an Unassigned Dinophyceae group, were more prevalent in Cluster 2. Based on the type and sign of associations that increased in Cluster 1, our findings indicate that mutualistic, competitive, or predatory relationships may have been more representative among microbes when conditions were less favorable in the estuary; however, such relationships require further exploration and validation in the field and lab. Coastal networks may also be driven by shifts in the abundance of certain taxonomic or functional groups. Sustained monitoring of microbial communities over environmental gradients, both spatial and temporal, is critical to predict microbial dynamics and biogeochemistry in future marine ecosystems.

Bacterial growth in multicellular aggregates leads to the emergence of complex life cycles

Facultative multicellular behaviors expand the metabolic capacity and physiological resilience of bacteria. Despite their ubiquity in nature, we lack an understanding of how these behaviors emerge from cellular-scale phenomena. Here, we show how the coupling between growth and resource gradient formation leads to the emergence of multicellular lifecycles in a marine bacterium. Under otherwise carbon-limited growth conditions, Vibrio splendidus 12B01 forms clonal multicellular groups to collectively harvest carbon from soluble polymers of the brown-algal polysaccharide alginate. As they grow, groups phenotypically differentiate into two spatially distinct sub-populations: a static "shell" surrounding a motile, carbon-storing "core." Differentiation of these two sub-populations coincides with the formation of a gradient in nitrogen-source availability within clusters. Additionally, we find that populations of cells containing a high proportion of carbon-storing individuals propagate and form new clusters more readily on alginate than do populations with few carbon-storing cells. Together, these results suggest that local metabolic activity and differential partitioning of resources leads to the emergence of reproductive cycles in a facultatively multicellular bacterium.

Diversity and Selection of Surface Marine Microbiomes in the Atlantic-Influenced Arctic

Arctic marine environments are experiencing rapid changes due to the polar amplification of global warming. These changes impact the habitat of the cold-adapted microbial communities, which underpin biogeochemical cycles and marine food webs. We comparatively investigated the differences in prokaryotic and microeukaryotic taxa between summer surface water microbiomes sampled along a latitudinal transect from the ice-free southern Barents Sea and into the sea-ice-covered Nansen Basin to disentangle the dominating community (ecological) selection processes driving phylogenetic diversity. The community structure and richness of each site-specific microbiome were assessed in relation to the physical and biogeochemical conditions of the environment. A strong homogeneous deterministic selection process was inferred across the entire sampling transect via a phylogenetic null modeling approach. The microbial species richness and diversity were not negatively influenced by northward decreasing temperature and salinity. The results also suggest that regional phytoplankton blooms are a major prevalent factor in governing the bacterial community structure. This study supports the consideration that strong homogeneous selection is imposed across these cold-water marine environments uniformly, regardless of geographic assignments within either the Nansen Basin or the Barents Sea.

Niche Partitioning of Labyrinthulomycete Protists Across Sharp Coastal Gradients and Their Putative Relationships With Bacteria and Fungi

While planktonic microbes play key roles in the coastal oceans, our understanding of heterotrophic microeukaryotes’ ecology, particularly their spatiotemporal patterns, drivers, and functions, remains incomplete. In this study, we focus on a ubiquitous marine fungus-like protistan group, the Labyrinthulomycetes, whose biomass can exceed that of bacterioplankton in coastal oceans but whose ecology is largely unknown. Using quantitative PCR and amplicon sequencing of their 18S rRNA genes, we examine their community variation in repeated five-station transects across the nearshore-to-offshore surface waters of North Carolina, United States. Their total 18S rRNA gene abundance and phylotype richness decrease significantly from the resource-rich nearshore to the oligotrophic offshore waters, but their Pielou’s community evenness appears to increase offshore. Similar to the bacteria and fungi, the Labyrinthulomycete communities are significantly structured by distance from shore, water temperature, and other environmental factors, suggesting potential niche partitioning. Nevertheless, only several Labyrinthulomycete phylotypes, which belong to aplanochytrids, thraustochytrids, or unclassified Labyrinthulomycetes, are prevalent and correlated with cohesive bacterial communities, while more phylotypes are patchy and often co-occur with fungi. Overall, these results complement previous time-series observations that resolve the Labyrinthulomycetes as persistent and short-blooming ecotypes with distinct seasonal preferences, further revealing their partitioning spatial patterns and multifaceted roles in coastal marine microbial food webs.

High-Frequency Variability of Bacterioplankton in Response to Environmental Drivers in Red Sea Coastal Waters

Autotrophic and heterotrophic bacterioplankton are essential to the biogeochemistry of tropical ecosystems. However, the processes that govern their dynamics are not well known. We provide here a high-frequency assessment of bacterial community dynamics and concurrent environmental factors in Red Sea coastal waters. Weekly sampling of surface samples during a full annual cycle at an enclosed station revealed high variability in ecological conditions, which reflected in changes of major bacterioplankton communities. Temperature varied between 23 and 34°C during the sampling period. Autotrophic (Synechococcus, 1.7–16.2 × 10⁴ cells mL⁻¹) and heterotrophic bacteria (1.6–4.3 × 10⁵ cells mL⁻¹) showed two maxima in abundance in spring and summer, while minima were found in winter and autumn. Heterotrophic cells with high nucleic acid content (HNA) peaked in July, but their contribution to the total cell counts (35–60%) did not show a clear seasonal pattern. Actively respiring cells (CTC+) contributed between 4 and 51% of the total number of heterotrophic bacteria, while live cells (with intact membrane) consistently accounted for over 90%. Sequenced 16S rRNA amplicons revealed a predominance of Proteobacteria in summer and autumn (>40%) and a smaller contribution in winter (21–24%), with members of the Alphaproteobacteria class dominating throughout the year. The contribution of the Flavobacteriaceae family was highest in winter (21%), while the Rhodobacteraceae contribution was lowest (6%). Temperature, chlorophyll-a, and dissolved organic carbon concentration were the environmental variables with the greatest effects on bacterial abundance and diversity patterns.

A comparative whole-genome approach identifies bacterial traits for marine microbial interactions

Microbial interactions shape the structure and function of microbial communities with profound consequences for biogeochemical cycles and ecosystem health. Yet, most interaction mechanisms are studied only in model systems and their prevalence is unknown. To systematically explore the functional and interaction potential of sequenced marine bacteria, we developed a trait-based approach, and applied it to 473 complete genomes (248 genera), representing a substantial fraction of marine microbial communities. We identified genome functional clusters (GFCs) which group bacterial taxa with common ecology and life history. Most GFCs revealed unique combinations of interaction traits, including the production of siderophores (10% of genomes), phytohormones (3-8%) and different B vitamins (57-70%). Specific GFCs, comprising Alpha- and Gammaproteobacteria, displayed more interaction traits than expected by chance, and are thus predicted to preferentially interact synergistically and/or antagonistically with bacteria and phytoplankton. Linked trait clusters (LTCs) identify traits that may have evolved to act together (e.g., secretion systems, nitrogen metabolism regulation and B vitamin transporters), providing testable hypotheses for complex mechanisms of microbial interactions. Our approach translates multidimensional genomic information into an atlas of marine bacteria and their putative functions, relevant for understanding the fundamental rules that govern community assembly and dynamics.

NetGAM: Using generalized additive models to improve the predictive power of ecological network analyses constructed using time-series data

Ecological network analyses are used to identify potential biotic interactions between microorganisms from species abundance data. These analyses are often carried out using time-series data; however, time-series networks have unique statistical challenges. Time-dependent species abundance data can lead to species co-occurrence patterns that are not a result of direct, biotic associations and may therefore result in inaccurate network predictions. Here, we describe a generalize additive model (GAM)-based data transformation that removes time-series signals from species abundance data prior to running network analyses. Validation of the transformation was carried out by generating mock, time-series datasets, with an underlying covariance structure, running network analyses on these datasets with and without our GAM transformation, and comparing the network outputs to the known covariance structure of the simulated data. The results revealed that seasonal abundance patterns substantially decreased the accuracy of the inferred networks. In addition, the GAM transformation increased the predictive power (F1 score) of inferred ecological networks on average and improved the ability of network inference methods to capture important features of network structure. This study underscores the importance of considering temporal features when carrying out network analyses and describes a simple, effective tool that can be used to improve results.

Diatom Biogeography, Temporal Dynamics, and Links to Bacterioplankton across Seven Oceanographic Time-Series Sites Spanning the Australian Continent

Diatom communities significantly influence ocean primary productivity and carbon cycling, but their spatial and temporal dynamics are highly heterogeneous and are governed by a complex diverse suite of abiotic and biotic factors. We examined the seasonal and biogeographical dynamics of diatom communities in Australian coastal waters using amplicon sequencing data (18S-16S rRNA gene) derived from a network of oceanographic time-series spanning the Australian continent. We demonstrate that diatom community composition in this region displays significant biogeography, with each site harbouring distinct community structures. Temperature and nutrients were identified as the key environmental contributors to differences in diatom communities at all sites, collectively explaining 21% of the variability observed in diatoms assemblages. However, specific groups of bacteria previously implicated in mutualistic ecological interactions with diatoms (Rhodobacteraceae, Flavobacteriaceae and Alteromonadaceae) also explained a further 4% of the spatial dynamics observed in diatom community structure. We also demonstrate that the two most temperate sites (Port Hacking and Maria Island) exhibited strong seasonality in diatom community and that at these sites, winter diatom communities co-occurred with higher proportion of Alteromonadaceae. In addition, we identified significant co-occurrence between specific diatom and bacterial amplicon sequence variants (ASVs), with members of the Roseobacter and Flavobacteria clades strongly correlated with some of the most abundant diatom genera (Skeletonema, Thalassiosira, and Cylindrotheca). We propose that some of these co-occurrences might be indicative of ecologically important interactions between diatoms and bacteria. Our analyses reveal that in addition to physico-chemical conditions (i.e., temperature, nutrients), the relative abundance of specific groups of bacteria appear to play an important role in shaping the spatial and temporal dynamics of marine diatom communities.

Resolving the structure of phage-bacteria interactions in the context of natural diversity

Microbial communities are shaped by viral predators. Yet, resolving which viruses (phages) and bacteria are interacting is a major challenge in the context of natural levels of microbial diversity. Thus, fundamental features of how phage-bacteria interactions are structured and evolve in the wild remain poorly resolved. Here we use large-scale isolation of environmental marine Vibrio bacteria and their phages to obtain estimates of strain-level phage predator loads, and use all-by-all host range assays to discover how phage and host genomic diversity shape interactions. We show that lytic interactions in environmental interaction networks (as observed in agar overlay) are sparse-with phage predator loads being low for most bacterial strains, and phages being host-strain-specific. Paradoxically, we also find that although overlap in killing is generally rare between tailed phages, recombination is common. Together, these results suggest that recombination during cryptic co-infections is an important mode of phage evolution in microbial communities. In the development of phages for bioengineering and therapeutics it is important to consider that nucleic acids of introduced phages may spread into local phage populations through recombination, and that the likelihood of transfer is not predictable based on lytic host range.

Seasonal niche differentiation among closely related marine bacteria

Bacteria display dynamic abundance fluctuations over time in marine environments, where they play key biogeochemical roles. Here, we characterized the seasonal dynamics of marine bacteria in a coastal oligotrophic time series station, tested how similar the temporal niche of closely related taxa is, and what are the environmental parameters modulating their seasonal abundance patterns. We further explored how conserved the niche is at higher taxonomic levels. The community presented recurrent patterns of seasonality for 297 out of 6825 amplicon sequence variants (ASVs), which constituted almost half of the total relative abundance (47%). For certain genera, niche similarity decreased as nucleotide divergence in the 16S rRNA gene increased, a pattern compatible with the selection of similar taxa through environmental filtering. Additionally, we observed evidence of seasonal differentiation within various genera as seen by the distinct seasonal patterns of closely related taxa. At broader taxonomic levels, coherent seasonal trends did not exist at the class level, while the order and family ranks depended on the patterns that existed at the genus level. This study identifies the coexistence of closely related taxa for some bacterial groups and seasonal differentiation for others in a coastal marine environment subjected to a strong seasonality.

Year-round dynamics of amplicon sequence variant communities differ among eukaryotes, Imitervirales, and prokaryotes in a coastal ecosystem

Coastal microbial communities are affected by seasonal environmental change, biotic interactions, and fluctuating nutrient availability. We investigated the seasonal dynamics of communities of eukaryotes, a major group of double-stranded DNA viruses that infect eukaryotes (order Imitervirales; phylum Nucleocytoviricota), and prokaryotes in the Uranouchi Inlet, Kochi, Japan. We performed metabarcoding using ribosomal RNA genes and viral polB genes as markers in 43 seawater samples collected over 20 months. Eukaryotes, prokaryotes, and Imitervirales communities characterized by the compositions of amplicon sequence variants (ASVs) showed synchronic seasonal cycles. However, the community dynamics showed intriguing differences in several aspects, such as the recovery rate after a year. We also showed that the differences in community dynamics were at least partially explained by differences in recurrence/persistence levels of individual ASVs among eukaryotes, prokaryotes, and Imitervirales. Prokaryotic ASVs were the most persistent, followed by eukaryotic ASVs and Imitervirales ASVs, which were the least persistent. We argue that the differences in the specificity of interactions (virus-eukaryote vs. prokaryote-eukaryote) as well as the niche breadth of community members were at the origin of the distinct community dynamics among eukaryotes, their viruses, and prokaryotes.

Composition change and decreased diversity of microbial eukaryotes in the coastal upwelling waters of South China Sea

Upwelling plays an important role in marine ecosystems and potentially reshapes microbial communities by enhanced dispersal and distinct environmental drivers. Relative to that of bacterioplankton, however, the response of eukaryotic microbes to upwelling is largely unknown. Here, we investigated the influence of coastal upwelling in South China Sea on the microbial eukaryotic communities. Unlike several folds of increase in the cell abundance of bacterioplankton in upwelling than non-upwelling stations at corresponding water layers, no significant difference was detected for the total microbial eukaryotic 18S rRNA gene abundance. Moreover, the microbial eukaryotes in the upwelling stations exhibited increasing 18S rRNA gene abundance from the surface to the deep, contrasting the vertical cell abundance pattern of the bacterioplankton; but their vertical abundance patterns were similar in non-upwelling stations. Importantly, the coastal upwelling significantly reduced the community evenness of the microbial eukaryotes and slightly reduced their Shannon diversity. Their community composition also varied obviously especially between the surface waters of upwelling and non-upwelling stations. Among the dominant supergroups, Alveolata was found to be less abundant while Stramenopiles, particularly thraustochytrids and diatoms, to be more abundant in the surface water of upwelling than non-upwelling stations. Temperature was identified as the most important factor of the microbial eukaryotic community composition, suggesting potential effects of the cold upwelling water masses on specific taxa. Overall, our results reveal significant and distinct impacts of coastal upwelling on the abundance, diversity, and community structure of microbial eukaryotes, filling the knowledge gap about the microbial responses to this important marine phenomenon.

Rapid evolutionary turnover of mobile genetic elements drives bacterial resistance to phages

[Figure: see text].

Heavy-tailed abundance distributions from stochastic Lotka-Volterra models

Microbial communities found in nature are composed of many rare species and few abundant ones, as reflected by their heavy-tailed abundance distributions. How a large number of species can coexist in those complex communities and why they are dominated by rare species is still not fully understood. We show how heavy-tailed distributions arise as an emergent property from large communities with many interacting species in population-level models. To do so, we rely on generalized Lotka-Volterra models for which we introduce a global maximal capacity. This maximal capacity accounts for the fact that communities are limited by available resources and space. In a parallel ad hoc approach, we obtain heavy-tailed abundance distributions from logistic models, without interactions, through specific distributions of the parameters. We expect both mechanisms, interactions between many species and specific parameter distributions, to be relevant to explain the observed heavy tails.

Successional dynamics and alternative stable states in a saline activated sludge microbial community over 9 years

BACKGROUND

Microbial communities in both natural and applied settings reliably carry out myriads of functions, yet how stable these taxonomically diverse assemblages can be and what causes them to transition between states remains poorly understood. We studied monthly activated sludge (AS) samples collected over 9 years from a full-scale wastewater treatment plant to answer how complex AS communities evolve in the long term and how the community functions change when there is a disturbance in operational parameters.

RESULTS

Here, we show that a microbial community in activated sludge (AS) system fluctuated around a stable average for 3 years but was then abruptly pushed into an alternative stable state by a simple transient disturbance (bleaching). While the taxonomic composition rapidly turned into a new state following the disturbance, the metabolic profile of the community and system performance remained remarkably stable. A total of 920 metagenome-assembled genomes (MAGs), representing approximately 70% of the community in the studied AS ecosystem, were recovered from the 97 monthly AS metagenomes. Comparative genomic analysis revealed an increased ability to aggregate in the cohorts of MAGs with correlated dynamics that are dominant after the bleaching event. Fine-scale analysis of dynamics also revealed cohorts that dominated during different periods and showed successional dynamics on seasonal and longer time scales due to temperature fluctuation and gradual changes in mean residence time in the reactor, respectively.

CONCLUSIONS

Our work highlights that communities can assume different stable states under highly similar environmental conditions and that a specific disturbance threshold may lead to a rapid shift in community composition. Video Abstract.

Reconstructing large interaction networks from empirical time series data

Reconstructing interactions from observational data is a critical need for investigating natural biological networks, wherein network dimensionality is usually high. However, these pose a challenge to existing methods that can quantify only small interaction networks. Here, we proposed a novel approach to reconstruct high-dimensional interaction Jacobian networks using empirical time series without specific model assumptions. This method, named "multiview distance regularised S-map," generalised the state space reconstruction to accommodate high dimensionality and overcome difficulties in quantifying massive interactions with limited data. When evaluating this method using time series generated from theoretical models involving hundreds of interacting species, estimated strengths of interaction Jacobians were in good agreement with theoretical expectations. Applying this method to a natural bacterial community helped identify important species from the interaction network and revealed mechanisms governing the dynamical stability of a bacterial community. The proposed method overcame the challenge of high dimensionality in large natural dynamical systems.

Fluctuation spectra of large random dynamical systems reveal hidden structure in ecological networks

Understanding the relationship between complexity and stability in large dynamical systems-such as ecosystems-remains a key open question in complexity theory which has inspired a rich body of work developed over more than fifty years. The vast majority of this theory addresses asymptotic linear stability around equilibrium points, but the idea of 'stability' in fact has other uses in the empirical ecological literature. The important notion of 'temporal stability' describes the character of fluctuations in population dynamics, driven by intrinsic or extrinsic noise. Here we apply tools from random matrix theory to the problem of temporal stability, deriving analytical predictions for the fluctuation spectra of complex ecological networks. We show that different network structures leave distinct signatures in the spectrum of fluctuations, and demonstrate the application of our theory to the analysis of ecological time-series data of plankton abundances.