**Preface to "Hotspots of Subterranean Biodiversity"**

Caves have long intrigued biologists because of the unique life forms they harbor—species without eyes or pigment, often with elongated, spindly appendages and other non-optic sensory hypertrophies. Caves and their inhabitants were once thought to be rare, but current research shows that neither is truly rare. Worldwide, there are well in excess of 100,000 caves and 5000 species specialized for caves and other aphotic, subterranean habitats. These species hold a special fascination for evolutionary biologists because of their bizarre morphology, and they present challenges to students of biodiversity because sampling cave fauna represents a considerable logistical challenge and because species ranges are typically very small due to the constraints of subterranean dispersal. More than 20 years ago, David Culver, together with the Slovenian speleobiologist, Boris Sket, assembled the first list of individual caves and karst wells (20 in total) with more than 20 specialized species. Since that time, this list of hotspot sites has grown to 22 caves and karst wells, each with 25 or more aquatic or terrestrial species (and in a few cases both) specialized for caves and wells.

In this Special Issue, detailed accounts are given of 14 of the 22 hotspot sites in Australia, Bermuda, Mexico, France, Romania, Slovenia, Spain, Sulawesi, and the United States. For each of these sites, a detailed taxonomic consideration is presented, as well as detailed description of the site, usually with a cave map. These descriptions should be of interest, not only to speleobiologists, but to anyone interested in biodiversity in general. Although the pattern is still incompletely known, the terrestrial pattern is a global one, with hotspots throughout the world. Aquatic hotspots are more constrained, somewhat concentrated in central Europe (six out of ten cases).

Our thanks go the publisher, especially Emma Li, the Managing Editor for *Diversity*, and to all the authors who contributed to this Special Issue—Hotspots of Subterranean Biodiversity—printed here.

> **Tanja Pipan, David C. Culver, Louis Deharveng** *Editors*

## *Editorial* **Hotspots of Subterranean Biodiversity**

#### **Tanja Pipan 1,\*, Louis Deharveng 2,\* and David C. Culver 3,\***


Received: 18 May 2020; Accepted: 21 May 2020; Published: 25 May 2020

**Abstract:** Worldwide, caves and groundwater habitats harbor thousands of species modified and limited to subterranean habitats in karst. Data are concentrated in Europe and USA, where a number of detailed analyses have been performed. Much less is known with respect to global patterns due to a lack of data. This special issue will focus on and discuss the global patterns of individual hotspot caves and groundwater habitats.

**Keywords:** α-diversity; biogeography; biospeleology; cave biology; caves; hotspots; invertebrates; subterranean biodiversity

#### **1. Introduction**

The denizens of caves and other habitats without light hold a special interest for students of biodiversity [1]. The fauna have bizarre and very distinctive morphology, including the lack of or reduction in eyes and pigment, elongated appendages, and often a ghost-like appearance [1–4]. How these species from disparate groups, including both vertebrates and invertebrates, evolved this convergent morphology remains a subject of intense research, both in the laboratory [5–7] and in the field [8,9]. Whatever the details, cave fauna are a magnificent example of convergent evolution and adaptation [10,11].

What is less well known is that there are thousands of species specialized for life in darkness, especially in caves. While speciation can occur as a result of isolation in caves, followed by subterranean dispersal and more isolation, opportunities for subterranean dispersal are restricted [4,12,13], and the specialized subterranean fauna is the result of many hundreds if not thousands of separate colonizations from the surface to the subsurface. These numbers are sufficient to expect predictable geographic patterns to emerge from their analysis [14].

While there are now several datasets on the georeferenced occurrence of specialized subterranean species from both Europe and the United States with over 2000 records [15–18], these records were hard-won and without parallel in other continents and regions. Brazil is the best-known region outside of Europe and the United States and the number of records is much smaller, and mostly for as yet undescribed species [19,20]. Exploration and access to caves is often formidable, and the number of caves is very large. For example, 45,000 caves were known from the USA in 1999 [21], and that number has grown considerably since that time. Of course, not all caves need to be sampled to obtain a good estimate of species richness for a region, but a rather large number do. This is because β-diversity dominates α-diversity [22,23]. One study of the European terrestrial fauna found that approximately 100 caves were needed to have accumulation curves of species numbers approach an asymptote [24].

The availability of large numbers of records for non-cave subterranean habitats is almost non-existent, with two exceptions. One is the analysis of epikarst communities in Slovenia [25], a habitat that is indirectly sampled through the continuous sampling of water dripping from cave ceilings [26]. The other is the larger scale sampling of porous aquifers in Europe as part of the PASCALIS project (Protocol for the ASsessment and Conservation of Aquatic Life in the Subsurface) [27]. There are also several studies that combine cave and non-cave subterranean aquatic diversity [28,29]. This sampling deficiency in subterranean habitats not directly accessible to humans has been named the "Racovitzan impediment" [30].

In spite of these impediments, we know a great deal about cave biodiversity patterns in Europe and, to a lesser extent, the United States. In Europe, there is a ridge of high species richness for both aquatic and terrestrial species at a latitude of approximately 45◦ N along the spine of the Pyrenees and through the Dinaric karst of Slovenia, Croatia, Serbia, and Montenegro [16,24]. Species range sizes follow Rapoport's Rule and increase with increasing latitude [16], and species turnover is high over short distances [16,23]. Explanations of patterns are highly scale-dependent. For example, the ridge of high species richness is largely concordant with areas of high productivity [24], but productivity is unimportant in determining hotspots within hotspots [15]. The explanations for regional differences have proven to be complex, often involving a combination of energy availability, spatial heterogeneity, and history, and the explanations vary from region to region, being spatially non-stationary [31].

