Metabarcoding and metabolome analyses of copepod grazing reveal feeding preference and linkage to metabolite classes in dynamic microbial plankton communities

Alpha-glucans from bacterial necromass indicate an intra-population loop within the marine carbon cycle

Phytoplankton blooms provoke bacterioplankton blooms, from which bacterial biomass (necromass) is released via increased zooplankton grazing and viral lysis. While bacterial consumption of algal biomass during blooms is well-studied, little is known about the concurrent recycling of these substantial amounts of bacterial necromass. We demonstrate that bacterial biomass, such as bacterial alpha-glucan storage polysaccharides, generated from the consumption of algal organic matter, is reused and thus itself a major bacterial carbon source in vitro and during a diatom-dominated bloom. We highlight conserved enzymes and binding proteins of dominant bloom-responder clades that are presumably involved in the recycling of bacterial alpha-glucan by members of the bacterial community. We furthermore demonstrate that the corresponding protein machineries can be specifically induced by extracted alpha-glucan-rich bacterial polysaccharide extracts. This recycling of bacterial necromass likely constitutes a large-scale intra-population energy conservation mechanism that keeps substantial amounts of carbon in a dedicated part of the microbial loop.

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.

Marine particle microbiomes during a spring diatom bloom contain active sulfate-reducing bacteria

Phytoplankton blooms fuel marine food webs with labile dissolved carbon and also lead to the formation of particulate organic matter composed of living and dead algal cells. These particles contribute to carbon sequestration and are sites of intense algal-bacterial interactions, providing diverse niches for microbes to thrive. We analyzed 16S and 18S ribosomal RNA gene amplicon sequences obtained from 51 time points and metaproteomes from 3 time points during a spring phytoplankton bloom in a shallow location (6-10 m depth) in the North Sea. Particulate fractions larger than 10 µm diameter were collected at near daily intervals between early March and late May in 2018. Network analysis identified two major modules representing bacteria co-occurring with diatoms and with dinoflagellates, respectively. The diatom network module included known sulfate-reducing Desulfobacterota as well as potentially sulfur-oxidizing Ectothiorhodospiraceae. Metaproteome analyses confirmed presence of key enzymes involved in dissimilatory sulfate reduction, a process known to occur in sinking particles at greater depths and in sediments. Our results indicate the presence of sufficiently anoxic niches in the particle fraction of an active phytoplankton bloom to sustain sulfate reduction, and an important role of benthic-pelagic coupling for microbiomes in shallow environments. Our findings may have implications for the understanding of algal-bacterial interactions and carbon export during blooms in shallow-water coastal areas.

Progress and challenges in exploring aquatic microbial communities using non-targeted metabolomics

Advances in bioanalytical technologies are constantly expanding our insights into complex ecosystems. Here, we highlight strategies and applications that make use of non-targeted metabolomics methods in aquatic chemical ecology research and discuss opportunities and remaining challenges of mass spectrometry-based methods to broaden our understanding of environmental systems.

Gut content metabarcoding of specialized feeders is not a replacement for environmental DNA assays of seawater in reef environments

In tropical marine ecosystems, the coral-based diet of benthic-feeding reef fishes provides a window into the composition and health of coral reefs. In this study, for the first time, we compare multi-assay metabarcoding sequences of environmental DNA (eDNA) isolated from seawater and partially digested gut items from an obligate corallivore butterflyfish (Chaetodon lunulatus) resident to coral reef sites in the South China Sea. We specifically tested the proportional and statistical overlap of the different approaches (seawater vs gut content metabarcoding) in characterizing eukaryotic community composition on coral reefs. Based on 18S and ITS2 sequence data, which differed in their taxonomic sensitivity, we found that gut content detections were only partially representative of the eukaryotic communities detected in the seawater based on low levels of taxonomic overlap (3 to 21%) and significant differences between the sampling approaches. Overall, our results indicate that dietary metabarcoding of specialized feeders can be complimentary to, but is no replacement for, more comprehensive environmental DNA assays of reef environments that might include the processing of different substrates (seawater, sediment, plankton) or traditional observational surveys. These molecular assays, in tandem, might be best suited to highly productive but cryptic oceanic environments (kelp forests, seagrass meadows) that contain an abundance of organisms that are often small, epiphytic, symbiotic, or cryptic.

