Reconstructing marine plankton food web interactions using DNA metabarcoding

Metabarcoding reveals potentially mixotrophic flagellates and picophytoplankton as key groups of phytoplankton in the Elbe estuary

In estuaries, phytoplankton are faced with strong environmental forcing (e.g. high turbidity, salinity gradients). Taxa that appear under such conditions may play a critical role in maintaining food webs and biological carbon pumping, but knowledge about estuarine biota remains limited. This is also the case in the Elbe estuary where the lower 70 km of the water body are largely unexplored. In the present study, we investigated the phytoplankton composition in the Elbe estuary via metabarcoding. Our aim was to identify key taxa in the unmonitored reaches of this ecosystem and compare our results from the monitored area with available microscopy data. Phytoplankton communities followed distinct seasonal and spatial patterns. Community composition was similar across methods. Contributions of key classes and genera were correlated to each other (p < 0.05) when obtained from reads and biovolume (R2 = 0.59 and 0.33, respectively). Centric diatoms (e.g. Stephanodiscus) were the dominant group - comprising on average 55 % of the reads and 66-69 % of the biovolume. However, results from metabarcoding imply that microscopy underestimates the prevalence of picophytoplankton and flagellates with a potential for mixotrophy (e.g. cryptophytes). This might be due to their small size and sensitivity to fixation agents. We argue that mixotrophic flagellates are ecologically relevant in the mid to lower estuary, where, e.g., high turbidity render living conditions rather unfavorable, and skills such as phagotrophy provide fundamental advantages. Nevertheless, further findings - e.g. important taxa missing from the metabarcoding dataset - emphasize potential limitations of this method and quantitative biases can result from varying numbers of gene copies in different taxa. Further research should address these methodological issues but also shed light on the causal relationship of taxa with the environmental conditions, also with respect to active mixotrophic behavior.

Environmental heterogeneity caused by large-scale cultivation of Pyropia haitanensis shapes multi-group biodiversity distribution in coastal areas

The response of marine biodiversity to mariculture has long been a research focus in marine ecology. However, the effects of seaweed cultivation on biological community assembly are poorly understood, especially in diverse communities with distinct ecological characteristics. In this study, we used environmental DNA metabarcoding to investigate the spatial distribution patterns of bacterial, protistan, and metazoan diversity, aiming to reveal the mechanisms of community assembly in the Pyropia haitanensis cultivation zone along the Fujian coast, China. We found that, compared with the biological communities in control zones, those in P. haitanensis cultivation zones exhibited stronger geographic distance-decay patterns and displayed more complex and stable network structures. Deterministic processes (environmental selection) played a more important role in the assembly of bacterial, protistan, and metazoan communities in P. haitanensis cultivation zones, especially metazoan communities. Variance partitioning analysis showed that environmental variables made greater contributions to the diversity of the three types of communities within the P. haitanensis cultivation zones than in the control zones. Partial least squares path modeling analysis identified nitrate‑nitrogen (NO3-N), pH, particulate organic carbon (POC), and dissolved organic carbon (DOC) as the key environmental variables affecting biodiversity. Overall, the environmental heterogeneity caused by the large-scale cultivation of P. haitanensis could be the crucial factor influencing the composition and structure of various biological communities. Our results highlight the importance of the responses of multi-group organisms to the cultivation of seaweed, and provide insights into the coexistence patterns of biodiversity at the spatial scale.

Brute force prey metabarcoding to explore the diets of small invertebrates

Prey metabarcoding has become a popular tool in molecular ecology for resolving trophic interactions at high resolution, from various sample types and animals. To date, most predator-prey studies of small-sized animals (<1 mm) have met the problem of overabundant predator DNA in dietary samples by adding blocking primers/peptide nucleic acids. These primers aim to limit the PCR amplification and detection of the predator DNA but may introduce bias to the prey composition identified by interacting with sequences that are similar to those of the predator. Here we demonstrate the use of an alternative method to explore the prey of small marine copepods using whole-body DNA extracts and deep, brute force metabarcoding of an 18S rDNA fragment. After processing and curating raw data from two sequencing runs of varying depths (0.4 and 5.4 billion raw reads), we isolated 1.3 and 52.2 million prey reads, with average depths of ~15,900 and ~120,000 prey reads per copepod individual, respectively. While data from both sequencing runs were sufficient to distinguish dietary compositions from disparate seasons, locations, and copepod species, greater sequencing depth led to better separation of clusters. As computation and sequencing are becoming ever more powerful and affordable, we expect the brute force approach to become a general standard for prey metabarcoding, as it offers a simple and affordable solution to consumers that is impractical to dissect or unknown to science.