Determining the global patterns of subterranean species richness has remained elusive. At one level, the possibilities of determining the patterns seem difficult at best and a long way in the future. If thousands of records are needed to describe and understand the European pattern, global studies with tens of thousands of records would be needed. However, there may be a shortcut at hand that at least holds promise of a global picture, albeit incomplete. It starts with the observation by Gibert and Deharveng [22] in their classic paper on subterranean biodiversity that regional diversity is a good predictor of local diversity and vice versa. This was buttressed by later findings that species accumulation curves rarely crossed, and thus the regional qualitative pattern could be captured by a relatively small number of samples [24,32]. Culver and Sket [33] took this idea to its logical extreme, and considered only caves and karst wells with the highest species richness, originally finding 20 sites with 20 or more species specialized for subterranean life. While the coverage of large numbers of caves in a relatively small area was (and is) limited to Europe and the United States, they reasoned that at least a few outstanding caves, extensively sampled, were known from most large karst areas. Since the publication of the first hotspot list in 2000, knowledge of the global cave fauna has grown exponentially. Species lists are available for several tropical countries [20,34,35], and a number of caves throughout the tropics and sub-tropics are now well sampled [36]. In a later update [37], the bar was raised to either 25 terrestrial or 25 aquatic species. In 2019, there were 24 examples of hotspot caves known [38]—16 from the temperate zone, 5 from the subtemperate zone, and 3 from the subtropics. The tropics have at least five caves with 20 or more species specialized for subterranean life, the original hotspot criterion [37]. The demonstration of hotspot caves in the tropics also raises some other issues. One is how to treat undescribed species, which is the case for the majority of tropical species [39]. If they are ignored, then tropical and subtropical caves will appear to be artificially depauperate. If they are fully counted, the likelihood that some of these species are not valid [40] is ignored. Of course, described species may also be wrongly thought to be limited to subterranean habitats. In addition, the tropical cave fauna often has a component specialized on guano but never found outside of caves. These guanobionts typically show less eye and pigment loss and less appendage elongation [38,39]. Should they be discounted as cave specialists? If so, then once again tropical and subtropical caves may appear to be artificially depauperate.

This special issue of *Diversity* aims to bring together information on hotspot caves and karst ground-watered habitats for in-depth analyses and comparisons. First and foremost, there will be a species list for each of the hotspot caves, information that is strangely unpublished for many of the hotspot caves. This is especially important given the controversy around the ecological status of cave species [41]. Deharveng and Bedos [41] pointed out that considerable confusion exists in the literature about the terms troglobiont—which should be used only for species not found outside of caves, irrespective of their morphology—and troglomorph [42,43], species with reduced eyes and pigment and elongated appendages. The two are not identical, a problem that arises not only with guanobionts but also with all species without troglomorphic features occurring in caves [44]. The special issue will also provide a physical setting for the caves and groundwater habitats, including their hydrogeological and environmental context, their use by humans, the nature of the karst in which they are situated, and the knowledge on nearby cave biodiversity. An interesting side note is that a number of hotspot caves are or have been commercial caves, including Mammoth Cave in Kentucky, USA; the Postojna-Planina Cave System in Slovenia; and Vjetrenica in Bosnia and Hercegovina. Meanwhile, some groundwater hotspots are industrially exploited for water consumption. Finally, there will also be a summary of what we know, the patterns of distribution, and future research directions. For more information, contact any of the editors.

**Conflicts of Interest:** The authors declare no conflict of interest.

#### **References**


© 2020 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).

**Ana Isabel Camacho 1,\* and Carlos Puch <sup>2</sup>**


**Abstract:** Ojo Guareña Natural Monument in Burgos (Spain) is an important and large karstic system. It consists of more than 110 km of surveyed galleries, and it has rich sources of organic material from the surface and permanent water circulation. It is the fourth largest cave system in the Iberian Peninsula, and one of the 10 largest in Europe. Ojo Guareña also ranks 23rd among the world's largest caves. To date, only volcanic caves in the Canary Islands, in which between 28 and 38 subterranean species occur, are considered subterranean diversity hotspots in Spain. Here, we provide the first list of subterranean taxa present in Ojo Guareñ, which is comprised of 54 taxa that includes 46 stygobiotic and eight troglobiotic species (some still unidentified at the species level), revealing Ojo Guareña as the largest known subterranean biodiversity hotspot in Spain and Portugal. In addition, we provide a list of an additional 48 taxa, 34 stygophiles and 14 troglophiles, found in the system, whose ecological status could change with detailed biological studies, which may change the number of strictly subterranean species present in the system. Indeed, at present, these numbers are provisional as they correspond to a small part of this sizeable cave system. The biodiversity of large areas of the system remains unknown as these areas have yet to be explored from the biological point of view. In addition, a large number of samples of both terrestrial and aquatic fauna are still under study by specialists. Furthermore, evidence of cryptic species within Bathynellacea (Crustacea) indicates an underestimation of biodiversity in the karstic system. Despite these limitations, the data available reveal the typical uneven distribution of subterranean aquatic fauna, and suggest that the great heterogeneity of the microhabitats in this wide and highly connected karstic extension led to the great richness of aquatic subterranean species.

**Keywords:** hotspot caves; Ojo Guareña natural monument; stygobionts; troglobionts; subterranean biodiversity

#### **1. Introduction**

Iberian Peninsula, together with Balearic Islands, is one of the regions in Europe with the greatest development of karst areas. Numerous caves, many of them large (more than 3 km of development and/or more than 300 m in depth), have been surveyed in the country [1]. Despite this, very little data on the subterranean fauna within these caves are available. Potential biospeleologists are discouraged by the difficulties of access and progression in the caves (e.g., large vertical shafts, meanders, narrow passages, etc.), the inaccessibility to man of much of the underground environment (flooded galleries, mesovoid shallow substratum-MSS-, hyporheic sites), and the low density of specimens in accessible populations (epikarst, especially). Thus, most samplings of these ecosystems have been limited and sporadic, and mainly focused on the study of specific taxonomic groups using capture techniques suited for them. Due to the great sampling deficiencies and the scarce biospeleological tradition in our country, there are very few inventories of the region's underground fauna beyond some lists of a few terrestrial or aquatic taxa in the more well-known caves. The only caves considered in the literature as subterranean diversity

**Citation:** Camacho, A.I.; Puch, C. Ojo Guareña: A Hotspot of Subterranean Biodiversity in Spain. *Diversity* **2021**, *13*, 199. https:// doi.org/10.3390/d13050199

Academic Editors: Tanja Pipan, David C. Culver and Louis Deharveng

Received: 31 March 2021 Accepted: 6 May 2021 Published: 8 May 2021

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

**Copyright:** © 2021 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https:// creativecommons.org/licenses/by/ 4.0/).

hotspots in Spain are four volcanic ones located in the Canary Islands: in these, between 28 and 36 troglobiotic species and 38 stygobiotic and interstitial species (volcanic anchialine habitats of Lanzarote) have been recorded [2,3]. Culver and Pipan (2013) also mention possible diversity hotspots in the Cantabrian Mountains in northern Spain, which contain the largest number of caves in the country [4]; however, more studies of this region are needed to confirm this hypothesis.

