Metabolic trade-offs constrain the cell size ratio in a nitrogen-fixing symbiosis

Biological dinitrogen (N2) fixation is a key metabolic process exclusively performed by prokaryotes, some of which are symbiotic with eukaryotes. Species of the marine haptophyte algae Braarudosphaera bigelowii harbor the N2-fixing endosymbiotic cyanobacteria UCYN-A, which might be evolving organelle-like characteristics. We found that the size ratio between UCYN-A and their hosts is strikingly conserved across sublineages/species, which is consistent with the size relationships of organelles in this symbiosis and other species. Metabolic modeling showed that this size relationship maximizes the coordinated growth rate based on trade-offs between resource acquisition and exchange. Our findings show that the size relationships of N2-fixing endosymbionts and organelles in unicellular eukaryotes are constrained by predictable metabolic underpinnings and that UCYN-A is, in many regards, functioning like a hypothetical N2-fixing organelle (or nitroplast).

Non-cyanobacterial diazotrophs: global diversity, distribution, ecophysiology, and activity in marine waters

Biological dinitrogen (N2) fixation supplies nitrogen to the oceans, supporting primary productivity, and is carried out by some bacteria and archaea referred to as diazotrophs. Cyanobacteria are conventionally considered to be the major contributors to marine N2 fixation, but non-cyanobacterial diazotrophs (NCDs) have been shown to be distributed throughout ocean ecosystems. However, the biogeochemical significance of marine NCDs has not been demonstrated. This review synthesizes multiple datasets, drawing from cultivation-independent molecular techniques and data from extensive oceanic expeditions, to provide a comprehensive view into the diversity, biogeography, ecophysiology, and activity of marine NCDs. A NCD nifH gene catalog was compiled containing sequences from both PCR-based and PCR-free methods, identifying taxa for future studies. NCD abundances from a novel database of NCD nifH-based abundances were colocalized with environmental data, unveiling distinct distributions and environmental drivers of individual taxa. Mechanisms that NCDs may use to fuel and regulate N2 fixation in response to oxygen and fixed nitrogen availability are discussed, based on a metabolic analysis of recently available Tara Oceans expedition data. The integration of multiple datasets provides a new perspective that enhances understanding of the biology, ecology, and biogeography of marine NCDs and provides tools and directions for future research.

Diazotrophs are overlooked contributors to carbon and nitrogen export to the deep ocean

Diazotrophs are widespread microorganisms that alleviate nitrogen limitation in 60% of our oceans, thereby regulating marine productivity. Yet, the group-specific contribution of diazotrophs to organic matter export has not been quantified, which so far has impeded an accurate assessment of their impact on the biological carbon pump. Here, we examine the fate of five groups of globally-distributed diazotrophs by using an original combination of mesopelagic particle sampling devices across the subtropical South Pacific Ocean. We demonstrate that cyanobacterial and non-cyanobacterial diazotrophs are exported down to 1000 m depth. Surprisingly, group-specific export turnover rates point to a more efficient export of small unicellular cyanobacterial diazotrophs (UCYN) relative to the larger and filamentous Trichodesmium. Phycoerythrin-containing UCYN-B and UCYN-C-like cells were recurrently found embedded in large (>50 µm) organic aggregates or organized into clusters of tens to hundreds of cells linked by an extracellular matrix, presumably facilitating their export. Beyond the South Pacific, our data are supported by analysis of the Tara Oceans metagenomes collected in other ocean basins, extending the scope of our results globally. We show that, when diazotrophs are found in the euphotic zone, they are also systematically present in mesopelagic waters, suggesting their transport to the deep ocean. We thus conclude that diazotrophs are a significant part of the carbon sequestered in the deep ocean and, therefore, they need to be accounted in regional and global estimates of export.

Primer design for the amplification of the ammonium transporter genes from the uncultured haptophyte algal species symbiotic with the marine nitrogen-fixing cyanobacterium UCYN-A1