Dynamics of Gut Bacteria Across Different Zooplankton Genera in the Baltic Sea

In aquatic ecosystems, zooplankton-associated bacteria potentially have a great impact on the structure of ecosystems and trophic networks by providing various metabolic pathways and altering the ecological niche of host species. To understand the composition and drivers of zooplankton gut microbiota, we investigated the associated microbial communities of four zooplankton genera from different seasons in the Baltic Sea using the 16S rRNA gene. Among the 143 ASVs (amplified sequence variants) observed belonging to heterotrophic bacteria, 28 ASVs were shared across all zooplankton hosts over the season, and these shared core ASVs represented more than 25% and up to 60% of relative abundance in zooplankton hosts but were present at low relative abundance in the filtered water. Zooplankton host identity had stronger effects on bacterial composition than seasonal variation, with the composition of gut bacterial communities showing host-specific clustering patterns. Although bacterial compositions and dominating core bacteria were different between zooplankton hosts, higher gut bacteria diversity and more bacteria contributing to the temporal variation were found in Temora and Pseudocalanus, compared to Acartia and Synchaeta. Diet diatom and filamentous cyanobacteria negatively correlated with gut bacteria diversity, but the difference in diet composition did not explain the dissimilarity of gut bacteria composition, suggesting a general effect of diet on the inner conditions in the zooplankton gut. Synchaeta maintained high stability of gut bacterial communities with unexpectedly low bacteria-bacteria interactions as compared to the copepods, indicating host-specific regulation traits. Our results suggest that the patterns of gut bacteria dynamics are host-specific and the variability of gut bacteria is not only related to host taxonomy but also related to host behavior and life history traits.

Profiling pathogenic protozoan and their functional pathways in wastewater using 18S rRNA and shotgun metagenomics

Despite extensive research, little is known about the composition of eukaryotic protists in environmental samples. This is due to low parasite concentrations, the complexity of parasite diversity, and a lack of suitable reference databases and standardized protocols. To bridge this knowledge gap, this study used 18S rRNA short amplicon and shotgun metagenomic sequencing approaches to profile protozoan microbial communities as well as their functional pathways in treated and untreated wastewater samples collected from different regions of South Africa. Results demonstrated that protozoan diversity (Shannon index P-value = 0.03) and taxonomic composition (PERMANOVA, P-value = 0.02) was mainly driven by the type of wastewater samples (treated & untreated) and geographic location. However, these WWTPs were also found to contain a core community of protozoan parasites. The untreated wastewater samples revealed a predominant presence of free-living, parasitic, and potentially pathogenic protists typically found in humans and animals, ranging from Alveolata (27 %) phylum (Apicomplexa and Ciliophora) to Excavata (3.88 %) (Discoba and Parasalia) and Amoebozoa (2.84 %) (Entamoeba and Acanthamoeba). Shotgun metagenomics analyses in a subset of the untreated wastewater samples confirmed the presence of public health-importance protozoa, including Cryptosporidium species (3.48 %), Entamoeba hystolitica (6.58 %), Blastocystis hominis (2.91 %), Naegleria gruberi (2.37 %), Toxoplasma gondii (1.98 %), Cyclospora cayetanensis (1.30 %), and Giardia intestinalis (0.31 %). Virulent gene families linked to pathogenic protozoa, such as serine/threonine protein phosphatase and mucin-desulfating sulfatase were identified. Additionally, enriched pathways included thiamine diphosphate biosynthesis III, heme biosynthesis, Methylerythritol 4-Phosphate Pathway, methyl erythritol phosphate (MEP), and pentose phosphate pathways. These findings suggest that protozoan pathogens may possess metabolic and growth potential within WWTPs, posing a severe risk of transmission to humans and animals if inadequately disinfected before release. This study provides a baseline for the future investigation of diverse protozoal communities in wastewater, which are of public health importance.

Uptake and trophic transfer of selenium into phytoplankton and zooplankton of the southern Baltic Sea

