Next Article in Journal
Essential Oils from Fruit and Vegetables, Aromatic Herbs, and Spices: Composition, Antioxidant, and Antimicrobial Activities
Previous Article in Journal
Abnormal Enhancement of Protein Disulfide Isomerase-like Activity of a Cyclic Diselenide Conjugated with a Basic Amino Acid by Inserting a Glycine Spacer
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Biology
Volume 10
Issue 11
10.3390/biology10111092
biology-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Testing the Influence of Incomplete DNA Barcode Libraries on Ecological Status Assessment of Mediterranean Transitional Waters

by
Maurizio Pinna
1,2,
Benedetta Saccomanno
1,
Gabriele Marini
1,2,
Francesco Zangaro
1,2,
Akbota Kabayeva
3,
Mina Khalaj
4,
Laura Shaimardan
3,
Simona D’Attis
1,
Eftychia Tzafesta
1 and
Valeria Specchia
1,*
1
Department of Biological and Environmental Sciences and Technologies, University of Salento, DiSTeBA, Via Monteroni 165, 73100 Lecce, Italy
2
Research Centre for Fisheries and Aquaculture of Aquatina di Frigole, DiSTeBA, University of Salento, 73100 Lecce, Italy
3
Department of Biotechnology, Al-Farabi Kazakh National University, Al-Farabi Avenue, Almaty 71-050040, Kazakhstan
4
Department of Basic Sciences, Sari Agricultural Sciences and Natural Resources University, 48 Sari, Iran
*
Author to whom correspondence should be addressed.
Biology 2021, 10(11), 1092; https://doi.org/10.3390/biology10111092
Submission received: 22 September 2021 / Revised: 14 October 2021 / Accepted: 20 October 2021 / Published: 25 October 2021
(This article belongs to the Section Conservation Biology and Biodiversity)
Browse Figures
Versions Notes

Abstract

:

Simple Summary

The biodiversity and ecological status assessment of transitional water ecosystems by benthic macroinvertebrates investigation could use DNA barcode tools for more rapid and efficient outputs. The principal limits of this application are the incompleteness of DNA barcode databases, the identification of optimal primers set, and the gap in the species sequences. The influence of the incompleteness of DNA barcode libraries on species diversity indices, ecological indicators, and ecological status assessment in transitional waters of the southeast Mediterranean were analysed, underlying the importance to implement DNA barcode libraries and to put an effort toward specific species at a local level.

Abstract

The ecological assessment of European aquatic ecosystems is regulated under the framework directives on strategy for water and marine environments. Benthic macroinvertebrates are the most used biological quality element for ecological assessment of rivers, coastal-marines, and transitional waters. The morphological identification of benthic macroinvertebrates is the current tool for their assessment. Recently, DNA-based tools have been proposed as effective alternatives. The main current limits of DNA-based applications include the incompleteness of species recorded in the DNA barcode reference libraries and the primers bias. Here, we analysed the influence of the incompleteness of DNA barcode databases on species diversity indices, ecological indicators, and ecological assessment in transitional waters of the southeast Mediterranean, taking into account the availability of commonly sequenced and deposited genomic regions for listed species. The ecological quality status assigned through the potential application of both approaches to the analysed transitional water ecosystems was different in 27% of sites. We also analysed the inter-specific genetic distances to evaluate the potential application of the DNA metabarcoding method. Overall, this work highlights the importance to expand the barcode databases and to analyse, at the regional level, the gaps in the DNA barcodes.

1. Introduction

A key global challenge in the 21st century is to maintain the supply of clean water and other aquatic ecosystem services without affecting the supporting biodiversity and ecosystem processes that underpin their sustainability. Accordingly, extensive national and international regulations have been adopted to protect aquatic ecosystems and water resources, including the Water Framework Directive (EC, 2000, Directive 2000/60/EC-WFD), the Marine Strategy Framework Directive (EC, 2008, Directive 2008/56/EC-MSFD), the Swiss Water Protection Ordinance (WPO, Swiss Federal Council 1998), the Clean Water Act (CWA, 1972) of the US Environmental Protection Agency, and the United Nations Convention on the Law of the Sea (UNCLOS, 1982). All of these regulations aim to protect aquatic ecosystems and restore altered ecosystems at least to the category of “good status”, which is defined as a condition slightly altered by anthropogenic activities [1]. The WFD indicates specific biological quality elements (BQEs), such as fish, benthic macroinvertebrates, and phytoplanktonic benthos to assess the ecological quality status (EQS). To evaluate these elements, the directive requires the development and validation of tools mainly based on taxonomic identification and composition, abundance, and species sensitivity [2,3,4,5,6,7,8]. They also include descriptors based on the individual body size [9,10,11,12], functional rates, such as the decomposition of dead organic matter [13,14,15,16,17,18,19,20,21], and ecosystem thermodynamics [22].
Among biological quality elements, the benthic macroinvertebrates are effective biological indicators that respond to environmental changes of aquatic ecosystems [23,24]. Their main characteristics as bio-indicators include poor mobility, high number of species and functional groups, long life-cycles, and their important role in aquatic trophic networks [25,26,27,28]. However, the main weaknesses in the use of benthic macroinvertebrates are the time-consuming process of taxonomic identification and possible misidentifications [29,30,31]. In recent years, DNA-based taxon identification has been emerging as a tool that could improve the accuracy and efficiency of biomonitoring, including the assessment of benthic macroinvertebrate communities [32,33,34,35,36,37,38]. The advantages of DNA-based tools include the identification of immature stages that lack recognizable features and the greater resolution at species level, especially for macroinvertebrates that present a high diversity, in comparison with other groups of species.
The DNA-based method uses specific genomic sequences or “barcodes” that uniquely identify the species [39,40,41]. For animals, the most commonly used and reliable barcode regions belong to the mitochondrial cytochrome c oxidase subunit I gene (COI) and rDNA genes [1,39,42,43,44]. Public database repositories, such as BOLD Systems and NCBI GenBank, currently hold a wide range of species based on the reference sequence. However, the use of DNA barcoding currently finds its weaknesses in the incompleteness of these databases and in the primers bias [35,45,46]. Concerning the gap in the taxonomic coverage of the DNA-reference libraries, an example is that only about 48% of benthic macroinvertebrate species included in the most recent AZTI’s Marine Biotic Index list (AMBI list) present at least one DNA barcode sequence in the reference libraries. Another example underlining the gap in the taxonomic coverage of the DNA-reference libraries is that, for the marine benthic macroinvertebrates of the AMBI list, the three most represented phyla (Annelida, Mollusca, and Arthropoda, which represent about 85% of the total species considered in the list) are moderately represented in the DNA barcode reference libraries (from 40 to 50%). On the contrary, a few groups constituted by Nemertea, Sipuncula, and Echinodermata are more represented (at least 65%) in the reference libraries [45]. Concerning the primers bias, the application of DNA metabarcoding ideally requires the identification of primer pairs that are useful for the detection of all taxa present in an environmental sample [35].
The aim of this research is to evaluate the current potential applicability of DNA-based tools in the assessment of the ecological status of transitional water ecosystems at a regional scale, taking into account the gaps in the barcode databases and the availability of useful primers sets. We used a checklist of benthic macroinvertebrates collected in two seasons from ten southeast Mediterranean transitional waters (Apulia, Italy) to compare and to correlate species diversity, ecological indicators, and assessment scores calculated for morphologically identified taxa and for taxa available in the reference libraries.

