Brief Research Report ARTICLE
Early Alert of Biological Risk in a Coastal Lagoon Through eDNA Metabarcoding
- 1Department of Functional Biology, University of Oviedo, Oviedo, Spain
- 2USR 3278-CRIOBE-CNRS-EPHE, Laboratoire d’excellence “CORAIL,” Universite de Perpignan-CBETM, Perpignan, France
- 3Laboratoire d’Excellence “CORAIL,” Centre de Recherche Insulaire et Observatoire de l’Environnement (CRIOBE), Papetoai, French Polynesia
Mediterranean coastal lagoons are environmentally important ecosystems whose conservation has been challenged due to anthropogenic impacts that promoted the expansion of non-indigenous and, sometimes, invasive species. Therefore, it is crucial to inventory biodiversity in these areas for the development of strategies of conservation and management. Classical methods used for biodiversity surveys and detection of non-native species may be unsuccessful for the detection and identification of species in early development stages such as cryptic, microscopic, elusive, and new coming species at low population density. The development of metabarcoding techniques in the last decade offers new opportunities for reliable biodiversity surveillance and facilitates early detection of nuisance species. The objective of this study was to analyze the species occurring in the protected coastal lagoon Canet-Saint Nazaire using a simple sampling protocol based on water samples and environmental DNA (eDNA) metabarcoding with a single barcode (cytochrome c oxidase subunit I [COI] gene). Two invasive species (Polydora cornuta and Acartia tonsa), two polychaete bioindicators of pollution (Hediste diversicolor and Capitella capitata), and one alga that produces harmful algal blooms were detected from only 6 L of water, indicating environmental degradation in the lagoon despite its protected status. These results demonstrate the importance of COI as single barcode together with eDNA as an ecological early warning system and suggest the need for environmental restoration in this lagoon.
Coastal lagoons occupy 13% of the coastal area worldwide (Barnes, 1980) and are among the marine habitats showing the highest biological productivity, by providing diverse habitat types for many species, nursery areas, and feeding grounds for marine and estuarine fishes (Pérez-Ruzafa et al., 2011). They are distinctive ecosystems because they are shallow coastal water bodies separated from the ocean by a barrier and connected intermittently to the sea (Kjerfve, 1994). They also support important fisheries and allow for intensive aquaculture exploitation (Cataudella et al., 2015). Despite their environmental and economic importance, and protected status, lagoons suffer from several threats derived from human activities such as the effects of climate change, pollution, eutrophication, and the introduction of non-indigenous species (NIS) (Reizopoulou et al., 1996; Chapman, 2012). Human activities, such as the increases in maritime traffic and the opening of the Suez Canal, have facilitated NIS introductions in the Mediterranean Sea (Katsanevakis et al., 2014). In coastal lagoons, NIS settlement is facilitated by their naturally stressful conditions, pollution, loss of native species, and intense shipping in the numerous nearby harbors (Ruiz et al., 2000; Crooks et al., 2011). NIS can become invasive and affect both native species and economic activities in the area (Galil, 2007).
Together with NIS, eutrophication plays an important role in the degradation of these estuarine systems. Inefficient wastewater management and intense agricultural activities in the catchment area of the lagoons can increase the concentration of nutrients (Carlier et al., 2008). The consequent eutrophic state facilitates the uncontrolled growth of organisms such as dinoflagellates, diatoms, and cyanobacteria that produce harmful algal blooms (HABs, also referred as red tides) (Anderson et al., 2002). The increasing algal bloom episodes is not only favored by eutrophication but also by climate change and the introduction of new strains by ballast water (Hallegraeff, 1993; Moore et al., 2008). HABs have several effects on the environment. The mere growth and following decay of organic material cause anoxic conditions, leading to the death of aquatic life. Some HABs also produce toxins that affect both marine life and human health even at low densities (Sellner et al., 2003). In humans, these toxins can cause different types of poisoning (diarrheal, neurotoxic, and paralytic). The most frequent are diarrheic shellfish poisoning (DSP), mainly due to toxic strains of Prorocentrum spp. and Dinophysis spp. (Yasumoto et al., 1980; Bravo et al., 2001), and paralytic shellfish poisoning (PSP) predominantly linked to toxic strains of Alexandrium spp. (Anderson et al., 2012).
The correct management and conservation of coastal areas requires efficient assessment of their biodiversity and detection and identification of species that may be of environmental concern, such as pollution indicators and NIS. For this purpose, species surveys are carried out periodically, generally consisting of manual sampling and visual identification of sampled specimens by taxonomical experts. Such a classical morphological analysis requires high sampling effort and might be inefficient in the detection of some species, for example those that are in early development stages, cryptic species, or microorganisms (Ficetola et al., 2008). For NIS, their early detection is crucial to prevent their establishment and dispersion because eradication or control are more efficient when the species are at low density soon after introduction (Gozlan et al., 2010). However, low population densities require greater survey efforts for detection, and many newly established populations may go unnoticed (Blanchet, 2012). Classical approaches for phytoplankton identification are often considered time consuming while requiring expertise in taxonomical identification. Identifying dinoflagellates represent a greater challenge due to high morphological similarities and a lack of unique characters between different species (Lin et al., 2009). DNA metabarcoding is a rapidly evolving method for assessing biodiversity that exceeds the limitations of traditional methods. It combines the use of environmental DNA (eDNA) and next-generation sequencing (NGS) allowing the detection of species from single cells in an environmental sample (e.g., soil or water) (Taberlet et al., 2012). Although eDNA metabarcoding was found to be less sensitive than a targeted monitoring approach for generating detailed specific distribution data (Furlan et al., 2015; Ardura et al., 2016; Bylemand et al., 2019), it has the ability to provide baseline information on biodiversity patterns (e.g., cryptic species and novel incursions of invasive species) (Blackman et al., 2017), capturing a great part of the aquatic community diversity (Zaiko et al., 2015). Borrell et al. (2017) have proposed the use of metabarcoding as an early alert method for the detection of invasive species in ports. They used two metabarcodes (cytochrome oxidase I and 18S rRNA genes) on 3 L of water from each sampling point. For the need of easy methods for alert of environmental disturbances, here, we will test a similar but simplified approach for the exploration of nuisance species, including invasive species and HABs.
We apply metabarcoding (eDNA) for the detection of aquatic nuisance organisms in the coastal lagoon Canet-Saint Nazaire, situated in the French Mediterranean coast. Previous studies developed in this lagoon using a rapid assessment of invertebrate species and barcoding (Ardura and Planes, 2017) suggested a high level of degradation and vulnerability of the lagoon. However, in that study, small species like HABs and organisms in early development stages could not be recognized. In the present investigation, eDNA was extracted from water samples from the canal that connects the lagoon to the open sea, and a region of the cytochrome c oxidase subunit I (COI) gene was used as a single metabarcode for species identification, since it has been used as a barcode of high-resolution power (Hebert et al., 2003a) and reference databases for aquatic organisms are more complete for this gene than other barcodes (Weigand et al., 2019). The principal aim of this study was focused on the detection of species that may reveal a risk for the conservation status and environmental health of the lagoon.
Materials and Methods
Sampling was conducted in the canal that connects the Canet-Saint Nazaire Lagoon to the open sea (Figure 1). This lagoon is a special protection area (SPA) for birds within the Natura 2000 network, and it is one of the many coastal lagoons spanning the French Mediterranean coastline. It is a semiclosed system, with its principal water sources coming from two rivers and streams and, to a lesser extent, a canal connected to the sea. This, along with seasonal variations in temperature, drought periods, and sudden and intense rainfalls, causes important changes in its depth and salinity that varies from 13.2 to 35.6 during the year (Vouvé et al., 2014). Besides this, the lagoon suffers from intense pollution, eutrophication, and sediment filling due to human activities in its catchment area (agriculture, water treatment, and tourism) and scarce water exchange with the sea (Carlier et al., 2008; Souchu et al., 2010). Within 4 km of the lagoon, Canet port is located, which is a popular location for recreational boating with space for 977 vessels up to 24 m in length (Portbooker website, 2018). The presence of this harbor, along with the lagoon’s brackish water, makes Canet-Saint Nazaire extremely vulnerable to introduction of non-native species (Paavola et al., 2005; Ardura and Planes, 2017).
Figure 1. Lagoon of Canet-Saint Nazaire (below) and port of Canet (above) indicating the sampling site (Google Maps, 2018).
Water samples of the lagoon were taken in the proximity of the canal (42°39′25.19″N 3° 1′39.30″E). Six liters of surface water was collected in three sterile bottles of 2 L (three replicates) on November–December 2016. All the water from each bottle was vacuum filtered through several filters of different material and pore size to collect all DNA from the sample without clogging the filters. Bottles are sterile, and cross-contamination is not expected because the sampling is done in unique point. First, the whole volume of each sample was filtered through 10-μm nylon filters; then, each sample was divided in four and filtered through 0.8-μm polyethersulfone (PES) filters. Finally, the whole volume of each sample was pooled together and filtered through 0.2-μm PES filters. In total, 18 filters (3 replicates × 6 filters each one) were obtained, which were thereafter preserved with 96% ethanol until eDNA extraction.
DNA and Bioinformatics Analysis
DNA was extracted from the filters using PowerWater® DNA Isolation Kit (MOBIO Laboratories, EE. UU.) following the manufacturer’s extraction protocol. In addition, ethanol was centrifuged, and the pellet, with precipitated DNA, was added in the extracting process. A negative control of pure water was added at this step to monitor contamination during the extraction process. In addition, bovine serum albumin (BSA) was added to the PCR reactions to increase PCR yields from low-purity templates and to avoid, as much as possible, the effect of inhibitors present in the water. Extracted DNA from all the filters was sent to Macrogen (Seoul, South Korea) were it was quantified and sequenced. The eDNA was quantified by a fluorescence-based method Victor 3 (Picogreen, Invitrogen). The modified universal COI primers mlCOIintF/jgHCO2198 (Leray et al., 2013) were used for PCR amplification of a fragment of ∼313 bp within the mitochondrial gene coding for the COI. The primers were modified to include Illumina sequencing adapters and sample-specific dual indices (i5 and i7) following the Illumina (2013) protocol in which the conditions for the amplicon PCR were changed for the ones described by Leray et al. (2013). The PCR reactions were undertaken by Macrogen, following the protocol Illumina (2013). After constructing the library, it was charged into the platform Illumina MiSeq v3 reagents that generates paired-end sequences (2 × 301). Adapters and indices were removed from the raw data along with reads <36 bp using Scythe and Buffalo software, respectively.
The trimmed data were received from Macrogen in Fastaq format and were further processed using the platform Qiime1 (Caporaso et al., 2010). All the sequence reads were assessed for quality by applying a Phred quality score threshold of 20 and were filtered by length (200 bp ≤ reads ≤400 bp) from the downstream analysis.
The paired-end reads from each sample were merged when they presented an overlapping region of at least 100 bp and <15% of differences within this region. This step is necessary to merge forward and reverse reads. Some errors often appear in the Illumina data as the merging is not always 100%, due to sequencing length and PCR errors. For this, one of the thresholds in the merging algorithm must be the percentage of coincidence within the merging region to avoid artifact sequences when merging them and to improve the following assignment against reference database.
For the taxonomic assignation, Basic Local Alignment Search Tool (BLAST) alignment was performed against NCBI database1 of COI sequences (obtained in 09/2017) filtered for environmental reads and using as threshold criteria: maximum E value = 1 × 10–50 and minimum percent identity = 97.0, which is enough in most of the cases for species identification from COI barcode (Hebert et al., 2003b). Some reads may have multiple BLAST hits at a 100% identity to different species, or the best BLAST hit have no species-level reference available. In those cases, the operational taxonomic unit (OTU) was assigned to a genus level.
OTU tables, a list of OTUs obtained for each sample and the number of sequences assigned to them, were constructed, clustering reads with a 100% identity between them and maintaining all assigned sequences including singletons to retain maximum sensitivity for species detection. The removal of singletons is usually employed to eliminate false positives as proposed by Scott et al. (2018), in the context of species survival, NIS early detection, or marine biosecurity surveillance; a false negative is most costly than a false positive (von Ammon et al., 2018). Sequences of organisms without relevance for the study (e.g., human, insects, terrestrial plants, etc.) were removed from the dataset (Figure 2). Generally, chimera formation is not extensive, although some chimeras can be found. In addition, one sequence per OTU was BLASTed manually to verify the reliability of the pipeline parameters and discard the presence of chimeras (Nilsson et al., 2012). The taxonomic information from the remaining OTUs was checked against the World Register of Marine Species2 and AlgaeBase3. The number of reads per OTU, as a proxy for species abundance, was used to generate OTU rarefaction curves using Vegan package in R software (Oksanen et al., 2013).
Figure 2. Rarefaction curve plotting the number of reads by the number of operational taxonomic units (OTUs) for the sample analyzed using high-throughput sequencing. Interpolated line in red. Graphics with the proportion of different taxonomical groups detected in the lagoon before and after expert check (eliminating the reads assigned to terrestrial and exclusive marine species).
The minimum concentration of the extracted eDNA was 0.076 ng/μl. A total of 581,786 raw paired-end reads were sequenced from positive PCR amplicons (average read length, 301 bp). Following initial filtering, a total of 328,220 reads (56.4%) were retained for downstream analysis (average read length, 269 bp). After merging, and a second quality filtering, 239,842 sequences remained (41.2%), of which 98,282 matched a correct reference sequence from the database (16.9%; average read length, 359 bp) (Table 1).
Table 1. Sequences matched a correct reference sequence from the GenBank database per filter, after merging and quality filtering.
The rarefaction curve of sequences indicated that the sample reached a plateau (Figure 2). This figure represents three 2-L superficial water samples, which were then grouped for the purpose of downstream analysis, hence the single plotted line. This indicates that, in terms of the number of sequences generated, our sequencing had enough depth to detect most of the possible species present within the sample.
After the application of expert check, eliminating the reads assigned to terrestrial and exclusive marine species, 51,624 sequences were left assigned to aquatic organisms at least at class level (Table 1). Of these reads, 10.5% were identified to a species level, while the majority was assigned to a genus or higher level. The presence of singletons may prevent reaching the asymptote in the curve of accumulation of species. In this case, the rare species (singletons and doublets) represent 11.1% (3 of the 27 OTU obtained), which did not impede to reach the asymptote, indicating that, in a number of sequences generated, the sample was adequately sized (Figure 2).
Regarding the biota profile captured from NGS after expert check, the assigned sequences were clustered in 27 different OTUs of which 17 belonged to the Class Dinophyceae (89.5% of the reads), while the rest was divided into four different classes by sequences number order: crustaceans, polychaetes, diatoms, and rotifers (Figure 2).
Considering only the sequences assigned at least at a genus level (28,029 in total), 25 taxa were identified, 10 of them at a species level. The genus with more reads was Protoperidinium (19.4%), followed by Scrippsiella (19.1%), and Levanderina (16.9%) (Table 1), pointing out at a clear dominance of Dinophyceae in the water samples analyzed not only in the number of taxa but also in abundance of reads. The following taxon in abundance of reads was the acorn barnacle Perforatus perforatus (8.3% of reads), then the whip mudworm Polydora cornuta (5.8%). Other two polychaetes, Capitella capitata and Hediste diversicolor, and another barnacle (Sacculina carcini) were represented by less reads: 3, 4, and 46, respectively. One rotifer accounting for 5% of the reads was identified: Testudinella clypeata. Diatoms were represented by only two species, Haslea crucigera and Cylindrotheca fusiformis, being together <0.04% of the reads. The copepod Acartia tonsa was also found at very low abundance in number of reads, which is only one read (Table 1).
One half of the species detected from metabarcoding in the analyzed lagoon can be considered nuisance species or indicators of bad environmental quality4 (Table 2). The three polychaetes are bioindicators of pollution. The polychaete P. cornuta and the copepod A. tonsa are considered invasive species in the Mediterranean Sea. Finally, the dinoflagellate Akashiwo sanguinea forms algal blooms (Table 2).
Moreover, seven of the dinoflagellate genera are known to contain numerous species that are considered HABs because they produce toxins and/or form algal blooms, being cataloged in the HAB list of Intergovernmental Oceanographic Commission (IOC)-Unesco (Moestrup et al., 2009). According to the Harmful Algal Events Database (HAEDAT)4, the genera Alexandrium (13 species) and Dynophysis (10 species) found in our study (Table 1) are responsible for most algal blooms in the area. The other five genera of dinoflagellates found in our study and listed as HABs were Karenia (10 species), Prorocentrum (12 species), Karlodinium (6 species), Protoceratium (1 species), and Azadinium (3 species) (Table 1).
The present study, based on a single metabarcode and 6 L of surface water from only one sampling point from one lagoon of 600 ha, did reveal bioindicators of three different environmental or biosecurity threats in Canet-Saint Nazaire lagoon: pollution (three bioindicator species), biological invasions (two species), and harmful algal blooms (one species). Despite extreme methodological simplicity, the resource-efficient methodology employed in this study was capable of obtaining important environmental information. Although it is known that DNA remains in the environment for prolonged periods especially in cold conditions and that individuals may not be present at the time of sampling, these results confirm the utility of COI as an eDNA metabarcoding tool for early alert of biological risks as proposed by Borrell et al. (2017) for biological invasions in ports, expanding the applications to pollution assessments and HABs.
The results indicate a poor ecological status for the lagoon Canet-Saint Nazaire. The three polychaete species were found to have been previously used as pollution indicators due to their ability to survive in highly polluted waters, with wide ranges of temperature and salinity, hypoxia, and eutrophication states (Surugiu, 2005; Dean, 2008; Maranho et al., 2014). P. cornuta is also considered an invasive species in the Mediterranean due to its ability to change the composition and abundance of native species, replacing native species like C. capitata (Cinar et al., 2005; Çinar et al., 2012; Zenetos et al., 2012). It was first discovered in the Mediterranean from polluted sediments in Valencia’s harbor (Spain) (Tena et al., 1991); however, it now has a wider distribution. Its larval phase can survive transport in ship ballast water, and adults may form dense settlements through fouling of ship hulls facilitating its spread (Radashevsky and Selifonova, 2013). On the other hand, the invasive F. enigmaticus that was found in a previous survey in this lagoon (Ardura and Planes, 2017) could have facilitated the introduction and establishment of P. cornuta in the lagoon as reported in other locations (Read and Gordon, 1991; Heiman et al., 2008).
Acartia tonsa deserves additional attention, since it is a neritic copepod commonly found in the western Atlantic, the Pacific, and the Indian Ocean, especially in estuarine and coastal waters. It was first described in the Mediterranean in 1985 (Gaudy and Viñas, 1985) and has been previously reported in the Gulf of Lion in Berre Lagoon (Cervetto et al., 1999). It is considered an invasive species due to its ability to replace native species of copepods and its propensity to spread to new areas through ballast water (David et al., 2007).
On the other hand, despite having a single sampling point with only 6 L of water and forming only a limited description of the diversity of species in the whole ecosystem; the data set shows a great diversity of dinoflagellates in the sampling point analyzed. Most reads (89.5%) and OTUs (60.7%) detected were assigned to the class Dinophycea. This could be explained by the sampling protocol, since the samples were taken from the lagoon surface where a greater number of these organism are expected. In any case, this is a signal of the lagoon’s poor ecological state and high eutrophication and nitrogen and phosphorus levels that facilitate the proliferation of these organisms (Carlier et al., 2008; Heisler et al., 2008). Moreover, several dinoflagellates can produce resistance cysts, allowing them to survive during harsher periods in the lagoon (Sellner et al., 2003). However, the lack of variability in the COI region used for the species identification does not allow differences between species or strains (native or not) from the same genus to be resolved, making it impossible to detect species that can cause harm to the ecosystem. Several species of the genera detected (e.g., Alexandrium, Dinophysis, Karenia, or Prorocentrum) are known to cause HABs that could harm the lagoon ecosystem (Sellner et al., 2003). A. sanguinea (only species in the genus Akashiwo) blooms have been reported to produce powerful surfactant-like proteins that can coat bird plumage and collapse their feathers, causing loss of waterproofing and thermal insulation and even mass mortality (Jones et al., 2017). This could have serious effects in the lagoon because it is located along one of Europe’s principal migration routes, and since March 2006, it has been classified as a SPA (Ardura and Planes, 2017). Further analysis with more specific identification tools would be necessary to assess the risks to the lagoon, but the data obtained in this study should be taken as warning signal. In Thau Lagoon, 100 km away from Canet, a toxic strain from Alexandrium catenella originating in the West Pacific was probably introduced from ballast water released in Séte’s harbor that is connected to the lagoon (Lilly et al., 2002).
The current study has some limitations, likely because the sampling was limited to only one point and 6 L of water taken from the surface; thus, benthic organisms are unlikely detected. One limitation could be a relatively modest number of reads recovered (0.5 million reads). The presence of inhibitors due to the polluted state of the lagoon cannot be discarded as an additional explanation for this, since pollution can inhibit the PCR amplification of eDNA (Jane et al., 2015). The lack of detection of some invasive macroscopic species detected by Ardura and Planes (2017) in the lagoon (the bivalves Abra sp. and Cerastoderma glaucum and the tubeworm Ficopomatus enigmaticus) could be explained by the different sampling methods employed in the two studies. The species C. glaucum occurs in the bottom and is very scarce in the lagoon (in the previous survey, more sampling points were analyzed and only one individual for was obtained). Abra sp. and F. enigmaticus, which were abundant in the lagoon, were not seen where the water samples were taken in the present study; thus, it is possible that their DNA was simply absent from our water samples. Additional explanations could be very degraded DNA (Barnes et al., 2014) and/or lack of universality of the primers that may not anneal equally well in all taxa (Wilcox et al., 2013). In Ardura and Planes (2017), two different barcodes were amplified, e.g., COI (Geller et al., 2013) and 16S rRNA (Palumbi, 1996), but high-quality sequences for F. enigmaticus and Abra sp. were obtained only for the 16S rRNA marker, while in the present study, only COI primers were employed.
The number of detected OTUs in metabarcoding samples tends to be higher than in the visual and barcoding analysis, and effectively, the number of species detected here was much higher than in the previous work by Ardura and Planes (2017). False positives in metabarcoding due to incorrect assignations during the bioinformatic analysis (Ficetola et al., 2015) could be reasonably discarded since all the OTUs detected in this study were BLASTed manually to avoid incorrect assignations. On the other hand, it is possible to detect species that are not alive in the sampling area but whose DNA could come from the sea or from neighboring rivers (Pochon et al., 2017). If it was the case, they would still be around and be able to reach the lagoon eventually. In general, and despite possible errors and the obvious limitations of the simplified sampling scheme, the metabarcoding technique used in this study allowed the successful completion of a general assessment of the lagoon’s environmental state from a list of nuisance and indicator species. Despite this, the data should be considered qualitative rather than empirical. Although eDNA concentration and number of sequences yielded from NGS are positively correlated with biomass or population density, the accurate conversion of genome abundances to cell numbers and estimates of absolute abundance are still imprecise (Kelly et al., 2014; Bonk et al., 2018). A greater number of reads assigned to a specific OTU do not necessarily imply a greater abundance (Thomas et al., 2016).
This study provides new results from a new case study that support the applicability of the COI barcode as a metabarcoding tool using eDNA (e.g., Thomsen and Willerslev, 2015; Zaiko et al., 2015; Borrell et al., 2017; Pochon et al., 2017; Forster et al., 2019), showing in this particular case its usefulness to detect HAB-causing organisms. However, although it is useful to assign some organisms at species level, such as the genus Akashiwo, in most cases, this assignment is not possible because this small fragment, 300 bp, is not informative enough.
The simple sampling scheme in this study, with only one-point sample and three replicates, provides a proof of concept for the use of single samples as an approximation of ecological status. It is important to remark that the study trials a simplified sampling method to remark the value that small sample numbers can have on remote locations for future studies: even one sample with three replicates is informative. However, the limitations associated to collect small sample numbers should not be forgotten, bringing attention to considerations that must be had when analyzing NGS data collected from one sample. In these cases, the data must be analyzed as proof of concept. Further research with additional samples (at different depths and situational periods) would be necessary to build a more complete survey of the study area. Moreover, larger COI fragment or different molecular markers and the enlargement and improvement of barcode databases (Ardura, 2018) are also necessary and will provide a better representation of the ecosystem’s biodiversity.
On the other hand, new methodological tools are been developed to improve biological surveys, such as the use of environmental RNA (eRNA) (Pochon et al., 2017). Compared to eDNA, eRNA degrades more rapidly in the marine environment (typically hours to days, Thomsen et al., 2012; Sassoubre et al., 2016) and is therefore considered a better proxy for detecting living biota. However, susceptibility of RNA makes it difficult to work with. The collection of RNA samples requires dedicated sampling protocols, more careful preservation, and storage. There is also additional processing time and costs associated with isolation and reverse transcription of RNA (Laroche et al., 2016), making it more expensive and challenging and thus a less attractive molecule to work with (Zaiko et al., 2018).
As a final remark, the presence of invasive species, bioindicators of contamination, and organisms that can cause damage to the ecosystem in this single sampling point suggests a high level of degradation and vulnerability of the lagoon. This is a protected area under the Natura 2000 network, and the results may question the effectiveness of the current conservation programs, indicating that COI and metabarcoding can play an important role in monitoring the progress of conservation efforts.
Data Availability Statement
The Fastaq files are available in GenBank, reference (PRJNA549294).
AA: conception and design. SP: sampling. MS-M and AA: laboratory analysis. MS-M: NGS analysis. EG-V and AA: statistical analysis. AA and EG-V: writing the article. Finally, all authors revised the manuscript.
This study has been supported by the Spanish Ministry of Economy and Competitiveness Grant CGL2016-79209-R and the project FUO-201-17 supported by the Centre National de la Recherche Scientifique: “eDNA in Mediterranean lagoons.” This is a contribution from the Marine Observatory of Asturias (OMA), Asturias University Institute of Biotechnology (IUBA), and the Spanish Research Group of Excellence ARENA. AA holds a Juan de la Cierva-reincorporation Grant from the Spanish Ministry.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The authors thank Martin Desmalades for helping with the sampling tasks.
- ^ https://www.ncbi.nlm.nih.gov/
- ^ http://www.marinespecies.org/
- ^ http://www.algaebase.org/
- ^ http://haedat.iode.org/ accessed on August 2019
Anderson, D. M., Alpermann, T. J., Cembella, A. D., Collos, Y., Masseret, E., and Montresor, M. (2012). The globally distributed genus Alexandrium: multifaceted roles in marine ecosystems and impacts on human health. Harmful Algae 14, 10–35. doi: 10.1016/j.hal.2011.10.012
Ardura, A. (2018). Species-specific markers for early detection of marine invertebrate invaders through eDNA methods: gaps and priorities in GenBank as database example. J. Nat. Conserv. 47, 51–57. doi: 10.1016/j.jnc.2018.11.005
Ardura, A., and Planes, S. (2017). Rapid assessment of non-indigenous species in the era of the eDNA barcoding: a mediterranean case study. Estuar.Coast. Shelf Sci. 188, 81–87. doi: 10.1016/j.ecss.2017.02.004
Ardura, A., Zaiko, A., Borrell, Y., Samuiloviene, A., and Garcia-Vazquez, E. (2016). Novel tools for early detection of a global aquatic invasive, the zebra mussel Dreissena polymorpha. Aquat. Conserv. 27, 165–176. doi: 10.1002/aqc.2655
Barnes, M. A., Turner, C. R., Jerde, C. L., Renshaw, M. A., Chadderton, W. L., and Lodge, D. M. (2014). Environmental conditions influence eDNA persistence in aquatic systems. Environ. Sci. Technol. 48, 1819–1827. doi: 10.1021/es404734p
Blackman, R. C., Constable, D., Hahn, C., Sheard, A., and Durkota, J. (2017). Detection of a new non-native freshwater species by DNA metabarcoding of environmental samples – first record of Gammarus fossarum in the UK. Aquat. Invasion. 12, 177–189. doi: 10.3391/ai.2017.12.2.06
Bonk, F., Popp, D., Harms, H., and Centler, F. (2018). PCR-based quantification of taxa-specific abundances in microbial communities: quantifying and avoiding common pitfalls. J. Microbiol. Methods 153, 139–147. doi: 10.1016/j.mimet.2018.09.015
Borrell, Y. J., Miralles, L., Do Huu, H., Mohammed-Geba, K., and Garcia-Vazquez, E. (2017). DNA in a bottle—Rapid metabarcoding survey for early alerts of invasive species in ports. PLoS One 12:e0183347. doi: 10.1371/journal.pone.0183347
Bravo, I., Fernández, M. L., Ramilo, I., and Martınez, A. (2001). Toxin composition of the toxic dinoflagellate Prorocentrum lima isolated from different locations along the Galician coast (NW Spain). Toxicon 39, 1537–1545. doi: 10.1016/s0041-0101(01)00126-x
Bylemand, J., Gleeson, D. M., Duncan, R. P., Hardy, C. M., and Furlan, E. M. (2019). A performance evaluation of targeted eDNA and eDNA metabarcoding analyses for freshwater fishes. Environ. DNA 1, 402–414. doi: 10.1002/edn3.41
Caporaso, J. G., Kuczynski, J., Stombaugh, J., Bittinger, K., Bushman, F. D., Costello, E. K., et al. (2010). QIIME allows analysis of high-throughput community sequencing data. Nat. Methods 7, 335–336. doi: 10.1038/nmeth.f.303
Carlier, A., Riera, P., Amouroux, J. M., Bodiou, J. Y., Desmalades, M., and Grémare, A. (2008). Food web structure of two Mediterranean lagoons under varying degree of eutrophication. J. Sea Res. 60, 264–275. doi: 10.1016/j.seares.2008.10.006
Cataudella, S., Crosetti, D., and Massa, F. (2015). “Mediterranean coastal lagoons: sustainable management and interactions among aquaculture, capture fisheries and the environment,” in General Fisheries Commission for the Mediterranean. Rome: FAO.
Cervetto, G., Gaudy, R., and Pagano, M. (1999). Influence of salinity on the distribution of Acartia tonsa (Copepoda, Calanoida). J. Exp. Mar. Biol Ecol. 239, 33–45. doi: 10.1016/s0022-0981(99)00023-4
Cinar, M. E., Ergen, Z., Dagli, E., and Petersen, M. E. (2005). Alien species of spionid polychaetes (Streblospio gynobranchiata and Polydora cornuta) in Izmir Bay, eastern Mediterranean. J. Mar. Biol. Assoc. UK. 85, 821–827. doi: 10.1017/s0025315405011768
Çinar, M. E., Katagan, T., Öztürk, B., Bakir, K., and Dagli, E. (2012). Spatio-temporal distributions of zoobenthos in soft substratum of Izmir Bay (Aegean Sea, eastern Mediterranean), with special emphasis on alien species and ecological quality status. J. Mar. Biol. Assoc. UK 92, 1457–1477. doi: 10.1017/s0025315412000264
David, M., Gollasch, S., Cabrini, M., Perkoviè, M., Bošnjak, D., and Virgilio, D. (2007). Results from the first ballast water sampling study in the Mediterranean Sea–the Port of Koper study. Mar. Pollut. Bull. 54, 53–65. doi: 10.1016/j.marpolbul.2006.08.041
Ficetola, G. F., Pansu, J., Bonin, A., Coissac, E., Giguet−Covex, C., De Barba, M., et al. (2015). Replication levels, false presences and the estimation of the presence/absence from eDNA metabarcoding data. Mol. Ecol. Resour. 15, 543–556. doi: 10.1111/1755-0998.12338
Forster, D., Filker, S., Kochems, R., Breiner, H.-W., Cordier, T., Pawlowski, J., et al. (2019). A comparison of different ciliate metabarcode genes as bioindicators for environmental impact assessments of salmon aquaculture. Eukaryot. Microbiol. 66, 294–308. doi: 10.1111/jeu.12670
Geller, J., Meyer, C., Parker, M., and Hawk, H. (2013). Redesign of PCR primers for mitochondrial cytochrome c oxidase subunit I for marine invertebrates and application in all−taxa biotic surveys. Mol. Ecol. Resour. 13, 851–861. doi: 10.1111/1755-0998.12138
Hebert, P. D., Ratnasingham, S., and de Waard, J. R. (2003b). Barcoding animal life: cytochrome c oxidase subunit 1 divergences among closely related species. Proc. R. Soc. Lon. B Biol. Sci. 270(Suppl. 1), S96–S99.
Heiman, K. W., Vidargas, N., and Micheli, F. (2008). Non-native habitat as home for non-native species: comparison of communities associated with invasive tubeworm and native oyster reefs. Aquati. Biol. 2, 47–56. doi: 10.3354/ab00034
Heisler, J., Glibert, P. M., Burkholder, J. M., Anderson, D. M., Cochlan, W., Dennison, W. C., et al. (2008). Eutrophication and harmful algal blooms: a scientific consensus. Harmful Algae 8, 3–13. doi: 10.1016/j.hal.2008.08.006
Jane, S. F., Wilcox, T. M., McKelvey, K. S., Young, M. K., Schwartz, M. K., Lowe, W. H., et al. (2015). Distance, flow and PCR inhibition: eDNA dynamics in two headwater streams. Mol. Ecol. Resour. 15, 216–227. doi: 10.1111/1755-0998.12285
Jones, T., Parrish, J. K., Punt, A. E., Trainer, V. L., Kudela, R., Lang, J., et al. (2017). Mass mortality of marine birds in the Northeast Pacific caused by Akashiwo sanguinea. Mar. Ecol. Progr. Ser. 579, 111–127. doi: 10.3354/meps12253
Katsanevakis, S., Coll, M., Piroddi, C., Steenbeek, J., Ben Rais Lasram, F., Zenetos, A., et al. (2014). Invading the mediterranean sea: biodiversity patterns shaped by human activities. Front. Mar. Sci. 1:32. doi: 10.3389/fmars.2014.00032
Laroche, O., Wood, S. A., Tremblay, L. A., Ellis, J. I., Lejzerowicz, F., Pawlowski, J., et al. (2016). First evaluation of foraminiferal metabarcoding for monitoring environmental impact from an offshore oil drilling site. Mar. Environ. Res. 120, 225–235. doi: 10.1016/j.marenvres.2016.08.009
Leray, M., Yang, J. Y., Meyer, C. P., Mills, S. C., Agudelo, N., Ranwez, V., et al. (2013). A new versatile primer set targeting a short fragment of the mitochondrial COI region for metabarcoding metazoan diversity: application for characterizing coral reef fish gut contents. Front. Zool. 10:34. doi: 10.1186/1742-9994-10-34
Lilly, E. L., Kulis, D. M., Gentien, P., and Anderson, D. M. (2002). Paralytic shellfish poisoning toxins in France linked to a human-introduced strain of Alexandrium catenella from the western Pacific: evidence from DNA and toxin analysis. J. Plankton Res. 24, 443–452. doi: 10.1093/plankt/24.5.443
Lin, S., Zhang, H., Hou, Y., Zhuang, Y., and Miranda, L. (2009). High-level diversity of dinoflagellates in the natural environment, revealed by assessment of mitochondrial cox1 and cob genes for dinoflagellate DNA barcoding. Appl. Environ. Microbiol. 75, 1279–1290. doi: 10.1128/AEM.01578-08
Maranho, L. A., Baena-Nogueras, R. M., Lara-Martín, P. A., DelValls, T. A., and Martín-Díaz, M. L. (2014). Bioavailability, oxidative stress, neurotoxicity and genotoxicity of pharmaceuticals bound to marine sediments. The use of the polychaete Hediste diversicolor as bioindicator species. Environ. Res. 134, 353–365. doi: 10.1016/j.envres.2014.08.014
Moestrup, Ø, Akselmann, R., Fraga, S., Hoppenrath, M., Iwataki, M., Komárek, J., et al. (eds) (2009). IOC-UNESCO Taxonomic Reference List of Harmful Micro Algae. Avaliable at: http://www.marinespecies.org/hab (accessed December 16, 2018).
Moore, S. K., Trainer, V. L., Mantua, N. J., Parker, M. S., Laws, E. A., Backer, L. C., et al. (2008). Impacts of climate variability and future climate change on harmful algal blooms and human health. Environ. Health 7(Suppl. 2):S4.
Nilsson, R. H., Tedersoo, L., Abarenkov, K., Ryberg, M., Kristiansson, E., Hartmann, M., et al. (2012). Five simple guidelines for establishing basic authenticity and reliability of newly generated fungal ITS sequences. MycoKeys 4, 37–63. doi: 10.3897/mycokeys.4.3606
Paavola, M., Olenin, S., and Leppäkoski, E. (2005). Are invasive species most successful in habitats of low native species richness across European brackish water seas? Estuar.Coast. Shelf Sci. 64, 738–750. doi: 10.1016/j.ecss.2005.03.021
Pérez-Ruzafa, A., Marcos, C., and Pérez-Ruzafa, I. M. (2011). Mediterranean coastal lagoons in an ecosystem and aquatic resources management context. Phys,. Chem. Earth Parts A/B/C 36, 160–166. doi: 10.1016/j.pce.2010.04.013
Pochon, X., Zaiko, A., Fletcher, L. M., Laroche, O., and Wood, S. A. (2017). Wanted dead or alive? Using metabarcoding of environmental DNA and RNA to distinguish living assemblages for biosecurity applications. PLoS One 12:e0187636. doi: 10.1371/journal.pone.0187636
Read, G. B., and Gordon, D. P. (1991). Adventive occurrence of the fouling serpulid Ficopomatus enigmaticus (Polychaeta) in New Zealand. N. Z. J Mar. Freshwa. Res. 25, 269–273. doi: 10.1080/00288330.1991.9516478
Ruiz, G. M., Fofonoff, P. W., Carlton, J. T., Wonham, M. J., and Hines, A. H. (2000). Invasion of coastal marine communities in North America: apparent patterns, processes, and biases. Annu. Rev. Ecol. Syst. 31, 481–531. doi: 10.1146/annurev.ecolsys.31.1.481
Sassoubre, L. M., Yamahara, K. M., Gardner, L. D., Block, B. A., and Boehm, A. B. (2016). Quantification of environmental DNA (eDNA) shedding and decay rates for three marine fish. Environ. Sci. Technol. 50, 10456–10464. doi: 10.1021/acs.est.6b03114
Scott, R., Zhan, A., Brown, E. A., Chain, F. J. J., Cristescu, M. E., Gras, R., et al. (2018). Optimization and performance testing of a sequence processing pipeline applied to detection of nonindigenous species. Evol. Appl. 11, 891–905. doi: 10.1111/eva.12604
Souchu, P., Bec, B., Smith, V. H., Laugier, T., Fiandrino, A., Benau, L., et al. (2010). Patterns in nutrient limitation and chlorophyll a along an anthropogenic eutrophication gradient in French Mediterranean coastal lagoons. Can. J. Fish Aquat. Sci. 67, 743–753. doi: 10.1139/f10-018
Surugiu, V. (2005). The use of polychaetes as indicators of eutrophication and organic enrichment of coastal waters: a study case–Romanian Black Sea Coast. Analele Ştiinţifice ale Universităţii “Al.I. Cuza”. Iaşi, s. Biol. Anim. 51, 55–62.
Taberlet, P., Coissac, E., Pompanon, F., Brochmann, C., and Willerslev, E. (2012). Towards next−generation biodiversity assessment using DNA metabarcoding. Mol. Ecol. 21, 2045–2050. doi: 10.1111/j.1365-294X.2012.05470.x
Tena, J., Capaccioni-Azzati, R., Porras, R., and Torres-Gavilá, F. J. (1991). Cuatro especies de poliquetos nuevas para las costas mediterráneas españolas en los sedimentos del antepuerto de Valencia. Miscel⋅ Zool. 15, 29–41.
Thomas, A. C., Deagle, B. E., Eveson, J. P., Harsch, C. H., and Trites, A. W. (2016). Quantitative DNA metabarcoding: improved estimates of species proportional biomass using correction factors derived from control material. Mol. Ecol. Resour. 16, 714–726. doi: 10.1111/1755-0998.12490
Thomsen, P. F., and Willerslev, E. (2015). Environmental DNA – An emerging tool in conservation for monitoring past and present biodiversity. Biol. Conserv. 183, 4–18. doi: 10.1016/j.biocon.2014.11.019
Thomsen, P. T., Kielgast, J., Iversen, L., Møller, P. R., Rasmussen, M., and Willerslev, E. (2012). Detection of a diverse marine fish fauna using environmental DNA from seawater samples. PLoS One 7:e41732. doi: 10.1371/journal.pone.0041732
von Ammon, U., Wood, S. A., Larocge, O., Zaiko, A., Tait, L., Lavery, S., et al. (2018). Combining morpho-taxonomy and metabarcoding enhances the detection of non-indigenous marine pests in biofouling communities. Sci. Rep. 8, 16290. doi: 10.1038/s41598-018-34541-1
Vouvé, F., Buscail, R., Aubert, D., Labadie, P., Chevreuil, M., Canal, C., et al. (2014). Bages-Sigean and Canet-St Nazaire lagoons (France): physico-chemical characteristics and contaminant concentrations (Cu, Cd, PCBs and PBDEs) as environmental quality of water and sediment. Environ. Sci. Pollut. Res. 21, 3005–3020. doi: 10.1007/s11356-013-2229-1
Weigand, A., Beermann, A. J., Èiampor, F., Costa, F. O., Csabai, Z., Duarte, S., et al. (2019). DNA barcode reference libraries for the monitoring of aquatic biota in Europe: Gap-analysis and recommendations for future work. Sci. Total Environ. 678, 499–524. doi: 10.1016/j.scitotenv.2019.04.247
Wilcox, T. M., McKelvey, K. S., Young, M. K., Jane, S. F., Lowe, W. H., Whiteley, A. R., et al. (2013). Robust detection of rare species using environmental DNA: the importance of primer specificity. PLoS One 8:e59520. doi: 10.1371/journal.pone.0059520
Yasumoto, T., Oshima, Y., Fukuyo, Y., Oguri, H., Igarashi, T., and Fujita, N. (1980). Identification of Dinophysis fortii as the causative organism of diarrhetic shellfish poisoning. Bull. Japan. Soc. Sci. Fish. 46, 1405–1411. doi: 10.2331/suisan.46.1405
Zaiko, A., Pochon, X., Garcia-Vazquez, E., Olenin, S., and Wood, S. A. (2018). Advantages and limitations of environmental DNA/RNA tools for marine biosecurity: management and surveillance of non-indigenous species. Front. Mar. Sci. 5:322. doi: 10.3389/fmars.2018.00322
Zaiko, A., Samuiloviene, A., Ardura, A., and Garcia-Vazquez, E. (2015). Metabarcoding approach for nonindigenous species surveillance in marine coastal waters. Mar. Pollut. Bull. 100, 53–59. doi: 10.1016/j.marpolbul.2015.09.030
Zenetos, A., Gofas, S., Morri, C., Rosso, A., Violanti, D., García Raso, J. E., et al. (2012). Alien species in the mediterranean Sea by 2012,” in a contribution to the application of european union’s marine strategy framework directive (MSFD), Part2. introduction trends and pathways. Mediterr. Mar. Sci. 13, 328–352. doi: 10.12681/mms.327
Keywords: COI “barcode,” next generation sequencing, Canet lagoon, Nature 2000 areas, conservation
Citation: Suarez-Menendez M, Planes S, Garcia-Vazquez E and Ardura A (2020) Early Alert of Biological Risk in a Coastal Lagoon Through eDNA Metabarcoding. Front. Ecol. Evol. 8:9. doi: 10.3389/fevo.2020.00009
Received: 03 October 2019; Accepted: 15 January 2020;
Published: 28 February 2020.
Edited by:David S. Thaler, Universität Basel, Switzerland
Reviewed by:Alejandro Trujillo Gonzalez, University of Canberra, Australia
Rebecca Barr Simmons, University of North Dakota, United States
Copyright © 2020 Suarez-Menendez, Planes, Garcia-Vazquez and Ardura. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Alba Ardura, email@example.com