The multiple symbiotic partnerships between closely related species of the haptophyte algae Braarudosphaera bigelowii and the nitrogen-fixing cyanobacteria Candidatus Atelocyanobacterium thalassa (UCYN-A) contribute importantly to the nitrogen and carbon cycles in vast areas of the ocean. The diversity of the eukaryotic 18S rDNA phylogenetic gene marker has helped to identify some of these symbiotic haptophyte species, yet we still lack a genetic marker to assess its diversity at a finer scale. One of such genes is the ammonium transporter (amt) gene, which encodes the protein that might be involved in the uptake of ammonium from UCYN-A in these symbiotic haptophytes. Here, we designed three specific PCR primer sets targeting the amt gene of the haptophyte species (A1-Host) symbiotic with the open ocean UCYN-A1 sublineage, and tested them in samples collected from open ocean and near-shore environments. Regardless of the primer pair used at Station ALOHA, which is where UCYN-A1 is the pre-dominant UCYN-A sublineage, the most abundant amt amplicon sequence variant (ASV) was taxonomically classified as A1-Host. In addition, two out of the three PCR primer sets revealed the existence of closely-related divergent haptophyte amt ASVs (>95% nucleotide identity). These divergent amt ASVs had higher relative abundances than the haptophyte typically associated with UCYN-A1 in the Bering Sea, or co-occurred with the previously identified A1-Host in the Coral Sea, suggesting the presence of new diversity of closely-related A1-Hosts in polar and temperate waters. Therefore, our study reveals an overlooked diversity of haptophytes species with distinct biogeographic distributions partnering with UCYN-A, and provides new primers that will help to gain new knowledge of the UCYN-A/haptophyte symbiosis.

Deep ocean metagenomes provide insight into the metabolic architecture of bathypelagic microbial communities

The deep sea, the largest ocean's compartment, drives planetary-scale biogeochemical cycling. Yet, the functional exploration of its microbial communities lags far behind other environments. Here we analyze 58 metagenomes from tropical and subtropical deep oceans to generate the Malaspina Gene Database. Free-living or particle-attached lifestyles drive functional differences in bathypelagic prokaryotic communities, regardless of their biogeography. Ammonia and CO oxidation pathways are enriched in the free-living microbial communities and dissimilatory nitrate reduction to ammonium and H2 oxidation pathways in the particle-attached, while the Calvin Benson-Bassham cycle is the most prevalent inorganic carbon fixation pathway in both size fractions. Reconstruction of the Malaspina Deep Metagenome-Assembled Genomes reveals unique non-cyanobacterial diazotrophic bacteria and chemolithoautotrophic prokaryotes. The widespread potential to grow both autotrophically and heterotrophically suggests that mixotrophy is an ecologically relevant trait in the deep ocean. These results expand our understanding of the functional microbial structure and metabolic capabilities of the largest Earth aquatic ecosystem.

Critical Role of Light in the Growth and Activity of the Marine N2-Fixing UCYN-A Symbiosis

The unicellular N₂-fixing cyanobacteria UCYN-A live in symbiosis with haptophytes in the Braarudosphaera bigelowii lineage. Maintaining N₂-fixing symbioses between two unicellular partners requires tight coordination of multiple biological processes including cell growth and division and, in the case of the UCYN-A symbiosis, N₂ fixation of the symbiont and photosynthesis of the host. In this system, it is thought that the host photosynthesis supports the high energetic cost of N₂ fixation, and both processes occur during the light period. However, information on this coordination is very limited and difficult to obtain because the UCYN-A symbiosis has yet to be available in culture. Natural populations containing the UCYN-A2 symbiosis were manipulated to explore the effects of alterations of regular light and dark periods and inhibition of host photosynthesis on N₂ fixation (single cell N₂ fixation rates), nifH gene transcription, and UCYN-A2 cell division (fluorescent in situ hybridization and nifH gene abundances). The results showed that the light period is critical for maintenance of regular patterns of gene expression, N₂ fixation and symbiont replication and cell division. This study suggests a crucial role for the host as a producer of fixed carbon, rather than light itself, in the regulation and implementation of these cellular processes in UCYN-A.

UCYN-A3, a newly characterized open ocean sublineage of the symbiotic N2 -fixing cyanobacterium Candidatus Atelocyanobacterium thalassa

The symbiotic unicellular cyanobacterium Candidatus Atelocyanobacterium thalassa (UCYN-A) is one of the most abundant and widespread nitrogen (N₂ )-fixing cyanobacteria in the ocean. Although it remains uncultivated, multiple sublineages have been detected based on partial nitrogenase (nifH) gene sequences, including the four most commonly detected sublineages UCYN-A1, UCYN-A2, UCYN-A3 and UCYN-A4. However, very little is known about UCYN-A3 beyond the nifH sequences from nifH gene diversity surveys. In this study, single cell sorting, DNA sequencing, qPCR and CARD-FISH assays revealed discrepancies involving the identification of sublineages, which led to new information on the diversity of the UCYN-A symbiosis. 16S rRNA and nifH gene sequencing on single sorted cells allowed us to identify the 16S rRNA gene of the uncharacterized UCYN-A3 sublineage. We designed new CARD-FISH probes that allowed us to distinguish and observe UCYN-A2 in a coastal location (SIO Pier; San Diego) and UCYN-A3 in an open ocean location (Station ALOHA; Hawaii). Moreover, we reconstructed about 13% of the UCYN-A3 genome from Tara Oceans metagenomic data. Finally, our findings unveil the UCYN-A3 symbiosis in open ocean waters suggesting that the different UCYN-A sublineages are distributed along different size fractions of the plankton defined by the cell-size ranges of their prymnesiophyte hosts.

