Development of eDNA methods for monitoring two stygobiotic species of the Dinaric Karst, Proteus anguinus and Congeria jalzici, using digital PCR Razvoj metod eDNA za monitoring dveh stigobiontov Dinarskega krasa, človeške ribice (Proteus anguinus) in Jalžićeve jamske školjke (Congeria jalzici), z uporabo digitalne PCR Špela GORIČKI*1, Primož PRESETNIK2, Uršula PROSENC-ZMRZLJAK3, Tajda GREDAR4, Matej BLATNIK5,6, Blaž KOGOVŠEK5,6, Oliver KOIT7, Cyril MAYAUD5,6, Sara STRAH8, Branko JALŽIĆ9, Gregor ALJANČIČ10, Dejan ŠTEBIH11, Andrej HUDOKLIN12, Rok KOŠIR3 1Scriptorium biologorum - Biološka pisarna d.o.o., Ulica Nikole Tesla 6, 9000 Murska Sobota, Slovenia; E-mail: goricki.spela@gmail.com 2Tolstojeva ulica 9b, SI-1000 Ljubljana, Slovenia; E-mail: primoz.presetnik@amis.net 3BIA Separations CRO, Labena d.o.o, Teslova ulica 30, SI-1000 Ljubljana, Slovenia; E-mails: ursula.prosenc@biaseparationscro.com, rok.kosir@biaseparationscro.com 4Department of Biology, Biotechnical Faculty, University of Ljubljana, Večna pot 111, SI-1000 Ljubljana, Slovenia; E-mail: tajda.gredar@gmail.com 5ZRC SAZU Karst Research Institute, Titov trg 2, SI-6230 Postojna, Slovenia; E-mails: mblatnik@zrc-sazu.si, blaz.kogovsek@zrc-sazu.si, cyril.mayaud@zrc-sazu.si 6UNESCO Chair on Karst Education, University of Nova Gorica, Glavni trg 8, SI-5271 Vipava, Slovenia 7Institute of Ecology, Tallinn University, Narva Rd 25, 10120 Tallinn, Estonia; E-mail: oliver.koit@tlu.ee 8University of Primorska, Titov trg 4, SI-6000 Koper, Slovenia; E-mail: sarastrah97@gmail.com 9Croatian Biospeleological Society, Demetrova 1, HR-10000 Zagreb, Croatia; E-mail: jalzicbranko@gmail.com 10Society for Cave Biology, Oldhamska cesta 8A, SI-4000 Kranj, Slovenia; E-mail: gregor.aljancic@guest.arnes.si 11National Institute of Biology, Večna pot 111, SI-1000 Ljubljana, Slovenia; E-mail: dejan.stebih@nib.si 12Institute of the Republic of Slovenia for Nature Conservation, Novo mesto Regional Unit, Adamičeva ulica 2, SI-8000 Novo mesto, Slovenia; E-mail: andrej.hudoklin@zrsvn.si The welfare of species and ecosystems are assessed by biological field surveys. Without knowing about the existence of particular species, we cannot protect them. Environmental DNA (eDNA) methods detect species at very low densities. When carefully validated and appropriately interpreted (see e.g. Cristescu & Herbert 2018), eDNA-based inventories improve biological surveys of ecosystems and are used to guide conservation efforts (e.g. Zaiko et al. 2018). Apart from establishing the presence of endangered or rare species, eDNA methods are also applied for monitoring the spread of alien invasive species and infectious agents such as fungal or bacterial spores, viruses and parasites. Analyses of total eDNA in a sample (i.e. environmental metagenome) are utilized to obtain information on the species community structure, the relative abundance of different trophic levels, and their interactions. However, a number of properties of eDNA molecules, such as differences in their release rates by different species during different life-history stages, their persistence under different environmental conditions and their variable detection rates using different techniques, are still insufficiently understood. Consequently, while eDNA methods are becoming wildly popular for establishing the presence of species in diverse environments, few studies to date have attempted to use eDNA quantification for determining population sizes and trends in these species (e.g. Doi et al. 2015a, 2017, Buxton et al. 2017, Chambert et al. 2018 and references therein). Recent records of the olm Proteus anguinus outside its historically known range, discovered through detection of its DNA dissolved in groundwater, introduced the eDNA methodology to the study and conservation in the cryptic subterranean environment (Gorički et al. 2017). Since then, we have been able to detect the presence of P. anguinus in karst groundwater by analyzing samples of water collected on the surface – at springs or in wells. An upgraded technology, droplet digital PCR (ddPCR), which is reportedly more sensitive, accurate and resistant to PCR inhibitors than classical qPCR (Doi et al. 2015b, Taylor et al. 2017, Baker et al. 2018), is being tested for direct quantification of eDNA molecules in groundwater. We are using the large, accessible and relatively well-characterized natural P. anguinus population inhabiting the Planinska jama Cave (southwestern Slovenia) as a model. The size of this population was recently estimated using a genetic mark-recapture method, although only the number of animals actually caught was reported (Zakšek & Trontelj 2017). In succession of the named survey, we sampled water at several sites along the subterranean Pivka River and recorded representative stream profiles and flow velocities at or near the sampling sites (Fig. 1). Thereupon we estimated the number of eDNA molecules specific to P. anguinus in three 100 m long sections of the stream (B, F and I of Zakšek & Trontelj 2017). Based on the data reported by Zakšek & Trontelj (2017), we also estimated the number of eDNA molecules per individual animal in the selected stream sections. Until now we have tested several sampling and filtration strategies and found that their success greatly depends on the amount of suspended particles in the water, which in turn reflects wet weather conditions during a yet undetermined period of time before sampling. On the other hand, during more stable, i.e. better hydrological conditions for the eDNA method, we were able to obtain comparable molecule counts in two consecutive years in all three river sections that were sampled twice. The relative number of eDNA copies per individual animal was the lowest in section F, where the highest number of animals had been observed, but where the slowest water flow was measured. This may indicate a higher rate of eDNA sedimentation (sinking), although we failed to detect more eDNA copies closer to the stream bottom than on the water surface. Alternatively, our result might reflect differences in eDNA release rates between animals from different sections. This explanation is even more plausible when we compare the other two sections with greater hydrologic similarity, yet we consistently detected more eDNA copies in section B than in I – both in absolute and relative terms – which apparently cannot be explained by transfer from the upstream sections alone. Our results indicate the need to investigate and characterize the Planinska jama Cave and its resident P. anguinus population in greater detail, as well as to measure the many variables in the eDNA model in a more controlled, experimental setting. Figure 1. Stream profiling of the Pivka channel in the Planina Cave (photo: C. Mayaud). Slika 1. Meritev pretoka na Pivškem rokavu Planinske jame (foto: C. Mayaud). In another line of eDNA research, the utility of ddPCR is being explored for improved detection of the much smaller and rare stygobiont, the cave clam Congeria jalzici. The clam is known from only one site in Slovenia and three sites in Croatia (Hudoklin & Ilenič 2012, Bilandžija et al. 2013). The spring of the Krupa River in southeastern Slovenia (Fig. 2) is another presumed site of C. jalzici, although only shells, but no living animals were found in it (Sket 1971 in Hudoklin & Ilenič 2012). This spring is the site of one of the greatest ecological disasters in the Dinaric Karst. Between 1962 and 1984, dozens of tons of pure polychlorinated biphenyls (PCBs) were dumped at waste disposal sites and in nearby dolines by a local condenser factory. The undegradable carcinogenic chemicals are still gradually released from the sediment into the groundwater and were detected in the tissues of P. anguinus residing in the aquifer (Pezdirc et al. 2011). We detected eDNA of C. jalzici in the Krupa spring water and thereby confirmed the presence of a living population there. In the future, we will survey additional sites in the area potentially inhabited by the mollusk to assess the presence and vulnerability of any newly discovered populations, and to monitor the spread of the alien invasive zebra mussel Dreissena polymorpha, which may soon become a new major threat to the Slovenian populations of C. jalzici. In conclusion, eDNA holds great potential for monitoring and conservation of fauna in inaccessible subterranean habitats. In the future, the eDNA methodology might be applied in the estimation of P. anguinus population sizes without having to see, mark or otherwise disturb the animals themselves. In parallel to eDNA assay development for various stygobiotic species of the Dinaric Karst, a groundwater-sample library is being created. This collection of samples will be available for future analysis of their species composition, once a reference sequence database for known Dinaric Karst species is created. It will serve as a source of information on any changes in species distribution over time, or for some species, as a record of their existence before they are lost forever. Figure 2. Water sampling at the Krupa spring (photo: P. Presetnik). Slika 2. Vzorčenje vode na izviru Krupe (foto: P. Presetnik). Acknowledgements We are grateful to Helena Bilandžija (Ruđer Bošković Institute, Zagreb, Croatia), Tjaša Lokovšek (ZRC SAZU Jovan Hadži Institute of Biology, Ljubljana, Slovenia), Franci Gabrovšek (ZRC SAZU Karst Research Institute, Postojna, Slovenia), Matjaž Kuntner (ZRC SAZU Jovan Hadži Institute of Biology, Ljubljana, Slovenia), Peter Trontelj (Department of Biology, University of Ljubljana, Slovenia), Marijan Govedič (Centre for Cartography of Fauna and Flora, Ljubljana, Slovenia), Barbara Kink (ZRSVN-IRSNC, Novo Mesto Regional Unit, Slovenia), William Jeffery, Jasmina Kotnik, Eva Pavlovič, Rudi Kraševec, Damjan Vinko, Klemen Kramar, Petra Kovač-Konrad, Jenny Barnjak, Larry Cohen and everyone who donated to the project Through a glass darkly: assessing population size of an endangered cave salamander from samples of spring and cave water (https://experiment.com/cavesalamander). The 2018 ddPCR part of the analysis was performed in collaboration and under sponsorship of Labena d.o.o. and their ddPCR Grant Challenge. References Baker C.S., Steel D., Nieukirk S., Klinck H. (2018): Environmental DNA (eDNA) from the wake of the whales: droplet digital PCR for detection and species identification. Front. Mar. Sci. 5: 133. Bilandžija H., Morton B., Podnar M., Ćetković H. (2013): Evolutionary history of relict Congeria (Bivalvia: Dreissenidae): unearthing the subterranean biodiversity of the Dinaric Karst. Front. Zool. 10: 5. Buxton A.S., Groombridge J.J., Zakaria N.B., Griffiths R.A. (2017): Seasonal variation in environmental DNA in relation to population size and environmental factors. Sci. Rep. 7: 46294. Chambert T., Pilliod, D.S., Goldberg C.S., Doi H., Takahara T. (2018): An analytical framework for estimating aquatic species density from environmental DNA. Ecol. Evol. 8: 3468-3477. Cristescu M.E., Hebert P.D.N. (2018): Uses and misuses of environmental DNA in biodiversity science and conservation. Annu. Rev. Ecol. Evol. Syst. 49: 209-230. Doi H., Uchii K., Takahara T., Matsuhashi S., Yamanaka H., Minamoto T. (2015a): Use of droplet digital PCR for estimation of fish abundance and biomass in environmental DNA surveys. PLoS ONE 10(3): e0122763. Doi H., Takahara T., Minamoto T., Matsuhashi S., Uchii K., Yamanaka H. (2015b): Droplet digital polymerase chain reaction (PCR) outperforms real-time PCR in the detection of environmental DNA from an invasive fish species. Environ. Sci. Technol. 49(9): 5601-5608. Doi H., Inui R., Akamatsu Y., Kanno K., Yamanaka H., Takahara T., Minamoto T. (2017): Environmental DNA analysis for estimating the abundance and biomass of stream .sh. Freshw. Biol. 62: 30-39. Gorički Š., Stanković D., Snoj A., Kuntner M., Jeffery W.R., Trontelj P., Pavićević M., Grizelj Z., Năpăruş-Aljančič M., Aljančič G. (2017): Environmental DNA in subterranean biology: range extension and taxonomic implications for Proteus. Sci. Rep. 7: 45054. Hudoklin A., Ilenič, T. (2012): Izvir jamske školjke. Dolenjski kras 6: 139-142. Pezdirc M., Heath E., Bizjak Mali L., Bulog B. (2011): PCB accumulation and tissue distribution in cave salamander (Proteus anguinus anguinus, Amphibia, Urodela) in the polluted karstic hinterland of the Krupa River, Slovenia. Chemosphere 84(7): 987-993. Taylor S.C., Laperriere G., Germain H. (2017): Droplet digital PCR versus qPCR for gene expression analysis with low abundant targets: from variable nonsense to publication quality data. Sci. Rep. 7: 2409. Zaiko A., Pochon X., Garcia-Vasquez E., Olenin S., 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. Zakšek V., Trontelj P. (2017): Conservation genetics of proteus in the Postojna-Planina Cave System. Nat. Slo. 19(1): 33-34.