Characterization of the Mycoremediation of n-Alkanes and Branched-Chain Alkanes by Filamentous Fungi from Oil-Polluted Soil Samples in Kazakhstan

For decades, researchers have focused on containing terrestrial oil pollution. The heterogeneity of soils, with immense microbial diversity, inspires them to transform pollutants and find cost-effective bioremediation methods. In this study, the mycoremediation potentials of five filamentous fungi isolated from polluted soils in Kazakhstan were investigated for their degradability of n-alkanes and branched-chain alkanes as sole carbon and energy sources. Dry weight estimation and gas chromatography-mass spectrometry (GC-MS) monitored the growth and the changes in the metabolic profile during degradation, respectively. Penicillium javanicum SBUG-M1741 and SBUG-M1742 oxidized medium-chain alkanes almost completely through mono- and di-terminal degradation. Pristane degradation by P. javanicum SBUG-M1741 was >95%, while its degradation with Purpureocillium lilacinum SBUG-M1751 was >90%. P. lilacinum SBUG-M1751 also exhibited the visible degradation potential of tetradecane and phytane, whereby in the transformation of phytane, both the mono- and di-terminal degradation pathways as well as α- and ß-oxidation steps could be described. Scedosporium boydii SBUG-M1749 used both mono- and di-terminal degradation pathways for n-alkanes, but with poor growth. Degradation of pristane by Fusarium oxysporum SBUG-M1747 followed the di-terminal oxidation mechanism, resulting in one dicarboxylic acid. These findings highlight the role of filamentous fungi in containing oil pollution and suggest possible degradation pathways.

Environmental changes associated with drying climate are expected to affect functional groups of pro- and microeukaryotes differently in temporary saline waters

Temporary ponds are among the most sensitive aquatic habitats to climate change. Their microbial communities have crucial roles in food webs and biogeochemical cycling, yet how their communities are assembled along environmental gradients is still understudied. This study aimed to reveal the environmental drivers of diversity (OTU-based richness, evenness, and phylogenetic diversity) and community composition from a network of saline temporary ponds, soda pans, in two consecutive spring seasons characterized by contrasting weather conditions. We used DNA-based molecular methods to investigate microbial community composition. We tested the effect of environmental variables on the diversity of prokaryotic (Bacteria, Cyanobacteria) and microeukaryotic functional groups (ciliates, heterotrophic flagellates and nanoflagellates, fungi, phytoplankton) within and across the years. Conductivity and the concentration of total suspended solids and phosphorus were the most important environmental variables affecting diversity patterns in all functional groups. Environmental conditions were harsher and they also had a stronger impact on community composition in the dry spring. Our results imply that these conditions, which are becoming more frequent with climate change, have a negative effect on microbial diversity in temporary saline ponds. This eventually might translate into community-level shifts across trophic groups with changing local conditions with implications for ecosystem functioning.

DNA extraction bias is more pronounced for microbial eukaryotes than for prokaryotes

DNA extraction and preservation bias is a recurring topic in DNA sequencing‐based microbial ecology. The different methodologies can lead to distinct outcomes, which has been demonstrated especially in studies investigating prokaryotic community composition. Eukaryotic microbes are ubiquitous, diverse, and increasingly a subject of investigation in addition to bacteria and archaea. However, little is known about how the choice of DNA preservation and extraction methodology impacts perceived eukaryotic community composition. In this study, we compared the effect of two DNA preservation methods and six DNA extraction methods on the community profiles of both eukaryotes and prokaryotes in phototrophic biofilms on seagrass (Zostera marina) leaves from the Baltic Sea. We found that, whereas both DNA preservation and extraction method caused significant bias in perceived community composition for both eukaryotes and prokaryotes, extraction bias was more pronounced for eukaryotes than for prokaryotes. In particular, soft‐bodied and hard‐shelled eukaryotes like nematodes and diatoms, respectively, were differentially abundant depending on the extraction method. We conclude that careful consideration of DNA preservation and extraction methodology is crucial to achieving representative community profiles of eukaryotes in marine biofilms and likely all other habitats containing diverse eukaryotic microbial communities.

