Phosphate acquisition genes in Prochlorococcus ecotypes: evidence for genome-wide adaptation

Bacteriophages carry auxiliary metabolic genes related to energy, sulfur and phosphorus metabolism during a harmful algal bloom in a freshwater lake

Cyanophages play an important role in nutrient cycling in lakes since they can modulate the metabolism of cyanobacteria. A proper understanding of the impact of cyanophage infection on the metabolism and ecology of cyanobacteria is critical during a complete cycle of harmful algal bloom (HAB). The ecology of cyanophages in marine environments has been well-delineated, but cyanophages in freshwater lakes remain less studied. Here, we studied the diversity of cyanophages and their impact on host ecology and metabolism through the succession of HAB in Utah Lake, which is a shallow eutrophic freshwater lake, in 2019. We collected water samples at three different periods from two locations in freshwater Utah Lake. The three sampling periods represented the pre-bloom, peak-bloom, and post-bloom events. We observed that the Utah Lake virome was dominated by families Myoviridae, Siphoviridae, and Podoviridae under the order Caudovirales. We detected photosystem-related genes, sulfur assimilation genes, and pho regulon (phosphorus metabolism) genes in genomes of predicted cyanophages. We were able to capture the changes in relative abundance and expression of functional genes in genomes of cyanophage at different stages of the bloom. We observed higher relative abundance and expression of cyanophage-encoded pho-regulon genes in the "pre-bloom" period. The higher expression of pho-regulon genes in P-limited ecosystem of Utah Lake indicated the possible contribution of cyanophage to enhance the fitness of the host cyanobacteria. Our study provides some insightful findings on the role of cyanophages in controlling the ecology and relative abundance of host cyanobacteria in freshwater lakes.

Expression and characterization of the complete cyanophage genome PP in the heterologous host Synechococcus elongatus PCC 7942

In this study, we successfully integrated the full-length genome of the cyanophage PP into the non-host cyanobacterium Synechococcus elongatus PCC 7942, facilitated by conjugation via Escherichia coli. To address the challenge posed by the toxic open reading frames (ORFs) of PP in E. coli, we first identified and characterized three toxic ORFs. The PP genome was subsequently rearranged, and the expression of these toxic ORFs was controlled using a tandem-induction switch system. The full-length PP genome was then successfully integrated into the genome of S. elongatus PCC 7942. Interestingly, the integration of the PP genome led to a reduction in photosynthesis and carbon fixation in S. elongatus PCC 7942, resembling the effects typically associated with cyanophage infection. Transcriptomic analysis showed that 32 of the 41 ORFs in the PP genome were actively transcribed in S. elongatus PCC 7942, significantly affecting energy metabolism and carbon fixation pathways. These effects were further confirmed by metabolomic analysis.

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.

Adaptive Responses of Cyanobacteria to Phosphate Limitation: A Focus on Marine Diazotrophs

Phosphorus is an essential component of numerous macromolecules and is vital for life. Its availability significantly influences primary production, particularly in oligotrophic environments. Marine diazotrophic cyanobacteria, which play key roles in biogeochemical cycles through nitrogen fixation (N2 fixation), have adapted to thrive in phosphate (Pi)-poor areas. However, the molecular mechanisms that facilitate their adaptation to such conditions remain incompletely understood. Bacteria have evolved various strategies to cope with Pi limitation, including detecting Pi availability, utilising high-affinity Pi transporters, and hydrolyzing dissolved organic phosphorus (DOP) with various enzymes. This review synthesises current knowledge regarding how cyanobacteria adapt to Pi scarcity, with particular emphasis on subtropical marine free-living diazotrophs and their ability to utilise diverse DOP molecules. Omics approaches, such as (meta)genomics and (meta)transcriptomics, reveal the resilience of marine diazotrophs in the face of Pi scarcity and highlight the need for further research into their molecular adaptive strategies. Adaptation to Pi limitation is often intertwined with the broader response of cyanobacteria to multiple limitations and stresses. This underscores the importance of understanding Pi adaptation to assess the ecological resilience of these crucial microorganisms in dynamic environments, particularly in the context of global climate change.

Enhanced carbon use efficiency and warming resistance of soil microorganisms under organic amendment

The frequency and intensity of extreme weather events, including rapid temperature fluctuations, are increasing because of climate change. Long-term fertilization practices have been observed to alter microbial physiology and community structure, thereby affecting soil carbon sequestration. However, the effects of warming on the carbon sequestration potential of soil microbes adapted to long-term fertilization remain poorly understood. In this study, we utilized 18O isotope labeling to assess microbial carbon use efficiency (CUE) and employed stable isotope probing (SIP) with 18O-H2O to identify growing taxa in response to temperature changes (5-35 °C). Organic amendment with manure or straw residue significantly increased microbial CUE by 86-181 % compared to unfertilized soils. The microorganisms inhabiting organic amended soils displayed greater resistance of microbial CUE to high temperatures (25-35 °C) compared to those inhabiting soils fertilized only with minerals. Microbial growth patterns determined by the classification of taxa into incorporators or non-incorporators based on 18O incorporation into DNA exhibited limited phylogenetic conservation in response to temperature changes. Microbial clusters were identified by grouping taxa with similar growth patterns across different temperatures. Organic amendments enriched microbial clusters associated with increased CUE, whereas clusters in unfertilized or mineral-only fertilized soils were linked to decreased CUE. Specifically, shifts in the composition of growing bacteria were correlated with enhanced microbial CUE, whereas modifications in the composition of growing fungi were associated with diminished CUE. Notably, the responses of microbial CUE to temperature fluctuations were primarily driven by changes in the bacterial composition. Overall, our findings demonstrate that organic amendments enhance soil microbial CUE and promote the enrichment of specific microbial clusters that are better equipped to cope with temperature changes. This study establishes a theoretical foundation for manipulating soil microbes to enhance carbon sequestration under global climate scenarios.

A distinct, high-affinity, alkaline phosphatase facilitates occupation of P-depleted environments by marine picocyanobacteria

Marine picocyanobacteria of the genera Prochlorococcus and Synechococcus, the two most abundant phototrophs on Earth, thrive in oligotrophic oceanic regions. While it is well known that specific lineages are exquisitely adapted to prevailing in situ light and temperature regimes, much less is known of the molecular machinery required to facilitate occupancy of these low-nutrient environments. Here, we describe a hitherto unknown alkaline phosphatase, Psip1, that has a substantially higher affinity for phosphomonoesters than other well-known phosphatases like PhoA, PhoX, or PhoD and is restricted to clade III Synechococcus and a subset of high light I-adapted Prochlorococcus strains, suggesting niche specificity. We demonstrate that Psip1 has undergone convergent evolution with PhoX, requiring both iron and calcium for activity and likely possessing identical key residues around the active site, despite generally very low sequence homology. Interrogation of metagenomes and transcriptomes from TARA oceans and an Atlantic Meridional transect shows that psip1 is abundant and highly expressed in picocyanobacterial populations from the Mediterranean Sea and north Atlantic gyre, regions well recognized to be phosphorus (P)-deplete. Together, this identifies psip1 as an important oligotrophy-specific gene for P recycling in these organisms. Furthermore, psip1 is not restricted to picocyanobacteria and is abundant and highly transcribed in some α-proteobacteria and eukaryotic algae, suggesting that such a high-affinity phosphatase is important across the microbial taxonomic world to occupy low-P environments.

Nitrogen uptake rates and phytoplankton composition across contrasting North Atlantic Ocean coastal regimes north and south of Cape Hatteras

Understanding nitrogen (N) uptake rates respect to nutrient availability and the biogeography of phytoplankton communities is crucial for untangling the complexities of marine ecosystems and the physical, biological, and chemical forces shaping them. In the summer of 2016, we conducted measurements of bulk microbial uptake rates for six 15N-labeled substrates: nitrate, nitrite, ammonium, urea, cyanate, and dissolve free amino acids across distinct marine provinces, including the continental shelf of the Mid-and South Atlantic Bights (MAB and SAB), the Slope Sea, and the Gulf Stream, marking the first instance of simultaneously measuring six different N uptake rates in this dynamic region. Total measured N uptake rates were lowest in the Gulf Stream followed by the SAB. Notably, the MAB exhibited significantly higher N uptake rates compared to the SAB, likely due to the excess levels of pre-existing phosphorus present in the MAB. Together, urea and nitrate uptake contributed approximately 50% of the total N uptake across the study region. Although cyanate uptake rates were consistently low, they accounted for up to 11% of the total measured N uptake at some Gulf Stream stations. Phytoplankton groups were identified based on specific pigment markers, revealing a dominance of diatoms in the shelf community, while Synechococcus, Prochlorococcus, and pico-eukaryotes dominated in oligotrophic Gulf Stream waters. The reported uptake rates in this study were mostly in agreement with previous studies conducted in coastal waters of the North Atlantic Ocean. This study suggests there are distinct regional patterns of N uptake in this physically dynamic region, correlating with nutrient availability and phytoplankton community composition. These findings contribute valuable insights into the intricate interplay of biological and chemical factors shaping N dynamics in disparate marine ecosystems.

