A marine heatwave drives significant shifts in pelagic microbiology.

Marine heatwaves (MHWs) cause disruption to marine ecosystems, deleteriously impacting macroflora and fauna. However, effects on microorganisms are relatively unknown despite ocean temperature being a major determinant of assemblage structure. Using data from thousands of Southern Hemisphere samples, we reveal that during an "unprecedented" 2015/16 Tasman Sea MHW, temperatures approached or surpassed the upper thermal boundary of many endemic taxa. Temperate microbial assemblages underwent a profound transition to niche states aligned with sites over 1000 km equatorward, adapting to higher temperatures and lower nutrient conditions bought on by the MHW. MHW conditions also modulate seasonal patterns of microbial diversity and support novel assemblage compositions. The most significant affects of MHWs on microbial assemblages occurred during warmer months, when temperatures exceeded the upper climatological bounds. Trends in microbial response across several MHWs in different locations suggest these are emergent properties of temperate ocean warming, which may facilitate monitoring, prediction and adaptation efforts.

Rhometa: Population recombination rate estimation from metagenomic read datasets.

Prokaryotic evolution is influenced by the exchange of genetic information between species through a process referred to as recombination. The rate of recombination is a useful measure for the adaptive capacity of a prokaryotic population. We introduce Rhometa (https://github.com/sid-krish/Rhometa), a new software package to determine recombination rates from shotgun sequencing reads of metagenomes. It extends the composite likelihood approach for population recombination rate estimation and enables the analysis of modern short-read datasets. We evaluated Rhometa over a broad range of sequencing depths and complexities, using simulated and real experimental short-read data aligned to external reference genomes. Rhometa offers a comprehensive solution for determining population recombination rates from contemporary metagenomic read datasets. Rhometa extends the capabilities of conventional sequence-based composite likelihood population recombination rate estimators to include modern aligned metagenomic read datasets with diverse sequencing depths, thereby enabling the effective application of these techniques and their high accuracy rates to the field of metagenomics. Using simulated datasets, we show that our method performs well, with its accuracy improving with increasing numbers of genomes. Rhometa was validated on a real S. pneumoniae transformation experiment, where we show that it obtains plausible estimates of the rate of recombination. Finally, the program was also run on ocean surface water metagenomic datasets, through which we demonstrate that the program works on uncultured metagenomic datasets.

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.

Defining marine bacterioplankton community assembly rules by contrasting the importance of environmental determinants and biotic interactions.

Bacterioplankton communities govern marine productivity and biogeochemical cycling, yet drivers of bacterioplankton assembly remain unclear. Here, we contrast the relative contribution of deterministic processes (environmental factors and biotic interactions) in driving temporal dynamics of bacterioplankton diversity at three different oceanographic time series locations, spanning 15° of latitude, which are each characterized by different environmental conditions and varying degrees of seasonality. Monthly surface samples (5.5 years) were analysed using 16S rRNA amplicon sequencing. The high- and mid-latitude sites of Maria Island and Port Hacking were characterized by high and intermediate levels of environmental heterogeneity, respectively, with both alpha diversity (72%; 24% of total variation) and beta diversity (32%; 30%) patterns within bacterioplankton assemblages explained by day length, ammonium, and mixed layer depth. In contrast, North Stradbroke Island, a sub-tropical location where environmental conditions are less variable, interspecific interactions were of increased importance in structuring bacterioplankton diversity (alpha: 33%; beta: 26%) with environment only contributing 11% and 13% to predicting diversity, respectively. Our results demonstrate that bacterioplankton diversity is the result of both deterministic environmental and biotic processes and that the importance of these different deterministic processes varies, potential in response to environmental heterogeneity.

Global Phylogeography of Marine Synechococcus in Coastal Areas Reveals Strong Community Shifts.

Marine Synechococcus comprise a numerically and ecologically prominent phytoplankton group, playing a major role in both carbon cycling and trophic networks in all oceanic regions except in the polar oceans. Despite their high abundance in coastal areas, our knowledge of Synechococcus communities in these environments is based on only a few local studies. Here, we use the global metagenome data set of the Ocean Sampling Day (June 21st, 2014) to get a snapshot of the taxonomic composition of coastal Synechococcus communities worldwide, by recruitment on a reference database of 141 picocyanobacterial genomes, representative of the whole Prochlorococcus, Synechococcus, and Cyanobium diversity. This allowed us to unravel drastic community shifts over small to medium scale gradients of environmental factors, in particular along European coasts. The combined analysis of the phylogeography of natural populations and the thermophysiological characterization of eight strains, representative of the four major Synechococcus lineages (clades I to IV), also brought novel insights about the differential niche partitioning of clades I and IV, which most often co-dominate the Synechococcus community in cold and temperate coastal areas. Altogether, this study reveals several important characteristics and specificities of the coastal communities of Synechococcus worldwide. IMPORTANCE Synechococcus is the second most abundant phytoplanktonic organism on Earth, and its wide genetic diversity allowed it to colonize all the oceans except for polar waters, with different clades colonizing distinct oceanic niches. In recent years, the use of global metagenomics data sets has greatly improved our knowledge of "who is where" by describing the distribution of Synechococcus clades or ecotypes in the open ocean. However, little is known about the global distribution of Synechococcus ecotypes in coastal areas, where Synechococcus is often the dominant phytoplanktonic organism. Here, we leverage the global Ocean Sampling Day metagenomics data set to describe Synechococcus community composition in coastal areas worldwide, revealing striking community shifts, in particular along the coasts of Europe. As temperature appears as an important driver of the community composition, we also characterize the thermal preferenda of 8 Synechococcus strains, bringing new insights into the adaptation to temperature of the dominant Synechococcus clades.

Correction: McLennan et al. Assessing the Use of Molecular Barcoding and qPCR for Investigating the Ecology of Prorocentrum minimum (Dinophyceae), a Harmful Algal Species. Microorganisms 2021, 9, 510.

The authors wish to make the following corrections to this paper [...].

Chemotaxis may assist marine heterotrophic bacterial diazotrophs to find microzones suitable for N2 fixation in the pelagic ocean.

Heterotrophic bacterial diazotrophs (HBDs) are ubiquitous in the pelagic ocean, where they have been predicted to carry out the anaerobic process of nitrogen fixation within low-oxygen microenvironments associated with marine pelagic particles. However, the mechanisms enabling particle colonization by HBDs are unknown. We hypothesized that HBDs use chemotaxis to locate and colonize suitable microenvironments, and showed that a cultivated marine HBD is chemotactic toward amino acids and phytoplankton-derived DOM. Using an in situ chemotaxis assay, we also discovered that diverse HBDs at a coastal site are motile and chemotactic toward DOM from various phytoplankton taxa and, indeed, that the proportion of diazotrophs was up to seven times higher among the motile fraction of the bacterial community compared to the bulk seawater community. Finally, three of four HBD isolates and 16 of 17 HBD metagenome assembled genomes, recovered from major ocean basins and locations along the Australian coast, each encoded >85% of proteins affiliated with the bacterial chemotaxis pathway. These results document the widespread capacity for chemotaxis in diverse and globally relevant marine HBDs. We suggest that HBDs could use chemotaxis to seek out and colonize low-oxygen microenvironments suitable for nitrogen fixation, such as those formed on marine particles. Chemotaxis in HBDs could therefore affect marine nitrogen and carbon biogeochemistry by facilitating nitrogen fixation within otherwise oxic waters, while also altering particle degradation and the efficiency of the biological pump.

The Microbiological Drivers of Temporally Dynamic Dimethylsulfoniopropionate Cycling Processes in Australian Coastal Shelf Waters.

The organic sulfur compounds dimethylsulfoniopropionate (DMSP) and dimethyl sulfoxide (DMSO) play major roles in the marine microbial food web and have substantial climatic importance as sources and sinks of dimethyl sulfide (DMS). Seasonal shifts in the abundance and diversity of the phytoplankton and bacteria that cycle DMSP are likely to impact marine DMS (O) (P) concentrations, but the dynamic nature of these microbial interactions is still poorly resolved. Here, we examined the relationships between microbial community dynamics with DMS (O) (P) concentrations during a 2-year oceanographic time series conducted on the east Australian coast. Heterogenous temporal patterns were apparent in chlorophyll a (chl a) and DMSP concentrations, but the relationship between these parameters varied over time, suggesting the phytoplankton and bacterial community composition were affecting the net DMSP concentrations through differential DMSP production and degradation. Significant increases in DMSP were regularly measured in spring blooms dominated by predicted high DMSP-producing lineages of phytoplankton (Heterocapsa, Prorocentrum, Alexandrium, and Micromonas), while spring blooms that were dominated by predicted low DMSP-producing phytoplankton (Thalassiosira) demonstrated negligible increases in DMSP concentrations. During elevated DMSP concentrations, a significant increase in the relative abundance of the key copiotrophic bacterial lineage Rhodobacterales was accompanied by a three-fold increase in the gene, encoding the first step of DMSP demethylation (dmdA). Significant temporal shifts in DMS concentrations were measured and were significantly correlated with both fractions (0.2-2 µm and >2 µm) of microbial DMSP lyase activity. Seasonal increases of the bacterial DMSP biosynthesis gene (dsyB) and the bacterial DMS oxidation gene (tmm) occurred during the spring-summer and coincided with peaks in DMSP and DMSO concentration, respectively. These findings, along with significant positive relationships between dsyB gene abundance and DMSP, and tmm gene abundance with DMSO, reinforce the significant role planktonic bacteria play in producing DMSP and DMSO in ocean surface waters. Our results highlight the highly dynamic nature and myriad of microbial interactions that govern sulfur cycling in coastal shelf waters and further underpin the importance of microbial ecology in mediating important marine biogeochemical processes.

Chemotaxis shapes the microscale organization of the ocean's microbiome.

The capacity of planktonic marine microorganisms to actively seek out and exploit microscale chemical hotspots has been widely theorized to affect ocean-basin scale biogeochemistry1-3, but has never been examined comprehensively in situ among natural microbial communities. Here, using a field-based microfluidic platform to quantify the behavioural responses of marine bacteria and archaea, we observed significant levels of chemotaxis towards microscale hotspots of phytoplankton-derived dissolved organic matter (DOM) at a coastal field site across multiple deployments, spanning several months. Microscale metagenomics revealed that a wide diversity of marine prokaryotes, spanning 27 bacterial and 2 archaeal phyla, displayed chemotaxis towards microscale patches of DOM derived from ten globally distributed phytoplankton species. The distinct DOM composition of each phytoplankton species attracted phylogenetically and functionally discrete populations of bacteria and archaea, with 54% of chemotactic prokaryotes displaying highly specific responses to the DOM derived from only one or two phytoplankton species. Prokaryotes exhibiting chemotaxis towards phytoplankton-derived compounds were significantly enriched in the capacity to transport and metabolize specific phytoplankton-derived chemicals, and displayed enrichment in functions conducive to symbiotic relationships, including genes involved in the production of siderophores, B vitamins and growth-promoting hormones. Our findings demonstrate that the swimming behaviour of natural prokaryotic assemblages is governed by specific chemical cues, which dictate important biogeochemical transformation processes and the establishment of ecological interactions that structure the base of the marine food web.

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.

Biogeography of Southern Ocean prokaryotes: a comparison of the Indian and Pacific sectors.

We investigated the Southern Ocean (SO) prokaryote community structure via zero-radius operational taxonomic unit (zOTU) libraries generated from 16S rRNA gene sequencing of 223 full water column profiles. Samples reveal the prokaryote diversity trend between discrete water masses across multiple depths and latitudes in Indian (71-99°E, summer) and Pacific (170-174°W, autumn-winter) sectors of the SO. At higher taxonomic levels (phylum-family) we observed water masses to harbour distinct communities across both sectors, but observed sectorial variations at lower taxonomic levels (genus-zOTU) and relative abundance shifts for key taxa such as Flavobacteria, SAR324/Marinimicrobia, Nitrosopumilus and Nitrosopelagicus at both epi- and bathy-abyssopelagic water masses. Common surface bacteria were abundant in several deep-water masses and vice-versa suggesting connectivity between surface and deep-water microbial assemblages. Bacteria from same-sector Antarctic Bottom Water samples showed patchy, high beta-diversity which did not correlate well with measured environmental parameters or geographical distance. Unconventional depth distribution patterns were observed for key archaeal groups: Crenarchaeota was found across all depths in the water column and persistent high relative abundances of common epipelagic archaeon Nitrosopelagicus was observed in deep-water masses. Our findings reveal substantial regional variability of SO prokaryote assemblages that we argue should be considered in wide-scale SO ecosystem microbial modelling.

Biogeographical and seasonal dynamics of the marine Roseobacter community and ecological links to DMSP-producing phytoplankton.

Ecological interactions between marine bacteria and phytoplankton play a pivotal role in governing the ocean's major biogeochemical cycles. Among these, members of the marine Roseobacter Group (MRG) can establish mutualistic relationships with phytoplankton that are, in part, maintained by exchanges of the organosulfur compound, dimethylsulfoniopropionate (DMSP). Yet most of what is known about these interactions has been derived from culture-based laboratory studies. To investigate temporal and spatial co-occurrence patterns between members of the MRG and DMSP-producing phytoplankton we analysed 16S and 18S rRNA gene amplicon sequence variants (ASVs) derived from 5 years of monthly samples from seven environmentally distinct Australian oceanographic time-series. The MRG and DMSP-producer communities often displayed contemporaneous seasonality, which was greater in subtropical and temperate environments compared to tropical environments. The relative abundance of both groups varied latitudinally, displaying a poleward increase, peaking (MRG at 33% of total bacteria, DMSP producers at 42% of eukaryotic phototrophs) during recurrent spring-summer phytoplankton blooms in the most temperate site (Maria Island, Tasmania). Network analysis identified 20,140 significant positive correlations between MRG ASVs and DMSP producers and revealed that MRGs exhibit significantly stronger correlations to high DMSP producers relative to other DMSP-degrading bacteria (Pelagibacter, SAR86 and Actinobacteria). By utilising the power of a continental network of oceanographic time-series, this study provides in situ confirmation of interactions found in laboratory studies and demonstrates that the ecological dynamics of an important group of marine bacteria are shaped by the production of an abundant and biogeochemically significant organosulfur compound.

Functional characterisation of substrate-binding proteins to address nutrient uptake in marine picocyanobacteria.

Marine cyanobacteria are key primary producers, contributing significantly to the microbial food web and biogeochemical cycles by releasing and importing many essential nutrients cycled through the environment. A subgroup of these, the picocyanobacteria (Synechococcus and Prochlorococcus), have colonised almost all marine ecosystems, covering a range of distinct light and temperature conditions, and nutrient profiles. The intra-clade diversities displayed by this monophyletic branch of cyanobacteria is indicative of their success across a broad range of environments. Part of this diversity is due to nutrient acquisition mechanisms, such as the use of high-affinity ATP-binding cassette (ABC) transporters to competitively acquire nutrients, particularly in oligotrophic (nutrient scarce) marine environments. The specificity of nutrient uptake in ABC transporters is primarily determined by the peripheral substrate-binding protein (SBP), a receptor protein that mediates ligand recognition and initiates translocation into the cell. The recent availability of large numbers of sequenced picocyanobacterial genomes indicates both Synechococcus and Prochlorococcus apportion >50% of their transport capacity to ABC transport systems. However, the low degree of sequence homology among the SBP family limits the reliability of functional assignments using sequence annotation and prediction tools. This review highlights the use of known SBP structural representatives for the uptake of key nutrient classes by cyanobacteria to compare with predicted SBP functionalities within sequenced marine picocyanobacteria genomes. This review shows the broad range of conserved biochemical functions of picocyanobacteria and the range of novel and hypothetical ABC transport systems that require further functional characterisation.

Uptake of Dimethylsulfoniopropionate (DMSP) by Natural Microbial Communities of the Great Barrier Reef (GBR), Australia.

Dimethylsulfoniopropionate (DMSP) is a key organic sulfur compound that is produced by many phytoplankton and macrophytes and is ubiquitous in marine environments. Following its release into the water column, DMSP is primarily metabolised by heterotrophic bacterioplankton, but recent evidence indicates that non-DMSP producing phytoplankton can also assimilate DMSP from the surrounding environment. In this study, we examined the uptake of DMSP by communities of bacteria and phytoplankton within the waters of the Great Barrier Reef (GBR), Australia. We incubated natural GBR seawater with DMSP and quantified the uptake of DMSP by different fractions of the microbial community (>8 µm, 3-8 µm, <3 µm). We also evaluated how microbial community composition and the abundances of DMSP degrading genes are influenced by elevated dissolved DMSP levels. Our results showed uptake and accumulation of DMSP in all size fractions of the microbial community, with the largest fraction (>8 µm) forming the dominant sink, increasing in particulate DMSP by 44-115% upon DMSP enrichment. Longer-term incubations showed however, that DMSP retention was short lived (<24 h) and microbial responses to DMSP enrichment differed depending on the community carbon and sulfur demand. The response of the microbial communities from inside the reef indicated a preference towards cleaving DMSP into the climatically active aerosol dimethyl sulfide (DMS), whereas communities from the outer reef were sulfur and carbon limited, resulting in more DMSP being utilised by the cells. Our results show that DMSP uptake is shared across members of the microbial community, highlighting larger phytoplankton taxa as potentially relevant DMSP reservoirs and provide new information on sulfur cycling as a function of community metabolism in deeper, oligotrophic GBR waters.

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.

Assessing the Use of Molecular Barcoding and qPCR for Investigating the Ecology of Prorocentrum minimum (Dinophyceae), a Harmful Algal Species.

Prorocentrum minimum is a species of marine dinoflagellate that occurs worldwide and can be responsible for harmful algal blooms (HABs). Some studies have reported it to produce tetrodotoxin; however, results have been inconsistent. qPCR and molecular barcoding (amplicon sequencing) using high-throughput sequencing have been increasingly applied to quantify HAB species for ecological analyses and monitoring. Here, we isolated a strain of P. minimum from eastern Australian waters, where it commonly occurs, and developed and validated a qPCR assay for this species based on a region of ITS rRNA in relation to abundance estimates from the cultured strain as determined using light microscopy. We used this tool to quantify and examine ecological drivers of P. minimum in Botany Bay, an estuary in southeast Australia, for over ~14 months in 2016-2017. We compared abundance estimates using qPCR with those obtained using molecular barcoding based on an 18S rRNA amplicon. There was a significant correlation between the abundance estimates from amplicon sequencing and qPCR, but the estimates from light microscopy were not significantly correlated, likely due to the counting method applied. Using amplicon sequencing, ~600 unique actual sequence variants (ASVs) were found, much larger than the known phytoplankton diversity from this region. P. minimum abundance in Botany Bay was found to be significantly associated with lower salinities and higher dissolved CO2 levels.

Cyanorak v2.1: a scalable information system dedicated to the visualization and expert curation of marine and brackish picocyanobacteria genomes.

Cyanorak v2.1 (http://www.sb-roscoff.fr/cyanorak) is an information system dedicated to visualizing, comparing and curating the genomes of Prochlorococcus, Synechococcus and Cyanobium, the most abundant photosynthetic microorganisms on Earth. The database encompasses sequences from 97 genomes, covering most of the wide genetic diversity known so far within these groups, and which were split into 25,834 clusters of likely orthologous groups (CLOGs). The user interface gives access to genomic characteristics, accession numbers as well as an interactive map showing strain isolation sites. The main entry to the database is through search for a term (gene name, product, etc.), resulting in a list of CLOGs and individual genes. Each CLOG benefits from a rich functional annotation including EggNOG, EC/K numbers, GO terms, TIGR Roles, custom-designed Cyanorak Roles as well as several protein motif predictions. Cyanorak also displays a phyletic profile, indicating the genotype and pigment type for each CLOG, and a genome viewer (Jbrowse) to visualize additional data on each genome such as predicted operons, genomic islands or transcriptomic data, when available. This information system also includes a BLAST search tool, comparative genomic context as well as various data export options. Altogether, Cyanorak v2.1 constitutes an invaluable, scalable tool for comparative genomics of ecologically relevant marine microorganisms.

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.

Microbial tropicalization driven by a strengthening western ocean boundary current.

Western boundary currents (WBCs) redistribute heat and oligotrophic seawater from the tropics to temperate latitudes, with several displaying substantial climate change-driven intensification over the last century. Strengthening WBCs have been implicated in the poleward range expansion of marine macroflora and fauna, however, the impacts on the structure and function of temperate microbial communities are largely unknown. Here we show that the major subtropical WBC of the South Pacific Ocean, the East Australian Current (EAC), transports microbial assemblages that maintain tropical and oligotrophic (k-strategist) signatures, to seasonally displace more copiotrophic (r-strategist) temperate microbial populations within temperate latitudes of the Tasman Sea. We identified specific characteristics of EAC microbial assemblages compared with non-EAC assemblages, including strain transitions within the SAR11 clade, enrichment of Prochlorococcus, predicted smaller genome sizes and shifts in the importance of several functional genes, including those associated with cyanobacterial photosynthesis, secondary metabolism and fatty acid and lipid transport. At a temperate time-series site in the Tasman Sea, we observed significant reductions in standing stocks of total carbon and chlorophyll a, and a shift towards smaller phytoplankton and carnivorous copepods, associated with the seasonal impact of the EAC microbial assemblage. In light of the substantial shifts in microbial assemblage structure and function associated with the EAC, we conclude that climate-driven expansions of WBCs will expand the range of tropical oligotrophic microbes, and potentially profoundly impact the trophic status of temperate waters.

Investigating the Diversity of Marine Bacteriophage in Contrasting Water Masses Associated with the East Australian Current (EAC) System.

Virus- and bacteriophage-induced mortality can have a significant impact on marine productivity and alter the flux of nutrients in marine microbial food-webs. Viral mediated horizontal gene transfer can also influence host fitness and community composition. However, there are very few studies of marine viral diversity in the Southern Hemisphere, which hampers our ability to fully understand the complex interplay of biotic and abiotic factors that shape microbial communities. We carried out the first genetic study of bacteriophage communities within a dynamic western boundary current (WBC) system, the east Australian current (EAC). Virus DNA sequences were extracted from 63 assembled metagenomes and six metaviromes obtained from various depths at 24 different locations. More than 1700 bacteriophage genomic fragments (>9 kbps) were recovered from the assembled sequences. Bacteriophage diversity displayed distinct depth and regional patterns. There were clear differences in the bacteriophage populations associated with the EAC and Tasman Sea euphotic zones, at both the taxonomic and functional level. In contrast, bathypelagic phages were similar across the two oceanic regions. These data provide the first characterisation of viral diversity across a dynamic western boundary current, which is an emerging model for studying the response of microbial communities to climate change.