DNA-metabarcoding reveals the importance of gelatinous zooplankton in the diet of Pandalus borealis, a keystone species in the Arctic

Information about the dietary composition of a species is crucial to understanding their position and role in the food web. Increasingly molecular approaches such as DNA-metabarcoding are used in studying trophic relations, not least because they may alleviate problems such as low taxonomic resolution or underestimation of digestible taxa in the diet. Here, we used DNA-metabarcoding with universal primers for cytochrome c oxidase I (COI), to study the diet composition of the Northern shrimp (Pandalus borealis), an Arctic keystone species with large socio-economic importance. Across locations, jellyfish and chaetognaths were the most important components in the diet of P. borealis, jointly accounting for 40-60% of the total read abundance. This dietary importance of gelatinous zooplankton contrasts sharply with published results based on SCA. At the same time, diet composition differed between fjord and shelf locations, pointing to different food webs supporting P. borealis in these two systems. Our study underscores the potential of molecular approaches to provide new insights into the diet of marine invertebrates that are difficult to obtain with traditional methods, and calls for a revision of the role of gelatinous zooplankton in the diet of the key Arctic species P. borealis, and in extension, Arctic food webs.

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.

Quantitative molecular detection of larval Atlantic herring (Clupea harengus) in stomach contents of Atlantic mackerel (Scomber scombrus) marks regions of predation pressure

Mortality rates in the early life-history stages of fishes are generally high yet identifying the causes remain unclear. Faltering recruitment rates of Atlantic herring (Clupea harengus) in the Norwegian Sea indicate a need to identify which mortality factors influence larval herring survival. Previous research suggests that increased predation pressure by Atlantic mackerel (Scomber scombrus) may contribute to the disconnect between spawning stock biomass and recruitment. To quantify the contribution of predation pressure by Atlantic mackerel to herring larval mortality, two research cruises were conducted within a probable “hot spot” (67–72° N) for intensified mackerel predation based on particle drift simulations. Mackerel stomach contents were analysed for herring larvae content using droplet digital polymerase chain reaction (ddPCR) with a quantitative molecular detection assay specific for herring. The ddPCR results demonstrate clear predation by mackerel on herring larvae and also suggest that the alternative use of visual examination may give misleading results. Our results show that mackerel should be considered a potentially important predator on herring larvae. The quantitative molecular assay presented here shows great promise as an efficient and specific tool to correctly identify and quantify predation pressure on early life-history stages of fishes.

Metabolomics-derived marker metabolites to characterize Phaeocystis pouchetii physiology in natural plankton communities

Phaeocystis pouchetii (Hariot) Lagerheim, 1893 regularly dominates phytoplankton blooms in higher latitudes spanning from the English Channel to the Arctic. Through zooplankton grazing and microbial activity, it is considered to be a key resource for the entire marine food web, but the actual relevance of biomass transfer to higher trophic levels is still under discussion. Cell physiology and algal nutritional state are suggested to be major factors controlling the observed variability in zooplankton grazing. However, no data have so far yielded insights into the metabolic state of Phaeocystis populations that would allow testing this hypothesis. Therefore, endometabolic markers of different growth phases were determined in laboratory batch cultures using comparative metabolomics and quantified in different phytoplankton blooms in the field. Metabolites, produced during exponential, early and late stationary growth of P. pouchetii, were profiled using gas chromatography-mass spectrometry. Then, metabolites were characterized that correlate with the growth phases using multivariate statistical analysis. Free amino acids characterized the exponential growth, whereas the early stationary phase was correlated with sugar alcohols, mono- and disaccharides. In the late stationary phase, free fatty acids, sterols and terpenes increased. These marker metabolites were then traced in Phaeocystis blooms during a cruise in the Barents Sea and North Norwegian fjords. About 50 endometabolites of P. pouchetii were detected in natural phytoplankton communities. Mannitol, scyllo-inositol, 24-methylcholesta-5,22-dien-3β-ol, and several free fatty acids were characteristic for Phaeocystis-dominated blooms but showed variability between them. Distinct metabolic profiles were detected in the nutrient-depleted community in the inner Porsangerfjord (< 0.5 µM NO₃^-, < 0.1 µM PO ₄ ^3- ), with high relative amounts of free mono- and disaccharides indicative for a limited culture. This study thereby shows how the variable physiology of phytoplankton can alter the metabolic landscape of entire plankton communities.

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.

