Francisco M. Cornejo-Castillo

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.

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.

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.

Global diversity and biogeography of deep-sea pelagic prokaryotes

The deep-sea is the largest biome of the biosphere, and contains more than half of the whole ocean’s microbes. Uncovering their general patterns of diversity and community structure at a global scale remains a great challenge, as only fragmentary information of deep-sea microbial diversity exists based on regional-scale studies. Here we report the first globally comprehensive survey of the prokaryotic communities inhabiting the bathypelagic ocean using high-throughput sequencing of the 16S rRNA gene. This work identifies the dominant prokaryotes in the pelagic deep ocean and reveals that 50% of the operational taxonomic units (OTUs) belong to previously unknown prokaryotic taxa, most of which are rare and appear in just a few samples. We show that whereas the local richness of communities is comparable to that observed in previous regional studies, the global pool of prokaryotic taxa detected is modest (~3600 OTUs), as a high proportion of OTUs are shared among samples. The water masses appear to act as clear drivers of the geographical distribution of both particle-attached and free-living prokaryotes. In addition, we show that the deep-oceanic basins in which the bathypelagic realm is divided contain different particle-attached (but not free-living) microbial communities. The combination of the aging of the water masses and a lack of complete dispersal are identified as the main drivers for this biogeographical pattern. All together, we identify the potential of the deep ocean as a reservoir of still unknown biological diversity with a higher degree of spatial complexity than hitherto considered.

Delineating ecologically significant taxonomic units from global patterns of marine picocyanobacteria

Significance Metagenomics has become an accessible approach to study complex microbial communities thanks to the advent of high-throughput sequencing technologies. However, molecular ecology studies often face interpretation issues, notably due to the lack of reliable reference databases for assigning reads to the correct taxa and use of fixed cutoffs to delineate taxonomic groups. Here, we considerably refined the phylogeography of marine picocyanobacteria, responsible for about 25% of global marine productivity, by recruiting reads targeting a high-resolution marker from Tara Oceans metagenomes. By clustering lineages based on their distribution patterns, we showed that there is significant diversity at a finer resolution than the currently defined “ecotypes,” a diversity that is tightly controlled by environmental cues. Prochlorococcus and Synechococcus are the two most abundant and widespread phytoplankton in the global ocean. To better understand the factors controlling their biogeography, a reference database of the high-resolution taxonomic marker petB, encoding cytochrome b6, was used to recruit reads out of 109 metagenomes from the Tara Oceans expedition. An unsuspected novel genetic diversity was unveiled within both genera, even for the most abundant and well-characterized clades, and 136 divergent petB sequences were successfully assembled from metagenomic reads, significantly enriching the reference database. We then defined Ecologically Significant Taxonomic Units (ESTUs)—that is, organisms belonging to the same clade and occupying a common oceanic niche. Three major ESTU assemblages were identified along the cruise transect for Prochlorococcus and eight for Synechococcus. Although Prochlorococcus HLIIIA and HLIVA ESTUs codominated in iron-depleted areas of the Pacific Ocean, CRD1 and the yet-to-be cultured EnvB were the prevalent Synechococcus clades in this area, with three different CRD1 and EnvB ESTUs occupying distinct ecological niches with regard to iron availability and temperature. Sharp community shifts were also observed over short geographic distances—for example, around the Marquesas Islands or between southern Indian and Atlantic Oceans—pointing to a tight correlation between ESTU assemblages and specific physico-chemical parameters. Together, this study demonstrates that there is a previously overlooked, ecologically meaningful, fine-scale diversity within some currently defined picocyanobacterial ecotypes, bringing novel insights into the ecology, diversity, and biology of the two most abundant phototrophs on Earth.

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.

Particle‐association lifestyle is a phylogenetically conserved trait in bathypelagic prokaryotes

The free‐living (FL) and particle‐attached (PA) marine microbial communities have repeatedly been proved to differ in their diversity and composition in the photic ocean and also recently in the bathypelagic ocean at a global scale. However, although high taxonomic ranks exhibit preferences for a PA or FL mode of life, it remains poorly understood whether two clear lifestyles do exist and how these are distributed across the prokaryotic phylogeny. We studied the FL (<0.8 μm) and PA (0.8–20 μm) prokaryotes at 30 stations distributed worldwide within the bathypelagic oceanic realm (2150–4000 m depth) using high‐throughput sequencing of the small subunit ribosomal RNA gene (16S rRNA). A high proportion of the bathypelagic prokaryotes were mostly found either attached to particles or freely in the surrounding water but rarely in both types of environments. In particular, this trait was deeply conserved through their phylogeny, suggesting that the deep‐ocean particles and the surrounding water constitute two highly distinct niches and that transitions from one to the other have been rare at an evolutionary timescale. As a consequence, PA and FL communities had clear alpha‐ and beta‐diversity differences that exceeded the global‐scale geographical variation. Our study organizes the bathypelagic prokaryotic diversity into a reasonable number of ecologically coherent taxa regarding their association with particles, a first step for understanding which are the microbes responsible for the processing of the dissolved and particulate pools of organic matter that have a very different biogeochemical role in the deep ocean.

Cyanobacterial symbionts diverged in the late Cretaceous towards lineage-specific nitrogen fixation factories in single-celled phytoplankton

The unicellular cyanobacterium UCYN-A, one of the major contributors to nitrogen fixation in the open ocean, lives in symbiosis with single-celled phytoplankton. UCYN-A includes several closely related lineages whose partner fidelity, genome-wide expression and time of evolutionary divergence remain to be resolved. Here we detect and distinguish UCYN-A1 and UCYN-A2 lineages in symbiosis with two distinct prymnesiophyte partners in the South Atlantic Ocean. Both symbiotic systems are lineage specific and differ in the number of UCYN-A cells involved. Our analyses infer a streamlined genome expression towards nitrogen fixation in both UCYN-A lineages. Comparative genomics reveal a strong purifying selection in UCYN-A1 and UCYN-A2 with a diversification process ∼91 Myr ago, in the late Cretaceous, after the low-nutrient regime period occurred during the Jurassic. These findings suggest that UCYN-A diversified in a co-evolutionary process, wherein their prymnesiophyte partners acted as a barrier driving an allopatric speciation of extant UCYN-A lineages.

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.

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 samples1–4, viral photosystem I (vPSI) genes have so far only been detected in environmental samples5,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 expedition7 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.Isolation of a cyanophage encoding photosystem I genes reveals that these are expressed during infection and inserted into host membranes, resulting in enhanced electron flow, and that phage carrying these genes are abundant in marine environments.