Single cell genomes of Prochlorococcus, Synechococcus, and sympatric microbes from diverse marine environments

A redox switch allows binding of Fe(II) and Fe(III) ions in the cyanobacterial iron-binding protein FutA from Prochlorococcus

The marine cyanobacterium Prochlorococcus is a main contributor to global photosynthesis, whilst being limited by iron availability. Cyanobacterial genomes generally encode two different types of FutA iron-binding proteins: periplasmic FutA2 ABC transporter subunits bind Fe(III), while cytosolic FutA1 binds Fe(II). Owing to their small size and their economized genome Prochlorococcus ecotypes typically possess a single futA gene. How the encoded FutA protein might bind different Fe oxidation states was previously unknown. Here, we use structural biology techniques at room temperature to probe the dynamic behavior of FutA. Neutron diffraction confirmed four negatively charged tyrosinates, that together with a neutral water molecule coordinate iron in trigonal bipyramidal geometry. Positioning of the positively charged Arg103 side chain in the second coordination shell yields an overall charge-neutral Fe(III) binding state in structures determined by neutron diffraction and serial femtosecond crystallography. Conventional rotation X-ray crystallography using a home source revealed X-ray-induced photoreduction of the iron center with observation of the Fe(II) binding state; here, an additional positioning of the Arg203 side chain in the second coordination shell maintained an overall charge neutral Fe(II) binding site. Dose series using serial synchrotron crystallography and an XFEL X-ray pump-probe approach capture the transition between Fe(III) and Fe(II) states, revealing how Arg203 operates as a switch to accommodate the different iron oxidation states. This switching ability of the Prochlorococcus FutA protein may reflect ecological adaptation by genome streamlining and loss of specialized FutA proteins.

Increased cyanophage infection at the bottom of the euphotic zone, especially in the fall

Picocyanobacteria contribute greatly to offshore primary production with cells extending through the deep euphotic zone. Literature indicates high viral infection of cyanobacteria in ocean transition zones. We postulate that the bottom of the euphotic zone is a transition zone, where communities transition from phototrophic to aphotic processes. We use single-copy core genes to examine cyanophage to cyanobacteria ratios in cellular metagenomes in the subtropical North Atlantic and Pacific. Cyanophage to cyanobacteria terL/rpoB ratios generally increase to >10 in the deep euphotic zone. As light levels decrease in the fall, Prochlorococcus in the deep euphotic zone experience reduced light levels. We find clear differences between spring (Geotraces GA02) and fall (GA03) in the North Atlantic, with terL/rpoB ratios increasing to >40 in the fall. When examining 23 months of the North Pacific Hawaii Ocean Timeseries, the depth of elevated cyanophage to cyanobacteria ratios in cellular metagenomes negatively correlated with surface photosynthetic radiation (PAR), particularly with the change in PAR, which reflected the season. In fall, all picocyanobacteria ecotypes were found at depths enriched with viruses, while in summer, only low light ecotypes were affected. Thus, we find high cyanophage infection both in the deep euphotic zone and during seasonal transitions.

Metagenome-assembled genomes reveal greatly expanded taxonomic and functional diversification of the abundant marine Roseobacter RCA cluster

BACKGROUND

The RCA (Roseobacter clade affiliated) cluster belongs to the family Roseobacteracea and represents a major Roseobacter lineage in temperate to polar oceans. Despite its prevalence and abundance, only a few genomes and one described species, Planktomarina temperata, exist. To gain more insights into our limited understanding of this cluster and its taxonomic and functional diversity and biogeography, we screened metagenomic datasets from the global oceans and reconstructed metagenome-assembled genomes (MAG) affiliated to this cluster.

RESULTS

The total of 82 MAGs, plus five genomes of isolates, reveal an unexpected diversity and novel insights into the genomic features, the functional diversity, and greatly refined biogeographic patterns of the RCA cluster. This cluster is subdivided into three genera: Planktomarina, Pseudoplanktomarina, and the most deeply branching Candidatus Paraplanktomarina. Six of the eight Planktomarina species have larger genome sizes (2.44-3.12 Mbp) and higher G + C contents (46.36-53.70%) than the four Pseudoplanktomarina species (2.26-2.72 Mbp, 42.22-43.72 G + C%). Cand. Paraplanktomarina is represented only by one species with a genome size of 2.40 Mbp and a G + C content of 45.85%. Three novel species of the genera Planktomarina and Pseudoplanktomarina are validly described according to the SeqCode nomenclature for prokaryotic genomes. Aerobic anoxygenic photosynthesis (AAP) is encoded in three Planktomarina species. Unexpectedly, proteorhodopsin (PR) is encoded in the other Planktomarina and all Pseudoplanktomarina species, suggesting that this light-driven proton pump is the most important mode of acquiring complementary energy of the RCA cluster. The Pseudoplanktomarina species exhibit differences in functional traits compared to Planktomarina species and adaptations to more resource-limited conditions. An assessment of the global biogeography of the different species greatly expands the range of occurrence and shows that the different species exhibit distinct biogeographic patterns. They partially reflect the genomic features of the species.

CONCLUSIONS

Our detailed MAG-based analyses shed new light on the diversification, environmental adaptation, and global biogeography of a major lineage of pelagic bacteria. The taxonomic delineation and validation by the SeqCode nomenclature of prominent genera and species of the RCA cluster may be a promising way for a refined taxonomic identification of major prokaryotic lineages and sublineages in marine and other prokaryotic communities assessed by metagenomics approaches. Video Abstract.

Synechococcus nitrogen gene loss in iron-limited ocean regions

Synechococcus are the most abundant cyanobacteria in high latitude regions and are responsible for an estimated 17% of annual marine net primary productivity. Despite their biogeochemical importance, Synechococcus populations have been unevenly sampled across the ocean, with most studies focused on low-latitude strains. In particular, the near absence of Synechococcus genomes from high-latitude, High Nutrient Low Chlorophyll (HNLC) regions leaves a gap in our knowledge of picocyanobacterial adaptations to iron limitation and their influence on carbon, nitrogen, and iron cycles. We examined Synechococcus populations from the subarctic North Pacific, a well-characterized HNLC region, with quantitative metagenomics. Assembly with short and long reads produced two near complete Synechococcus metagenome-assembled genomes (MAGs). Quantitative metagenome-derived abundances of these populations matched well with flow cytometry counts, and the Synechococcus MAGs were estimated to comprise >99% of the Synechococcus at Station P. Whereas the Station P Synechococcus MAGs contained multiple genes for adaptation to iron limitation, both genomes lacked genes for uptake and assimilation of nitrate and nitrite, suggesting a dependence on ammonium, urea, and other forms of recycled nitrogen leading to reduced iron requirements. A global analysis of Synechococcus nitrate reductase abundance in the TARA Oceans dataset found nitrate assimilation genes are also lower in other HNLC regions. We propose that nitrate and nitrite assimilation gene loss in Synechococcus may represent an adaptation to severe iron limitation in high-latitude regions where ammonium availability is higher. Our findings have implications for models that quantify the contribution of cyanobacteria to primary production and subsequent carbon export.