High CO2 concentration and iron availability determine the metabolic inventory in an Emiliania huxleyi-dominated phytoplankton community

Ocean acidification (OA), a consequence of anthropogenic carbon dioxide (CO₂ ) emissions, strongly impacts marine ecosystems. OA also influences iron (Fe) solubility, affecting biogeochemical and ecological processes. We investigated the interactive effects of CO₂ and Fe availability on the metabolome response of a natural phytoplankton community. Using mesocosms we exposed phytoplankton to ambient (390 μatm) or future CO₂ levels predicted for the year 2100 (900 μatm), combined with ambient (4.5 nM) or high (12 nM) dissolved iron (dFe). By integrating over the whole phytoplankton community, we assigned functional changes based on altered metabolite concentrations. Our study revealed the complexity of phytoplankton metabolism. Metabolic profiles showed three stages in response to treatments and phytoplankton dynamics. Metabolome changes were related to the plankton group contributing respective metabolites, explaining bloom decline and community succession. CO₂ and Fe affected metabolic profiles. Most saccharides, fatty acids, amino acids and many sterols significantly correlated with the high dFe treatment at ambient pCO₂ . High CO₂ lowered the abundance of many metabolites irrespective of Fe. However, sugar alcohols accumulated, indicating potential stress. We demonstrate that not only altered species composition but also changes in the metabolic landscape affecting the plankton community may change as a consequence of future high-CO₂ oceans.

How do microbiota associated with an invasive seaweed vary across scales?

Communities are shaped by scale dependent processes. To study the diversity and variation of microbial communities across scales, the invasive and widespread seaweed Agarophyton vermiculophyllum presents a unique opportunity. We characterized pro- and eukaryotic communities associated with this holobiont across its known distribution range, which stretches over the northern hemisphere. Our data reveal that community composition and diversity in the holobiont vary at local but also larger geographic scales. While processes acting at the local scale (i.e., within population) are the main structuring drivers of associated microbial communities, changes in community composition also depend on processes acting at larger geographic scales. Interestingly, the largest analysed scale (i.e., native and non-native ranges) explained variation in the prevalence of predicted functional groups, which could suggest a functional shift in microbiota occurred over the course of the invasion process. While high variability in microbiota at the local scale supports A. vermiculophyllum to be a generalist host, we also identified a number of core taxa. These geographically independent holobiont members imply that cointroduction of specific microbiota may have additionally promoted the invasion process.

Long-Term Rewetting of Three Formerly Drained Peatlands Drives Congruent Compositional Changes in Pro- and Eukaryotic Soil Microbiomes through Environmental Filtering

Drained peatlands are significant sources of the greenhouse gas (GHG) carbon dioxide. Rewetting is a proven strategy used to protect carbon stocks; however, it can lead to increased emissions of the potent GHG methane. The response to rewetting of soil microbiomes as drivers of these processes is poorly understood, as are the biotic and abiotic factors that control community composition. We analyzed the pro- and eukaryotic microbiomes of three contrasting pairs of minerotrophic fens subject to decade-long drainage and subsequent long-term rewetting. Abiotic soil properties including moisture, dissolved organic matter, methane fluxes, and ecosystem respiration rates were also determined. The composition of the microbiomes was fen-type-specific, but all rewetted sites showed higher abundances of anaerobic taxa compared to drained sites. Based on multi-variate statistics and network analyses, we identified soil moisture as a major driver of community composition. Furthermore, salinity drove the separation between coastal and freshwater fen communities. Methanogens were more than 10-fold more abundant in rewetted than in drained sites, while their abundance was lowest in the coastal fen, likely due to competition with sulfate reducers. The microbiome compositions were reflected in methane fluxes from the sites. Our results shed light on the factors that structure fen microbiomes via environmental filtering.