Selenium (Se) is an essential trace element and displays a narrow range of concentration between essentiality and toxicity. Se plays an important role in ameliorating mercury toxicity in organisms. Despite this there are only a few reports concerning Se concentration in plankton, the first link in the trophic chain that determines the uptake and transfer of Se to subsequent trophic levels. This paper aimed to determine Se concentration in water, phytoplankton, and zooplankton in the Baltic Sea, and factors affecting Se absorption from the environment and its transfer to higher trophic levels. Sea water and plankton samples were collected from the Gulf of Gdańsk during 5 cruises (4 seasons: 2019-2022) at 4 research stations. An additional cruise was undertaken in July 2020 in the open waters of the southern Baltic Sea. The median Se concentrations in the Gulf of Gdańsk was 0.25 μg·dm-3. While the median of Se concentration in phytoplankton was 1.11 μg·g-1 and in zooplankton was 1.25 μg·g-1. The biomass of organisms in the phytoplankton and zooplankton in the Gulf of Gdańsk showed an important role in shaping Se concentration. Seasonal trends in Se concentration in zooplankton could be the result of taxa composition changes, changes to dietary intake of Se, changes in growth dilution, or potentially some combination of factors. The highest biomagnification rate occurred in the summer. In contrast, in autumn and winter, when plankton biomass was dominated by the ciliate species Mesodinium rubrum, the highest Se concentration in plankton was measured. Further scientific studies are needed into the active biocomponents of the Se concentration process, including Se speciation, to more fully understand the dynamics of Se concentrations in the pelagic food webs of this and other freshwater and marine systems.

Biogeographic patterns of micro-eukaryotic generalists and specialists and their effects on regional α-diversity at inter-oceanic scale

Inter-oceanic scale studies allow us to understand the global spread of micro-organisms in marine ecosystems. In this study, micro-eukaryotic communities in marine surface sediment were collected from tropical to Arctic sites. We found that micro-eukaryotic generalists had much higher intraspecific variation than specialists which allow them to distribute more widely through higher spatiotemporal asynchrony and complementary niche preferences among conspecific taxa. Moreover, comparing to the host-associated protozoa and small metazoa, the algae and free-living protozoa with higher intraspecific variation allow them to have wider distribution ranges. Species abundance also played an important role in driving the distribution ranges of generalists and specialists. The generalists had important effects on regional α-diversity even at an inter-oceanic scale which led to the micro-eukaryotic species richness in polar sites to be mainly influenced by the regional generalists but not the local specialists. In particular, more than 97% of algal species in polar sites were shared with the tropical and subtropical sites (including toxic dinoflagellate). Overall, our study suggests that the effects of global change and human activities on the vulnerable high latitude habitats may lead to biotic homogenization for the whole microbial community (not only the dispersal of some harmful algae) through the potential long-distance spread of generalists.

Role of Syndiniales parasites in depth-specific networks and carbon flux in the oligotrophic ocean

Microbial associations that result in phytoplankton mortality are important for carbon transport in the ocean. This includes parasitism, which in microbial food webs is dominated by the marine alveolate group, Syndiniales. Parasites are expected to contribute to carbon recycling via host lysis; however, knowledge on host dynamics and correlation to carbon export remain unclear and limit the inclusion of parasitism in biogeochemical models. We analyzed a 4-year 18S rRNA gene metabarcoding dataset (2016-19), performing network analysis for 12 discrete depths (1-1000 m) to determine Syndiniales-host associations in the seasonally oligotrophic Sargasso Sea. Analogous water column and sediment trap data were included to define environmental drivers of Syndiniales and their correlation with particulate carbon flux (150 m). Syndiniales accounted for 48-74% of network edges, most often associated with Dinophyceae and Arthropoda (mainly copepods) at the surface and Rhizaria (Polycystinea, Acantharea, and RAD-B) in the aphotic zone. Syndiniales were the only eukaryote group to be significantly (and negatively) correlated with particulate carbon flux, indicating their contribution to flux attenuation via remineralization. Examination of Syndiniales amplicons revealed a range of depth patterns, including specific ecological niches and vertical connection among a subset (19%) of the community, the latter implying sinking of parasites (infected hosts or spores) on particles. Our findings elevate the critical role of Syndiniales in marine microbial systems and reveal their potential use as biomarkers for carbon export.

DNA metabarcoding reveals impact of local recruitment, dispersal, and hydroperiod on assembly of a zooplankton metacommunity