A compendium of bacterial and archaeal single-cell amplified genomes from oxygen deficient marine waters

Oxygen-deficient marine waters referred to as oxygen minimum zones (OMZs) or anoxic marine zones (AMZs) are common oceanographic features. They host both cosmopolitan and endemic microorganisms adapted to low oxygen conditions. Microbial metabolic interactions within OMZs and AMZs drive coupled biogeochemical cycles resulting in nitrogen loss and climate active trace gas production and consumption. Global warming is causing oxygen-deficient waters to expand and intensify. Therefore, studies focused on microbial communities inhabiting oxygen-deficient regions are necessary to both monitor and model the impacts of climate change on marine ecosystem functions and services. Here we present a compendium of 5,129 single-cell amplified genomes (SAGs) from marine environments encompassing representative OMZ and AMZ geochemical profiles. Of these, 3,570 SAGs have been sequenced to different levels of completion, providing a strain-resolved perspective on the genomic content and potential metabolic interactions within OMZ and AMZ microbiomes. Hierarchical clustering confirmed that samples from similar oxygen concentrations and geographic regions also had analogous taxonomic compositions, providing a coherent framework for comparative community analysis.

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.

Regulatory and metabolic adaptations in the nitrogen assimilation of marine picocyanobacteria

Prochlorococcus and Synechococcus are the two most abundant photosynthetic organisms on Earth, with a strong influence on the biogeochemical carbon and nitrogen cycles. Early reports demonstrated the streamlining of regulatory mechanisms in nitrogen metabolism and the removal of genes not strictly essential. The availability of a large series of genomes, and the utilization of latest generation molecular techniques have allowed elucidating the main mechanisms developed by marine picocyanobacteria to adapt to the environments where they thrive, with a particular interest in the strains inhabiting oligotrophic oceans. Given that nitrogen is often limited in those environments, a series of studies have explored the strategies utilized by Prochlorococcus and Synechococcus to exploit the low concentrations of nitrogen-containing molecules available in large areas of the oceans. These strategies include the reduction in the GC and the cellular protein contents; the utilization of truncated proteins; a reduced average amount of N in the proteome; the development of metabolic mechanisms to perceive and utilize nanomolar nitrate concentrations; and the reduced responsiveness of key molecular regulatory systems such as NtcA to 2-oxoglutarate. These findings are in sharp contrast with the large body of knowledge obtained in freshwater cyanobacteria. We will outline the main discoveries, stressing their relevance to the ecological success of these important microorganisms.

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.

Synthetic biology in marine cyanobacteria: Advances and challenges

The current economic and environmental context requests an accelerating development of sustainable alternatives for the production of various target compounds. Biological processes offer viable solutions and have gained renewed interest in the recent years. For example, photosynthetic chassis organisms are particularly promising for bioprocesses, as they do not require biomass-derived carbon sources and contribute to atmospheric CO2 fixation, therefore supporting climate change mitigation. Marine cyanobacteria are of particular interest for biotechnology applications, thanks to their rich diversity, their robustness to environmental changes, and their metabolic capabilities with potential for therapeutics and chemicals production without requiring freshwater. The additional cyanobacterial properties, such as efficient photosynthesis, are also highly beneficial for biotechnological processes. Due to their capabilities, research efforts have developed several genetic tools for direct metabolic engineering applications. While progress toward a robust genetic toolkit is continuously achieved, further work is still needed to routinely modify these species and unlock their full potential for industrial applications. In contrast to the understudied marine cyanobacteria, genetic engineering and synthetic biology in freshwater cyanobacteria are currently more advanced with a variety of tools already optimized. This mini-review will explore the opportunities provided by marine cyanobacteria for a greener future. A short discussion will cover the advances and challenges regarding genetic engineering and synthetic biology in marine cyanobacteria, followed by a parallel with freshwater cyanobacteria and their current genetic availability to guide the prospect for marine species.

Bacterioplankton communities reveal horizontal and vertical influence of an Island Mass Effect

Coral reefs are highly productive ecosystems with distinct biogeochemistry and biology nestled within unproductive oligotrophic gyres. Coral reef islands have often been associated with a nearshore enhancement in phytoplankton, a phenomenon known as the Island Mass Effect (IME). Despite being documented more than 60 years ago, much remains unknown about the extent and drivers of IMEs. Here we utilized 16S rRNA gene metabarcoding as a biological tracer to elucidate horizontal and vertical influence of an IME around the islands of Mo'orea and Tahiti, French Polynesia. We show that those nearshore oceanic stations with elevated chlorophyll a included bacterioplankton found in high abundance in the reef environment, suggesting advection of reef water is the source of altered nearshore biogeochemistry. We also observed communities in the nearshore deep chlorophyll maximum (DCM) with enhanced abundances of upper euphotic bacterioplankton that correlated with intrusions of low-density, O₂ rich water, suggesting island influence extends into the DCM.

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.

New cyanobacterial genus Argonema is hidding in soil crusts around the world

Cyanobacteria are crucial primary producers in soil and soil crusts. However, their biodiversity in these habitats remains poorly understood, especially in the tropical and polar regions. We employed whole genome sequencing, morphology, and ecology to describe a novel cyanobacterial genus Argonema isolated from Antarctica. Extreme environments are renowned for their relatively high number of endemic species, but whether cyanobacteria are endemic or not is open to much current debate. To determine if a cyanobacterial lineage is endemic is a time consuming, elaborate, and expensive global sampling effort. Thus, we propose an approach that will help to overcome the limits of the sampling effort and better understand the global distribution of cyanobacterial clades. We employed a Sequencing Read Archive, which provides a rich source of data from thousands of environmental samples. We developed a framework for a characterization of the global distribution of any microbial species using Sequencing Read Archive. Using this approach, we found that Argonema is actually cosmopolitan in arid regions. It provides further evidence that endemic microbial taxa are likely to be much rarer than expected.

Picocyanobacterial Synechococcus in marine ecosystem: Insights from genetic diversity, global distribution, and potential function