The potential of sedimentary ancient DNA for reconstructing past sea ice evolution

Sea ice is a crucial component of the Arctic climate system, yet the tools to document the evolution of sea ice conditions on historical and geological time scales are few and have limitations. Such records are essential for documenting and understanding the natural variations in Arctic sea ice extent. Here we explore sedimentary ancient DNA (aDNA), as a novel tool that unlocks and exploits the genetic (eukaryote) biodiversity preserved in marine sediments specifically for past sea ice reconstructions. Although use of sedimentary aDNA in paleoceanographic and paleoclimatic studies is still in its infancy, we use here metabarcoding and single-species quantitative DNA detection methods to document the sea ice conditions in a Greenland Sea marine sediment core. Metabarcoding has allowed identifying biodiversity changes in the geological record back to almost ~100,000 years ago that were related to changing sea ice conditions. Detailed bioinformatic analyses on the metabarcoding data revealed several sea-ice-associated taxa, most of which previously unknown from the fossil record. Finally, we quantitatively traced one known sea ice dinoflagellate in the sediment core. We show that aDNA can be recovered from deep-ocean sediments with generally oxic bottom waters and that past sea ice conditions can be documented beyond instrumental time scales. Our results corroborate sea ice reconstructions made by traditional tools, and thus demonstrate the potential of sedimentary aDNA, focusing primarily on microbial eukaryotes, as a new tool to better understand sea ice evolution in the climate system.

Metabarcoding Insights Into the Trophic Behavior and Identity of Intertidal Benthic Foraminifera

Foraminifera are ubiquitous marine protists with an important role in the benthic carbon cycle. However, morphological observations often fail to resolve their exact taxonomic placement and there is a lack of field studies on their particular trophic preferences. Here, we propose the application of metabarcoding as a tool for the elucidation of the in situ feeding behavior of benthic foraminifera, while also allowing the correct taxonomic assignment of the feeder, using the V9 region of the 18S (small subunit; SSU) rRNA gene. Living foraminiferal specimens were collected from two intertidal mudflats of the Wadden Sea and DNA was extracted from foraminiferal individuals and from the surrounding sediments. Molecular analysis allowed us to confirm that our foraminiferal specimens belong to three genetic types: Ammonia sp. T6, Elphidium sp. S5 and Haynesina sp. S16. Foraminiferal intracellular eukaryote communities reflected to an extent those of the surrounding sediments but at different relative abundances. Unlike sediment eukaryote communities, which were largely determined by the sampling site, foraminiferal intracellular eukaryote communities were driven by foraminiferal species, followed by sediment depth. Our data suggests that Ammonia sp. T6 can predate on metazoan classes, whereas Elphidium sp. S5 and Haynesina sp. S16 are more likely to ingest diatoms. These observations, alongside the use of metabarcoding in similar ecological studies, significantly contribute to our overall understanding of the ecological roles of these protists in intertidal benthic environments and their position and function in the benthic food webs.

The microbiota of technosols resembles that of a nearby forest soil three years after their establishment

Technosols can be used to rehabilitate degraded land and reuse wastes. Ideally, these newly formed soils should also fulfil the main soil functions. In this study, initially, we characterized the physicochemical and microbial properties of different formulations and their ingredients (i.e., dirt from a waste recovery plant, recycled bentonite, sewage sludge). When these technosols were then used for the rehabilitation of a quarry, the evolution of such properties was monitored for three consecutive years. Physicochemical and microbial properties were compared to those of a reference soil from a nearby forest. Diversity and composition of prokaryotes and eukaryotes were determined using 16S and 18S rRNA amplicon sequencing. Three years after establishment, as much as 78.8% and 63.9% of the prokaryotic and eukaryotic orders, respectively, were shared between the technosols and the reference forest soil. Although technosols initially showed lower values of CO₂ emission, compaction and functional diversity (Biolog EcoPlates™), at the end of the study these values were similar to those observed in the reference forest soil. It was concluded that the microbiota of the studied technosols resembles that of the nearby forest soil after just three years of establishment.

