Coral-associated nitrogen fixation rates and diazotrophic diversity on a nutrient-replete equatorial reef

The role of diazotrophs in coral physiology and reef biogeochemistry remains poorly understood, in part because N2 fixation rates and diazotrophic community composition have only been jointly analyzed in the tissue of one tropical coral species. We performed field-based 15N2 tracer incubations during nutrient-replete conditions to measure diazotroph-derived nitrogen (DDN) assimilation into three species of scleractinian coral (Pocillopora acuta, Goniopora columna, Platygyra sinensis). Using multi-marker metabarcoding (16S rRNA, nifH, 18S rRNA), we analyzed DNA- and RNA-based communities in coral tissue and skeleton. Despite low N2 fixation rates, DDN assimilation supplied up to 6% of the holobiont's N demand. Active coral-associated diazotrophs were chiefly Cluster I (aerobes or facultative anaerobes), suggesting that oxygen may control coral-associated diazotrophy. Highest N2 fixation rates were observed in the endolithic community (0.20 µg N cm-2 per day). While the diazotrophic community was similar between the tissue and skeleton, RNA:DNA ratios indicate potential differences in relative diazotrophic activity between these compartments. In Pocillopora, DDN was found in endolithic, host, and symbiont compartments, while diazotrophic nifH sequences were only observed in the endolithic layer, suggesting a possible DDN exchange between the endolithic community and the overlying coral tissue. Our findings demonstrate that coral-associated diazotrophy is significant, even in nutrient-rich waters, and suggest that endolithic microbes are major contributors to coral nitrogen cycling on reefs.

Comparison of Deep-Sea Picoeukaryotic Composition Estimated from the V4 and V9 Regions of 18S rRNA Gene with a Focus on the Hadal Zone of the Mariana Trench

Diversity of microbial eukaryotes is estimated largely based on sequencing analysis of the hypervariable regions of 18S rRNA genes. But the use of different regions of 18S rRNA genes as molecular markers may generate bias in diversity estimation. Here, we compared the differences between the two most widely used markers, V4 and V9 regions of the 18S rRNA gene, in describing the diversity of epipelagic, bathypelagic, and hadal picoeukaryotes in the Challenger Deep of the Mariana Trench, which is a unique and little explored environment. Generally, the V9 region identified more OTUs in deeper waters than V4, while the V4 region provided greater Shannon diversity than V9. In the epipelagic zone, where Alveolata was the dominant group, picoeukaryotic community compositions identified by V4 and V9 markers are similar at different taxonomic levels. However, in the deep waters, the results of the two datasets show clear differences. These differences were mainly contributed by Retaria, Fungi, and Bicosoecida. The primer targeting the V9 region has an advantage in amplifying Bicosoecids in the bathypelagic and hadal zone of the Mariana Trench, and its high abundance in V9 dataset pointed out the possibility of Bicosoecids as a dominant group in this environment. Chrysophyceae, Fungi, MALV-I, and Retaria were identified as the dominant picoeukaryotes in the bathypelagic and hadal zone and potentially play important roles in deep-sea microbial food webs and biogeochemical cycling by their phagotrophic, saprotrophic, and parasitic life styles. Overall, the use of different markers of 18S rRNA gene allows a better assessment and understanding of the picoeukaryotic diversity in deep-sea environments.

In situ Prokaryotic and Eukaryotic Communities on Microplastic Particles in a Small Headwater Stream in Germany

The ubiquitous use of plastic products in our daily life is often accompanied by improper disposal. The first interactions of plastics with organisms in the environment occur by overgrowth or biofilm formation on the particle surface, which can facilitate the ingestion by animals. In order to elucidate the colonization of plastic particles by prokaryotic and eukaryotic microorganisms in situ, we investigated microbial communities in biofilms on four different polymer types and on mineral particles in a small headwater stream 500 m downstream of a wastewater treatment plant in Germany. Microplastic and mineral particles were exposed to the free-flowing water for 4 weeks in spring and in summer. The microbial composition of the developing biofilm was analyzed by 16S and 18S amplicon sequencing. Despite the expected seasonal differences in the microbial composition of pro- and eukaryotic communities, we repeatedly observed polymer type-specific differentiation in both seasons. The order of polymer type-specific prokaryotic and eukaryotic community distances calculated by Robust Aitchison principal component analysis (PCA) was the same in spring and summer samples. However, the magnitude of the distance differed considerably between polymer types. Prokaryotic communities on polyethylene particles exhibited the most considerable difference to other particles in summer, while eukaryotic communities on polypropylene particles showed the most considerable difference to other spring samples. The most contributing bacterial taxa to the polyethylene-specific differentiation belong to the Planctomycetales, Saccharimonadales, Bryobacterales, uncultured Acidiomicrobia, and Gemmatimonadales. The most remarkable differences in eukaryotic microorganism abundances could be observed in several distinct groups of Ciliophora (ciliates) and Chlorophytes (green algae). Prediction of community functions from taxonomic abundances revealed differences between spring and summer, and – to a lesser extent – also between polymer types and mineral surfaces. Our results show that different microplastic particles were colonized by different biofilm communities. These findings may be used for advanced experimental designs to investigate the role of microorganisms on the fate of microplastic particles in freshwater ecosystems.

Environmental DNA metabarcoding for monitoring metazoan biodiversity in Antarctic nearshore ecosystems

Antarctic benthic ecosystems support high biodiversity but their characterization is limited to a few well-studied areas, due to the extreme environment and remoteness making access and sampling difficult. Our aim was to compare water and sediment as sources of environmental DNA (eDNA) to better characterise Antarctic benthic communities and further develop practical approaches for DNA-based biodiversity assessment in remote environments. We used a cytochrome c oxidase subunit I (COI) metabarcoding approach to characterise metazoan communities in 26 nearshore sites across 12 locations in the Vestfold Hills (East Antarctica) based on DNA extracted from either sediment cores or filtered seawater. We detected a total of 99 metazoan species from 12 phyla across 26 sites, with similar numbers of species detected in sediment and water eDNA samples. However, significantly different communities were detected in the two sample types at sites where both were collected (i.e., where paired samples were available). For example, nematodes and echinoderms were more likely to be detected exclusively in sediment and water eDNA samples, respectively. eDNA from water and sediment core samples are complementary sample types, with epifauna more likely to be detected in water column samples and infauna in sediment. More reference DNA sequences are needed for infauna/meiofauna to increase the proportion of sequences and number of taxa that can be identified. Developing a better understanding of the temporal and spatial dynamics of eDNA at low temperatures would also aid interpretation of eDNA signals from polar environments. Our results provide a preliminary scan of benthic metazoan communities in the Vestfold Hills, with additional markers required to provide a comprehensive biodiversity survey. However, our study demonstrates the choice of sample type for eDNA studies of benthic ecosystems (sediment, water or both) needs to be carefully considered in light of the research or monitoring question of interest.

Freshwater protists: unveiling the unexplored in a large floodplain system

Protists play a fundamental role in all ecosystems, but we are still far from estimating the total diversity of many lineages, in particular in highly diverse environments, such as freshwater. Here, we survey the protist diversity of the Paraná River using metabarcoding, and we applied an approach that includes sequence similarity and phylogeny to evaluate the degree of genetic novelty of the protists' communities against the sequences described in the reference database PR² . We observed that ~28% of the amplicon sequence variants were classified as novel according to their similarity with sequences from the reference database; most of them were related to heterotrophic groups traditionally overlooked in freshwater systems. This lack of knowledge extended to those groups within the green algae (Archaeplastida) that are well documented such as Mamiellophyceae, and also to the less studied Pedinophyceae, for which we found sequences representing novel deep-branching clusters. Among the groups with potential novel protists, Bicosoecida (Stramenopiles) were the best represented, followed by Codosiga (Opisthokonta), and the Perkinsea (Alveolata). This illustrates the lack of knowledge on freshwater planktonic protists and also the need for isolation and/or cultivation of new organisms to better understand their role in ecosystem functioning.

Dimethyl sulfide mediates microbial predator-prey interactions between zooplankton and algae in the ocean