2. Materials and Methods

2.1. Benthic Macroinvertebrates Database

For this study, the species checklist of the Apulia Regional Environmental Protection Agency (ARPA-Puglia, Italy) published in 2011 was acquired, consisting of a checklist of benthic macroinvertebrate species from the transitional waters of the Apulia region in the southeast of Italy. An overall number of ten transitional water ecosystems was considered (Figure 1), namely: Laguna di Lesina (1); Lago Varano (2); Vasche Evaporanti, Lago Salpi (3); Torre Guaceto (4); Punta della Contessa (5); Cesine (6); Alimini Grande (7); Baia di Porto Cesareo (8); Mar Piccolo-Primo Seno (9); Mar Piccolo-Secondo Seno (10). A total of fifteen sampling sites were investigated, three in the Laguna di Lesina (1a, 1b, 1c), three in the Lago di Varano (2a, 2b, 2c), two in the Alimini Grande (7a, 7b), and one in each of the other transitional water ecosystems. The data analysed refer to the samplings carried out during the fall of 2010 and the spring of 2011. Each species was defined by its presence or absence in the sampled transitional aquatic ecosystems (Supplementary Table S1). Taxon names of each reported species were verified using both the European platform EU-NOMEN (http://www.eu-nomen.eu) (accessed on 15 March 2021) and the worldwide platform WORMS (http://www.marinespecies.org) (accessed on 15 March 2021). The data matrix was updated to the most recent and accepted taxonomy. DNA barcode libraries (BOLD systems, http://www.boldsystems.org (accessed on 22 March 2021); Ranasignham and Hebert 2007; and NCBI GenBank, https://www.ncbi.nlm.nih.gov/genbank (accessed on 22 March 2021)) were examined in order to produce a reduced database containing all the barcoded species that were needed for our analysis.

2.2. Descriptors, Ecological Indicators, and Correlations between Morphological and Reduced Databases

The species richness and the Shannon diversity index were calculated for both databases. Two ecological indicators, AMBI and M-AMBI, were estimated using the AZTI-AMBI software v6.0 (http://ambi.azti.es, access on 24 March 2021), which was used to assess the ecological quality status of transitional water ecosystems by classifying them into five quality categories (High, Good, Moderate, Poor, and Bad). Correlation tests between ecological indicators of morphological and molecular data were carried out to compare the results.

2.3. Species Delimitation Analysis

COI sequences of the taxa listed in Table S1 were downloaded from BOLD Systems and NCBI using PrimerMiner 0.3b [47,48,49]. We only omitted taxa where the species was not specified or was not present in the database, giving a total of 38 taxa. For each of them, we aligned the sequences of all the accessions available in the databases to generate one consensus sequence, with the only exceptions of Nereis falsa, Eunice vittate, and Actinia fragacea where only one accession was found. These consensus sequences were then aligned to construct a maximum likelihood tree (ML) in RAxml (Randomized Axelerated Maximum Likelihood), using the general time reversible + gamma (GTR+G) model [50] and a neighbor-joining (NJ) bootstrap method with the Kimura-2-parameter model in MEGAX. Both trees were then used as constraint trees for species delimitation with bPTP, the Bayesian implementation of the Poisson Tree Processes (PTP) model with a 100,000 MCMC generation and a 1% burn-in [51].

3. Results

3.1. Presence of Target Species in DNA Barcode Libraries

We examined a total number of 82 benthic macroinvertebrate species within six phyla from ten Apulian transitional water sites (Figure 1 and Table S1) that were reported in the ARPA-Puglia checklist. Among these, only 53 species, corresponding to the 64% of the benthic macroinvertebrates checklist, presented at least one DNA barcode sequence within the BOLD or NCBI GenBank reference libraries. For the remaining 29 species, (36%), a DNA barcode was not available in the reference libraries. The phyla that show a relevant gap in the reference libraries were Mollusca (45%), Annelida (28%), and Arthropoda (21%). These species were also characterized by a reduced occurrence and low density. In particular, the highest number of species without a DNA barcode was found in two of the sampled transitional water ecosystems: Lago di Varano and Mar Piccolo-Secondo Seno.

3.2. Species Richness

The species richness, being the simplest descriptor of the biodiversity of the ecological communities, which can also be assessed from the metabarcoding data [1], was calculated for morphological and reduced databases, as reported in Figure 2. In the fall, the percentage of species with DNA barcodes ranged from 40% in Alimini Grande (sampling site 7a) to 100% in Torre Guaceto and Punta della Contessa (sampling sites 4 and 5, respectively); 13 of the 15 percentage values obtained by comparing the morphological and the reduced databases were lower than 75%. In the spring, the percentage of species with a DNA barcode ranged from 43% in Mar Piccolo-Secondo Seno (sampling site 10) to 100% in Vasche Evaporanti Lago Salpi, Torre Guaceto, and Punta della Contessa (sampling sites 3, 4 and 5, respectively); 8 of the 15 percentage values obtained comparing the morphological and the reduced databases were lower than 75%. In each season, species-richness values based on both the morphological and the sequence data were directly and significantly correlated (p < 0.001; Figure 3).

3.3. Shannon Diversity Index

The Shannon diversity index takes into account the distribution of individuals among the species and was calculated for both databases. The values obtained from the morphological database were higher than those obtained from the reduced database according to the higher number of species recognized in the first one (Figure 4). Only in a few sampling sites, the index was similar or higher in the reduced database than in the morphological one. In both seasons, the Shannon diversity index values from both databases were directly and significantly correlated (p < 0.001; Figure 5).

3.4. AMBI

According to the morphological database, the AMBI indicator was lower in the fall season than in the spring season. The lowest AMBI values in the fall were recorded in Torre Guaceto (sampling site 4), whereas in the spring, the lowest values were recorded in Alimini Grande (sampling site 7b) (Figure 6). The maximum AMBI index was recorded in the Lago di Varano (sampling site 2c) in the fall and at Punta della Contessa in the spring (sampling site 5) (Figure 6). The reduced database showed higher AMBI values than the morphological ones in both seasons (Figure 6). In each season, the AMBI values calculated from morphological and reduced databases were directly and significantly correlated (p < 0.001; Figure 7).

3.5. M-AMBI

The M-AMBI indicator values were higher in the morphological database than in the reduced database in both seasons and similar between databases (Figure 8). In each season, the M-AMBI indicators from both databases were directly and significantly correlated (p < 0.001; Figure 9).

3.6. Comparison of Ecological Quality Status Assignment to the Transitional Water Ecosystems Using Morphological and DNA Barcode Databases

Using M-AMBI values, we assigned an ecological quality class to each sampling site, for each season (fall and spring), and for each reference database (morphological and DNA barcode database) (Table 1). After 30 comparisons between the morphological and barcode data, eight sites presented quality class divergences. Therefore, in 27% of the sites, the classification of the EQS obtained through the traditional approach differs from that assigned through a potential application of the DNA barcode approach. It is important to underline that the distance between the quality classifications assigned through both approaches differs on a single quality class. In addition, among the eight cases of divergence, the classification result based on the barcode database exceeds only in an instance (sampling site: 9): the one obtained from the morphological approach. This analysis and the linear relationship of the results from both approaches support the great potential of the application of DNA metabarcoding in biomonitoring programs, confirming the importance of populating the databases with barcode sequences.

3.7. Primers and DNA-Barcoded Region Analysis

A fundamental condition to perform a metabarcoding analysis is to identify a primer pair or primer set that anneals and amplifies the DNA of all of the species represented in the environmental sample. We checked the information of the primers used for all the listed and barcoded species in the BOLD Systems database. A primer pair is reported for 39 of the 53 barcoded species (74%), and, among these, 32 (82%) were amplified using LCO1490/HCO2198 or degenerated/modified LCO1490/HCO2198 primers (Supplementary Table S2). For the other 14 barcoded species, we verified the alignment by CLUSTALW multi-alignment with degenerated LCO1490/HCO2198 primers, and we obtained a match for eight species. This analysis outlines the importance of updating the databases with the primers used for the amplification of the DNA-barcoded region and the importance of the local/regional analysis to set the metabarcoding experiments in relation to the expected specific pool of species.

3.8. Interspecific Genetic Distance Analysis

To deepen the efficiency of the application of DNA metabarcoding, we analysed the sequences of the species present in the list to establish the interspecific genetic distances. Comparing the species delimitation results between the maximum likelihood and the neighbor-joining method (Figure 10, Tables S3 and S4), we can see that the first separates 92% of the species taken into account and exhibits Bayesian support values above 0.7 in 41%, while NJ separates 87% of the species with Bayesian support values above 0.7% in 27%. Both methods clustered together Actinia equina with Actinia fragacea and Musculista senhousia with Mytilaster minimus, but they were able to distinguish between Dardanus arrosor and Dardanus calidus, and Gibbula umbilicaris and Gibbula varia. Interestingly, Cerastoderma edule and Cerastoderma glaucum are two separate species in NJ (support value = 1.000) while they are together in ML (support value = 0.506) and, in turn, Branchiomma bairdii and Branchiomma boholense are one species in NJ (0.504) and two separate ones in ML (0.342). In addition, Lekanesphaera hookeri clusters with Sphaeroma serratum in NJ (0.281) and with Cymodoce truncate in ML (0.336).

4. Discussion

The strength of genetic identification methods to assess species biodiversity depends decisively on the completeness of the reference sequence databases. Taxa lacking barcodes in the databases cannot be identified through DNA-based approaches. Many studies show a discrepancy between molecular and morphological datasets, both in terms of species presence and abundance [45,46,47,48].
In the case study proposed in this work, the three phyla with the highest species richness (Mollusca, Arthropoda, and Annelida), covering 93% of total species in the list, present 58%, 73%, and 65% of barcoded species, respectively. The results of the analysis of DNA barcode reference libraries confirm the presence of gaps in the sampled species data coverage, comparing the Apulia Regional Environmental Protection Agency checklist with the sequences available in the BOLD Systems and GenBank, which has revealed a gap in DNA barcode sequences coverage of 42.3% [46]. Additionally, Leite et al. highlighted a lack of representative barcodes for many marine macroinvertebrate species belonging to dominant faunal groups [48].
Our results, in terms of ecological descriptors (species richness and Shannon index diversity), highlighted how the current process of identification and analysis of sampled species biodiversity through the innovative DNA barcoding approach underestimates the real biodiversity of transitional water ecosystems. Also, this trend is amplified in the transitional aquatic ecosystems where the species richness, calculated through the morphological approach, is higher [28].
In addition, 27% of the analysed transitional water ecosystems differed in the ecological quality status assigned through the potential application of both approaches. This confirms the importance of expanding barcode databases and defining useful primers sets for molecular identification at the regional level.
As proof of concept, we tried to discriminate the species present in this study solely by analysing the genetic variation of publicly available COI sequences. We took a phylogeny-based species delimitation approach by comparing two different methods: neighbor-joining and maximum likelihood. The first method finds the best tree in the dataset with a clustering algorithm, while the second considers a set of all the possible trees and selects the best based on the highest log-likelihood (lnL) tree [52]. Both methods produced a similar number of species with NJ separating 33 out of 38 species and ML separating 34 out of 38 species, but they differently designated some of the clusters. With NJ being a rapid and computationally less demanding method, it is often used to quickly separate species based on pairwise distances [53], but its accuracy decreases when applied to a large dataset, short sequences, or an unequal rate of substitution [54,55]. For these reasons, ML is increasingly used in barcoding studies [56,57,58,59,60,61] and can give a deeper resolution in terms of molecular evolution, particularly when considering a heterogeneous or previously uncharacterized dataset. This is also shown in this study where the ML approach shows higher supporting values in species delimitation. However, it is important to underline that the present analysis was obtained by building a consensus sequence from a different number of sequenced accessions found in several parts of the world. This may explain why, even though the species delimitation mostly matches our expectation, many of the reported species delimitation have low support values. In addition, molecular barcoding is often based on the amplification of a small fragment of a conserved gene, which may underestimate genetic variation and divergence time [62], and thus, incorrectly assign different species to the same cluster. We, therefore, demonstrated the importance of not only generating NGS data at a local/regional level, but also of populating reference libraries with several barcoding genes in order to fine-tune and effectively apply molecular biomonitoring.

5. Conclusions

DNA-based methods for assessing an ecosystem’s health are an innovative tool but require further standardization and improvement processes. This work, focused on the transitional aquatic ecosystems of the southeast Mediterranean, underlines both the importance of upgrading global DNA barcode databases and the validity of local databases for a more suitable identification of primer sets. In addition, the interspecific genetic distances analysis is relevant to verify the percentage of species potentially identified in a DNA metabarcoding study; the species with similar sequences require the analysis of other gene markers or species-specific primers for amplification of the barcoded genomic region.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/biology10111092/s1, Table S1: List of taxa identified in the transitional aquatic ecosystems of Apulia (South-East Italy) for biomonitoring purposes (ARPA, 2011), Table S2: List of barcoded species and primers pair reported in the BOLD Systems database, Table S3: bPTP species delimitation based on a maximum likelihood constraint tree generated in RAxml., Table S4: bPTP species delimitation based on a Neighbor/joining constraint tree generated in MEGAX.

Author Contributions

Conceptualization, M.P. and V.S.; methodology, M.P. and V.S.; formal analysis, M.P., B.S. and G.M.; data curation, M.P., B.S., G.M., F.Z., S.D., E.T., A.K., M.K. and L.S.; writing original draft preparation, M.P. and V.S.; writing-review and editing, M.P., V.S., B.S. and F.Z.; supervision, M.P. and V.S.; funding acquisition M.P. and V.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the European Community (Interreg ADRION Programme 2014–2020), ImPrEco project N.450 awarded to M. Pinna by the Italian Ministry of University and Research. Funding of Basic Research Activities (FFABR) 2017 was awarded to M. Pinna and V. Specchia; the APC was funded by FFABR 2017 and was awarded to V. Specchia.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors are grateful to the Inter-Asia Project (CUP F82I17000000008) for supporting A. Akhmetkaliyeva and L. Shaimardan with a fellowship. The authors thank F. Leese for the useful comments and suggestions on the preliminary draft, as well as to those who provided during the meetings of the EU Cost Action CA15219 DNAqua-Net on “Developing new genetic tools for bioassessment of aquatic ecosystems in Europe”.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Pawlowski, J.; Kelly-Quinn, M.; Altermatt, F.; Apothéloz-Perret-Gentil, L.; Beja, P.; Boggero, A.; Borja, A.; Bouchez, A.; Cordier, T.; Domaizon, I.; et al. The future of biotic indices in the ecogenomic era: Integrating (e)DNA metabarcoding in biological assessment of aquatic ecosystems. Sci. Total Environ. 2018, 637–638, 1295–1310. [Google Scholar] [CrossRef]
  2. Borja, A.; Franco, J.; Pérez, V. A Marine Biotic Index to Establish the Ecological Quality of Soft-Bottom Benthos Within European Estuarine and Coastal Environments. Mar. Pollut. Bull. 2000, 40, 1100–1114. [Google Scholar] [CrossRef]
  3. Simboura, N.; Zenetos, A. Benthic indicators to use in Ecological Quality classification of Mediterranean soft bottom marine ecosystems, including a new Biotic Index. Mediterr. Mar. Sci. 2002, 3, 77. [Google Scholar] [CrossRef] [Green Version]
  4. Bellino, A.; Mangano, M.C.; Baldantoni, D.; Russell, B.D.; Mannino, A.M.; Mazzola, A.; Vizzini, S.; Sarà, G. Seasonal patterns of biodiversity in Mediterranean coastal lagoons. Divers. Distrib. 2019, 25, 1512–1526. [Google Scholar] [CrossRef] [Green Version]
  5. Orfanidis, S.; Papathanasiou, V.; Sabetta, L.; Pinna, M.; Gigi, V.; Gounaris, S.; Tsiagga, E.; Nakou, K.; Theodosiou, T.H. Benthic macrophyte communities as bioindicators of transitional and coastal waters: Relevant approaches and tools. Transit. Waters Bull. 2007, 1, 45–49. [Google Scholar]
  6. Orfanidis, S.; Pinna, M.; Sabetta, L.; Stamatis, N.; Nakou, K. Variation of structural and functional metrics in macrophyte communities within two habitats of eastern Mediterranean coastal lagoons: Natural versus human effects. Aquat. Conserv. Mar. Freshw. Ecosyst. 2008, 18, S45–S61. [Google Scholar] [CrossRef]
  7. Ponti, M.; Vadrucci, M.R.; Orfanidis, S.; Pinna, M. Biotic indices for ecological status of transitional water ecosystems. Transit. Waters Bull. 2009, 3, 32–90. [Google Scholar]
  8. Pinto, R.; Patriício, J.; Baeta, A.; Fath, B.D.; Neto, J.M.; Marques, J.C. Review and evaluation of estuarine biotic indices to assess benthic condition. Ecol. Indic. 2009, 9, 1–25. [Google Scholar] [CrossRef] [Green Version]
  9. Elbrecht, V.; Vamos, E.E.; Meissner, K.; Aroviita, J.; Leese, F. Assessing strengths and weaknesses of DNA metabarcoding-based macroinvertebrate identification for routine stream monitoring. Methods Ecol. Evol. 2017, 8, 1265–1275. [Google Scholar] [CrossRef] [Green Version]
  10. Reizopoulou, S.; Nicolaidou, A. Index of size distribution (ISD): A method of quality assessment for coastal lagoons. Hydrobiologia 2007, 577, 141–149. [Google Scholar] [CrossRef]
  11. Logez, M.; Pont, D. Development of metrics based on fish body size and species traits to assess European coldwater streams. Ecol. Indic. 2011, 11, 1204–1215. [Google Scholar] [CrossRef]
  12. Lucena-Moya, P.; Pardo, I. An invertebrate multimetric index to classify the ecological status of small coastal lagoons in the Mediterranean ecoregion (MIBIIN). Mar. Freshw. Res. 2012, 63, 801–814. [Google Scholar] [CrossRef]
  13. Gessner, M.O.; Chauvet, E. A case for using litter breakdown to assess functional stream integrity. Ecol. Appl. 2002, 12, 498–510. [Google Scholar] [CrossRef]
  14. Pinna, M.; Sangiorgio, F.; Fonnesu, A.; Basset, A. Spatial variation of plant detritus processing in a Mediterranean river type: The case of River Tirso basin (Sardinia, Italy). J. Environ. Sci. 2003, 15, 227–240. [Google Scholar]
  15. Kominoski, J.S.; Pringle, C.M.; Ball, B.A.; Bradford, M.A.; Coleman, D.C.; Hall, D.B.; Hunter, M.D. Nonadditive effects of leaf litter species diversity on breakdown dynamics in a detritus-based stream. Ecology 2007, 88, 1167–1176. [Google Scholar] [CrossRef] [PubMed]
  16. Wallace, J.B.; Webster, J.R.; Cuffney, T.F. Stream detritus dynamics: Regulation by invertebrate consumers. Oecologia 1982, 53, 197–200. [Google Scholar] [CrossRef]
  17. Wissinger, S.A.; Perchik, M.E.; Klemmer, A.J. Role of animal detritivores in the breakdown of emergent plant detritus in temporary ponds. Freshw. Sci. 2018, 3, 826–835. [Google Scholar] [CrossRef]
  18. Di Sabatino, A.; Cicolani, B.; Miccoli, F.P.; Cristiano, G. Plant detritus origin and microbial–detritivore interactions affect leaf litter breakdown in a Central Apennine (Italy) cold spring. Aquat. Ecol. 2020, 54, 495–504. [Google Scholar] [CrossRef]
  19. Wu, W.; Huang, H.; Biber, P.; Bethel, M. Litter decomposition of Spartina alterniflora and Juncus roemerianus: Implications of climate change in salt marshes. J. Coast. Res. 2017, 33, 372–384. [Google Scholar] [CrossRef]
  20. Bozinovic, F.; Munoz, J.L.P.; Cruz-Neto, A.P. Intraspecific variability in the basal metabolic rate: Testing the food habits hypothesis. Physiol. Biochem. Zool. 2007, 80, 452–460. [Google Scholar] [CrossRef] [PubMed]
  21. Di Sabatino, A.; Cristiano, G.; Pinna, M.; Lombardo, P.; Miccoli, F.P.; Marini, G.; Vignini, P.; Cicolani, B. Structure, functional organization and biological traits of macroinvertebrate assemblages from leaf-bags and benthic samples in a third-order stream of Central Apennines (Italy). Ecol. Indic. 2014, 46, 84–91. [Google Scholar] [CrossRef]
  22. Silow, E.A.; In-Hye, O. Aquatic ecosystem assessment using exergy. Ecol. Indic. 2004, 4, 189–198. [Google Scholar] [CrossRef]
  23. Rosenberg, D.M.; Resh, V.H. Introduction to Freshwater Biomonitoring and Benthic Macroinvertebrates. In Freshwater Biomonitoring and Benthic Macroinvertebrates; Rosenberg, D.M., Resh, V.H., Eds.; Chapman/Hall: New York, NY, USA, 1993; pp. 1–9. [Google Scholar]
  24. Basset, A.; Barbone, E.; Rosati, I.; Vignes, F.; Breber, P.; Specchiulli, A.; D’Adamo, R.; Renzi, M.; Focardi, S.; Ungaro, N.; et al. Resistance and resilience of ecosystem descriptors and properties to dystrophic events: A study case in a Mediterranean lagoon. Transit. Waters Bull. 2013, 7, 1–22. [Google Scholar]
  25. Dauer, D.M. Biological criteria, environmental health and estuarine macrobenthic community structure. Mar. Pollut. Bull. 1993, 26, 249–257. [Google Scholar] [CrossRef]
  26. Petchey, O.L.; Belgrano, A. Body-size distributions and size-spectra: Universal indicators of ecological status? Biol. Lett. 2010, 6, 434–437. [Google Scholar] [CrossRef] [Green Version]
  27. Sharma, K.K.; Chowdhary, S. Macroinvertebrate assemblages as biological indicators of pollution in a Central Himalayan River, Tawi (JK). Int. J. Biodivers. Conserv. 2011, 3, 167–174. [Google Scholar] [CrossRef]
  28. Pinna, M.; Janzen, S.; Franco, A.; Specchia, V.; Marini, G. Role of habitats and sampling techniques on macroinvertebrate descriptors and ecological indicators: An experiment in a protected Mediterranean lagoon. Ecol. Indic. 2017, 83, 495–503. [Google Scholar] [CrossRef]
  29. Barbone, E.; Rosati, I.; Pinna, M.; Basset, A. Taxonomic and dimensional structure of benthic macroinvertebrate guilds in the Margherita di Savoia Salt Pans (Italy). Transit. Waters Bull. 2007, 1, 21–31. [Google Scholar] [CrossRef]
  30. Pawlowski, M.B.; Branstrator, D.K.; Hrabik, R. Major shift in the phenology of crustacean biomass in western Lake Superior associated with temperature anomaly. J. Great Lakes Res. 2018, 44, 788–797. [Google Scholar] [CrossRef]
  31. Aylagas, E.; Borja, A.; Pochon, X.; Zaiko, A.; Keeley, N.; Bruce, K.; Hong, P.; Ruiz, G.M.; Stein, E.D.; Theroux, S.; et al. Translational Molecular Ecology in practice: Linking DNA-based methods to actionable marine environmental management. Sci. Total Environ. 2020, 744, 140780. [Google Scholar] [CrossRef]
  32. Hajibabaei, M.; Shokralla, S.; Zhou, X.; Singer, G.A.C.; Baird, D.J. Environmental barcoding: A next-generation sequencing approach for biomonitoring applications using river benthos. PLoS ONE 2011, 6, e17497. [Google Scholar] [CrossRef] [Green Version]
  33. Carew, M.E.; Pettigrove, V.J.; Metzeling, L.; Hoffmann, A.A. Environmental monitoring using next generation sequencing: Rapid identification of macroinvertebrate bioindicator species. Front. Zool. 2013, 10, 45. [Google Scholar] [CrossRef] [Green Version]
  34. Elbrecht, V.; Leese, F. Can DNA-Based Ecosystem Assessments Quantify Species Abundance? Testing Primer Bias and Biomass-Sequence Relationships with an Innovative Metabarcoding Protocol. PLoS ONE 2015, 10, e0130324. [Google Scholar] [CrossRef] [PubMed]
  35. Elbrecht, V.; Leese, F. Validation and development of COI metabarcoding primers for freshwater macroinvertebrate bioassessment. Front. Environ. Sci. 2017, 5, 11. [Google Scholar] [CrossRef] [Green Version]
  36. Leese, F.; Altermatt, F.; Bouchez, A.; Ekrem, T.; Hering, D.; Meissner, K.; Mergen, P.; Pawlowski, J.; Piggott, J.J.; Rimet, F.; et al. DNAqua-Net: Developing new genetic tools for bioassessment and monitoring of aquatic ecosystems in Europe. Res. Ideas Outcomes 2016, 2, e11321. [Google Scholar] [CrossRef] [Green Version]
  37. Hering, D.; Borja, A.; Jones, J.I.; Pont, D.; Boets, P.; Bouchez, A.; Bruce, K.; Drakare, S.; Hänfling, B.; Kahlert, M.; et al. Implementation options for DNA-based identification into ecological status assessment under the European water framework directive. Water Res. 2018, 138, 192–205. [Google Scholar] [CrossRef]
  38. Beentjes, K.K.; Speksnijder, A.G.C.L.; Schilthuizen, M.; Schaub, B.E.; Van Der Hoorn, B.B. The influence of macroinvertebrate abundance on the assessment of freshwater quality in The Netherlands. Metabarcoding Metagenom. 2018, 2, e26744. [Google Scholar] [CrossRef] [Green Version]
  39. Hebert, P.D.; Cywinska, A.; Ball, S.L.; de Waard, J.R. Biologcal identifications through DNA barcodes. Proc. Biol. Sci. 2003, 270, 313–321. [Google Scholar] [CrossRef] [Green Version]
  40. Sweeney, B.W.; Battle, J.M.; Jackson, J.K.; Dapkey, T. Can DNA barcodes of stream macroinvertebrates improve descriptions of community structure and water quality? J. N. Am. Benthol. Soc. 2011, 30, 195–216. [Google Scholar] [CrossRef]
  41. Kress, W.J.; García-Robledo, C.; Uriarte, M.; Erickson, D.L. DNA barcodes for ecology, evolution, and conservation. Trends Ecol. Evol. 2015, 30, 25–35. [Google Scholar] [CrossRef] [PubMed]
  42. Aylagas, E.; Borja, A.; Rodriguez-Ezpeleta, N. Environmental status assessment using DNA metabarcoding: Towards a genetics based Marine Biotic Index (gAMBI). PLoS ONE 2014, 9, e90529. [Google Scholar] [CrossRef] [Green Version]
  43. Leray, M.; Knowlton, N. DNA barcoding and metabarcoding of standardized samples reveal patterns of marine benthic diversity. Proc. Natl. Acad. Sci. USA 2015, 112, 2076–2081. [Google Scholar] [CrossRef] [Green Version]
  44. Aylagas, E.; Borja, Á.; Muxika, I.; Rodríguez-Ezpeleta, N. Adapting metabarcoding-based benthic biomonitoring into routine marine ecological status assessment networks. Ecol. Indic. 2018, 95, 194–202. [Google Scholar] [CrossRef]
  45. Weigand, H.; Beermann, A.J.; Čiampor, F.; Costa, F.O.; Csabai, Z.; Duarte, S.; Geigerg, M.F.; Grabowski, M.; Rimet, F.; Rulik, B.; et al. DNA barcode reference libraries for the monitoring of aquatic biota in Europe: Gap-analysis and recommendations for future work. Sci. Total Environ. 2019, 678, 499–524. [Google Scholar] [CrossRef] [PubMed]
  46. Specchia, V.; Tzafesta, E.; Marini, G.; Scarcella, S.; D’Attis, S.; Pinna, M. Gap Analysis for DNA Barcode Reference Libraries for Aquatic Macroinvertebrate Species in the Apulia Region (Southeast of Italy). J. Mar. Sci. Eng. 2020, 8, 538. [Google Scholar] [CrossRef]
  47. Kelly, R.P.; Closek, C.J.; O’Donnell, J.L.; Kralj, J.E.; Shelton, A.O.; Samhouri, J.F. Genetic and manual survey methods yield different and complementary views of an ecosystem. Front. Mar. Sci. 2017, 3, 283. [Google Scholar] [CrossRef] [Green Version]
  48. Leite, B.; Vieira, P.; Teixeira, M.; Lobo-Arteaga, J.; Hollatz, C.; Borges, L.; Duarte, S.; Troncoso, J.; Costa, F. Gap-analysis and annotated reference library for supporting macroinvertebrate metabarcoding in Atlantic Iberia. Reg. Stud. Mar. Sci. 2020, 36, 101307. [Google Scholar] [CrossRef]
  49. Elbrecht, V.; Leese, F. PrimerMiner: An r package for development and in silico validation of DNA metabarcoding primers. Methods Ecol. Evol. 2016, 8, 622–626. [Google Scholar] [CrossRef]
  50. Kozlov, A.M.; Darriba, D.; Flouri, T.; Morel, B.; Stamatakis, A. RAxML-NG: A fast, scalable, and user-friendly tool for maximum likelihood phylogenetic inference. Bioinformatics 2019, 35, 4453–4455. [Google Scholar] [CrossRef] [Green Version]
  51. Zhang, J.; Kapli, P.; Pavlidis, P.; Stamatakis, A. A general species delimitation method with applications to phylogenetic placements. Bioinformatics 2013, 29, 2869–2876. [Google Scholar] [CrossRef] [Green Version]
  52. White, B.P.; Pilgrim, E.M.; Boykin, L.M.; Stein, E.D.; Mazor, R.D. Comparison of four species-delimitation methods applied to a DNA barcode data set of insect larvae for use in routine bioassessment. Freshw. Sci. 2014, 33, 338–348. [Google Scholar] [CrossRef]
  53. Saitou, N.; Nei, M. The neighbor-joining method: A new method for reconstructing phylogenetic trees. Mol. Biol. Evol. 1987, 4, 406–425. [Google Scholar] [CrossRef] [PubMed]
  54. Tamura, K.; Nei, M.; Kumar, S. Prospects for inferring very large phylogenies by using the neighbor-joining method. Proc. Natl. Acad. Sci. USA 2004, 101, 11030–11035. [Google Scholar] [CrossRef] [Green Version]
  55. Zharkikh, A.; Li, W.H. Statistical properties of bootstrap estimation of phylogenetic variability from nucleotide sequences: II. Four taxa without a molecular clock. J. Mol. Evol. 1992, 35, 356–366. [Google Scholar] [CrossRef]
  56. Cerca, J.; Meyer, C.; Purschke, G.; Struck, T.H. Delimitation of cryptic species drastically reduces the geographical ranges of marine interstitial ghost-worms (Stygocapitella; Annelida, Sedentaria). Mol. Phylogenet. Evol. 2019, 143, 106663. [Google Scholar] [CrossRef]
  57. Jungcurt, T.; Bocksberger, G.; Brüggemann, W.; Berberich, T. Phylogenetic relationship and molecular taxonomy of African grasses of the genus Panicum inferred from four chloroplast DNA-barcodes and nuclear gene sequences. J. Plant Res. 2012, 126, 363–371. [Google Scholar] [CrossRef]
  58. Lee, T.R.C.; Anderson, S.J.; Tran-Nguyen, L.T.T.; Sallam, N.; Le Ru, B.P.; Conlong, D.; Powell, K.; Ward, A.; Mitchell, A. Towards a global DNA barcode reference library for quarantine identifications of lepidopteran stemborers, with an emphasis on sugarcane pests. Sci. Rep. 2019, 9, 7039. [Google Scholar] [CrossRef]
  59. Pinto, M.; Ojala-Barbour, R.; Brito, J.; Menchaca, A.; Carvalho, A.L.; Weksler, M.; Amato, G.; Lee, T. Rodents of the eastern and western slopes of the Tropical Andes: Phylogenetic and taxonomic insights using DNA barcodes. Therya 2018, 9, 15–27. [Google Scholar] [CrossRef] [Green Version]
  60. Peña, V.; De Clerck, O.; Afonso-Carrillo, J.; Ballesteros, E.; Bárbara, I.; Barreiro, R.; Le Gall, L. An integrative systematic approach to species diversity and distribution in the genusMesophyllum(Corallinales, Rhodophyta) in Atlantic and Mediterranean Europe. Eur. J. Phycol. 2014, 50, 20–36. [Google Scholar] [CrossRef]
  61. Silva, F.; Wiedenbrug, S. Integrating DNA barcodes and morphology for species delimitation in the Corynoneura group (Diptera: Chironomidae: Orthocladiinae). Bull. Entomol. Res. 2014, 104, 65–78. [Google Scholar] [CrossRef] [PubMed]
  62. Liu, J.; Jiang, J.; Song, S.; Tornabene, L.; Chabarria, R.; Naylor, G.J.P.; Li, C. Multilocus DNA barcoding–Species Identification with Multilocus Data. Sci. Rep. 2017, 7, 16601. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Biology 10 01092 g001 550
Figure 1. Map of the transitional aquatic ecosystems sampled in the Apulia region (South-East Italy). The numbers on the map correspond to the following ecosystems and sampling sites of each transitional water ecosystem: Laguna di Lesina-1a, 1b, 1c; Lago Varano-2a, 2b, 2c; Vasche Evaporanti, Lago Salpi-3; Torre Guaceto-4; Punta della Contessa-5; Cesine-6; Alimini Grande-7a, 7b; Baia di Porto Cesareo-8; Mar Piccolo-Primo Seno-9; Mar Piccolo-Secondo Seno-10.
Figure 1. Map of the transitional aquatic ecosystems sampled in the Apulia region (South-East Italy). The numbers on the map correspond to the following ecosystems and sampling sites of each transitional water ecosystem: Laguna di Lesina-1a, 1b, 1c; Lago Varano-2a, 2b, 2c; Vasche Evaporanti, Lago Salpi-3; Torre Guaceto-4; Punta della Contessa-5; Cesine-6; Alimini Grande-7a, 7b; Baia di Porto Cesareo-8; Mar Piccolo-Primo Seno-9; Mar Piccolo-Secondo Seno-10.
Biology 10 01092 g001
Biology 10 01092 g002a 550Biology 10 01092 g002b 550
Figure 2. Variation of species richness in the sampling sites of the Apulia transitional waters where the benthic macroinvertebrates were collected during the fall of 2010 and the spring of 2011. Vertical bars indicate the standard error.
Figure 2. Variation of species richness in the sampling sites of the Apulia transitional waters where the benthic macroinvertebrates were collected during the fall of 2010 and the spring of 2011. Vertical bars indicate the standard error.
Biology 10 01092 g002aBiology 10 01092 g002b
Biology 10 01092 g003 550
Figure 3. Correlation analysis between the species richness values calculated in the fall and spring sampling periods from both morphological and reduced database.
Figure 3. Correlation analysis between the species richness values calculated in the fall and spring sampling periods from both morphological and reduced database.
Biology 10 01092 g003
Biology 10 01092 g004 550
Figure 4. Variation of Shannon index in the sampling sites of the Apulia region transitional waters where the benthic macroinvertebrates were collected during the fall of 2010 and the spring of 2011. Vertical bars indicate the standard error.
Figure 4. Variation of Shannon index in the sampling sites of the Apulia region transitional waters where the benthic macroinvertebrates were collected during the fall of 2010 and the spring of 2011. Vertical bars indicate the standard error.
Biology 10 01092 g004
Biology 10 01092 g005 550
Figure 5. Correlation analysis between the Shannon index diversity values calculated in the fall and spring sampling periods from the morphological and reduced databases.
Figure 5. Correlation analysis between the Shannon index diversity values calculated in the fall and spring sampling periods from the morphological and reduced databases.
Biology 10 01092 g005
Biology 10 01092 g006a 550Biology 10 01092 g006b 550
Figure 6. Variation of AMBI in the sampling sites of the Apulia region transitional waters where the benthic macroinvertebrates were collected during the fall of 2010 and the spring of 2011. Vertical bars indicate the standard error.
Figure 6. Variation of AMBI in the sampling sites of the Apulia region transitional waters where the benthic macroinvertebrates were collected during the fall of 2010 and the spring of 2011. Vertical bars indicate the standard error.
Biology 10 01092 g006aBiology 10 01092 g006b
Biology 10 01092 g007 550
Figure 7. Correlation analysis between AMBI indicator values calculated in the fall and spring sampling period from the morphological and reduced databases.
Figure 7. Correlation analysis between AMBI indicator values calculated in the fall and spring sampling period from the morphological and reduced databases.
Biology 10 01092 g007
Biology 10 01092 g008 550
Figure 8. Variation of M-AMBI in the in the sampling sites of the Apulia region transitional waters where the benthic macroinvertebrates were collected during the fall of 2010 and the spring of 2011. Vertical bars indicate the standard error.
Figure 8. Variation of M-AMBI in the in the sampling sites of the Apulia region transitional waters where the benthic macroinvertebrates were collected during the fall of 2010 and the spring of 2011. Vertical bars indicate the standard error.
Biology 10 01092 g008
Biology 10 01092 g009 550
Figure 9. Correlation analysis between M-AMBI ecological indicator values calculated in the fall and spring sampling period from the morphological and reduced databases.
Figure 9. Correlation analysis between M-AMBI ecological indicator values calculated in the fall and spring sampling period from the morphological and reduced databases.
Biology 10 01092 g009
Biology 10 01092 g010 550
Figure 10. Phylogenetic analysis of consensus COI sequences available in BOLD Systems. Tree generated under the maximum likelihood criterion using the general time reversible and gamma model in RAxML. The tree is drawn to scale with branch lengths measured in number of substitutions per site.
Figure 10. Phylogenetic analysis of consensus COI sequences available in BOLD Systems. Tree generated under the maximum likelihood criterion using the general time reversible and gamma model in RAxML. The tree is drawn to scale with branch lengths measured in number of substitutions per site.
Biology 10 01092 g010
Table
Table 1. Ecological quality classes of M-AMBI calculated for each lagoon/station of Apulian transitional water ecosystems, as shown in Figure 1, in the fall and spring, using both the morphological and DNA-barcode data.
Table 1. Ecological quality classes of M-AMBI calculated for each lagoon/station of Apulian transitional water ecosystems, as shown in Figure 1, in the fall and spring, using both the morphological and DNA-barcode data.
Lagoon/Sampling SiteFallSpring
Morphological DatabaseReduced DatabaseMorphological DatabaseReduced Database
1aModeratePoorPoorPoor
1bPoorBadModerateModerate
1cModerateModeratePoorPoor
2aPoorPoorGoodGood
2bModeratePoorPoorPoor
2cBadBadGoodModerate
3BadBadBadBad
4PoorPoorBadBad
5BadBadBadBad
6BadBadBadBad
7aBadBadBadBad
7bPoorPoorModerateModerate
8ModerateBadModerateModerate
9PoorModerateModeratePoor
10HighHighGoodModerate
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Pinna, M.; Saccomanno, B.; Marini, G.; Zangaro, F.; Kabayeva, A.; Khalaj, M.; Shaimardan, L.; D’Attis, S.; Tzafesta, E.; Specchia, V. Testing the Influence of Incomplete DNA Barcode Libraries on Ecological Status Assessment of Mediterranean Transitional Waters. Biology 2021, 10, 1092. https://doi.org/10.3390/biology10111092

AMA Style

Pinna M, Saccomanno B, Marini G, Zangaro F, Kabayeva A, Khalaj M, Shaimardan L, D’Attis S, Tzafesta E, Specchia V. Testing the Influence of Incomplete DNA Barcode Libraries on Ecological Status Assessment of Mediterranean Transitional Waters. Biology. 2021; 10(11):1092. https://doi.org/10.3390/biology10111092

Chicago/Turabian Style

Pinna, Maurizio, Benedetta Saccomanno, Gabriele Marini, Francesco Zangaro, Akbota Kabayeva, Mina Khalaj, Laura Shaimardan, Simona D’Attis, Eftychia Tzafesta, and Valeria Specchia. 2021. "Testing the Influence of Incomplete DNA Barcode Libraries on Ecological Status Assessment of Mediterranean Transitional Waters" Biology 10, no. 11: 1092. https://doi.org/10.3390/biology10111092

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop