Host preferences of coexisting Perkinsea parasitoids during coastal dinoflagellate blooms.

Parasites in aquatic systems are highly diverse and ubiquitous. In marine environments, parasite-host interactions contribute substantially to shaping microbial communities, but their nature and complexity remain poorly understood. In this study, we examined the relationship between Perkinsea parasitoids and bloom-forming dinoflagellate species. Our aim was to determine whether parasite-host species interactions are specific and whether the diversity and distribution of parasitoids are shaped by their dinoflagellate hosts. Several locations along the Catalan coast (NW Mediterranean Sea) were sampled during the blooms of five dinoflagellate species and the diversity of Perkinsea was determined by combining cultivation-based methods with metabarcoding of the V4 region of 18S rDNA. Most known species of Parviluciferaceae, and others not yet described, were detected, some of them coexisting in the same coastal location, and with a wide distribution. The specific parasite-host interactions determined for each of the studied blooms demonstrated the host preferences exhibited by parasitoids in nature. The dominance of a species within the parasitoid community is driven by the presence and abundances of its preferred host(s). The absence of parasitoid species, often associated with a low abundance of their preferred hosts, suggested that high infection rates are reached only under conditions that favour parasitoid propagation, especially dinoflagellate blooms.

Temperature regulation of marine heterotrophic prokaryotes increases latitudinally as a breach between bottom-up and top-down controls.

Planktonic heterotrophic prokaryotes make up the largest living biomass and process most organic matter in the ocean. Determining when and where the biomass and activity of heterotrophic prokaryotes are controlled by resource availability (bottom-up), predation and viral lysis (top-down) or temperature will help in future carbon cycling predictions. We conducted an extensive survey across subtropical and tropical waters of the Atlantic, Indian and Pacific Oceans during the Malaspina 2010 Global Circumnavigation Expedition and assessed indices for these three types of controls at 109 stations (mostly from the surface to 4,000 m depth). Temperature control was approached by the apparent activation energy in eV (ranging from 0.46 to 3.41), bottom-up control by the slope of the log-log relationship between biomass and production rate (ranging from -0.12 to 1.09) and top-down control by an index that considers the relative abundances of heterotrophic nanoflagellates and viruses (ranging from 0.82 to 4.83). We conclude that temperature becomes dominant (i.e. activation energy >1.5 eV) within a narrow window of intermediate values of bottom-up (0.3-0.6) and top-down 0.8-1.2) controls. A pervasive latitudinal pattern of decreasing temperature regulation towards the Equator, regardless of the oceanic basin, suggests that the impact of global warming on marine microbes and their biogeochemical function will be more intense at higher latitudes. Our analysis predicts that 1°C ocean warming will result in increased biomass of heterotrophic prokaryoplankton only in waters with <26°C of mean annual surface temperature.

Marine protist diversity in European coastal waters and sediments as revealed by high-throughput sequencing.

Although protists are critical components of marine ecosystems, they are still poorly characterized. Here we analysed the taxonomic diversity of planktonic and benthic protist communities collected in six distant European coastal sites. Environmental deoxyribonucleic acid (DNA) and ribonucleic acid (RNA) from three size fractions (pico-, nano- and micro/mesoplankton), as well as from dissolved DNA and surface sediments were used as templates for tag pyrosequencing of the V4 region of the 18S ribosomal DNA. Beta-diversity analyses split the protist community structure into three main clusters: picoplankton-nanoplankton-dissolved DNA, micro/mesoplankton and sediments. Within each cluster, protist communities from the same site and time clustered together, while communities from the same site but different seasons were unrelated. Both DNA and RNA-based surveys provided similar relative abundances for most class-level taxonomic groups. Yet, particular groups were overrepresented in one of the two templates, such as marine alveolates (MALV)-I and MALV-II that were much more abundant in DNA surveys. Overall, the groups displaying the highest relative contribution were Dinophyceae, Diatomea, Ciliophora and Acantharia. Also, well represented were Mamiellophyceae, Cryptomonadales, marine alveolates and marine stramenopiles in the picoplankton, and Monadofilosa and basal Fungi in sediments. Our extensive and systematic sequencing of geographically separated sites provides the most comprehensive molecular description of coastal marine protist diversity to date.

A fungi hotspot deep in the ocean: explaining the presence of Gjaerumia minor in equatorial Pacific bathypelagic waters.

