Abstract
Dytiscus latissimus (Coleoptera Dytiscidae) is an endangered diving beetle throughout its range. It is one of the two species of Dytiscidae listed in Annex II of the Habitats Directive, IUCN red list and in many national level legislations and therefore strictly protected. The conservation of endangered species first of all requires an assessment of their population size. Until now, a method has not been developed for estimating the size of D. latissimus populations. The article summarizes the results of two studies carried out independently in Germany and Latvia. Both studies were carried out in one water body used recapture method but with a different spatial placement of traps, which, according to our data, is an important factor in population estimation. We evaluated Jolly-Seber and Schnabel approaches of estimating aquatic beetle's populations and found that confidence intervals obtained by different methods in our research do not differ significantly, but combination of both models provide the most accurate estimates of population dynamics. As part of the study, we concluded that the populations of Dytiscus latissimus are relatively closed, so we accept that the Schnabel estimate shows more accurate data. By fixing the places of capture of each individual, it was found that females live mainly locally, and males actively move within the water body. This aspect indicates the advantage of the spatial placement of traps compared to the use of transects. The results of our study show a significantly higher number of both captured and recaptured males Such a sex ratio may indicate both a greater activity of males and differences in the sex ratio in the population. The study confirmed that environmental changes, such as the water level in a water body, can also significantly affect the result of a population assessment. In the frame of D. latissimus monitoring, to obtain an objective estimation of the species population size we recommend using four traps for each 100 m of water body shoreline with 4–8 censuses, dependently on the recapture rate.
Similar content being viewed by others
Introduction
Dytiscidae is the most diverse beetle family that inhabits freshwater environments. The family contains more than 4200 known species that are divided in 10 or 11 tribes with predicted total species richness of around 5400 species1,2,3,4. Species of the Dytiscidae family inhabit a broad range of habitats from large lakes to seasonally drying ponds, peatlands and salt marshes, mountain and lowland ponds which can be oligotrophic to eutrophic5,6,7. Despite the fact that there is an increase in taxonomic, systematic and phylogeographic studies dytiscids are still an understudied group in terms of their ecology and biology8. One of such understudied species that attracts more attention due to its inclusion in several legislative acts in Europe is Dytiscus latissimus, with 44 mm body length one of the largest diving beetle species in the world (Figure B). D. latissimus is western Palearctic species with relatively large distribution range, occurs in most central European countries from France and Italy in the south, Sweden and Finland in the north, ranging eastwards to Western Siberia9. It inhabits a wide variety of freshwater environments such as shorelines of dystrophic and mesotrophic lakes that are densely vegetated by Carex and Equisetum6 or surrounded by flooded Phragmites belts10,11. In the south of its distribution, D. latissimus prefers to inhabit more oligotrophic lakes and in north it occurs mainly in the most nutrient-rich lakes12. Three of the most important factors necessary for D. latissimus are (1) the permanency of the water body that covers at least 1 ha and is at least 1 m deep (2) the specific plants along the shoreline necessary for oviposition (Caltha palustris, Carex acuta, Carex pseudocyperus, Carex rostrata, Menyanthes trifoliata) and (3) sufficient food supply for larvae, which is mainly formed by caddisflies Limnephilidae (Trichoptera) larvae6,12,13,14.
The populations of D. latissimus have been declining over most of its range since the early twentieth century, especially in the western part, while Latvia has the largest number of confirmed localities in the Baltic states6,15. In Germany D. latissimus population trends indicate a rapid decline since the 1960s16. At the beginning of the twentieth century the species occurred in most of the territories in Germany but after 2000 it has been observed only in five out of 16 federal states of the country, with stable populations only in the federal states of Brandenburg, Mecklenburg-Vorpommern and Bavaria10,11,17,18.
D. latissimus is included in the European Habitats Directive 82/43/EEC19 as well as in the Bern Convention Annex II, and it is also protected by legislation in many EU member states. The putative causes for the population decline are understudied and include drainage of lakes and agricultural land, agricultural intensification which contaminates freshwater systems with fertilizers and pesticides, urbanization, industrialization, deforestation, climate change and invasions of alien species20 or by increasing the abundance of more competitive species that prefer higher mean annual temperatures (e.g. the diving beetle Cybister lateralimarginalis)11. These factors may have an indirect effect on D. latissimus, contributing to the decline of the caddisfly population and thus the beetle larvae lack prey21.
Most of the research on D. latissimus is about distribution areas in different European countries, which has provided the basic knowledge about preferred habitats and possible factors of decline22,23,24,25,26. Studies on feeding preferences of D. latissimus suggest that during larval development they prefer Limnephilidae larvae, lack of which may influence the habitability of a water body13,21. Majority of the available data on the species currently is simple reports on presence/absence. Very little is known research on the population dynamics, rates of dispersion and spatial distribution within natural habitats has been done27. Given that estimating population sizes is a key in animal ecology, evolution, and conservation biology28,29,30,31 and that the assessment of the D. latissimus populations must be carried out according to European Union legislation, it is essential to develop an effective method. Due to a hidden lifestyle it is impossible to count all individuals in the population thus, the Capture-Mark-Recapture (CMR) method has been widely used by field biologists and ecologists to investigate the dynamics of biological populations27,32,33,34,35. There are various statistical models for analysing CMR data35,36,37,38, and methodologies for zoological field work38,39,40. To develop effective methods for specific groups of organisms, it is necessary to conduct studies that evaluate various approaches and consider the ecological characteristics of the target species, resulting in the most reliable data on the population size.
This is the first comprehensive study of the population biology of the critically endangered diving beetle D. latissimus, in which population size, density, age distribution and life expectancy were assessed. To evaluate the CMR approach in assessing the D. latissimus population size and the patterns of spatial distribution, we compared two independent extensive field surveys in different parts of its distribution range in the Northeast of Germany, in the nature reserve “Rothes Moor bei Wesenberg”41 and in the East of Latvia, in the nature protected territory "Rāznas Nacionālais Parks". In the analysis of the collected CMR data, we also compared two approaches, the Schnabel42 and Jolly Seber method43,44. Both studies used the same recapture method but differed in the spatial placement of traps. The aims of this study are to evaluate different approaches to estimating aquatic beetle populations using recapture method and describe movement patterns. Moreover, we provide recommendations for future monitoring programs.
Methods
Study area
The study area included two water bodies in which traps were set: The Lake Rothemoorsee (13.0253275°E 53.2774556°N) situated in the Northeast of Germany, in the nature reserve “Rothes Moor bei Wesenberg” and Lake Glušonoks (27,3,599,552°E 56,3,814,767°N) situated in the East of Latvia, in the nature protected area "Rāznas Nacionālais Parks" (Fig. 1).
Lake Rothemoorsee (Figs. 2a, 3a) is a slightly dystrophic to mesotrophic lake deficient in lime, surrounded by large bogs and marshes of intermediate type with former peat layer in the bogs. Several species of sedges (Carex.spp.) are common in the shallow water zone of the lake; most parts of the lakeside strip are overgrown with dense stands of the common reed Phragmites australis. Near the lakeside the white lily (Nymphaea alba) is growing inside the stands of the common reed. In more open water the yellow water lily (Nuphar lutea) covers most part of the water surface. The lake is inhabited by predatory fish. Pike and perch were regularly caught in traps during the survey. Special studies of the ichthyofauna in the lake were not conducted, therefore, determining the population density of predatory species is not possible. according to the authors of the study, it can be considered that the population density of predatory fish in this lake is quite low. The Lake Glušonoks corresponds to the habitat included in the first Annex of the Habitats Directive: Natural eutrophic lakes with Magnopotamion or Hydrocharition type vegetation19. It is a brown-water (poly-humic) lake with diverse vegetation. Several species of sedges (Carex spp.) are common in the coastal zone of the lake; part of the coastal strip is overgrown with common reeds. Near the shore, the water soldiers (Stratiotes aloides) and the yellow water lily constitute most of the surface coverage. No predatory fish were found in the lake during the study period, so the presence of predatory fish during this period is unlikely. However, in the course of later research, one specimen of pike was found in the lake. We assume that it got into the lake after the research, but the existence of a very small population of pike in the lake is not excluded.
The average depth of Lake Rothemoorsee is 0.4 to 1.5 m and the Glušonoks is 1 to 1.5 m., The Lake Rothemoorsee has an area approximately twice the size of the Glušonoks 5.5 ha and 2.4 ha respectively.
Trap construction and placement in Rothemoorsee
In the Lake Rothemoorsee nine line transects were selected along shallow shores with dense macrophyte growths (Fig. 2a). The transects were calibrated and marked on site. The traps within the transects were placed 3 to 5 m apart. Due to fluctuations in the water level, the number of traps varied between the capture sessions and study years, the exact locations of the traps were slightly adjusted. The maximum number of traps used on the nine transects was 51, but minimum 21; the traps were removed in November. Based on the assessment of lakeside the water level was measured in ordinal scale with five categories: very low (water level significantly below the edge of the shore), low (water level slightly below the edge of the shore), moderate (water level in the same line with the edge of the shore), high (water level slightly above the edge of the shore) and very high (water level significantly above the edge of the shore).
Large (55 × 23 × 23 cm) traps, made of plastic mesh and a metal frame, were baited with bloody pig liver (Fig. 3d). Traps had self-made plastic funnel entrances at both ends with a funnel opening of 30 mm. Two floats were attached along the sides of the traps so that the top was kept above the water level and the captured individuals could breathe. The traps were removed 48 h after being set up. In the shallow water zone, they were located at the bottom, further and deeper in the lake they were floating partly beneath the water surface.
Trap construction and placement in Lake Glušonoks
The surface of Lake Glušonoks (Fig. 4a) was divided into 25 × 25 m quadrants (Fig. 2b). Two paired traps were set in the centre of each quadrant one slightly floating on top of the surface and one on the lakebed. Considering that both traps formed a single trapping unit we did not distinguish whether beetle was caught in the bottom or top trap in the analysis. We used cylindrical baitfish traps covered with netting. The square lattices of the net measured 1 × 1 cm. Trap height was 55 cm, diameter 32 cm. Each trap had circular entrance with 12 cm diameter. Sealed foam tubes were attached to the top traps so they would float and have oxygen for the beetles (Fig. 4d). Traps set on the lakebed were covered with waterproof fabric which was filled with atmospheric air when being submerged (Fig. 4b). In order to keep the bottom traps on the lakebed we attached two 2 kg weights to them. Both traps were connected by a rope at both ends of the cylinder which in order to anchor the trap near the water surface from drifting outside of the quadrant. We used fresh cattle liver as bait. Traps were removed 48 h after being set up.
Timing of data collection
Data acquisition in the Lake Rothemoorsee was carried out in three periods, during which the traps were set up once a month: May–November 2009 (eight capture sessions, 20–39 traps per session), March–November 2010 (10 sessions, 41–51 traps per session), April–November 2011 (nine sessions, 38–46 traps per session, in November only 21 traps). Field surveys in Lake Glušonoks were done from April 2018 to November 2018 with eight surveying trips with approximately one-month intervals. In total 71 trap systems were used for each survey.
Capture-Mark-Recapture (CMR) technique
The CMR methodology in both studies followed similar principles, assuring that so the implementation did not differ significantly. Each captured D. latissimus was marked with a unique, permanent number directly as the traps were emptied. The collection site, sex, time date, and notes, which may include additional information about the age and morphological characteristics of the specimen were documented for each marked individual. The specimens were released back into the same place where they were captured. The slight differences in the CMR methodology between the countries are described below.
In the Lake Rothemoorsee, the marking took place by engraving the pronotum (Fig. 3c), as tested by Mölle45 on a large number of Dytiscus marginalis. The engraving was performed with great caution with the help of a dental burr without damaging the animal. The marking was a continuous individual number for each year for each sex. The marking year was denoted by a special symbol. Freshly metamorphosed specimens with soft exoskeletons were not marked because the risk of injury was too high.
In the Lake Glušonoks Beetles were marked by fixing labels to the elytra with cyanoacrylate “Superglue”, following a modified method described by Davy-Bowker46. After removing the trap, each unlabelled specimen was marked. Each individual was then placed in a separate clean container without water. Shortly afterwards the specimen was dried using a paper towel and the surface of the prothorax was cleaned using a paper towel soaked in 96% ethanol. We used commercially available labels used to mark bee queens. These are made of plastic and are approximately 2 mm in diameter. Each contains three pieces of information: a unique two-digit number, colour of background and colour of the numbers. Once ethanol evaporated from the prothorax surface a label was prepared by applying a small amount of cyanoacrylate glue to the bottom. Label was applied on the pronotum of the beetle (Fig. 4c) and care was taken to avoid any of the glue getting on the sutures. After about five minutes of letting the glue dry, the specimen was placed back into a plastic container for further polymerization of the glue for additional 10 min. When the glue was completely polymerized the specimen was released in the same place where it was collected.
Statistical analysis
Firstly, descriptive statistics were calculated: (1) total number of trap checks, beetles per trap, captures, marked and first captured, and recaptured of marked specimens every year (Table 1); (2) the centre point (median) and spread of captured and recaptured individuals (interquartile range) (Table 2); (3) the Schnabel estimate for closed populations (Table 2) and the Jolly-Seber model for open populations with a 95% confidence intervals (Fig. 5).
Secondly, to test statistical hypotheses the inferential statistics were calculated. To compare the number of captured and recaptured individuals between water levels the Kruskal – Wallis H test with Bonferroni post hoc correction was conducted for Germany (Fig. 6).
To determine the strength and direction of the association between water level and captured and recaptured individuals the Spearman's correlation was calculated. To develop recommendations for the minimum required number of traps for monitoring, we took advantage of the shortcoming of this study—the different number of traps. We used the confidence interval of population size as an indicator of the number of traps in order to find the optimal ratio of traps and the population size indicator, i.e., point where the confidence interval has the smallest size at the minimum number of traps.
As this observation study is replicated across seasons introducing random variable, the Generalized Linear Mixed Model (GLMM) was performed to identify factors influencing Dytiscus occurrence in natural habitat. Because the response variable (captured Dytiscus individuals) comes from the count distribution with overdispersion, the negative binomial distribution with the logarithmic link was applied. The sex ratio and the recapture individuals were used as explanatory variables. All possible models including the null model and interactions were calculated. To select the best model, the Akaike Information Criterion (AIC) was used. The GLMM was carried out separately for the countries because the number of traps in Latvia was the same throughout the study, but in Germany it varies between years. The number of traps was modelled as an offset. Statistical data analysis was conducted using the packages ‘fishmethods’47, marked, lmer and RMark implemented in R 4.0.4.48.
Results
During three field seasons, 1193 specimens were captured in the Lake Rothemoorsee, and 239 were captured in the Lake Glušonoks during one season. Of all marked specimens, 181 were recaptured in the Lake Rothemoorsee and 63 in the Lake Glušonoks. Details of number of trap checks, beetles per trap, capture recapture events are appended in Table 1. Population estimates using the Schnabel approach are summarized in the Table 2, and using the Jolly-Seber (Fig. 5).
The results of GLMM showed that null model better explain captured individuals indicating no effect of explanatory variables (Supplementary Table S1). Significant fluctuations in lake water levels were observed during the study in the Lake Rothemoorsee. Spearman's correlation test showed a moderately significant negative association between the number of captured individuals and the water level in the lake (rs = − 0.486, n = 25, p = 0.014), indicating that the decline in the number of captured beetles could be correlated with the increased level of water in the lake; there was, however, no association between recaptured individuals and the level of water in the lake (p > 0.05). The Kruskal – Wallis H test showed significant difference of captured (X2(2) = 8.08, p = 0.018) and recaptured (X2(2) = 7.96, p = 0.019) individuals between water levels, where post hoc test revealed significant difference only between moderate and high water level (p < 0.05) (Fig. 6).
Based on data from recaptured individuals, the study estimated the distance between the points of the capture of individuals. The Mann -Whitney U test showed that the mean distance covered by recaptured specimens differed significantly between Latvia and Germany U = 1956.5, Z = − 3.754, p < 0.001, r = 0.354 (medium effect size) (Table 3).
The multi-year trapping in the Lake Rothemoorsee allowed us to make inferences of the maximum lifespan of D. latissimus in the wild. One male (marked in October 2008) was found in April 2011. Therefore, it is reported here for the first time that the adults of D. latissimus can reach an age up to three years, – even in the wild.
Discussion
We compared different approaches in estimating D. latissimus populations in two lakes in Latvia and Germany. Both studies were based on two main methods: trapping beetles with modified fish traps and the CMR method. These methods had been used in similar studies, on other organisms, however not on D.latissimus. Obtaining an accurate quantification of population size is often of prime importance in ecology and conservation biology49,50. CMR techniques have been widely used and specialized to study organisms throughout the field of biology51. This method is most suitable for obtaining reliable information on the number of individuals in populations of animal species with a hidden lifestyle. It has for example been successfully applied to estimate population size in the hermit beetle Osmoderma eremita52. In the present study, two approaches were used for CMR data analysis: the Schnabel method based on closed population model and the Jolly-Seber method based on open population model43,44,53. Combination of both models provide the most accurate estimates of population dynamics54.
The Schnabel method has been recognized as the most accurate in estimating freshwater crab54, and the Freshwater Anomura Aegla longirostri55 populations. The method provided results with relatively narrow confidence intervals, expectedly reflecting the real size of Orthoptera, Diptera and Hymenoptera populations38. The Jolly-Seber method is less efficient in assessing insect populations with potentially low recapture rate possibility, due to short-living stage of imago, lacking territorial behaviour and relatively high mobility38. Comparing the confidence intervals obtained by different methods in our research, the results do not differ significantly. This indicates that the populations of Dytiscus latissimus are relatively closed, most of the individuals lives in a limited space, and the proven life span of individuals reaches three years, which makes it possible to effectively use the CMR method. The significant difference lies in the population size estimates, Jolly-Seber estimator showing significantly lower population size results. Similar results, where the Jolly-Seber estimator showed significantly lower population size, were obtained when evaluating the population of Psophus stridulus (Orthoptera)38. It should be noted that this result is related to the low recapture rates, which is not true in our case. At the same time, for example, in Lake Glušonoks, based on the Jolly-Seber estimator data, more than 70% of the population was captured, which is objectively too optimistic. Based on these findings, we accept that the Schnabel estimate shows more accurate data, however, there are several additional factors to consider, which are presented in the next part of the discussion.
Methods of CMR analysis work particularly well for univoltine insects with clear spatial and temporal population boundaries, characterizing many habitat specialists. Dytiscus latissimus is a suitable species for use of the CMR method. In calculating the population size of a diving beetle species, a range of factors can influence the result. The first is that most dytiscid species disperse from one waterbody to another, primarily using aerial flight12. This process may be related to oogenesis flight syndrome56, characterizing female insect that migrate from their birthplaces before reproduction. Beetles’ migration can also be caused by changing conditions in the waterbody. Most previous studies were based on detecting flying insects, using light or flight traps57,58. However, using the light trap method it is impossible to draw conclusions about the rhythm of dispersal flight; thus, the CMR is the more effective method to determine flight intensity of aquatic insects59. For this purpose, the CMR method has been used on diving beetles by Davy-Bowker (2002), but that study did not estimate, population size60. Considering that no studies have been conducted to evaluate the intensity of dispersion of D. latissimus specimens, it is impossible to determine the effect of this factor on the results of population assessment. Another factor that can affect the estimate of the population size is life expectancy. Based on our results, D. latissimus can reach an age of up to three years. Accordingly, CMR allows obtaining additional estimates of mortality, life expectancy and maximum lifespan of individuals of the target species, which are key parameters characterising the dynamics of wildlife populations61. In case of D. latissimus several generations of beetles are included in the population estimation process.
According to the assessment of the abundance of D. latissimus populations in the Lake Rothemooresee and Lake Glušonoks, the number of individuals varies significantly. These differences are plausible, and can be explained by the facts that the Lake Glušonoks is just half of the size than Rothemooresee in Germany. Moreover, the shore regions of the the Lake Rothemoorsee are more richly structured and offer better protection for the larvae against predators (fishes) with their broad Phragmites australis reed beds that are flooded all year round. The analysis of various factors influencing the population size was not carried out within the framework of this study, which is laborious and impracticable in the framework of monitoring activities. The main difference between the methods applied in the Lake Rothemoorsee and the Lake Glušonoks was the placement of traps in the waterbody. In the Lake Rothemoorsee, the transect method was used and the traps covered only parts of the lake´s coastline, while in the Lake Glušonoks the grid method was used and the traps were evenly distributed throughout the lake. In both studies, the traps were placed in fixed locations, but given that the German study was conducted over several years, the results were influenced by the water level in the lake. For this reason, at different stages of the study, the number of traps differed.
Statistical methods were used to compare CMR data from both studies. A measure of spread (such as min–max, IQR, standard deviation etc.) shows how well measures of central tendency (such as mean, median or mode) represents the data62. The narrow range of IQR and min–max measures in the Latvian data (Table 1) indicates that the mean value describes the data better, and that the estimate of population size is more robust here than the corresponding results from the Lake Rothemoorsee. The range of 95% confidence interval of the Schnabel estimate from the Lake Rothemoorsee shows great variation among years of investigation. This suggests either a large effect of water level on population size, or a problem in the methodological approach. Special attention should be given to the results of 2009 where the range of 95% confidence interval was very high, indicating low predictive power of the Schnabel estimate. Finally, the Schnabel estimates for 2009 and 2011 differ by factor of 2, raising a question about the importance of the number and spatial arrangement of traps in the application of the CMR method. Certainly, the influence of water level may affect the population size, but the slight variations observed in Lake Rothemoorsee cannot be considered as the main factor affecting the population size of D. latissimus; therefore, we concluded that the changes in the number of individuals in different study periods are related to the spatial placement of traps.
Great effort has been done in the identification of ecological and methodological factors affecting detectability and estimation of the population size of various animal groups63,64,65,66,67. Since our work is pioneering the population estimates of D. latissimus many possibly confounding factors are yet to be understood. For instance, it still remains unknowns, since it remains unknown how effective the traps used in the study are and whether all the beetles attracted by the bait actually swim into the trap, and what is the risk of the beetles escaping from the trap.
Elevated wind speed has been identified as a factor influencing the recapture of Culex tarsalis67. Although diving beetles are primarily aquatic, wind speed can also affect these organisms. Wind induced waves affect the movement of both nutrients and various organisms in the water bodies, including diving beetles, which may affect recapturing success. Climatic conditions such as weather and atmospheric pressure can also have an effect.
All the above factors may bias population estimation. Thus, for a more accurate population assessment, that accounts for the influence of all factors, comprehensive models should be created, for example, in the68 or CAPTURE program. Future work may require Bayesian methods as they are increasingly used to analyse ecological data69 including spatial capture–recapture data sets69.
According to the results of the study in Lake Rothemoorsee, we have prepared a suggestion for counting D. latissimus, in the frame of EU monitoring. In order to be able to use the survey results in lakes of different sizes, we related the required number of traps to the length of the shoreline of the current body of water. The shoreline of Lake Rothemoorsee is 1200 m long, and according to our estimates the optimal number of traps for a water body with this shoreline length comprises 45–50 traps. Conversion of this result to 100 m of shoreline gives a ratio of 4 traps per 100 m, which can be used to plan the number of traps to monitor D. latissimus. This number of traps will provide the optimal ratio of data variation and the effort expended on setting the traps. In turn, the minimum number of traps should not be less than 35, as this greatly increases the range of the confidence interval. Estimators based on closed population model provided little fluctuating results with relatively narrow confidence interval after 4–8 censuses, dependently on the recapture rate39.
It seems that the number of traps per month plays a more significant role than the total number of traps per year. For example, it will be more effective to set 35 traps for 3 consecutive months (together 105 traps per year) than 21 traps for 5 consecutive months (together 105 traps per year). As water level plays a significant role in estimating population size, three consecutive years with less than 35 traps per month is not enough to achieve good monitoring results.
Conclusions
Dytiscus latissimus is a species of high conservation value and the need to estimate its abundance is specified in European Union legislation. Prior to the present study, only a single publication on a small relict population in the Netherlands27 provided data for calculating the population size of D. latissimus. Our study provides a larger comprehensive framework for further studies aimed at monitoring and conservation of D. latissimus populations. This publication discusses two approaches to assessing populations of this species. In both studies, sufficient numbers of individuals were captured and recaptured to estimate population sizes, so the method used is appropriate for estimating populations of species. The comparison of data spread between years and countries suggests that setting shoreline traps only leads to overestimation, especially in lakes with an unstable hydrological regime, and needs to be improved to determine the population size more accurately. The absolute population size was determined with the method used in Latvia, which is more practical for smaller lakes, whereas in the larger lake in north-eastern Germany the activity density was determined.
Data availability
Currently data are available from the corresponding author on reasonable request. The datasets analysed during the current study are not publicly available yet because this data will be used to prepare at least one additional publication after which data will be publicly available.
References
Miller, K. On the phylogeny of the Dytiscidae (Insecta: Coleoptera) with emphasis on the morphology of the female reproductive system. Insect Syst. Evol. 32, 45–89 (2001).
Nilsson-Ortmann, V. & Nilsson, A. N. Using taxonomic revision data to estimate the global species richness and characteristics of undescribed species of diving beetles (Coleoptera: Dytiscidae). Biodiv. Inf. 7, 1–16. https://doi.org/10.17161/bi.v7i1.3631 (2010).
Hájek, J. Family Dytiscidae Leach, 1815, pp. 844–914. In: Löbl I. and Löbl D. (eds.): Catalogue of Palaearctic Coleoptera. Volume 1. Revised and updated edition. Archostemata Myxophaga – Adephaga. Leiden/Boston: Brill, XXXIV + 1443 pp. (2017).
Nilsson, A.N. & Hájek, J. A World Catalogue of the Family Dytiscidae, or the Diving Beetles (Coleoptera, Adephaga). Version 1.I.2021. Distributed as a PDF file via Internet (2021). Available from: http://www.waterbeetles.eu (accessed 26 April 2021).
Foster, G.N. Species Status No. 1: A review of the scarce and threatened Coleoptera of Great Britain—Part 3: Water beetles of Great Britain. (2010) http://www.jncc.gov.uk/
Vahruševs, V. & Kalninš, M. Broadest diver Dytiscus latissimus Linnaeus, 1758 (Coleoptera: Dytiscidae) in the Baltic States: A rare or little-known species. Zool. Ecol. 23(3), 203–216. https://doi.org/10.1080/21658005.2013.811906 (2013).
Gioria, M. Habitats (2014), pp. 307–362, DOI: https://doi.org/10.1007/978-94-017-9109-0_7.
Yee, D. Ecology, Systematics, and the Natural History of Predaceous Diving Beetles (Coleoptera: Dytiscidae) (2014), pp. 468.
Nilsson, A.N. Catalogue of Palearctic Dytiscidae (Coleoptera). Internet version 2010–01–01, p. 49 (2010).
Hendrich, L., Müller, R., Schmidt, G. & Frase, T. D. Breitrandkäfer Dytiscus latissimus (Linnaeus, 1768) in Brandenburg–Wiederfund nach über 20 Jahren sowie eine kritische Betrachtung historischer Fundmeldungen und Sammlungsdaten. Naturschutz Landsch. Brandenbg. 21(3), 120–126 (2012).
Hendrich, L., Schmidt, G., Müller, R., Frase, T. Distribution and habitat requirements of the large diving beetle Dytiscus latissimus Linnaeus, 1758 [“Der Breitrand”] (Coleoptera, Dytiscidae) in north-east Germany. Jahrestagung der DGL/GFÖ ausgerichtet von der Universität Potsdam und dem Institut für Gewässerökologie und Binnenfischerei (IGB), 9–13. September 2013 in Potsdam, Abstract Volume: 200 (2013).
Nilsson, A.N., Holmen, M. The Aquatic Adephaga (Coleoptera) of the Fennoscandia and Denmark. Ii. II-Dytiscidea. Brill, Dytiscidae pp. 192 (1995).
Johansson, A. & Nilsson, A. Dytiscus latissimus and D. circumcinctus (Coleoptera, Dytiscidae) larvae as predators on three case-making caddis larvae. Hydrobiologia 248, 201–213. https://doi.org/10.1007/bf00006148 (1992).
Blunck, H. Zur Kenntnis des ‘Breitrands’ Dytiscus latissimus L. und seiner Junglarve. Zool. Anz. 57, 157–168 (1923).
Holmen, M. Fredede insekter i Danmark. III: Biller knyttet til vand. Entomol. Medd. 61(4), 117–134 (1993).
Hendrich, L. Mythos Breitrand—vom Leben und „leisen Sterben“ des zweitgrößten Schwimmkäfers der welt (Dytiscidae: Dytiscus latissimus Linnaeus, 1758). Nachr. Bayer. Entomol. 60(1/2), 2–9 (2011).
Hendrich, L. & Balke, M. Verbreitung, Habitatbindung, Gefährdung und mögliche Schutzmaßnahmen der FFH-Arten Dytiscus latissimus Linnaeus, 1758 (Der Breitrand) und Graphoderus bilineatus (De Geer, 1774) in Deutschland (Coleoptera: Dytiscidae). Insecta 6, 98–114 (2000).
Hendrich, L., Balke, M. Graphoderus bilineatus (De Geer, 1774) (Coleoptera: Dytiscidae). In: Petersen, B., Ellwanger, G., Biewald, G., Hauke, U., Ludwig, G., Pretsher, P., Schroder, E., Ssymank, A. (Hrsg.): Das europäische Schutzgebietssystem Natura 2000. Ökologie und Verbreitung von Arten der FFH-Richtlinie in Deutschland. Band 1: Pflanzen und Wirbellose. Schriftenreihe für Landschaftspflege und Naturschutz. Bonn—Bad Godesberg, 69(1), 388–396 (2005).
Commission of the European Communities. Council directive 92/43/EEC of 21 May 1992 on the conservation of natural habitats and of wild fauna and flora. As amended by the Accession Act of Austria, Finland and Sweden (1995) and the Accession Act of the Czech Republic, the Republic of Estonia, the Republic of Cyprus, the Republic of Latvia, the Republic of Lithuania, the Republic of Hungary, the Republic of Malta, the Republic of Poland, the Republic of Slovenia and the Slovak Republic (2003). Off. J. Eur. Union L. Brussels. European Commission 1992/95/2003.
Foster, G., Bilton, D. The conservation of predaceous diving beetles: knowns, unknowns and anecdotes In: Yee DA (ed) Ecology, Systematics, and the Natural History of Predaceous Diving Beetles (Coleoptera: Dytiscidae). (Springer, London, 2014).
Scholten, I., Kleef, H., van Dijk, G., Brouwer, J. & Verberk, W. Larval development, metabolism and diet are possible key factors explaining the decline of the threatened Dytiscus latissimus. Insect Conserv. Div. 11, 565–577. https://doi.org/10.1111/icad.12294 (2018).
Cuppen, J.G.M., van Dijk, G., Koese, B.,Vorst, O. De Breede Geelgebrande Waterroofkever Dytiscus latissimus in Zuidwest. Drenthe. (European Invertebrate Survey (EIS) Nederland, 2006) p. 59.
Reemer, M., Cuppen, J.G.M., van Dijk, G., Koese, B., Vorst, O. De Brede Geelgerqande Waterroofkever Dytiscus latissimus in Nederland. (European Invertebrate Survey (EIS) Nederland, 2008), p. 33.
Przewoźny, M. & Lubecki, K. Nowe stanowiska Dytiscus latissimus Linnaeus, 1758 i Graphoderus bilineatus (DeGeer, 1774) (Coleoptera: Dytiscidae) w Polsce. Wiad. Entomol. 30, 261–263 (2011).
Scheers, K. Dytiscus latissimus Linnaeus, 1758 and Graphoderus bilineatus (Degeer, 1774) in Belgium: a detailed account of the known records (Coleoptera: Dytiscidae). Bull. Soc. R. Entomol. Bull K. Belg. Ver. Voor Entomol. 151, 34–39 (2015).
Markowski, P. & Buczyński, P. Rediscovery of Dytiscus latissimus Linnaeus, 1758 (Coleoptera: Dytiscidae) on the southern edge of its distribution area in Central and Eastern Poland. Pol. J. Entomol. 89, 81–90. https://doi.org/10.5604/01.3001.0014.1547 (2020).
Koese, B., Cuppen, J., van Dijk, G., Vorst, O. Populatieschatting van de Geelgerande Waterroofkever Dytiscus latissimus in Nederland. (European Invertebrate Survey (EIS) Nederland, Leiden, 2010), p. 33.
Norton-Griffiths, M. Counting Animals: Revised Second Edition. Handbook No. 1. Serengeti Ecological Monitoring Programme. African Wildlife Leadership Foundation. pp. 140 (1978).
Davis, D. E. & Winstead, R. L. Estimating the numbers of wildlife populations. In Wildlife Management Techniques Manual 4th edn (ed. Schemnitz, S. D.) 221–245 (The Wildlife Society, 1980).
Van Lavieren, L. P. Wildlife Management in the Tropics with Special Emphasis on South East Asia: A Guidebook to the Warden. Part 1 Introduction, Taking Fielnotes and Wildlife Census Methods. (School of Environmental Conservation Management, 1982).
Grimm, A., Gruber, B. & Henle, K. Reliability of different mark-recapture methods for population size estimation tested against reference population sizes constructed from field data. PLoS ONE 9(6), e98840 (2014).
Chao, A. Estimating the population size for capture-recapture data with unequal catchability. Biometrics 43(4), 783–791 (1987).
Eberhardt, L. L. Population estimates from recapture frequencies. J. Wildl. Manag. 33, 28–39 (1969).
Otis, D. L., Burnham, K. P., White, G. C. & Anderson, D. R. Statistical interence from capture data on closed animal populations. Wildl. Monogr. 62, 1–135 (1978).
Amstrup, S. C., McDonald, T. L. & Manly, B. F. J. Handbook of Capture-recapture Analyses 336 (Princeton University Press, 2006).
Schwarz, C.J., Seber, G.A. Estimating animal abundance. Rev. III 14, 427–456 (1999)
Williams, B. K., Nichols, J. D. & Conroy, M. J. Analysis and Management of Animal Populations 817 (Academic Press, 2002).
Budrys, E., Budrienë, A. & Pakalniškis, S. Population size assessment using mark-release-recapture of 12 species of orthoptera, diptera and hymenoptera: A comparison of methods. Latv. Entomol. 41, 32–43 (2004).
Romairone, J., Jime’nez, J., Luque-Larena, J. J. & Mougeot, F. Spatial capture-recapture design and modelling for the study of small mammals. PLoS ONE 13(6), e0198766 (2018).
Tourani, M. A review of spatial capture-recapture: Ecological insights, limitations, and 651 prospects. Ecol. Evol. 12, e8468 (2022).
Schmidt, G., Hendrich, L. Abundance and population dynamics of Dytiscus latissimus Linnaeus, 1758 (Coleoptera, Dytiscidae) in north-east Germany, a species listed in the EU Species and Habitats Directive. Jahrestagung der DGL/GFÖ ausgerichtet von der Universität Potsdam und dem Institut für Gewässerökologie und Binnenfischerei (IGB), 9–13. 2013 in Potsdam, Abstract Volume: 188 (2013).
Schnabel, Z. E. The estimation of the total fish population of a lake. Am. Math. Mon. 45, 348–352 (1938).
Jolly, G. M. Mark-recapture models with parameters constant in time. Biometrics 38, 301–321 (1982).
Seber, G.A.F. The Estimation of Animal Abundance and Related Parameters, 2nd edn. Macmillan Publishing Co. Inc. (1982), pp. 654.
Mölle, J. Zur Bedeutung von Amphibienlarven für die Populationsentwicklung des Gemeinen Gelbrandkäfers Dytiscus marginalis, L. 1758. Dissertation. Rheinischen Friedrich-Wilhelms-Universität Bonn (2001), p. 211.
Davy-Bowker, J. A mark and recapture study of water beetles (Coleoptera: Dytiscidae) in a group of semi-permanent and temporary ponds. Aquat. Ecol. 36, 435–446. https://doi.org/10.1023/a:1016550127986 (2002).
Nelson, G. A. Fishery Science Methods and Models. R package version 1.11–2 (2021). Available at https://cran.r-project.org/web/packages/fishmethods/index.html (accessed 25 October 2021).
R Development Core Team R: A Language and Environment for Statistical Computing. Vienna, Austria, R Foundation for Statistical Computing (2021). Available at http://www.R-project.org (accessed 25 October 2021).
Sutherland, W. J. Ecological Census Techniques, A Handbook, 2nd edn. (Cambridge University Press, 2006), pp. 432.
Turlure, C., Pe’er, G., Baguette, M. & Schtickzelle, N. A simplified mark-release-recapture protocol to improve the cost effectiveness of repeated population size quantification. Methods Ecol. Evol. 9, 645–656. https://doi.org/10.1111/2041-210X.12900 (2018).
White, A., Minch, R., Bidder, L. & Gaff, H. A simple, inexpensive method for mark-recapture of Ixodid ticks. J. Insect Sci. 20(6), 9. https://doi.org/10.1093/jisesa/ieaa052 (2020).
Lancia, R. A., Nichols, J. D. & Pollock, K. H. Estimating the number of animals in wildlife populations. In Research and Management Techniques for Wildlife and Habitats 5th edn (ed. Bookhout, T. A.) 215–253 (The Wildlife Society, 1994).
Ranius, T. & Hedin, J. The dispersal rate of a beetle, Osmoderma eremita, living in tree hollows. Oecologia 126, 363–370. https://doi.org/10.1007/s004420000529 (2001).
Pollock, K. H. A capture-recapture design robust to unequal probability of capture. J. Wildl. Manag. 46, 752–757 (1982).
Carvalho, D. D., Collins, P. A. & De Bonis, C. J. The mark recapture method applied to population estimates of a freshwater crab on an alluvial plain. Mar. Freshw. Res. 64, 317–323 (2013).
Baumart, J. S., Cogo, G. B., Morales, F. E. C. & Santos, S. Population size of Aegla longirostri Bond- Buckup and Buckup, 1994 (Crustacea, Decapoda, Anomura): Comparison of methods withthe mark-recapture technique in closed population. Nauplius 26, e2018016-8 (2018).
Johnson, C. G. Migration and Dispersal of Insects by Flight. (Methuen and Co Ltd., 1969), pp. 763.
Zalom, F. G., Grigarick, A. A. & Way, M. O. Diel flight periodicities of some Dytiscidae (Coleoptera) associated with California rice paddies. Ecol. Entomol. 5, 183–187. https://doi.org/10.1111/j.1365-2311.1980.tb01140.x (1980).
Lundkvist, E., Landin, J. & Karlsson, F. Dispersing diving beetles (Dytiscidae) in agricultural and urban landscapes in south-eastern Sweden. Ann. Zool. Fenn. 39(2), 109–123 (2002).
Boda, P. & Csabai, Z. When do beetles and bugs fly? A unified scheme for describing seasonal flight behaviour of highly dispersing primary aquatic insects. Hydrobiologia 703, 133–147 (2013).
Mayne, B., Berry, O. & Jarman, S. Redefining life expectancy and maximum lifespan for wildlife management. Austral. Ecol. 45, 855–857 (2020).
Field, A. Discovering Statistics using IBM SPSS Statistics 4th edn. (SAGE Publications, 2013).
Bröder, L. et al. Optimization of capture–recapture monitoring of elusive species illustrated with a threatened grasshopper. Conserv. Biol. 34(3), 743–753. https://doi.org/10.1111/cobi.13449 (2020).
Burns, P. A., McCall, C., Rowe, K. C., Parrott, M. L. & Phillips, B. L. Accounting for detectability and abundance in survey design for a declining species. Divers. Distrib. 25(10), 1655–1665. https://doi.org/10.1111/ddi.12966 (2019).
Crawford, J. A., Dreslik, M. J., Baker, S. J., Phillips, C. A. & Peterman, W. E. Factors affecting the detection of an imperiled and cryptic species. Diversity 12, 177 (2020).
Willson, J. D., Winne, C. T. & Todd, B. D. Ecological and methodological factors affecting detectability and population estimation in elusive species. J. Wildl. Manag. 75, 36–45 (2011).
Reisen, W. K., Lothrop, H. D. & Lothrop, B. Factors infuencing the outcome of mark-release- recapture studies with Culex tarsalis (Diptera: Culicidae). J. Med. Entomol. 40, 820–829 (2003).
Cooch, E., White, G. A Gentle Introduction (2nd edition) Program MARK. Analysis of data from marked individuals (2001). http://www.phidot.org/software/mark/docs/book.
Dey, S., Delampady, M. & Gopalaswamy, A. M. Bayesian model selection for spatial capture–recapture models. Ecol. Evol. 9, 11569–11583 (2019).
Acknowledgements
The authors express their gratitude to Dr. Garth Foster for helpful comments. We wish to thank all the colleagues, for their technical assistance in the field work of this study.
Funding
Work in Latvia was supported by the LVAF project Nr. 1/27/2020–8/27/2020 "Improvement and approbation of the methodology for monitoring and estimating the population size of three EU protected invertebrate species—Dytiscus latissimus, Graphoderus bilineatus and Hirudo medicinalis". Work in Germany was supported by the "University of Applied Science Neubrandenburg and the Ministry of Climate Protection, Agriculture and Environment Mecklenburg-Vorpommern". The second author would like to say a special word of thanks to Prof. Dr. Grünwald (†) for his invaluable encouragement and academic stimulus.
Author information
Authors and Affiliations
Contributions
M.B., K.S., L.H. and G.S. wrote the main manuscript text. M.N. prepared cartographical materials. M.Z. analysed data. K.A, K.K, M.B, U.V, G.S. and L.H collected data. All authors reviewed the manuscript.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Publisher's note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Balalaikins, M., Schmidt, G., Aksjuta, K. et al. The first comprehensive population size estimations for the highly endangered largest diving beetle Dytiscus latissimus in Europe. Sci Rep 13, 9715 (2023). https://doi.org/10.1038/s41598-023-36242-w
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-023-36242-w
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.