Phytoplankton are key components of the oceanic carbon and sulfur cycles¹. During bloom events, some species can emit large amounts of the organosulfur volatile dimethyl sulfide (DMS) into the ocean and consequently the atmosphere, where it can modulate aerosol formation and affect climate^2,3. In aquatic environments, DMS plays an important role as a chemical signal mediating diverse trophic interactions. Yet, its role in microbial predator-prey interactions remains elusive with contradicting evidence for its role in either algal chemical defence or in the chemo-attraction of grazers to prey cells^4,5. Here we investigated the signalling role of DMS during zooplankton-algae interactions by genetic and biochemical manipulation of the algal DMS-generating enzyme dimethylsulfoniopropionate lyase (DL) in the bloom-forming alga Emiliania huxleyi⁶. We inhibited DL activity in E. huxleyi cells in vivo using the selective DL-inhibitor 2-bromo-3-(dimethylsulfonio)-propionate⁷ and overexpressed the DL-encoding gene in the model diatom Thalassiosira pseudonana. We showed that algal DL activity did not serve as an anti-grazing chemical defence but paradoxically enhanced predation by the grazer Oxyrrhis marina and other microzooplankton and mesozooplankton, including ciliates and copepods. Consumption of algal prey with induced DL activity also promoted O. marina growth. Overall, our results demonstrate that DMS-mediated grazing may be ecologically important and prevalent during prey-predator dynamics in aquatic ecosystems. The role of algal DMS revealed here, acting as an eat-me signal for grazers, raises fundamental questions regarding the retention of its biosynthetic enzyme through the evolution of dominant bloom-forming phytoplankton in the ocean.

pr2-primers: An 18S rRNA primer database for protists