Application of sewage sludge to agricultural soil increases the abundance of antibiotic resistance genes without altering the composition of prokaryotic communities

The application of sewage sludge as soil amendment is a common agricultural practice. However, wastewater treatment plants, sewage sludge and sewage sludge-amended soils have been reported as hotspots for the appearance and dissemination of antibiotic resistance, driven, among other factors, by selection pressure exerted by co-exposure to antibiotics and heavy metals. To address this threat to environmental and human health, soil samples from a long-term (24 years) field experiment, carried out to study the impact of thermally dried and anaerobically digested sewage sludge (at different doses and frequencies of application) on agricultural soil quality, were investigated for the presence of genes encoding antibiotic resistance (ARGs) and mobile genetic elements (MGEs). Sewage sludge-induced changes in specific soil physicochemical and microbial properties, as indicators of soil quality, were also investigated. The application of sewage sludge increased the total concentration of copper and zinc in amended soils, but without affecting the bioavailability of these metals, possibly due to the high values of soil pH and organic matter content. Soil microbal quality, as reflected by the value of the Soil Quality Index, was higher in sewage sludge-amended soils. Similarly, the application of sewage sludge increased soil microbial activity and biomass, as well as the abundance of ARGs and MGE genes, posing a risk of dissemination of antibiotic resistance. In contrast, the composition of soil prokaryotic communities was not significantly altered by the application of sewage sludge. We found correlation between soil Cu and Zn concentrations and the abundance of ARGs and MGE genes. It was concluded that sewage sludge-derived amendments must be properly treated and managed if they are to be applied to agricultural soil.

DNA metabarcoding unveils multiscale trophic variation in a widespread coastal opportunist

A thorough understanding of ecological networks relies on comprehensive information on trophic relationships among species. Since unpicking the diet of many organisms is unattainable using traditional morphology-based approaches, the application of high-throughput sequencing methods represents a rapid and powerful way forward. Here, we assessed the application of DNA metabarcoding with nearly universal primers for the mitochondrial marker cytochrome c oxidase I in defining the trophic ecology of adult brown shrimp, Crangon crangon, in six European estuaries. The exact trophic role of this abundant and widespread coastal benthic species is somewhat controversial, while information on geographical variation remains scant. Results revealed a highly opportunistic behaviour. Shrimp stomach contents contained hundreds of taxa (>1,000 molecular operational taxonomic units), of which 291 were identified as distinct species, belonging to 35 phyla. Only twenty ascertained species had a mean relative abundance of more than 0.5%. Predominant species included other abundant coastal and estuarine taxa, including the shore crab Carcinus maenas and the amphipod Corophium volutator. Jacobs' selectivity index estimates based on DNA extracted from both shrimp stomachs and sediment samples were used to assess the shrimp's trophic niche indicating a generalist diet, dominated by crustaceans, polychaetes and fish. Spatial variation in diet composition, at regional and local scales, confirmed the highly flexible nature of this trophic opportunist. Furthermore, the detection of a prevalent, possibly endoparasitic fungus (Purpureocillium lilacinum) in the shrimp's stomach demonstrates the wide range of questions that can be addressed using metabarcoding, towards a more robust reconstruction of ecological networks.

Metabarcoding of shrimp stomach content: Harnessing a natural sampler for fish biodiversity monitoring

Given their positioning and biological productivity, estuaries have long represented key providers of ecosystem services and consequently remain under remarkable pressure from numerous forms of anthropogenic impact. The monitoring of fish communities in space and time is one of the most widespread and established approaches to assess the ecological status of estuaries and other coastal habitats, but traditional fish surveys are invasive, costly, labour intensive and highly selective. Recently, the application of metabarcoding techniques, on either sediment or aqueous environmental DNA, has rapidly gained popularity. Here, we evaluate the application of a novel, high‐throughput DNA‐based monitoring tool to assess fish diversity, based on the analysis of the gut contents of a generalist predator/scavenger, the European brown shrimp, Crangon crangon. Sediment and shrimp samples were collected from eight European estuaries, and DNA metabarcoding (using both 12S and COI markers) was carried out to infer fish assemblage composition. We detected 32 teleost species (16 and 20, for 12S and COI, respectively). Twice as many species were recovered using metabarcoding than by traditional net surveys. By comparing and interweaving trophic, environmental DNA and traditional survey‐based techniques, we show that the DNA‐assisted gut content analysis of a ubiquitous, easily accessible, generalist species may serve as a powerful, rapid and cost‐effective tool for large‐scale, routine estuarine biodiversity monitoring.