The Ojo Guareña Natural Monument (OGNM) karstic system is situated on the southern rim of the Cantabrian Mountains. It lies within a hypothetical zone (between ca 42◦ and 46◦ N) in Europe characterized as having "high biodiversity for terrestrial cave fauna" [5]. According to Gulden [6], OGNM is the fourth largest cave system in the Iberian Peninsula, and one of the 10 largest in Europe; it also ranks 23rd on the list of the world's largest caves. With more than 110 km of surveyed galleries, it is a vast cave system that has rich sources of organic material from the surface and also permanent groundwater flow (phreatic water). It is one of the few places in Spain where the terrestrial and aquatic environments have been sampled repeatedly and extensively [7]. However, we are far from knowing its true diversity. All of the "shortfalls" impeding our knowledge of biodiversity that are commonly mentioned at the global level (e.g., Linnean, Wallacean, Darwinian, Prestonian) [8] also apply at local scales. A more recently described shortfall, the "Racovitzan impediment" [9], is of particular importance in the case of the subterranean biodiversity of the OGNM: what remains unexplored, from a biological point of view, cannot be known.

The most accessible (and thus the best sampled) area of the OGNM extends from the main entrance at Palomera Cave to Museo de Cera Hall (covering about 3 km of passages). Moreover, this area contains sumps from the Guareña River and water that filters in from the Trema River. However, many of the specimens collected in this area and throughout the monument remain unidentified at the species level (taxonomic shortfall). Cryptic species of crustaceans have also been discovered in the area [10], indicating that the data presented here represent an underestimation of the diversity of the subterranean fauna. In addition, the biological composition of many kilometers of underground passages has yet to be surveyed, and the accumulation curves of stygobitic species [7,11] of the best sampled areas do not reach saturation, demonstrating our incomplete knowledge of the true biodiversity of this cave system's fauna.

In addition to its biological assemblages, the OG endokarst has unique geological, geomorphological, hydrological, paleoclimatic, archaeological, and paleontological characteristics that make it worthy of protection. It has a continuous and non-fragmented surface that is subjected to natural evolution with little human intervention and that can maintain the physical and biological characteristics of the system, ensuring the continued functioning of natural processes. Moreover, the absence of resource extraction by agricultural, forestry, hydraulic or mining (natural or artificial) activities reduce the likelihood that the functioning and resources of the ecosystems within the system, and the aesthetics of the landscape, would be significantly altered.

#### **2. A Bit of History**

The main cave of the system (OG), through its main entrances (Palomera Cave and Dolencias Shaft), has been known since the 1940s [4]. Although the cave did not initially arouse great interest among the scientific community—beyond simple sporadic visits and specific samplings—, interest later increased as the importance of its biological species [12–15], archaeological sites [16], and hydrogeological and geomorphological peculiarities became more evident. Due to these features, various public administrations began to establish general protection regimes for the karst system; eventually, in 1996, it was incorporated as a natural monument in the Network of Protected Natural Spaces in Spain and, more recently, in the European Ecological Network Natura 2000 [7]. In addition, given the uniqueness of its geology, it was included in a publication on the "Points of Geological Interest of the Cantabrian Coast" (1983). However, despite these acknowledgments, very

little scientific study has been conducted in the system. In fact, until just this century, faunal sampling of the system had been sparse, and the findings almost anecdotal.

Here, we provide some details of past explorations of this natural monument in order to put into context the knowledge that exists of this vast karst system. The presence of numerous archaeological remains from the Paleolithic and Neolithic, and from the Iron and Middle ages, testifies that OG has been known since ancient times. However, geologists did not become interested in it until 1933 [17]. In 1956, the Edelweiss Cave Club (GEE), from Burgos, began to explore the large galleries that converge at the main entrance (development = 8 km). In 1958, with the collaboration of various speleological groups (Spanish and foreigners), the "OG 58" campaign was held, resulting in 8900 m being surveyed and produced the first plan of the cave. A complete topographic survey of the system continues to this day, and in parallel, exploration inside the main cave and the rest of the caves intensifies. By 1986 [18], surveys of new sectors, only accessible when the water is low, together with the union of other caves, considerably increased the known development of the complex to 89,071 m. Currently, 110 km in an area encompassing 15 km<sup>2</sup> have been explored and surveyed.

Biological sampling of OG was occasional during the second half of the twentieth century, though it is worth highlighting the initiative of Professor Eugenio Ortiz, director of the Museo Nacional de Ciencias Naturales de Madrid (MNCN), who, from 1968 to 1975, tried to establish an underground laboratory at the system (that never materialized) and an associated biological station to be sponsored by the Consejo Superior de Investigaciones Científicas (CSIC). Among the aquatic faunal samplings carried out during those years, Ortiz observed the presence of numerous specimens of the isopod *Stenasellus virei buchneri* Stammer, 1936 (Figure 1) in a pond near the Palomera entrance. He also collected specimens of Bathynellacea and some other groups; however, they were never studied or, at least, the results were never published. Other samples collected by Ortiz, especially terrestrial ones, were later studied by specialists of different fauna groups (e.g., oligochaetes, opilions, mites, myriapods, and beetles). Later (1975), T. Antón, a member of the Edelweiss Cave Club, and X. Bellés collected and studied some other terrestrial fauna, which led to several publications [19–26]. In 1984, Notenboom [27,28] reported on some aquatic faunal samples from one of the resurgences of the system (La Torcona). Camacho and Puch also conducted sporadic samplings of aquatic fauna in the cave during the 1980s and 1990s. Thanks to these works, which were performed without institutional support, the cave's fauna began to be inventoried [29]. As of 1993, a total of 81 taxa were known in the cave, of which 71 were identified to the species level (63 terrestrial and 8 aquatic). Of these, only a dozen of the well-identified species could be considered as strictly subterranean fauna [29]. From 2002 to 2004, OGNM was one of four areas under study as part of a project on European stygofauna (PASCALIS) [11]. The project areas included various caves, springs, wells, unsaturated zones (hyporheic system of epigean rivers) and porous aquifers. Finally, the Junta de Castilla y León and the MNCN of the CSIC signed agreements to carry out faunistic studies of the cave (from 2002 to 2004 and 2006 to 2009) prior to its limited opening to tourism. Small parts of the terrestrial [30] and aquatic ecosystems [31] were studied prior to any potential effects that tourism could have on the faunal composition. The results of all of these investigations have contributed to the lists of species constituting the main objective of this publication.