Metabarcoding of microbial eukaryotes (collectively known as protists) has developed tremendously in the last decade, almost solely relying on the 18S rRNA gene. As microbial eukaryotes are extremely diverse, many primers and primer pairs have been developed. To cover a relevant and representative fraction of the protist community in a given study system, an informed primer choice is necessary, as no primer pair can target all protists equally well. As such, a smart primer choice is very difficult even for experts and there are very few online resources available to list existing primers. We built a database listing 285 primers and 83 unique primer pairs that have been used for eukaryotic 18S rRNA gene metabarcoding. In silico performance of primer pairs was tested against two sequence databases: PR² version 4.12.0 for eukaryotes and a subset of silva version 132 for bacteria and archaea. We developed an R-based web application enabling browsing of the database, visualization of the taxonomic distribution of the amplified sequences with the number of mismatches, and testing any user-defined primer or primer set (https://app.pr2-primers.org). Taxonomic specificity of primer pairs, amplicon size and location of mismatches can also be determined. We identified universal primer sets that matched the largest number of sequences and analysed the specificity of some primer sets designed to target certain groups. This tool enables guided primer choices that will help a wide range of researchers to include protists as part of their investigations.

Resolving cryptic species complexes in marine protists: phylogenetic haplotype networks meet global DNA metabarcoding datasets

Marine protists have traditionally been assumed to be lowly diverse and cosmopolitan. Yet, several recent studies have shown that many protist species actually consist of cryptic complexes of species whose members are often restricted to particular biogeographic regions. Nonetheless, detection of cryptic species is usually hampered by sampling coverage and application of methods (e.g. phylogenetic trees) that are not well suited to identify relatively recent divergence and ongoing gene flow. In this paper, we show how these issues can be overcome by inferring phylogenetic haplotype networks from global metabarcoding datasets. We use the Chaetoceros curvisetus (Bacillariophyta) species complex as study case. Using two complementary metabarcoding datasets (Ocean Sampling Day and Tara Oceans), we equally resolve the cryptic complex in terms of number of inferred species. We detect new hypothetical species in both datasets. Gene flow between most of species is absent, but no barcoding gap exists. Some species have restricted distribution patterns whereas others are widely distributed. Closely related taxa occupy contrasting biogeographic regions, suggesting that geographic and ecological differentiation drive speciation. In conclusion, we show the potential of the analysis of metabarcoding data with evolutionary approaches for systematic and phylogeographic studies of marine protists.

Testing cophylogeny between coral reef invertebrates and their bacterial and archaeal symbionts

Marine invertebrates harbour a complex suite of bacterial and archaeal symbionts, a subset of which are probably linked to host health and homeostasis. Within a complex microbiome it can be difficult to tease apart beneficial or parasitic symbionts from nonessential commensal or transient microorganisms; however, one approach is to detect strong cophylogenetic patterns between microbial lineages and their respective hosts. We employed the Procrustean approach to cophylogeny (PACo) on 16S rRNA gene derived microbial community profiles paired with COI, 18S rRNA and ITS1 host phylogenies. Second, we undertook a network analysis to identify groups of microbes that were co-occurring within our host species. Across 12 coral, 10 octocoral and five sponge species, each host group and their core microbiota (50% prevalence within host species replicates) had a significant fit to the cophylogenetic model. Independent assessment of each microbial genus and family found that bacteria and archaea affiliated to Endozoicomonadaceae, Spirochaetaceae and Nitrosopumilaceae have the strongest cophylogenetic signals. Further, local Moran's I measure of spatial autocorrelation identified 14 ASVs, including Endozoicomonadaceae and Spirochaetaceae, whose distributions were significantly clustered by host phylogeny. Four co-occurring subnetworks were identified, each of which was dominant in a different host group. Endozoicomonadaceae and Spirochaetaceae ASVs were abundant among the subnetworks, particularly one subnetwork that was exclusively comprised of these two bacterial families and dominated the octocoral microbiota. Our results disentangle key microbial interactions that occur within complex microbiomes and reveal long-standing, essential microbial symbioses in coral reef invertebrates.

A Comparison of DNA Metabarcoding and Microscopy Methodologies for the Study of Aquatic Microbial Eukaryotes

The procedures and methodologies employed to study microbial eukaryotic plankton have been thoroughly discussed. Two main schools exist—one insisting on classic microscopy methodologies and the other supporting modern high-throughput sequencing (DNA metabarcoding). However, few studies have attempted to combine both these approaches; most studies implement one method while ignoring the other. This work aims to contribute to this discussion and examine the advantages and disadvantages of each methodology by comparing marine plankton community results from microscopy and DNA metabarcoding. The results obtained by the two methodologies do not vary significantly for Bacillariophyta, although they do for Dinoflagellata and Ciliophora. The lower the taxonomic level, the higher the inconsistency between the two methodologies for all the studied groups. Considering the different characteristics of microscopy-based identification and DNA metabarcoding, this work underlines that each method should be chosen depending on the aims of the study. DNA metabarcoding provides a better estimate of the taxonomic richness of an ecosystem while microscopy provides more accurate quantitative results regarding abundance and biomass. In any case, the combined use of the two methods, if properly standardized, can provide much more reliable and accurate results for the study of marine microbial eukaryotes.

Seasonal plankton succession is in accordance with phycotoxin occurrence in Disko Bay, West Greenland

Harmful algal blooms (HABs) are occurring more frequently in the world's oceans, probably as a consequence of climate change. HABs have not been considered a serious concern in the Arctic, even though the Arctic warms faster than any other region. While phycotoxins and toxin-producing phytoplankton have been found in Arctic waters on several occasions, there is a lack of information on seasonal succession of species and whether the occurrence of harmful species correlates with the presence of their respective phycotoxins. Hence, there is no baseline to assess future changes of HABs in this area. Here, we investigated two periods, from winter to spring and from the spring bloom until summer, in Disko Bay, West Greenland and followed the succession of toxins and their producers using metabarcoding, as well as analyses of particulate and dissolved toxins. We observed a typical seasonal succession with a spring bloom dominated by diatoms, followed by dinoflagellates in summer, with the two most important potentially toxic taxa found being Pseudo-nitzschia spp. and Alexandrium ostenfeldii. The Pseudo-nitzschia spp. peak correlated with a clear increase in particulate domoic acid, reaching 0.05 pg/L. Presence of Alexandrium ostenfeldii could be linked to an increase in spirolides, up to 56.4 pg/L in the particulate phase. Generally, the majority of detected dissolved toxins followed the succession pattern of the particulate toxins with a delay in time. Our results further show that Arctic waters are a suitable habitat for various toxin producers and that the strong seasonality of this environment is reflected by changing abundances of different toxins that pose a potential threat to the ecosystem and its beneficiaries.

Diversity of dinoflagellate assemblages in coastal temperate and offshore tropical waters of Australia

BACKGROUND

Dinoflagellates are a ubiquitous and ecologically important component of marine phytoplankton communities, with particularly notable species including those associated with harmful algal blooms (HABs) and those that bioluminesce. High-throughput sequencing offers a novel approach compared to traditional microscopy for determining species assemblages and distributions of dinoflagellates, which are poorly known especially in Australian waters.

RESULTS

We assessed the composition of dinoflagellate assemblages in two Australian locations: coastal temperate Port Phillip Bay and offshore tropical waters of Davies Reef (Great Barrier Reef). These locations differ in certain environmental parameters reflecting latitude as well as possible anthropogenic influences. Molecular taxonomic assessment revealed more species than traditional microscopy, and it showed statistically significant differences in dinoflagellate assemblages between locations. Bioluminescent species and known associates of HABs were present at both sites. Dinoflagellates in both areas were mainly represented by the order Gymnodiniales (66%-82% of total sequence reads). In the warm waters of Davies Reef, Gymnodiniales were equally represented by the two superclades, Gymnodiniales sensu stricto (33%) and Gyrodinium (34%). In contrast, in cooler waters of Port Phillip Bay, Gymnodiniales was mainly represented by Gyrodinium (82%). In both locations, bioluminescent dinoflagellates represented up to 0.24% of the total sequence reads, with Protoperidinium the most abundant genus. HAB-related species, mainly represented by Gyrodinium, were more abundant in Port Phillip Bay (up to 47%) than at Davies Reef (28%), potentially reflecting anthropogenic influence from highly populated and industrial areas surrounding the bay. The entire assemblage of dinoflagellates, as well as the subsets of HAB and bioluminescent species, were strongly correlated with water quality parameters (R2 = 0.56-0.92). Significant predictors differed between the subsets: HAB assemblages were explained by salinity, temperature, dissolved oxygen, and total dissolved solids; whereas, bioluminescent assemblages were explained only by salinity and dissolved oxygen, and had greater variability.

CONCLUSION

High-throughput sequencing and genotyping revealed greater diversity of dinoflagellate assemblages than previously known in both subtropical and temperate Australian waters. Significant correlations of assemblage structure with environmental variables suggest the potential for explaining the distribution and composition of both HAB species and bioluminescent species.

Annual phytoplankton dynamics in coastal waters from Fildes Bay, Western Antarctic Peninsula

Year-round reports of phytoplankton dynamics in the West Antarctic Peninsula are rare and mainly limited to microscopy and/or pigment-based studies. We analyzed the phytoplankton community from coastal waters of Fildes Bay in the West Antarctic Peninsula between January 2014 and 2015 using metabarcoding of the nuclear and plastidial 18/16S rRNA gene from both size-fractionated and flow cytometry sorted samples. Overall 14 classes of photosynthetic eukaryotes were present in our samples with the following dominating: Bacillariophyta (diatoms), Pelagophyceae and Dictyochophyceae for division Ochrophyta, Mamiellophyceae and Pyramimonadophyceae for division Chlorophyta, Haptophyta and Cryptophyta. Each metabarcoding approach yielded a different image of the phytoplankton community with for example Prymnesiophyceae more prevalent in plastidial metabarcodes and Mamiellophyceae in nuclear ones. Diatoms were dominant in the larger size fractions and during summer, while Prymnesiophyceae and Cryptophyceae were dominant in colder seasons. Pelagophyceae were particularly abundant towards the end of autumn (May). In addition of Micromonas polaris and Micromonas sp. clade B3, both previously reported in Arctic waters, we detected a new Micromonas 18S rRNA sequence signature, close to, but clearly distinct from M. polaris, which potentially represents a new clade specific of the Antarctic. These results highlight the need for complementary strategies as well as the importance of year-round monitoring for a comprehensive description of phytoplankton communities in Antarctic coastal waters.

Empirical evidence for concerted evolution in the 18S rDNA region of the planktonic diatom genus Chaetoceros

Concerted evolution is a process of homogenisation of repetitive sequences within a genome through unequal crossing over and gene conversion. This homogenisation is never fully achieved because mutations always create new variants. Classically, concerted evolution has been detected as "noise" in electropherograms and these variants have been characterised through cloning and sequencing of subsamples of amplified products. However, this approach limits the number of detectable variants and provides no information about the abundance of each variant. In this study, we investigated concerted evolution by using environmental time-series metabarcoding data, single strain high-throughput sequencing (HTS) and a collection of Sanger reference barcode sequences. We used six species of the marine planktonic diatom genus Chaetoceros as study system. Abundance plots obtained from environmental metabarcoding and single strain HTS showed the presence of a haplotype far more abundant than all the others (the "dominant" haplotype) and identical to the reference sequences of that species obtained with Sanger sequencing. This distribution fitted best with Zipf's law among the rank abundance/ dominance models tested. Furthermore, in each strain 99% of reads showed a similarity of 99% with the dominant haplotype, confirming the efficiency of the homogenisation mechanism of concerted evolution. We also demonstrated that minor haplotypes found in the environmental samples are not only technical artefacts, but mostly intragenomic variation generated by incomplete homogenisation. Finally, we showed that concerted evolution can be visualised inferring phylogenetic networks from environmental data. In conclusion, our study provides an important contribution to the understanding of concerted evolution and to the interpretation of DNA barcoding and metabarcoding data based on multigene family markers.

Heterotrophic Foraminifera Capable of Inorganic Nitrogen Assimilation

Nitrogen availability often limits biological productivity in marine systems, where inorganic nitrogen, such as ammonium is assimilated into the food web by bacteria and photoautotrophic eukaryotes. Recently, ammonium assimilation was observed in kleptoplast-containing protists of the phylum foraminifera, possibly via the glutamine synthetase/glutamate synthase (GS/GOGAT) assimilation pathway imported with the kleptoplasts. However, it is not known if the ubiquitous and diverse heterotrophic protists have an innate ability for ammonium assimilation. Using stable isotope incubations (¹⁵N-ammonium and ¹³C-bicarbonate) and combining transmission electron microscopy (TEM) with quantitative nanoscale secondary ion mass spectrometry (NanoSIMS) imaging, we investigated the uptake and assimilation of dissolved inorganic ammonium by two heterotrophic foraminifera; a non-kleptoplastic benthic species, Ammonia sp., and a planktonic species, Globigerina bulloides. These species are heterotrophic and not capable of photosynthesis. Accordingly, they did not assimilate ¹³C-bicarbonate. However, both species assimilated dissolved ¹⁵N-ammonium and incorporated it into organelles of direct importance for ontogenetic growth and development of the cell. These observations demonstrate that at least some heterotrophic protists have an innate cellular mechanism for inorganic ammonium assimilation, highlighting a newly discovered pathway for dissolved inorganic nitrogen (DIN) assimilation within the marine microbial loop.

Functional Seasonality of Free-Living and Particle-Associated Prokaryotic Communities in the Coastal Adriatic Sea

Marine snow is an important habitat for microbes, characterized by chemical and physical properties contrasting those of the ambient water. The higher nutrient concentrations in marine snow lead to compositional differences between the ambient water and the marine snow-associated prokaryotic community. Whether these compositional differences vary due to seasonal environmental changes, however, remains unclear. Thus, we investigated the seasonal patterns of the free-living and marine snow-associated microbial community composition and their functional potential in the northern Adriatic Sea. Our data revealed seasonal patterns in both, the free-living and marine snow-associated prokaryotes. The two assemblages were more similar to each other in spring and fall than in winter and summer. The taxonomic distinctness resulted in a contrasting functional potential. Motility and adaptations to low temperature in winter and partly anaerobic metabolism in summer characterized the marine snow-associated prokaryotes. Free-living prokaryotes were enriched in genes indicative for functions related to phosphorus limitation in winter and in genes tentatively supplementing heterotrophic growth with proteorhodopsins and CO-oxidation in summer. Taken together, the results suggest a strong influence of environmental parameters on both free-living and marine snow-associated prokaryotic communities in spring and fall leading to higher similarity between the communities, while the marine snow habitat in winter and summer leads to a specific prokaryotic community in marine snow in these two seasons.

High-throughput sequencing of 5S-IGS in oaks: Exploring intragenomic variation and algorithms to recognize target species in pure and mixed samples

Measuring biological diversity is a crucial but difficult undertaking, as exemplified in oaks where complex patterns of morphological, ecological, biogeographical and genetic differentiation collide with traditional taxonomy, which measures biodiversity in number of species (or higher taxa). In this pilot study, we generated high-throughput sequencing amplicon data of the intergenic spacer of the 5S nuclear ribosomal DNA cistron (5S-IGS) in oaks, using six mock samples that differ in geographical origin, species composition and pool complexity. The potential of the marker for automated genotaxonomy applications was assessed using a reference data set of 1,770 5S-IGS cloned sequences, covering the entire taxonomic breadth and distribution range of western Eurasian Quercus, and applying similarity (blast) and evolutionary approaches (maximum-likelihood trees and Evolutionary Placement Algorithm). Both methods performed equally well, allowing correct identification of species in sections Ilex and Cerris in the pure and mixed samples, and main lineages shared by species of sect. Quercus. Application of different cut-off thresholds revealed that medium- to high-abundance (>10 or 25) sequences suffice for a net species identification of samples containing one or a few individuals. Lower thresholds identify phylogenetic correspondence with all target species in highly mixed samples (analogous to environmental bulk samples) and include rare variants pointing towards reticulation, incomplete lineage sorting, pseudogenic 5S units and in situ (natural) contamination. Our pipeline is highly promising for future assessments of intraspecific and interpopulation diversity, and of the genetic resources of natural ecosystems, which are fundamental to empower fast and solid biodiversity conservation programmes worldwide.

Capturing open ocean biodiversity: Comparing environmental DNA metabarcoding to the continuous plankton recorder

Environmental DNA (eDNA) metabarcoding is emerging as a novel, objective tool for monitoring marine metazoan biodiversity. Zooplankton biodiversity in the vast open ocean is currently monitored through continuous plankton recorder (CPR) surveys, using ship-based bulk plankton sampling and morphological identification. We assessed whether eDNA metabarcoding (2 L filtered seawater) could capture similar Southern Ocean zooplankton biodiversity as conventional CPR bulk sampling (~1,500 L filtered seawater per CPR sample). We directly compared eDNA metabarcoding with (a) conventional morphological CPR sampling and (b) bulk DNA metabarcoding of CPR collected plankton (two transects for each comparison, 40 and 44 paired samples, respectively). A metazoan-targeted cytochrome c oxidase I (COI) marker was used to characterize species-level diversity. In the 2 L seawater eDNA samples, this marker amplified large amounts of non-metazoan picoplanktonic algae, but eDNA metabarcoding still detected up to 1.6 times more zooplankton species than morphologically analysed bulk CPR samples. COI metabarcoding of bulk DNA samples mostly avoided nonmetazoan amplifications and recovered more zooplankton species than eDNA metabarcoding. However, eDNA metabarcoding detected roughly two thirds of metazoan species and identified similar taxa contributing to community differentiation across the subtropical front separating transects. We observed a diurnal pattern in eDNA data for copepods which perform diel vertical migrations, indicating a surprisingly short temporal eDNA signal. Compared to COI, a eukaryote-targeted 18S ribosomal RNA marker detected a higher proportion, but lower diversity, of metazoans in eDNA. With refinement and standardization of methodology, eDNA metabarcoding could become an efficient tool for monitoring open ocean biodiversity.

Diverse coral reef invertebrates exhibit patterns of phylosymbiosis

Microbiome assemblages of plants and animals often show a degree of correlation with host phylogeny; an eco-evolutionary pattern known as phylosymbiosis. Using 16S rRNA gene sequencing to profile the microbiome, paired with COI, 18S rRNA and ITS1 host phylogenies, phylosymbiosis was investigated in four groups of coral reef invertebrates (scleractinian corals, octocorals, sponges and ascidians). We tested three commonly used metrics to evaluate the extent of phylosymbiosis: (a) intraspecific versus interspecific microbiome variation, (b) topological comparisons between host phylogeny and hierarchical clustering (dendrogram) of host-associated microbial communities, and (c) correlation of host phylogenetic distance with microbial community dissimilarity. In all instances, intraspecific variation in microbiome composition was significantly lower than interspecific variation. Similarly, topological congruency between host phylogeny and the associated microbial dendrogram was more significant than would be expected by chance across all groups, except when using unweighted UniFrac distance (compared with weighted UniFrac and Bray-Curtis dissimilarity). Interestingly, all but the ascidians showed a significant positive correlation between host phylogenetic distance and associated microbial dissimilarity. Our findings provide new perspectives on the diverse nature of marine phylosymbioses and the complex roles of the microbiome in the evolution of marine invertebrates.

Do plankton reflect the environmental quality status? The case of a post-industrial Mediterranean Bay

While the effects of industrial contamination in coastal areas may persist for years in benthos communities, plankton should not show permanent impairments because of their high spatial dynamics, fast turnover times and pronounced seasonality. To test this hypothesis, in 2019 we conducted five surveys in the Bay of Pozzuoli (Gulf of Naples, Mediterranean Sea), in front of a dismissed steel factory and in the adjacent inshore coastal waters. High seasonal variability was observed for bacteria, phytoplankton and mesozooplankton, whereas plankton spatial gradients were relatively smooth during each survey. Plankton biomass and diversity did not reveal any effects of past industrial activities not even at the innermost stations of the Bay, which however showed some signals of present anthropogenic pressure. Hydrodynamic and morphological features likely play a prominent role in maintaining a relatively good status of the plankton of the Bay, which hints at the relevance of coastal circulation and meteorological dynamics to revitalize areas impacted by human activities.

Unveiling protist diversity associated with the Pacific oyster Crassostrea gigas using blocking and excluding primers

BACKGROUND

Microbiome of macroorganisms might directly or indirectly influence host development and homeostasis. Many studies focused on the diversity and distribution of prokaryotes within these assemblages, but the eukaryotic microbial compartment remains underexplored so far.

RESULTS

To tackle this issue, we compared blocking and excluding primers to analyze microeukaryotic communities associated with Crassostrea gigas oysters. High-throughput sequencing of 18S rRNA genes variable loops revealed that excluding primers performed better by not amplifying oyster DNA, whereas the blocking primer did not totally prevent host contaminations. However, blocking and excluding primers showed similar pattern of alpha and beta diversities when protist communities were sequenced using metabarcoding. Alveolata, Stramenopiles and Archaeplastida were the main protist phyla associated with oysters. In particular, Codonellopsis, Cyclotella, Gymnodinium, Polarella, Trichodina, and Woloszynskia were the dominant genera. The potential pathogen Alexandrium was also found in high abundances within some samples.

CONCLUSIONS

Our study revealed the main protist taxa within oysters as well as the occurrence of potential oyster pathogens. These new primer sets are promising tools to better understand oyster homeostasis and disease development, such as the Pacific Oyster Mortality Syndrome (POMS) targeting juveniles.

Performance of the melting seawater-ice elution method on the metabarcoding characterization of benthic protist communities

Massive amplicon sequencing approaches to characterize the diversity of microbial eukaryotes in sediments are scarce and controls about the effects introduced by different methods to recover DNA are lacking. In this study, we compare the performance of the melting seawater-ice elution method on the characterization of benthic protist communities by 18S rRNA gene metabarcoding with results obtained by direct cell lysis and DNA purification from sediments. Even though the most abundant operational taxonomic units were recovered by both methods, eluted samples yielded higher richness than samples undergoing direct lysis. Both treatments allowed recovering the same taxonomic groups, although we observed significant differences in terms of relative abundance for some of them. Dinoflagellata and Ciliophora strongly dominated the community in eluted samples (> 80% reads). In directly lysed samples, they only represented 37%, while groups like Fungi and Ochrophytes were highly represented (> 20% reads respectively). Our results show that the elution process yields a higher protist richness estimation, most likely as a result of the higher sample volume used to recover organisms as compared to commonly used volumes for direct benthic DNA purification. Motile groups, like dinoflagellates and ciliates, are logically more enriched during the elution process.

Temporal Variability and Ecological Interactions of Parasitic Marine Syndiniales in Coastal Protist Communities

Syndiniales are a ubiquitous group of protist parasites that infect and kill a wide range of hosts, including harmful bloom-forming dinoflagellates. Despite the importance of parasitism as an agent of plankton mortality, parasite-host dynamics remain poorly understood, especially over time, hindering the inclusion of parasitism in food web and ecosystem models. For a full year in the Skidaway River Estuary (Georgia), we employed weekly 18S rRNA sampling and co-occurrence network analysis to characterize temporal parasite-host infection dynamics of Syndiniales. Over the year, Syndiniales exhibited strong temporal variability, with higher relative abundance from June to October (7 to 28%) than other months in the year (0.01% to 6%). Nonmetric dimensional scaling of Syndiniales composition revealed tight clustering in June to October that coincided with elevated temperatures (23 to 31°C), though in general, abiotic factors poorly explained composition (canonical correspondence analysis [CCA] and partial least-squares [PLS]) and were less important in the network than biotic relationships. Syndiniales amplicon sequence variants (ASVs) were well represented in the co-occurrence network (20% of edges) and had significant positive associations (Spearman r > 0.7), inferred to be putative parasite-host relationships, with known dinoflagellate hosts (e.g., Akashiwo and Gymnodinium) and other protist groups (e.g., ciliates, radiolarians, and diatoms). Positive associations rarely involved a single Syndiniales and dinoflagellate species, implying flexible parasite-host infection dynamics. These findings provide insight into the temporal dynamics of Syndiniales over a full year and reinforce the importance of single-celled parasites in driving plankton population dynamics. Further empirical work is needed to confirm network interactions and to incorporate parasitism within the context of ecosystem models.IMPORTANCE Protist parasites in the marine alveolate group, Syndiniales, have been observed within infected plankton host cells for decades, and recently, global-scale efforts (Tara Ocean exploration) have confirmed their importance within microbial communities. Yet, protist parasites remain enigmatic, particularly with respect to their temporal dynamics and parasite-host interactions. We employed weekly 18S amplicon surveys over a full year in a coastal estuary, revealing strong temporal shifts in Syndiniales parasites, with highest relative abundance during warmer summer to fall months. Though influenced by temperature, Syndiniales population dynamics were also driven by a high frequency of biological interactions with other protist groups, as determined through co-occurrence network analysis. Parasitic interactions implied by the network highlighted a range of confirmed (dinoflagellates) and putative (diatoms) interactions and suggests parasites may be less selective in their preferred hosts. Understanding parasite-host dynamics over space and time will improve our ability to include parasitism as a loss term in microbial food web models.

Unraveling the Diversity of Eukaryotic Microplankton in a Large and Deep Perialpine Lake Using a High Throughput Sequencing Approach

The structure of microbial communities, microalgae, heterotrophic protozoa and fungi contributes to characterize food webs and productivity and, from an anthropogenic point of view, the qualitative characteristics of water bodies. Traditionally, in freshwater environments many investigations have been directed to the study of pelagic microalgae (“phytoplankton”) and periphyton (i.e., photosynthetic and mixotrophic protists) through the use of light microscopy (LM). While the number of studies on bacterioplankton communities have shown a substantial increase after the advent of high-throughput sequencing (HTS) approaches, the study of the composition, structure, and spatio-temporal patterns of microbial eukaryotes in freshwater environments was much less widespread. Moreover, the understanding of the correspondence between the relative phytoplankton abundances estimated by HTS and LM is still incomplete. Taking into account these limitations, this study examined the biodiversity and seasonality of the community of eukaryotic microplankton in the epilimnetic layer of a large and deep perialpine lake (Lake Garda) using HTS. The analyses were carried out at monthly frequency during 2014 and 2015. The results highlighted the existence of a rich and well diversified community and the presence of numerous phytoplankton taxa that were never identified by LM in previous investigations. Furthermore, the relative abundances of phytoplankton estimated by HTS and LM showed a significant relationship at different taxonomic ranks. In the 2 years of investigation, the temporal development of the whole micro-eukaryotic community showed a clear non-random and comparable distribution pattern, with the main taxonomic groups coherently distributed in the individual seasons. In perspective, the results obtained in this study highlight the importance of HTS approaches in assessing biodiversity and the relative importance of the main protist groups along environmental gradients, including those caused by anthropogenic impacts (e.g., eutrophication and climate change).

Human access impacts biodiversity of microscopic animals in sandy beaches

Whereas most work to understand impacts of humans on biodiversity on coastal areas has focused on large, conspicuous organisms, we highlight effects of tourist access on the diversity of microscopic marine animals (meiofauna). We used a DNA metabarcoding approach with an iterative and phylogeny-based approach for the taxonomic assignment of meiofauna and relate diversity patterns to the numbers of tourists accessing sandy beaches on an otherwise un-impacted island National Park. Tourist frequentation, independently of differences in sediment granulometry, beach length, and other potential confounding factors, affected meiofaunal diversity in the shallow “swash” zone right at the mean water mark; the impacts declined with water depth (up to 2 m). The indicated negative effect on meiofauna may have a consequence on all the biota including the higher trophic levels. Thus, we claim that it is important to consider restricting access to beaches in touristic areas, in order to preserve biodiversity.

Martínez et al. use DNA metabarcoding and a phylogeny-based approach to demonstrate the effects of tourist access on meiofauna diversity of beaches in Asinara National Park. Their results show that tourist frequentation decreases meiofaunal diversity at the shallow “swash” zone, and can be used to inform tourist access and management of beaches.

Comparative analyses of the V4 and V9 regions of 18S rDNA for the extant eukaryotic community using the Illumina platform

Illumina sequencing is a representative tool for understanding the massive diversity of microbial eukaryotes in natural ecosystems. Here, we investigated the eukaryotic community in a pond (salinity of 2–4) on Dokdo (island) in the East Sea, Korea, using Illumina sequencing with primer sets for the V4 and V9 regions of 18S rDNA from 2016 to 2018 for the first time. Totally, 1,413 operational taxonomic units (OTUs) and 915 OTUs were detected using the V9 and V4 primer sets, respectively. Taxonomic analyses of these OTUs revealed that although the V4 primer set failed to describe the extant diversity for some major sub-division groups, the V9 primer set represented their diversity. Moreover, the rare taxa with <1% of total reads were exclusively detected using V9 primer set. Hence, the diversity of the eukaryotic community can vary depending on the choice of primers. The Illumina sequencing data of the V9 region of 18S rDNA may be advantageous for estimating the richness of the eukaryotic community including a rare biosphere, whereas the simultaneous application of two biomarkers may be suitable for understanding the molecular phylogenetic relationships. We strongly recommend both biomarkers be used to assess the diversity and phylogenetic relationship within the eukaryotic community in natural samples.

Species detection and delineation in the marine planktonic diatoms Chaetoceros and Bacteriastrum through metabarcoding: making biological sense of haplotype diversity

High-throughput sequencing (HTS) metabarcoding is commonly applied to assess phytoplankton diversity. Usually, haplotypes are grouped into operational taxonomic units (OTUs) through clustering, whereby the resulting number of OTUs depends on chosen similarity thresholds. We applied, instead, a phylogenetic approach to infer taxa among 18S rDNA V4-metabarcode haplotypes gathered from 48 time-series samples using the marine planktonic diatoms Chaetoceros and Bacteriastrum as test case. The 73 recovered taxa comprised both solitary haplotypes and polytomies, the latter composed each of a highly abundant, dominant haplotype and one to several minor, peripheral haplotypes. The solitary and dominant haplotypes usually matched reference sequences, enabling species assignation of taxa. We hypothesise that the super-abundance of reads in dominant haplotypes results from the homogenization effect of concerted evolution. Reads of populous peripheral haplotypes and dominant haplotypes show comparable distribution patterns over the sample dates, suggesting that they are part of the same population. Many taxa revealed marked seasonality, with closely related ones generally showing distinct periodicity, whereas others occur year-round. Phylogenies inferred from metabarcode haplotypes enable delineation of biologically meaningful taxa, whereas OTUs resulting from clustering algorithms often deviate markedly from such taxa.

Variability and Community Composition of Marine Unicellular Eukaryote Assemblages in a Eutrophic Mediterranean Urban Coastal Area with Marked Plankton Blooms and Red Tides

The Thessaloniki Bay is a eutrophic coastal area which has been characterized in recent years by frequent and intense phytoplankton blooms and red tides. The aim of the study was to investigate the underexplored diversity of marine unicellular eukaryotes in four different sampling sites in Thessaloniki Bay during a year of plankton blooms, red tides, and mucilage aggregates. High-Throughput Sequencing (HTS) was applied in extracted DNA from weekly water samples targeting the 18S rRNA gene. In almost all samples, phytoplankton blooms and/or red tides and mucilage aggregates were observed. The metabarcoding analysis has detected the known unicellular eukaryotic groups frequently observed in the Bay, dominated by Bacillariophyta and Dinoflagellata, and revealed taxonomic groups previously undetected in the study area (MALVs, MAST, and Cercozoa). The dominant OTUs were closely related to species known to participate in red tides, harmful blooms, and mucilage aggregates. Other OTUs, present also during the blooms in low abundance (number of reads), were closely related to known harmful species, suggesting the occurrence of rare taxa with potential negative impacts on human health not detectable with classical microscopy. Overall, the unicellular eukaryote assemblages showed temporal patterns rather than small-scale spatial separation responding to the variability of physical and chemical factors.

A metagenomic assessment of microbial eukaryotic diversity in the global ocean

Surveying microbial diversity and function is accomplished by combining complementary molecular tools. Among them, metagenomics is a PCR free approach that contains all genetic information from microbial assemblages and is today performed at a relatively large scale and reasonable cost, mostly based on very short reads. Here, we investigated the potential of metagenomics to provide taxonomic reports of marine microbial eukaryotes. We prepared a curated database with reference sequences of the V4 region of 18S rDNA clustered at 97% similarity and used this database to extract and classify metagenomic reads. More than half of them were unambiguously affiliated to a unique reference whilst the rest could be assigned to a given taxonomic group. The overall diversity reported by metagenomics was similar to that obtained by amplicon sequencing of the V4 and V9 regions of the 18S rRNA gene, although either one or both of these amplicon surveys performed poorly for groups like Excavata, Amoebozoa, Fungi and Haptophyta. We then studied the diversity of picoeukaryotes and nanoeukaryotes using 91 metagenomes from surface down to bathypelagic layers in different oceans, unveiling a clear taxonomic separation between size fractions and depth layers. Finally, we retrieved long rDNA sequences from assembled metagenomes that improved phylogenetic reconstructions of particular groups. Overall, this study shows metagenomics as an excellent resource for taxonomic exploration of marine microbial eukaryotes.

Cryptic species in the parasitic Amoebophrya species complex revealed by a polyphasic approach

As critical primary producers and recyclers of organic matter, the diversity of marine protists has been extensively explored by high-throughput barcode sequencing. However, classification of short metabarcoding sequences into traditional taxonomic units is not trivial, especially for lineages mainly known by their genetic fingerprints. This is the case for the widespread Amoebophrya ceratii species complex, parasites of their dinoflagellate congeners. We used genetic and phenotypic characters, applied to 119 Amoebophrya individuals sampled from the same geographic area, to construct practical guidelines for species delineation that could be applied in DNA/RNA based diversity analyses. Based on the internal transcribed spacer (ITS) regions, ITS2 compensatory base changes (CBC) and genome k-mer comparisons, we unambiguously defined eight cryptic species among closely related ribotypes that differed by less than 97% sequence identity in their SSU rDNA. We then followed the genetic signatures of these parasitic species during a three-year survey of Alexandrium minutum blooms. We showed that these cryptic Amoebophrya species co-occurred and shared the same ecological niche. We also observed a maximal ecological fitness for parasites having narrow to intermediate host ranges, reflecting a high cost for infecting a broader host range. This study suggests that a complete taxonomic revision of these parasitic dinoflagellates is long overdue to understand their diversity and ecological role in the marine plankton.

Habitat Heterogeneity and Connectivity: Effects on the Planktonic Protist Community Structure at Two Adjacent Coastal Sites (the Lagoon and the Gulf of Venice, Northern Adriatic Sea, Italy) Revealed by Metabarcoding

The Lagoon of Venice (LoV) and the Gulf of Venice (GoV), two adjacent coastal Long Term Ecological Research (LTER) sites in the northern Adriatic Sea, represent a transitional/marine coupled ecosystem under the influence of regional and local factors. In this study, these sites were sampled on four dates from April 2016 to February 2017 for environmental DNA and relevant abiotic variables, aiming to assess the relative importance of habitat heterogeneity and connectivity in structuring the protist community. High Throughput Sequencing of V4-18S rRNA gene from 56 samples collected at seven stations produced ca 6 million reads, grouped into 7,336 Operational Taxonomic Units (OTUs) at 97% similarity, which were affiliated to protists belonging to 34 taxonomic groups. The whole community was dominated by Bacillariophyta, especially in spring-summer in the LoV, and by Dinophyta, mainly in the GoV. Ciliophora, Syndiniales, and Cryptophyceae were the next more abundant groups. The community structure varied across the seasons and was different in the two ecosystems, which shared 96% of the reads but showed a high proportion of OTUs distributed preferentially in one of the two sites (specialists) and a different partitioning of trophic categories. GoV specialists were mainly Dinophyceae (>56%), followed by Syndiniales and Bacillariophyta, while the LoV specialists were distributed among several groups, including Bacillariophyta, Syndiniales, Ciliophora, Cryptophyceae, and Trebouxiophyceae. The main abiotic drivers of the differences between protist communities were salinity and temperature, which however explained a minor part of the variance (17%), pointing at a higher relevance of biotic factors and inter-taxa relationships. This was more evident in the LoV, where the network analysis highlighted a higher number of OTUs' connections than in the GoV. Overall, the metabarcoding approach allowed to depict the composition of the whole protist community in the lagoon and adjacent coastal waters with high resolution, revealing many taxa so far not reported in the area. In addition, despite no clear barrier to dispersal processes, differences in the relative abundance and temporal variability of local protist communities indicate that environmental heterogeneity, in these adjacent and connected ecosystems, can be strong enough to allow for ecological segregation.

Temporal and spatial dynamics of Bacteria, Archaea and protists in equatorial coastal waters

Singapore, an equatorial island in South East Asia, is influenced by a bi-annual reversal of wind directions which defines two monsoon seasons. We characterized the dynamics of the microbial communities of Singapore coastal waters by collecting monthly samples between February 2017 and July 2018 at four sites located across two straits with different trophic status, and sequencing the V6-V8 region of the small sub-unit ribosomal RNA gene (rRNA gene) of Bacteria, Archaea, and Eukaryota. Johor Strait, which is subjected to wider environmental fluctuations from anthropogenic activities, presented a higher abundance of copiotrophic microbes, including Cellvibrionales and Rhodobacterales. The mesotrophic Singapore Strait, where the seasonal variability is caused by changes in the oceanographic conditions, harboured a higher proportion of typically marine microbe groups such as Synechococcales, Nitrosupumilales, SAR11, SAR86, Marine Group II Archaea and Radiolaria. In addition, we observed seasonal variability of the microbial communities in the Singapore Strait, which was possibly influenced by the alternating monsoon regime, while no seasonal pattern was detected in the Johor Strait.

Metabarcoding successfully tracks temporal changes in eukaryotic communities in coastal sediments

Metabarcoding is a method that combines high-throughput DNA sequencing and DNA-based identification. Previously, this method has been successfully used to target spatial variation of eukaryote communities in marine sediments, however, the temporal changes in these communities remain understudied. Here, we follow the temporal changes of the eukaryote communities in Baltic Sea surface sediments collected from two coastal localities during three seasons of two consecutive years. Our study reveals that the structure of the sediment eukaryotic ecosystem was primarily driven by annual and seasonal changes in prevailing environmental conditions, whereas spatial variation was a less significant factor in explaining the variance in eukaryotic communities over time. Therefore, our data suggests that shifts in regional climate regime or large-scale changes in the environment are the overdriving factors in shaping the coastal eukaryotic sediment ecosystems rather than small-scale changes in local environmental conditions or heterogeneity in ecosystem structure. More studies targeting temporal changes are needed to further understand the long-term trends in ecosystem stability and response to climate change. Furthermore, this work contributes to the recent efforts in developing metabarcoding applications for environmental biomonitoring, proving a comprehensive option for traditional monitoring approaches.

Diversity and Composition of Pelagic Prokaryotic and Protist Communities in a Thin Arctic Sea-Ice Regime

One of the most prominent manifestations of climate change is the changing Arctic sea-ice regime with a reduction in the summer sea-ice extent and a shift from thicker, perennial multiyear ice towards thinner, first-year ice. These changes in the physical environment are likely to impact microbial communities, a key component of Arctic marine food webs and biogeochemical cycles. During the Norwegian young sea ICE expedition (N-ICE2015) north of Svalbard, seawater samples were collected at the surface (5 m), subsurface (20 or 50 m), and mesopelagic (250 m) depths on 9 March, 27 April, and 16 June 2015. In addition, several physical and biogeochemical data were recorded to contextualize the collected microbial communities. Through the massively parallel sequencing of the small subunit ribosomal RNA amplicon and metagenomic data, this work allows studying the Arctic's microbial community structure during the late winter to early summer transition. Results showed that, at compositional level, Alpha- (30.7%) and Gammaproteobacteria (28.6%) are the most frequent taxa across the prokaryotic N-ICE2015 collection, and also the most phylogenetically diverse. Winter to early summer trends were quite evident since there was a high relative abundance of thaumarchaeotes in the under-ice water column in late winter while this group was nearly absent during early summer. Moreover, the emergence of Flavobacteria and the SAR92 clade in early summer might be associated with the degradation of a spring bloom of Phaeocystis. High relative abundance of hydrocarbonoclastic bacteria, particularly Alcanivorax (54.3%) and Marinobacter (6.3%), was also found. Richness showed different patterns along the depth gradient for prokaryotic (highest at mesopelagic depth) and protistan communities (higher at subsurface depths). The microbial N-ICE2015 collection analyzed in the present study provides comprehensive new knowledge about the pelagic microbiota below drifting Arctic sea-ice. The higher microbial diversity found in late winter/early spring communities reinforces the need to continue with further studies to properly characterize the winter microbial communities under the pack-ice.

Development of an autonomous biosampler to capture in situ aquatic microbiomes

The importance of planktonic microbial communities is well acknowledged, since they are fundamental for several natural processes of aquatic ecosystems. Microorganisms naturally control the flux of nutrients, and also degrade and recycle anthropogenic organic and inorganic contaminants. Nevertheless, climate change effects and/or the runoff of nutrients/pollutants can affect the equilibrium of natural microbial communities influencing the occurrence of microbial pathogens and/or microbial toxin producers, which can compromise ecosystem environmental status. Therefore, improved microbial plankton monitoring is essential to better understand how these communities respond to environmental shifts. The study of marine microbial communities typically involves highly cost and time-consuming sampling procedures, which can limit the frequency of sampling and data availability. In this context, we developed and validated an in situ autonomous biosampler (IS-ABS) able to collect/concentrate in situ planktonic communities of different size fractions (targeting prokaryotes and unicellular eukaryotes) for posterior genomic, metagenomic, and/or transcriptomic analysis at a home laboratory. The IS-ABS field prototype is a small size and compact system able to operate up to 150 m depth. Water is pumped by a micropump (TCS MG2000) through a hydraulic circuit that allows in situ filtration of environmental water in one or more Sterivex filters placed in a filter cartridge. The IS-ABS also includes an application to program sampling definitions, allowing pre-setting configuration of the sampling. The efficiency of the IS-ABS was tested against traditional laboratory filtration standardized protocols. Results showed a good performance in terms of DNA recovery, as well as prokaryotic (16S rDNA) and eukaryotic (18S rDNA) community diversity analysis, using either methodologies. The IS-ABS automates the process of collecting environmental DNA, and is suitable for integration in water observation systems, what will contribute to substantially increase biological surveillances. Also, the use of highly sensitive genomic approaches allows a further study of the diversity and functions of whole or specific microbial communities.

Time series and beyond: multifaceted plankton research at a marine Mediterranean LTER site

Plankton are a pivotal component of the diversity and functioning of coastal marine ecosystems. A long time-series of observations is the best tool to trace their patterns and variability over multiple scales, ultimately providing a sound foundation for assessing, modelling and predicting the effects of anthropogenic and natural environmental changes on pelagic communities. At the same time, a long time-series constitutes a formidable asset for different kinds of research on specific questions that emerge from the observations, whereby the results of these complementary studies provide precious interpretative tools that augment the informative value of the data collected. In this paper, we review more than 140 studies that have been developed around a Mediterranean plankton time series gathered in the Gulf of Naples at the station LTER-MC since 1984. These studies have addressed different topics concerning marine plankton, which have included: i) seasonal patterns and trends; ii) taxonomic diversity, with a focus on key or harmful algal species and the discovery of many new taxa; iii) molecular diversity of selected species, groups of species or the whole planktonic community; iv) life cycles of several phyto- and zooplankton species; and v) interactions among species through trophic relationships, parasites and viruses. Overall, the products of this research demonstrate the great value of time series besides the record of fluctuations and trends, and highlight their primary role in the development of the scientific knowledge of plankton much beyond the local scale.