Substrate specificity and ecological significance of PstS homologs in phosphorus uptake in marine Synechococcus sp. WH8102

Phosphorus, a vital macronutrient, often limits primary productivity in marine environments. Marine Synechococcus strains, including WH8102, rely on high-affinity phosphate-binding proteins (PstS) to scavenge inorganic phosphate in oligotrophic oceans. However, WH8102 possesses three distinct PstS homologs whose substrate specificity and ecological roles are unclear. The three PstS homologs were heterologously expressed and purified to investigate their substrate specificity and binding kinetics. Our study revealed that all three PstS homologs exhibited a high degree of specificity for phosphate but differed in phosphate binding affinities. Notably, PstS1b displayed nearly 10-fold higher binding affinity (KD = 0.44 µM) compared to PstS1a (KD = 3.3 μM) and PstS2 (KD = 4.3 μM). Structural modeling suggested a single amino acid variation in the binding pocket of PstS1b (threonine instead of serine in PstS1a and PstS2) likely contributed to its higher Pi affinity. Genome context data, together with the protein biophysical data, suggest distinct ecological roles for the three PstS homologs. We propose that PstS1b may be involved in scavenging inorganic phosphorus in oligotrophic conditions and that PstS1a may be involved in transporting recycled phosphate derived from organic phosphate cleavage. The role of PstS2 is less clear, but it may be involved in phosphate uptake when environmental phosphate concentrations are transiently higher. The conservation of three distinct PstS homologs in Synechococcus clade III strains likely reflects distinct adaptations for P acquisition under varying oligotrophic conditions.IMPORTANCEPhosphorus is an essential macronutrient that plays a key role in marine primary productivity and biogeochemistry. However, intense competition for bioavailable phosphorus in the marine environment limits growth and productivity of ecologically important cyanobacteria. In oligotrophic oceans, marine Synechococcus strains, like WH8102, utilize high-affinity phosphate-binding proteins (PstS) to scavenge inorganic phosphate. However, WH8102 possesses three distinct PstS homologs, with unclear substrate specificity and ecological roles, creating a knowledge gap in understanding phosphorus acquisition mechanisms in picocyanobacteria. Through genomic, functional, biophysical, and structural analysis, our study unravels the ecological functions of these homologs. Our findings enhance our understanding of cyanobacterial nutritional uptake strategies and shed light on the crucial role of these conserved nutrient uptake systems in adaptation to specific niches, which ultimately underpins the success of marine Synechococcus across a diverse array of marine ecosystems.

Exploring virus-host-environment interactions in a chemotrophic-based underground estuary

BACKGROUND

Viruses play important roles in modulating microbial communities and influencing global biogeochemistry. There is now growing interest in characterising their ecological roles across diverse biomes. However, little is known about viral ecology in low-nutrient, chemotrophic-based environments. In such ecosystems, virus-driven manipulation of nutrient cycles might have profound impacts across trophic levels. In particular, anchialine environments, which are low-energy underground estuaries sustained by chemotrophic processes, represent ideal model systems to study novel virus-host-environment interactions.

RESULTS

Here, we employ metagenomic sequencing to investigate the viral community in Bundera Sinkhole, an anchialine ecosystem rich in endemic species supported by microbial chemosynthesis. We find that the viruses are highly novel, with less than 2% representing described viruses, and are hugely abundant, making up as much as 12% of microbial intracellular DNA. These highly abundant viruses largely infect important prokaryotic taxa that drive key metabolic processes in the sinkhole. Further, the abundance of viral auxiliary metabolic genes (AMGs) involved in nucleotide and protein synthesis was strongly correlated with declines in environmental phosphate and sulphate concentrations. These AMGs encoded key enzymes needed to produce sulphur-containing amino acids, and phosphorus metabolic enzymes involved in purine and pyrimidine nucleotide synthesis. We hypothesise that this correlation is either due to selection of these AMGs under low phosphate and sulphate concentrations, highlighting the dynamic interactions between viruses, their hosts, and the environment; or, that these AMGs are driving increased viral nucleotide and protein synthesis via manipulation of host phosphorus and sulphur metabolism, consequently driving nutrient depletion in the surrounding water.

CONCLUSION

This study represents the first metagenomic investigation of viruses in anchialine ecosystems, and provides new hypotheses and insights into virus-host-environment interactions in such 'dark', low-energy environments. This is particularly important since anchialine ecosystems are characterised by diverse endemic species, both in their microbial and faunal assemblages, which are primarily supported by microbial chemosynthesis. Thus, virus-host-environment interactions could have profound effects cascading through all trophic levels.

Coordinated transcriptional response to environmental stress by a Synechococcus virus

Viruses are a major control on populations of microbes. Often, their virulence is examined in controlled laboratory conditions. Yet, in nature, environmental conditions lead to changes in host physiology and fitness that may impart both costs and benefits on viral success. Phosphorus (P) is a major abiotic control on the marine cyanobacterium Synechococcus. Some viruses infecting Synechococcus have acquired, from their host, a gene encoding a P substrate binding protein (PstS), thought to improve virus replication under phosphate starvation. Yet, pstS is uncommon among cyanobacterial viruses. Thus, we asked how infections with viruses lacking PstS are affected by P scarcity. We show that the production of infectious virus particles of such viruses is reduced in low P conditions. However, this reduction in progeny is not caused by impaired phage genome replication, thought to be a major sink for cellular phosphate. Instead, transcriptomic analysis showed that under low P conditions, a PstS-lacking cyanophage increased the expression of a specific gene set that included mazG, hli2, and gp43 encoding a pyrophosphatase, a high-light inducible protein and DNA polymerase, respectively. Moreover, several of the upregulated genes were controlled by the host's phoBR two-component system. We hypothesize that recycling and polymerization of nucleotides liberates free phosphate and thus allows viral morphogenesis, albeit at lower rates than when phosphate is replete or when phages encode pstS. Altogether, our data show how phage genomes, lacking obvious P-stress-related genes, have evolved to exploit their host's environmental sensing mechanisms to coordinate their own gene expression in response to resource limitation.

Phosphate-related genomic islands as drivers of environmental adaptation in the streamlined marine alphaproteobacterial HIMB59

IMPORTANCE

These results shed light on the evolutionary strategies of microbes with streamlined genomes to adapt and survive in the oligotrophic conditions that dominate the surface waters of the global ocean. At the individual level, these microbes have been subjected to evolutionary constraints that have led to a more efficient use of nutrients, removing non-essential genes named as "streamlining theory." However, at the population level, they conserve a highly diverse gene pool in flexible genomic islands resulting in polyclonal populations on the same genomic background as an evolutionary response to environmental pressures. Localization of these islands at equivalent positions in the genome facilitates horizontal transfer between clonal lineages. This high level of environmental genomic heterogeneity could explain their cosmopolitan distribution. In the case of the order HIMB59 within the class Alphaproteobacteria, two factors exert evolutionary pressure and determine this intraspecific diversity: phages and the concentration of P in the environment.

Global scale phylogeography of functional traits and microdiversity in Prochlorococcus

Prochlorococcus is the most numerically abundant photosynthetic organism in the surface ocean. The Prochlorococcus high-light and warm-water adapted ecotype (HLII) is comprised of extensive microdiversity, but specific functional differences between microdiverse sub-clades remain elusive. Here we characterized both functional and phylogenetic diversity within the HLII ecotype using Bio-GO-SHIP metagenomes. We found widespread variation in gene frequency connected to local environmental conditions. Metagenome-assembled marker genes and genomes revealed a globally distributed novel HLII haplotype defined by adaptation to chronically low P conditions (HLII-P). Environmental correlation analysis revealed different factors were driving gene abundances verses phylogenetic differences. An analysis of cultured HLII genomes and metagenome-assembled genomes revealed a subclade within HLII, which corresponded to the novel HLII-P haplotype. This work represents the first global assessment of the HLII ecotype's phylogeography and corresponding functional differences. These findings together expand our understanding of how microdiversity structures functional differences and reveals the importance of nutrients as drivers of microdiversity in Prochlorococcus.

Intra-population genomic diversity of the bloom-forming cyanobacterium, Aphanizomenon gracile, at low spatial scale

Cyanobacteria are oxygenic photosynthetic bacteria that perform a substantial part of the global primary production. Some species are responsible for catastrophic environmental events, called blooms, which have become increasingly common in lakes and freshwater bodies as a consequence of global changes. Genotypic diversity is considered essential for marine cyanobacterial population, allowing it to cope with spatio-temporal environmental variations and to adapt to specific micro-niches in the ecosystem. This aspect is underestimated in the study of bloom development, however, and given little notice in studies of the ecology of harmful cyanobacteria. Here we compared the genomes of four strains of Aphanizomenon gracile, a species of filamentous toxinogenic cyanobacteria (Nostocales) found worldwide in fresh and brackish water. Millimeter-sized fascicles were isolated from a single water sample and have been maintained in culture since 2010. A comparative study revealed extensive heterogeneity in gene contents, despite similar genome size and high similarity indices. These variations were mainly associated with mobile genetic elements and biosynthetic gene clusters. For some of the latter, metabolomic analysis confirmed the production of related secondary metabolites, such as cyanotoxins and carotenoids, which are thought to play a fundamental role in the cyanobacterial fitness. Altogether, these results demonstrated that an A. gracile bloom could be a highly diverse population at low spatial scale and raised questions about potential exchanges of essential metabolites between individuals.

Marine picocyanobacterial PhnD1 shows specificity for various phosphorus sources but likely represents a constitutive inorganic phosphate transporter

Despite being fundamental to multiple biological processes, phosphorus (P) availability in marine environments is often growth-limiting, with generally low surface concentrations. Picocyanobacteria strains encode a putative ABC-type phosphite/phosphate/phosphonate transporter, phnDCE, thought to provide access to an alternative phosphorus pool. This, however, is paradoxical given most picocyanobacterial strains lack known phosphite degradation or carbon-phosphate lyase pathway to utilise alternate phosphorus pools. To understand the function of the PhnDCE transport system and its ecological consequences, we characterised the PhnD1 binding proteins from four distinct marine Synechococcus isolates (CC9311, CC9605, MITS9220, and WH8102). We show the Synechococcus PhnD1 proteins selectively bind phosphorus compounds with a stronger affinity for phosphite than for phosphate or methyl phosphonate. However, based on our comprehensive ligand screening and growth experiments showing Synechococcus strains WH8102 and MITS9220 cannot utilise phosphite or methylphosphonate as a sole phosphorus source, we hypothesise that the picocyanobacterial PhnDCE transporter is a constitutively expressed, medium-affinity phosphate transporter, and the measured affinity of PhnD1 to phosphite or methyl phosphonate is fortuitous. Our MITS9220_PhnD1 structure explains the comparatively lower affinity of picocyanobacterial PhnD1 for phosphate, resulting from a more limited H-bond network. We propose two possible physiological roles for PhnD1. First, it could function in phospholipid recycling, working together with the predicted phospholipase, TesA, and alkaline phosphatase. Second, by having multiple transporters for P (PhnDCE and Pst), picocyanobacteria could balance the need for rapid transport during transient episodes of higher P availability in the environment, with the need for efficient P utilisation in typical phosphate-deplete conditions.

Differential global distribution of marine picocyanobacteria gene clusters reveals distinct niche-related adaptive strategies

The ever-increasing number of available microbial genomes and metagenomes provides new opportunities to investigate the links between niche partitioning and genome evolution in the ocean, especially for the abundant and ubiquitous marine picocyanobacteria Prochlorococcus and Synechococcus. Here, by combining metagenome analyses of the Tara Oceans dataset with comparative genomics, including phyletic patterns and genomic context of individual genes from 256 reference genomes, we show that picocyanobacterial communities thriving in different niches possess distinct gene repertoires. We also identify clusters of adjacent genes that display specific distribution patterns in the field (eCAGs) and are thus potentially involved in the same metabolic pathway and may have a key role in niche adaptation. Several eCAGs are likely involved in the uptake or incorporation of complex organic forms of nutrients, such as guanidine, cyanate, cyanide, pyrimidine, or phosphonates, which might be either directly used by cells, for example for the biosynthesis of proteins or DNA, or degraded to inorganic nitrogen and/or phosphorus forms. We also highlight the enrichment of eCAGs involved in polysaccharide capsule biosynthesis in Synechococcus populations thriving in both nitrogen- and phosphorus-depleted areas vs. low-iron (Fe) regions, suggesting that the complexes they encode may be too energy-consuming for picocyanobacteria thriving in the latter areas. In contrast, Prochlorococcus populations thriving in Fe-depleted areas specifically possess an alternative respiratory terminal oxidase, potentially involved in the reduction of Fe(III) to Fe(II). Altogether, this study provides insights into how phytoplankton communities populate oceanic ecosystems, which is relevant to understanding their capacity to respond to ongoing climate change.

Global scale phylogeography of functional traits and microdiversity in Prochlorococcus

Prochlorococcus is the most numerically abundant photosynthetic organism in the surface ocean. The Prochlorococcus high-light and warm-water adapted ecotype (HLII) is comprised of extensive microdiversity, but specific functional differences between microdiverse sub-clades remain elusive. Here we characterized both functional and phylogenetic diversity within the HLII ecotype using Bio-GO-SHIP metagenomes. We found widespread variation in gene frequency connected to local environmental conditions. Metagenomically assembled marker genes and genomes revealed a globally distributed novel HLII haplotype defined by adaptation to chronically low P conditions (HLII-P). Environmental correlation analysis revealed different factors were driving gene abundances verses phylogenetic differences. An analysis of cultured HLII genomes and metagenomically assembled genomes revealed a subclade within HLII, which corresponded to the novel HLII-P haplotype. This work represents the first global assessment of the HLII ecotype’s phylogeography and corresponding functional differences. These findings together expand our understanding of how microdiversity structures functional differences and reveals the importance of nutrients as drivers of microdiversity in Prochlorococcus .

Novel integrative elements and genomic plasticity in ocean ecosystems

Horizontal gene transfer accelerates microbial evolution. The marine picocyanobacterium Prochlorococcus exhibits high genomic plasticity, yet the underlying mechanisms are elusive. Here, we report a novel family of DNA transposons-"tycheposons"-some of which are viral satellites while others carry cargo, such as nutrient-acquisition genes, which shape the genetic variability in this globally abundant genus. Tycheposons share distinctive mobile-lifecycle-linked hallmark genes, including a deep-branching site-specific tyrosine recombinase. Their excision and integration at tRNA genes appear to drive the remodeling of genomic islands-key reservoirs for flexible genes in bacteria. In a selection experiment, tycheposons harboring a nitrate assimilation cassette were dynamically gained and lost, thereby promoting chromosomal rearrangements and host adaptation. Vesicles and phage particles harvested from seawater are enriched in tycheposons, providing a means for their dispersal in the wild. Similar elements are found in microbes co-occurring with Prochlorococcus, suggesting a common mechanism for microbial diversification in the vast oligotrophic oceans.

Seasonal bacterial niche structures and chemolithoautotrophic ecotypes in a North Atlantic fjord

Quantifying the temporal change of bacterial communities is essential to understanding how both natural and anthropogenic pressures impact the functions of coastal marine ecosystems. Here we use weekly microbial DNA sampling across four years to show that bacterial phyla have distinct seasonal niches, with a richness peak in winter (i.e., an inverse relationship with daylength). Our results suggest that seasonal fluctuations, rather than the kinetic energy or resource hypotheses, dominated the pattern of bacterial diversity. These findings supplement those from global analyses which lack temporal replication and present few data from winter months in polar and temperate regions. Centered log-ratio transformed data provided new insights into the seasonal niche partitioning of conditionally rare phyla, such as Modulibacteria, Verrucomicrobiota, Synergistota, Deinococcota, and Fermentibacterota. These patterns could not be identified using the standard practice of ASV generation followed by rarefaction. Our study provides evidence that five globally relevant ecotypes of chemolithoautotrophic bacteria from the SUP05 lineage comprise a significant functional group with varying seasonal dominance patterns in the Bedford Basin.

