Transition from a mixotrophic/heterotrophic protist community during the dark winter to a photoautotrophic spring community in surface waters of Disko Bay, Greenland

Unicellular eukaryotic plankton communities (protists) are the major basis of the marine food web. The spring bloom is especially important, because of its high biomass. However, it is poorly described how the protist community composition in Arctic surface waters develops from winter to spring. We show that mixotrophic and parasitic organisms are prominent in the dark winter period. The transition period toward the spring bloom event was characterized by a high relative abundance of mixotrophic dinoflagellates, while centric diatoms and the haptophyte Phaeocystis pouchetii dominated the successive phototrophic spring bloom event during the study. The data shows a continuous community shift from winter to spring, and not just a dormant spring community waiting for the right environmental conditions. The spring bloom initiation commenced while sea ice was still scattering and absorbing the sunlight, inhibiting its penetration into the water column. The initial increase in fluorescence was detected relatively deep in the water column at ~55 m depth at the halocline, at which the photosynthetic cells accumulated, while a thick layer of snow and sea ice was still obstructing sunlight penetration of the surface water. This suggests that water column stratification and a complex interplay of abiotic factors eventually promote the spring bloom initiation.

From genes to toxins: Profiling Prymnesium parvum during a riverine harmful algal bloom

Blooms of Prymnesium parvum, a unicellular alga globally distributed in marine and brackish environments, frequently result in massive fish kills due to the production of toxins called prymnesins by this haptophyte. In August 2022, a harmful algal bloom (HAB) of this species occurred in the lower Oder River (Poland and Germany), which caused mass mortalities of fish and other organisms. This HAB was linked to low discharge of the Oder and mining activities that caused a significant increase in salinity. In this context, we report on the molecular detection and screening of this haptophyte and its toxins in environmental samples and clonal cultures derived thereof. Both conventional PCR and droplet digital PCR assays reliably detected P. parvum in environmental samples. eDNA metabarcoding using the V4 region of the 18S rRNA gene revealed a single Prymnesium sequence variant, but failed to identify it to species level. Four clonal cultures established from environmental samples were unambiguously identified as P. parvum by molecular phylogenetics (near full-length 18S rRNA gene) and light microscopy. Phylogenetic analysis (ITS1-5.8S-ITS2 marker region) placed the cultured phylotype within a clade containing other P. parvum strains known to produce B-type prymnesins. Toxin-screening of the cultures using liquid chromatography-electrospray ionization - time of flight mass spectrometry identified B-type prymnesins, which were also detected in extracts of filter residues from water samples of the Oder collected during the HAB. Overall, our investigation provides a detailed characterization of P. parvum, including their prymnesins, during this HAB in the Oder River, contributing valuable insights into this ecological disaster. In addition, the droplet digital PCR assay established here will be useful for future monitoring of low levels of P. parvum on the Oder River or any other salt-impacted and brackish water bodies.

Marine Fungal Diversity and Dynamics in the Gulf of Trieste (Northern Adriatic Sea)

Fungi contribute to different important ecological processes, including decomposition of organic matter and nutrient cycling, but in the marine environment the main factors influencing their diversity and dynamics at the spatial and temporal levels are still largely unclear. In this study, we performed DNA metabarcoding on seawater sampled monthly over a year and a half in the Gulf of Trieste (northern Adriatic Sea), targeting the internal transcribed spacer (ITS) and the 18S rRNA gene regions. The fungal communities were diverse, very dynamic, and belonged predominantly to marine taxa. Samples could be clustered in two groups, mainly based on the high (> 30%) or low relative proportion of the ascomycetes Parengyodontium album, which emerged as a key taxon in this area. Dissolved and particulate organic C:N ratio played important roles in shaping the mycoplankton assemblages, suggesting that differently bioavailable organic matter pools may be utilized by different consortia. The proportion of fungal over total reads was 31% for ITS and 0.7% for 18S. ITS had the highest taxonomic resolution but low power to detect early divergent fungal lineages. Our results on composition, distribution, and environmental drivers extended our knowledge of the structure and function of the mycobiome of coastal waters.

A full-length 18S ribosomal DNA metabarcoding approach for determining protist community diversity using Nanopore sequencing

Protist diversity studies are frequently conducted using DNA metabarcoding methods. Currently, most studies have utilized short read sequences to assess protist diversity. One limitation of using short read sequences is the low resolution of the markers. For better taxonomic resolution longer sequences of the 18S rDNA are required because the full-length has both conserved and hypervariable regions. In this study, a new primer pair combination was used to amplify the full-length 18S rDNA and its efficacy was validated with a test community and then validated with field samples. Full-length sequences obtained with the Nanopore MinION for protist diversity from field samples were compared with Illumina MiSeq V4 and V8-V9 short reads. Sequences generated from the high-throughput sequencers are Amplicon Sequence Variants (ASVs). Metabarcoding results show high congruency among the long reads and short reads in taxonomic annotation at the major taxonomic group level; however, not all taxa could be successfully detected from sequences. Based on the criteria of ≥95% similarity and ≥1000 bp query length, 298 genera were identified by all markers in the field samples, 250 (84%) were detected by 18S, while only 226 (76%) by V4 and 213 (71%) by V8-V9. Of the total 85 dinoflagellate genera observed, 19 genera were not defined by 18S dinoflagellate ASVs compared to only three among the total 52 diatom genera. The discrepancy in this resolution is due to the lack of taxonomically available 18S reference sequences in particular for dinoflagellates. Overall, this preliminary investigation demonstrates that application of the full-length 18S rDNA approach can be successful in field studies.

Spring protistan communities in response to warming in the northeastern East China Sea

The northeastern East China Sea is a highly dynamic marine ecosystem influenced by seasonally varying water mass properties. However, despite being among the world's fastest-warming ocean, there has been limited investigation into the impacts of warming on protistan communities. We collected seawater from two stations (E42 and E46) with different natural protist communities and environmental attributes to investigate the acclimation of the two communities to artificially elevated temperatures (ambient T, +2, and +4 °C). Nutrient and Chl-a conditions reflected oceanographic differences, providing insights into protistan community dynamics. Notably, small-sized autotrophic protists prevailed in the phosphate-deficient E42 community, with mid-incubation heterotrophic conversions. Higher temperatures exacerbated the effects of the P deficiency on the E42 community. While the proportions of Bacillariophyta increased only in the nutrient-balanced E46 communities, those of mixotrophic dinoflagellates increased with elevated temperature, regardless of P deficiency, suggesting that mixotrophy likely aids adaptation in changing marine environments. In summary, the findings of this microcosm study illuminate the potential modulation of spring protistan communities in the northeastern East China Sea under anticipated future warming.

Temporal and spatial variability of constitutive mixotroph abundance and proportion

Mixotrophic plankton can comprise a substantial portion of the plankton community compared to phytoplankton and zooplankton. However, there is a gap in the understanding of conditions that influence mixotroph prevalence and activity in situ because current methods often over- or underestimate mixotroph abundance. A labeled prey-tracer method was utilized to identify active mixotrophs present at two locations in a temperate estuary over a year. The tracer method was combined with light microscopy data to estimate active mixotroph abundance and proportion. This study estimated that actively grazing mixotrophic taxa were more abundant in the spring and autumn compared to summer. Dinoflagellates typically dominated the mixotrophic taxa except during autumn at the low salinity location when cryptophytes dominated. Further analysis suggested that active mixotroph abundances might not be only regulated by environmental conditions favorable to mixotrophy but, instead, environmental conditions favorable to different mixotrophs utilization of phagotrophy. By focusing on mixotrophic taxa that were identified to be actively grazing at time of sampling, this study provided a more nuanced estimation of mixotroph abundance, increasing the understanding of how mixotrophic abundance and proportion in situ are influenced by the planktonic community composition and environmental factors.