Protist Diversity and Seasonal Dynamics in Skagerrak Plankton Communities as Revealed by Metabarcoding and Microscopy

Protist community composition and seasonal dynamics are of major importance for the production of higher trophic levels, such as zooplankton and fish. Our aim was to reveal how the protist community in the Skagerrak changes through the seasons by combining high-throughput sequencing and microscopy of plankton collected monthly over two years. The V4 region of the 18S rRNA gene was amplified by eukaryote universal primers from the total RNA/cDNA. We found a strong seasonal variation in protist composition and proportional abundances, and a difference between two depths within the euphotic zone. Highest protist richness was found in late summer-early autumn, and lowest in winter. Temperature was the abiotic factor explaining most of the variation in diversity. Dinoflagellates was the most abundant and diverse group followed by ciliates and diatoms. We found about 70 new taxa recorded for the first time in the Skagerrak. The seasonal pattern in relative read abundance of major phytoplankton groups was well in accordance with microscopical biovolumes. This is the first metabarcoding study of the protist plankton community of all taxonomic groups and through seasons in the Skagerrak, which may serve as a baseline for future surveys to reveal effects of climate and environmental changes.

High contribution of Rhizaria (Radiolaria) to vertical export in the California Current Ecosystem revealed by DNA metabarcoding

Passive sinking of particulate organic matter (POM) is the main mechanism through which the biological pump transports surface primary production to the ocean interior. However, the contribution and variability of different biological sources to vertical export is not fully understood. Here, we use DNA metabarcoding of the 18S rRNA gene and particle interceptor traps (PITs) to characterize the taxonomic composition of particles sinking out of the photic layer in the California Current Ecosystem (CCE), a productive system with high export potential. The PITs included formalin-fixed and 'live' traps to investigate eukaryotic communities involved in the export and remineralization of sinking particles. Sequences affiliated with Radiolaria dominated the eukaryotic assemblage in fixed traps (90%), with Dinophyta and Metazoa making minor contributions. The prominence of Radiolaria decreased drastically in live traps, possibly due to selective consumption by copepods, heterotrophic nanoflagellates, and phaeodarians that were heavily enriched in these traps. These patterns were consistent across the water masses surveyed extending from the coast to offshore, despite major differences in productivity and trophic structure of the epipelagic plankton community. Our findings identify Radiolaria as major actors in export fluxes in the CCE.

