Ground-dwelling spider assemblages in contrasting habitats in the central South African Grassland Biome

In Africa, the Grassland Biome is largely limited to the central plateau of South Africa, Lesotho and parts of Swaziland (Mucina & Rutherford 2006). It is characterised by extremely high plant biodiversity, second only to that of the Fynbos Biome (Low & Rebelo 1996). Grasslands can be defined as a single-layered herbaceous plant community, with a few woody plants, which are usually restricted to specific habitats, including drainage lines and rocky hilltops (Carbutt et al. 2011). The Grassland Biome is one of the most transformed biomes in South Africa and is under continuous threat from cultivation, overgrazing and urban expansion (Bredenkamp & Van Rooyen 1996). Only an estimated 2.04% to 2.80% of this biome is formally conserved (Carbutt et al. 2011; O’Connor & Kuyler 2005), and therefore, effective management and conservation of private land is critical to protect its highly endemic fauna and flora (Wessels et al. 2003).


Introduction
In Africa, the Grassland Biome is largely limited to the central plateau of South Africa, Lesotho and parts of Swaziland (Mucina & Rutherford 2006).It is characterised by extremely high plant biodiversity, second only to that of the Fynbos Biome (Low & Rebelo 1996).Grasslands can be defined as a single-layered herbaceous plant community, with a few woody plants, which are usually restricted to specific habitats, including drainage lines and rocky hilltops (Carbutt et al. 2011).The Grassland Biome is one of the most transformed biomes in South Africa and is under continuous threat from cultivation, overgrazing and urban expansion (Bredenkamp & Van Rooyen 1996).Lotz, Seaman & Kok 1991;Van den Berg & Dippenaar-Schoeman 1991), plant-dwelling (Dippenaar-Schoeman, Hamer & Haddad 2011;Fourie et al. 2013;Haddad 2005;Neethling & Haddad 2013), litter-dwelling (Butler & Haddad 2011) and termitophilous assemblages (Haddad & Dippenaar-Schoeman 2002, 2006).Different species are adapted to particular microhabitats within grasslands, either living on the ground, on grasses or on foliage of woody plants, with few species abundant in multiple strata (Haddad et al. 2013).Consequently, there is considerable scope for research on spider biodiversity, ecology and biology in this unique biome.
Only three of the aforementioned studies have investigated finer scale differences in spider assemblages associated with different plant species, notably litter-dwellers (Butler & Haddad 2011) and foliage-dwellers (Fourie et al. 2013;Neethling & Haddad 2013).However, assemblage structure in contrasting vegetation communities has so far only been investigated for grass-dwellers in structurally variable grasslands (Fourie et al. 2013).
All the pitfall trapping surveys listed above were conducted in open grasslands with sparse or absent woody vegetation, and spider assemblages in shrublands and woodlands remain largely unknown in this biome (Butler & Haddad 2011).The aims of this study were (1) to sample grounddwelling spider assemblages in different plant communities (predominantly shrublands) to determine possible habitat associations of spider species; (2) assess how habitats affected the activity density and species richness of spiders; (3) determine seasonality of ground-dwelling spider assemblages, and (4) determine whether indicator species could be identified for any of the sampled habitats.Further, this study aims to add to the current knowledge base of ground-dwelling spider biodiversity in South Africa, for which relatively little information is currently available (Dippenaar-Schoeman et al. 2015;Janion-Scheepers et al. 2016).

