Search

Over the past decades the sampling of environmental DNA (eDNA) — encompassing the DNA of all organisms present in an environmental sample1 — has emerged as a technique for biodiversity monitoring and discovery in a diversity of environments. Avoiding the physical collection and identification of biota, this approach is praised for its independence of taxonomic expertise and has changed the way biologists study biodiversity. However, a common result in eDNA studies is the finding of unexpected taxa which are often removed by conservative bioinformatic filters or disregarded, since the authors are uncertain about the result and rarely have the interest, time, skills, and/or resources to return to the field and confirm with actual specimens2. Here, we report a case in which an eDNA discovery led to the physical localization of a member of the Micrognathozoa (Figure 1B) — a rare group of limnic micrometazoans, and the animal phylum to be discovered last3, which is the sister group to rotifers4,5. To this day, Micrognathozoa still comprises only a single named species from Greenland and a few additional disparate places

Although there is high interest in realizing the great potential of applying eDNA methods to marine monitoring in the Northeast Atlantic, the implementation process has proven challenging. Researchers and stakeholders lack a common understanding about how eDNA methods can be used as well as a basis for discussing the various challenges and ambitions of both parties. The aim of the UNIFIeD project was to assess the current status of implementation; identifying these challenges and ambitions and provide a basis for enhanced communication between researchers and stakeholders. This report describes results from surveys addressing the key issues of eDNA application, standardization and implementation. In addition, the discussions and consensus from the following workshop is presented.

The Southern Ocean, in particular the southwest Atlantic sector, is experiencing rapid environmental changes. A long-term trend of density changes of key pelagic species has been noted over the last decades: Antarctic krill populations are declining whilst salps are on the rise and shifting their distribution poleward. A similar poleward expansion is anticipated for a third key player, the hyperiid amphipod crustacean Themisto gaudichaudii, leading to an increasing overlap of the distributions of these three species. Due to major knowledge gaps in the ecology, and genetic connectivity of T. gaudichaudii, the likelihood of this shift and its consequences for the pelagic food web structure remain largely unexplored. In this context, Themisto’s genetic and trophic connectivity as well as thermal response were investigated with state-of-the-art molecular methods. Phylogeographic analyses showed genetic homogeneity between localities in the Southern Ocean and Atlantic waters combined with high degree of phenotypic plasticity enabling different lineages to thrive in regions further south. Diet analyses using DNA metabarcoding were applied to characterize regional variation in diet. These analyses showed a diet predominantly composed of krill, in particular in the Antarctic Peninsula region, showing that Themisto’s poleward range expansion can further impact the already declining krill stocks. It also unexpectedly revealed ctenophores to be an important prey, despite their reputation as “trophic dead-end”. Transcriptome analyses were used to study the thermal response of Themisto individuals from different geographic populations that were experimentally exposed to heat and cold treatments. The analysis of differentially expressed genes showed that genetic lineages differ in thermal tolerances. It also revealed a wide range of molecular mechanisms in Themisto amphipods to cope with thermal stress. These findings contribute to better predict the impact of climate-driven rang

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 relationships, 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 stomach content analysis. 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 underlines 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.

Ecosystem biomonitoring is crucial for proper management of natural communities during the Anthropocene era. With the advent of new sequencing technologies, DNA metabarcoding has been proposed as a game-changing tool for biomonitoring. In this Thesis we plead for the use of metabarcoding of a highly variable marker to infer not only the interspecies but also the intraspecies variability to assess both biogeographic, at the species level, and metaphylogeographic patterns, at the haplotype level. We focused on highly complex hard-substratum benthic littoral communities. The term "Metaphylogeography", coined in this Thesis, refers to the study of phylogeographic patterns of many species at the same time using metabarcoding data. However, as of the start of this Thesis, only a few studies had tested the metabarcoding method to directly characterize the whole eukaryotic community in highly diverse benthic ecosystems. This required to set up and calibrate methods for these communities as a prior step. We first evaluated both the sampling methods and the bioinformatic pipelines. We assessed the viability of detecting the environmental DNA released from the benthic community into the adjacent water layer using metabarcoding of COI with highly degenerated primers targeting the whole eukaryotic community. We sampled water from 0 to 20m from shallow rocky benthic communities and compared the DNA signal with the results obtained from metabarcoding directly the benthic communities by traditional quadrat sampling. We also designed a pipeline combining clustering and denoising methods to treat metabarcoding data of COI. We considered the entropy of each codon position of this coding fragment both to improve the detection of spurious sequences and to calibrate the best performing parameters of the software used. In addition, we created our own denoising program, DnoisE, to incorporate information on the codon position. This new code and parameter calibration were required as the co