Novel diversity within marine Mamiellophyceae (Chlorophyta) unveiled by metabarcoding

Mamiellophyceae (unicellular green algae) are a key phytoplankton group in coastal waters. Although extensively studied over the last 20 years, the overall oceanic distribution of the major species/clades is still poorly known. To address this problem, we analyzed the 2014 Ocean Sampling Day (OSD) metabarcoding dataset providing sequences from the V4 hypervariable region of the 18S rRNA gene for 157 samples collected at 143 mostly coastal stations. Mamiellophyceae were found at nearly all OSD stations and represented 55% of the green microalgae (Chlorophyta) reads. We performed phylogenetic analyses of unique OSD metabarcodes (amplicon single variants, ASVs) and GenBank reference sequences from cultures and from the environment, focusing on the four most represented genera: Ostreococcus (45% of the Mamiellophyceae reads), Micromonas (34%), Bathycoccus (10%) and Mantoniella (8.7%). These analyses uncovered novel diversity within each genus except Bathycoccus. In Ostreococcus, a new clade (E) was the second most represented clade after Ostreococcus "lucimarinus". Micromonas could be separated into nine clades, exceeding the six species and candidate species already described. Finally, we found two new environmental clades within Mantoniella. Each Mamiellophyceae clade had a specific distribution in the OSD dataset suggesting that they are adapted to different ecological niches.