A marine heatwave drives significant shifts in pelagic microbiology

Marine heatwaves (MHWs) cause disruption to marine ecosystems, deleteriously impacting macroflora and fauna. However, effects on microorganisms are relatively unknown despite ocean temperature being a major determinant of assemblage structure. Using data from thousands of Southern Hemisphere samples, we reveal that during an "unprecedented" 2015/16 Tasman Sea MHW, temperatures approached or surpassed the upper thermal boundary of many endemic taxa. Temperate microbial assemblages underwent a profound transition to niche states aligned with sites over 1000 km equatorward, adapting to higher temperatures and lower nutrient conditions bought on by the MHW. MHW conditions also modulate seasonal patterns of microbial diversity and support novel assemblage compositions. The most significant affects of MHWs on microbial assemblages occurred during warmer months, when temperatures exceeded the upper climatological bounds. Trends in microbial response across several MHWs in different locations suggest these are emergent properties of temperate ocean warming, which may facilitate monitoring, prediction and adaptation efforts.

Constraining the oxygen requirements for modern microbial eukaryote diversity

Eukaryotes originated prior to the establishment of modern marine oxygen (O2) levels. According to the body fossil and lipid biomarker records, modern (crown) microbial eukaryote lineages began diversifying in the ocean no later than ~800 Ma. While it has long been predicted that increasing atmospheric O2 levels facilitated the early diversification of microbial eukaryotes, the O2 levels needed to permit this diversification remain unconstrained. Using time-resolved geochemical parameter and gene sequence information from a model marine oxygen minimum zone spanning a range of dissolved O2 levels and redox states, we show that microbial eukaryote taxonomic richness and phylogenetic diversity remain the same until O2 declines to around 2 to 3% of present atmospheric levels, below which these diversity metrics become significantly reduced. Our observations suggest that increasing O2 would have only directly promoted early crown-eukaryote diversity if atmospheric O2 was below 2 to 3% of modern levels when crown-eukaryotes originated and then later met or surpassed this range as crown-eukaryotes diversified. If atmospheric O2 was already consistently at or above 2 to 3% of modern levels by the time that crown-eukaryotes originated, then the subsequent diversification of modern microbial eukaryotes was not directly driven by atmospheric oxygenation.

Evaluation of sequential filtration and centrifugation to capture environmental DNA and survey microbial eukaryotic communities in aquatic environments

Sequential membrane filtration of water samples is commonly used to monitor the diversity of aquatic microbial eukaryotes. This capture method is efficient to focus on specific taxonomic groups within a size fraction, but it is time-consuming. Centrifugation, often used to collect microorganisms from pure culture, could be seen as an alternative to capture microbial eukaryotic communities from environmental samples. Here, we compared the two capture methods to assess diversity and ecological patterns of eukaryotic communities in the Thau lagoon, France. Water samples were taken twice a month over a full year and sequential filtration targeting the picoplankton (0.2-3 μm) and larger organisms (>3 μm) was used in parallel to centrifugation. The microbial eukaryotic community in the samples was described using an environmental DNA approach targeting the V4 region of the 18S rRNA gene. The most abundant divisions in the filtration fractions and the centrifugation pellet were Dinoflagellata, Metazoa, Ochrophyta, Cryptophyta. Chlorophyta were dominant in the centrifugation pellet and the picoplankton fraction but not in the larger fraction. Diversity indices and structuring patterns of the community in the two size fractions and the centrifugation pellet were comparable. Twenty amplicon sequence variants were significantly differentially abundant between the two size fractions and the centrifugation pellet, and their temporal patterns of abundance in the two fractions combined were similar to those obtained with centrifugation. Overall, centrifugation led to similar ecological conclusions as the two filtrated fractions combined, thus making it an attractive time-efficient alternative to sequential filtration.

Covariation between microeukaryotes and bacteria associated with Planorbidae snails

Background

Microbial communities associated with macroorganisms might affect host physiology and homeostasis. Bacteria are well studied in this context, but the diversity of microeukaryotes, as well as covariations with bacterial communities, remains almost unknown.

Methods

To study microeukaryotic communities associated with Planorbidae snails, we developed a blocking primer to reduce amplification of host DNA during metabarcoding analyses. Analyses of alpha and beta diversities were computed to describe microeukaryotes and bacteria using metabarcoding of 18S and 16S rRNA genes, respectively.

Results

Only three phyla (Amoebozoa, Opisthokonta and Alveolata) were dominant for microeukaryotes. Bacteria were more diverse with five dominant phyla (Proteobacteria, Bacteroidetes, Tenericutes, Planctomycetes and Actinobacteria). The composition of microeukaryotes and bacteria were correlated for the Biomphalaria glabrata species, but not for Planorbarius metidjensis. Network analysis highlighted clusters of covarying taxa. Among them, several links might reflect top-down control of bacterial populations by microeukaryotes, but also possible competition between microeukaryotes having opposite distributions (Lobosa and Ichthyosporea). The role of these taxa remains unknown, but we believe that the blocking primer developed herein offers new possibilities to study the hidden diversity of microeukaryotes within snail microbiota, and to shed light on their underestimated interactions with bacteria and hosts.

Mapping harmful microalgal species by eDNA monitoring: A large-scale survey across the southwestern South China Sea

A large-scale sampling was undertaken during a research cruise across the South China Sea in August 2016, covering an area of about 100,000 km2 to investigate the molecular diversity and distributions of micro-eukaryotic protists, with a focus on the potentially harmful microalgal (HAB) species along the east coast of Peninsular Malaysia. Environmental DNAs from 30 stations were extracted and DNA metabarcoding targeting the V4 and V9 markers in the 18S rDNA was performed. Many protistan molecular units, including previously unreported HAB taxa, were discovered for the first time in the water. Our findings also revealed interesting spatial distribution patterns, with a marked signal of compositional turnover between latitudinal regimes of water masses, where dinophytes and diatom compositions were among the most strongly enhanced at the fronts, leading to distinct niches. Our results further confirmed the widespread distribution of HAB species, such as the toxigenic Alexandrium tamiyavaichii and Pseudo-nitzschia species, and the fish-killing Margalefidinium polykrikoides and Karlodinium veneficum. The molecular information obtained from this study provides an updated HAB species inventory and a toolset that could facilitate existing HAB monitoring schemes in the region to better inform management decisions.

Microclimate shapes the phylosymbiosis of rodent gut microbiota in Jordan's Great Rift Valley