In this study, we provide the first lists of stygobiotic and troglobiotic (and stygophilic and troglophilic) taxa present in OGNM, which were compiled on the basis of the results of the aforementioned investigations. Although only a small fraction of this vast karst system has been surveyed, and some gaps in knowledge remain, we show that it is a hotspot of subterranean biodiversity with great biological potential.

**Figure 1.** Isopod *Stenasellus virei buchneri*. Photo ESEM, MNCN, Madrid.

#### **3. Materials and Methods**

The study area (OGNM) is located in the Sotoscueva region (Burgos, Northern Spain), within the southernmost limit of the Cantabrian Mountains. OGNM is a paradigmatic example of a complete karst system, with a water absorption zone in the northwestern and central parts, an extensive network of galleries of more than 110 km in the main cave and evacuation points located at the southeastern limit of the system. It develops in limestones and dolomites of the Upper Coniacian (about 130 m of thickness), which lean on clay limestones and marls of the Middle Coniacian–Turonian (Upper Cretaceous) that act as an impermeable substrate (Figure 2). The surface landscape is arranged in the form of a geosyncline with a WNW–ESE alignment, giving rise to the characteristic slopes that stand out in the northern part of the area [4,18,32]. The rivers Guareña, to the north, and Trema, to the west, course underground when they encounter limestones and dolomites and resurface, after a long, partially known route, in La Torcona, a cave resurgence near the confluence of the Trema River with La Hoz Stream (Figure 3).

**Figure 2.** Geological section approximately North-South, corresponding to line A–B in Figure 3.

**Figure 3.** Location of caves, underground galleries, flows, and drainages on a digital elevation model of the OGNM terrain.

The main cave (OG) has 13 entrances: Palomera Cave, Dolencias Shaft, Huesos Shaft, Rizuelos Shaft, Cornejo Cave, Cuatro Pisos Cave, San Bernabé caves, Moro Cave, Trema river sinkholes, La Mina Cave, Guareña river sinkhole (the "Ojo" or Eye), Villallana Cave, and Torcona cave spring (main resurgence of the system) (Table S1) (Figure 3). The last two are connected to the main cave by long flooded passages. The general structure is comprised of a maze of galleries, mainly horizontal, that are arranged along the W–E and N–S axes. In the western part of the cave, the galleries are mainly oriented in a NE–SW direction. Other caves of the system that are not yet connected to the main cave include Kubía (surveyed length 550 m), Último Sumidero (350 m), Prado Vargas (130 m), Kaite (585 m), Covanería (320 m), and Jaime (650 m). A little further away from the main cave lies Las Yeguas Cave (1900 m). In 1981, the Trema river sinkholes (on the right bank of the river) and La Mina Cave (on the left bank) became connected through two conduits that traverse underneath the bottom of the Trema river canyon. The hypogeal course of this river is only known in the space between two sumps (300 m). In Cornejo Cave, the progression of the underground stream stops in a semi-flooded squeeze.

In short, the current ensemble of galleries developed in a 15 km2 area (Figure 3). Of these, 30 km are concentrated in an area of just 1 km2, underneath Alto de San Bernabé Hill, just above and to the south of the Guareña river sinkhole. The galleries are distributed in six largely overlapping levels, creating a three-dimensional maze known as "the West Daedalus". The upper levels are almost clogged inactive passages. The intermediate ones present a seasonal hydrological activity that can be very important, depending on the weather. The lower levels are the permanent beds of the Guareña and Trema rivers, to the north and west, respectively, and of Villamartín Stream, to the northwest, feeding a large underground aquifer installed around the axis of the syncline. After a long, partially known path, these water sources reappear in La Torcona, a cave near the eastern limit of the massif, and El Torcón, an inaccessible flooded spring. Hydrological recharge occurs by infiltration of rainwater and snowmelt through permeable formations, as well as by feeding from small hillside springs and by karst sinks that induce flow losses along epigeal channels. The most spectacular cases of recharging in the entire unit are found at the Guareña river sinkhole, at the base of the San Bernabé escarpment, through a penetrable sinkhole ("el Ojo"), and at the Trema river sinkhole in its own channel, shortly before the

village of Cornejo. The entire endorheic basin represents an area of about 27 km2. After a relatively rapid underground flow (1.5 km/day in low water conditions and 4 to 5 km/day in flooding), the flow drains towards the Trema River, on the southeastern limit of the karst. Although much is known about the behavior of the underground waters that come from the Guareña River, very little is known about those that come from the Trema River.

#### *3.1. Sampling and Studies of the Fauna*

As mentioned above, the OGNM was only occasionally sampled prior to the first decade of the twenty-first century, when the most intensive and systematic samplings, to date, occurred. The taxa lists compiled here are derived from previous studies of the material obtained from all known samplings, though with a bias towards the more recent ones. Below, we briefly mention some of the sampling techniques used, as well as the extent of sampling.

#### 3.1.1. Samples of Aquatic Fauna

Although we do not know the number or specific locations of the samples collected or the specific sampling techniques used by Ortiz during his exploration of the cave, we do have detailed information on the systematic samplings carried out during this century. Between 2002 and 2009, 344 samples of aquatic fauna were collected at OGNM [31]. Sampling was carried out in the main cave (OG), plus 12 other caves and 12 springs, and in the hyporheic environments of the epigean rivers (Table S1):

(a) Main Cave (OG), about 3 km of main gallery (Palomera to Museo de Cera). A total of 59 sites were sampled (see Table S1 for site names): OG01 to OG16, periodically sampled (methods described below), and OG17 to OG59, sampled once or twice using hand nets (Figure 4). A total of 244 samples was studied.

**Figure 4.** Part of the survey of the main cave, indicating the sampling sites.

(b) Other caves in the system that were sampled on one or two occasions: Rizuelos, Kaite 2, Jaime, La Mina, Prado Vargas, García, San Bernabé, San Miguel river sinkhole, Las Llanas, Cornejo, Racino, and Redonda. A total of 31 samples were studied.

(c) Springs sampled only once: Torcona, Pozo del Infierno, Cornejo, La Mea, Villa, Salce, Cubío, La Calzada, Jordana, Avellanos, Mazo 1, and Mazo 2. A total of 24 samples were studied.

(d) Hyporheic and interstitial environments of the surrounding rivers and streams that are part of the OGNM aquifer, sampled on one or two occasions and at various points: Guareña, Ulemas, Trema, Trueba, Engaña, Nela, and La Hoz. A total of 43 samples were studied.