Siderophores as an iron source for picocyanobacteria in deep chlorophyll maximum layers of the oligotrophic ocean

Prochlorococcus and Synechococcus are the most abundant photosynthesizing organisms in the oceans. Gene content variation among picocyanobacterial populations in separate ocean basins often mirrors the selective pressures imposed by the region's distinct biogeochemistry. By pairing genomic datasets with trace metal concentrations from across the global ocean, we show that the genomic capacity for siderophore-mediated iron uptake is widespread in Synechococcus and low-light adapted Prochlorococcus populations from deep chlorophyll maximum layers of iron-depleted regions of the oligotrophic Pacific and S. Atlantic oceans: Prochlorococcus siderophore consumers were absent in the N. Atlantic ocean (higher new iron flux) but constituted up to half of all Prochlorococcus genomes from metagenomes in the N. Pacific (lower new iron flux). Picocyanobacterial siderophore consumers, like many other bacteria with this trait, also lack siderophore biosynthesis genes indicating that they scavenge exogenous siderophores from seawater. Statistical modeling suggests that the capacity for siderophore uptake is endemic to remote ocean regions where atmospheric iron fluxes are the smallest, especially at deep chlorophyll maximum and primary nitrite maximum layers. We argue that abundant siderophore consumers at these two common oceanographic features could be a symptom of wider community iron stress, consistent with prior hypotheses. Our results provide a clear example of iron as a selective force driving the evolution of marine picocyanobacteria.

Biochemical and structural characterization of the cyanophage-encoded phosphate binding protein: implications for enhanced phosphate uptake of infected cyanobacteria

To acquire phosphorus, cyanobacteria use the typical bacterial ABC-type phosphate transporter, which is composed of a periplasmic high-affinity phosphate-binding protein PstS and a channel formed by two transmembrane proteins PstC and PstA. A putative pstS gene was identified in the genomes of cyanophages that infect the unicellular marine cyanobacteria Prochlorococcus and Synechococcus. However, it has not been determined whether the cyanophage PstS protein is functional during infection to enhance the phosphate uptake rate of host cells. Here we showed that the cyanophage P-SSM2 PstS protein was abundant in the infected Prochlorococcus NATL2A cells and the host phosphate uptake rate was enhanced after infection. This is consistent with our biochemical and structural analyses showing that the phage PstS protein is indeed a high-affinity phosphate-binding protein. We further modeled the complex structure of phage PstS with host PstCA and revealed three putative interfaces that may facilitate the formation of a chimeric ABC transporter. Our results provide insights into the molecular mechanism by which cyanophages enhance the phosphate uptake rate of cyanobacteria. Phosphate acquisition by infected bacteria can increase the phosphorus contents of released cellular debris and virus particles, which together constitute a significant proportion of the marine dissolved organic phosphorus pool. This article is protected by copyright. All rights reserved.

Patterns and ecological drivers of viral communities in acid mine drainage sediments across Southern China

Recent advances in environmental genomics have provided unprecedented opportunities for the investigation of viruses in natural settings. Yet, our knowledge of viral biogeographic patterns and the corresponding drivers is still limited. Here, we perform metagenomic deep sequencing on 90 acid mine drainage (AMD) sediments sampled across Southern China and examine the biogeography of viruses in this extreme environment. The results demonstrate that prokaryotic communities dictate viral taxonomic and functional diversity, abundance and structure, whereas other factors especially latitude and mean annual temperature also impact viral populations and functions. In silico predictions highlight lineage-specific virus-host abundance ratios and richness-dependent virus-host interaction structure. Further functional analyses reveal important roles of environmental conditions and horizontal gene transfers in shaping viral auxiliary metabolic genes potentially involved in phosphorus assimilation. Our findings underscore the importance of both abiotic and biotic factors in predicting the taxonomic and functional biogeographic dynamics of viruses in the AMD sediments.

The biogeography of viral communities in extreme environments remains understudied. Here, the authors use metagenomic sequencing on 90 acid mine drainage sediments sampled across Southern China, showing the predominant effects of prokaryotic communities and the influence of environmental variables on viral taxonomy and function.

Distribution of soil viruses across China and their potential role in phosphorous metabolism

BACKGROUND

Viruses are the most abundant biological entities on the planet and drive biogeochemical cycling on a global scale. Our understanding of biogeography of soil viruses and their ecological functions lags significantly behind that of Bacteria and Fungi. Here, a viromic approach was used to investigate the distribution and ecological functions of viruses from 19 soils across China.

RESULTS

Soil viral community were clustered more significantly by geographical location than type of soil (agricultural and natural). Three clusters of viral communities were identified from North, Southeast and Southwest regions; these clusters differentiated using taxonomic composition and were mainly driven by geographic location and climate factors. A total of 972 viral populations (vOTUs) were detected spanning 23 viral families from the 19 viromes. Phylogenetic analyses of the phoH gene showed a remarkable diversity and the distribution of viral phoH genes was more dependent on the environment. Notably, five proteins involved in phosphorus (P) metabolism-related nucleotide synthesis functions, including dUTPase, MazG, PhoH, Thymidylate synthase complementing protein (Thy1), and Ribonucleoside reductase (RNR), were mainly identified in agricultural soils.

CONCLUSIONS

The present work revealed that soil viral communities were distributed across China according to geographical location and climate factors. In addition, P metabolism genes encoded by these viruses probably drive the synthesis of nucleotides for their own genomes inside bacterial hosts, thereby affecting P cycling in the soil ecosystems.

The ecological assembly of bacterial communities in Antarctic wetlands varies across levels of phylogenetic resolution

As functional traits are conserved at different phylogenetic depths, the ability to detect community assembly processes can be conditional on the phylogenetic resolution; yet most previous work quantifying their influence has focused on a single level of phylogenetic resolution. Here, we have studied the ecological assembly of bacterial communities from an Antarctic wetland complex, applying null models across different levels of phylogenetic resolution (i.e. clustering ASVs into OTUs with decreasing sequence identity thresholds). We found that the relative influence of the community assembly processes varies with phylogenetic resolution. More specifically, selection processes seem to impose stronger influence at finer (100% sequence similarity ASV) than at coarser (99%–97% sequence similarity OTUs) resolution. We identified environmental features related with the ecological processes and propose a conceptual model for the bacterial community assembly in this Antarctic ecosystem. Briefly, eco‐evolutionary processes appear to be leading to different but very closely related ASVs in lotic, lentic and terrestrial environments. In all, this study shows that assessing community assembly processes at different phylogenetic resolutions is key to improve our understanding of microbial ecology. More importantly, a failure to detect selection processes at coarser phylogenetic resolution does not imply the absence of such processes at finer resolutions.

Taxon-Specific Shifts in Bacterial and Archaeal Transcription of Dissolved Organic Matter Cycling Genes in a Stratified Fjord

A considerable fraction of organic matter derived from photosynthesis in the euphotic zone settles into the ocean’s interior and, as it progresses, is degraded by diverse microbial consortia that utilize a suite of extracellular enzymes and membrane transporters. Still, the molecular details that regulate carbon cycling across depths remain little explored. As stratification in fjords has made them attractive models to explore patterns in biological oceanography, we here analyzed bacterial and archaeal transcription in samples from five depth layers in the Gullmar Fjord, Sweden. Transcriptional variation over depth correlated with gradients in chlorophyll a and nutrient concentrations. Differences in transcription between sampling dates (summer and early autumn) were strongly correlated with ammonium concentrations, which potentially was linked with a stronger influence of (micro-)zooplankton grazing in summer. Transcriptional investment in carbohydrate-active enzymes (CAZymes) decreased with depth and shifted toward peptidases, partly a result of elevated CAZyme transcription by Flavobacteriales, Cellvibrionales, and Synechococcales at 2 to 25 m and a dominance of peptidase transcription by Alteromonadales and Rhodobacterales from 50 m down. In particular, CAZymes for chitin, laminarin, and glycogen were important. High levels of transcription of ammonium transporter genes by Thaumarchaeota at depth (up to 18% of total transcription), along with the genes for ammonia oxidation and CO₂ fixation, indicated that chemolithoautotrophy contributed to the carbon flux in the fjord. The taxon-specific expression of functional genes for processing of the marine pool of dissolved organic matter and inorganic nutrients across depths emphasizes the importance of different microbial foraging mechanisms over spatiotemporal scales for shaping biogeochemical cycles.