Marine Synechococcus, a main group of picocyanobacteria, has been ubiquitously observed across the global oceans. Synechococcus exhibits high phylogenetical and phenotypical diversity, and horizontal gene transfer makes its genetic evolution much more intricate. With the development of measurement technologies and analysis methods, the genomic information and niche partition of each Synechococcus lineage tend to be precisely described, but the global analysis is still lacking. Therefore, it is necessary to summarize existing studies and integrate published data to gain a comprehensive understanding of Synechococcus on genetic variation, niche division, and potential functions. In this review, the maximum likelihood trees are constructed based on existing sequence data, including both phylogenetic and pigmentary gene markers. The global distribution characteristics of abundance, lineages, and pigment types are concluded through pooled analysis of more than 700 samples obtained from approximately 50 scientific research cruises. The potential functions of Synechococcus are explored in element cycles and biological interactions. Future work on Synechococcus is suggested to focus on not only elucidating the nature of Synechococcus biodiversity but also demonstrating its interactions with the ecosystem by combining bioinformatics and macroscopic isotope-labeled environmental parameters.

Toward quantifying the adaptive role of bacterial pangenomes during environmental perturbations

Metagenomic surveys have revealed that natural microbial communities are predominantly composed of sequence-discrete, species-like populations but the genetic and/or ecological processes that maintain such populations remain speculative, limiting our understanding of population speciation and adaptation to perturbations. To address this knowledge gap, we sequenced 112 Salinibacter ruber isolates and 12 companion metagenomes from four adjacent saltern ponds in Mallorca, Spain that were experimentally manipulated to dramatically alter salinity and light intensity, the two major drivers of this ecosystem. Our analyses showed that the pangenome of the local Sal. ruber population is open and similar in size (~15,000 genes) to that of randomly sampled Escherichia coli genomes. While most of the accessory (noncore) genes were isolate-specific and showed low in situ abundances based on the metagenomes compared to the core genes, indicating that they were functionally unimportant and/or transient, 3.5% of them became abundant when salinity (but not light) conditions changed and encoded for functions related to osmoregulation. Nonetheless, the ecological advantage of these genes, while significant, was apparently not strong enough to purge diversity within the population. Collectively, our results provide an explanation for how this immense intrapopulation gene diversity is maintained, which has implications for the prokaryotic species concept.

Viruses affect picocyanobacterial abundance and biogeography in the North Pacific Ocean

The photosynthetic picocyanobacteria Prochlorococcus and Synechococcus are models for dissecting how ecological niches are defined by environmental conditions, but how interactions with bacteriophages affect picocyanobacterial biogeography in open ocean biomes has rarely been assessed. We applied single-virus and single-cell infection approaches to quantify cyanophage abundance and infected picocyanobacteria in 87 surface water samples from five transects that traversed approximately 2,200 km in the North Pacific Ocean on three cruises, with a duration of 2-4 weeks, between 2015 and 2017. We detected a 550-km-wide hotspot of cyanophages and virus-infected picocyanobacteria in the transition zone between the North Pacific Subtropical and Subpolar gyres that was present in each transect. Notably, the hotspot occurred at a consistent temperature and displayed distinct cyanophage-lineage composition on all transects. On two of these transects, the levels of infection in the hotspot were estimated to be sufficient to substantially limit the geographical range of Prochlorococcus. Coincident with the detection of high levels of virally infected picocyanobacteria, we measured an increase of 10-100-fold in the Synechococcus populations in samples that are usually dominated by Prochlorococcus. We developed a multiple regression model of cyanophages, temperature and chlorophyll concentrations that inferred that the hotspot extended across the North Pacific Ocean, creating a biological boundary between gyres, with the potential to release organic matter comparable to that of the sevenfold-larger North Pacific Subtropical Gyre. Our results highlight the probable impact of viruses on large-scale phytoplankton biogeography and biogeochemistry in distinct regions of the oceans.

A Novel and Ubiquitous Marine Methylophage Provides Insights into Viral-Host Coevolution and Possible Host-Range Expansion in Streamlined Marine Heterotrophic Bacteria

The methylotrophic OM43 clade are Gammaproteobacteria that comprise some of the smallest free-living cells known and have highly streamlined genomes. OM43 represents an important microbial link between marine primary production and remineralization of carbon back to the atmosphere. Bacteriophages shape microbial communities and are major drivers of mortality and global marine biogeochemistry. Recent cultivation efforts have brought the first viruses infecting members of the OM43 clade into culture. Here, we characterize a novel myophage infecting OM43 called Melnitz. Melnitz was isolated independently from water samples from a subtropical ocean gyre (Sargasso Sea) and temperate coastal (Western English Channel) systems. Metagenomic recruitment from global ocean viromes confirmed that Melnitz is globally ubiquitous, congruent with patterns of host abundance. Bacteria with streamlined genomes such as OM43 and the globally dominant SAR11 clade use riboswitches as an efficient method to regulate metabolism. Melnitz encodes a two-piece tmRNA (ssrA), controlled by a glutamine riboswitch, providing evidence that riboswitch use also occurs for regulation during phage infection of streamlined heterotrophs. Virally encoded tRNAs and ssrA found in Melnitz were phylogenetically more closely related to those found within the alphaproteobacterial SAR11 clade and their associated myophages than those within their gammaproteobacterial hosts. This suggests the possibility of an ancestral host transition event between SAR11 and OM43. Melnitz and a related myophage that infects SAR11 were unable to infect hosts of the SAR11 and OM43, respectively, suggesting host transition rather than a broadening of host range.

IMPORTANCE Isolation and cultivation of viruses are the foundations on which the mechanistic understanding of virus-host interactions and parameterization of bioinformatic tools for viral ecology are based. This study isolated and characterized the first myophage known to infect the OM43 clade, expanding our knowledge of this understudied group of microbes. The nearly identical genomes of four strains of Melnitz isolated from different marine provinces and the global abundance estimations from metagenomic data suggest that this viral population is globally ubiquitous. Genome analysis revealed several unusual features in Melnitz and related genomes recovered from viromes, such as a curli operon and virally encoded tmRNA controlled by a glutamine riboswitch, neither of which are found in the host. Further phylogenetic analysis of shared genes indicates that this group of viruses infecting the gammaproteobacterial OM43 shares a recent common ancestor with viruses infecting the abundant alphaproteobacterial SAR11 clade. Host ranges are affected by compatible cell surface receptors, successful circumvention of superinfection exclusion systems, and the presence of required accessory proteins, which typically limits phages to singular narrow groups of closely related bacterial hosts. This study provides intriguing evidence that for streamlined heterotrophic bacteria, virus-host transitioning may not be necessarily restricted to phylogenetically related hosts but is a function of shared physical and biochemical properties of the cell.