(e) Wells sampled only once: Torcón and Villabascones. Only two samples were studied.

Most samplings were restricted to the epikarstic zone of the caves [31]. However, sampling also occurred in some other areas including a couple of underground rivers in the main cave (Erizos and Guareña), the phreatic zone in Aburrimiento gallery (a long southernmost passage), and the resurgences of La Torcona and Villallana.

Eight sampling methods were used [33,34]: (1) traps with bait in free water (ponds and gours) and located within the substrate on riverbanks; (2) direct captures with a manual aspirator; (3) drip collection; (4) the Karaman-Chappuis method (on sandy shores of underground pools and epi- and hypogeal rivers); (5) removing the substrate and filtering the water with a hand net (ponds, gours, potholes) (Figure 5); (6) Bou-Rouch pumping (in La Torcona and in epigean rivers); (7) kicking (in the benthos of the Guareña and Erizos river streams); and (8) Cvetkov nets or other types of phreatobiological nets (in deep lakes and wells). All sampling nets (0.100 mm mesh size) were designed ad-hoc in order to adapt them for the different habitats and types of fauna studied. With the mesh size used, most adult forms of crustaceans and interstitial fauna can be collected.

**Figure 5.** Sampling in "Sala Edelweiss" OG09. Photo C. Puch.

The samples collected between 2002 and 2004 and between 2007 and 2009 were processed in the MNCN laboratories in Madrid under the supervision of A.I. Camacho [31]. As a great deal of the material collected (mainly ostracods, amphipods and copepods) are still being studied and identified by specialists, some of the entries on the stygobiotic species list have only been identified to the family or genus level. In most cases, samples were identified using classical (morphological) taxonomic methods. For samples of Bathynellacea, molecular analyses of nuclear (18S) and mitochondrial (COI) gene fragments were also carried out for taxonomic purposes. For this, we used DNA extraction and PCR amplification techniques previously developed for this type of fauna [10,35,36].

#### 3.1.2. Samples of Terrestrial Fauna

The sampling techniques used in the earlier sporadic sampling visits to the cave are not known. We assume that most samplings were made de visu, a low yield technique, hence the scarcity of available material and data on the few troglobiotic species known at that time. From 2002 to 2004, samples of terrestrial fauna were collected using three techniques:

(1) Baited traps, after being set, were collected periodically (30, 60, 90, and 120 days after setting) in four selected spots in OG, along the main route from Palomera to Museo de Cera, over a distance of about 3 km. Samples collected outside the main cave (OG), from the access doline to the Palomera entrance, [29] were also included in the study.

(2) Sediment samples were collected at four locations in OG at six time points throughout the year, in an attempt to collect cave-dwelling fauna that are not attracted to the baited traps. Sediment samples of the access doline to Palomera, taken over four occasions, were also processed.

(3) Manual aspirators were used to capture samples de visu on every visit to OG.

The collected material was processed in the laboratories of the Department of Animal Biology I of the Faculty of Biology at the Complutense University of Madrid, under the supervision of Pérez-Zaballos [30]. Fauna were extracted from the sediment using the Berlese technique, and taxa were separated into individual groups. In a few cases, specimens could be identified only to the order or family level, and most specimens remain under study by taxonomists specializing in the different terrestrial groups. Some of the material has been deposited in the collections housed at the MNCN (Madrid, Spain).

#### **4. Results**

Overall, only partial results have been obtained from the samplings: most of the material collected, particularly the terrestrial specimens, are still under study, with their identification to a specific level still pending.

A total of 299 taxa has been identified in OGNM (Table 1), of which 112 are terrestrial and 187 are aquatic [30,31]. A total of 54 of these taxa constitutes subterranean fauna: 46 stygobiotics and 8 troglobiotics. However, another 48 taxa comprising 34 stygophilics and 14 troglophilics (Table S2) were also identified. In many cases, the ecological status of cave animals is difficult to determine [37] and therefore we feel their inclusion is warranted. In the future, the surface/subterranean status of some of these species could change as more information about their ecology and biology comes to light.

**Table 1.** Number of taxa known in the Ojo Guareña Natural Monument.


#### **Table 1.** *Cont.*


#### *4.1. Aquatic Fauna*

Table 1 shows the distribution of the 187 aquatic taxa in OGNM by faunal group. Crustacea is the largest group, with 78 species, 30 of which are stygobiotics (some without a specific name). Groups with the largest species numbers include Copepoda, Ostracoda and Bathynellacea. Only five of the 49 aquatic species of Oligochaeta are stygobionts. In the case of the Acari, nine of the 20 aquatic species are stygobiotics. The two remaining stygobiotics species are Mollusca (two gastropod species). Similar to the case of the stygobiotics, most of the 34 stygophilic species belong to the Copepoda (12 species), Ostracoda (9), and Oligochaeta (8) (Table 1).

Although this work focuses on the subterranean species present in OGNM, it is important to highlight the repeated presence of the up to 187 taxa in the samples, and their unequal distributions and abundances among the different subterranean aquatic habitats and sampling points. As an important energetic resource of the ecosystem, this aquatic fauna may be able to provide information on the diversity and disparity of populations in different localities [31].

Below, we provide a list of stygobiont taxa found in OGNM. All of the samples have been reviewed by expert taxonomists; thus, we consider the list as valid, without specifying whether the "sp." notation refers to a new species, an already known species, or, possibly, a cryptic species. The specific site(s) from which each taxa was collected in OGNM is also included.

List of stygobiont species found in the Ojo Guareña National Monument, and their specific sites (sampling sites in the main cave, through the Palomera entrance: OG01–OG59, Figure 4).

#### Oligochaeta

*Gianius navarroi* Rodriguez and Achurra, 2010. OG10 and OG14. *Delaya navarrensis* (Delay, 1973). OG21. *Stylodrilus mariae* Achurra, Rodriguez and Erséus, 2015. OG07 and OG08. Trichodrilus tenuis Habre, 1960. OG14. *Trichodrilus* sp. 1. OG09.

#### Mollusca

*Palaospeum septentrionale* (Rolan and Ramos, 1995). Sima Jaime Cave, De la Hoz Creek, and Torcona Spring.

*Spiralix (Burgosia) burgensis* Boeters, 2003. Sima Jaime Cave, Trema River, and Trueba River.

#### Hydrachnidia

