Next Article in Journal
A Comparison of DNA Metabarcoding and Microscopy Methodologies for the Study of Aquatic Microbial Eukaryotes
Previous Article in Journal
Planarians, a Neglected Component of Biodiversity in Groundwaters
Previous Article in Special Issue
Updated Distribution of the Mysid Antromysis cenotensis (Crustacea: Peracarida), a Protected Key Species in Yucatan Peninsula Cenotes
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Primer on Spider Assemblages in Levantine Caves: The Neglected Subterranean Habitats of the Levant—A Biodiversity Mine

1
The National Natural History Collections, Edmond J. Safra Campus, Givat Ram, The Hebrew University of Jerusalem, Jerusalem 9190401, Israel
2
The Department of Ecology, Evolution and Behavior, Edmond J. Safra Campus, Givat Ram, The Hebrew University of Jerusalem, Jerusalem 9190401, Israel
3
Mitrani Department of Desert Ecology, Swiss Institute for Dryland Environmental and Energy Research, Blaustein Institutes for Desert Research, Sede Boqer Campus, Ben-Gurion University of the Negev, Midreshet Ben-Gurion 8499000, Israel
*
Author to whom correspondence should be addressed.
Diversity 2021, 13(5), 179; https://doi.org/10.3390/d13050179
Submission received: 25 March 2021 / Revised: 18 April 2021 / Accepted: 20 April 2021 / Published: 23 April 2021
(This article belongs to the Special Issue Cave Communities: From the Surface Border to the Deep Darkness)

Abstract

:
Caves share unique conditions that have led to convergent adaptations of cave-dwelling animals. In addition, local factors act as filters on regional species-pools to shape the assemblage composition of local caves. Surveys of 35 Levantine caves, distributed along a climate gradient from the mesic in the north of Israel to hyper-arid areas in the south of Israel, were conducted to test the effect of cave characteristics, location, climate, bat presence, and guano level on the spider assemblage. We found 62 spider species and assigned four species as troglobites, 28 as troglophiles, and 30 as accidentals. Precipitation, elevation, latitude, minimum temperature, and guano levels significantly affected the composition of cave-dwelling spider assemblages. Caves situated in the Mediterranean region had higher species richness and abundance, as well as more troglobite and troglophile arachnids. These discoveries contribute to the knowledge of the local arachnofauna and are important for the conservation of cave ecosystems. By comparing spider assemblages of Levantine caves to European caves, we identified gaps in the taxonomic research, focusing our efforts on spider families that may have additional cryptic or yet to be described cave-dwelling spider species. Our faunistic surveys are crucial stages for understanding the evolutionary and ecological mechanisms of arachnid speciation in Levantine caves.

Graphical Abstract

1. Introduction

Different caves, similar patterns. Subterranean habitats (caves and other hypogean habitats) around the world can be found in different climates, rock formations, and biogeographical regions. Additionally, they can be formed by various means such as volcanic, glacial, mechanical, and erosion/solution processes. However, the majority of these subterranean habitats share unique abiotic conditions such as a limitation of light, stable and narrow range of temperature, and high relative humidity [1,2,3]. The peculiar abiotic conditions found in subterranean habitats, together with specific regional and local factors, determine species richness and assemblage composition of a particular cave. Processes at the regional scale include geological and climatic events together with historical biogeography, dispersal, extinction, and speciation, which shape the regional species pool [4,5,6]. At the local scale, ecological interactions and local abiotic conditions act as filters on the regional species pool to shape each cave assemblage composition [7,8,9,10].
Depending on the number of cave openings, their size and location, as well as the size and depth of the cave, a gradient of light intensity, climatic conditions and nutrients can be found, creating up to five defined zones: entrance, twilight, transition, deep, and air-stagnant zones [1,11,12,13]. The abiotic conditions at the deep and air-stagnant zones of caves are less affected by seasonality and daily cycles and therefore are more stable environments in comparison to the cave entrance, twilight, and transition zones [14,15,16]. There is some evidence that different cave zones are inhabited by species with different environmental requirements and constitute discrete assemblages [14].
Connectivity of epigean and hypogean habitats also affect assemblage composition inside the cave via colonization and dispersal within and between epigean and hyopgean systems [17]. While it is suggested that species richness of troglobites (organisms obligated to life in caves) is explained by historical biogeography, species richness of troglophiles (organisms with strong affinity to caves) and accidental or occasional visitors in caves may be explained by local ecological factors [18]. Regions of the world that have rich regional species pools coupled with diverse cave formations and a range of abiotic conditions are expected to have rich assemblages of cave-dwelling visitors and resident species (both troglophiles and troglobites).

1.1. The Levant and Its Caves

The Levant is a diverse biogeographic unit formed by the northeastern African and northwestern Arabian plates and the eastern Mediterranean Levantine basin [19]. Its geographical location at the junction of three continents (Europe, Asia, and Africa), and diversity of habitats and climate zones (Mediterranean, steppe, and arid) lay the foundation for a diverse regional species pool.
The Levant has a diversity of caves that differ in microclimate, age and the type of their rock substrate [20]. The different rock substrates, with intensive karstification in the north of Israel [21], and a climatic gradient from a mesic alpine climate at the north of Israel (Mt. Hermon) to the arid Negev desert in the south of Israel, resulted in a gradient of Karst features with very few or no karst features in the southern Negev desert and Arava. Furthermore, both epigenic and hypogenic caves (formed above or below the water table, respectively) [22,23] can be found in the Levant in basaltic bedrock [24], sandstone, limestone, dolostone, marl, chalk, chert [25,26,27,28], and salt rocks [29]. While many of the caves in the Levant are natural, some are manmade, such as burial caves, and some chambers of natural caves were formed or enlarged as part of a secondary use by man [30].

1.2. Spider and Other Arachnid Assemblages in Caves

It is very common to find troglophile arachnids in the entrance of Levantine caves, as in caves of other regions of the world. Arachnids (and among them spiders) are often numerically dominant in caves and are considered dominant predators in many cave foodwebs [31,32,33,34]. Of the 11 extant terrestrial arachnid orders, five orders (i.e., Araneae, Opiliones, Palpigradi, Pseudoscorpiones, and Scorpiones) and the polyphyletic sub-class Acari were reported to have troglobite species around the world [35,36,37,38], while troglobite species of three orders (i.e., Amblypygi [39], Ricinulei, and Schizomida) were reported only from subtropical and tropical regions [35].
The arachnofauna of caves in the Levant is poorly known in comparison with European caves. In the present study, we investigated the spider assemblages, and specifically the presence of accidentals, troglophile, and troglobite species in eastern Levantine caves (Israel and Palestine) in relation to geographic location, physical characteristics, and ecological zonation (entrance, twilight, and deep) of the caves, as well as the presence of bats and guano in the caves. The data were collected in field surveys between 2013 to 2015. Several new family records for Israel and species new to science discovered in the course of this field survey were reported elsewhere [38,39,40,41], including a distribution model for a common Levantine troglophile pholcid spider: Artema nephilit Aharon, Huber, and Gavish-Regev (2017) [17], while some are reported here for the first time. Here, we used the data collected in the same ecological field survey to describe the spider assemblages in 35 caves along a rainfall gradient and analyzed the environmental factors shaping these assemblages.

2. Materials and Methods

2.1. Field Surveys