Host phylogeny and the environment play vital roles in shaping animal microbiomes. However, the effects of these variables on the diversity and richness of the gut microbiome in different bioclimatic zones remain underexplored. In this study, we investigated the effects of host phylogeny and bioclimatic zone on the diversity and composition of the gut microbiota of two heterospecific rodent species, the spiny mouse Acomys cahirinus and the house mouse Mus musculus, in three bioclimatic zones of the African Great Rift Valley (GRV). We confirmed host phylogeny using the D-loop sequencing method and analyzed the influence of host phylogeny and bioclimatic zone parameters on the rodent gut microbiome using high-throughput amplicon sequencing of 16S rRNA gene fragments. Phylogenetic analysis supported the morphological identification of the rodents and revealed a marked genetic difference between the two heterospecific species. We found that bioclimatic zone had a significant effect on the gut microbiota composition while host phylogeny did not. Microbial alpha diversity of heterospecific hosts was highest in the Mediterranean forest bioclimatic zone, followed by the Irano-Turanian shrubland, and was lowest in the Sudanian savanna tropical zone. The beta diversity of the two rodent species showed significant differences across the Mediterranean, Irano-Turanian, and Sudanian regions. The phyla Firmicutes and Bacteroidetes were highly abundant, and Deferribacterota, Cyanobacteria and Proteobacteria were also prominent. Amplicon sequence variants (ASVs) were identified that were unique to the Sudanian bioclimatic zone. The core microbiota families recovered in this study were consistent among heterospecific hosts. However, diversity decreased in conspecific host populations found at lower altitudes in Sudanian bioclimatic zone. The composition of the gut microbiota is linked to the adaptation of the host to its environment, and this study underscores the importance of incorporating climatic factors such as elevation and ambient temperature, in empirical microbiome research and is the first to describe the rodent gut microbiome from the GRV.

A metabarcode based (species) inventory of the northern Adriatic phytoplankton

Background

The northern Adriatic is characterised as the coldest and most productive marine area of the Mediterranean, which is due to high nutrient levels introduced by river discharges, the largest of which is the Italian Po River (at the same time also the largest freshwater input into the Mediterranean). The northern Adriatic is a very shallow marine ecosystem with ocean current patterns that result in long retention times of plankton in the area. The northern Adriatic phytoplankton biodiversity and abundance are well-studied, through many scientific and long-term monitoring reports. These datasets were based on phytoplankton morphological traits traditionally obtained with light microscopy. The most recent comprehensive eastern Adriatic phytoplankton checklist was published more than 20 years ago and is still valuable today. Since phytoplankton taxonomy and systematics are constantly being reviewed (partly also due to new molecular methods of species identification that complement classical methodologies), checklists need to be updated and complemented. Today, metabarcoding of molecular markers gains more and more importance in biodiversity research and monitoring. Here, we report the use of high throughput sequencing methods to re-examine taxonomic richness and provide updated knowledge of phytoplankton diversity in the eastern northern Adriatic to complement the standardised light microscopy method.

New information

This study aimed to report an up-to-date list of the phytoplankton taxonomic richness and phylogenetic relationships in the eastern northern Adriatic, based on sequence variability of barcoding genes resolved with advanced molecular tools, namely metabarcoding. Here, metabarcoding is used to complement standardised light microscopy to advance conventional monitoring and research of phytoplankton communities for the purpose of assessing biodiversity and the status of the marine environments. Monthly two-year net sampling targeted six phytoplankton groups including Bacillariophyceae (diatoms) and Chrysophyceae (golden algae) belonging to Ochrophyta, Dinophyceae (dinoflagellates), Cryptophyceae (cryptophytes), Haptophyta (mostly coccolithophorids) and Chlorophyta with Prasinophyceae (prasinophytes) and Chlorophyceae (protist green algae). Generated sequence data were taxonomically assigned and redistributed in two kingdoms, five classes, 32 orders, 49 families and 67 genera. The most diverse group were dinoflagellates, comprising of 34 found genera (48.3%), following by diatoms with 23 (35.4%) and coccolithophorids with three genera (4.0%). In terms of genetic diversity, results were a bit different: a great majority of sequences with one nucleotide tolerance (ASVs, Amplicon sequence variants) assigned to species or genus level were dinoflagellates (83.8%), 13.7% diatoms and 1.6% Chlorophyta, respectively. Although many taxa have not been detected that have been considered as common in this area, metabarcoding revealed five diatoms and 20 dinoflagellate genera that were not reported in previous checklists, along with a few species from other targeted groups that have been reported previously. We here describe the first comprehensive 18S metabarcode inventory for the northern Adriatic Sea.

Intragenomic diversity of the V9 hypervariable domain in eukaryotes has little effect on metabarcoding

Metabarcoding revolutionized our understanding of diversity and ecology of microorganisms in different habitats. However, it is also associated with several inherent biases, one of which is associated with intragenomic diversity of a molecular barcode. Here, we compare intragenomic variability of the V9 region of the 18S rRNA gene in 19 eukaryotic phyla abundant in marine plankton. The level of intragenomic variability is comparable across all the phyla, and in most genomes and transcriptomes one V9 sequence and one OTU is predominant. However, most of the variability observed at the barcode level is probably caused by sequencing errors and is mitigated by using a denoising tool, DADA2. The SWARM algorithm commonly used in metabarcoding studies is not optimal for collapsing genuine and erroneous sequences into a single OTU, leading to an overestimation of diversity in metabarcoding data. For an unknown reason, SWARM inflates diversity of eupelagonemids more than that of other eukaryotes.

Identification of the heterotrophic nanoflagellate Bilabrum latius in the southern Adriatic (Mediterranean Sea)

Heterotrophic flagellates (HF) represent an important protist group in marine ecosystem functioning. Characterised by high taxonomic diversity, identification and classification of HF is often difficult using classical methods of light microscopy (LM). Complementing LM with molecular methods, such as barcoding, enables reliable taxonomic identification of even small size nanoflagellates that share similar or unnoticeable morphological features. The order Bicosoecida is a group of heterotrophic nanoflagellates that are important part of protist plankton and benthic communities of the world oceans. Recently, on the basis of high-resolution light microscopy and barcoding, a new bicosoecid genus, Bilabrum, was described with B. latius sp. as a type species. Our study reports on identification of B. latius from co-culture with the diatom species Chaetoceros affinis isolated from fresh plankton samples collected in the southern Adriatic. This detection of the Adriatic B.latius represents first record of this species outside itś up to now known and described habitat (deep-sea sediment of the South - East Atlantic Ocean) and in diatom co-culture.

Genetic Markers for Metabarcoding of Freshwater Microalgae: Review

The metabarcoding methods for studying the diversity of freshwater microalgae and routine biomonitoring are actively used in modern research. A lot of experience has been accumulated already, and many methodological questions have been solved (such as the influence of the methods and time of sample conservation, DNA extraction and bioinformatical processing). The reproducibility of the method has been tested and confirmed. However, one of the main problems-choosing a genetic marker for the study-still lacks a clear answer. We analyzed 70 publications and found out that studies on eukaryotic freshwater microalgae use 12 markers (different nuclear regions 18S and ITS and plastids rbcL, 23S and 16S). Each marker has its peculiarities; they amplify differently and have various levels of efficiency (variability) in different groups of algae. The V4 and V9 18S and rbcL regions are used most often. We concentrated especially on the studies that compare the results of using different markers and microscopy. We summarize the data on the primers for each region and on how the choice of a marker affects the taxonomic composition of a community.

Community dynamics and co-occurrence relationships of pelagic ciliates and their potential prey at a coastal and an offshore station in the ultra-oligotrophic Eastern Mediterranean Sea