*Axonopsis (Paraxonopsis) vietsi* Motas and Tanasachi, 1947. Trueba River. *Frontipodopsis reticulatifrons* Szalay, 1945. OG08, OG10, and Engaña River. *Frontipodopsis* sp. Lasía River. *Kongsbergia* sp. Trueba and Lasía rivers. *Neoacarus hibernicus* Halbert, 1944. Nela River. *Protzia* sp. Trueba River. Pseudotorrenticola sp. Nela River. *Sperchonopsis* sp. Lasía River. Limnohalacarida (Halacaridae)

*Soldanellonyx chappuisi* Walter, 1917. OG: 07–10, 12, 14, 15, 17, 19, 30, 31, 34, 36–38, and 50; caves: Rizuelos, García, San Miguel, Cornejo and Redonda; springs: Torcona and Calzada; rivers: Trema, Trueba, Engaña, Nela, and De la Hoz and Torcon Well.

#### Copepoda, Cyclopoida

*Acanthocyclops* cf. *venustus* (Norman and Scott, 1906). Engaña and Nela rivers.

*Acanthocyclops* n. sp. OG: 10, 31, 34, 37, and 38.

*Graeteriella* (*G*.) *unisetigera* (Graeter, 1910). OG02 and OG07.

*Speocyclops infernus* (Kiefer, 1930). OG: 01, 08, 09, 14, 16–19, 21, 31, 38, and Huesos 2 and 4.

*Speocyclops sebastianus* Kiefer, 1937. Trema River.

*Speocyclops spelaeus* Kiefer, 1937. Caves: García, San Bernabé, Las Llanas, and Racino. Copepoda, Harpacticoida

*Bryocamptus* (*Rheocamptus*) *pyrenaicus* (Chappuis, 1923). Caves: Kaite 2, García, Cornejo, Racino, and Redonda; Salce Spring.

*Paracamptus* sp. SA. Nela River.

*Parastenocaris* sp. OG: 07–09, 12, 14, 16, 17, 20, 21, 31–38, 41–43, 45, 46, 51–54, and Huesos 2; Villa Spring and De la Hoz and De la Cueva creeks.

#### Ostracoda

Candoninae cf. gen. sp. 1 "Rounded". OG: 07, 09, 12, 17, 19, 21, 39; spring: Torcona, Cornejo and Calzada; Nela River.

Candoninae "Trapezoid" 1. OG: 09, 19, 21; caves: Jaime, La Mina, Las Llanas; Guareña River.

Candoninae cf. gen. sp. 2 "Trapezoid" 2. OG09.

Candoninae cf. gen. sp. 2 "Trapezoid" 3. OG: 1, 18, 37, 38.

Candoninae "Triangular" 1. OG: 15, 30–32 and Trema River.

Candoninae "Triangular" 2. OG16 and Sima Jaime Cave.

Candoninae "Triangular" 3. Torcona Spring and Trema River.

*Cypria* sp. OG: 7, 9, 10, 35, 40, 42, 43, 51, 55–59, Huesos 2 and 3; San Miguel Cave; springs: Cornejo, Villa, Salce, Cubio, and Jordana; rivers: Guareña and Lasía.

#### Amphipoda

*Haploginglymus* sp. 1. OG38; Torcona Spring; river: Trema, Trueba, and Nela; Torcon well.

*Niphargus* sp. OG37 and 38; Cubio Spring and Trueba River. *Pseudoniphargus burgensis* Notenboom, 1986. Torcona Spring.

*Pseudoniphargus* n. sp. 1. Torcona Spring.

#### Isopoda

*Cantabroniscus primitivus* Vandel, 1965. OG21.

*Stenasellus virei buchneri* (Stammer, 1936). OG: 02, 07, 09, 10, 14, 15, 17, 21, 39, 40, 44, 53, 55, and 57; caves: García, San Bernabé, San Miguel, Las Llanas, and Redonda and Salce Spring.

#### Bathynellacea

*Iberobathynella burgalensis* Camacho, 2005. OG53.

*Iberobathynella cornejoensis* Camacho, 2005. Redonda and Racino caves and Trema River. *Iberobathynella guarenensis* Camacho, 2003. OG57 ("Erizos River").

*Vejdovskybathynella edelweiss* Camacho, 2007. OG01, 09, 16, 17, 38, 50, 57, and Huesos 3 and 4; Racino and Mina caves and Cubio Spring

*Vejdovskybathynella* n. sp. 1 (cryptic species "*edelweiss*"). OG57 ("Erizos River").

*Vejdovskybathynella* n. sp. 2 (cryptic species "*edelweiss*"). Redonda and Sima Jaime caves.

Bathynellidae n. gen. n. sp. OG09, 14, Torcona Spring and Prado Vargas Cave.

In OG, 25 stygobiotic species have been found across the 59 sampling sites, with many of the species being found in OG09 (12 species), followed by OG07, OG08, and OG10 (11 species each) and OG14 (seven species) (Figure 4).

Of the five stygobiotic species of Oligochaeta found in OG, *Gianius navarroi* described from OG [38], has only been found in two small ponds, OG10 and OG14. The two species of Mollusca live in other caves, springs and rivers but not in OG. The eight stygobiotics of Hydrachnidia live in the rivers, though one, *Frontipodopsis reticulatifrons*, has also been found in OG8 and OG10, two small epikarstic ponds that are filled with water from diffuse drips. The stygobiotic species of Limnohalacaridae is widely distributed throughout the principal cave (in many of the sites), and in all of the other types of habitats (multiple springs and rivers, and one of the two wells).

Crustaceans, a major animal group found in groundwater worldwide [2], are very well represented in OGNM by 78 taxa, 30 of which are stygobionts (Table 1). Copepods are present in almost all ponds and pools in OG. Of the 37 copepod species found, nine are stygobiotics that live in several caves and rivers and in a couple of springs. Twentytwo identified species of ostracods are distributed around OGNM. Of these, eight are stygobiotics that inhabit all the different types of sites sampled. Species of three of the five genera of Amphipoda found at the system are stigobiotics. They live mainly in springs and rivers, and also in two sites of OG. Of the two stygobiotic isopod species

found, *Cantabroniscus primitivus* has only been found in the "siete lagos" area (OG21). The other species, *Stenasellus virei buchneri* (Figure 1), has been found in many sites in OG and in some of the other caves and springs. In addition to these isopod species, we found *Proasellus* cf *ortizi*, (Table S2), which we conservatively consider a stygophilic, despite its cave-like morphology. This species has been consistently found in OG12, although at varying abundances, and in three springs (Avellanos, Mea, and Jordana).