A myovirus encoding both photosystem I and II proteins enhances cyclic electron flow in infected Prochlorococcus cells

Cyanobacteria are important contributors to primary production in the open oceans. Over the past decade, various photosynthesis-related genes have been found in viruses that infect cyanobacteria (cyanophages). Although photosystem II (PSII) genes are common in both cultured cyanophages and environmental samples ^1-4 , viral photosystem I (vPSI) genes have so far only been detected in environmental samples ^5,6 . Here, we have used a targeted strategy to isolate a cyanophage from the tropical Pacific Ocean that carries a PSI gene cassette with seven distinct PSI genes (psaJF, C, A, B, K, E, D) as well as two PSII genes (psbA, D). This cyanophage, P-TIM68, belongs to the T4-like myoviruses, has a prolate capsid, a long contractile tail and infects Prochlorococcus sp. strain MIT9515. Phage photosynthesis genes from both photosystems are expressed during infection, and the resultant proteins are incorporated into membranes of the infected host. Moreover, photosynthetic capacity in the cell is maintained throughout the infection cycle with enhancement of cyclic electron flow around PSI. Analysis of metagenomic data from the Tara Oceans expedition ⁷ shows that phages carrying PSI gene cassettes are abundant in the tropical Pacific Ocean, composing up to 28% of T4-like cyanomyophages. They are also present in the tropical Indian and Atlantic Oceans. P-TIM68 populations, specifically, compose on average 22% of the PSI-gene-cassette carrying phages. Our results suggest that cyanophages carrying PSI and PSII genes are likely to maintain and even manipulate photosynthesis during infection of their Prochlorococcus hosts in the tropical oceans.

Viral to metazoan marine plankton nucleotide sequences from the Tara Oceans expedition

A unique collection of oceanic samples was gathered by the Tara Oceans expeditions (2009–2013), targeting plankton organisms ranging from viruses to metazoans, and providing rich environmental context measurements. Thanks to recent advances in the field of genomics, extensive sequencing has been performed for a deep genomic analysis of this huge collection of samples. A strategy based on different approaches, such as metabarcoding, metagenomics, single-cell genomics and metatranscriptomics, has been chosen for analysis of size-fractionated plankton communities. Here, we provide detailed procedures applied for genomic data generation, from nucleic acids extraction to sequence production, and we describe registries of genomics datasets available at the European Nucleotide Archive (ENA, www.ebi.ac.uk/ena). The association of these metadata to the experimental procedures applied for their generation will help the scientific community to access these data and facilitate their analysis. This paper complements other efforts to provide a full description of experiments and open science resources generated from the Tara Oceans project, further extending their value for the study of the world’s planktonic ecosystems.

Exploring Microdiversity in Novel Kordia sp. (Bacteroidetes) with Proteorhodopsin from the Tropical Indian Ocean via Single Amplified Genomes

Marine Bacteroidetes constitute a very abundant bacterioplankton group in the oceans that plays a key role in recycling particulate organic matter and includes several photoheterotrophic members containing proteorhodopsin. Relatively few marine Bacteroidetes species have been described and, moreover, they correspond to cultured isolates, which in most cases do not represent the actual abundant or ecologically relevant microorganisms in the natural environment. In this study, we explored the microdiversity of 98 Single Amplified Genomes (SAGs) retrieved from the surface waters of the underexplored North Indian Ocean, whose most closely related isolate is Kordia algicida OT-1. Using Multi Locus Sequencing Analysis (MLSA) we found no microdiversity in the tested conserved phylogenetic markers (16S rRNA and 23S rRNA genes), the fast-evolving Internal Transcribed Spacer and the functional markers proteorhodopsin and the beta-subunit of RNA polymerase. Furthermore, we carried out a Fragment Recruitment Analysis (FRA) with marine metagenomes to learn about the distribution and dynamics of this microorganism in different locations, depths and size fractions. This analysis indicated that this taxon belongs to the rare biosphere, showing its highest abundance after upwelling-induced phytoplankton blooms and sinking to the deep ocean with large organic matter particles. This uncultured Kordia lineage likely represents a novel Kordia species (Kordia sp. CFSAG39SUR) that contains the proteorhodopsin gene and has a widespread spatial and vertical distribution. The combination of SAGs and MLSA makes a valuable approach to infer putative ecological roles of uncultured abundant microorganisms.