Ciliates have been recognized as one of the major components of the microbial food web, especially in ultra-oligotrophic waters, such as the Eastern Mediterranean Sea, where nutrients are scarce and the microbial community is dominated by pico- and nano-sized organisms. For this reason, ciliates play an important role in these ecosystems since they are the main planktonic grazers. Regardless the importance of these organisms, little is known about the community structure of heterotrophic and mixotrophic ciliates and how they are associated to their potential prey. In this study, we used 18S V4 rRNA gene metabarcoding to analyze ciliate community dynamics and how the relationship with potential prey changes according to different seasons and depths. Samples were collected seasonally at two stations of the Eastern Mediterranean Sea (HCB: coastal, M3A: offshore) from the surface and deep chlorophyll maximum (DCM) layers. The ciliate community structure varied across depths in HCB and across seasons in M3A, and the network analysis showed that in both stations, mixotrophic oligotrichs were positively associated with diatoms and showed few negative associations with ASVs annotated as marine Stramenopiles (MAST). On the other hand, heterotrophic tintinnids showed negative relationships in both HCB and M3A stations, mostly with Ochrophyta and Chlorophyta. These results showed, in first place that, although the two stations are close to each other, the ciliate dynamics differed between them. Moreover, mixotrophic and heterotrophic ciliates may have different ecological niches since mixotrophic ciliates may be more selective compared to heterotrophic species regarding their prey. These findings are the first glimpse into an understanding of the dynamics between heterotrophic and mixotrophic ciliates and their role in microbial assemblages and dynamics of ultra-oligotrophic environments.

Trajectories of freshwater microbial genomics and greenhouse gas saturation upon glacial retreat

Due to climate warming, ice sheets around the world are losing mass, contributing to changes across terrestrial landscapes on decadal time spans. However, landscape repercussions on climate are poorly constrained mostly due to limited knowledge on microbial responses to deglaciation. Here, we reveal the genomic succession from chemolithotrophy to photo- and heterotrophy and increases in methane supersaturation in freshwater lakes upon glacial retreat. Arctic lakes at Svalbard also revealed strong microbial signatures form nutrient fertilization by birds. Although methanotrophs were present and increased along lake chronosequences, methane consumption rates were low even in supersaturated systems. Nitrous oxide oversaturation and genomic information suggest active nitrogen cycling across the entire deglaciated landscape, and in the high Arctic, increasing bird populations serve as major modulators at many sites. Our findings show diverse microbial succession patterns, and trajectories in carbon and nitrogen cycle processes representing a positive feedback loop of deglaciation on climate warming.

The distribution and diversity of eukaryotic phytoplankton in the Icelandic marine environment

Phytoplankton play a crucial role in the marine food web and are sensitive indicators of environmental change. Iceland is at the center of a contrasting hydrography, with cold Arctic water coming in from the north and warmer Atlantic water from the south, making this geographical location very sensitive to climate change. We used DNA metabarcoding to determine the biogeography of phytoplankton in this area of accelerating change. Seawater samples were collected in spring (2012-2018), summer (2017) and winter (2018) together with corresponding physico-chemical metadata around Iceland. Amplicon sequencing of the V4 region of the 18S rRNA gene indicates that eukaryotic phytoplankton community composition is different between the northern and southern water masses, with some genera completely absent from Polar Water masses. Emiliania was more dominant in the Atlantic-influenced waters and in summer, and Phaeocystis was more dominant in the colder, northern waters and in winter. The Chlorophyta picophytoplankton genus, Micromonas, was similarly dominant to the dominant diatom genus, Chaetoceros. This study presents an extensive dataset which can be linked with other 18s rRNA datasets for further investigation into the diversity and biogeography of marine protists in the North Atlantic.

Linking biodiversity and ecological function through extensive microeukaryotic movement across different habitats in six urban parks

Highly diverse but divergent microeukaryotes dwell in all types of habitats in urban park ecosystems. Extensive microbial migration occurs between both terrestrial and aquatic habitats. Microbial movement is beneficial to the maintenance of biodiversity and the exchange of functional guilds.

Impact of photobioreactor design on microalgae-bacteria communities grown on wastewater: Differences between thin-layer cascade and thin-layer raceway ponds

Thin-layer (TL) photobioreactors (PBRs) are characterised by high productivity. However, their use is limited to lab/pilot-scale, and a deeper level of characterisation is needed to reach industrial scale and test the resistance of multiple microalgae. Here, the performance and composition of eight microalgal communities cultivated in the two main TLs design (thin-layer cascade (TLC) and thin-layer raceway pond (RW)) were investigated through Illumina sequencing. Chlorella vulgaris showed robustness in both designs and often acted as an "invasive" species. Inoculum and reactor type brought variability. Eukaryotic microalgae inocula led to a more robust and stable community (higher similarity), however, RWs were characterised by a higher variability and did not favour the eukaryotic microalgae. The only cyanobacterial inoculum, Nostoc piscinale, was maintained, however the community was variable between designs. The reactor design had an effect on the N cycle with the TLC and RW configurations, enhancing nitrification and denitrification respectively.

Evaluation of the rbcL marker for metabarcoding of marine diatoms and inference of population structure of selected genera

Diatoms are one of the most important phytoplankton groups in the world's oceans. There are responsible for up to 40% of the photosynthetic activity in the Ocean, and they play an important role in the silicon and carbon cycles by decoupling carbon from atmospheric interactions through sinking and export. These processes are strongly influenced by the taxonomic composition of diatom assemblages. Traditionally, these have been assessed using microscopy, which in some cases is not reliable or reproducible. Next-generation sequencing enabled us to study diversity in a high-throughput manner and uncover new distribution patterns and diversity. However, phylogenetic markers used for this purpose, such as various 18S rDNA regions, are often insufficient because they cannot distinguish between some taxa. In this work, we demonstrate the performance of the chloroplast-encoded rbcL marker for metabarcoding marine diatoms compared to microscopy and 18S-V9 metabarcoding using a series of monthly samples from the Gulf of Trieste (GoT), northern Adriatic Sea. We demonstrate that rbcL is able to detect more taxa compared to 18S-V9 metabarcoding or microscopy, while the overall structure of the diatom assemblage was comparable to the other two methods with some variations, that were taxon dependent. In total, 6 new genera and 22 new diatom species for the study region were identified. We were able to spot misidentification of genera obtained with microscopy such as Pseudo-nitzschia galaxiae, which was mistaken for Cylindrotheca closterium, as well as genera that were completely overlooked, such as Minidiscus and several genera from the Cymatosiraceae family. Furthermore, on the example of two well-studied genera in the region, namely Chaetoceros and particularly Pseudo-nitzschia, we show how the rbcL method can be used to infer even deeper phylogenetic and ecologically significant differences at the species population level. Despite a very thorough community analysis obtained by rbcL the incompleteness of reference databases was still evident, and we shed light on possible improvements. Our work has further implications for studies dealing with taxa distribution and population structure, as well as carbon and silica flux models and networks.

Salp blooms drive strong increases in passive carbon export in the Southern Ocean

The Southern Ocean contributes substantially to the global biological carbon pump (BCP). Salps in the Southern Ocean, in particular Salpa thompsoni, are important grazers that produce large, fast-sinking fecal pellets. Here, we quantify the salp bloom impacts on microbial dynamics and the BCP, by contrasting locations differing in salp bloom presence/absence. Salp blooms coincide with phytoplankton dominated by diatoms or prymnesiophytes, depending on water mass characteristics. Their grazing is comparable to microzooplankton during their early bloom, resulting in a decrease of ~1/3 of primary production, and negative phytoplankton rates of change are associated with all salp locations. Particle export in salp waters is always higher, ranging 2- to 8- fold (average 5-fold), compared to non-salp locations, exporting up to 46% of primary production out of the euphotic zone. BCP efficiency increases from 5 to 28% in salp areas, which is among the highest recorded in the global ocean.