Bathynellacea, a strictly stygobiotic group of crustaceans, is represented in OGNM by seven species belonging to two families: Parabathynellidae and Bathynellidae (Figures 6 and 7). Four of the species, three belonging to Parabathynellidae (*Iberobathynella* Schminke, 1973 genus) and one to Bathynellidae (*Vejdovskybathynella* Serban and Leclerc, 1984 genus), have already been formally described [12–14,39]. The second genus found belonging to Bathynelidae is pending description. On the basis of molecular data, we discovered two cryptic species of *Vejdovskybathynella* within the family Bathynellidae [40]. These species are currently under further study. To our knowledge, no other karst area in the world has as many known species of Bathynellacea living together. The specimens of Bathynellidae found in Erizos River and Redonda Cave were initially identified as *Vejdovskybathynella edelweiss* (Figure 6), which is known from seven sites in OG, and from other caves of the system (Huesos, Racino, La Mina, and Cubio Spring). It is the most widely distributed species in OGNM. However, the specimens from Erizos River show a genetic divergence in the COI gene of 14.8% with the nominal species, and 17% with the specimens from Redonda Cave, which, in turn, show a 15.2% divergence with the nominal species [40]. These results imply that the specimens from Erizos and from Redonda each constitute a distinct species, and that there are actually three morphologically similar species of *Vejdovskybathynella* in the sampled sites. We expect this case will not be an isolated event, suggesting the biodiversity of the system is far greater than has been estimated. The distribution of the bathynellaceans is patchy. Species of both families have been found together in Redonda Cave and Erizos River (the only site from which *I. guarenensis* has been described). *Vejdovskybathynella edelweiss* and the new genus coexist in OG09. *Iberobathynella burgalensis* (Figure 7) only lives in gours in OG53. *Iberobathynella cornejoensis*, which has not been found in the main cave, inhabits a couple of small caves of the complex including Redonda (near the Trema sinkholes, where it has also been found) and Racino (to the west).

**Figure 6.** Habitus of *Vejdovskybathynella edelweiss* Camacho, 2007.

**Figure 7.** Habitus of *Iberobathynella burgalensis* Camacho, 2005.

#### *4.2. Terrestrial Fauna*

The terrestrial fauna in OGNM is poor compared with those in other systems: only 8 troglobiotic and 14 troglophilic species have been identified (Table S2). This is more likely due to a lack of sampling and study rather than to an absence of species. Given that we still have some samples waiting to be identified to the species level, we expect this number to increase.

Between 1968 and 1975, Ortiz sampled the terrestrial and aquatic fauna in Palomera. The material collected from two clayey areas of Palomera and from zones with organic residues throughout OG were subsequently studied by specialists [22–25] and gave rise to two monographs, one on terrestrial oligochaetes (eight species but only two troglobionts) [41] and another on oribatid mites (43 species but no troglobionts; three were new to Spain, and two were new to science) [42].

Other previous studies of terrestrial fauna in OG include those by Bellés and Antón [18]; Prieto and Gómez [43], who studied some terrestrial molluscs collected by the Niphargus Cave Club. Demange and Serra [44] described the only myriapod known from OG. Rambla cites the presence of a troglobiont opilion in OG [45] that Prieto and Zubiaga later studied [46]. Español [47] studied 10 coleopterans collected by Ortiz and found six species, one of them new to science (*Trechus ortizi)*. The presence of some of these species in OG were confirmed through the samplings carried out between 2002 and 2003; however, most of the specimens from these later campaigns have yet to be identified by specialists.

The list of well-identified troglobiotic species is short, and specific information on the sites in which they have been found is lacking.

List of troglobiotic species found in the main cave of the Ojo Guareña National Monument:

Oligochaeta:

*Aporrectodea rosea troglodyta* (Álvarez, 1971).

*Orodrilus paradoxoides* (Álvarez, 1971).

Mollusca:

*Zospeum schaufussi* Frauenfeld, 1862.

Arachnida, Opiliones:

*Litocampa zaldivarae* Sendra, Salgado and Monedero, 2003.

Myriapoda, Chilopoda:

*Lithobius deroutae sexusbispiniger* Demange and Serra, 1978.

Collembola:

*Verhoeffiella* cf *hispanicus* (Bonet, 1931).

Coleoptera:

*Speocharis sharpi* Escalera, 1898. *Trechus ortizi* Español, 1970.

#### **5. Discussion and Conclusions**

Only a small fraction of the terrestrial ecosystem of OGNM has been studied: knowledge of the species that occur in the complex stems from sporadic samplings conducted during the second half of the twentieth century and between 2002 and 2003, with many of the specimens collected still unstudied. Therefore, a few data available about the troglobiotic fauna in the system makes any discussion irrelevant. By contrast, the abundance of information on the subterranean aquatic fauna of OGNM makes their discussion important, despite the fact that only a small part of the system's aquatic environment has been studied. Over a total of six years (2002–2004 and 2006–2009), a series of selected sites from the epikarst of the main gallery had been regularly sampled, with many other areas of the vast subterranean territory sporadically sampled.

The information available for the stygofauna goes beyond that provided by a still photo taken during sampling. Repeated samplings at different times of the annual cycle and over successive years (thus considering changes in environmental conditions or variability of drought-flood cycles) are necessary to find a stygofauna that is formed, in general, by "rare" species (i.e., those that are scarce or difficult to locate) that live in small populations in places with little to no access.

The obtained data reveal the typical patchy distribution of subterranean aquatic fauna [48]. The great heterogeneity of microhabitats, in a highly connected, vast continuous karstic extension, could explain the high richness of subterranean species observed in OGNM, as has been demonstrated in other cases [49,50]. Habitat heterogeneity and the spatiotemporal dynamics of this subterranean aquatic ecosystem are also evidenced by the great variation of species and communities observed, even over short periods of time and at a local scale of a few meters [31], a result common for subterranean environments [51–53].

Despite numerous samplings, covering the different microhabitats (gours, drips, ponds, lakes, and the interstitial medium of hypogean streams), large sections of OGNM remain unexplored, from both a topological (e.g., permanently flooded sections) and biological perspective. The inability to sample the entirety of the karst system implies an underestimation of its species richness.