Diversity and evolution of pigment types in marine Synechococcus cyanobacteria

Synechococcus cyanobacteria are ubiquitous and abundant in the marine environment and contribute to an estimated 16% of the ocean net primary productivity. Their light-harvesting complexes, called phycobilisomes (PBS), are composed of a conserved allophycocyanin core, from which radiates six to eight rods with variable phycobiliprotein and chromophore content. This variability allows Synechococcus cells to optimally exploit the wide variety of spectral niches existing in marine ecosystems. Seven distinct pigment types or subtypes have been identified so far in this taxon based on the phycobiliprotein composition and/or the proportion of the different chromophores in PBS rods. Most genes involved in their biosynthesis and regulation are located in a dedicated genomic region called the PBS rod region. Here, we examine the variability of gene content and organization of this genomic region in a large set of sequenced isolates and natural populations of Synechococcus representative of all known pigment types. All regions start with a tRNA-Phe_GAA and some possess mobile elements for DNA integration and site-specific recombination, suggesting that their genomic variability relies in part on a “tycheposon”-like mechanism. Comparison of the phylogenies obtained for PBS and core genes revealed that the evolutionary history of PBS rod genes differs from the core genome and is characterized by the co-existence of different alleles and frequent allelic exchange. We propose a scenario for the evolution of the different pigment types and highlight the importance of incomplete lineage sorting in maintaining a wide diversity of pigment types in different Synechococcus lineages despite multiple speciation events.

Beyond Soil-Dwelling Actinobacteria: Fantastic Antibiotics and Where to Find Them

Bacterial secondary metabolites represent an invaluable source of bioactive molecules for the pharmaceutical and agrochemical industries. Although screening campaigns for the discovery of new compounds have traditionally been strongly biased towards the study of soil-dwelling Actinobacteria, the current antibiotic resistance and discovery crisis has brought a considerable amount of attention to the study of previously neglected bacterial sources of secondary metabolites. The development and application of new screening, sequencing, genetic manipulation, cultivation and bioinformatic techniques have revealed several other groups of bacteria as producers of striking chemical novelty. Biosynthetic machineries evolved from independent taxonomic origins and under completely different ecological requirements and selective pressures are responsible for these structural innovations. In this review, we summarize the most important discoveries related to secondary metabolites from alternative bacterial sources, trying to provide the reader with a broad perspective on how technical novelties have facilitated the access to the bacterial metabolic dark matter.

Genetic engineering of marine cyanophages reveals integration but not lysogeny in T7-like cyanophages

Marine cyanobacteria of the genera Synechococcus and Prochlorococcus are the most abundant photosynthetic organisms on earth, spanning vast regions of the oceans and contributing significantly to global primary production. Their viruses (cyanophages) greatly influence cyanobacterial ecology and evolution. Although many cyanophage genomes have been sequenced, insight into the functional role of cyanophage genes is limited by the lack of a cyanophage genetic engineering system. Here, we describe a simple, generalizable method for genetic engineering of cyanophages from multiple families, that we named REEP for REcombination, Enrichment and PCR screening. This method enables direct investigation of key cyanophage genes, and its simplicity makes it adaptable to other ecologically relevant host-virus systems. T7-like cyanophages often carry integrase genes and attachment sites, yet exhibit lytic infection dynamics. Here, using REEP, we investigated their ability to integrate and maintain a lysogenic life cycle. We found that these cyanophages integrate into the host genome and that the integrase and attachment site are required for integration. However, stable lysogens did not form. The frequency of integration was found to be low in both lab cultures and the oceans. These findings suggest that T7-like cyanophage integration is transient and is not part of a classical lysogenic cycle.

Basin-scale biogeography of marine phytoplankton reflects cellular-scale optimization of metabolism and physiology

Extensive microdiversity within Prochlorococcus, the most abundant marine cyanobacterium, occurs at scales from a single droplet of seawater to ocean basins. To interpret the structuring role of variations in genetic potential, as well as metabolic and physiological acclimation, we developed a mechanistic constraint-based modeling framework that incorporates the full suite of genes, proteins, metabolic reactions, pigments, and biochemical compositions of 69 sequenced isolates spanning the Prochlorococcus pangenome. Optimizing each strain to the local, observed physical and chemical environment along an Atlantic Ocean transect, we predicted variations in strain-specific patterns of growth rate, metabolic configuration, and physiological state, defining subtle niche subspaces directly attributable to differences in their encoded metabolic potential. Predicted growth rates covaried with observed ecotype abundances, affirming their significance as a measure of fitness and inferring a nonlinear density dependence of mortality. Our study demonstrates the potential to interpret global-scale ecosystem organization in terms of cellular-scale processes.

A Genomic Perspective Across Earth's Microbiomes Reveals That Genome Size in Archaea and Bacteria Is Linked to Ecosystem Type and Trophic Strategy

Our view of genome size in Archaea and Bacteria has remained skewed as the data has been dominated by genomes of microorganisms that have been cultivated under laboratory settings. However, the continuous effort to catalog Earth's microbiomes, specifically propelled by recent extensive work on uncultivated microorganisms, provides an opportunity to revise our perspective on genome size distribution. We present a meta-analysis that includes 26,101 representative genomes from 3 published genomic databases; metagenomic assembled genomes (MAGs) from GEMs and stratfreshDB, and isolates from GTDB. Aquatic and host-associated microbial genomes present on average the smallest estimated genome sizes (3.1 and 3.0 Mbp, respectively). These are followed by terrestrial microbial genomes (average 3.7 Mbp), and genomes from isolated microorganisms (average 4.3 Mbp). On the one hand, aquatic and host-associated ecosystems present smaller genomes sizes in genera of phyla with genome sizes above 3 Mbp. On the other hand, estimated genome size in phyla with genomes under 3 Mbp showed no difference between ecosystems. Moreover, we observed that when using 95% average nucleotide identity (ANI) as an estimator for genetic units, only 3% of MAGs cluster together with genomes from isolated microorganisms. Although there are potential methodological limitations when assembling and binning MAGs, we found that in genome clusters containing both environmental MAGs and isolate genomes, MAGs were estimated only an average 3.7% smaller than isolate genomes. Even when assembly and binning methods introduce biases, estimated genome size of MAGs and isolates are very similar. Finally, to better understand the ecological drivers of genome size, we discuss on the known and the overlooked factors that influence genome size in different ecosystems, phylogenetic groups, and trophic strategies.

MOSGA 2: Comparative genomics and validation tools