Understanding the environmental impact on the assembly of local communities in relation to their spatial and temporal connectivity is still a challenge in metacommunity ecology. This study aims to unravel underlying metacommunity processes and environmental factors that result in observed zooplankton communities. Unlike most metacommunity studies, we jointly examine active and dormant zooplankton communities using a DNA metabarcoding approach to overcome limitations of morphological species identification. We applied two-fragment (COI and 18S) metabarcoding to monitor communities of 24 kettle holes over a two-year period to unravel (i) spatial and temporal connectivity of the communities, (ii) environmental factors influencing local communities, and (iii) dominant underlying metacommunity processes in this system. We found a strong separation of zooplankton communities from kettle holes of different hydroperiods (degree of permanency) throughout the season, while the community composition within single kettle holes did not differ between years. Species richness was primarily dependent on pH and permanency, while species diversity (Shannon Index) was influenced by kettle hole location. Community composition was impacted by kettle hole size and surrounding field crops. Environmental processes dominated temporal and spatial processes. Sediment communities showed a different composition compared to water samples but did not differ between ephemeral and permanent kettle holes. Our results suggest that communities are mainly structured by environmental filtering based on pH, kettle hole size, surrounding field crops, and permanency. Environmental filtering based on specific conditions in individual kettle holes seems to be the dominant process in community assembly in the studied zooplankton metacommunity.

Does environmental DNA reflect the actual phytoplankton diversity in the aquatic environment? Case study of marine mucilage in the Sea of Marmara

The present study was designed to assess the effectiveness of the eDNA metabarcoding approach to determine the phytoplankton composition in the marine environment with a special focus on mucilage episodes in the Sea of Marmara. For this purpose, the samples were collected from 5 different sites located in the Sea of Marmara and the northern Aegean Sea during the mucilage episode in June 2021. The phytoplankton diversity was analyzed morphologically and by 18S rRNA gene amplicon sequencing, and the dataset of both methods was compared, accordingly. The results showed significant differences between methods in terms of composition and the abundance of the phytoplankton groups. While Miozoa was the most abundant group by metabarcoding, light microscopy (LM) indicated a dominance of Bacillariophyta. Katablepharidophyta was found at lower abundances by the metabarcoding (representing < 1% of the community); the members of this phylum were not observed by a microscope. At the lower taxonomic levels, Chaetoceros was the only genus detected in all samples by both methods. Additionally, while mucilage-forming Gonyaulax fragilis, Cylindrotheca closterium, and Thalassiosira rotula were detected to species-level by LM, metabarcoding was able to determine these organisms at the genus level. On the other hand, the genus Arcocellulus was found in all metabarcoding datasets and not detected by microscopy. The results indicated that metabarcoding can detect a greater number of genera and reveal taxa that were overlooked by light microscopy but to develop a complete picture of phytoplankton diversity in the sample, microscopical observations still are in need.

DNA metabarcoding highlights cyanobacteria as the main source of primary production in a pelagic food web model

Models that estimate rates of energy flow in complex food webs often fail to account for species-specific prey selectivity of diverse consumer guilds. While DNA metabarcoding is increasingly used for dietary studies, methodological biases have limited its application for food web modeling. Here, we used data from dietary metabarcoding studies of zooplankton to calculate prey selectivity indices and assess energy fluxes in a pelagic resource-consumer network. We show that food web dynamics are influenced by prey selectivity and temporal match-mismatch in growth cycles and that cyanobacteria are the main source of primary production in the investigated coastal pelagic food web. The latter challenges the common assumption that cyanobacteria are not supporting food web productivity, a result that is increasingly relevant as global warming promotes cyanobacteria dominance. While this study provides a method for how DNA metabarcoding can be used to quantify energy fluxes in a marine food web, the approach presented here can easily be extended to other ecosystems.

Phytoplankton and Microzooplankton Community Structure and Assembly Mechanisms in Northwestern Pacific Ocean Estuaries with Environmental Heterogeneity and Geographic Segregation

