Applicability of diatom metabarcoding in the ecological status assessment of Hungarian lotic and soda pan habitats

Diatoms, eukaryotic algae with silica cell wall have been proved to be reliable bioindicators and are applied in ecological status assessment of aquatic ecosystems using indices calculated from pollution sensitivity and indicator values of species. Traditional method of identification of diatom species is based on morphological features of frustule that requires in-depth taxonomical knowledge and expertise. Metabarcoding that combines barcodes and high-throughput sequencing offers a promising alternative. In this pilot study we tested the applicability of metabarcoding of benthic diatom assemblages for ecological status assessment of water bodies in Hungary comparing its performance to that of morphology-based identification of species. Samples from 78 rivers and streams arranged along a trophic gradient and 14 soda pans representing unique types of habitats were investigated. Sequences of the barcode region of rbcL gene suggested by Vasselon et al. (2017) were acquired revealing 1135 diatom amplicon sequence variants (ASVs) of which 54% could be assigned at species level using Diat.barcode (Rimet et al. 2019) as reference database. Morphological investigation detected more species and intraspecific taxa in the lotic and lentic samples (413 and 78, respectively) than metabarcoding (190 and 75). Correspondence between taxa detected based on morphology and DNA sequences was relatively low (on average 24% in lotic and 26% in lentic samples) but considerably increased (on average 66% and 56%, respectively) when taking into account only the taxa reaching higher than 5% relative abundance in microscopy counting. The differences mainly derived from that a considerable portion (59% and 32%, respectively) of the morphologically identified taxa were not recorded in the reference database. Community structure of samples from both running and standing waters based on microscopy and DNA sequence analysis showed significant correlation as revealed using Mantel test. For rivers and streams diatom index IPS (Specific Pollution Sensitivity Index, (Coste in Cemagref, 1982) for soda pans indices H (Ziemann et al., 2001) and IBD (Indice Biologique Diatomées, Lenoir and Coste, 1996) were calculated from the results of the two aspects. Strong correlation was found between morphology- and sequence-based indices and in about half of the samples, the ecological status class obtained with the two methods coincided. Our results suggest that metabarcoding inheres a great opportunity and could be successfully applied in benthic diatom-based ecological status assessment of Hungarian waters after the reference database is supplemented with taxa frequently occurring in these habitats. References Cemagref, 1982. Etude des méthodes biologiques quantitativesd’appréciation de la qualité des eaux, Rapport Division Qualité des Eaux Lyon. Agence financière de Bassin Rhône-Méditerranée-Corse, Pierre-Bénite, France. Lenoir, A. , Coste, M. , 1996. Development of a practical diatom index of overall water quality applicable to the French National Water Board Network. Presented at the International symposium, Volksbildungsheim Grilhof Vill, AUT, 17-19 September 1995, Universität Innsbruck, pp. 29–43. Rimet, F. , Gusev, E. , Kahlert, M. , Kelly, M.G. , Kulikovskiy, M. , Maltsev, Y. , Mann, D.G. , Pfannkuchen, M. , Trobajo, R. , Vasselon, V. , Zimmermann, J., Bouchez, A. , 2019. Diat.barcode, an open-access curated barcode library for diatoms. Sci. Rep. 9, 15116. Vasselon, V. , Rimet, F. , Tapolczai, K. , Bouchez, A. , 2017. Assessing ecological status with diatoms DNA metabarcoding: Scaling-up on a WFD monitoring network (Mayotte island, France). Ecol. Indic. 82, 1–12. Ziemann, H. , Kies, L. , Schulz, C. -J., 2001. Desalinization of running waters: III. Changes in the structure of diatom assemblages caused by a decreasing salt load and changing ion spectra in the river Wipper (Thuringia, Germany). Limnologica 31, 257–280.