Due to the highly growing number of available genomic information, the need for accessible and easy-to-use analysis tools is increasing. To facilitate eukaryotic genome annotations, we created MOSGA. In this work, we show how MOSGA 2 is developed by including several advanced analyses for genomic data. Since the genomic data quality greatly impacts the annotation quality, we included multiple tools to validate and ensure high-quality user-submitted genome assemblies. Moreover, thanks to the integration of comparative genomics methods, users can benefit from a broader genomic view by analyzing multiple genomic data sets simultaneously. Further, we demonstrate the new functionalities of MOSGA 2 by different use-cases and practical examples. MOSGA 2 extends the already established application to the quality control of the genomic data and integrates and analyzes multiple genomes in a larger context, e.g., by phylogenetics.

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

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

Dysbiosis in marine aquaculture revealed through microbiome analysis: reverse ecology for environmental sustainability

The increasing demand for products for human consumption is leading to the fast-growing expansion of numerous food sectors such as marine aquaculture (mariculture). However, excessive input of nutrients and pollutants modifies marine ecosystems. Here, we applied a metagenomic approach to investigate these perturbations in samples from marine farms of gilthead seabream cultures. Results revealed dysbiosis and functional imbalance within the net cage with a unique structure, with little interference with samples from the fish microbiota or those collected far away from the coast. Remarkably, below the cage the prokaryotic community was highly similar to the marine microbiome of photic offshore samples. We recovered 48 novel metagenome-assembled genomes. Metagenomic recruitment revealed a significant change in the microbial community which was dominated by several Proteobacteria orders (Sphingomonadales, Pseudomonadales, Caudobacterales and Rhizobiales). Genomic potential for bioremediation processes, including nitrate removal through aerobic denitrification, and degradation of aromatic compounds and other toxic products were enriched in these microbes. The detrimental side effects were the increased number of antimicrobial resistance genes and the presence of potentially emergent pathogens. Knowledge of this metabolic diversity and the microbes involved in ecological balance recovery can be used to reduce the environmental impact of these practices.

Opening up the Single-Cell Toolbox for Microbial Natural Products Research

The diverse microbes that produce natural products represent an important source of novel therapeutics, drug leads, and scientific tools. However, the vast majority have not been grown in axenic culture and are members of complex communities. While meta-'omic methods such as metagenomics, -transcriptomics, and -proteomics reveal collective molecular features of this "microbial dark matter", the study of individual microbiome members can be challenging. To address these limits, a number of techniques with single-bacterial resolution have been developed in the last decade and a half. While several of these are embraced by microbial ecologists, there has been less use by researchers interested in mining microbes for natural products. In this review, we discuss the available and emerging techniques for targeted single-cell analysis with a particular focus on applications to the discovery and study of natural products.

Pili allow dominant marine cyanobacteria to avoid sinking and evade predation

How oligotrophic marine cyanobacteria position themselves in the water column is currently unknown. The current paradigm is that these organisms avoid sinking due to their reduced size and passive drift within currents. Here, we show that one in four picocyanobacteria encode a type IV pilus which allows these organisms to increase drag and remain suspended at optimal positions in the water column, as well as evade predation by grazers. The evolution of this sophisticated floatation mechanism in these purely planktonic streamlined microorganisms has important implications for our current understanding of microbial distribution in the oceans and predator-prey interactions which ultimately will need incorporating into future models of marine carbon flux dynamics.

Uncultivated Viral Populations Dominate Estuarine Viromes on the Spatiotemporal Scale

Viruses are ubiquitous and abundant in the oceans, and viral metagenomes (viromes) have been investigated extensively via several large-scale ocean sequencing projects. However, there have not been any systematic viromic studies in estuaries. Here, we investigated the viromes of the Delaware Bay and Chesapeake Bay, two Mid-Atlantic estuaries. Deep sequencing generated a total of 48,190 assembled viral sequences (>5 kb) and 26,487 viral populations (9,204 virus clusters and 17,845 singletons), including 319 circular viral contigs between 7.5 kb and 161.8 kb. Unknown viruses represented the vast majority of the dominant populations, while the composition of known viruses, such as pelagiphage and cyanophage, appeared to be relatively consistent across a wide range of salinity gradients and in different seasons. A difference between estuarine and ocean viromes was reflected by the proportions of Myoviridae, Podoviridae, Siphoviridae, Phycodnaviridae, and a few well-studied virus representatives. The difference in viral community between the Delaware Bay and Chesapeake Bay is significantly more pronounced than the difference caused by temperature or salinity, indicating strong local profiles caused by the unique ecology of each estuary. Interestingly, a viral contig similar to phages infecting Acinetobacter baumannii (“Iraqibacter”) was found to be highly abundant in the Delaware Bay but not in the Chesapeake Bay, the source of which is yet to be identified. Highly abundant viruses in both estuaries have close hits to viral sequences derived from the marine single-cell genomes or long-read single-molecule sequencing, suggesting that important viruses are still waiting to be discovered in the estuarine environment.

IMPORTANCE This is the first systematic study about spatial and temporal variation of virioplankton communities in estuaries using deep metagenomics sequencing. It is among the highest-quality viromic data sets to date, showing remarkably consistent sequencing depth and quality across samples. Our results indicate that there exists a large pool of abundant and diverse viruses in estuaries that have not yet been cultivated, their genomes only available thanks to single-cell genomics or single-molecule sequencing, demonstrating the importance of these methods for viral discovery. The spatiotemporal pattern of these abundant uncultivated viruses is more variable than that of cultured viruses. Despite strong environmental gradients, season and location had surprisingly little impact on the viral community within an estuary, but we saw a significant distinction between the two estuaries and also between estuarine and open ocean viromes.

Frequency of mispackaging of Prochlorococcus DNA by cyanophage

Prochlorococcus cells are the numerically dominant phototrophs in the open ocean. Cyanophages that infect them are a notable fraction of the total viral population in the euphotic zone, and, as vehicles of horizontal gene transfer, appear to drive their evolution. Here we examine the propensity of three cyanophages-a podovirus, a siphovirus, and a myovirus-to mispackage host DNA in their capsids while infecting Prochlorococcus, the first step in phage-mediated horizontal gene transfer. We find the mispackaging frequencies are distinctly different among the three phages. Myoviruses mispackage host DNA at low and seemingly fixed frequencies, while podo- and siphoviruses vary in their mispackaging frequencies by orders of magnitude depending on growth light intensity. We link this difference to the concentration of intracellular reactive oxygen species and protein synthesis rates, both parameters increasing in response to higher light intensity. Based on our findings, we propose a model of mispackaging frequency determined by the imbalance between the production of capsids and the number of phage genome copies during infection: when protein synthesis rate increase to levels that the phage cannot regulate, they lead to an accumulation of empty capsids, in turn triggering more frequent host DNA mispackaging errors.

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.