Arachnids were collected from 35 caves located in Israel and Palestine (West Bank) (Figure 1) from three cave ecological zones: the cave entrance (inside the cave in the vicinity of the entrance); the twilight zone (in the intermediate part of the cave when it was applicable); and the inner dark zone (when it was applicable), as well as outside each cave (near the cave entrance). The caves are distributed along the climatic gradient from the Mediterranean (mesic) climate in the north and center of Israel (12 caves for each region, among them 6 caves in Palestine, Figure 1 and Figure 2A–F), to the arid and hyper-arid climate in the south of Israel (11 caves, Figure 1 and Figure 2G–I). Data on the caves is summarized in Appendix A.
In each cave ecological zone, temperature and illumination were measured. The temperature was measured using PicoLite 16-K and a single-trip USB temperature logger (FOURTEC), with measurements taken once an hour for 74–77 days. The illumination was recorded at the time of each sampling using an ExTech 401025 Lux Light Meter. The light meter was positioned on the ground until the reading stabilized for a minimum of 1 min. Measurements inside caves ranged between 0 to 420 lux, while measurements outside caves ranged between 60 to 70,000 lux (Supplementary Table S1). Temperature and illuminance were also measured outside the cave. Caves were assigned to three cave size categories: large (13 caves, more than 50 m with twilight and dark zones), medium (12 caves, 11 to 50 m with twilight zone, lacking a dark zone), and small (10 caves, 10 m or shorter, lacking twilight and dark zones) (Figure 1). Cave size and length were estimated from cave maps when available (Israel Cave Research Center) or in the field, and representing the distance from the cave opening toward the darkest region of the cave we could reach. Elevation and geological data were provided by the GIS (geographic information system) center at the Hebrew University of Jerusalem. Annual average precipitation data was taken from the Israel Meteorological Service (https://ims.gov.il/en/ClimateAtlas, accessed on 1 February 2021) from 1980–2010, which came from the closest meteorological station to each cave (Appendix A).
Each of the 35 caves was sampled according to a specific protocol twice during 2014 (6 March–6 April 2014 and 22 May–22 June 2014). Due to differences in cave morphology (including microhabitat, fractal shape of the substrates, size, and volume) we standardized our sampling effort by time. Therefore, our protocol included a 20-min thorough visual search by one of three experienced arachnologists in 3 to 10 m long sectors using headlamps and UV lights in each cave ecological zones. In the first visit to each cave, most of the arachnid observed were collected by hand for further identification in the lab. In further visits to each cave, we thought some species populations were more sensitive or common, thus finding it possible to identify a species level in situ. Such data were recorded and not collected.
We visited several additional caves outside of the formal 2014 cave survey and collected additional arachnids. On these occasions, only the locality and date were documented for the arachnids. Therefore, we refer to the arachnid species collected in these caves only in the general paragraph of the results and discussion, but not in the analysis of the survey data.

2.2. Statistical Analysis

Generalized linear models (GLMs) were used to test the effect of the geographical region (north, center, and south of Israel) on spider richness (Poisson distribution and log link) and abundance (normal distribution and log link) using the entire dataset of troglobite, troglophile, and accidental spider species.
Using a set of GLMs with binomial distribution (logit link), we tested separately for each spider guild [42] the probability that a spider species observed in a cave belongs to the particular spider guild, and whether the probability changes in relation to the chosen environmental variables. Geographic region can encompass different environmental variables. In Israel, there is a strong north to south gradient of several correlated climatic variables such as precipitation, air humidity, evaporation, and air temperature. We added the presence of bats in caves as an additional explanatory variable as bat activity changes the habitat for other cave dwellers, such as spiders. To further test for such effects, we included an alternative model in which, instead of bat presence, we used an estimation of guano level on the cave floor from our field assessment, i.e., none, low (when less than 15% of the cave floor had guano deposition), medium (when 15% to 40% of the cave floor had guano deposition), or high (when more than 40% of the cave floor had guano deposition).
Two multivariate analyses methods were performed in order to better understand the factors affecting the assemblage composition of troglobite and troglophile species: hierarchical clustering, using R v.4.0.3. [43] and a direct CCA ordination (canonical correspondence analysis), using Canoco [44,45]. For the cluster analysis, pairwise Pearson correlation between the cave communities was determined using the ‘cor’ function with the ‘pairwise.complete.obs’ method, and the output hierarchically clustered as a Euclidean distance matrix with an ‘average’ linkage using the ‘hclust’ function. The resultant clustering was plotted using the ‘heatmap.2′ function of the ‘gplots’ package [46], with scaling by row. For the CCA ordination, we used the unrestricted Monte Carlo permutation tests (4999 runs) and forward selection, testing 9 explanatory environmental variables: geographic location (latitude and longitude, continuous), elevation (continuous), cave estimated length (continuous), minimum temperature in the cave during the survey (continuous), average annual precipitation (continuous), rock category of the cave (6 categories: basalt, carbonate rocks, chalk, marlstone, salt, and sandstone), minimum age of the rock (continuous), bat inhabitance (three categories: no bats, insectivorous, frugivorous), and guano level (four categories: no guano, low level, medium level, and high level of guano).
For both multivariate analyses, only 34 of the 35 caves were used, as ‘Ammude ‘Amram (a small cave, 34) had two Solifugae individuals, but no spiders were found in our two visits. We restricted the statistical analyses to troglobite and troglophile species (Nspecies = 32), omitting the accidental species (Nspecies = 30). The latter were epigeic species that occurred in low abundance and were scattered among the different caves.

3. Results

3.1. Cave Arachnofauna from Field Surveys and Additional Visits to Caves

A total of 1132 arachnids were found during the surveys, comprised mainly of species of the order Araneae (1054 individuals). However, representatives of additional arachnids were also found (Figure 3), i.e., Acari (65 individuals, among them two soft-tick species (Argasidae) Argas vespertilionis (Latreille, 1796) (Figure 3A) and Ornithodoros tholozani Laboulbène and Mégnin, 1882 (Figure 3B), as well as hard-ticks (Ixodidae)); Pseudoscorpiones (four individuals, awaiting identification, Figure 3C,D); Amblypygi (although many individuals were observed, only three were collected: two Charinidae species Charinus ioanniticus (Kritscher, 1959) and Charinus israelensis Miranda, Aharon, Gavish-Regev, Giupponi, and Wizen, 2016, Figure 3E) [39]; Scorpiones (although more individuals were observed, we were able to catch only three juveniles (Buthidae)); Solifugae (two juveniles; identification was not possible) and one Opiliones of the species Mediostoma haasi (Roewer, 1953) (Nemastomatidae).
In other visits to several additional caves outside of the 2014 survey, we found one additional Opiliones species, the troglobitic eyeless Haasus naasane Aharon et al., 2019 (Pyramidopidae, Figure 3F) [38], an Opilioacaridae species (in preparation, Figure 3G, [47]), Palpigradi species (in preparation, Figure 3H, [47]), and one spider family that was not reported before from Israel: Nesticidae (the invasive species Eidmannella pallida (Emerton, 1875) and an additional unidentified species, both reported here for the first time).
The 1054 spiders found during the 2014 survey in 34 of the 35 caves were identified to 62 species and morphospecies in 38 genera and 22 families (Supplementary Table S2). As in many faunistic surveys, the Levantine cave spider rank-abundance curve was skewed, with nine species that have 33–152 individuals each and account for 72% of the overall abundance (762 individuals out of 1054), 29 singleton, and doubleton species (42 individuals out of 1054). Moreover, 24 species have 3–28 individuals each (249 individuals out of 1054) (Figure 4). Of the 62 species, 46 are either species with valid taxonomic status (26 species), are in the process of description (six species), or need further taxonomic study (14 species: eight linyphiids, two theridiids (Steatoda Sundevall, 1833), two pholcids (Pholcus Walckenaer, 1805 and Spermophorides Wunderlich, 1992), one theraphosid (Chaetopelma Ausserer, 1871), and one oonopid (Megaoonops Saaristo, 2007)). The rest (16 morphospecies) were not adults and could be identified only to the genus or family level based on morphology (Supplementary Table S2).

3.2. Cave Geographic Region, Bat Presence and Spider Richness and Abundance

The geographical region of the caves significantly affected spider richness (Nspecies = 62; Nindividuals = 1054, χ 2 2 = 27.41 ,   p < 0.001 ) . Specifically, spider richness in the southern region was significantly lower than in the other geographical regions monitored (Figure 5A, Table 1). Bat presence did not significantly affect spider richness ( χ 1 2 = 0.002 ,   p = 0.968 ) . Similarly, a non-significant effect was found when the estimate of guano level was used (guano: χ 2 2 = 3.845 ,   p = 0.279 ). Spider abundance followed a similar trend as richness, showing a significant effect of geographic region on abundance ( χ 2 2 = 7.127 ,   p = 0.028 ), with fewer individuals detected in southern caves (Figure 5B, Table 1), and no significant effect of bat presence ( χ 1 2 = 1.884 ,   p = 0.170 ) ). The similar model with bat guano level was also not significant (geographical region: χ 2 2 = 6.328 ,   p = 0.487 ; guano level: χ 2 2 = 2.438 ,   p = 0.487 ).

3.3. Troglobites and Troglophiles

We assigned each of the 62 spider species in our survey to one of three categories: troglobite, troglophile, and accidental, according to their distribution and the known use of caves as a habitat in Israel and Europe [33]. We also categorized them based on troglomorphic phenotype. This resulted in a list of 32 troglophile and troglobite spider species inhabiting caves in our region and 30 accidental spider species (Supplementary Table S2). Only three spider families included troglobite species: Agelenidae, with two eyeless species found during visits to several additional caves outside of the 2014 survey, and two eye-reduced species in the genus Tegenaria Latreille, 1804; Dysderidae with one eyeless species in the genus Harpactea Bristowe, 1939, which may be a species complex based on our preliminary morphological study; and Leptonetidae, with one species Cataleptoneta edentula Denis, 1955 (Figure 6 and Figure 7) [41].
In total, 11 spider families in Israel and Palestine included troglophile species (Figure 8 and Figure 9). Pholcidae has 8 troglophile species and is family with the largest number of troglophile species in Israel and Palestine (Figure 7). It includes some abundant species found in caves along the north–south climatic gradient (Figure 8): Artema nephilit Aharon, Huber and Gavish-Regev 2017 (Figure 8A) [9,40], Hoplopholcus cecconii Kulczyński, 1908 (Figure 8B), Holocnemus pluchei Scopoli, 1763 (Figure 8C), and a species in the genus Pholcus Walckenaer, 1805 (Figure 8D). However, the three most abundant troglophile species are Tegenaria pagana C. L. Koch, 1840 (Agelenidae, Figure 8E), Filistata insidiatrix (Forskål, 1775) (Filistatidae, Figure 8F), and Loxosceles rufescens (Dufour, 1820) (Sicaridae, Figure 8G), respectively (Figure 4 and Figure 8). We found additional troglophile species (Figure 4, Figure 7, Figure 8 and Figure 9) in the families Theridiidae (Steatoda triangulosa (Walckenaer, 1802), Figure 8H), Linyphiidae (six or seven species of the sub-family Micronetinae, Figure 9A–D), Theraphosidae (Chaetopelma Ausserer, 1871 species, Figure 9E), Phyxelididae (Phyxelida anatolica Griswold, 1990, Figure 9F, [41]), Sparassidae (Heteropoda variegata (Simon, 1874), Figure 9G), Dysderidae (Harpactea Bristowe, 1939, Figure 9H), and Nesticidae. Supplementary Table S2 includes the full list of the spider species with their localities and distribution.