16S rRNA gene metabarcoding and TEM reveals different ecological strategies within the genus Neogloboquadrina (planktonic foraminifer)

Uncovering the complexities of trophic and metabolic interactions among microorganisms is essential for the understanding of marine biogeochemical cycling and modelling climate-driven ecosystem shifts. High-throughput DNA sequencing methods provide valuable tools for examining these complex interactions, although this remains challenging, as many microorganisms are difficult to isolate, identify and culture. We use two species of planktonic foraminifera from the climatically susceptible, palaeoceanographically important genus Neogloboquadrina, as ideal test microorganisms for the application of 16S rRNA gene metabarcoding. Neogloboquadrina dutertrei and Neogloboquadrina incompta were collected from the California Current and subjected to either 16S rRNA gene metabarcoding, fluorescence microscopy, or transmission electron microscopy (TEM) to investigate their species-specific trophic interactions and potential symbiotic associations. 53–99% of 16S rRNA gene sequences recovered from two specimens of N. dutertrei were assigned to a single operational taxonomic unit (OTU) from a chloroplast of the phylum Stramenopile. TEM observations confirmed the presence of numerous intact coccoid algae within the host cell, consistent with algal symbionts. Based on sequence data and observed ultrastructure, we taxonomically assign the putative algal symbionts to Pelagophyceae and not Chrysophyceae, as previously reported in this species. In addition, our data shows that N. dutertrei feeds on protists within particulate organic matter (POM), but not on bacteria as a major food source. In total contrast, of OTUs recovered from three N. incompta specimens, 83–95% were assigned to bacterial classes Alteromonadales and Vibrionales of the order Gammaproteobacteria. TEM demonstrates that these bacteria are a food source, not putative symbionts. Contrary to the current view that non-spinose foraminifera are predominantly herbivorous, neither N. dutertrei nor N. incompta contained significant numbers of phytoplankton OTUs. We present an alternative view of their trophic interactions and discuss these results within the context of modelling global planktonic foraminiferal abundances in response to high-latitude climate change.

Eelgrass Leaf Surface Microbiomes Are Locally Variable and Highly Correlated with Epibiotic Eukaryotes

Eelgrass (Zostera marina) is a marine foundation species essential for coastal ecosystem services around the northern hemisphere. Like all macroscopic organisms, it possesses a microbiome (here defined as an associated prokaryotic community) which may play critical roles in modulating the interaction of eelgrass with its environment. For example, its leaf surface microbiome could inhibit or attract eukaryotic epibionts which may overgrow the eelgrass leading to reduced primary productivity and subsequent eelgrass meadow decline. We used amplicon sequencing of the 16S and 18S rRNA genes of prokaryotes and eukaryotes to assess the leaf surface microbiome (prokaryotes) as well as eukaryotic epibionts in- and outside lagoons on the German Baltic Sea coast. Prokaryote microbiomes varied substantially both between sites inside lagoons and between open coastal and lagoon sites. Water depth, leaf area and biofilm chlorophyll a concentration explained a large amount of variation in both prokaryotic and eukaryotic community composition. The prokaryotic microbiome and eukaryotic epibiont communities were highly correlated, and network analysis revealed disproportionate co-occurrence between a limited number of eukaryotic taxa and several bacterial taxa. This suggests that eelgrass leaf surfaces are home to a mosaic of microbiomes of several epibiotic eukaryotes, in addition to the microbiome of the eelgrass itself. Our findings thereby underline that eukaryotic diversity should be taken into account in order to explain prokaryotic microbiome assembly and dynamics in aquatic environments.