A single-cell polony method reveals low levels of infected Prochlorococcus in oligotrophic waters despite high cyanophage abundances

Long-term stability of picocyanobacteria in the open oceans is maintained by a balance between synchronous division and death on daily timescales. Viruses are considered a major source of microbial mortality, however, current methods to measure infection have significant methodological limitations. Here we describe a method that pairs flow-cytometric sorting with a PCR-based polony technique to simultaneously screen thousands of taxonomically resolved individual cells for intracellular virus DNA, enabling sensitive, high-throughput, and direct quantification of infection by different virus lineages. Under controlled conditions with picocyanobacteria-cyanophage models, the method detected infection throughout the lytic cycle and discriminated between varying infection levels. In North Pacific subtropical surface waters, the method revealed that only a small percentage of Prochlorococcus (0.35–1.6%) were infected, predominantly by T4-like cyanophages, and that infection oscillated 2-fold in phase with the diel cycle. This corresponds to 0.35–4.8% of Prochlorococcus mortality daily. Cyanophages were 2–4-fold more abundant than Prochlorococcus, indicating that most encounters did not result in infection and suggesting infection is mitigated via host resistance, reduced phage infectivity and inefficient adsorption. This method will enable quantification of infection for key microbial taxa across oceanic regimes and will help determine the extent that viruses shape microbial communities and ecosystem level processes.

GToTree: a user-friendly workflow for phylogenomics

Summary

Genome-level evolutionary inference (i.e. phylogenomics) is becoming an increasingly essential step in many biologists’ work. Accordingly, there are several tools available for the major steps in a phylogenomics workflow. But for the biologist whose main focus is not bioinformatics, much of the computational work required—such as accessing genomic data on large scales, integrating genomes from different file formats, performing required filtering, stitching different tools together etc.—can be prohibitive. Here I introduce GToTree, a command-line tool that can take any combination of fasta files, GenBank files and/or NCBI assembly accessions as input and outputs an alignment file, estimates of genome completeness and redundancy, and a phylogenomic tree based on a specified single-copy gene (SCG) set. Although GToTree can work with any custom hidden Markov Models (HMMs), also included are 13 newly generated SCG-set HMMs for different lineages and levels of resolution, built based on searches of ∼12 000 bacterial and archaeal high-quality genomes. GToTree aims to give more researchers the capability to make phylogenomic trees.

Availability and implementation

GToTree is open-source and freely available for download from: github.com/AstrobioMike/GToTree. It is implemented primarily in bash with helper scripts written in python.

Supplementary information

Supplementary data are available at Bioinformatics online.

Mixotrophy in marine picocyanobacteria: use of organic compounds by Prochlorococcus and Synechococcus

Marine picocyanobacteria of the Prochlorococcus and Synechococcus genera have been longtime considered as autotrophic organisms. However, compelling evidence published over the last 15 years shows that these organisms can use different organic compounds containing key elements to survive in oligotrophic oceans, such as N (amino acids, amino sugars), S (dimethylsulfoniopropionate, DMSP), or P (ATP). Furthermore, marine picocyanobacteria can also take up glucose and use it as a source of carbon and energy, despite the fact that this compound is devoid of limiting elements and can also be synthesized by using standard metabolic pathways. This review will outline the main findings suggesting mixotrophy in the marine picocyanobacteria Prochlorococcus and Synechococcus, and its ecological relevance for these important primary producers.

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.

GEOTRACES: Accelerating Research on the Marine Biogeochemical Cycles of Trace Elements and Their Isotopes

The biogeochemical cycles of trace elements and their isotopes (TEIs) constitute an active area of oceanographic research due to their role as essential nutrients for marine organisms and their use as tracers of oceanographic processes. Selected TEIs also provide diagnostic information about the physical, geological, and chemical processes that supply or remove solutes in the ocean. Many of these same TEIs provide information about ocean conditions in the past, as their imprint on marine sediments can be interpreted to reflect changes in ocean circulation, biological productivity, the ocean carbon cycle, and more. Other TEIs have been introduced as the result of human activities and are considered contaminants. The development and implementation of contamination-free methods for collecting and analyzing samples for TEIs revolutionized marine chemistry, revealing trace element distributions with oceanographically consistent features and new insights about the processes regulating them. Despite these advances, the volume and geographic coverage of high-quality TEI data by the end of the twentieth century were insufficient to constrain their global biogeochemical cycles. To accelerate progress in this field of research, marine geochemists developed a coordinated international effort to systematically study the marine biogeochemical cycles of TEIs-the GEOTRACES program. Following a decade of planning and implementation, GEOTRACES launched its main field effort in 2010. This review, roughly midway through the field program, summarizes the steps involved in designing the program, its management structure, and selected findings.

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.

Ecogenomics of the SAR11 clade

Members of the SAR11 clade, despite their high abundance, are often poorly represented by metagenome‐assembled genomes. This fact has hampered our knowledge about their ecology and genetic diversity. Here we examined 175 SAR11 genomes, including 47 new single‐amplified genomes. The presence of the first genomes associated with subclade IV suggests that, in the same way as subclade V, they might be outside the proposed Pelagibacterales order. An expanded phylogenomic classification together with patterns of metagenomic recruitment at a global scale have allowed us to define new ecogenomic units of classification (genomospecies), appearing at different, and sometimes restricted, metagenomic data sets. We detected greater microdiversity across the water column at a single location than in samples collected from similar depth across the global ocean, suggesting little influence of biogeography. In addition, pangenome analysis revealed that the flexible genome was essential to shape genomospecies distribution. In one genomospecies preferentially found within the Mediterranean, a set of genes involved in phosphonate utilization was detected. While another, with a more cosmopolitan distribution, was unique in having an aerobic purine degradation pathway. Together, these results provide a glimpse of the enormous genomic diversity within this clade at a finer resolution than the currently defined clades.

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

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

Charting the Complexity of the Marine Microbiome through Single-Cell Genomics

Marine bacteria and archaea play key roles in global biogeochemistry. To improve our understanding of this complex microbiome, we employed single-cell genomics and a randomized, hypothesis-agnostic cell selection strategy to recover 12,715 partial genomes from the tropical and subtropical euphotic ocean. A substantial fraction of known prokaryoplankton coding potential was recovered from a single, 0.4 mL ocean sample, which indicates that genomic information disperses effectively across the globe. Yet, we found each genome to be unique, implying limited clonality within prokaryoplankton populations. Light harvesting and secondary metabolite biosynthetic pathways were numerous across lineages, highlighting the value of single-cell genomics to advance the identification of ecological roles and biotechnology potential of uncultured microbial groups. This genome collection enabled functional annotation and genus-level taxonomic assignments for >80% of individual metagenome reads from the tropical and subtropical surface ocean, thus offering a model to improve reference genome databases for complex microbiomes.