A plant parasite associated with the white haze disease in apples, the Basidiomycota Gjaerumia minor, has been found in most samples of the global bathypelagic ocean. An analysis of environmental 18S rDNA sequences on 12 vertical profiles of the Malaspina 2010 expedition shows that the relative abundance of this cultured species increases with depth while its distribution is remarkably different between the deep waters of the Pacific and Atlantic oceans, being present in higher concentrations in the former. This is evident from sequence analysis and a microscopic survey with a species-specific newly designed TSA-FISH probe. Several hints point to the hypothesis that G. minor is transported to the deep ocean attached to particles, and the absence of G. minor in bathypelagic Atlantic waters could then be explained by the absence of this organism in surface waters of the equatorial Atlantic. The good correlation of G. minor biomass with Apparent Oxygen Utilization, recalcitrant carbon and free-living prokaryotic biomass in South Pacific waters, together with the identification of the observed cells as yeasts and not as resting spores (teliospores), point to the possibility that once arrived at deep layers this species keeps on growing and thriving.

Marine picoplankton metagenomes and MAGs from eleven vertical profiles obtained by the Malaspina Expedition.

The Ocean microbiome has a crucial role in Earth's biogeochemical cycles. During the last decade, global cruises such as Tara Oceans and the Malaspina Expedition have expanded our understanding of the diversity and genetic repertoire of marine microbes. Nevertheless, there are still knowledge gaps regarding their diversity patterns throughout depth gradients ranging from the surface to the deep ocean. Here we present a dataset of 76 microbial metagenomes (MProfile) of the picoplankton size fraction (0.2-3.0 µm) collected in 11 vertical profiles covering contrasting ocean regions sampled during the Malaspina Expedition circumnavigation (7 depths, from surface to 4,000 m deep). The MProfile dataset produced 1.66 Tbp of raw DNA sequences from which we derived: 17.4 million genes clustered at 95% sequence similarity (M-GeneDB-VP), 2,672 metagenome-assembled genomes (MAGs) of Archaea and Bacteria (Malaspina-VP-MAGs), and over 100,000 viral genomic sequences. This dataset will be a valuable resource for exploring the functional and taxonomic connectivity between the photic and bathypelagic tropical and sub-tropical ocean, while increasing our general knowledge of the Ocean microbiome.

Cryo-electron microscopy of extracellular vesicles associated with the marine toxic dinoflagellate Alexandrium minutum.

Extracellular Vesicles (EVs) are likely an important strategy of transport and communication in marine microbial community. Their isolation and characterization from axenic culture of microbial eukaryotes represents a technological challenge not fully solved. Here, for the first time, we isolated EVs from a near-axenic culture of the toxic dinoflagellate Alexandrium minutum. Pictures of the isolated vesicles were done with Cryo TEM (Cryogenic Transmission Electron Microscopy). Based on their morphotype the EVs were clustered in five major groups (rounded, rounded electron-dense, lumen electron-dense, double and irregular) and each EV was measured resulting in an average size of 0.36 µm of diameter. Taking in account that in prokaryotes it has been demonstrated that EVs play an important role in the mechanism of toxicity, this descriptive work aims to be the first step to study the possible role of EVs in the toxicity of dinoflagellates.

Automated flow cytometry as a tool to obtain a fine-grain picture of marine prokaryote community structure along an entire oceanographic cruise.

On a standard oceanographic cruise, flow cytometry data are usually collected sparsely through a bottle-based sampling and with stations separated by kilometers leading to a fragmented view of the ecosystem; to improve the resolution of the datasets produced by this technique here it is proposed the application of an automatic method of sampling and staining. The system used consists of a flow-cytometer (Accuri-C6) connected to an automated continuous sampler (OC-300) that collects samples of marine surface waters every 15 min. We tested this system for five days during a brief Mediterranean cruise with the aim of estimating the abundance, relative size and phenotypic diversity of prokaryotes. Seawater was taken by a faucet linked to an inlet pump (ca. 5 m depth). Once the sample was taken, the Oncyt-300 stained it and sent it to the flow cytometer. A total of 366 samples were collected, effectively achieving a fine-grained scale view of microbial community composition both through space and time. A significative positive relationship was found comparing data obtained with the automatic method and 10 samples collected from the faucet but processed with the standard protocol. Abundance values retrieved varied from 3.56·105 cell mL-1 in the coastal area till 6.87 105 cell mL-1 in open waters, exceptional values were reached in the harbor area where abundances peaked to 1.28 106 cell mL-1. The measured features (abundance and size) were associated with metadata (temperature, salinity, conductivity) also taken in continuous, of which conductivity was the one that better explained the variability of abundance. A full 24 h measurement cycle was performed resulting in slightly higher median bacterial abundances values during daylight hours compared to night. Alpha diversity, calculated using computational cytometry techniques, showed a higher value in the coastal area above 41° of latitude and had a strong inverse relationship with both salinity and conductivity. This is the first time to our knowledge that the OC-300 is directly applied to the marine environment during an oceanographic cruise; due to its high-resolution, this set-up shows great potential both to cover large sampling areas, and to monitor day-night cycles in situ.