Global distribution and vertical patterns of a prymnesiophyte-cyanobacteria obligate symbiosis

A marine symbiosis has been recently discovered between prymnesiophyte species and the unicellular diazotrophic cyanobacterium UCYN-A. At least two different UCYN-A phylotypes exist, the clade UCYN-A1 in symbiosis with an uncultured small prymnesiophyte and the clade UCYN-A2 in symbiosis with the larger Braarudosphaera bigelowii. We targeted the prymnesiophyte-UCYN-A1 symbiosis by double CARD-FISH (catalyzed reporter deposition-fluorescence in situ hybridization) and analyzed its abundance in surface samples from the MALASPINA circumnavigation expedition. Our use of a specific probe for the prymnesiophyte partner allowed us to verify that this algal species virtually always carried the UCYN-A symbiont, indicating that the association was also obligate for the host. The prymnesiophyte-UCYN-A1 symbiosis was detected in all ocean basins, displaying a patchy distribution with abundances (up to 500 cells ml(-1)) that could vary orders of magnitude. Additional vertical profiles taken at the NE Atlantic showed that this symbiosis occupied the upper water column and disappeared towards the Deep Chlorophyll Maximum, where the biomass of the prymnesiophyte assemblage peaked. Moreover, sequences of both prymnesiophyte partners were searched within a large 18S rDNA metabarcoding data set from the Tara-Oceans expedition around the world. This sequence-based analysis supported the patchy distribution of the UCYN-A1 host observed by CARD-FISH and highlighted an unexpected homogeneous distribution (at low relative abundance) of B. bigelowii in the open ocean. Our results demonstrate that partners are always in symbiosis in nature and show contrasted ecological patterns of the two related lineages.

Closing the gaps on the viral photosystem-I psaDCAB gene organization

Marine photosynthesis is largely driven by cyanobacteria, namely Synechococcus and Prochlorococcus. Genes encoding for photosystem (PS) I and II reaction centre proteins are found in cyanophages and are believed to increase their fitness. Two viral PSI gene arrangements are known, psaJF→C→A→B→K→E→D and psaD→C→A→B. The shared genes between these gene cassettes and their encoded proteins are distinguished by %G + C and protein sequence respectively. The data on the psaD→C→A→B gene organization were reported from only two partial gene cassettes coming from Global Ocean Sampling stations in the Pacific and Indian oceans. Now we have extended our search to 370 marine stations from six metagenomic projects. Genes corresponding to both PSI gene arrangements were detected in the Pacific, Indian and Atlantic oceans, confined to a strip along the equator (30°N and 30°S). In addition, we found that the predicted structure of the viral PsaA protein from the psaD→C→A→B organization contains a lumenal loop conserved in PsaA proteins from Synechococcus, but is completely absent in viral PsaA proteins from the psaJF→C→A→B→K→E→D gene organization and most Prochlorococcus strains. This may indicate a co-evolutionary scenario where cyanophages containing either of these gene organizations infect cyanobacterial ecotypes biogeographically restricted to the 30°N and 30°S equatorial strip.

Ocean plankton. Determinants of community structure in the global plankton interactome

Species interaction networks are shaped by abiotic and biotic factors. Here, as part of the Tara Oceans project, we studied the photic zone interactome using environmental factors and organismal abundance profiles and found that environmental factors are incomplete predictors of community structure. We found associations across plankton functional types and phylogenetic groups to be nonrandomly distributed on the network and driven by both local and global patterns. We identified interactions among grazers, primary producers, viruses, and (mainly parasitic) symbionts and validated network-generated hypotheses using microscopy to confirm symbiotic relationships. We have thus provided a resource to support further research on ocean food webs and integrating biological components into ocean models.

Ocean plankton. Structure and function of the global ocean microbiome

Microbes are dominant drivers of biogeochemical processes, yet drawing a global picture of functional diversity, microbial community structure, and their ecological determinants remains a grand challenge. We analyzed 7.2 terabases of metagenomic data from 243 Tara Oceans samples from 68 locations in epipelagic and mesopelagic waters across the globe to generate an ocean microbial reference gene catalog with >40 million nonredundant, mostly novel sequences from viruses, prokaryotes, and picoeukaryotes. Using 139 prokaryote-enriched samples, containing >35,000 species, we show vertical stratification with epipelagic community composition mostly driven by temperature rather than other environmental factors or geography. We identify ocean microbial core functionality and reveal that >73% of its abundance is shared with the human gut microbiome despite the physicochemical differences between these two ecosystems.