Viral infection switches the balance between bacterial and eukaryotic recyclers of organic matter during coccolithophore blooms

Algal blooms are hotspots of marine primary production and play central roles in microbial ecology and global elemental cycling. Upon demise of the bloom, organic carbon is partly respired and partly transferred to either higher trophic levels, bacterial biomass production or sinking. Viral infection can lead to bloom termination, but its impact on the fate of carbon remains largely unquantified. Here, we characterize the interplay between viral infection and the composition of a bloom-associated microbiome and consequently the evolving biogeochemical landscape, by conducting a large-scale mesocosm experiment where we monitor seven induced coccolithophore blooms. The blooms show different degrees of viral infection and reveal that only high levels of viral infection are followed by significant shifts in the composition of free-living bacterial and eukaryotic assemblages. Intriguingly, upon viral infection the biomass of eukaryotic heterotrophs (thraustochytrids) rivals that of bacteria as potential recyclers of organic matter. By combining modeling and quantification of active viral infection at a single-cell resolution, we estimate that viral infection causes a 2-4 fold increase in per-cell rates of extracellular carbon release in the form of acidic polysaccharides and particulate inorganic carbon, two major contributors to carbon sinking into the deep ocean. These results reveal the impact of viral infection on the fate of carbon through microbial recyclers of organic matter in large-scale coccolithophore blooms.

Metabarcoding reveals high genetic diversity of harmful algae in the coastal waters of Texas, Gulf of Mexico

Environmental-DNA (eDNA) for metabarcoding is a rapid and effective means to investigate microplankton community composition and species diversity. The objective of this study was to examine the genetic diversity of the phytoplankton community in the Gulf of Mexico, with particular emphasis on harmful algal bloom species. Samples were collected at stations along the coast of Texas in September-October 2017 that were inundated by low salinity waters in the aftermath of Hurricane Harvey. Metabarcodes were generated from the eDNA targeting both the V4 and V8-V9 regions of the 18S rDNA gene. Evaluation of the metabarcodes revealed an unexpectedly high number of harmful algal species during this short period, including five that had not been documented in this region previously. A total of 36 harmful algal species could be differentiated based on V4 and V8-V9 metabarcode markers. Using a phylogenetic approach, the taxonomic resolution of each marker differed and not all species could be differentiated using solely one marker. The V4 region resolved species within some genera (e.g., Heterocapsa), while the V8-V9 marker was necessary to resolve species within other genera (e.g., Chattonella). In other cases, species differentiation within a genus required a combination of both markers (e.g., Prorocentrum, Karenia), or another marker will be needed to resolve all species (e.g., Alexandrium, Dinophysis). We conclude that no single marker can delineate all species, so it is recommended HAB monitoring programs use more than one marker. Overall, the observed diversity of HAB species along the Texas coast using metabarcoding exceeded reports from other parts of the world.

Macroevolutionary patterns in intragenomic rDNA variability among planktonic foraminifera

Ribosomal intragenomic variability in prokaryotes and eukaryotes is a genomic feature commonly studied for its inflationary impact on molecular diversity assessments. However, the evolutionary mechanisms and distribution of this phenomenon within a microbial group are rarely explored. Here, we investigate the intragenomic variability in 33 species of planktonic foraminifera, calcifying marine protists, by inspecting 2,403 partial SSU sequences obtained from single-cell clone libraries. Our analyses show that polymorphisms are common among planktonic foraminifera species, but the number of polymorphic sites significantly differs among clades. With our molecular simulations, we could assess that most of these mutations are located in paired regions that do not affect the secondary structure of the SSU fragment. Finally, by mapping the number of polymorphic sites on the phylogeny of the clades, we were able to discuss the evolution and potential sources of intragenomic variability in planktonic foraminifera, linking this trait to the distinctive nuclear and genomic dynamics of this microbial group.

Centrate as a sustainable growth medium: Impact on microalgal inocula and bacterial communities in tubular photobioreactor cultivation systems

Centrate is a low-cost alternative to synthetic fertilizers for microalgal cultivation, reducing environmental burdens and remediation costs. Adapted microalgae need to be selected and characterised to maximise biomass production and depuration efficiency. Here, the performance and composition of six microalgal communities cultivated both on synthetic media and centrate within semi-open tubular photobioreactors were investigated through Illumina sequencing. Biomass grown on centrate, exposed to a high concentration of ammonium, showed a higher quantity of nitrogen (5.6% dry weight) than the biomass grown on the synthetic media nitrate (3.9% dry weight). Eukaryotic inocula were replaced by other microalgae while cyanobacterial inocula were maintained. Communities were generally similar for the same inoculum between media, however, inoculation with cyanobacteria led to variability within the eukaryotic community. Where communities differed, centrate resulted in a higher richness and diversity. The higher nitrogen of centrate possibly led to higher abundance of genes coding for N metabolism enzymes.

metaPR2 : A database of eukaryotic 18S rRNA metabarcodes with an emphasis on protists