Phytoplankton and microzooplankton are crucial players in marine ecosystems and first responders to environmental changes, but their community structures and how they are shaped by environmental conditions have rarely been studied simultaneously. In this study, we conducted an eDNA metabarcoding sequencing combined with multiple statistical methods to simultaneously analyze the phytoplankton and microzooplankton in Liaohe (LH) and Yalujiang (YLJ) estuaries. The major objective was to examine how plankton community structure and assembly mechanism may differ between two estuaries with similar latitudinal position and climate but geographical segregation and differential level of urbanization (more in LH). Clear differences in diversity and composition of phytoplankton and microzooplankton communities between LH and YLJ estuaries were observed. Richness of phytoplankton was significantly higher in LH than YLJ, while richness of microzooplankton was higher in YLJ. The magnitude of intrahabitat variations in phytoplankton communities was significantly stronger than that of microzooplankton. Some phytoplankton and microzooplankton taxa also showed interhabitat differences in their relative abundances. Phytoplankton showed a stronger geographic distance-decay of similarity than microzooplankton, while significant environmental distance-decay of similarity in microzooplankton was found in the less urbanized YLJ estuary. Community assembly of phytoplankton was, based on the neutral community models, driven primarily by stochastic processes, while deterministic processes contributed more for microzooplankton. Furthermore, we detected wider habitat niche breadths and stronger dispersal abilities in phytoplankton than in microzooplankton. These results suggest that passive dispersal shapes the phytoplankton community whereas environmental selection shapes the microzooplankton community. IMPORTANCE Understanding the underlying mechanisms shaping a metacommunity is useful to management for improving the ecosystem function. The research presented in the manuscript mainly tried to address the effects of habitat geography and environmental conditions on the phytoplankton and microzooplankton communities, and the underlying mechanisms of community assembly in temperate estuaries. In order to achieve this purpose, we developed a metabarcoding sequencing method based on 18S rRNA gene. The phytoplankton and microzooplankton communities from two estuaries with similar latitude and climatic conditions but obvious geographical segregation and significant environmental heterogeneity were investigated. The results of our study could lay a solid foundation for ascertaining phytoplankton and microzooplankton communities in estuaries with obvious environmental heterogeneity and geographic segregation and mechanisms underlying community assembly.

Opportunistic vs selective feeding strategies of zooplankton under changing environmental conditions

The plankton community consists of diverse interacting species. The estimation of species interactions in nature is challenging. There is limited knowledge on how plankton interactions are influenced by environmental conditions because of limited understanding of zooplankton feeding strategies and factors affecting trophic interactions. In this study, we used DNA-metabarcoding to investigate trophic interactions in mesozooplankton predators and the influence of prey availability on their feeding behavior. We found that mesozooplankton feeding strategies vary within species across an environmental gradient. Some species, such as Temora longicornis consistently used a selective strategy, while diets of Centropages hamatus and Acartia spp. varied between stations, showing a trophic plasticity with the prey community. We found a dominance of Synechococcales reads in Temora's gut content and a high prey diversity for the cladoceran Evadne nordmanni. Our study shows the wide range of prey species that supports mesozooplankton community and helps to understand the spatial and temporal complexity of plankton species interactions and discriminate the selectivity ability of four zooplankton key species. Due to the central role of plankton in marine waters, a better comprehension of the spatiotemporal variability in species interactions helps to estimate fluxes to benthic and pelagic predators.

The predator problem and PCR primers in molecular dietary analysis: swamped or silenced; depth or breadth?

Dietary metabarcoding has vastly improved our ability to analyse the diets of animals, but it is hampered by a plethora of technical limitations including potentially reduced data output due to the disproportionate amplification of the DNA of the focal predator, here termed 'the predator problem'. We review the various methods commonly used to overcome this problem, from deeper sequencing to exclusion of predator DNA during PCR, and how they may interfere with increasingly common multi-predator-taxon studies. We suggest that multi-primer approaches with an emphasis on achieving both depth and breadth of prey detections may overcome the issue to some extent, although multi-taxon studies require further consideration, as highlighted by an empirical example. We also review several alternative methods for reducing the prevalence of predator DNA that are conceptually promising but require additional empirical examination. The predator problem is a key constraint on molecular dietary analyses but, through this synthesis, we hope to guide researchers in overcoming this in an effective and pragmatic way.

Niche partitioning between planktivorous fish in the pelagic Baltic Sea assessed by DNA metabarcoding, qPCR and microscopy

Marine communities undergo rapid changes related to human-induced ecosystem pressures. The Baltic Sea pelagic food web has experienced several regime shifts during the past century, resulting in a system where competition between the dominant planktivorous mesopredatory clupeid fish species herring (Clupea harengus) and sprat (Sprattus sprattus) and the rapidly increasing stickleback (Gasterosteus aculeatus) population is assumed to be high. Here, we investigate diet overlap between these three planktivorous fishes in the Baltic Sea, utilizing DNA metabarcoding on the 18S rRNA gene and the COI gene, targeted qPCR, and microscopy. Our results show niche differentiation between clupeids and stickleback, and highlight that rotifers play an important role in this pattern, as a resource that is not being used by the clupeids nor by other zooplankton in spring. We further show that all the diet assessment methods used in this study are consistent, but also that DNA metabarcoding describes the plankton-fish link at the highest taxonomic resolution. This study suggests that rotifers and other understudied soft-bodied prey may have an important function in the pelagic food web and that the growing population of pelagic stickleback may be supported by the open feeding niche offered by the rotifers.