3.4. Spider Foraging Guilds

The spiders found in the cave survey represent six foraging guilds [42]: sheet-web weavers, sensing-web weavers, space-web weavers, orb-web weavers, ambush hunters, and other hunters (Supplementary Table S2). Three out of four troglobite species and 24 out of the 28 troglophile species are web-builders (sheet-web: Agelenidae, Linyphiidae, Phyxelididae; space-web: Leptonetidae, Pholcidae, Theridiidae; sensing-web: Filistatidae, Theraphosidae), and only one troglobite and four troglophiles are hunters (Dysderidae, Sicariidae, Sparassidae).
GLM analyses of the geographic region and bat presence or guano level showed that neither factor in the two alternative models was significant (Table 2).

3.5. Levantine Cave Spider Assemblages: Similarities and Environmental Variables from Field Survey

We used only the 32 troglobite and troglophile spider species to further analyze the spider assemblages in this cave survey (34 caves). We identified five groups of caves that cluster together based on the correlation between their species assemblages (Figure 10). The most noticeable clusters are based each on one of five dominant (T. pagana, F. insidiatrix, and A. nephilit) and sub-dominant (S. triangulosa and H. cecconii) troglophile species. One cluster was distinctive from the others by most sharing a single species, A. nephilit. This cluster includes caves in arid climates, all of the Dead Sea caves (caves 23, 26–29), one cave from the Negev desert (cave 33), and one cave from the Arava (cave 35). Another cluster was based on F. insidiatrix (with either L. rufescens or H. pluchei), with three upper Galilee caves (caves 3, 5, 8) and four caves from the Judean desert, Galilee, and Karmel (caves 9–10, 20, 22), respectively. All caves in another cluster (which could be divided to three sub-clusters) had T. pagana in them (caves 7, 12–13, 15–18, 21). The fourth cluster was based on S. triangulosa and includes two caves from the Negev desert (caves 30–31) and a cave from the Karmel (cave 11). Two caves from the upper Galilee (caves 1–2) form the fifth cluster based on H. cecconii.
A CCA ordination showed five significant variables explaining 91.4% cumulative percentage variance of the species–environment relation: precipitation (F-ratio = 2.52, p = 0.0002), elevation (F-ratio = 2.37, p = 0.0002), latitude (F-ratio = 2.03, p = 0.0032), minimum temperature (F-ratio = 1.79, p = 0.0302), and high guano level (F-ratio = 1.87, p = 0.0266) (Figure 11, Table 3). Additional CCA ordination analyses were done using all 62 species (including the accidental species): (1) with 34 caves and the three different ecological zones within the caves (67 samples all together), and (2) with 34 caves and without cave ecological zones, and found similar results. In analysis (1), precipitation, elevation, latitude, minimum temperature bats, and guano level were significant, as were some of the rock categories. Moreover, in analysis (2), precipitation, latitude, and guano level were significant, but not elevation and minimum temperature (see supporting information Figure S1, Tables S3 and S4).

4. Discussion

4.1. Levantine Cave Arachnofauna

We discovered many unique arthropods, including troglobite [38,39] and troglophile [40,41] arachnids and species endemic to our region. Representatives of the orders Amblypygi, Palpigradi, and Pseudoscorpiones were found only in mesic caves in the Mediterranean region, while Opilioacaridae were found in caves in the semi-desert region. Spiders were the most abundant and diverse arachnid order in the caves studied. Currently, 53 spider families, 297 genera, and 758 species are known from Israel ([48], unpublished data). Of them six troglobite species in three families (Agelenidae (two from the survey and two from caves that were not included in the 2014 survey), Dysderidae, and Leptonetidae), and 34 troglophile species in 11 spider families (Figure 7).

4.2. Troglobites and Troglophiles

Based on our results, 20% of the families known from Israel and Palestine include troglobite and/or troglophile species, compared to 34% in Europe (22 out of 64 families reported from Europe [33,49]). From our current data, only about 5% of the spider species known from Israel and Palestine are cave-dwellers (40 species) with approximately 2% endemics. In 2018, Mammola et al. reported a total of 486 cave-dwelling spider species in Europe (195 troglobite and 291 troglophile species) [33], with 90% of them considered endemics of single countries. In the well-studied cave arachnofauna of Slovenia [50,51], for example, 30% of families (out of 43 families) include troglobite or troglophile species, representing about 11% of species known in Slovenia (N = 753). Some of the families with troglobite and troglophile species in Slovenia [33,50] also occurred in our cave survey in Israel and Palestine (Ageleindae, Dysderidae, Leptonetidae, Lyniphiidae, Nesticidae, Pholcidae, and Theridiidae). Linyphiidae had the greatest number of cave-dwelling spider species in Slovenia, with 48 cave-dwelling species out of 221 linyphiid species in total, followed by Dysderidae with seven cave-dwelling species (22 species in total), Agelenidae with seven cave-dwelling species (25 species in total), and Pholcidae with five cave-dwelling species (five species in total). Our survey uncovered a potentially high species richness of troglophiles and troglobites in the same four families as in Slovenia, but with a different order of richness: Agelenidae, Pholcidae, Linyphiidae, and Dysderidae.
Preference of shaded habitats is a common feature of the four most abundant troglophile spider species in our caves. Placing those four species on a cave-affinity continuum provides a better understanding of their distribution in caves. Loxosceles rufescens is an opportunistic synanthropic species common in the Mediterranean basin in houses and other manmade habitats (and introduced to other areas around the world) [52]. Although it can be found in shaded natural habitats such as caves and under stones [33], and is very common in caves as reported above, it is not a cave specialist. Artema nephilit and F. insidiatrix have a higher affinity to caves than to other habitats (Supplementary Table S2, personal observation). Although A. nephilit can be found also under boulders, in crevices, and in basements, F. insidiatrix can be found in other natural shaded habitats and they each have preference for cave entrances [9,33,49]. Of these four troglophile species, T. pagana has the highest affinity to caves and can be found in cave entrance and twilight zones in large numbers. Similar to some other troglophiles, T. pagana could also be found in other suitable natural shaded habitats, but it is much more common in caves than in any other shaded habitat in Israel and Palestine (Supplementary Table S2, personal observation, Aharon et al., in preparation).

4.3. Foraging Guilds

Among the cave-dwelling spiders we found, web-builders are more species-rich than hunters (27 vs. five species, respectively; Supplementary Table S2). The four most abundant troglophile spider species in our caves represent four out of the six foraging guilds found: sheet-web weavers (T. pagana), sensing-web weavers (F. insidiatrix), ambush hunters (L. rufescens), and space-web weavers (A. nephilit). Higher species richness of web-builders vs. hunters was also found in caves in the Iberian peninsula [53], while ambush hunters and sensing-web weavers were absent from these caves [32,53]. A single species of each of these guilds was abundant in the Levantine caves (L. rufescens and F. insidiatrix, respectively). The troglobite species of the Levantine caves are represented mainly by sheet-web weavers (Agelenidae), one space-web weaver (Leptonetidae), and one hunter species (or species complex) (Dysderidae). Although the spider family Dysderidae has species that specialize on isopods [42], the diet preferences of the dysderids in caves is unknown [53]. Specifically, we lack diet information on the eyeless and eye-bearing Harpactea, eye-bearing Dysdera, and Dasumia crassipalpis (Simon, 1882) found in the Levantine caves, and therefore we assigned them to the guild of ‘other hunters’.

4.4. Species-Pool, Regions, and Bat Inhabitance