Marked changes in diversity and relative activity of picoeukaryotes with depth in the world ocean.

Microbial eukaryotes are key components of the ocean plankton. Yet, our understanding of their community composition and activity in different water layers of the ocean is limited, particularly for picoeukaryotes (0.2-3 µm cell size). Here, we examined the picoeukaryotic communities inhabiting different vertical zones of the tropical and subtropical global ocean: surface, deep chlorophyll maximum, mesopelagic (including the deep scattering layer and oxygen minimum zones), and bathypelagic. Communities were analysed by high-tthroughput sequencing of the 18S rRNA gene (V4 region) as represented by DNA (community structure) and RNA (metabolism), followed by delineation of Operational Taxonomic Units (OTUs) at 99% similarity. We found a stratification of the picoeukaryotic communities along the water column, with assemblages corresponding to the sunlit and dark ocean. Specific taxonomic groups either increased (e.g., Chrysophyceae or Bicosoecida) or decreased (e.g., Dinoflagellata or MAST-3) in abundance with depth. We used the rRNA:rDNA ratio of each OTU as a proxy of metabolic activity. The highest relative activity was found in the mesopelagic layer for most taxonomic groups, and the lowest in the bathypelagic. Altogether, we characterize the change in community structure and metabolic activity of picoeukaryotes with depth in the global ocean, suggesting a hotspot of activity in the mesopelagic.

Seasonal impact of grazing, viral mortality, resource availability and light on the group-specific growth rates of coastal Mediterranean bacterioplankton.

Estimation of prokaryotic growth rates is critical to understand the ecological role and contribution of different microbes to marine biogeochemical cycles. However, there is a general lack of knowledge on what factors control the growth rates of different prokaryotic groups and how these vary between sites and along seasons at a given site. We carried out several manipulation experiments during the four astronomical seasons in the coastal NW Mediterranean in order to evaluate the impact of grazing, viral mortality, resource competition and light on the growth and loss rates of prokaryotes. Gross and net growth rates of different bacterioplankton groups targeted by group-specific CARD-FISH probes and infrared microscopy (for aerobic anoxygenic phototrophs, AAP), were calculated from changes in cell abundances. Maximal group-specific growth rates were achieved when both predation pressure and nutrient limitation were experimentally minimized, while only a minimal effect of viral pressure on growth rates was observed; nevertheless, the response to predation removal was more remarkable in winter, when the bacterial community was not subjected to nutrient limitation. Although all groups showed increases in their growth rates when resource competition as well as grazers and viral pressure were reduced, Alteromonadaceae consistently presented the highest rates in all seasons. The response to light availability was generally weaker than that to the other factors, but it was variable between seasons. In summer and spring, the growth rates of AAP were stimulated by light whereas the growth of the SAR11 clade (likely containing proteorhodopsin) was enhanced by light in all seasons. Overall, our results set thresholds on bacterioplankton group-specific growth and mortality rates and contribute to estimate the seasonally changing contribution of various bacterioplankton groups to the function of microbial communities. Our results also indicate that the least abundant groups display the highest growth rates, contributing to the recycling of organic matter to a much greater extent than what their abundances alone would predict.

Large-scale ocean connectivity and planktonic body size.

Global patterns of planktonic diversity are mainly determined by the dispersal of propagules with ocean currents. However, the role that abundance and body size play in determining spatial patterns of diversity remains unclear. Here we analyse spatial community structure - ß-diversity - for several planktonic and nektonic organisms from prokaryotes to small mesopelagic fishes collected during the Malaspina 2010 Expedition. ß-diversity was compared to surface ocean transit times derived from a global circulation model, revealing a significant negative relationship that is stronger than environmental differences. Estimated dispersal scales for different groups show a negative correlation with body size, where less abundant large-bodied communities have significantly shorter dispersal scales and larger species spatial turnover rates than more abundant small-bodied plankton. Our results confirm that the dispersal scale of planktonic and micro-nektonic organisms is determined by local abundance, which scales with body size, ultimately setting global spatial patterns of diversity.

Unveiling the role and life strategies of viruses from the surface to the dark ocean.