Annotated 18S and 28S rDNA reference sequences of taxa in the planktonic diatom family Chaetocerotaceae

The species-rich diatom family Chaetocerotaceae is common in the coastal marine phytoplankton worldwide where it is responsible for a substantial part of the primary production. Despite its relevance for the global cycling of carbon and silica, many species are still described only morphologically, and numerous specimens do not fit any described taxa. Nowadays, studies to assess plankton biodiversity deploy high throughput sequencing metabarcoding of the 18S rDNA V4 region, but to translate the gathered metabarcodes into biologically meaningful taxa, there is a need for reference barcodes. However, 18S reference barcodes for this important family are still relatively scarce. We provide 18S rDNA and partial 28S rDNA reference sequences of 443 morphologically characterized chaetocerotacean strains. We gathered 164 of the 216 18S sequences and 244 of the 413 28S sequences of strains from the Gulf of Naples, Atlantic France, and Chile. Inferred phylogenies showed 84 terminal taxa in seven principal clades. Two of these clades included terminal taxa whose rDNA sequences contained spliceosomal and Group IC1 introns. Regarding the commonly used metabarcode markers in planktonic diversity studies, all terminal taxa can be discriminated with the 18S V4 hypervariable region; its primers fit their targets in all but two species, and the V4-tree topology is similar to that of the 18S. Hence V4-metabarcodes of unknown Chaetocerotaceae are assignable to the family. Regarding the V9 hypervariable region, most terminal taxa can be discriminated, but several contain introns in their primer targets. Moreover, poor phylogenetic resolution of the V9 region affects placement of metabarcodes of putative but unknown chaetocerotacean taxa, and hence, uncertainty in taxonomic assignment, even of higher taxa.