Contributions of single-cell genomics to our understanding of planktonic marine archaea

Single-cell genomics has transformed many fields of biology, marine microbiology included. Here, we consider the impact of single-cell genomics on a specific group of marine microbes-the planktonic marine archaea. Despite single-cell enabled discoveries of novel metabolic function in the marine thaumarchaea, population-level investigations are hindered by an overall lower than expected recovery of thaumarchaea in single-cell studies. Metagenome-assembled genomes have so far been a more useful method for accessing genome-resolved insights into the Marine Group II euryarchaea. Future progress in the application of single-cell genomics to archaeal biology in the ocean would benefit from more targeted sorting approaches, and a more systematic investigation of potential biases against archaea in single-cell workflows including cell lysis, genome amplification and genome screening. This article is part of a discussion meeting issue 'Single cell ecology'.

Targeted metagenomic recovery of four divergent viruses reveals shared and distinctive characteristics of giant viruses of marine eukaryotes

Giant viruses have remarkable genomic repertoires—blurring the line with cellular life—and act as top–down controls of eukaryotic plankton. However, to date only six cultured giant virus genomes are available from the pelagic ocean. We used at-sea flow cytometry with staining and sorting designed to target wild predatory eukaryotes, followed by DNA sequencing and assembly, to recover novel giant viruses from the Pacific Ocean. We retrieved four ‘PacV’ partial genomes that range from 421 to 1605 Kb, with 13 contigs on average, including the largest marine viral genomic assembly reported to date. Phylogenetic analyses indicate that three of the new viruses span a clade with deep-branching members of giant Mimiviridae, incorporating the Cafeteria roenbergensis virus, the uncultivated terrestrial Faunusvirus, one PacV from a choanoflagellate and two PacV with unclear hosts. The fourth virus, oPacV-421, is phylogenetically related to viruses that infect haptophyte algae. About half the predicted proteins in each PacV have no matches in NCBI nr (e-value < 10⁻⁵), totalling 1735 previously unknown proteins; the closest affiliations of the other proteins were evenly distributed across eukaryotes, prokaryotes and viruses of eukaryotes. The PacVs encode many translational proteins and two encode eukaryotic-like proteins from the Rh family of the ammonium transporter superfamily, likely influencing the uptake of nitrogen during infection. cPacV-1605 encodes a microbial viral rhodopsin (VirR) and the biosynthesis pathway for the required chromophore, the second finding of a choanoflagellate-associated virus that encodes these genes. In co-collected metatranscriptomes, 85% of cPacV-1605 genes were expressed, with capsids, heat shock proteins and proteases among the most highly expressed. Based on orthologue presence–absence patterns across the PacVs and other eukaryotic viruses, we posit the observed viral groupings are connected to host lifestyles as heterotrophs or phototrophs.

This article is part of a discussion meeting issue ‘Single cell ecology’.

Metagenomic Analysis of the Diversity of DNA Viruses in the Surface and Deep Sea of the South China Sea

A metagenomic analysis of the viral community from five surface and five deep sea water (>2000 m below the surface, mbs) samples collected from the central basin of the South China Sea and adjacent Northwest Pacific Ocean during July-August 2017 was conducted herein. We builded up a South China Sea DNA virome (SCSV) dataset of 29,967 viral Operational Taxonomic Units (vOTUs), which is comparable to the viral populations from the original Tara Ocean and Malaspina expeditions. The most abundant and widespread viral populations were from the uncultivated viruses annotated from the viral metagenomics. Only 74 and 37 vOTUs have similarity with the reported genomes from the cultivated viruses and the single-virus genomics, respectively. The community structures of deep sea viromes in the SCSV were generally different from the surface viromes. The carbon flux and nutrients (PO₄ and NOx) were related to the surface and deep sea viromes in the SCSV, respectively. In the SCSV, the annotated vOTUs could be affiliated to the cultivated viruses mainly including Pelagibacter (SAR11) phage HTVC010P, Prochlorococcus phages (P-GSP1, P-SSM4, and P-TIM68), Cyanophages (MED4-184 and MED4-117) and Mycobacterium phages (Sparky and Squirty). It indicated that phage infection to the SAR11 cluster may occur ubiquitously and has significant impacts on bathypelagic SAR11 communities in the deep sea. Meanwhile, as Prochlorococcus is prominently distributed in the euphotic ocean, the existence of their potential phages in the deep sea suggested the sedimentation mechanism might contribute to the formation of the deep sea viromes. Intriguingly, the presence of Mycobacterium phages only in the deep sea viromes, suggests inhabitance of endemic viral populations in the deep sea viromes in the SCSV. This study provided an insight of the viral community in the South China Sea and for the first time uncovered the deep sea viral diversity in the central basin of the South China Sea.

The structure of a highly-conserved picocyanobacterial protein reveals a Tudor domain with an RNA-binding function

Cyanobacteria of the Prochlorococcus and marine Synechococcus genera are the most abundant photosynthetic microbes in the ocean. Intriguingly, the genomes of these bacteria are strongly divergent even within each genus, both in gene content and at the amino acid level of the encoded proteins. One striking exception to this is a 62-amino-acid protein, termed Prochlorococcus/ Synechococcus hyper-conserved protein (PSHCP). PSHCP is not only found in all sequenced Prochlorococcus and marine Synechococcus genomes, but it is also nearly 100% identical in its amino acid sequence across all sampled genomes. Such universal distribution and sequence conservation suggest an essential cellular role of PSHCP in these bacteria. However, its function is unknown. Here, we used NMR spectroscopy to determine its structure, finding that 53 of the 62 amino acids in PSHCP form a Tudor domain, whereas the remainder of the protein is disordered. NMR titration experiments revealed that PSHCP has only a weak affinity for DNA, but an 18.5-fold higher affinity for tRNA, hinting at an involvement of PSHCP in translation. Isothermal titration calorimetry experiments further revealed that PSHCP also binds single-stranded, double-stranded, and hairpin RNAs. These results provide the first insight into the structure and function of PSHCP, suggesting that PSHCP appears to be an RNA-binding protein that can recognize a broad array of RNA molecules.

Protein-level assembly increases protein sequence recovery from metagenomic samples manyfold