In recent years, metabarcoding has become the method of choice for investigating the composition and assembly of microbial eukaryotic communities. The number of environmental data sets published has increased very rapidly. Although unprocessed sequence files are often publicly available, processed data, in particular clustered sequences, are rarely available in a usable format. Clustered sequences are reported as operational taxonomic units (OTUs) with different similarity levels or more recently as amplicon sequence variants (ASVs). This hampers comparative studies between different environments and data sets, for example examining the biogeographical patterns of specific groups/species, as well analysing the genetic microdiversity within these groups. Here, we present a newly-assembled database of processed 18S rRNA metabarcodes that are annotated with the PR² reference sequence database. This database, called metaPR² , contains 41 data sets corresponding to more than 4000 samples and 90,000 ASVs. The database, which is accessible through both a web-based interface (https://shiny.metapr2.org) and an R package, should prove very useful to all researchers working on protist diversity in a variety of systems.

The Diversity Patterns of Rare to Abundant Microbial Eukaryotes Across a Broad Range of Salinities in a Solar Saltern

Solar salterns are excellent artificial systems for examining species diversity and succession along salinity gradients. Here, the eukaryotic community in surface water of a Korean solar saltern (30 to 380 practical salinity units) was investigated from April 2019 to October 2020 using Illumina sequencing targeting the V4 and V9 regions of 18S rDNA. A total of 926 operational taxonomic units (OTUs) and 1,999 OTUs were obtained with the V4 and V9 regions, respectively. Notably, most of the OTUs were microbial eukaryotes, and the high-abundance groups (> 5% relative abundance (RA), Alveolata, Stramenopila, Archaeplastida, and Opisthokonta) usually accounted for > 90% of the total cumulative read counts and > 80% of all OTUs. Moreover, the high-abundance Alveolata (larger forms) and Stramenopila (smaller forms) groups displayed a significant inverse relationship, probably due to predator-prey interactions. Most of the low-abundance (0.1-5% RA) and rare (< 0.1% RA) groups remained small portion during the field surveys. Taxonomic novelty (at < 90% sequence identity) was high in the Amoebozoa, Cryptista, Haptista, Rhizaria, and Stramenopila groups (69.8% of all novel OTUs), suggesting the presence of a large number of hidden species in hypersaline environments. Remarkably, the high-abundance groups had little overlap with the other groups, implying the weakness of rare-to-prevalent community dynamics. The low-abundance Discoba group alone temporarily became the high-abundance group, suggesting that it is an opportunistic group. Overall, the composition and diversity of the eukaryotic community in hypersaline environments may be persistently stabilized, despite diverse disturbance events.

Spatial and biological oceanographic insights into the massive fish-killing bloom of the haptophyte Chrysochromulina leadbeateri in northern Norway

A bloom of the fish-killing haptophyte Chrysochromulina leadbeateri in northern Norway during May and June 2019 was the most harmful algal event ever recorded in the region, causing massive mortalities of farmed salmon. Accordingly, oceanographic and biodiversity aspects of the bloom were studied in unprecedented detail, based on metabarcoding and physico-chemical and biotic factors related with the dynamics and distribution of the bloom. Light- and electron-microscopical observations of nanoplankton samples from diverse locations confirmed that C. leadbeateri was dominant in the bloom and the primary cause of associated fish mortalities. Cell counts by light microscopy and flow cytometry were obtained throughout the regional bloom within and adjacent to five fjord systems. Metabarcoding sequences of the V4 region of the 18S rRNA gene from field material collected during the bloom and a cultured isolate from offshore of Tromsøy island confirmed the species identification. Sequences from three genetic markers (18S, 28S rRNA gene and ITS region) verified the close if not identical genetic similarity to C. leadbeateri from a previous massive fish-killing bloom in 1991 in northern Norway. The distribution and cell abundance of C. leadbeateri and related Chrysochromulina species in the recent incident were tracked by integrating observations from metabarcoding sequences of the V4 region of the 18S rRNA gene. Metabarcoding revealed at least 14 distinct Chrysochromulina variants, including putative cryptic species. C. leadbeateri was by far the most abundant of these species, but with high intraspecific genetic variability. Highest cell abundance of up to 2.7 × 10⁷ cells L ^- 1 of C. leadbeateri was found in Balsfjorden; the high cell densities were associated with stratification near the pycnocline (at ca. 12 m depth) within the fjord. The cell abundance of C. leadbeateri showed positive correlations with temperature, negative correlation with salinity, and a slightly positive correlation with ambient phosphate and nitrate concentrations. The spatio-temporal succession of the C. leadbeateri bloom suggests independent initiation from existing pre-bloom populations in local zones, perhaps sustained and supplemented over time by northeastward advection of the bloom from the fjords.

DNA sequence and taxonomic gap analyses to quantify the coverage of aquatic cyanobacteria and eukaryotic microalgae in reference databases: Results of a survey in the Alpine region

The taxonomic identification of organisms based on the amplification of specific genetic markers (metabarcoding) implicitly requires adequate discriminatory information and taxonomic coverage of environmental DNA sequences in taxonomic databases. These requirements were quantitatively examined by comparing the determination of cyanobacteria and microalgae obtained by metabarcoding and light microscopy. We used planktic and biofilm samples collected in 37 lakes and 22 rivers across the Alpine region. We focused on two of the most used and best represented genetic markers in the reference databases, namely the 16S rRNA and 18S rRNA genes. A sequence gap analysis using blastn showed that, in the identity range of 99-100%, approximately 30% (plankton) and 60% (biofilm) of the sequences did not find any close counterpart in the reference databases (NCBI GenBank). Similarly, a taxonomic gap analysis showed that approximately 50% of the cyanobacterial and eukaryotic microalgal species identified by light microscopy were not represented in the reference databases. In both cases, the magnitude of the gaps differed between the major taxonomic groups. Even considering the species determined under the microscope and represented in the reference databases, 22% and 26% were still not included in the results obtained by the blastn at percentage levels of identity ≥95% and ≥97%, respectively. The main causes were the absence of matching sequences due to amplification and/or sequencing failure and potential misidentification in the microscopy step. Our results quantitatively demonstrated that in metabarcoding the main obstacles in the classification of 16S rRNA and 18S rRNA sequences and interpretation of high-throughput sequencing biomonitoring data were due to the existence of important gaps in the taxonomic completeness of the reference databases and the short length of reads. The study focused on the Alpine region, but the extent of the gaps could be much greater in other less investigated geographic areas.

Towards quantitative metabarcoding of eukaryotic plankton: an approach to improve 18S rRNA gene copy number bias

Plankton metabarcoding is increasingly implemented in marine ecosystem assessments and is more cost-efficient and less time-consuming than monitoring based on microscopy (morphological). 18S rRNA gene is the most widely used marker for groups’ and species’ detection and classification within marine eukaryotic microorganisms. These datasets have commonly relied on the acquisition of organismal abundances directly from the number of DNA sequences (i.e. reads). Besides the inherent technical biases in metabarcoding, the largely varying 18S rRNA gene copy numbers (GCN) among marine protists (ranging from tens to thousands) is one of the most important biological biases for species quantification. In this work, we present a gene copy number correction factor (CF) for four marine planktonic groups: Bacillariophyta, Dinoflagellata, Ciliophora miscellaneous and flagellated cells. On the basis of the theoretical assumption that ‘1 read’ is equivalent to ‘1 GCN’, we used the GCN median values per plankton group to calculate the corrected cell number and biomass relative abundances. The species-specific absolute GCN per cell were obtained from various studies published in the literature. We contributed to the development of a species-specific 18S rRNA GCN database proposed by previous authors. To assess the efficiency of the correction factor we compared the metabarcoding, morphological and corrected relative abundances (in cell number and biomass) of 15 surface water samples collected in the Belgian Coastal Zone. Results showed that the application of the correction factor over metabarcoding results enables us to significantly improve the estimates of cell abundances for Dinoflagellata, Ciliophora and flagellated cells, but not for Bacillariophyta. This is likely to due to large biovolume plasticity in diatoms not corresponding to genome size and gene copy numbers. C-biomass relative abundance estimations directly from amplicon reads were only improved for Dinoflagellata and Ciliophora. The method is still facing biases related to the low number of species GCN assessed. Nevertheless, the increase of species in the GCN database may lead to the refinement of the proposed correction factor.

Surface microbiota of Mediterranean loggerhead sea turtles unraveled by 16S and 18S amplicon sequencing

The loggerhead sea turtle is considered a keystone species with a major ecological role in Mediterranean marine environment. As is the case with other wild reptiles, their outer microbiome is rarely studied. Although there are several studies on sea turtle’s macro-epibionts and endo-microbiota, there has been little research on epibiotic microbiota associated with turtle skin and carapace. Therefore we aimed to provide the identification of combined epibiotic eukaryotic, bacterial and archaeal microbiota on Mediterranean loggerhead sea turtles. In this study, we sampled skins and carapaces of 26 loggerheads from the Mediterranean Sea during 2018 and 2019. To investigate the overall microbial diversity and composition, amplicon sequencing of 16S and 18S rRNA genes was performed. We found that the Mediterranean loggerhead sea turtle epibiotic microbiota is a reservoir of a vast variety of microbial species. Microbial communities mostly varied by different locations and seas, while within bacterial communities’ significant difference was observed between sampled body sites (carapace vs. skin). In terms of relative abundance, Proteobacteria and Bacteroidota were the most represented phyla within prokaryotes, while Alveolata and Stramenopiles thrived among eukaryotes. This study, besides providing a first survey of microbial eukaryotes on loggerheads via metabarcoding, identifies fine differences within both bacterial and eukaryotic microbial communities that seem to reflect the host anatomy and habitat. Multi-domain epi-microbiome surveys provide additional layers of information that are complementary with previous morphological studies and enable better understanding of the biology and ecology of these vulnerable marine reptiles.

Composition and spatial-temporal dynamics of phytoplankton community shaped by environmental selection and interactions in the Jiaozhou Bay

The Jiaozhou Bay as a model marine ecosystem in China has been intensively investigated over the last 90 years. However, detailed phytoplankton community composition, spatial-temporal dynamics, and its assembly mechanism were still unclear. To address these, we systematically examined the composition and spatial-temporal dynamics of phytoplankton in the Jiaozhou Bay through high-throughput sequencing of 18S rDNA V4. Analysis of 468 samples from 12 sampling sites over one full year revealed much higher phytoplankton diversity than previous reports, and strong seasonal succession patterns. Some phytoplankton also showed spatial variations, although the phytoplankton community didn't show significant distance-decay pattern. Environmental factors (especially temperature), species-species interrelationships and unique resting stages were uncovered to be the main contributors instead of stochastic process in shaping the phytoplankton community assembly. The overwhelming positive correlations between phytoplankton and other protists suggested that coevolution might be critical in this marine ecosystem. Complementary distributions of different amplicon sequence variants (ASVs) of same genera, such as Skeletonema marinoi (ASV_2) and Skeletonema tropicum (ASV_263) of the genus Skeletonema, suggested that phytoplankton have evolved differentially to exploit a wide range of ecological niches. This study laid a solid foundation for asertaining phytoplankton composition and spatial-temporal dynamics in temperate seawaters and mechanisms underlying phytoplankton community assembly, allowing in-depth studies of marine ecology.

The Microbiological Drivers of Temporally Dynamic Dimethylsulfoniopropionate Cycling Processes in Australian Coastal Shelf Waters

The organic sulfur compounds dimethylsulfoniopropionate (DMSP) and dimethyl sulfoxide (DMSO) play major roles in the marine microbial food web and have substantial climatic importance as sources and sinks of dimethyl sulfide (DMS). Seasonal shifts in the abundance and diversity of the phytoplankton and bacteria that cycle DMSP are likely to impact marine DMS (O) (P) concentrations, but the dynamic nature of these microbial interactions is still poorly resolved. Here, we examined the relationships between microbial community dynamics with DMS (O) (P) concentrations during a 2-year oceanographic time series conducted on the east Australian coast. Heterogenous temporal patterns were apparent in chlorophyll a (chl a) and DMSP concentrations, but the relationship between these parameters varied over time, suggesting the phytoplankton and bacterial community composition were affecting the net DMSP concentrations through differential DMSP production and degradation. Significant increases in DMSP were regularly measured in spring blooms dominated by predicted high DMSP-producing lineages of phytoplankton (Heterocapsa, Prorocentrum, Alexandrium, and Micromonas), while spring blooms that were dominated by predicted low DMSP-producing phytoplankton (Thalassiosira) demonstrated negligible increases in DMSP concentrations. During elevated DMSP concentrations, a significant increase in the relative abundance of the key copiotrophic bacterial lineage Rhodobacterales was accompanied by a three-fold increase in the gene, encoding the first step of DMSP demethylation (dmdA). Significant temporal shifts in DMS concentrations were measured and were significantly correlated with both fractions (0.2-2 μm and >2 μm) of microbial DMSP lyase activity. Seasonal increases of the bacterial DMSP biosynthesis gene (dsyB) and the bacterial DMS oxidation gene (tmm) occurred during the spring-summer and coincided with peaks in DMSP and DMSO concentration, respectively. These findings, along with significant positive relationships between dsyB gene abundance and DMSP, and tmm gene abundance with DMSO, reinforce the significant role planktonic bacteria play in producing DMSP and DMSO in ocean surface waters. Our results highlight the highly dynamic nature and myriad of microbial interactions that govern sulfur cycling in coastal shelf waters and further underpin the importance of microbial ecology in mediating important marine biogeochemical processes.

Seasonal dynamics of marine protist communities in tidally mixed coastal waters

Major seasonal community reorganizations and associated biomass variations are landmarks of plankton ecology. However, the processes of plankton community turnover rates have not been fully elucidated so far. Here, we analyse patterns of planktonic protist community succession in temperate latitudes, based on quantitative taxonomic data from both microscopy counts (cells >10 μm) and ribosomal DNA metabarcoding (size fraction >3 μm, 18S rRNA gene) from plankton samples collected bimonthly over 8 years (2009–2016) at the SOMLIT‐Astan station (Roscoff, Western English Channel). Based on morphology, diatoms were clearly the dominating group all year round and over the study period. Metabarcoding uncovered a wider diversity spectrum and revealed the prevalence of Dinophyceae and diatoms but also of Cryptophyta, Chlorophyta, Cercozoa, Syndiniales and Ciliophora in terms of read counts and or richness. The use of morphological and molecular analyses in combination allowed improving the taxonomic resolution and to identify the sequence of the dominant species and OTUs (18S V4 rDNA‐derived taxa) that drive annual plankton successions. We detected that some of these dominant OTUs were benthic as a result of the intense tidal mixing typical of the French coasts in the English Channel. Our analysis of the temporal structure of community changes point to a strong seasonality and resilience. The temporal structure of environmental variables (especially Photosynthetic Active Radiation, temperature and macronutrients) and temporal structures generated by species life cycles and or species interactions, are key drivers of the observed cyclic annual plankton turnover.

Diversity of Microbial Eukaryotes Along the West Antarctic Peninsula in Austral Spring

During a cruise from October to November 2019, along the West Antarctic Peninsula, between 64.32 and 68.37°S, we assessed the diversity and composition of the active microbial eukaryotic community within three size fractions: micro- (> 20 μm), nano- (20-5 μm), and pico-size fractions (5-0.2 μm). The communities and the environmental parameters displayed latitudinal gradients, and we observed a strong similarity in the microbial eukaryotic communities as well as the environmental parameters between the sub-surface and the deep chlorophyll maximum (DCM) depths. Chlorophyll concentrations were low, and the mixed layer was shallow for most of the 17 stations sampled. The richness of the microplankton was higher in Marguerite Bay (our southernmost stations), compared to more northern stations, while the diversity for the nano- and pico-plankton was relatively stable across latitude. The microplankton communities were dominated by autotrophs, mostly diatoms, while mixotrophs (phototrophs-consuming bacteria and kleptoplastidic ciliates, mostly alveolates, and cryptophytes) were the most abundant and active members of the nano- and picoplankton communities. While phototrophy was the dominant trophic mode, heterotrophy (mixotrophy, phagotrophy, and parasitism) tended to increase southward. The samples from Marguerite Bay showed a distinct community with a high diversity of nanoplankton predators, including spirotrich ciliates, and dinoflagellates, while cryptophytes were observed elsewhere. Some lineages were significantly related-either positively or negatively-to ice coverage (e.g., positive for Pelagophyceae, negative for Spirotrichea) and temperature (e.g., positive for Cryptophyceae, negative for Spirotrichea). This suggests that climate changes will have a strong impact on the microbial eukaryotic community.

Short- and long-read metabarcoding of the eukaryotic rRNA operon: evaluation of primers and comparison to shotgun metagenomics sequencing

High-throughput sequencing for analysis of microbial diversity has evolved vastly over the last decade. Currently the go-to method for studying microbial eukaryotes is short-read metabarcoding of variable regions of the 18S rRNA gene with <500 bp amplicons. However, there is a growing interest in applying long-read sequencing of amplicons covering the rRNA operon for improving taxonomic resolution. For both methods, the choice of primers is crucial. It determines if community members are covered, if they can be identified at a satisfactory taxonomic level, and if the obtained community profile is representative. Here, we designed new primers targeting 18S and 28S rRNA based on 177,934 and 21,072 database sequences, respectively. The primers were evaluated in silico along with published primers on reference sequence databases and marine metagenomics datasets. We further evaluated a subset of the primers for short- and long-read sequencing on environmental samples in vitro and compared the obtained community profile with primer-unbiased metagenomic sequencing. Of the short-read pairs, a new V6-V8 pair and the V4_Balzano pair used with a simplified PCR protocol provided good results in silico and in vitro. Fewer differences were observed between the long-read primer pairs. The long-read amplicons and ITS1 alone provided higher taxonomic resolution than V4. Together, our results represent a reference and guide for selection of robust primers for research on and environmental monitoring of microbial eukaryotes.

Photoperiod-driven rhythms reveal multi-decadal stability of phytoplankton communities in a highly fluctuating coastal environment

Phytoplankton play a pivotal role in global biogeochemical and trophic processes and provide essential ecosystem services. However, there is still no broad consensus on how and to what extent their community composition responds to environmental variability. Here, high-frequency oceanographic and biological data collected over more than 25 years in a coastal Mediterranean site are used to shed light on the temporal patterns of phytoplankton species and assemblages in their environmental context. Because of the proximity to the coast and due to large-scale variations, environmental conditions showed variability on the short and long-term scales. Nonetheless, an impressive regularity characterised the annual occurrence of phytoplankton species and their assemblages, which translated into their remarkable stability over decades. Photoperiod was the dominant factor related to community turnover and replacement, which points at a possible endogenous regulation of biological processes associated with species-specific phenological patterns, in analogy with terrestrial plants. These results highlight the considerable stability and resistance of phytoplankton communities in response to different environmental pressures, which contrast the view of these organisms as passively undergoing changes that occur at different temporal scales in their habitat, and show how, under certain conditions, biological processes may prevail over environmental forcing.

Biodiversity exploration in autumn using environmental DNA in the South China sea

The South China Sea (SCS) is an important part of the Indo-Pacific convergence zone, with high biodiversity and abundant marine resources. Traditional methods are primarily used to monitor biodiversity. However, a few studies have used environmental DNA (eDNA) metabarcoding to research the assemblage structure of the SCS. This study used eDNA metabarcoding to survey the SCS assemblage and its relationship with environmental factors over a month-long time-series (August 30th to September 30th, 2020) of seawater samples from the central part of the SCS (9°-20°86' N, 113°-118°47' E). 32 stations were divided into six groups (A, B, C, D, E, F) according to longitude. We collected water samples, extracted eDNA, and amplified 18S rRNA gene V4 region (18S V4), 18S rRNA gene V9 region (18S V9), and 12S rRNA gene (12S). Krona diagrams were used to show species composition. We identified 192 phytoplankton, 104 invertebrate, and 61 fish species from 18S V4, 18S V9, and 12S, respectively. Generally, the three assemblage structures exhibited an increase in species diversity with increasing longitude. Group E had the highest fish diversity. Groups F and C had the highest phytoplankton and invertebrate diversity, respectively. Canonical correspondence analysis showed that four factors (chlorophyll a, depth, salinity, and temperature) were correlated with assemblage structure. Chlorophyll a was the main environmental factor that affected fish, phytoplankton, and invertebrate assemblage structures; salinity was strongly correlated with fish and invertebrate assemblage structures; temperature was a key factor that impacted fish and invertebrate assemblage structures; and depth was strongly correlated with invertebrate assemblage structure. Our results revealed that eDNA metabarcoding is a powerful tool for improving detection rate and using multiple markers is an effective approach for monitoring biodiversity. This study provided information that can be used to enhance biodiversity protection efforts in the SCS.

Biogeographical and seasonal dynamics of the marine Roseobacter community and ecological links to DMSP-producing phytoplankton

Ecological interactions between marine bacteria and phytoplankton play a pivotal role in governing the ocean's major biogeochemical cycles. Among these, members of the marine Roseobacter Group (MRG) can establish mutualistic relationships with phytoplankton that are, in part, maintained by exchanges of the organosulfur compound, dimethylsulfoniopropionate (DMSP). Yet most of what is known about these interactions has been derived from culture-based laboratory studies. To investigate temporal and spatial co-occurrence patterns between members of the MRG and DMSP-producing phytoplankton we analysed 16S and 18S rRNA gene amplicon sequence variants (ASVs) derived from 5 years of monthly samples from seven environmentally distinct Australian oceanographic time-series. The MRG and DMSP-producer communities often displayed contemporaneous seasonality, which was greater in subtropical and temperate environments compared to tropical environments. The relative abundance of both groups varied latitudinally, displaying a poleward increase, peaking (MRG at 33% of total bacteria, DMSP producers at 42% of eukaryotic phototrophs) during recurrent spring-summer phytoplankton blooms in the most temperate site (Maria Island, Tasmania). Network analysis identified 20,140 significant positive correlations between MRG ASVs and DMSP producers and revealed that MRGs exhibit significantly stronger correlations to high DMSP producers relative to other DMSP-degrading bacteria (Pelagibacter, SAR86 and Actinobacteria). By utilising the power of a continental network of oceanographic time-series, this study provides in situ confirmation of interactions found in laboratory studies and demonstrates that the ecological dynamics of an important group of marine bacteria are shaped by the production of an abundant and biogeochemically significant organosulfur compound.

Plankton diversity in Anthropocene: Shipping vs. aquaculture along the eastern Adriatic coast assessed through DNA metabarcoding

Coastal ecosystems globally are exposed to the most pervasive anthropogenic activities, caused by a suite of human infrastructure and enterprises such as shipping ports, aquaculture facilities, fishing, and tourism. These anthropogenic activities may lead to changes in ecosystem biodiversity, followed by loss of ecosystem functioning and services. Shipping industry and aquaculture have also been recognized as the main vectors for introduction of marine non-indigenous species (NIS) worldwide. In this study, we used DNA metabarcoding-based methods to investigate plankton biodiversity under varying anthropogenic pressures (shipping and bivalve aquaculture) along the eastern Adriatic coast (the northernmost part of the Mediterranean Sea). Our comparative assessment revealed similar community structures among investigated coastal locations (Northern, Central and Southern Adriatic). When the whole plankton communities were considered, they did not differ significantly between port and aquaculture sites. However, the proportion of the unique zOTUs in the port samples was remarkably higher than that in aquaculture sites (40.5% vs 8.2%), indicating that port areas may receive higher abundance and species richness of NIS than aquaculture sites. Further important difference between the two types of anthropogenically impacted habitats was a high abundance of three notorious invaders - M. leidyi, M. gigas, and H. elegans in late summer at the aquaculture site in Northern Adriatic. Therefore, the plankton community of the area is under pressure not only from aquaculture activities, but also establishment of NIS. Port areas are probably under greater introduction pressure from NIS, but aquaculture sites may experience greater community changes due to their establishment.

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.