Study area
Amanzi Private Game Reserve (APGR) is located about 80 km north-east of Bloemfontein in the central Free State (Figure 1a) and falls within the summer rainfall region of central South Africa, with an average of approximately 475 mm of rainfall received annually (Butler 2017).The area experiences hot summers, with day temperatures sometimes between 35 °C and 40 °C (averaging above 30 °C), and cold winters, with night temperatures frequently below freezing and day temperatures usually ranging between 15 °C and 20 °C (Butler 2017).
The study area is located in the Grassland Biome (Rutherford & Westfall 1994), with the vegetation of the surrounding area being described as Dry Sandy Highveld Grassland (Bredenkamp & Van Rooyen 1996) or Vaal-Vet Sandy Grassland (Mucina & Rutherford 2006).The vegetation in APGR is, however, more representative of Winburg Grassy Shrubland, which occurs in a series of larger patches from Trompsburg through Bloemfontein and Winburg to Ventersburg (Mucina & Rutherford 2006).
The landscape of this vegetation type consists of isolated hills, slopes and escarpments, creating habitats ranging from open grassland to shrubland (Mucina & Rutherford 2006).
A comprehensive vegetation survey was conducted in the western section of APGR, and details of the community composition in this part of the reserve are provided in Butler (2017).
Based on the results of the vegetation survey, spiders were collected in eight different habitats representing four plant communities (Table 1; Figure 1b).In the Buddelja saligna-Searsia burchellii community, two different subcommunities were sampled, one of which included three habitats associated with the northern and southern slopes and plateau of a hill (Buddleja saligna-Searsia burchellii-Olea europaea africana subcommunity), and the second including two habitats, one dominated by Searsia burchellii and the other by Tarchonanthus camphoratus (Buddleja saligna-Searsia burchellii-Vachellia karroo subcommunity).The remaining three habitats sampled were the Themeda triandra-Digitaria eriantha community, which was dominated in the woody layer by Vachellia karroo; cultivated D. eriantha pastures (Digitaria eriantha-Cynodon dactylon community); and an area in and around a freshwater dam (Persicaria lapathifolia-Panicum coloratum community).The latter two communities had little or no plants contributing to the woody layer (Table 1).

Spider sampling
Two sites were sampled in each habitat, with at least 100 m separating them to avoid pseudoreplication.Five pitfall traps were placed 5 m apart in a straight line at each site.A soil auger was used for drilling holes, and plastic buckets 10 cm in diameter were used as pitfall traps.Ethylene glycol (100 mL) was added to each pitfall trap to preserve terrestrial arthropods for later identification.Collected material was removed and the pitfall traps refilled at the end of each month, from the start of September 2012 (early spring) to the end of August 2013 (end of winter).The sampled material was sorted in the laboratory and all arachnids were extracted from the samples and preserved in 70% ethanol.Following identification and tallying of adult spiders, material was deposited in the National Collection of Arachnida at the ARC-Plant Protection Research in Pretoria, South Africa.
Although spiders were sampled for a full year, several sampling months represented incomplete sampling efforts as a result of flooding of the pitfall traps at some sites following heavy rains, or damage caused to pitfall traps by large herbivorous mammals or vervet monkeys.To ensure that sampling effort was equal between habitats, we opted to provide data for 1 month in each season for which all sites had a complete sampling effort: early spring (September 2012), mid-summer (January 2013), mid-autumn (April 2013) and mid-winter (July 2013).Only adult spiders were included in the analysis.

Statistical analysis
We calculated the estimated species richness for each habitat using the equation: where the number of species represented by a single individual (i.e.singletons) and two individuals only (i.e.doubletons) are represented by F 1 and F 2 , respectively (Chao 1984).Chao1 is an estimator calculated using the available abundance data and is a function of the ratio between the singletons and doubletons in the data.With an increase in the number of samples, an accumulation curve reaches an asymptote when all species in the community are represented by at least two individuals.
We calculated sampling completeness as the ratio of the observed species richness (S obs ) and the Chao1-estimated species richness (S Chao1 ) (e.g.Cardoso et al. 2008;Scharff et al. 2003;Sørensen, Coddington & Scharff 2002).The herbaceous layer of this plant community was dominated by Themeda triandra, Sporobolus fimbriatus and Digitaria eriantha, with the woody layer dominated by the deciduous tree Vachellia karroo Chao and Jost (2012) proposed the use of coverage-based rarefaction and extrapolation to assess community richness and sampling effort.They defined sample coverage as the proportion of the total number of individuals in a community that belong to the species represented in the sample.Subtracting the sample coverage from unity gives the proportion of the community belonging to as yet unsampled species, which they referred to as the 'coverage deficit'.This can be inferred as the likelihood that a new, previously unsampled species will be found if the sample was increased by one individual (Chao & Jost 2012).Coverage for each habitat and for the total spider assemblage was calculated using the following equation: where n represents the number of individuals in the sample and f 1 and f 2 represent the number of singleton and doubleton species, respectively.Chao and Lee (1992) proposed that an estimated coverage value should be at least 50%, that is, 0.5.