The Effects of Water Level Fluctuation on Zooplankton Communities in Shahu Lake Based on DNA Metabarcoding and Morphological Methods

Background: The water level of Poyang Lake (China) fluctuates seasonally. Shahu Lake, a smaller body of water connected to Poyang Lake during the wet season, is separated in the dry season. Due to a special fishing method termed ‘lake enclosed in autumn’, the water level is lowered and reaches its lowest point in January, which is <0.5 m deep in the middle of the lake. Our research investigated the effect of water level changes on the zooplankton community composition in Shahu Lake. Methods: We used both DNA metabarcoding method (MBC) (18S rRNA gene V4 region) and morphological method (MOI) to track the zooplankton community structure over four seasons in Shahu Lake (China). Results: Totals of 90 and 98 species of zooplankton were detected by MOI and MBC, respectively, with rotifers being the main zooplankton component. The α-diversity index of both methods increased from spring to summer and decreased from summer to autumn, reaching the lowest value in winter. NMDS and a cluster analysis showed that all zooplankton communities detected by MOI and MBC were significantly separated by season. The zooplankton community in winter was separated from that of the other three seasons, but the summer and autumn communities were more similar. Conclusions: Changes in the water level had significant effects on the zooplankton community composition. We found that MBC was more able to detect the differences in the zooplankton composition than MOI. MBC also had more advantages in copepod recognition. In our study, 37 species of copepods were detected by MBC, but only 11 species were detected by MOI. We concluded that MBC should be used to research the seasonal variations of zooplankton.

Neustonic microplastics and zooplankton in coastal waters of Cabrera Marine Protected Area (Western Mediterranean Sea)

The high abundance of microplastics in marine environments is becoming a growing threat for Marine Protected Areas (MPAs). Recent studies have identified microplastics (MP) as having multiple effects on biota, and it is of special interest to highlight their effects on zooplanktonic organisms. These organisms play a key role on local food web structures and there is a need to better understand the accumulation of MPs and associated contaminants within the food web. The present study addresses MP distribution and composition within Cabrera MPA as well as its effects on local zooplanktonic community composition. Neustonic microplastic and zooplankton abundance and composition were determined in the shallow coastal waters of Cabrera MPA during the summer season of 2019. Samples were taken using a sea surface manta-net, with a mesh size of 335 μm. MPs were found in all manta trawl samples for a total of 7047 MP with an overall mean abundance (± SD) of 3.52 (± 8.81) items/m³, a value higher than those reported for the majority of Western Mediterranean Sea basins. Qualitative and semi-quantitative analyses of the local zooplankton community suggest that its composition could be sensitive to MP abundances. Although no correlation was found between overall mean abundance of MP and zooplankton, a positive correlation was found between MP abundance and the abundance of the planktonic stage of the foraminifer Rosalina globularis (Tretomphalus). This species seems to be able to use MP for its dispersion. On the other hand, a negative correlation between MP abundance and Copepoda abundance was detected. This work confirms that the coastal zooplankton community composition within the MPA of Cabrera are not exempt from MP pollution and suggests important and complex interactions between MP and zooplankton organisms in coastal waters.

Metabarcoding analysis of trophic sources and linkages in the plankton community of the Kuroshio and neighboring waters

Trophic sources and pathways supporting early life stages are crucial for survival of forage fishes recruiting around the oligotrophic and unproductive Kuroshio. However, information is limited for the Kuroshio planktonic food web and its trophodynamics because of its high biodiversity. Here, we explore trophic sources and linkages in the Kuroshio plankton community using metabarcoding analysis of gut-content DNA for 22 mesozooplankton groups. The major prey was dinoflagellates and calanoids for omnivorous groups, and calanoids and gelatinous organisms for carnivorous groups. Larvaceans and hydrozoans were the most frequently appeared prey for both omnivores and carnivores, whereas they were minor constituents of the available prey in water samples. Although calanoids overlapped as major prey items for both omnivores and carnivores because they were the most available, contributions from phytoplankton and gelatinous prey differed among taxonomic groups. Further analysis of the metabarcoding data showed that in addition to omnivorous copepods like calanoids, gelatinous groups like larvaceans and hydrozoans were important hubs in the planktonic food web with their multiple trophic linkages to many components. These findings suggest that gelatinous organisms are important as supplementary prey and provide evidence of niche segregation on trophic sources among mesozooplankton groups in the Kuroshio.