Viruses are a key component of marine ecosystems, but the assessment of their global role in regulating microbial communities and the flux of carbon is precluded by a paucity of data, particularly in the deep ocean. We assessed patterns in viral abundance and production and the role of viral lysis as a driver of prokaryote mortality, from surface to bathypelagic layers, across the tropical and subtropical oceans. Viral abundance showed significant differences between oceans in the epipelagic and mesopelagic, but not in the bathypelagic, and decreased with depth, with an average power-law scaling exponent of -1.03 km-1 from an average of 7.76 × 106 viruses ml-1 in the epipelagic to 0.62 × 106 viruses ml-1 in the bathypelagic layer with an average integrated (0 to 4000 m) viral stock of about 0.004 to 0.044 g C m-2, half of which is found below 775 m. Lysogenic viral production was higher than lytic viral production in surface waters, whereas the opposite was found in the bathypelagic, where prokaryotic mortality due to viruses was estimated to be 60 times higher than grazing. Free viruses had turnover times of 0.1 days in the bathypelagic, revealing that viruses in the bathypelagic are highly dynamic. On the basis of the rates of lysed prokaryotic cells, we estimated that viruses release 145 Gt C year-1 in the global tropical and subtropical oceans. The active viral processes reported here demonstrate the importance of viruses in the production of dissolved organic carbon in the dark ocean, a major pathway in carbon cycling.

Large variability of bathypelagic microbial eukaryotic communities across the world's oceans.

In this work, we study the diversity of bathypelagic microbial eukaryotes (0.8-20 µm) in the global ocean. Seawater samples from 3000 to 4000 m depth from 27 stations in the Atlantic, Pacific and Indian Oceans were analyzed by pyrosequencing the V4 region of the 18S ribosomal DNA. The relative abundance of the most abundant operational taxonomic units agreed with the results of a parallel metagenomic analysis, suggesting limited PCR biases in the tag approach. Although rarefaction curves for single stations were seldom saturated, the global analysis of all sequences together suggested an adequate recovery of bathypelagic diversity. Community composition presented a large variability among samples, which was poorly explained by linear geographic distance. In fact, the similarity between communities was better explained by water mass composition (26% of the variability) and the ratio in cell abundance between prokaryotes and microbial eukaryotes (21%). Deep diversity appeared dominated by four taxonomic groups (Collodaria, Chrysophytes, Basidiomycota and MALV-II) appearing in different proportions in each sample. Novel diversity amounted to 1% of the pyrotags and was lower than expected. Our study represents an essential step in the investigation of bathypelagic microbial eukaryotes, indicating dominating taxonomic groups and suggesting idiosyncratic assemblages in distinct oceanic regions.

Benthic protists: the under-charted majority.

Marine protist diversity inventories have largely focused on planktonic environments, while benthic protists have received relatively little attention. We therefore hypothesize that current diversity surveys have only skimmed the surface of protist diversity in marine sediments, which may harbor greater diversity than planktonic environments. We tested this by analyzing sequences of the hypervariable V4 18S rRNA from benthic and planktonic protist communities sampled in European coastal regions. Despite a similar number of OTUs in both realms, richness estimations indicated that we recovered at least 70% of the diversity in planktonic protist communities, but only 33% in benthic communities. There was also little overlap of OTUs between planktonic and benthic communities, as well as between separate benthic communities. We argue that these patterns reflect the heterogeneity and diversity of benthic habitats. A comparison of all OTUs against the Protist Ribosomal Reference database showed that a higher proportion of benthic than planktonic protist diversity is missing from public databases; similar results were obtained by comparing all OTUs against environmental references from NCBI's Short Read Archive. We suggest that the benthic realm may therefore be the world's largest reservoir of marine protist diversity, with most taxa at present undescribed.

Global abundance of planktonic heterotrophic protists in the deep ocean.

The dark ocean is one of the largest biomes on Earth, with critical roles in organic matter remineralization and global carbon sequestration. Despite its recognized importance, little is known about some key microbial players, such as the community of heterotrophic protists (HP), which are likely the main consumers of prokaryotic biomass. To investigate this microbial component at a global scale, we determined their abundance and biomass in deepwater column samples from the Malaspina 2010 circumnavigation using a combination of epifluorescence microscopy and flow cytometry. HP were ubiquitously found at all depths investigated down to 4000 m. HP abundances decreased with depth, from an average of 72±19 cells ml(-1) in mesopelagic waters down to 11±1 cells ml(-1) in bathypelagic waters, whereas their total biomass decreased from 280±46 to 50±14 pg C ml(-1). The parameters that better explained the variance of HP abundance were depth and prokaryote abundance, and to lesser extent oxygen concentration. The generally good correlation with prokaryotic abundance suggested active grazing of HP on prokaryotes. On a finer scale, the prokaryote:HP abundance ratio varied at a regional scale, and sites with the highest ratios exhibited a larger contribution of fungi molecular signal. Our study is a step forward towards determining the relationship between HP and their environment, unveiling their importance as players in the dark ocean's microbial food web.