Diatom diversity through HTS-metabarcoding in coastal European seas

Diatoms constitute a diverse lineage of unicellular organisms abundant and ecologically important in aquatic ecosystems. Compared to other protists, their biology and taxonomy are well-studied, offering the opportunity to combine traditional approaches and new technologies. We examined a dataset of diatom 18S rRNA- and rDNA- (V4 region) reads from different plankton size-fractions and sediments from six European coastal marine sites, with the aim of identifying peculiarities and commonalities with respect to the whole protistan community. Almost all metabarcodes (99.6%) were assigned to known genera (121) and species (236), the most abundant of which were those already known from classic studies and coincided with those seen in light microscopy. rDNA and rRNA showed comparable patterns for the dominant taxa, but rRNA revealed a much higher diversity particularly in the sediment communities. Peculiar to diatoms is a tight bentho-pelagic coupling, with many benthic or planktonic species colonizing both water column and sediments and the dominance of planktonic species in both habitats. Overall metabarcoding results reflected the marked specificity of diatoms compared to other protistan groups in terms of morphological and ecological characteristics, at the same time confirming their great potential in the description of protist communities.

Temporal variability of bacterioplankton is habitat driven

Temporal dynamics of bacterioplankton are rarely investigated for multiple habitats and years within individual lakes, limiting our understanding of the variability of bacterioplankton community (BC) composition with respect to environmental factors. We assessed the BC composition of a littoral and two pelagic habitats (euphotic zone and hypolimnion) of Lake Tovel monthly from April 2014 to May 2017 by high-throughput sequencing of the V3-V4 hypervariable region of the 16S rRNA gene. The three habitats differed in temperature, light, oxygen and hydrology. In particular, the littoral was the most hydrologically unstable because it receives most of the lake inflow, the hypolimnion was the most stable because of its hydrologically sheltered position, and the pelagic euphotic habitat was intermediate. Consequently, we hypothesized different temporal patterns of BC composition for all three habitats according to their environmental differences. We applied PERMANOVA, nonmetric multidimensional scaling and source-sink analysis to characterize BC composition. Overall, BCs were different among habitats with the littoral showing the highest variability and the hypolimnion the highest stability. The BC of rainy 2014 was distinct from the BCs of other years irrespective of the habitats considered. Seasonal differences in BCs were limited to spring, probably linked to meltwater inflow and mixing. Thus, temporal effects related to year and season were linked to the hydrological gradient of habitats. We suggest that despite potential within-lake dispersal of bacterioplankton by water flow and mixing, local environmental conditions played a major role in Lake Tovel, fostering distinct BCs in the three habitats.