Inventory completeness was analysed in Paleontological
Statistics (PAST) version 2.07 using the sample rarefaction function, which implements the 'Mao tau' analytical procedure, with standard errors indicated as 95% confidence intervals on the resulting graphs.We produced curves for each of the habitats, as well as the for the whole spider assemblage.
Using PAST versions 2.07 and 3.06 (Hammer, Harper & Ryan 2001), we calculated whether the abundance and species richness of spiders differed between the eight habitats and also whether these parameters varied seasonally, using linear multivariate regression.In each analysis, habitats and samples were used as independent variables.Because count data follows a Poisson distribution, values were log-transformed prior to analysis to approach a normal distribution.
We then determined whether the two sites from each habitat sampled similar spider assemblages by performing a twodimensional non-metric multidimensional scaling (NMDS), based on a Bray-Curtis similarity matrix.Ideally, the stress value of an NMDS should be lower than 0.2; otherwise, the resulting diagram needs to be interpreted with caution (Clarke 1993).Further, we performed a cluster analysis using the cluster function and unweighted pair-group average algorithm (Clarke & Warwick 2001).In both analyses, the pooled data from each of the 16 sampling sites were used.Further, a permutational multivariate analysis of variance (PERMANOVA) was performed to test for differences in assemblages between the eight habitats, using the Bray-Curtis distance measure and 10 000 permutations.These analyses were carried out in PAST version 2.07.
Lastly, we attempted to identify indicator spider species, which are considered to be characteristic of a particular habitat.Indicator values were obtained by multiplying a species' relative abundance in a particular habitat, expressed as a percentage of its total abundance, with its relative frequency of occurrence in that particular habitat, that is, proportion of samples in which a species was collected (Dufrene & Legendre 1997).Thus, a species' specificity (narrow association with a particular habitat) and fidelity (frequency of occurrence in that habitat) is expressed as a percentage that can be compared with other species in the sampled habitats (Dufrene & Legendre 1997).A high indicator value illustrates a high affiliation of a species to a particular habitat, with a suitable benchmark of 70% being suggested (e.g.Haddad et al. 2010;McGeoch, Van Rensburg & Botes 2002;Van Rensburg et al. 1999).

Ethical consideration
Permission to collect arachnids in the Free State province was obtained from the Free State Department of Economic Development, Tourism and Environmental Affairs.