Aquatic Organisms Research with DNA Barcodes

Since their inception, DNA barcodes have become a powerful tool for understanding the biodiversity and biology of aquatic species, with multiple applications in diverse fields such as food security, fisheries, environmental DNA, conservation, and exotic species detection. Nevertheless, most aquatic ecosystems, from marine to freshwater, are understudied, with many species disappearing due to environmental stress, mostly caused by human activities. Here we highlight the progress that has been made in studying aquatic organisms with DNA barcodes, and encourage its further development in assisting sustainable use of aquatic resources and conservation.

Dumpster diving for diatom plastid 16S rRNA genes

High throughput sequencing is improving the efficiency of monitoring diatoms, which inhabit and support aquatic ecosystems across the globe. In this study, we explored the potential of a standard V4 515F-806RB primer pair in recovering diatom plastid 16S rRNA sequences. We used PhytoREF to classify the 16S reads from our freshwater biofilm field sampling from three stream segments across two streams in south-eastern Australia and retrieved diatom community data from other, publicly deposited, Australian 16S amplicon datasets. When these diatom operational taxonomic units (OTUs) were traced using the default RDPII and NCBI databases, 68% were characterized as uncultured cyanobacteria. We analysed the 16S rRNA sequences from 72 stream biofilm samples, separated the chloroplast OTUs, and classified them using the PhytoREF database. After filtering the reads attributed to Bacillariophyta (relative abundance >1%), 71 diatom OTUs comprising more than 90% of the diatom reads in each stream biofilm sample were identified. Beta-diversity analyses demonstrated significantly different diatom assemblages and discrimination among river segments. To further test the approach, the diatom OTUs from our biofilm sampling were used as reference sequences to identify diatom reads from other Australian 16S rRNA datasets in the NCBI-SRA database. Across the three selected public datasets, 67 of our 71 diatom OTUs were detected in other Australian ecosystems. Our results show that diatom plastid 16S rRNA genes are readily amplified with existing 515F-806RB primer sets. Therefore, the volume of existing 16S rRNA amplicon datasets initially generated for microbial community profiling can also be used to detect, characterize, and map diatom distribution to inform phylogeny and ecological health assessments, and can be extended into a range of ecological and industrial applications. To our knowledge, this study represents the first attempt to classify freshwater samples using this approach and the first application of PhytoREF in Australia.

DNA metabarcoding reveals trophic niche diversity of micro and mesozooplankton species

Alternative pathways of energy transfer guarantee the functionality and productivity in marine food webs that experience strong seasonality. Nevertheless, the complexity of zooplankton interactions is rarely considered in trophic studies because of the lack of detailed information about feeding interactions in nature. In this study, we used DNA metabarcoding to highlight the diversity of trophic niches in a wide range of micro- and mesozooplankton, including ciliates, rotifers, cladocerans, copepods and their prey, by sequencing 16- and 18S rRNA genes. Our study demonstrates that the zooplankton trophic niche partitioning goes beyond both phylogeny and size and reinforces the importance of diversity in resource use for stabilizing food web efficiency by allowing for several different pathways of energy transfer. We further highlight that small, rarely studied zooplankton (rotifers and ciliates) fill an important role in the Baltic Sea pelagic primary production pathways and the potential of ciliates, rotifers and crustaceans in the utilization of filamentous and picocyanobacteria within the pelagic food web. The approach used in this study is a suitable entry point to ecosystem-wide food web modelling considering species-specific resource use of key consumers.

To denoise or to cluster, that is not the question: optimizing pipelines for COI metabarcoding and metaphylogeography

BACKGROUND

The recent blooming of metabarcoding applications to biodiversity studies comes with some relevant methodological debates. One such issue concerns the treatment of reads by denoising or by clustering methods, which have been wrongly presented as alternatives. It has also been suggested that denoised sequence variants should replace clusters as the basic unit of metabarcoding analyses, missing the fact that sequence clusters are a proxy for species-level entities, the basic unit in biodiversity studies. We argue here that methods developed and tested for ribosomal markers have been uncritically applied to highly variable markers such as cytochrome oxidase I (COI) without conceptual or operational (e.g., parameter setting) adjustment. COI has a naturally high intraspecies variability that should be assessed and reported, as it is a source of highly valuable information. We contend that denoising and clustering are not alternatives. Rather, they are complementary and both should be used together in COI metabarcoding pipelines.