The cave-dwelling spiders of the Levant represent a subset of the regional species-pool, with a maximum of 58 species (the 62 that were found in caves without the four troglobite species) that are known or potentially can be found in epigean habitats. Only two spider families, Leptonetidae and Nesticidae, were found only in hypogean habitats in Israel, while 29 spider families were found only in epigean habitats in Israel, and 22 spider families were found both in hypogean and epigean habitats in Israel (including troglobites, troglophiles and accidentals). This distribution is similar to in agreement with findings in other places of the world [32,53].
Caves in the southern region of Israel had relatively low species richness, with no troglobite arachnids and few troglophile spider species. This desert region, covering about 55% of Israel, harbors many epigean spider species that are adapted to the desert climate (at least 246 according to Zonstein and Marusik [48]) and are rarely found in caves. By contrast, all troglobite arachnids and most troglophile spiders were found in caves located in the Mediterranean mesic climate zone, mainly in the north of Israel but also in central Israel and Palestine (Figure 1, Figure 2 and Figure 7A,B). This could be explained by both the higher precipitation (see results and Figure 11) and intensive karstification in the north of Israel [21], i.e., more caves with regional mesic climate may support higher pre-adapted cave-dwelling species. Our analysis found that geographic region is a significant factor effecting both species richness and abundance in caves, and that precipitation and latitude correlated and had a significant effect on the assemblages (Figure 11).
Our preliminary taxonomic study revealed phenotypic variation between populations of different caves in Tegenaria, Harpactea, and the micronetine spiders (Linyphiidae). For some of these spiders, morphological variation was found to be lower within cave populations in comparison to between different caves, and these unique morphological characters are sufficient for describing new species. For other spiders, molecular methods were used to assess cryptic species and possible species complexes. Therefore, our future taxonomic efforts will focus on describing the troglobitic and troglophilic species new to science that have mainly been found in caves in the Mediterranean climate region.
One cave in the arid region, Zavoa cave, hosted an exceptional number of spider species (eight species) and individuals (56 individuals) (Figure 1, Figure 2G and Figure 11, Appendix A). Zavoa cave is situated in the southern Judean desert, with an average of 100 mm rainfall a year, and harbors several cave-dwelling species that may be endemics. One additional cave that was not part of our cave survey, A’rak Na‘asane, in the northern Judean desert, is home to several troglobite and troglophile endemic arachnids [38,54]. Both Zavoa and A’rak Na‘asane are large caves located in a region that experienced aridification during the last glacial interval [55]. Both caves are inhabited by bat colonies. A’rak Na‘asane has a high level of guano, while the Zavoa cave has less guano due to the currently lower numbers of bats. We visited several additional caves in the Judean desert, outside the 2014 survey, and found opilioacriformes and other troglophile arachnids as well as evidence of former use of the caves by bats (dry old guano), but it seems that the bats no longer occupy these caves (unpublished data).
Trajano and de Carvalho (2017) suggested that troglophile species richness could be explained by local ecological factors [18]. Guano is probably an important energy source for cave arthropods and may have a bottom-up effect on the cave food-web, as was shown in caves in Brazil [56,57]. One of our hypotheses was that the level of guano in the cave will positively affect spider richness and abundance, but we did not find this to be significant, but we did find a significant effect of the guano level on the spider assemblage composition in caves (Figure 11, Figure S1, Tables S2 and S3). To test the effect of guano on the Levantine cave spider assemblages in more depth, more information is needed on the ecology of spiders in these caves.

4.5. Levantine Cave Spider Assemblages

We showed that Levantine cave spider assemblages are diverse and are affected by the specific geographic location and its climatic characteristics, as well as by the presence of bats and guano level. While dry caves in the desert climate are dominated by pholcid and theridiid spiders, humid (but not wet) caves in the mesic Mediterranean climate region are dominated by agelenid, filistatid, sicariid, and linyphiid spiders. Dysderid, leptonetid, nesticid, and theraphosid spiders were found only in the mesic Mediterranean climate caves. We could not find effects of the estimated age of the cave, its estimated size, nor the cave ecological zone on spider assemblages. In one analysis (see supporting information), we found a significant effect of the surrounding rock type on the spider assemblage, but most of our caves are located in carbonate rocks (27 of 35) and more samples from basaltic, chalk, marlstone, salt, and sandstone caves are needed in order to have a more balanced analysis (see Appendix A). Additionally, we cannot separate between the geographic location of the cave and its surrounding rock, as some rocks are found solely in one region.

5. Conclusions

To conclude, Levantine cave spider assemblages are diverse, with higher species richness and abundance in caves located in the north of Israel. They tend to have high levels of bat guano and high humidity. We suggest that future research on Levantine cave arachnids focus on caves in the mesic Mediterranean climate region, caves with bats, and caves in arid climate regions that have chambers with high humidity. Our faunistic cave survey is the first and crucial stage in understanding the evolutionary and ecological mechanisms of speciation of arachnids in Levantine caves. Our discoveries contribute to the knowledge of the local arachnofauna in general and are important for conservation of these cave ecosystems. By comparing spider assemblages of Levantine caves to European caves, we identified gaps in the taxonomic research, and expect to focus our efforts on spider families that may have additional cryptic or undescribed cave-dwelling spider species: Agelenidae, Dysderidae, Linyphiidae, and Pholcidae.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/d13050179/s1, Figure S1: CCA ordination graph of the first and second axes testing 34 caves with ecological zones (total 67 samples), and 62 spider species. The significant explanatory variables (precipitation, elevation, latitude, no bats, frugivorous bats and insectivorous bats, low, medium and high guano levels, carbonate and salt rocks) are plotted on the graph, as well as the cave, ecological zone, and species names. Cave sizes are represented by large, medium, and small triangles, respectively; troglobite species by black points; troglophile species by grey points; accidental by empty points; nominal explanatory variable by a black star; continuous explanatory variables by arrows with dashed line, while non-dashed line arrows are used to connect between species name and its centroid; Table S1: LUX values for caves and caves ecological sectors (first visit; second visit). “No“ denote that the ecological sector was missing due to cave size, “Missing” denote a measurement that was not taken. * measurement was taken just before sunset; Table S2: Full list of the spider species with their localities, distribution and guilds; Table S3: CCA ordination results (67 samples), and 62 spider species); Table S4: CCA values for significant environmental variables

Author Contributions

Conceptualization: E.G.-R., S.A. and Y.L.; methodology: E.G.-R., S.A., M.S. and Y.L.; analysis: E.G.-R. and M.S.; investigation: S.A., I.A.S. and E.G.-R.; resources: E.G.-R., Y.L. and M.S.; data curation: E.G.-R. and S.A.; writing—original draft preparation: E.G.-R.; writing—review and editing: E.G.-R., S.A., I.A.S., M.S. and Y.L.; project administration: E.G.-R.; visualization, S.A., I.A.S. and E.G.-R.; supervision, E.G.-R., M.S. and Y.L.; funding acquisition: E.G.-R. and Y.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Israel Taxonomy Initiative (ITI) biodiversity survey grant to E.G.R. and Y.L., and ITI and Ben-Gurion University fellowships to S.A.

Institutional Review Board Statement

Not appliable.

Data Availability Statement

The data presented in this study are available in Appendix A and Supplementary Tables S1 and S2.

Acknowledgments

The research was conducted with collection permits 2013/40027, 2013/40085 and 2014/40313 from the Israel National Parks Authority. We thank J. Cuff for assistance in conducting the cluster analysis in R; M. Isaia and S. Mammola for discussions on the arachnid cave survey; B.A. Huber, I.L.F. Magalhaes, F. Ballarin, M. Rezac, Z. Ganem, and M. Arnedo for assistance in identification and taxonomical collaboration; A. Frumkin, B. Langford and E. Cohen (Israel Cave Research Center) for assistance with cave data and in the field; and E. Lunghi and two anonymous reviewers for comments on an earlier version of the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