Green microalgae in marine coastal waters: The Ocean Sampling Day (OSD) dataset

The ecology and distribution of green phytoplankton (Chlorophyta) in the ocean is poorly known because most studies have focused on groups with large cell size such as diatoms or dinoflagellates that are easily recognized by traditional techniques such as microscopy. The Ocean Sampling Day (OSD) project sampled surface waters quasi-simultaneously at 141 marine locations, mostly in coastal waters. The analysis of the 18S V4 region OSD metabarcoding dataset reveals that Chlorophyta are ubiquitous and can be locally dominant in coastal waters. Chlorophyta represented 29% of the global photosynthetic reads (Dinoflagellates excluded) and their contribution was especially high at oligotrophic stations (up to 94%) and along the European Atlantic coast. Mamiellophyceae dominated most coastal stations. At some coastal stations, they were replaced by Chlorodendrophyceae, Ulvophyceae, Trebouxiophyceae or Chlorophyceae as the dominating group, while oligotrophic stations were dominated either by Chloropicophyceae or the uncultured prasinophytes clade IX. Several Chlorophyta classes showed preferences in terms of nitrate concentration, distance to the coast, temperature and salinity. For example, Chlorophyceae preferred cold and low salinity coastal waters, and prasinophytes clade IX warm, high salinity, oligotrophic oceanic waters.

Systematic, continental scale temporal monitoring of marine pelagic microbiota by the Australian Marine Microbial Biodiversity Initiative

Sustained observations of microbial dynamics are rare, especially in southern hemisphere waters. The Australian Marine Microbial Biodiversity Initiative (AMMBI) provides methodologically standardized, continental scale, temporal phylogenetic amplicon sequencing data describing Bacteria, Archaea and microbial Eukarya assemblages. Sequence data is linked to extensive physical, biological and chemical oceanographic contextual information. Samples are collected monthly to seasonally from multiple depths at seven sites: Darwin Harbour (Northern Territory), Yongala (Queensland), North Stradbroke Island (Queensland), Port Hacking (New South Wales), Maria Island (Tasmania), Kangaroo Island (South Australia), Rottnest Island (Western Australia). These sites span ~30° of latitude and ~38° longitude, range from tropical to cold temperate zones, and are influenced by both local and globally significant oceanographic and climatic features. All sequence datasets are provided in both raw and processed fashion. Currently 952 samples are publically available for bacteria and archaea which include 88,951,761 bacterial (72,435 unique) and 70,463,079 archaeal (24,205 unique) 16 S rRNA v1-3 gene sequences, and 388 samples are available for eukaryotes which include 39,801,050 (78,463 unique) 18 S rRNA v4 gene sequences.