Life-style and genome structure of marine Pseudoalteromonas siphovirus B8b isolated from the northwestern Mediterranean Sea

Marine viruses (phages) alter bacterial diversity and evolution with impacts on marine biogeochemical cycles, and yet few well-developed model systems limit opportunities for hypothesis testing. Here we isolate phage B8b from the Mediterranean Sea using Pseudoalteromonas sp. QC-44 as a host and characterize it using myriad techniques. Morphologically, phage B8b was classified as a member of the Siphoviridae family. One-step growth analyses showed that this siphovirus had a latent period of 70 min and released 172 new viral particles per cell. Host range analysis against 89 bacterial host strains revealed that phage B8b infected 3 Pseudoalteromonas strains (52 tested, >99.9% 16S rRNA gene nucleotide identity) and 1 non-Pseudoaltermonas strain belonging to Alteromonas sp. (37 strains from 6 genera tested), which helps bound the phylogenetic distance possible in a phage-mediated horizontal gene transfer event. The Pseudoalteromonas phage B8b genome size was 42.7 kb, with clear structural and replication modules where the former were delineated leveraging identification of 16 structural genes by virion structural proteomics, only 4 of which had any similarity to known structural proteins. In nature, this phage was common in coastal marine environments in both photic and aphotic layers (found in 26.5% of available viral metagenomes), but not abundant in any sample (average per sample abundance was 0.65% of the reads). Together these data improve our understanding of siphoviruses in nature, and provide foundational information for a new ‘rare virosphere’ phage–host model system.

Validation of a new catalysed reporter deposition-fluorescence in situ hybridization probe for the accurate quantification of marine Bacteroidetes populations

Catalysed reporter deposition-fluorescence in situ hybridization (CARD-FISH) is a powerful approach to quantify bacterial taxa. In this study, we compare the performance of the widely used Bacteroidetes CF319a probe with the new CF968 probe. In silico analyses and tests with isolates demonstrate that CF319a hybridizes with non-Bacteroidetes sequences from the Rhodobacteraceae and Alteromonadaceae families. We test the probes' accuracy in 37 globally distributed marine samples and over two consecutive years at the Blanes Bay Microbial Observatory (NW Mediterranean). We also compared the CARD-FISH data with the Bacteroidetes 16S rRNA gene sequences retrieved from 27 marine metagenomes from the TARA Oceans expedition. We find no significant differences in abundances between both approaches, although CF319a targeted some unspecific sequences and both probes displayed different abundances of specific Bacteroidetes phylotypes. Our results demonstrate that quantitative estimations by using both probes are significantly different in certain oceanographic regions (Mediterranean Sea, Red Sea and Arabian Sea) and that CF968 shows seasonality within marine Bacteroidetes, notably large differences between summer and winter that is overlooked by CF319a. We propose CF968 as an alternative to CF319a for targeting the whole Bacteroidetes phylum since it has better coverage, greater specificity and overall better quantifies marine Bacteroidetes.

Metagenomic 16S rDNA Illumina tags are a powerful alternative to amplicon sequencing to explore diversity and structure of microbial communities

Sequencing of 16S rDNA polymerase chain reaction (PCR) amplicons is the most common approach for investigating environmental prokaryotic diversity, despite the known biases introduced during PCR. Here we show that 16S rDNA fragments derived from Illumina-sequenced environmental metagenomes (mi tags) are a powerful alternative to 16S rDNA amplicons for investigating the taxonomic diversity and structure of prokaryotic communities. As part of the Tara Oceans global expedition, marine plankton was sampled in three locations, resulting in 29 subsamples for which metagenomes were produced by shotgun Illumina sequencing (ca. 700 Gb). For comparative analyses, a subset of samples was also selected for Roche-454 sequencing using both shotgun (m454 tags; 13 metagenomes, ca. 2.4 Gb) and 16S rDNA amplicon (454 tags; ca. 0.075 Gb) approaches. Our results indicate that by overcoming PCR biases related to amplification and primer mismatch, mi tags may provide more realistic estimates of community richness and evenness than amplicon 454 tags. In addition, mi tags can capture expected beta diversity patterns. Using mi tags is now economically feasible given the dramatic reduction in high-throughput sequencing costs, having the advantage of retrieving simultaneously both taxonomic (Bacteria, Archaea and Eukarya) and functional information from the same microbial community.