DistrictGeographic RegionCaveNEElevationClimate TypePrecipitation (mm)Cave
Minimum Temp.
(E, T, D)
Lithology CategoryLithology Geological AgeCave SizeLength EstimateBatGuano LevelSpider Abundance, Richness
NorthUpper
Galilee
1 Shetula (P)33.087335.3169690Mediterenaean840.411.5, 14.5, 15.5Carbonate rocksCenomanianLarge150InsectivorousLow21, 6
2 Yir’on (cave 1) (P)33.067935.4665528716.416.5, 13, 20EoceneLarge150InsectivorousHigh49, 5
3 Yir’on (cave 2) (P)33.067235.4672541716.48Small10withoutZero26, 8
4 Pelekh (P)32.932435.238488648.512.5TuronianSmall5withoutZero16, 5
5 Yonim (P)32.923635.2168216648.512.5, 12.5Medium50withoutZero61, 14
Golan6 Susita (P)32.779335.657770382.521,21BasaltPlioceneMedium50InsectivorousLow31, 5
Lower
Galilee
7 Berniki (cave 3) (P)32.777535.5401−10245517Carbonate rocksTuronianSmall5withoutZero23, 6
8 Berniki (cave 1) (P)32.777535.5401−10245519.5, 17, 20.5 Large480InsectivorousHigh53, 6
9 Berniki (cave 2) (P)32.776835.5413−16645512.5, 12.5Medium15withoutZero86, 11
Karmel10 Oren (P)32.714434.974973611.216.5, 15.5Medium36withoutZero55, 5
11 Horvat Raqqit (P)32.712835.0123355611.214.5CenomanianSmall3withoutZero8, 5
12 Ezba’ (P)32.711834.9747120611.214.5, 15.5, 18TuronianLarge52FrugivorousHigh58, 7
Central Israel and PalestineNorthern Samaria13 Sal’it (P)32.245435.0456254620.417, 18, 19.5CenomanianLarge108withoutZero26, 6
Central Jordan Valley14 Andartat HaBiqa’ * (E)32.052435.4589−184Hyper-arid19.919.5, 25.5TuronianMedium23.5InsectivorousMedium31, 10
Western
Samaria
15 Oah (P)32.005334.9722123Mediterenaean56916Small10withoutZero19, 7
16 Bet ‘A’rif (P)32.002634.96429556915.5, 18, 18Medium45withoutZero42, 7
17 Tinshemet (P)31.999434.968110056915.5Medium43FrugivorousHigh51, 6
HaShfela18 Haruva (P)31.913334.9607180537.814.5, 16, 18ChalkConiacianCampanianLarge100withoutZero43, 9
Northern Judean Desert19 Perat Southern Slope
(cave 4) * (PE)
31.833435.3054295Semi-arid25015Carbonate rocksTuronianSmall10withoutZero17, 6
20 Perat ‘Inbal (cave 1) * (PE)31.833235.301931425017, 16.5Medium40withoutZero37, 9
21 Perat Ro’im
(cave 2) * (PE)
31.832535.31323825013.5, 16.5Medium25InsectivorousLow45, 9
22 Perat (cave 3) * (PE)31.832135.308326825013Small10withoutZero19, 7
Nortern Dead-Sea Area23 Qumeran * (E)31.755635.459−308Arid95.821, 21CenomanianMedium20withoutZero30, 4
Judean Mountain24 Te’omim (P)31.726235.0217375Mediterenaean50913, 13.5, 14.5Large134FrugivorousHigh45, 5
SouthSouthern Judean
Desert
25 Zavoa’ (PE)31.208635.2311495Arid10022, 23, 23.5Carbonate rocksTuronianLarge600UnknownLow56, 8
Southern Dead-Sea Area26 Arubotayim (E)31.101635.39−348Hyper-arid41.120, 19.5, 19.5SaltPlioceneLarge60withoutZero5, 1
27 Sedom (E)31.087235.3958−38141.120, 19.5, 19.5Large1800withoutZero28, 1
28 Malcham (E)31.076535.3971−38041.120, 19.5, 19.5Large10,000withoutZero11, 1
29 Ne’ot HaKikkar (Nezirim Burial cave) (E)30.991135.3465−33341.120.5MarlstoneQuaternarySmall4withoutZero5, 2
Negev
Desert
30 Telalim (PE)30.973434.7929482Arid91.417, 15Carbonate rocksTuronianMedium27withoutZero7, 2
31 Ashalim (PE)30.943434.739140491.417.5, 19.5, 17Large540InsectivorousLow13, 1
32 Besor (PE)30.941534.696135691.412, 16.5Medium35withoutZero26, 5
33 ‘Avedat (Nezirim cave) (PE)30.794134.77260193.311EoceneSmall5withoutZero9, 2
‘Arava D
esert
34 ‘Ammude ‘Amram (cave 2) (PE)29.651834.9337288Hyper-arid22.521SandstoneCambrianSmall5withoutZero0, 0
35 ‘Ammude ‘Amram (cave 1) (PE)29.651534.933629322.521, 24, 24Large80InsectivorousLow2, 2
List of the 35 caves sampled in this study (north to south) and their environmental variables. Localities in Israel and Palestine (West Bank) and transliterated names of the localities follow the “Israel Touring Map” (1:250,000) and “List of Settlements”, published by the Israel Survey, Ministry of Labor. Geographic coordinates are given in WGS84 (decimal degrees). Cave size estimates were corrected from Mammola et al., 2019 [9]: Bet ‘Arif was changed from large to medium category, Arubotayim was changed from medium to large category, Tinshemet was changed from small to medium, and Andartat HaBiqa’ was added to medium in the current analysis.
(1)
Localities in Palestine (West Bank) are marked by asterisk (*). Letters in parentheses after cave name indicate the zoogeographical region: (P)—Palaearctic; (E)—Ethiopian; (PE)—Palaeoeremic (after Por, 1975).
(2)
Precipitation data is taken from Israel Meteorological Service (https://ims.gov.il/en/ClimateAtlas, 1 February 2021) for the average annual mean for 1980–2010, from the closest meteorological station for each cave.
(3)
Minimum Temperature is the average of the minimum temperature measured by us for each cave for two months (between the first and second visit to each cave) during the survey.
(4)
Length was estimated from cave maps when available, or in field, and represent the distance from the cave opening toward the darker region of the cave that we could reach.

References

  1. Poulson, T.L.; White, W.B. The Cave Environment. Science 1969, 165, 971–981. [Google Scholar] [CrossRef]
  2. Howarth, F.G. Ecology of Cave Arthropods. Annu. Rev. Entomol. 1983, 28, 365–389. [Google Scholar] [CrossRef]
  3. Lauritzen, S.-E. Physiography of the Caves. In Cave Ecology; Moldovan, O.T., Kováč, Ľ., Halse, S., Eds.; Springer International Publishing: Cham, Switzerland, 2018; pp. 7–21. ISBN 978-3-319-98852-8. [Google Scholar]
  4. Zobel, M.; Otto, R.; Laanisto, L.; Naranjo-Cigala, A.; Pärtel, M.; Fernández-Palacios, J.M. The Formation of Species Pools: Historical Habitat Abundance Affects Current Local Diversity: Species Pool Formation and Current Diversity. Global Ecol. Biogeogr. 2011, 20, 251–259. [Google Scholar] [CrossRef]
  5. Trivellone, V.; Bougeard, S.; Giavi, S.; Krebs, P.; Balseiro, D.; Dray, S.; Moretti, M. Factors Shaping Community Assemblages and Species Co-Occurrence of Different Trophic Levels. Ecol. Evol. 2017, 7, 4745–4754. [Google Scholar] [CrossRef] [Green Version]
  6. Willis, K.J.; Whittaker, R.J. Species Diversity: Scale Matters. Science 2002, 295, 1245–1248. [Google Scholar] [CrossRef]
  7. Gibert, J.; Deharveng, L. Subterranean Ecosystems: A Truncated Functional Biodiversity. BioScience 2002, 52, 473. [Google Scholar] [CrossRef]
  8. Mammola, S.; Arnedo, M.A.; Pantini, P.; Piano, E.; Chiappetta, N.; Isaia, M. Ecological Speciation in Darkness? Spatial Niche Partitioning in Sibling Subterranean Spiders (Araneae: Linyphiidae: Troglohyphantes). Invert. Syst. 2018, 32, 1069. [Google Scholar] [CrossRef] [Green Version]
  9. Mammola, S.; Aharon, S.; Seifan, M.; Lubin, Y.; Gavish-Regev, E. Exploring the Interplay Between Local and Regional Drivers of Distribution of a Subterranean Organism. Diversity 2019, 11, 119. [Google Scholar] [CrossRef] [Green Version]
  10. Schneider, K.; Christman, M.C.; Fagan, W.F. The Influence of Resource Subsidies on Cave Invertebrates: Results from an Ecosystem-Level Manipulation Experiment. Ecology 2011, 92, 765–776. [Google Scholar] [CrossRef]
  11. Mammola, S. Finding Answers in the Dark: Caves as Models in Ecology Fifty Years after Poulson and White. Ecography 2019, 42, 1331–1351. [Google Scholar] [CrossRef] [Green Version]
  12. Howarth, F.G. High-Stress Subterranean Habitats and Evolutionary Change in Cave-Inhabiting Arthropods. Am. Nat. 1993, 142, S65–S77. [Google Scholar] [CrossRef]
  13. Manenti, R.; Lunghi, E.; Ficetola, G.F. The Distribution of Cave Twilight-Zone Spiders Depends on Microclimatic Features and Trophic Supply. Invertebr. Biol. 2015, 134, 242–251. [Google Scholar] [CrossRef]
  14. Tobin, B.W.; Hutchins, B.T.; Schwartz, B.F. Spatial and Temporal Changes in Invertebrate Assemblage Structure from the Entrance to Deep-Cave Zone of a Temperate Marble Cave. Int. J. Speleobiol. 2013, 42, 203–221. [Google Scholar] [CrossRef] [Green Version]
  15. Lunghi, E.; Manenti, R.; Ficetola, G.F. Cave Features, Seasonality and Subterranean Distribution of Non-Obligate Cave Dwellers. PeerJ 2017, 5, e3169. [Google Scholar] [CrossRef]
  16. Lunghi, E.; Manenti, R.; Ficetola, G.F. Seasonal Variation in Microhabitat of Salamanders: Environmental Variation or Shift of Habitat Selection? PeerJ 2015, 3, e1122. [Google Scholar] [CrossRef] [Green Version]
  17. Culver, D.C.; Pipan, T. The Biology of Caves and Other Subterranean Habitats, 2nd ed.; Oxford University Press: Oxford, UK; New York, NY, USA, 2019; ISBN 978-0-19-882076-5. [Google Scholar]
  18. Trajano, E.; de Carvalho, M.R. Towards a Biologically Meaningful Classification of Subterranean Organisms: A Critical Analysis of the Schiner-Racovitza System from a Historical Perspective, Difficulties of Its Application and Implications for Conservation. Subterr. Biol. 2017, 22, 1–26. [Google Scholar] [CrossRef]
  19. Avni, Y. Tectonic and Physiographic Settings of the Levant. In Quaternary of the Levant: Environments, Climate Change, and Humans; Bar-Yosef, O., Enzel, Y., Eds.; Cambridge University Press: Cambridge, UK, 2017; pp. 3–16. ISBN 978-1-107-09046-0. [Google Scholar]
  20. Frumkin, A. Quaternary Evolution of Caves and Associated Palaeoenvironments of the Southern Levant. In Quaternary of the Levant: Environments, Climate Change, and Humans; Bar-Yosef, O., Enzel, Y., Eds.; Cambridge University Press: Cambridge, UK, 2017; pp. 135–144. ISBN 978-1-107-09046-0. [Google Scholar]
  21. Frumkin, A. Interaction between Karst, Water and Agriculture over the Climatic Gradient of Israel. Int. J. Speleol. 1999, 28, 99–110. [Google Scholar] [CrossRef]
  22. Audra, P.; Palmer, A. The Pattern of Caves: Controls of Epigenic Speleogenesis. Geomorphol. Relief Process. Environ. 2011, 17, 359–378. [Google Scholar] [CrossRef]
  23. Frumkin, A.; Langford, B.; Lisker, S.; Amrani, A. Hypogenic Karst at the Arabian Platform Margins: Implications for Far-Field Groundwater Systems. GSA Bull. 2017, 129, 1636–1659. [Google Scholar] [CrossRef]
  24. Frumkin, A.; Naor, R. Formation and Modification of Pit Craters–Example from the Golan Volcanic Plateau, Southern Levant. Z. Geomorphol. 2019, 62, 163–181. [Google Scholar] [CrossRef]
  25. Sivan, D.; Sisma-Ventura, G.; Greenbaum, N.; Bialik, O.M.; Williams, F.H.; Tamisiea, M.E.; Rohling, E.J.; Frumkin, A.; Avnaim-Katav, S.; Shtienberg, G.; et al. Eastern Mediterranean Sea Levels through the Last Interglacial from a Coastal-Marine Sequence in Northern Israel. Quat. Sci. Rev. 2016, 145, 204–225. [Google Scholar] [CrossRef]
  26. Frumkin, A.; Fischhendler, I. Morphometry and Distribution of Isolated Caves as a Guide for Phreatic and Confined Paleohydrological Conditions. Geomorphology 2005, 67, 457–471. [Google Scholar] [CrossRef]
  27. Frumkin, A.; Ben-Tov, D.; Buslov, V.; Langford, B.; Valdman, E.; Yaaran, S. The Southernmost Authigenic Karst Plateau of the Levant. In Proceedings of the 17th International Congress of Speleology, Sydney, Australia, 23–29 July 2017; pp. 13–16. [Google Scholar]
  28. Frumkin, A.; Langford, B.; Porrat, R. The Judean Desert–The major hypogene cave region of the southern Levant. In Hypogene Karst Regions and Caves of the World; Springer: Berlin/Heidelberg, Germany, 2017; pp. 463–477. [Google Scholar]
  29. Frumkin, A. Morphology and Development of Salt Caves. Natl. Speleol. Soc. Bull. 1994, 56, 82–95. [Google Scholar]
  30. Frumkin, A.; Bar-Matthews, M.; Davidovich, U.; Langford, B.; Porat, R.; Ullman, M.; Zissu, B. In-Situ Dating of Ancient Quarries and the Source of Flowstone (‘Calcite-Alabaster’) Artifacts in the Southern Levant. J. Archaeol. Sci. 2014, 41, 749–758. [Google Scholar] [CrossRef]
  31. Resende, L.P.A.; Bichuette, M.E. Sharing the Space: Coexistence among Terrestrial Predators in Neotropical Caves. J. Nat. Hist. 2016, 50, 2107–2128. [Google Scholar] [CrossRef]
  32. Mammola, S.; Isaia, M. Spiders in Caves. Proc. R. Soc. B. 2017, 284. [Google Scholar] [CrossRef] [Green Version]
  33. Mammola, S.; Cardoso, P.; Ribera, C.; Pavlek, M.; Isaia, M. A Synthesis on Cave-Dwelling Spiders in Europe. J. Zool. Syst. Evol. Res. 2018, 56, 301–316. [Google Scholar] [CrossRef]
  34. Parimuchová, A.; Dušátková, L.P.; Kováč, Ľ.; Macháčková, T.; Slabý, O.; Pekár, S. The Food Web in a Subterranean Ecosystem Is Driven by Intraguild Predation. Sci. Rep. 2021, 11, 4994. [Google Scholar] [CrossRef]
  35. Reddell, J.R. Spiders and Related Groups. In Encyclopedia Caves; Culver, D.C., White, W.B., Eds.; Elsevier Academic Press: Burlington, MA, USA, 2005; pp. 554–564. [Google Scholar]
  36. Levy, G. The First Troglobite Scorpion from Israel and a New Chactoid Family (Arachnida: Scorpiones). Zool. Middle East 2007, 40, 91–96. [Google Scholar] [CrossRef]
  37. Karaman, I. A New Euscorpius Species (Scorpiones: Euscorpiidae) from a Dinaric Cave–the First Record of Troglobite Scorpion in European Fauna. Biol. Serbica 2020, 42, 14–31. [Google Scholar] [CrossRef]
  38. Aharon, S.; Ballesteros, J.A.; Crawford, A.R.; Friske, K.; Gainett, G.; Langford, B.; Santibáñez-López, C.E.; Ya’aran, S.; Gavish-Regev, E.; Sharma, P.P. The Anatomy of an Unstable Node: A Levantine Relict Precipitates Phylogenomic Dissolution of Higher-Level Relationships of the Armoured Harvestmen (Arachnida: Opiliones: Laniatores). Invert. Syst. 2019. [Google Scholar] [CrossRef]
  39. Miranda, G.S.; Aharon, S.; Gavish-Regev, E.; Giupponi, A.P.L.; Wizen, G. A New Species of Charinus Simon, 1892 (Arachnida: Amblypygi: Charinidae) from Israel and New Records of C. Ioanniticus (Kritscher, 1959). Eur. J. Taxon. 2016. [Google Scholar] [CrossRef] [Green Version]
  40. Aharon, S.; Huber, B.A.; Gavish-Regev, E. Daddy-Long-Leg Giants: Revision of the Spider Genus Artema Walckenaer, 1837 (Araneae, Pholcidae). Eur. J. Taxon. 2017. [Google Scholar] [CrossRef] [Green Version]
  41. Gavish-Regev, E.; Aharon, S.; Armiach, I.; Lubin, Y. Cave Survey Yields a New Spider Family Record for Israel. Aramit 2016, 51, 39–42. [Google Scholar] [CrossRef]
  42. Cardoso, P.; Pekár, S.; Jocqué, R.; Coddington, J.A. Global Patterns of Guild Composition and Functional Diversity of Spiders. PLoS ONE 2011, 6, e21710. [Google Scholar] [CrossRef] [Green Version]
  43. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2021. [Google Scholar]
  44. Ter Braak, C.; Šmilauer, P. CANOCO Reference Manual and User’s Guide to Canoco for Windows. Software for Canonical Community Ordination (Version 4.5); Centre for Biometry: Wageningen, The Netherlands, 2002. [Google Scholar]
  45. Lepš, J.; Šmilauer, P. Multivariate Analysis of Ecological Data Using CANOCO, 1st ed.; Cambridge University Press: Cambridge, UK, 2003; ISBN 978-0-521-81409-6. [Google Scholar]
  46. Warnes, G.R.; Bolker, B.; Bonebakker, L.; Gentleman, R.; Huber, W.; Liaw, A.; Lumley, T.; Maechler, M.; Magnusson, A.; Moeller, M.; et al. Gplots: Various R Programming Tools for Plotting Data. R Package. 2020. Available online: https://github.com/talgalili/gplots (accessed on 18 February 2021).
  47. Ballesteros, J.A.; Santibáñez López, C.E.; Kováč, Ľ.; Gavish-Regev, E.; Sharma, P.P. Ordered Phylogenomic Subsampling Enables Diagnosis of Systematic Errors in the Placement of the Enigmatic Arachnid Order Palpigradi. Proc. R. Soc. B. 2019, 286. [Google Scholar] [CrossRef] [PubMed]
  48. Zonstein, S.; Marusik, Y.M. Checklist of the Spiders (Araneae) of Israel. Zootaxa 2013, 3671, 1. [Google Scholar] [CrossRef]
  49. Nentwig, W.; Blick, T.; Bosmans, R.; Gloor, D.; Hänggi, A.; Kropf, C. Spiders of Europe. Version 02.2021. Available online: https://araneae.nmbe.ch/ (accessed on 18 February 2021).
  50. Kostanjšek, R.; Kuntner, M. BioPortal–Katalog–Araneae Sloveniae. Available online: http://www.bioportal.si/katalog/araneae.php (accessed on 18 February 2021).
  51. Kostanjšek, R.; Kuntner, M. Araneae Sloveniae: A National Spider Species Checklist. ZooKeys 2015, 474, 1–91. [Google Scholar] [CrossRef] [Green Version]
  52. Massa, M.; Planas, E.; Ribera, C. The Mediterranean as a Melting Pot: Phylogeography of Loxosceles Rufescens (Sicariidae) in the Mediterranean Basin. PLoS ONE 2018, 13, e0210093. [Google Scholar] [CrossRef]
  53. Cardoso, P. Diversity and Community Assembly Patterns of Epigean vs. Troglobiont Spiders in the Iberian Peninsula. Int. J. Speleol. 2012, 41, 83–94. [Google Scholar] [CrossRef]
  54. Frumkin, A.; Aharon, S.; Davidovich, U.; Langford, B.; Negev, Y.; Ullman, M.; Vaks, A.; Ya‘aran, S.; Zissu, B. Old and Recent Processes in a Warm and Humid Desert Hypogene Cave: ‘A’Rak Na‘asane, Israel. Int. J. Speleol. 2018, 47, 307–321. [Google Scholar] [CrossRef]
  55. Vaks, A.; Bar-Matthews, M.; Ayalon, A.; Schilman, B.; Gilmour, M.; Hawkesworth, C.J.; Frumkin, A.; Kaufman, A.; Matthews, A. Paleoclimate Reconstruction Based on the Timing of Speleothem Growth and Oxygen and Carbon Isotope Composition in a Cave Located in the Rain Shadow in Israel. Quat. Res. 2003, 59, 182–193. [Google Scholar] [CrossRef]
  56. Ferreira, R.L.; Martins, R.P. Trophic Structure and Natural History of Bat Guano Invertebrate Communities, with Special Reference to Brazilian Caves. Trop. Zool. 1999, 12, 231–252. [Google Scholar] [CrossRef]
  57. Ferreira, R.; Martins, R.; Prous, X. Structure of Bat Guano Communities in a Dry Brazilian Cave. Trop. Zool. 2007, 20, 55–74. [Google Scholar]
Figure 1. Distribution of the 35 caves. Large caves are represented by pink dots, medium caves by blue dots, and small caves by green dots. Numbers denote a specific cave (see Appendix A).
Figure 1. Distribution of the 35 caves. Large caves are represented by pink dots, medium caves by blue dots, and small caves by green dots. Numbers denote a specific cave (see Appendix A).
Diversity 13 00179 g001
Figure 2. Caves and habitats from the survey (cave numbers in parentheses). (AC): Northern Israel, Mediterranean climate. (A): The area near the opening of Yonim cave (upper Galilee, 5). (B): The north-facing slope of Oren wadi where Ezba’ cave (12) is located (taken from Oren cave, Karmel). (C): The south-facing slope of Oren wadi where Oren cave (Karmel, 10) is located. (DF): Central Israel and Palestine, Mediterranean and semi-arid climates. (D): The area near the opening of Teomim cave (Judean Mountains, 24). (E): The area near the opening of Haruva cave (HaShfela, 18). (F): The north-facing slope of Perat wadi where ‘Perat Inbal cave (cave 1, northern Judean desert, Palestine, 20) is located. (GI): Southern Israel, arid climates. (G): The east-facing slope where Zavoa cave is located (southern Judean desert, 25). (H): The east-facing slope where Arubotayim cave is located (southern Dead-sea region, 26). (I): The area where Ammude Amram caves are located (Arava valley, 34–35). Pictures: (AC,FI) Shlomi Aharon, (D,E) Igor Armiach Steinpress.
Figure 2. Caves and habitats from the survey (cave numbers in parentheses). (AC): Northern Israel, Mediterranean climate. (A): The area near the opening of Yonim cave (upper Galilee, 5). (B): The north-facing slope of Oren wadi where Ezba’ cave (12) is located (taken from Oren cave, Karmel). (C): The south-facing slope of Oren wadi where Oren cave (Karmel, 10) is located. (DF): Central Israel and Palestine, Mediterranean and semi-arid climates. (D): The area near the opening of Teomim cave (Judean Mountains, 24). (E): The area near the opening of Haruva cave (HaShfela, 18). (F): The north-facing slope of Perat wadi where ‘Perat Inbal cave (cave 1, northern Judean desert, Palestine, 20) is located. (GI): Southern Israel, arid climates. (G): The east-facing slope where Zavoa cave is located (southern Judean desert, 25). (H): The east-facing slope where Arubotayim cave is located (southern Dead-sea region, 26). (I): The area where Ammude Amram caves are located (Arava valley, 34–35). Pictures: (AC,FI) Shlomi Aharon, (D,E) Igor Armiach Steinpress.
Diversity 13 00179 g002
Figure 3. Arachnids found during cave surveys. (A): Argas vespertilionis (Latreille, 1796), (Acari, Argasidae); (B): Ornithodoros tholozani Laboulbène and Mégnin, 1882, (Acari, Argasidae); (C,D): Pseudoscorpiones; (E): Charinus israelensis Miranda, Aharon, Gavish-Regev, Giupponi, and Wizen, 2016, (Amblypygi, Charinidae); (F): Haasus naasane Aharon et al., 2019 (Opiliones, Pyramidopidae); (G): Opilioacaridae (Acari); (H): Eukoenenia Borner, 1901 (Palpigradi, Eukoeneniidae). Pictures: Shlomi Aharon.
Figure 3. Arachnids found during cave surveys. (A): Argas vespertilionis (Latreille, 1796), (Acari, Argasidae); (B): Ornithodoros tholozani Laboulbène and Mégnin, 1882, (Acari, Argasidae); (C,D): Pseudoscorpiones; (E): Charinus israelensis Miranda, Aharon, Gavish-Regev, Giupponi, and Wizen, 2016, (Amblypygi, Charinidae); (F): Haasus naasane Aharon et al., 2019 (Opiliones, Pyramidopidae); (G): Opilioacaridae (Acari); (H): Eukoenenia Borner, 1901 (Palpigradi, Eukoeneniidae). Pictures: Shlomi Aharon.
Diversity 13 00179 g003
Figure 4. Rank-abundance graph of 62 spider species from 34 Levantine caves (‘Ammude ‘Amram small cave (34) had no observed spiders during our two visits).
Figure 4. Rank-abundance graph of 62 spider species from 34 Levantine caves (‘Ammude ‘Amram small cave (34) had no observed spiders during our two visits).
Diversity 13 00179 g004
Figure 5. Box plot of (A) species richness and (B) abundance of troglobite, troglophile, and accidental spiders in 35 caves at the north, center, and south of Israel (minimum, first quartile, median, third quartile, and maximum).
Figure 5. Box plot of (A) species richness and (B) abundance of troglobite, troglophile, and accidental spiders in 35 caves at the north, center, and south of Israel (minimum, first quartile, median, third quartile, and maximum).
Diversity 13 00179 g005
Figure 6. Troglobite spiders found during cave surveys (cave numbers in parentheses). (A,B): Tegenaria Latreille, 1804 (Agelenidae). (A): Tegenaria from Yir’on large cave (upper Galilee, 2); (B): Tegenaria from Teomim large cave (Judean mountains, 24); (C) Harpactea Bristowe, 1939 (Dysderidae) eyeless species from Shetula large cave (upper Galilee, 1); (D): Cataleptoneta edentula Denis, 1955 (Leptonetidae) from Ezba’ large cave (Karmel ridge, 12). Pictures: Shlomi Aharon.
Figure 6. Troglobite spiders found during cave surveys (cave numbers in parentheses). (A,B): Tegenaria Latreille, 1804 (Agelenidae). (A): Tegenaria from Yir’on large cave (upper Galilee, 2); (B): Tegenaria from Teomim large cave (Judean mountains, 24); (C) Harpactea Bristowe, 1939 (Dysderidae) eyeless species from Shetula large cave (upper Galilee, 1); (D): Cataleptoneta edentula Denis, 1955 (Leptonetidae) from Ezba’ large cave (Karmel ridge, 12). Pictures: Shlomi Aharon.
Diversity 13 00179 g006
Figure 7. Proportion and number of troglobite, troglophile, and epigean spider species per family in Israel and Palestine (53 families and 758 species including undescribed species, description in progress). Families that have only epigean species were pooled in one column.
Figure 7. Proportion and number of troglobite, troglophile, and epigean spider species per family in Israel and Palestine (53 families and 758 species including undescribed species, description in progress). Families that have only epigean species were pooled in one column.
Diversity 13 00179 g007
Figure 8. Eight of the nine most abundant troglophile spiders found during cave surveys (cave numbers in parentheses). (AD): Pholcidae: (A): Artema nephilit Aharon, Huber and Gavish-Regev 2017 from Oren medium cave (Karmel ridge, 10); (B): Hoplopholcus cecconii Kulczyński, 1908 from Yir’on large cave (upper Galilee, 2); (C): Holocnemus pluchei Scopoli, 1763 from a cave in Perat wadi (northern Judean desert, 19–22); (D): Pholcus Walckenaer, 1805 from ‘Inbal cave 1 in Perat wadi (northern Judean desert); (E): Tegenaria pagana C. L. Koch, 1840 (Agelenidae) from Teomim cave (Judean mountains); (F): Filistata insidiatrix (Forskål, 1775) from ‘Inbal cave 1 in Perat wadi (northern Judean desert, 20); (G): Loxosceles rufescens (Dufour, 1820) (Sicaridae) from Berniki caves (lower Galilee, 7–9); (H): Steatoda triangulosa (Walckenaer, 1802) (Theridiidae) from Tinshemet cave (western Samaria, 17). Pictures: (AE,G,H): Shlomi Aharon, (F): Igor Armiach Steinptress.
Figure 8. Eight of the nine most abundant troglophile spiders found during cave surveys (cave numbers in parentheses). (AD): Pholcidae: (A): Artema nephilit Aharon, Huber and Gavish-Regev 2017 from Oren medium cave (Karmel ridge, 10); (B): Hoplopholcus cecconii Kulczyński, 1908 from Yir’on large cave (upper Galilee, 2); (C): Holocnemus pluchei Scopoli, 1763 from a cave in Perat wadi (northern Judean desert, 19–22); (D): Pholcus Walckenaer, 1805 from ‘Inbal cave 1 in Perat wadi (northern Judean desert); (E): Tegenaria pagana C. L. Koch, 1840 (Agelenidae) from Teomim cave (Judean mountains); (F): Filistata insidiatrix (Forskål, 1775) from ‘Inbal cave 1 in Perat wadi (northern Judean desert, 20); (G): Loxosceles rufescens (Dufour, 1820) (Sicaridae) from Berniki caves (lower Galilee, 7–9); (H): Steatoda triangulosa (Walckenaer, 1802) (Theridiidae) from Tinshemet cave (western Samaria, 17). Pictures: (AE,G,H): Shlomi Aharon, (F): Igor Armiach Steinptress.
Diversity 13 00179 g008
Figure 9. Additional troglophile spiders found during cave surveys. (AD): Linyphiidae, Micronetinae: (A): Male micronetine sp. 1 form Besor medium cave (Negev desert, 32); (B): Male micronetine sp. 6 from Teomim large cave (Judean mountains, 24); (C): Male micronetine sp. 1? from Bet ‘Arif medium cave (western Samaria, 16); (D): Male micronetine from ‘Avedat cave (Negev desert, 33); (E): Chaetopelma Ausserer, 1871 (Theraphosidae) from Modi’in cave (Judea, not part of the survey); (F): Phyxelida anatolica Griswold, 1990 (Phyxelididae) from Haruva cave (HaShfela, 18); (G): Heteropoda variegata (Simon, 1874) (Sparassidae) from Yir’on large cave (upper Galilee, 2); (H): Harpactea Bristowe, 1939 (Dysderidae) from Ezba’ large cave (Karmel ridge, 12). Pictures: (AD,FH): Shlomi Aharon, (E): Igor Armiach Steinptress.
Figure 9. Additional troglophile spiders found during cave surveys. (AD): Linyphiidae, Micronetinae: (A): Male micronetine sp. 1 form Besor medium cave (Negev desert, 32); (B): Male micronetine sp. 6 from Teomim large cave (Judean mountains, 24); (C): Male micronetine sp. 1? from Bet ‘Arif medium cave (western Samaria, 16); (D): Male micronetine from ‘Avedat cave (Negev desert, 33); (E): Chaetopelma Ausserer, 1871 (Theraphosidae) from Modi’in cave (Judea, not part of the survey); (F): Phyxelida anatolica Griswold, 1990 (Phyxelididae) from Haruva cave (HaShfela, 18); (G): Heteropoda variegata (Simon, 1874) (Sparassidae) from Yir’on large cave (upper Galilee, 2); (H): Harpactea Bristowe, 1939 (Dysderidae) from Ezba’ large cave (Karmel ridge, 12). Pictures: (AD,FH): Shlomi Aharon, (E): Igor Armiach Steinptress.
Diversity 13 00179 g009
Figure 10. Cluster heatmap for 34 caves (rows) and 32 troglobite and troglophile species (columns) (‘Ammude’ Amram (cave 34) had no observed spiders during our two visits). Caves are clustered using hierarchical clustering from pairwise Pearson correlation. Darker colors signify higher correlation. Numbers before cave names as in Figure 1 and Appendix A.
Figure 10. Cluster heatmap for 34 caves (rows) and 32 troglobite and troglophile species (columns) (‘Ammude’ Amram (cave 34) had no observed spiders during our two visits). Caves are clustered using hierarchical clustering from pairwise Pearson correlation. Darker colors signify higher correlation. Numbers before cave names as in Figure 1 and Appendix A.
Diversity 13 00179 g010
Figure 11. CCA ordination graph of the first and second axes testing 34 caves and 32 troglobite/troglophile species. The significant explanatory variables (precipitation, elevation, latitude, minimum temperature, and high guano level) are plotted on the graph, as well as the cave and species names. Cave sizes are represented by large, medium, and small triangles, respectively, troglobite species by black points, troglophile species by grey points, nominal explanatory variable by a black star, continuous explanatory variables by arrows with the dashed line. The non-dashed line arrows are used to connect between species name and its centroid.
Figure 11. CCA ordination graph of the first and second axes testing 34 caves and 32 troglobite/troglophile species. The significant explanatory variables (precipitation, elevation, latitude, minimum temperature, and high guano level) are plotted on the graph, as well as the cave and species names. Cave sizes are represented by large, medium, and small triangles, respectively, troglobite species by black points, troglophile species by grey points, nominal explanatory variable by a black star, continuous explanatory variables by arrows with the dashed line. The non-dashed line arrows are used to connect between species name and its centroid.
Diversity 13 00179 g011
Table 1. Spider species richness and abundance (troglobite, troglophile, and accidentals) in caves by region (35 caves). Numbers denoted a specific cave (see Appendix A).
Table 1. Spider species richness and abundance (troglobite, troglophile, and accidentals) in caves by region (35 caves). Numbers denoted a specific cave (see Appendix A).
South (11 Caves)Center (12 Caves)North (12 Caves)
Spider species richness (minimum, median, maximum; mean)0 1, 2, 8 2; 2.274 3, 7, 10 4; 7.085 5, 6, 14 6; 6.91
Spider abundance (minimum, median, maximum; mean)0 1, 9, 56 2; 14.7317 7, 34, 51 8; 33.758 5, 40, 86 9; 40.58
1 ‘Ammude ‘Amram (small cave, 34); 2 Zavoa’ (25); 3 Qumeran (23); 4 Andartat HaBiqa’ (14); 5 Raqqit (11); 6 Yonim (5); 7 Perat cave # 4 (19); 8 Tinshemet (17); 9 Berniki (medium cave, 9).
Table 2. GLM analyses of the effect of environmental variables on the probability to detect specific spider guilds in the caves. Wald χ2 and p values are shown for the two models. Neither geographic region nor bat presence (or guano level) was a significant explanatory variable for any of the guilds.
Table 2. GLM analyses of the effect of environmental variables on the probability to detect specific spider guilds in the caves. Wald χ2 and p values are shown for the two models. Neither geographic region nor bat presence (or guano level) was a significant explanatory variable for any of the guilds.
Ambush HuntersSpecialists and Other HuntersSensing WebSheet or Space WebOrb Web
DfWald χ 2 p-ValueWald χ 2 p-ValueWald χ 2 p-ValueWald χ 2 p-ValueWald χ 2 p-Value
Model 1Geographic region22.1440.3422.9630.2270.7060.7034.4400.1091.7330.420
Bat presence10.0240.8770.3940.5300.0180.8920.0650.7980.0010.979
Model 2Geographic region22.0470.3594.1460.1260.3140.8553.7760.1511.7290.421
Guano amount30.5860.9003.0210.3880.6090.8940.2770.9640.3390.953
Table 3. CCA ordination results.
Table 3. CCA ordination results.
EigenvaluesSpecies-Environment
Correlations
Cumulative Percentage Variance Sum of All
Eigenvalues
Sum of All
Canonical
Eigenvalues
of Species Dataof Species—Environment Relation
Axis 10.6420.9178.9317.1952.070
Axis 20.5890.83117.159.5
Axis 30.4480.83923.381.1
Axis 40.2140.72526.391.4
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Gavish-Regev, E.; Aharon, S.; Armiach Steinpress, I.; Seifan, M.; Lubin, Y. A Primer on Spider Assemblages in Levantine Caves: The Neglected Subterranean Habitats of the Levant—A Biodiversity Mine. Diversity 2021, 13, 179. https://doi.org/10.3390/d13050179

AMA Style

Gavish-Regev E, Aharon S, Armiach Steinpress I, Seifan M, Lubin Y. A Primer on Spider Assemblages in Levantine Caves: The Neglected Subterranean Habitats of the Levant—A Biodiversity Mine. Diversity. 2021; 13(5):179. https://doi.org/10.3390/d13050179

Chicago/Turabian Style

Gavish-Regev, Efrat, Shlomi Aharon, Igor Armiach Steinpress, Merav Seifan, and Yael Lubin. 2021. "A Primer on Spider Assemblages in Levantine Caves: The Neglected Subterranean Habitats of the Levant—A Biodiversity Mine" Diversity 13, no. 5: 179. https://doi.org/10.3390/d13050179

APA Style

Gavish-Regev, E., Aharon, S., Armiach Steinpress, I., Seifan, M., & Lubin, Y. (2021). A Primer on Spider Assemblages in Levantine Caves: The Neglected Subterranean Habitats of the Levant—A Biodiversity Mine. Diversity, 13(5), 179. https://doi.org/10.3390/d13050179

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

Article Metrics

Back to TopTop