Patterns of rare and abundant marine microbial eukaryotes.

BACKGROUND: Biological communities are normally composed of a few abundant and many rare species. This pattern is particularly prominent in microbial communities, in which most constituent taxa are usually extremely rare. Although abundant and rare subcommunities may present intrinsic characteristics that could be crucial for understanding community dynamics and ecosystem functioning, microbiologists normally do not differentiate between them. Here, we investigate abundant and rare subcommunities of marine microbial eukaryotes, a crucial group of organisms that remains among the least-explored biodiversity components of the biosphere. We surveyed surface waters of six separate coastal locations in Europe, independently considering the picoplankton, nanoplankton, and microplankton/mesoplankton organismal size fractions. RESULTS: Deep Illumina sequencing of the 18S rRNA indicated that the abundant regional community was mostly structured by organismal size fraction, whereas the rare regional community was mainly structured by geographic origin. However, some abundant and rare taxa presented similar biogeography, pointing to spatiotemporal structure in the rare microeukaryote biosphere. Abundant and rare subcommunities presented regular proportions across samples, indicating similar species-abundance distributions despite taxonomic compositional variation. Several taxa were abundant in one location and rare in other locations, suggesting large oscillations in abundance. The substantial amount of metabolically active lineages found in the rare biosphere suggests that this subcommunity constitutes a diversity reservoir that can respond rapidly to environmental change. CONCLUSIONS: We propose that marine planktonic microeukaryote assemblages incorporate dynamic and metabolically active abundant and rare subcommunities, with contrasting structuring patterns but fairly regular proportions, across space and time.

General patterns of diversity in major marine microeukaryote lineages.

Microeukaryotes have vital roles for the functioning of marine ecosystems, but still some general characteristics of their current diversity and phylogeny remain unclear. Here we investigated both aspects in major oceanic microeukaryote lineages using 18S rDNA (V4-V5 hypervariable regions) sequences from public databases that derive from various marine environmental surveys. A very carefully and manually curated dataset of 8291 Sanger sequences was generated and subsequently split into 65 taxonomic groups (roughly to Class level based on KeyDNATools) prior to downstream analyses. First, we calculated genetic distances and clustered sequences into Operational Taxonomic Units (OTUs) using different distance cut-off levels. We found that most taxonomic groups had a maximum pairwise genetic distance of 0.25. Second, we used phylogenetic trees to study general evolutionary patterns. These trees confirmed our taxonomic classification and served to run Lineage Through Time (LTT) plots. LTT results indicated different cladogenesis dynamics across groups, with some displaying an early diversification and others a more recent one. Overall, our study provides an improved description of the microeukaryote diversity in the oceans in terms of genetic differentiation within groups as well as in the general phylogenetic structure. These results will be important to interpret the large amount of sequence data that is currently generated by High Throughput Sequencing technologies.

Diversity patterns and activity of uncultured marine heterotrophic flagellates unveiled with pyrosequencing.

Flagellated heterotrophic microeukaryotes have key roles for the functioning of marine ecosystems as they channel large amounts of organic carbon to the upper trophic levels and control the population sizes of bacteria and archaea. Still, we know very little on the diversity patterns of most groups constituting this evolutionary heterogeneous assemblage. Here, we investigate 11 groups of uncultured flagellates known as MArine STramenopiles (MASTs). MASTs are ecologically very important and branch at the base of stramenopiles. We explored the diversity patterns of MASTs using pyrosequencing (18S rDNA) in coastal European waters. We found that MAST groups range from highly to lowly diversified. Pyrosequencing (hereafter '454') allowed us to approach to the limits of taxonomic diversity for all MAST groups, which varied in one order of magnitude (tens to hundreds) in terms of operational taxonomic units (98% similarity). We did not evidence large differences in activity, as indicated by ratios of DNA:RNA-reads. Most groups were strictly planktonic, although we found some groups that were active in sediments and even in anoxic waters. The proportion of reads per size fraction indicated that most groups were composed of very small cells (∼2-5 µm). In addition, phylogenetically different assemblages appeared to be present in different size fractions, depths and geographic zones. Thus, MAST diversity seems to be highly partitioned in spatial scales. Altogether, our results shed light on these ecologically very important but poorly known groups of uncultured marine flagellates.