IMPORTANCE It is generally recognized that stratification in the ocean strongly influences both the community composition and the distribution of ecological functions of microbial communities, which in turn are expected to shape the biogeochemical cycling of essential elements over depth. Here, we used metatranscriptomics analysis to infer molecular detail on the distribution of gene systems central to the utilization of organic matter in a stratified marine system. We thereby uncovered that pronounced shifts in the transcription of genes encoding CAZymes, peptidases, and membrane transporters occurred over depth among key prokaryotic orders. This implies that sequential utilization and transformation of organic matter through the water column is a key feature that ultimately influences the efficiency of the biological carbon pump.

Novel PhoH-encoding vibriophages with lytic activity against environmental Vibrio strains

Cholera is a devastating diarrheal disease that accounts for more than 10% of children's lives worldwide, but its treatment is hampered by a rise in antibiotic resistance. One promising alternative to antibiotic therapy is the use of bacteriophages to treat antibiotic-resistant cholera infections, and control Vibrio cholera in clinical cases and in the environment, respectively. Here, we report four novel, closely related environmental myoviruses, VP4, VP6, VP18, and VP24, which we isolated from two environmental toxigenic Vibrio cholerae strains from river Kuja and Usenge beach in Kenya. High-throughput sequencing followed by bioinformatics analysis indicated that the genomes of the four bacteriophages have closely related sequences, with sizes of 148,180 bp, 148,181 bp, 148,179 bp, and 148,179 bp, and a G + C content of 36.4%. The four genomes carry the phoH gene, which is overrepresented in marine cyanophages. The isolated phages displayed a lytic activity against 15 environmental, as well as one clinical, Vibrio cholerae strains. Thus, these novel lytic vibriophages represent potential biocontrol candidates for water decontamination against pathogenic Vibrio cholerae and ought to be considered for future studies of phage therapy.

N-acyl-homoserine-lactones signaling as a critical control point for phosphorus entrapment by multi-species microbial aggregates

Quorum sensing (QS) has been extensively studied in pure stains of microorganisms, but the ecological roles of QS in multi-species microbial aggregates are poorly understood due to the aggregates' heterogeneity and complexity, in particular the phosphorus (P) entrapment, a key aspect of element cycling. Using periphytic biofilm as a microbial-aggregate model, we addressed how QS signaling via N-acyl-homoserine-lactones (AHLs) regulated P entrapment. The most-abundant AHLs detected were C8-HSL, 3OC8-HSL, and C12-HSL, are the primary regulator of P entrapment in the periphytic biofilm. QS signaling-AHL is a beneficial molecule for bacterial growth in periphytic biofilm and the addition of these three AHLs optimized polyphosphate accumulating organisms (PAOs) community. Growth promotion was accompanied by up-regulation of pyrimidine, purine and energy metabolism. Both intra- and extra-cellular P entrapment were enhanced in the addition of AHLs. AHLs increased extracellular polymeric substances (EPS) production to drive extracellular P entrapment, via up-regulating amino acids biosynthesis and amino sugar/nucleotide sugar metabolism. Also, AHLs improved intracellular P entrapment potential by regulating genes involved in inorganic-P accumulation (ppk, ppx) and P uptake and transport (pit, pstSCAB). This proof-of-concept evidence about how QS signaling regulates P entrapment by microbial aggregates paves the way for managing QS to enhance P removal by microbial aggregates in aquatic environments.

Metagenomic analysis reveals global-scale patterns of ocean nutrient limitation

Nutrient supply regulates the activity of phytoplankton, but the global biogeography of nutrient limitation and co-limitation is poorly understood. Prochlorococcus adapt to local environments by gene gains and losses, and we used genomic changes as an indicator of adaptation to nutrient stress. We collected metagenomes from all major ocean regions as part of the Global Ocean Ship-based Hydrographic Investigations Program (Bio-GO-SHIP) and quantified shifts in genes involved in nitrogen, phosphorus, and iron assimilation. We found regional transitions in stress type and severity as well as widespread co-stress. Prochlorococcus stress genes, bottle experiments, and Earth system model predictions were correlated. We propose that the biogeography of multinutrient stress is stoichiometrically linked by controls on nitrogen fixation. Our omics-based description of phytoplankton resource use provides a nuanced and highly resolved description of nutrient stress in the global ocean.

Metabolic pathways inferred from a bacterial marker gene illuminate ecological changes across South Pacific frontal boundaries

Global oceanographic monitoring initiatives originally measured abiotic essential ocean variables but are currently incorporating biological and metagenomic sampling programs. There is, however, a large knowledge gap on how to infer bacterial functions, the information sought by biogeochemists, ecologists, and modelers, from the bacterial taxonomic information (produced by bacterial marker gene surveys). Here, we provide a correlative understanding of how a bacterial marker gene (16S rRNA) can be used to infer latitudinal trends for metabolic pathways in global monitoring campaigns. From a transect spanning 7000 km in the South Pacific Ocean we infer ten metabolic pathways from 16S rRNA gene sequences and 11 corresponding metagenome samples, which relate to metabolic processes of primary productivity, temperature-regulated thermodynamic effects, coping strategies for nutrient limitation, energy metabolism, and organic matter degradation. This study demonstrates that low-cost, high-throughput bacterial marker gene data, can be used to infer shifts in the metabolic strategies at the community scale.

Phosphorus Release and Regeneration Following Laboratory Lysis of Bacterial Cells

The availability of phosphorus limits primary production in large regions of the oceans, and marine microbes use a variety of strategies to overcome this limitation. One strategy is the production of alkaline phosphatase (APase), which allows hydrolysis of larger dissolved organic phosphorus (DOP) compounds in the periplasm or at the cell surface for transport of orthophosphate into the cell. Cell lysis, driven by grazing and viral infection, releases phosphorus-containing cell components, along with active enzymes that could persist after lysis. The importance of this continued enzymatic activity for orthophosphate regeneration is unknown. We used three model bacteria – Escherichia coli K-12 MG1655, Synechococcus sp. WH7803, and Prochlorococcus sp. MED4 – to assess the impact of continued APase activity after cell lysis, via lysozyme treatment, on orthophosphate regeneration. Direct release of orthophosphate scaled with cell size and was reduced under phosphate-starved conditions where APase activity continued for days after lysis. All lysate incubations showed post-lysis orthophosphate generation suggesting phosphatases other than APase maintain activity. Rates of DOP hydrolysis and orthophosphate remineralization varied post-lysis among strains. Escherichia coli K-12 MG1655 rates of remineralization were 0.6 and 1.2 amol cell^–1hr^–1 under deplete and replete conditions; Synechococcus WH7803 lysates ranged from 0.04 up to 0.3 amol cell^–1hr^–1 during phosphorus deplete and replete conditions, respectively, while in Prochlorococcus MED4 lysates, rates were stable at 0.001 amol cell^–1hr^–1 in both conditions. The range of rates of hydrolysis and regeneration underscores the taxonomic and biochemical variability in the process of nutrient regeneration and further highlights the complexity of quantitatively resolving the major fluxes within the microbial loop.

Distribution of Chromophytic Phytoplankton in the Eddy-Induced Upwelling Region of the West Pacific Ocean Revealed Using rbcL Genes

Marine chromophytic phytoplankton are a diverse group of algae and contribute significantly to the total oceanic primary production. However, the spatial distribution of chromophytic phytoplankton is understudied in the West Pacific Ocean (WPO). In this study, we have investigated the community structure and spatial distribution of chromophytic phytoplankton using RuBisCO genes (Form ID rbcL). Our results showed that Haptophyceae, Pelagophyceae, Cyanophyceae, Xanthophyceae, and Bacillariophyceae were the dominant groups. Further, chromophytic phytoplankton can be distinguished between upwelling and non-upwelling zones of the WPO. Surface and 75 m depths of a non-upwelling area were dominated by Prochlorococcus strains, whereas chromophytic phytoplankton were homogenously distributed at the surface layer in the upwelling zone. Meanwhile, Pelagomonas-like sequences were dominant at DCM (75 m) and 150 m depths of the upwelling zone. Non-metric multidimensional scaling (NMDS) analysis did not differentiate between chromophytic phytoplankton in the upwelling and non-upwelling areas, however, it showed clear trends of them at different depths. Further, redundancy analysis (RDA) showed the influence of physicochemical parameters on the distribution of chromophytic phytoplankton. Along with phosphate (p < 0.01), temperature and other dissolved nutrients were important in driving community structure. The upwelling zone was impacted by a decrease in temperature, salinity, and re-supplement of nutrients, where Pelagomonas-like sequences outnumbered other chromophytic groups presented.