DNA metabarcoding is broadly used in biodiversity studies encompassing a wide range of organisms. Erroneous amplicons, generated during amplification and sequencing procedures, constitute one of the major sources of concern for the interpretation of metabarcoding results. Several denoising programs have been implemented to detect and eliminate these errors. However, almost all denoising software currently available has been designed to process non-coding ribosomal sequences, most notably prokaryotic 16S rDNA. The growing number of metabarcoding studies using coding markers such as COI or RuBisCO demands a re-assessment and calibration of denoising algorithms. Here we present DnoisE, the first denoising program designed to detect erroneous reads and merge them with the correct ones using information from the natural variability (entropy) associated to each codon position in coding barcodes. We have developed an open-source software using a modified version of the UNOISE algorithm. DnoisE implements different merging procedures as options, and can incorporate codon entropy information either retrieved from the data or supplied by the user. In addition, the algorithm of DnoisE is parallelizable, greatly reducing runtimes on computer clusters. Our program also allows different input file formats, so it can be readily incorporated into existing metabarcoding pipelines.

Traditional morphological methods for species identification are highly time consuming, especially for small organisms, such as Foraminifera, a group of shell-building microbial eukaryotes. To analyze large amounts of samples more efficiently, species identification methods have extended to molecular tools in the last few decades. Although a wide range of phyla have good markers available, for Foraminifera only one hypervariable marker from the ribosomal region (18S) is widely used. Recently a new mitochondrial marker cytochrome oxidase subunit 1 (COI) has been sequenced. Here we investigate whether this marker has a higher potential for species identification compared to the ribosomal marker. We explore the genetic variability of both the 18S and COI markers in 22 benthic foraminiferal morphospecies (orders Miliolida and Rotaliida). Using single-cell DNA, the genetic variability within specimens (intra) and between specimens (inter) of each species was assessed using next-generation sequencing. Amplification success rate was twice as high for COI (151/200 specimens) than for 18S (73/200 specimens). The COI marker showed greatly decreased intra- and inter-specimen variability compared to 18S in six out of seven selected species. The 18S phylogenetic reconstruction fails to adequately cluster multiple species together in contrast to COI. Additionally, the COI marker helped recognize misclassified specimens difficult to morphologically identify to the species level. Integrative taxonomy, combining morphological and molecular characteristics, provides a robust picture of the foraminiferal species diversity. Finally, we suggest the use of a set of sequences (two or more) to describe species showing intra-genomic variability additionally to using multiple markers. Our findings highlight the potential of the newly discovered mitochondrial marker for molecular species identification and metabarcoding purposes.

In a recent paper, “Environmental DNA: What's behind the term? Clarifying the terminology and recommendations for its future use in biomonitoring,” Pawlowski et al. argue that the term eDNA should be used to refer to the pool of DNA isolated from environmental samples, as opposed to only extra-organismal DNA from macro-organisms. We agree with this view. However, we are concerned that their proposed two-level terminology specifying sampling environment and targeted taxa is overly simplistic and might hinder rather than improve clear communication about environmental DNA and its use in biomonitoring. This terminology is based on categories that are often difficult to assign and uninformative, and it overlooks a fundamental distinction within eDNA: the type of DNA (organismal or extra-organismal) from which ecological interpretations are derived.