The subterranean aquatic fauna in most karst areas in Europe has been fairly well studied, with around 1800 stygobiotic species known to date [54]. In Spain, slightly more than 200 are known from very few areas (representing 11% of all known species in Europe) [11]. In OGNM, of the 187 aquatic taxa found (see Table 1), 46 are stygobiotics, accounting for almost a quarter of all stygobiotic species known in our country. However, many of these taxa have been identified to only the genus level (e.g., for some taxa of Oligochaeta, Copepoda, Amphipoda, and Acari), or the subfamily level (e.g., Candoninae (Ostracoda)), and many specimens remain unidentified, suggesting the number of stygobiotic species will likely increase as more studies are published. We do not yet know the extent of endemism or the number of new species in OGNM, as many specimens from all of the major groups remain under study and have not yet been identified to a specific level.

Crustaceans constitute a major animal group in groundwater [2]: around 3400 species are known worldwide, 1200 of which are known to occur in Europe [54]. We identified 30 stygobiotic crustacean species in OGNM. Remarkably, we found seven species of

Bathynellacea in OGNM. This is a unique case as in no other karst system has this many species been found in coexistence. Moreover, the seven species are distributed among three genera: two known ones (*Iberobathynella* Schminke, 1973, and *Vejdovskybathynella* Serban and Leclerc, 1984, each with three species) and a third (Bathynellidae n. gen. n. sp.) that is new to science. Specimens of these species have been found, within the main cave, in five pools and the Erizos River, and in five of the other caves, two springs, and the Trema River. The two cryptic species of *Vejdovskybathynella*, *Vejdovskybathynella* sp. 1 and *Vejdovskybathynella* sp. 2, have also been found in two nearby caves outside OGNM, Río Chico, and Imunía [55]. The other five species of Bathynellacea are endemic to OGNM. A possible scenario is that the ancestral populations of these species, which were likely widely distributed in the area, became isolated due to droughts and vicariance. Although Bathynellacea lack active forms of dispersal (they do not have swimming larvae, like other crustaceans) and are generally not very mobile, a later (passive) dispersal event may have led to their genetic isolation and consequent speciation, resulting in the unusually high number of species coexisting, though in a patchy way, in a relatively small area. However, phylogeographic studies of these species are needed to confirm any of these hypotheses.

The largest number of stygobiotic species found in OGNM occur in the main cave (OG), with 25 species. A large number of these species is found at several sites: 12 species in OG09; 11 in OG07, OG08, and OG10; and seven in OG14. In addition, in general, the species were repeatedly found in all samplings, although in different proportions [31]. The 25 stygobiotics, along with the seven known troglobiotics, make OG a subterranean biodiversity hotspot within the wider OGNM hotspot (with its 54 subterranean species).

The number of stygobiont species in OGMN is striking. In no other area of Spain and Portugal has such a level of subterranean biodiversity been reported (see Table 1), even in the four hotspot caves in the Canary Islands (Felipe Reventón, Viento, Sobrado, and La Corona lava tube), which have between 28 and 36 troglobiotic species the first three and 38 stygobiotic species La Corona [2,3].

To improve our knowledge of the diversity of a region, it is essential to have large inventories of well-identified, georeferenced species in order to conduct standardized geospatial analyses that can be compared among regions and at different spatial scales [2,56]. Although the available data clearly demonstrate that the subterranean biodiversity in OGNM is relatively high, especially compared with that of other caves in the world [57,58]. The balance of sites and specific samples studied shows the relevance that the small known fraction can have compared with all that remains to be explored and studied from a biological perspective.

Among the common shortfalls of biodiversity data, the "Racovitzan impediment" is of great importance in the case of the OGNM. The sampling deficiency is large: there are numerous habitats not directly accessible to humans that, therefore, have not been sampled to date. Thus, we lack information about the species they may harbor. In addition, the lag between sampling efforts and the taxonomic identification of all specimens makes it difficult to obtain accurate data on the overall biodiversity of the system and its distribution. Although we are likely underestimating the real diversity of subterranean fauna in OGNM due to these issues, we confirm the prediction made by Culver and Pipan [2]: a subterranean hotspot is indeed present in northern Spain. Here, we provide the first list of stygobiotic and troglobiotic taxa present in OGNM, as well as demonstrate that, despite some knowledge gaps, this complex is a hotspot of subterranean biodiversity. Moreover, the diversity observed to date for OGNM implies its high rank among hotspot caves in not only Europe, but the world.

Future lines of research should focus on resolving the ecological status of the 48 taxa provisionally considered as stygophilics (34) or troglophilics (14) (see Table S2); discovering potential new cryptic species among the fauna through molecular analyses, such as those used to reveal cryptic species in the families Bathynellidae [55] and Parabathynellidae [59]; and identifying all remaining material to the lowest taxonomic level possible. These studies, along with a greater sampling effort, will all contribute to increase our knowledge of the subterranean biodiversity of this exceptionally conserved karst ecosystem.

**Supplementary Materials:** The following materials are available online at https://www.mdpi. com/article/10.3390/d13050199/s1: Figure S1, *Stenasellus virei buchneri* (Stammer, 1936) in OG14; Table S1, Georeferenced list of sampled sites on the Ojo Guareña Natural Monument; Table S2, List of stygophile and troglophile species and the sites where they have been found in Ojo Guareña Natural Monument.

**Author Contributions:** Both authors participated in all phases of the realization and writing of the paper. Both authors have read and agreed to the published version of the manuscript.

**Funding:** This work was supported by project PID2019-110243GB-100 (Ministerio de Ciencia e Innovación/FEDER).

**Institutional Review Board Statement:** Not applicable.

**Informed Consent Statement:** Not applicable.

**Data Availability Statement:** All data are available from the authors upon request.

**Acknowledgments:** We express our gratitude to J. Fernandez and B. Sanchez-Chillón for reviewing their bibliographic and specimen collections, respectively, for this study. We thank the Spanish PAS-CALIS team, as well as Fortunato Lázaro, Ana M. de Juan, and Jesús Robador for their participation in the samplings. Thanks to Pilar Rodríguez and Sanda Lepure for reviewing and correcting the list of oligochaetes and copepods, respectively. Thanks also to Consuelo Temiño and Beatriz Cabezas for their help from 2002–2010 under the agreements of the Junta de Castilla y León and the CSIC, and to the MNCN and the CSIC for their support in carrying out this work. And lastly, thanks to Miguel Rioseras for providing us with the surveys and cartographic data and Melinda Modrell, who helped us with the English translations.

**Conflicts of Interest:** The authors declare no conflict of interest.

#### **References**