The global explosion of eukaryotic algae: The potential role of phosphorus?

There arose one of the most important ecological transitions in Earth's history approximately 750 million years ago during the middle Neoproterozoic Era (1000 to 541 million years ago, Ma). Biomarker evidence suggests that around this time there was a rapid shift from a predominantly bacterial-dominated world to more complex ecosystems governed by eukaryotic primary productivity. The resulting 'Rise of the algae' led to dramatically altered food webs that were much more efficient in terms of nutrient and energy transfer. Yet, what triggered this ecological shift? In this study we examined the theory that it was the alleviation of phosphorus (P) deficiency that gave eukaryotic alga the prime opportunity to flourish. We performed laboratory experiments on the cyanobacterium Synechocystis salina and the eukaryotic algae Tetraselmis suecica and examined their ability to compete for phosphorus. Both these organisms co-occur in modern European coastal waters and are not known to have any allelopathic capabilities. The strains were cultured in mono and mixed cultures in chemostats across a range of dissolved inorganic phosphorus (DIP) concentrations to reflect modern and ancient oceanic conditions of 2 μM P and 0.2 μM P, respectively. Our results show that the cyanobacteria outcompete the algae at the low input (0.2 μM P) treatment, yet the eukaryotic algae were not completely excluded and remained a constant background component in the mixed-culture experiments. Also, despite their relatively large cell size, the algae T. suecica had a high affinity for DIP. With DIP input concentrations resembling modern-day levels (2 μM), the eukaryotic algae could effectively compete against the cyanobacteria in terms of total biomass production. These results suggest that the availability of phosphorus could have influenced the global expansion of eukaryotic algae. However, P limitation does not seem to explain the complete absence of eukaryotic algae in the biomarker record before ca. 750 Ma.

Unlocking the Genomic Taxonomy of the Prochlorococcus Collective

Prochlorococcus is the most abundant photosynthetic prokaryote on our planet. The extensive ecological literature on the Prochlorococcus collective (PC) is based on the assumption that it comprises one single genus comprising the species Prochlorococcus marinus, containing itself a collective of ecotypes. Ecologists adopt the distributed genome hypothesis of an open pan-genome to explain the observed genomic diversity and evolution patterns of the ecotypes within PC. Novel genomic data for the PC prompted us to revisit this group, applying the current methods used in genomic taxonomy. As a result, we were able to distinguish the five genera: Prochlorococcus, Eurycolium, Prolificoccus, Thaumococcus, and Riococcus. The novel genera have distinct genomic and ecological attributes.

Cyanophage host-derived genes reflect contrasting selective pressures with depth in the oxic and anoxic water column of the Eastern Tropical North Pacific

Cyanophages encode host-derived genes that may increase their fitness. We examined the relative abundance of 18 host-derived cyanophages genes in metagenomes and viromes along depth profiles from the Eastern Tropical North Pacific Oxygen Deficient Zone (ETNP ODZ) where Prochlorococcus dominates a secondary chlorophyll maximum within the ODZ. Cyanophages at the oxic primary chlorophyll maximum encoded genes related to light and phosphate stress (psbA, psbD and pstS in T4-like and psbA in T7-like), but the proportion of cyanophage with these genes decreased with depth. The proportion of cyanophage with purine biosynthesis genes increased with depth in T4-like, but not T7-like cyanophages. No additional host-derived genes were found in deep T7-like cyanophages, suggesting that T4-like and T7-like cyanophages have different host-derived gene acquisition strategies, possibly linked to their different genome packaging mechanisms. In contrast to the ETNP, in the oxic North Atlantic T4-like cyanophages encoded psbA and pstS throughout the euphotic zone. Differences in pstS between the ETNP and the North Atlantic stations were consistent with differences in phosphate concentrations in those regimes. We suggest that the low proportion of cyanophage with psbA within the ODZ reflects the stably stratified low-light conditions occupied by their hosts, a Prochlorococcus ecotype endemic to ODZs.

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.

Evolutionary Mechanisms of Long-Term Genome Diversification Associated With Niche Partitioning in Marine Picocyanobacteria

Marine picocyanobacteria of the genera Prochlorococcus and Synechococcus are the most abundant photosynthetic organisms on Earth, an ecological success thought to be linked to the differential partitioning of distinct ecotypes into specific ecological niches. However, the underlying processes that governed the diversification of these microorganisms and the appearance of niche-related phenotypic traits are just starting to be elucidated. Here, by comparing 81 genomes, including 34 new Synechococcus, we explored the evolutionary processes that shaped the genomic diversity of picocyanobacteria. Time-calibration of a core-protein tree showed that gene gain/loss occurred at an unexpectedly low rate between the different lineages, with for instance 5.6 genes gained per million years (My) for the major Synechococcus lineage (sub-cluster 5.1), among which only 0.71/My have been fixed in the long term. Gene content comparisons revealed a number of candidates involved in nutrient adaptation, a large proportion of which are located in genomic islands shared between either closely or more distantly related strains, as identified using an original network construction approach. Interestingly, strains representative of the different ecotypes co-occurring in phosphorus-depleted waters (Synechococcus clades III, WPC1, and sub-cluster 5.3) were shown to display different adaptation strategies to this limitation. In contrast, we found few genes potentially involved in adaptation to temperature when comparing cold and warm thermotypes. Indeed, comparison of core protein sequences highlighted variants specific to cold thermotypes, notably involved in carotenoid biosynthesis and the oxidative stress response, revealing that long-term adaptation to thermal niches relies on amino acid substitutions rather than on gene content variation. Altogether, this study not only deciphers the respective roles of gene gains/losses and sequence variation but also uncovers numerous gene candidates likely involved in niche partitioning of two key members of the marine phytoplankton.

Phosphorus Feast and Famine in Cyanobacteria: Is Luxury Uptake of the Nutrient Just a Consequence of Acclimation to Its Shortage?

To cope with fluctuating phosphorus (P) availability, cyanobacteria developed diverse acclimations, including luxury P uptake (LPU)-taking up P in excess of the current metabolic demand. LPU is underexplored, despite its importance for nutrient-driven rearrangements in aquatic ecosystems. We studied the LPU after the refeeding of P-deprived cyanobacterium Nostoc sp. PCC 7118 with inorganic phosphate (Pi), including the kinetics of Pi uptake, turnover of polyphosphate, cell ultrastructure, and gene expression. The P-deprived cells deployed acclimations to P shortage (reduction of photosynthetic apparatus and mobilization of cell P reserves). The P-starved cells capable of LPU exhibited a biphasic kinetic of the Pi uptake and polyphosphate formation. The first (fast) phase (1-2 h after Pi refeeding) occurred independently of light and temperature. It was accompanied by a transient accumulation of polyphosphate, still upregulated genes encoding high-affinity Pi transporters, and an ATP-dependent polyphosphate kinase. During the second (slow) phase, recovery from P starvation was accompanied by the downregulation of these genes. Our study revealed no specific acclimation to ample P conditions in Nostoc sp. PCC 7118. We conclude that the observed LPU phenomenon does not likely result from the activation of a mechanism specific for ample P conditions. On the contrary, it stems from slow disengagement of the low-P responses after the abrupt transition from low-P to ample P conditions.

Spatial patterns of microbial communities across surface waters of the Great Barrier Reef

Microorganisms are fundamental drivers of biogeochemical cycling, though their contribution to coral reef ecosystem functioning is poorly understood. Here, we infer predictors of bacterioplankton community dynamics across surface-waters of the Great Barrier Reef (GBR) through a meta-analysis, combining microbial with environmental data from the eReefs platform. Nutrient dynamics and temperature explained 41.4% of inter-seasonal and cross-shelf variation in bacterial assemblages. Bacterial families OCS155, Cryomorphaceae, Flavobacteriaceae, Synechococcaceae and Rhodobacteraceae dominated inshore reefs and their relative abundances positively correlated with nutrient loads. In contrast, Prochlorococcaceae negatively correlated with nutrients and became increasingly dominant towards outershelf reefs. Cyanobacteria in Prochlorococcaceae and Synechococcaceae families occupy complementary cross-shelf biogeochemical niches; their abundance ratios representing a potential indicator of GBR nutrient levels. One Flavobacteriaceae-affiliated taxa was putatively identified as diagnostic for ecosystem degradation. Establishing microbial observatories along GBR environmental gradients will facilitate robust assessments of microbial contributions to reef health and inform tipping-points in reef condition.