RESULTS

Using a COI dataset from benthic marine communities, we compared two denoising procedures (based on the UNOISE3 and the DADA2 algorithms), set suitable parameters for denoising and clustering, and applied these steps in different orders. Our results indicated that the UNOISE3 algorithm preserved a higher intra-cluster variability. We introduce the program DnoisE to implement the UNOISE3 algorithm taking into account the natural variability (measured as entropy) of each codon position in protein-coding genes.  This correction increased the number of sequences retained by 88%. The order of the steps (denoising and clustering) had little influence on the final outcome.

CONCLUSIONS

We highlight the need for combining denoising and clustering, with adequate choice of stringency parameters, in COI metabarcoding. We present a program that uses the coding properties of this marker to improve the denoising step. We recommend researchers to report their results in terms of both denoised sequences (a proxy for haplotypes) and clusters formed (a proxy for species), and to avoid collapsing the sequences of the latter into a single representative. This will allow studies at the cluster (ideally equating species-level diversity) and at the intra-cluster level, and will ease additivity and comparability between studies.

Ciguatera Mini Review: 21st Century Environmental Challenges and the Interdisciplinary Research Efforts Rising to Meet Them

Globally, the livelihoods of over a billion people are affected by changes to marine ecosystems, both structurally and systematically. Resources and ecosystem services, provided by the marine environment, contribute nutrition, income, and health benefits for communities. One threat to these securities is ciguatera poisoning; worldwide, the most commonly reported non-bacterial seafood-related illness. Ciguatera is caused by the consumption of (primarily) finfish contaminated with ciguatoxins, potent neurotoxins produced by benthic single-cell microalgae. When consumed, ciguatoxins are biotransformed and can bioaccumulate throughout the food-web via complex pathways. Ciguatera-derived food insecurity is particularly extreme for small island-nations, where fear of intoxication can lead to fishing restrictions by region, species, or size. Exacerbating these complexities are anthropogenic or natural changes occurring in global marine habitats, e.g., climate change, greenhouse-gas induced physical oceanic changes, overfishing, invasive species, and even the international seafood trade. Here we provide an overview of the challenges and opportunities of the 21st century regarding the many facets of ciguatera, including the complex nature of this illness, the biological/environmental factors affecting the causative organisms, their toxins, vectors, detection methods, human-health oriented responses, and ultimately an outlook towards the future. Ciguatera research efforts face many social and environmental challenges this century. However, several future-oriented goals are within reach, including digital solutions for seafood supply chains, identifying novel compounds and methods with the potential for advanced diagnostics, treatments, and prediction capabilities. The advances described herein provide confidence that the tools are now available to answer many of the remaining questions surrounding ciguatera and therefore protection measures can become more accurate and routine.

Reconstructing marine plankton food web interactions using DNA metabarcoding

Knowledge of zooplankton in situ diet is critical for accurate assessment of marine ecosystem function and structure, but due to methodological constraints, there is still a limited understanding of ecological networks in marine ecosystems. Here, we used DNA-metabarcoding to study trophic interactions, with the aim to unveil the natural diet of zooplankton species under temporal variation of food resources. Several target consumers, including copepods and cladocerans, were investigated by sequencing 16S rRNA and 18S rRNA genes to identify prokaryote and eukaryote potential prey present in their guts. During the spring phytoplankton bloom, we found a dominance of diatom and dinoflagellate trophic links to copepods. During the summer period, zooplankton including cladocerans showed a more diverse diet dominated by cyanobacteria and heterotrophic prey. Our study suggests that copepods present trophic plasticity, changing their natural diet over seasons, and adapting their feeding strategies to the available prey spectrum, with some species being more selective. We did not find a large overlap of prey consumed by copepods and cladocerans, based on prey diversity found in their guts, suggesting that they occupy different roles in the trophic web. This study represents the first molecular approach to investigate several zooplankton-prey associations under seasonal variation, and highlights how, unlike other techniques, the diversity coverage is high when using DNA, allowing the possibility to detect a wide range of trophic interactions in plankton communities.