Results
A total of 2982 adult spiders were collected, representing 129 species and 33 families (Table 2; Appendix 1).Ammoxenidae was the dominant family (n = 1218, 40.85%), followed by Gnaphosidae (n = 634, 21.26%), Zodariidae (n = 322, 10.80%) Gnaphosidae was the most species-rich family (32 spp., 24.81%), followed by Salticidae (17 spp., 13.18%), Lycosidae (15 spp., 11.63%) and Zodariidae (8 spp., 6.20%).Total species richness was quite similar between habitats, ranging between 34 and 48 species.However, Chao1-estimated species richness varied considerably, between 58 and 106 species per habitat, with the total ground-dwelling assemblage in the sampled habitats estimated at 167 species (Table 3).Although coverage values were above 0.85 for all the habitats, sample completion was much more variable (0.41-0.76), with 0.77 for the total assemblage (Table 3).This pattern was confirmed by the sample rarefaction curve for the whole spider assemblage (Figure 2a), which approached an asymptote but did not level out, indicating that the majority of the ground-dwelling species in the sampled habitats had been collected.However, none of the habitats' rarefaction curves approached an asymptote (Figures 2b-i), indicating that a considerable portion of the fauna of each was still unsampled.Because open grasslands were not sampled in this study, it could be expected that the total species richness at APGR may further exceed this projected value of 167 species.
Nearly half of all the spiders (n = 1480, 49.63%) were collected from the two cultivated D. eriantha pasture sites; indeed, spider numbers in this habitat were exactly five times higher than the habitat with the second highest activity densities, Searsia burchellii closed evergreen shrubland (n = 296).However, linear multivariate regression showed no significant differences between habitats in log-transformed spider activity densities (F 2,61 = 1.125, p = 0.3313) or species richness (F 2,61 = 0.3105, p = 0.7343).
Linear multivariate regression showed that there was significant seasonality across all sites in log-transformed spider activity densities (F 2,61 = 8.388, p = 0.0006), with summer and spring activity densities clearly much higher than those of the colder seasons (Figure 3a).Species richness showed a similar pattern (F 2,61 = 13.72,p < 0.0001), although spring and summer species richness was very similar (Figure 3b), although markedly higher than autumn and winter species richness.
Spider assemblages showed some interesting patterns.Not surprisingly, the assemblages of the cultivated D. eriantha pasture sites were the most distinct and showed the greatest similarity to one another compared to all the other habitat site-pairs.Among the remaining seven habitats, only the three habitats associated with the hill had their paired sites grouping together (Figures 4a, b).There was considerable overlap in the assemblages of the lower lying shrubland types and the dam, with most site-pairs not grouping close together (Figure 4b).These results were supported by the PERMANOVA analysis, which showed highly significant differences between the spider assemblages in the sampled habitats (pseudo-F = 2.614, p < 0.0001).Particularly, the pairwise post hoc comparisons showed significant differences in assemblages between the cultivated pastures (p < 0.0007), plateau (p < 0.0128) and southern slope (p < 0.0105) and all the other habitats (Table 4).For the other habitats, several other paired habitats were also significantly different from each other, although the assemblages of the freshwater dam and Searsia burchellii closed evergreen shrubland had the fewest significant paired values, confirming their assemblage overlap with other habitats (Table 4). Of

Discussion
Previous studies on ground-dwelling spiders in the Grassland Biome have focused on open grasslands, while the shrubland faunas have received little attention (Butler & Haddad 2011).
In this study, the first focused on shrubland habitats and cultivated pastures in the grasslands of central South Africa, the hill-associated habitats and cultivated D. eriantha pastures had the most distinct assemblages, while those of the lower lying shrublands and freshwater dam showed considerable overlap and often lacked clear distinction.The sample completion values varied between 0.41 and 0.76 for each habitat, which suggests that further sampling is necessary for a better representation of the species richness of those habitats with sample completion < 0.5.However, had the data from the other months sampled, but not included in this analysis, been incorporated into this study, then it is quite likely that this threshold value would have been exceeded in all the habitats.
Our results indicate that the aspect of hillside habitats has a strong influence in shaping assemblages.This could be because of the northern slopes getting more exposure to direct sunlight compared to the southern slopes.Southern slopes have denser vegetation cover in the woody layer, with especially Olea europaea africana providing a lot more shade.Assemblages of the plateau are also very unique, which can possibly be attributed to the differences in the herbaceous layer and slope of the habitat.We would propose conducting a large-scale survey of various hills in central South Africa to determine whether each aspect of these hills (including east and west that were not sampled here) contains a distinct assemblage and whether this pattern varies geographically.
The assemblages of cultivated D. eriantha pastures were very unique, largely because of the absence of woody plants.
Although studies in open grasslands in South Africa have shown very contrasting patterns of family dominance   (Haddad et al. 2015;Jansen et al. 2013;Lotz et al. 1991;Van den Berg & Dippenaar-Schoeman 1991), the dominance of Ammoxenidae in the cultivated D. eriantha pastures studied here is quite extreme (66.15%).This could be attributed to the high abundance of Hodotermes mossambicus (Hagen, 1853) termites in the pastures, which are the sole prey of A. amphalodes (Petráková et al. 2015).Similarly, the high activity of Ranops sp. in this habitat (13.24%) can be attributed to high activity densities of their mimetic model and possible prey ant, Anoplolepis custodiens (F.Smith, 1858) (Haddad 2012).Interestingly, these were the only two spider species with indicator values above 70.0%,suggesting that the other species sampled had more general habitat preferences, were more strongly seasonal in occurrence or were too scarce to serve as meaningful indicators of particular habitats.
The ten species of mygalomorph trapdoor spiders collected in four months' sampling (~9600 pitfall trap-days) is quite remarkable.This is higher than the eight species collected during nearly 22 000 pitfall trap-days' sampling in open grasslands in the Erfenis Dam Nature Reserve (Haddad et al. 2015)  was not collected in this study but has been recorded from the other two sites (Fourie et al. 2011;Haddad et al. 2015).
In addition, an unidentified Harpactira sp.(Theraphosidae) has also been collected from burrows and/or at night at all three localities, but has yet to be sampled using pitfall trapping.
Two faunistic records are of particular interest: the treetrapdoor genus Moggridgea (Migidae) is recorded from the Free State province for the first time (Dippenaar-Schoeman 2002;Dippenaar-Schoeman et al. 2010;Griswold 1987), with five male spiders sampled in very contrasting habitats (Appendix 1).This suggests that this species may occur in shrublands throughout central South Africa.This study also yielded the first records of Opiliones from the central Free State (Assamiidae: Polycoryphus asper Loman, 1902); all previous records of harvestmen in the province are from the eastern or southern fringes (Lotz 2002).Previously, this species was only known from the Kogelbeen caves in Northern Cape, Port Elizabeth in the Eastern Cape and on the Namibia-Angola border, and is thus widespread although scarce (Lotz 2009).Although only two specimens were collected in pitfall traps from the southern slope sites in this study, more than 20 additional specimens were collected from beneath large rocks on various hills at Amanzi, all on the southern slopes.Further studies are needed elsewhere in central South Africa to clarify the distribution and microhabitat preferences of these two arachnids.