Prochlorococcus Cells Rely on Microbial Interactions Rather than on Chlorotic Resting Stages To Survive Long-Term Nutrient Starvation

The ability of microorganisms to withstand long periods of nutrient starvation is key to their survival and success under highly fluctuating conditions that are common in nature. Therefore, one would expect this trait to be prevalent among organisms in the nutrient-poor open ocean. Here, we show that this is not the case for Prochlorococcus, a globally abundant and ecologically important marine cyanobacterium. Instead, Prochlorococcus relies on co-occurring heterotrophic bacteria to survive extended phases of nutrient and light starvation. Our results highlight the power of microbial interactions to drive major biogeochemical cycles in the ocean and elsewhere with consequences at the global scale.

Many microorganisms produce resting cells with very low metabolic activity that allow them to survive phases of prolonged nutrient or energy stress. In cyanobacteria and some eukaryotic phytoplankton, the production of resting stages is accompanied by a loss of photosynthetic pigments, a process termed chlorosis. Here, we show that a chlorosis-like process occurs under multiple stress conditions in axenic laboratory cultures of Prochlorococcus, the dominant phytoplankton linage in large regions of the oligotrophic ocean and a global key player in ocean biogeochemical cycles. In Prochlorococcus strain MIT9313, chlorotic cells show reduced metabolic activity, measured as C and N uptake by Nanoscale secondary ion mass spectrometry (NanoSIMS). However, unlike many other cyanobacteria, chlorotic Prochlorococcus cells are not viable and do not regrow under axenic conditions when transferred to new media. Nevertheless, cocultures with a heterotrophic bacterium, Alteromonas macleodii HOT1A3, allowed Prochlorococcus to survive nutrient starvation for months. We propose that reliance on co-occurring heterotrophic bacteria, rather than the ability to survive extended starvation as resting cells, underlies the ecological success of Prochlorococcus.

Alternative strategies of nutrient acquisition and energy conservation map to the biogeography of marine ammonia-oxidizing archaea

Ammonia-oxidizing archaea (AOA) are among the most abundant and ubiquitous microorganisms in the ocean, exerting primary control on nitrification and nitrogen oxides emission. Although united by a common physiology of chemoautotrophic growth on ammonia, a corresponding high genomic and habitat variability suggests tremendous adaptive capacity. Here, we compared 44 diverse AOA genomes, 37 from species cultivated from samples collected across diverse geographic locations and seven assembled from metagenomic sequences from the mesopelagic to hadopelagic zones of the deep ocean. Comparative analysis identified seven major marine AOA genotypic groups having gene content correlated with their distinctive biogeographies. Phosphorus and ammonia availabilities as well as hydrostatic pressure were identified as selective forces driving marine AOA genotypic and gene content variability in different oceanic regions. Notably, AOA methylphosphonate biosynthetic genes span diverse oceanic provinces, reinforcing their importance for methane production in the ocean. Together, our combined comparative physiological, genomic, and metagenomic analyses provide a comprehensive view of the biogeography of globally abundant AOA and their adaptive radiation into a vast range of marine and terrestrial habitats.

Proteomic Profiling of Emiliania huxleyi Using a Three-Dimensional Separation Method Combined with Tandem Mass Spectrometry

Emiliania huxleyi is one of the most abundant marine planktons, and it has a crucial feature in the carbon cycle. However, proteomic analyses of Emiliania huxleyi have not been done extensively. In this study, a three-dimensional liquid chromatography (3D-LC) system consisting of strong cation exchange, high- and low-pH reversed-phase liquid chromatography was established for in-depth proteomic profiling of Emiliania huxleyi. From tryptic proteome digest, 70 fractions were generated and analyzed using liquid chromatography-tandem mass spectrometry. In total, more than 84,000 unique peptides and 10,000 proteins groups were identified with a false discovery rate of ≤0.01. The physicochemical properties of the identified peptides were evaluated. Using ClueGO, approximately 700 gene ontology terms and 15 pathways were defined from the identified protein groups with p-value ≤0.05, covering a wide range of biological processes, cellular components, and molecular functions. Many biological processes associated with CO₂ fixation, photosynthesis, biosynthesis, and metabolic process were identified. Various molecular functions relating to protein binding and enzyme activities were also found. The 3D-LC strategy is a powerful approach for comparative proteomic studies on Emiliania huxleyi to reveal changes in its protein level and related mechanism.

Trait-based approach to bacterial growth efficiency

Bacterial growth efficiency (BGE) is the proportion of assimilated carbon that is converted into biomass and reflects the balance between growth and energetic demands. Often measured as an aggregate property of the community, BGE is highly variable within and across ecosystems. To understand this variation, we first identified how species identity and resource type affect BGE using 20 bacterial isolates belonging to the phylum Proteobacteria that were enriched from north temperate lakes. Using a trait-based approach that incorporated genomic and phenotypic information, we characterized the metabolism of each isolate and tested for predicted trade-offs between growth rate and efficiency. A substantial amount of variation in BGE could be explained at broad (i.e., order, 20%) and fine (i.e., strain, 58%) taxonomic levels. While resource type was a relatively weak predictor across species, it explained >60% of the variation in BGE within a given species. A metabolic trade-off (between maximum growth rate and efficiency) and genomic features revealed that BGE may be a species-specific metabolic property. Our study suggests that genomic and phylogenetic information may help predict aggregate microbial community functions like BGE and the fate of carbon in ecosystems.

New insights into the distributions of nitrogen fixation and diazotrophs revealed by high-resolution sensing and sampling methods

Nitrogen availability limits marine productivity across large ocean regions. Diazotrophs can supply new nitrogen to the marine environment via nitrogen (N₂) fixation, relieving nitrogen limitation. The distributions of diazotrophs and N₂ fixation have been hypothesized to be generally controlled by temperature, phosphorus, and iron availability in the global ocean. However, even in the North Atlantic where most research on diazotrophs and N₂ fixation has taken place, environmental controls remain contentious. Here we measure diazotroph composition, abundance, and activity at high resolution using newly developed underway sampling and sensing techniques. We capture a diazotrophic community shift from Trichodesmium to UCYN-A between the oligotrophic, warm (25-29 °C) Sargasso Sea and relatively nutrient-enriched, cold (13-24 °C) subpolar and eastern American coastal waters. Meanwhile, N₂ fixation rates measured in this study are among the highest ever recorded globally and show significant increase with phosphorus availability across the transition from the Gulf Stream into subpolar and coastal waters despite colder temperatures and higher nitrate concentrations. Transcriptional patterns in both Trichodesmium and UCYN-A indicate phosphorus stress in the subtropical gyre. Over this iron-replete transect spanning the western North Atlantic, our results suggest that temperature is the major factor controlling the diazotrophic community structure while phosphorous drives N₂ fixation rates. Overall, the occurrence of record-high UCYN-A abundance and peak N₂ fixation rates in the cold coastal region where nitrate concentrations are highest (~200 nM) challenges current paradigms on what drives the distribution of diazotrophs and N₂ fixation.

Accumulation of ambient phosphate into the periplasm of marine bacteria is proton motive force dependent

Bacteria acquire phosphate (P_i) by maintaining a periplasmic concentration below environmental levels. We recently described an extracellular P_i buffer which appears to counteract the gradient required for P_i diffusion. Here, we demonstrate that various treatments to outer membrane (OM) constituents do not affect the buffered P_i because bacteria accumulate P_i in the periplasm, from which it can be removed hypo-osmotically. The periplasmic P_i can be gradually imported into the cytoplasm by ATP-powered transport, however, the proton motive force (PMF) is not required to keep P_i in the periplasm. In contrast, the accumulation of P_i into the periplasm across the OM is PMF-dependent and can be enhanced by light energy. Because the conventional mechanism of P_i-specific transport cannot explain P_i accumulation in the periplasm we propose that periplasmic P_i anions pair with chemiosmotic cations of the PMF and millions of accumulated P_i pairs could influence the periplasmic osmolarity of marine bacteria.

The ubiquitous oceanic bacteria harbour an external phosphate buffer for modulating phosphate (P_i) uptake. Here, using both oceanic SAR11, Prochlorococcus and Synechococcus strains as a model, the authors show that the P_i buffer accumulation in the periplasm is proton motive force-dependent and can be enhanced by light energy.

Linking regional shifts in microbial genome adaptation with surface ocean biogeochemistry