Marine biodiversity underpins ecosystem health and societal well-being. Preservation of biodiversity hotspots is a global challenge. Molecular tools, like DNA barcoding and metabarcoding, hold great potential for biodiversity monitoring, possibly outperforming more traditional taxonomic methods. However, metabarcoding-based biodiversity assessments are limited by the availability of sequences in barcoding reference databases; a lack thereof results in high percentages of unassigned sequences. In this study we (i) present the current status of known vs. barcoded marine species at a global scale based on online taxonomic and genetic databases; and (ii) compare the current status with data from ten years ago. Then we analyzed occurrence data of marine animal species from five Large Marine Ecosystems (LMEs) classified as biodiversity hotspots, to identify any consistent disparities in COI barcoding coverage between geographic regions and at phylum level. Barcoding coverage varied among LMEs (from 36.8% to 62.4% COI-barcoded species) and phyla (from 4.8% to 74.7% COI-barcoded species), with Porifera, Bryozoa and Platyhelminthes being highly underrepresented, compared to Chordata, Arthropoda and Mollusca. We demonstrate that although barcoded marine species increased from 9.5% to 14.2% since the last assessment in 2011, about 15,000 (corresponding to 7.8% increase) new species were described from 2011 to 2021. The next ten years will thus be crucial to enroll concrete collaborative measures and long term initiatives (e.g., Horizon 2030, Ocean Decade) to populate barcoding libraries for the marine realm.

Global environmental changes may have a profound impact on ecosystems. In this context, it is crucial to gather biological and ecological information of the main species in marine communities to predict and mitigate potential effects of shifts in their distribution, abundance, and interactions. Using genotyping by sequencing (GBS), we assessed the genetic structure of a keystone species in the Mediterranean shallow littoral ecosystems, the black sea urchin Arbacia lixula. This bioengineer species can shape their communities due to its grazing activity and it is experiencing an ongoing expansion with increasing temperatures. The population genomic analyses on 5,241 loci sequenced in 240 individuals from 11 Mediterranean sampled populations revealed that all populations were diverse and showed significant departure from equilibrium. Albeit genetic differentiation was in general shallow, a significant break separated the western and eastern Mediterranean populations, a break not detected in previous studies with less resolutive markers. Notably, no clear effect of the Almería-Oran front, an important break in the Atlanto-Mediterranean transition, could be detected among the western basin populations, where only a slight differentiation of the two northernmost populations was found. Despite the generally low levels of genetic differentiation found, we identified candidate regions for local adaptation by combining different genomic analysis with environmental data. Salinity, rather than temperature, seemed to be an important driver of genetic structure in A. lixula. Overall, from a population genomics standpoint, there is ample scope for A. lixula to continue thriving and adapting in the warming Mediterranean.

In a recent paper, “Environmental DNA: What's behind the term? Clarifying the terminology and recommendations for its future use in biomonitoring,” Pawlowski et al. argue that the term eDNA should be used to refer to the pool of DNA isolated from environmental samples, as opposed to only extra-organismal DNA from macro-organisms. We agree with this view. However, we are concerned that their proposed two-level terminology specifying sampling environment and targeted taxa is overly simplistic and might hinder rather than improve clear communication about environmental DNA and its use in biomonitoring. This terminology is based on categories that are often difficult to assign and uninformative, and it overlooks a fundamental distinction within eDNA: the type of DNA (organismal or extra-organismal) from which ecological interpretations are derived.

Reproductive rates of copepods are temperature-dependent, but poorly known for small copepods at low temperatures, hindering the predictions of population dynamics and secondary production in high-latitude ecosystems. We investigated egg hatching rates, hatching success and egg production of the small copepods Oithona similis and Microsetella norvegica (sac spawners) and Microcalanus pusillus (broadcast spawner) between March and August. Incubations were performed at ecologically relevant temperatures between 1.3 and 13.2°C, and egg production rates were calculated. All egg hatching rates were positively correlated to temperature, although with large species-specific differences. At the lowest temperatures, M. pusillus eggs hatched within 4 days, whereas the eggs from sac spawners took 3–8 weeks to hatch. The egg hatching success was ≤25% for M. pusillus, >75% for O. similis and variable for M. norvegica. The maximum weight-specific egg production rate (μg C μg−1 C d−1) of M. pusillus was higher (0.22) than O. similis (0.12) and M. norvegica (0.06). M. norvegica reproduction peaked at 6–8°C, the prevailing in situ temperatures during its reproductive period. The difference in reproductive rates indicates species-specific thermal plasticity for the three copepods, which could have implications for present and future population dynamics of the species in arctic fjords.