Conclusion
This

FIGURE 1 :
FIGURE 1: (a) Location of the Amanzi Private Game Reserve, Brandfort district, South Africa (star), north-east of Bloemfontein (circle) and (b) vegetation map of the western section of the Amanzi Private Game Reserve.

FIGURE 3 :FIGURE 4 :
FIGURE 3: Seasonal variation in mean ± SD of (a) activity densities and (b) species richness of ground-dwelling spiders across 16 sites (eight habitats) sampled in the Amanzi Private Game Reserve.

TABLE 2 :
Family activity density and species richness of ground-dwelling spiders sampled by pitfall trapping from eight habitats in the Amanzi Private Game Reserve.

TABLE 3 :
Summary of actual and estimated species richness of ground-dwelling spiders sampled by pitfall trapping from eight habitats in the Amanzi Private Game Reserve.Searsia burchellii closed evergreen shrubland; C, cultivated Digitaria eriantha pastures; D, freshwater dam; N, Olea europaea closed evergreen shrubland on northern slopes of hill; P, O. europaea closed evergreen shrubland on plateau of hill; S, O. europaea closed evergreen shrubland on southern slope of hill; T, Tarchonanthus camphoratus closed evergreen shrubland; V, Vachellia karroo open deciduous shrubland.

TABLE 4 :
Post hoc p-values calculated from permutational multivariate analysis of variance analysis of spider assemblages associated with eight habitats in the Amanzi Private Game Reserve.
Note: Values in bold indicate significant differences between habitats.Habitat abbreviations: B, Searsia burchellii closed evergreen shrubland; C, cultivated Digitaria eriantha pastures; D, freshwater dam; N, Olea europaea closed evergreen shrubland on northern slopes of hill; P, O. europaea closed evergreen shrubland on plateau of hill; S, O. europaea closed evergreen shrubland on southern slope of hill; T, Tarchonanthus camphoratus closed evergreen shrubland; V, Vachellia karroo open deciduous shrubland.
and the five species collected in open grassland during 36 500 trap-days' sampling in the Free State National Botanical Gardens (J.A. Neethling & C.R. Haddad [University of the Free State] unpubl., August 2010 to May 2011), both within a radius of 70 km from Amanzi.This is probably because of the more structurally and topographically variable habitats sampled in the current study compared to open grasslands sampled in the latter two studies.The inclusion of open grassland habitats in the current study would likely have increased the trapdoor spider diversity at this site.For example, Calommata meridionalisFourie, Haddad and Jocqué, 2011 (Atypidae)

TABLE 1 -
A1: Ground-dwelling spiders collected by pitfall trapping from eight habitats in the Amanzi Private Game Reserve.