Linking 'omics measurements with biogeochemical cycles is a widespread challenge in microbial community ecology. Here, we propose applying genomic adaptation as 'biosensors' for microbial investments to overcome nutrient stress. We then integrate this genomic information with a trait-based model to predict regional shifts in the elemental composition of marine plankton communities. We evaluated this approach using metagenomic and particulate organic matter samples from the Atlantic, Indian and Pacific Oceans. We find that our genome-based trait model significantly improves our prediction of particulate C : P (carbon : phosphorus) across ocean regions. Furthermore, we detect previously unrecognized ocean areas of iron, nitrogen and phosphorus stress. In many ecosystems, it can be very challenging to quantify microbial stress. Thus, a carefully calibrated genomic approach could become a widespread tool for understanding microbial responses to environmental changes and the biogeochemical outcomes. This article is part of the theme issue 'Conceptual challenges in microbial community ecology'.

Phytoplankton transcriptomic and physiological responses to fixed nitrogen in the California current system

Marine phytoplankton are responsible for approximately half of photosynthesis on Earth. However, their ability to drive ocean productivity depends on critical nutrients, especially bioavailable nitrogen (N) which is scarce over vast areas of the ocean. Phytoplankton differ in their preferences for N substrates as well as uptake efficiencies and minimal N requirements relative to other critical nutrients, including iron (Fe) and phosphorus. In this study, we used the MicroTOOLs high-resolution environmental microarray to examine transcriptomic responses of phytoplankton communities in the California Current System (CCS) transition zone to added urea, ammonium, nitrate, and also Fe in the late summer when N depletion is common. Transcript level changes of photosynthetic, carbon fixation, and nutrient stress genes indicated relief of N limitation in many strains of Prochlorococcus, Synechococcus, and eukaryotic phytoplankton. The transcriptomic responses helped explain shifts in physiological and growth responses observed later. All three phytoplankton groups had increased transcript levels of photosynthesis and/or carbon fixation genes in response to all N substrates. However, only Prochlorococcus had decreased transcript levels of N stress genes and grew substantially, specifically after urea and ammonium additions, suggesting that Prochlorococcus outcompeted other community members in these treatments. Diatom transcript levels of carbon fixation genes increased in response to Fe but not to Fe with N which might have favored phytoplankton that were co-limited by N and Fe. Moreover, transcription patterns of closely related strains indicated variability in N utilization, including nitrate utilization by some high-light adapted Prochlorococcus. Finally, up-regulation of urea transporter genes by both Prochlorococcus and Synechococcus in response to filtered deep water suggested a regulatory mechanism other than classic control via the global N regulator NtcA. This study indicated that co-existing phytoplankton strains experience distinct nutrient stresses in the transition zone of the CCS, an understudied region where oligotrophic and coastal communities naturally mix.

Selection Is a Significant Driver of Gene Gain and Loss in the Pangenome of the Bacterial Genus Sulfurovum in Geographically Distinct Deep-Sea Hydrothermal Vents

Microbial genomes have highly variable gene content, and the evolutionary history of microbial populations is shaped by gene gain and loss mediated by horizontal gene transfer and selection. To evaluate the influence of selection on gene content variation in hydrothermal vent microbial populations, we examined 22 metagenome-assembled genomes (MAGs) (70 to 97% complete) from the ubiquitous vent Epsilonbacteraeota genus Sulfurovum that were recovered from two deep-sea hydrothermal vent regions, Axial Seamount in the northeastern Pacific Ocean (13 MAGs) and the Mid-Cayman Rise in the Caribbean Sea (9 MAGs). Genes involved in housekeeping functions were highly conserved across Sulfurovum lineages. However, genes involved in environment-specific functions, and in particular phosphate regulation, were found mostly in Sulfurovum genomes from the Mid-Cayman Rise in the low-phosphate Atlantic Ocean environment, suggesting that nutrient limitation is an important selective pressure for these bacteria. Furthermore, genes that were rare within the pangenome were more likely to undergo positive selection than genes that were highly conserved in the pangenome, and they also appeared to have experienced gene-specific sweeps. Our results suggest that selection is a significant driver of gene gain and loss for dominant microbial lineages in hydrothermal vents and highlight the importance of factors like nutrient limitation in driving microbial adaptation and evolution.IMPORTANCE Microbes can alter their gene content through the gain and loss of genes. However, there is some debate as to whether natural selection or neutral processes play a stronger role in molding the gene content of microbial genomes. In this study, we examined variation in gene content for the Epsilonbacteraeota genus Sulfurovum from deep-sea hydrothermal vents, which are dynamic habitats known for extensive horizontal gene transfer within microbial populations. Our results show that natural selection is a strong driver of Sulfurovum gene content and that nutrient limitation in particular has shaped the Sulfurovum genome, leading to differences in gene content between ocean basins. Our results also suggest that recently acquired genes undergo stronger selection than genes that were acquired in the more distant past. Overall, our results highlight the importance of natural selection in driving the evolution of microbial populations in these dynamic habitats.

A minimum set of regulators to thrive in the ocean

Marine cyanobacteria of the genus Prochlorococcus thrive in high cell numbers throughout the euphotic zones of the world's subtropical and tropical oligotrophic oceans, making them some of the most ecologically relevant photosynthetic microorganisms on Earth. The ecological success of these free-living phototrophs suggests that they are equipped with a regulatory system competent to address many different stress situations. However, Prochlorococcus genomes are compact and streamlined, with the majority encoding only five different sigma factors, five to six two-component systems and eight types of other transcriptional regulators. Here, we summarize the existing information about the functions of these protein regulators, about transcriptomic responses to defined stress conditions, and discuss the current knowledge about riboswitches, RNA-based regulation and the roles of certain metabolites as co-regulators. We focus on the best-studied isolate, Prochlorococcus MED4, but extend to other strains and ecotypes when appropriate, and we include some information gained from metagenomic and metatranscriptomic analyses.

Toward a genetic system in the marine cyanobacterium Prochlorococcus

As the smallest and most abundant primary producer in the oceans, the cyanobacterium Prochlorococcus is of interest to diverse branches of science. For the past 30 years, research on this minimal phototroph has led to a growing understanding of biological organization across multiple scales, from the genome to the global ocean ecosystem. Progress in understanding drivers of its diversity and ecology, as well as molecular mechanisms underpinning its streamlined simplicity, has been hampered by the inability to manipulate these cells genetically. Multiple attempts have been made to develop an efficient genetic transformation method for Prochlorococcus over the years; all have been unsuccessful to date, despite some success with their close relative, Synechococcus. To avoid the pursuit of unproductive paths, we report here what has not worked in our hands, as well as our progress developing a method to screen the most efficient electroporation parameters for optimal DNA delivery into Prochlorococcus cells. We also report a novel protocol for obtaining axenic colonies and a new method for differentiating live and dead cells. The electroporation method can be used to optimize DNA delivery into any bacterium, making it a useful tool for advancing transformation systems in other genetically recalcitrant microorganisms.

Genomic mosaicism underlies the adaptation of marine Synechococcus ecotypes to distinct oceanic iron niches

Phytoplankton are limited by iron (Fe) in ~40% of the world's oceans including high-nutrient low-chlorophyll (HNLC) regions. While low-Fe adaptation has been well-studied in large eukaryotic diatoms, less is known for small, prokaryotic marine picocyanobacteria. This study reveals key physiological and genomic differences underlying Fe adaptation in marine picocyanobacteria. HNLC ecotype CRD1 strains have greater physiological tolerance to low Fe congruent with their expanded repertoire of Fe transporter, storage and regulatory genes compared to other ecotypes. From metagenomic analysis, genes encoding ferritin, flavodoxin, Fe transporters and siderophore uptake genes were more abundant in low-Fe waters, mirroring paradigms of low-Fe adaptation in diatoms. Distinct Fe-related gene repertories of HNLC ecotypes CRD1 and CRD2 also highlight how coexisting ecotypes have evolved independent approaches to life in low-Fe habitats. Synechococcus and Prochlorococcus HNLC ecotypes likewise exhibit independent, genome-wide reductions of predicted Fe-requiring genes. HNLC ecotype CRD1 interestingly was most similar to coastal ecotype I in Fe physiology and Fe-related gene content, suggesting populations from these different biomes experience similar Fe-selective conditions. This work supports an improved perspective that phytoplankton are shaped by more nuanced Fe niches in the oceans than previously implied from mostly binary comparisons of low- versus high-Fe habitats and populations.