The open-source de novo protein-level assembler, Plass ( https://plass.mmseqs.com ), assembles six-frame-translated sequencing reads into protein sequences. It recovers 2-10 times more protein sequences from complex metagenomes and can assemble huge datasets. We assembled two redundancy-filtered reference protein catalogs, 2 billion sequences from 640 soil samples (soil reference protein catalog) and 292 million sequences from 775 marine eukaryotic metatranscriptomes (marine eukaryotic reference catalog), the largest free collections of protein sequences.

Methylphosphonate Oxidation in Prochlorococcus Strain MIT9301 Supports Phosphate Acquisition, Formate Excretion, and Carbon Assimilation into Purines

Until recently, MPn was only known to be degraded in the environment by the bacterial carbon-phosphorus (CP) lyase pathway, a reaction that releases the greenhouse gas methane. The identification of a formate-yielding MPn oxidative pathway in the marine planctomycete Gimesia maris (S. R. Gama, M. Vogt, T. Kalina, K. Hupp, et al., ACS Chem Biol 14:735–741, 2019, https://doi.org/10.1021/acschembio.9b00024) and the presence of this pathway in Prochlorococcus indicate that this compound can follow an alternative fate in the environment while providing a valuable source of P to organisms. In the ocean, where MPn is a major component of dissolved organic matter, the oxidation of MPn to formate by Prochlorococcus may direct the flow of this one-carbon compound to carbon dioxide or assimilation into biomass, thus limiting the production of methane.

The marine unicellular cyanobacterium Prochlorococcus is an abundant primary producer and widespread inhabitant of the photic layer in tropical and subtropical marine ecosystems, where the inorganic nutrients required for growth are limiting. In this study, we demonstrate that Prochlorococcus high-light strain MIT9301, an isolate from the phosphate-depleted subtropical North Atlantic Ocean, can oxidize methylphosphonate (MPn) and hydroxymethylphosphonate (HMPn), two phosphonate compounds present in marine dissolved organic matter, to obtain phosphorus. The oxidation of these phosphonates releases the methyl group as formate, which is both excreted and assimilated into purines in RNA and DNA. Genes encoding the predicted phosphonate oxidative pathway of MIT9301 were predominantly present in Prochlorococcus genomes from parts of the North Atlantic Ocean where phosphate availability is typically low, suggesting that phosphonate oxidation is an ecosystem-specific adaptation of some Prochlorococcus populations to cope with phosphate scarcity.

IMPORTANCE Until recently, MPn was only known to be degraded in the environment by the bacterial carbon-phosphorus (CP) lyase pathway, a reaction that releases the greenhouse gas methane. The identification of a formate-yielding MPn oxidative pathway in the marine planctomycete Gimesia maris (S. R. Gama, M. Vogt, T. Kalina, K. Hupp, et al., ACS Chem Biol 14:735–741, 2019, https://doi.org/10.1021/acschembio.9b00024) and the presence of this pathway in Prochlorococcus indicate that this compound can follow an alternative fate in the environment while providing a valuable source of P to organisms. In the ocean, where MPn is a major component of dissolved organic matter, the oxidation of MPn to formate by Prochlorococcus may direct the flow of this one-carbon compound to carbon dioxide or assimilation into biomass, thus limiting the production of methane.

Droplet Digital PCR for Estimating Absolute Abundances of Widespread Pelagibacter Viruses

Absolute abundances of prokaryotes are typically determined by FISH. Due to the lack of a universal conserved gene among all viruses, metagenomic fragment recruitment is commonly used to estimate the relative viral abundance. However, the paucity of absolute virus abundance data hinders our ability to fully understand how viruses drive global microbial populations. The cosmopolitan marine Pelagibacter ubique is host for the highly widespread HTVC010P pelagiphage isolate and the extremely abundant uncultured virus vSAG 37-F6 recently discovered by single-virus genomics. Here we applied droplet digital PCR (ddPCR) to calculate the absolute abundance of these pelagiphage genotypes in the Mediterranean Sea and the Gulf of Maine. Abundances were between 360 and 8,510 virus mL-1 and 1,270–14,400 virus mL-1 for vSAG 37-F6 and HTVC010P, respectively. Illumina PCR-amplicon sequencing corroborated the absence of ddPCR non-specific amplifications for vSAG 37-F6, but showed an overestimation of 6% for HTVC010P from off-targets, genetically unrelated viruses. Absolute abundances of both pelagiphages, two of the most abundance marine viruses, suggest a large viral pelagiphage diversity in marine environments, and show the efficiency and power of ddPCR to disentangle the structure of marine viral communities. Results also highlight the need for a standardized workflow to obtain accurate quantification that allows cross data comparison.

Emergence of trait variability through the lens of nitrogen assimilation in Prochlorococcus

Intraspecific trait variability has important consequences for the function and stability of marine ecosystems. Here we examine variation in the ability to use nitrate across hundreds of Prochlorococcus genomes to better understand the modes of evolution influencing intraspecific allocation of ecologically important functions. Nitrate assimilation genes are absent in basal lineages but occur at an intermediate frequency that is randomly distributed within recently emerged clades. The distribution of nitrate assimilation genes within clades appears largely governed by vertical inheritance, gene loss, and homologous recombination. By mapping this process onto a model of Prochlorococcus’ macroevolution, we propose that niche-constructing adaptive radiations and subsequent niche partitioning set the stage for loss of nitrate assimilation genes from basal lineages as they specialized to lower light levels. Retention of these genes in recently emerged lineages has likely been facilitated by selection as they sequentially partitioned into niches where nitrate assimilation conferred a fitness benefit.

Single-cell genomics uncover Pelagibacter as the putative host of the extremely abundant uncultured 37-F6 viral population in the ocean

The identification of relevant virus-host pairs that globally account for a large pool of carbon and nutrients in the ocean is paramount to build accurate ecological models. A previous work using single-virus genomics led to the discovery of the uncultured single-virus vSAG 37-F6, originally sorted from the Mediterranean Sea (Blanes Bay Microbial Observatory), that represents one of the most abundant dsDNA viral population in the marine surface virosphere. Here, from same sampling site, we report that a Pelagibacter single-cell contained a viral member of vSAG 37-F6 population, by means of PCR screening of sorted, genome-amplified single cells with vSAG 37-F6-specific primers and whole-genome sequencing. Furthermore, viruses from this population were also found in three other Pelagibacter single cells from the South Pacific and Atlantic oceans. These new uncultured pelagiphages were genetically different from the previously characterized pelagiphage isolates. Data showed that the uncultured vSAG 37-F6 population represents the Pelagibacter phages that inhabit the sunlit ocean better, and contains a vast unrecognized microdiversity.