Effects of urbanisation and landscape heterogeneity mediated by feeding guild and body size in a community of coprophilous beetles

Although the impacts of urbanisation on biodiversity are well studied, the precise response of some invertebrate groups remains poorly known. Dung-associated beetles are little studied in an urban context, especially in temperate regions. We considered how landscape heterogeneity, assessed at three spatial scales (250, 500 and 1000 m radius), mediates the community composition of coprophilous beetles on a broad urban gradient. Beetles were sampled using simple dung-baited traps, placed at 48 sites stratified across three distance bands around a large urban centre in England. The most urban sites hosted the lowest abundance of saprophagous beetles, with a lower mean body length relative to the least urban sites. Predicted overall species richness and the richness of saprophagous species were also lowest at the most urban sites. Ordination analyses followed by variation partitioning revealed that landscape heterogeneity across the urban gradient explained a small but significant proportion of community composition. Heterogeneity data for a 500-m radius around each site provided the best fit with beetle community data. Larger saprophagous species were associated with lower amounts of manmade surface and improved grassland. Some individual species, particularly predators, appeared to be positively associated with urban or urban fringe sites. This study is probably the first to examine the response of the whole coprophilous beetle community to urbanisation. Our results suggest that the response of this community to urbanisation matches expectations based on other taxonomic groups, whilst emphasising the complex nature of this response, with some smaller-bodied species potentially benefitting from urbanisation.


Introduction
Urbanisation is among the foremost threats to biodiversity, with an increasing proportion of the global human population living in cities (United Nations 2014). Yet it is this high population density that makes urban green spaces a key point of interaction between humans and wildlife, driving the need for researching biodiversity in urban areas (Niemelä 1999). Key to this is understanding how functional traits drive community assembly. Urban species communities may be filtered by traits, at a regional scale (Croci et al. 2008a;Vallet et al. 2010), with a consequential reduction in phylogenetic diversity (Knapp et al. 2008). At a local level (i.e. along a single urban gradient) the community composition may be more randomly structured (Magura et al. 2018) and not necessarily filtered by the same traits as at a regional scalepotentially threatening ecosystem function (Croci et al. 2008a;Magura et al. 2018).
Historically, much urban biodiversity research concerned birds, which undergo significant community homogenisation and restructuring along gradients of urbanisation (Baker et al. 2010;Gagné and Fahrig 2011). An increasing number of studies address invertebrates in urban areas (Jones and Leather 2012) and show that the response to urbanisation is far from uniform, varying both within and between taxonomic groups (McIntyre et al. 2001;Gibb and Hochuli 2002;Angold et al. 2006;Egerer et al. 2017) and mediated by factors such as body size (Magura et al. 2006), feeding guild (Hochuli et al. 2004;Magura et al. 2013), degree of specialism (Gaublomme et al. 2008) and mobility (Angold et al. 2006;Snep et al. 2006;Delgado de la Flor et al. 2017).
These non-uniformly negative impacts mean that urban greenspaces can host a diverse range of species (Angold et al. 2006), and ecosystem function is not necessarily impaired by urbanisation beyond its fragmentation of habitats (Wolf and Gibbs 2004). Bee diversity can be particularly high in urban green spaces (Lowenstein et al. 2014;Baldock et al. 2015;Banaszak-Cibicka et al. 2018), suggesting that urban areas may contribute to maintaining pollination services across the wider landscape (Theodorou et al. 2017). Pest regulation may also be supported by the relatively high diversity of open-habitat spiders (Lövei et al. 2019). This highlights that urban areas are an important part of the contemporary landscape mosaic and their 'spillover' effect into agricultural landscape or neighbouring semi-natural habitats, whether positive or negative, needs to be addressed. Some argue that by hosting novel ecosystems (Hobbs et al. 2009;Kowarik 2011) urban areas actually boost overall biodiversity at a broad landscape scale, e.g., regional or national (Sattler et al. 2011). Species communities in a habitat patch may even be modified by the presence of urban areas in the surrounding landscape (Neumann et al. 2016a). To assess the full effect of landscape heterogeneity on biodiversity, including how urban areas interact with other land cover types, ideally a measure of gamma diversity would be obtained, i.e., the species diversity of the whole landscape. However, studies that address this are not particularly common (Duflot et al. 2014(Duflot et al. , 2017. Though not a complete substitute, sampling insects on ecotones using a common attractant, such as baited traps or flowers (e.g. Foster et al. 2019), may provide a snapshot of landscape biodiversity that can be assessed at a large number of sites relatively easily. Such a method can be carried out in urban sites as well as rural, where a habitat-focussed experimental design is not always feasible. It may also facilitate focussing on a particular functional group, which may be a more meaningful measure for assessing the value of urban greenspaces than the more random assemblages of species collected by passive trapping methods (Gagic et al. 2015;Pinho et al. 2016).
Scarabaeoidea are a well-studied group, especially the true dung beetles in the families Geotrupidae and Scarabaeidae subfamilies Aphodiinae and Scarabainae (henceforth 'dung beetles'). Dung beetles are important ecosystem-service providers in pastoral agricultural systems (Nichols et al. 2008;Manning et al. 2016) and often used as indicators in tropical ecosystems (Nichols et al. 2007), where they are threatened by habitat fragmentation including urbanisation (Korasaki et al. 2013). They are usually sampled using a variety of dungbaited trapping techniques or by direct searching of dung. Very few studies consider dung beetles in an urban context (Ramírez-Restrepo and Halffter 2016). Not all dung-feeding Scarabaeoidea are obligate dung feeders (Gittings and Giller 1997) and some will feed on dog dung, even preferentially so for species that prefer omnivore to herbivore dung (Cave 2004;Carpaneto et al. 2005). As such, a stable community of dung beetles might persist in urban green spaces that provide either a ready supply of dung or other sources of rotting organic material for species with broadly saprophagous larvae. In addition, other beetles besides Scarabs are coprophilous, attracted to dung either as a direct food source or to predate on the eggs and larvae of other insects (Hanski and Cambefort 1991). Many of these species belong to the Staphylinidae, which despite potentially being a good indicator group (Bohac 1999;Vásquez-Vélez et al. 2010) are less often the focus of urban or landscape studies (though see Delgado de la Flor et al. 2017;Magura et al. 2013;Vergnes et al. 2014).
In this study, we use simple baited traps to sample the coprophilous beetle community in and around a large urban area in southern England. By stratifying sampling sites into three urban distance bands, we investigated the effect of proximity to urban land cover. We expected the most rural sites to contain more resources for coprophilous beetles (due to higher density of both wild mammals (Baker and Harris 2007) and livestock) and therefore test the hypotheses that i) coprophilous beetle abundance and ii) species richness are highest at rural sites and lowest at urban sites. As urbanisation filters out larger species and habitat specialists in some other groups of beetles (Magura et al. 2010;Merckx et al. 2018), we also test the hypotheses that iii) beetles collected at urban sites will have a lower average body size than at rural sites and iv) community composition differs between urban distance bands. We examine the response of predatory and saprophagous species separately, anticipating that predators may be more able to exploit alternative habitat resources to dung. Finally, we explore the effect of landscape heterogeneity on the beetle community using ordination analyses, identifying the key landscape elements and the spatial scales at which they most strongly influence community composition in coprophilous beetles.

Study design
Forty-eight study sites were selected on woodland edges or the perimeters of partly wooded urban green spaces (edges >50 m in length). Whilst sampling may therefore be biased towards edge-preferring (or edge-tolerating) species, we considered this appropriate as the highly heterogeneous urban sites provided very little interior habitat. Sampling along edges was also designed to capture species preferring both open and shaded habitats, rather than confine the species pool under consideration to a single habitat. Overall abundance and diversity can be higher at habitat edges (Magura et al. 2017) though some dung beetle species may avoid ecotones (Martínez-falcón et al. 2018).
Study sites were selected in three groups based on their distance from Reading, a town at the centre of a large conurbation (area 83.7km 2 , population c. 340,000 (Brinkhoff 2018)) in southern England (lat. 51.45,. A circular area with a 15 km radius originating at the town centre was divided into eight sections following compass intervals, e.g., the area between bearings north and northwest. In each of the eight sections six sites were selected. Two of the six were 'Urban' sites, defined as being less than 200 m from the main Reading conurbation, as mapped in Land Cover 2007 (Morton et al. 2011). Two were 'Fringe' sites, between 200 m and 1000 m away from either the Reading conurbation or other urban areas >10 ha in extent. Finally, two 'Rural' sites were selected, all of which were > 1000 m away from any urban area > 10 ha. Urban, Fringe and Rural are henceforth referred to as urban distance bands. All selected sites were at least 1 km from their nearest neighbour to maximise independence between beetle populations. During fieldwork, site locations occasionally had to be changed to find a suitable place to set traps. Sites were subsequently recategorised after fieldwork was complete, again according to the distance bands above but using the precise field survey locations and updated urban land cover data from Land Cover 2015 (Rowland et al. 2017). This resulted in a final sample of 15 'Urban' sites, 16 'Fringe' sites and 17 'Rural' sites ( Fig. 1).

Beetle sampling
Due to the high likelihood of trap disturbance in the urban sites, we assessed the dung-attracted beetle community using simple baited traps deployed for a short time period. Traps consisted of a 140 mm wide, 20 mm deep Petri dish filled to the brim with horse dung, which is attractive to a wide variety of species (Mroczyński and Komosiński 2014). Dung was collected early in the morning from an enclosed stable in order to reduce the potential for colonisation by beetles ahead of its use for sampling (Krell 2007). Between collection and use for trapping, the dung was kept inside in a covered bucket to keep it moist and exclude insects. Blends of dung collected on three different days were used to fill the traps in order to reduce variability in the age and consistency of dung provided between sampling sites. Any dung older than three days was discarded and not used for sampling since dry dung is less attractive to dung beetles (Aschenborn et al. 1989;Hanski and Cambefort 1991) and possibly other coprophilous beetles.
Sampling sites were visited twice each between 19th May 2015 and 25th July 2015 in dry, calm conditions. Two randomly selected sites in each urban distance band were visited on each day of fieldwork. On the day, sites were visited in a random order to further reduce collinearity between weather or time of day and urban distance. Four traps were set out in a line at each sampling location, following a woodland edge or other semi-natural habitat boundary. Traps were placed on the ground at least 10 m apart. Presence/absence of livestock at the sampling location was also recorded during sampling as dung availability is known to modify the efficacy of dungbaited traps (Finn et al. 1998). After one hour, traps were sealed using the petri dish lid and returned to the laboratory.
Trap contents were emptied into a bucket of water to float out any adult beetles present, which were removed using a small sieve. Beetles were stored in 70% ethanol prior to identification. Beetles were identified to species following the relevant identification literature (Lott 2009;Lott and Anderson 2011;Duff 2012;Dung Beetle UK Mapping Project 2018). Members of the family Ptilidae and subfamily Staphylinidae: Aleocharinae were counted a maximum of once each towards the species total for each site but not identified to species for use in community analysis. We used double-blind identification for quality control. Where necessary, specimens were dissected to confirm identification and compared to reference material in the Oxford University Museum of Natural History (Specimens for this study are held in the Centre for Wildlife Assessment and Conservation at the University of Reading and available to view on request). Species were split into larval feeding guilds based on expert information held in the PANTHEON database (Webb et al. 2018). The majority fell into one of two groups, either saprophagous (feeding on rotting organic material including dung) or predatory. They were also assigned an average body length from information provided in the identification resources for each family.

Spatial scale, landscape variables
Dung beetles in pasture habitat generally remain within the field in which they emerged (Roslin 2000), though dispersal of individuals is leptokurtic, with some 1% of marked beetles observed by Roslin (2000) found 1 km away. If the appropriate spatial scale at which to examine the effects of landscape heterogeneity can be drawn from recorded dispersal flight distances, evidence from other beetle families is mixed. Irmler et al. (2010) reported very few saproxylic beetle species in flight intercept traps >80 m from woodland, whereas Doležal et al. (2016) observed a maximum movement in a bark beetle of over 1 km. The harlequin ladybird Harmonia axyridis may be capable of regular dispersal flights >18 km (Jeffries et al. 2013), so clearly some beetles are unlikely to be impacted by habitat fragmentation. Most flighted rove beetles are probably highly mobile, potentially dispersing long distances assisted by weather patterns. To account for this uncertainty, landscape variables were extracted for three different scales. These were 250-m, 500-m and 1-km buffers around the centre of the sampling sites. Landscape variables are split into two conceptual categories: landscape composition, which is the amount of different land cover types, and landscape configuration, which describes the spatial arrangement of those types, i.e., heterogeneity and fragmentation.

Landscape composition
The topography of land cover patches for landscape variables was taken from OS Mastermap (Ordnance Survey 2015) and analysed in ArcGIS 10.4 (ESRI 2016). All patches intersecting 2000-m buffers around the sampling points were classified into one of eight landscape composition variables, in total comprising a complete landscape mosaic analysis. All buildings, structures, paved roads and unsurfaced tracks were classified as MANMADE. Patches listed in Mastermap as 'Make = Multiple' (a mixture of manmade and natural surfaces) were classified as GARDEN. Semi-natural habitat patches were classified as either WOODED (nonconiferous trees, coniferous trees, coppice), SCRUB (scrub, scattered coniferous trees, scattered non coniferous trees, orchard) or OPEN (freshwater marsh, heath, heather grassland, neutral grassland, acid grassland and calcareous grassland). Linear patches of habitat alongside railway lines were classified as OPEN. Linear patches alongside roads were also classified as OPEN unless the underlying land cover type in LCM2015 (Rowland et al. 2017) was urban, in which case they were classified as improved or amenity grassland (IMPGRASS). Freshwater types including rivers, lakes, streams and ditches were classified as WATER. Remaining unclassified polygons with no attributes listed in OS Mastermap were assigned to a category with reference to Natural England priority habitats data (Natural England 2018) and Land Cover Map 2015 (Rowland et al. 2017), categorising agricultural land as either improved grassland (IMPGRASS) or arable/horticulture (ARABLE).

Landscape configuration
A variable estimating the availability of livestock dung (DUNGAV) in the surrounding landscape (Webb et al. 2010) was created by multiplying the area of each grassland patch (improved and semi-natural) by a 'livestock index'. The index was derived from livestock headage data (DEFRA 2010) at a 5 km 2 grid cell resolution. Vales for each cell were rescaled to 0-1 where 0 = no livestock in the corresponding 5-km 2 grid Fig. 1 (a) The location of 48 sampling sites in the vicinity of Reading, UK, separated into three types according to their distance from urban areas >10 ha. Urban areas in grey shading, blue shading denotes rivers and other bodies of freshwater. Landscape composition for a 250 m and 500 m radius around the sites (black circles) is shown for (b) an urban site, (c) a 'fringe' (periurban) site and (d) a rural site. © Crown copyright and database rights 2019 Ordnance Survey (100025252) cell and 1 = the maximum value for any cell intersecting with the study area.
To calculate mean patch size for the semi-natural habitat categories WOOD, SCRUB and OPEN, any patches closer than 20 m were functionally considered to represent a single patch (Neumann et al. 2016b) and aggregated into one. Most beetles attracted to dung-baited traps are likely to arrive by flying, so gaps of this size are unlikely to present a significant barrier to dispersal. The mean size of aggregated patches that intersected the landscape buffers was then used to create the variables WOPATCH, SCPATCH and OPPATCH. Road verges were not included in the patch size calculation for OPPATCH, so this variable refers only to more extensive semi-natural grassland. As a further measure of habitat fragmentation, the mean distance to each of the three semi-natural habitats within landscape buffers was calculated, forming the variables WOODDIS, SCRUBDIS and OPENDIS. These were based on the original non-aggregated habitat patches, again excluding road verges for OPENDIS.
Finally, mean garden size (GARSIZE) was included as an alternative measure of urbanisation. Garden size tends to be smaller in core urban areas and larger in suburbs (Smith et al. 2009). Larger gardens may be more likely to contain a diverse mix of habitats, with a range of sources of decaying vegetation such as compost heaps that may provide habitat for saprophagous beetles that are also attracted to dung. A summary of all landscape composition and configuration variables and their mean values in each urban distance band is given in Table 1.

Statistical analysis
All analysis was carried out in R 3.4 (R Core Team 2017). Data from the four traps at each sampling site were pooled. Data from the two visits were considered separately when testing for differences in abundance and species number. For species richness estimates and all community analyses, data from the two visits were pooled into a single value for each site.
Landscape variables, abundance, species richness and mean body length Differences in the values of i) landscape variables, ii) beetle abundance, iii) species number and iv) mean body length between urban distance bands were assessed using Kruskal-Wallis tests (adjusted for ties), followed by Dunn's test for multiple comparisons with Bonferroni corrected p values. These were carried out using the R package 'dunn.test'.
For abundance and species richness, 12 initial Kruskal-Wallis tests were performed, examining the abundance and species richness in visit 1 and 2 of i) all beetles ii) those with saprophagous larvae and iii) those with predatory larvae. As numbers of dung beetles appeared to strongly influence the results for saprophagous species, the abundance of these species alone during each visit were examined in two further tests. Estimated species richness for each distance band was calculated using the function 'Specpool' in the R package 'Vegan' (Oksanen et al. 2016), examining predicted species richness ± standard error using the Jacknife, and Bootstrap methods (Smith and van Belle 1984;Palmer 1990). Estimates were made for all beetles and for those with saprophagous larvae only.

Community analysis: Distance categories
All community analyses were carried out using the R package Vegan (Oksanen et al. 2016). Permutational analysis of variance was performed using the function Adonis (with 999 permutations) in order to assess whether beetle communities differed between urban distance bands. Beetle community data were standardised using the alternative Gower measure (Anderson et al. 2006) and log base 2 transformed such that a doubling in abundance was weighted the same as a change in composition (we were equally interested in both potential modifications of the species community). Analysis of multivariate homogeneity of group dispersions using the function Betadisper was used to check whether any observed difference in community structure between distance bands might in fact be attributable to within group variation. Tests were carried out for all beetles and for each of the main feeding guilds (saprophagous or predatory) only.

Community analysis: Landscape heterogeneity
Only species occurring in four or more sites were included in the community data used for ordination analyses. The community data matrix was log transformed (log 2 (x) + 1 for x > 0) (Anderson et al. 2006). Key landscape composition and configuration variables were identified via canonical correspondence analysis (CCA), followed by a forward selection procedure. Landscape composition variables were log transformed to reduce multicollinearity between predictors and ensure the influence of small but potentially important components of the landscape was not overlooked (Neumann et al. 2016a).
Six initial models were run in total, one for all landscape composition and one for all landscape configuration variables at each of the spatial scales (250 m, 500 m and 1000 m buffers). Collinearity between predictors was tested using Variance Inflation Factors (VIF) (Neter et al. 1996) and correlation coefficients. Terms with high VIF were removed sequentially until all terms in the global model had VIF <5.0 and correlation coefficients were all <0.6 (Radford and Bennett 2007;Neumann et al. 2016b).
Forward selection to obtain a reduced model for composition or configuration at the three scales was carried out only after a significant global test of all variables (p < 0.05). Terms were added sequentially to the model based on their predicted contribution to adjusted-R 2 (Peres- Neto et al. 2006), and retained only if significant at p < 0.05 (Blanchet et al. 2008). Model building stopped before the adjusted-R 2 of the reduced model exceeded that of the global model to avoid over-fitting of explanatory variables (Blanchet et al. 2008). Probability values for the global model and individual terms retained in the reduced model were estimated using Monte Carlo tests with 9999 permutation via the function Anova.cca.
A final model for each of the three spatial scales was obtained by adding all significant composition and configuration variables as predictors and removing any that were collinear by referring again to VIF and correlation coefficients. If any terms in the resulting model were still non-significant, they were considered likely to be explaining the same environmental gradient and a further forward selection procedure (without the adjusted R 2 stopping rule) was carried out to obtain a more parsimonious model. Finally, this procedure was repeated with the best fitting variables from each scale entered into a single model to obtain a final set of landscape variables that best described the beetle community composition.

Spatial and temporal autocorrelation, variation partitioning
The effects of temporal autocorrelation and livestock presence at a site were also assessed using CCA following the model selection steps above. Dates of the first and second round of sampling were entered as two continuous variables; this controls for peaks in species activity as CCA assumes a unimodal, not linear relationship. Livestock presence was a factor with two levels indicating simple presence/absence.
We tested for spatial autocorrelation by calculating Principle Coordinates of Neighbourhood Matrices based on the XY coordinates of the sampling sites. All PCNM were used as predictors in a CCA followed by a further forward selection procedure to identify significant gradients of spatial autocorrelation. The interaction between space (represented by significant PCNM), landscape (key landscape variables) and date/livestock presence was then explored using variation partitioning. This finds the unique fraction of variation explained by each set of constraints, based on their contribution to adjusted-R 2 (estimated for CCA by permutation following Peres-Neto et al. 2006). Probability values for the unique effects of space, landscape or time/dung were estimated by using partial canonical correspondence analysis to partition out the other two sets of variables, followed by Monte Carlo tests with 9999 permutations.

Beetle sampling
A total of 2636 individual beetles were collected from the 48 sites, with representatives of 34 species and a further two groups not identified to species (Staphylinidae: Aleocharinae, Ptilidae). The 19 species recorded from four or more sites are listed in Table 2. Of these, seven belonged to the saprophagous larval guild (of which some also feed on dung as adults); the rest all had predatory larvae and Heleophorus brevipalpus are also predators as adults. Five of the top 10 most abundant species, including the two most frequently recorded, belonged to the subfamily Oxytelinae (Staphylinidae), which are often among the most abundant rove beetles in dung (Caballero and León-Cortés 2012;Mroczyński and Komosiński 2014;Yamamoto et al. 2014). Eight species of dung beetle were recorded, with four among the most widespread species. The most common, Volinus sticticus, is widespread in woodland and scrub habitats (Dung Beetle UK Mapping Project 2018).

Abundance
There was a significant difference between distance bands in visit 1 for all species with saprophagous larvae (K-W χ 2 = 11.660, p < 0.01) but not for all beetles together (K-W χ 2 = 5.074, p = 0.08). For saprophagous species, numbers at Urban sites were significantly lower than at Rural sites (z = 3.414, p = 0.001) but the difference between Fringe and Urban sites was not significant (z = 1.718, p = 0.130;) (Fig. 2a). There was no difference in the abundance of predators between distance bands (K-W χ 2 = 3.21, p = 0.20) (Fig. 2c).

Species richness
Neither the overall species richness (K-W χ 2 = 3.638, p = 0.16) nor the richness of saprophagous species (Fig. 2e; K-W χ 2 = 4.939, p = 0.08) was significantly different between urban bands during the first visit, though the number of dung beetle species was (K-W χ 2 = 8.513, p = 0.01). More dung beetle species were recorded at Rural sites than Urban sites (z = 2.906, p = 0.006). There was no difference in the species richness of predators ( Fig. 2g; K-W χ 2 = 1.652, p = 0.44). The same pattern occurred during the second visit, with no difference between bands for total species number (K-W χ 2 = 1.475, p = 0.48), saprophagous species ( Fig. 2f; K-W χ 2 = 1.605, p = 0.45) or predators (K-W χ 2 = 0.862, p = 0.65). Between the first and second visit patterns of species richness in predators reversed, with the highest values in Rural sites in visit 1 (Fig. 2g) and Urban sites in visit 2 (Fig. 2h). Both estimates of species richness for urban bands predicted that Rural sites are the most species rich overall, higher than Fringe, with Urban the most species poor (Table 3), though the standard error ranges overlapped in each case. The same predicted pattern of species richness occurred for saprophagous species. The bootstrap estimate predicted that Rural sites host 50% more saprophagous species than Urban sites, with no overlap in standard error range between Rural and Urban sites. Conversely, Rural sites are predicted to have the lowest richness of predatory species with the highest richness in Fringe sites (Table 3).

Mean body length
There was a significant difference in mean body length between distance bands in the first visit (K-W χ 2 = 6.033, p = 0.05). Beetles at Rural sites were on average very slightly larger than those at Urban sites (Fig. 3a, z = 2.393, p = 0.025). In the second visit, there was no difference in average body length between distance bands (K-W χ 2 = 1.826, p = 0.40), driven by greater variation among Fringe and Rural sites (Fig. 3b).

Community dissimilarity
In the permutational analysis of variance, communities were significantly dissimilar between urban band (F = 1.931, p = 0.032), which explained 8.7% of the dissimilarity between sites. There was no significant difference in within band homogeneity (F = 0.539, p = 0.587). When Fringe sites were excluded, the difference between Urban and Rural sites alone explained slightly more variation (9.1%; F = 2.889, p = 0.009). When only saprophagous species were examined, urban distance band explained 15.0% of variation between sites (F = 3.174, p = 0.013), whereas for predators urban distance band explained only 4.5% of variation between sites and was not significant (F = 1.027, p = 0.397).

Landscape variables
Cover of manmade surfaces, garden and arable land were significantly different between urban distance bands at all measured spatial scales. Woodland cover, dung availability, woodland patch size, distance between woodland patches and average garden size were all different between distance bands when measured at 500 m or 1 km. Improved grassland was only significantly different at 1 km. The highest woodland cover and patch size was around Rural sites. Manmade surface and garden cover was much higher in Urban than in Rural or Fringe sites, which were not significantly different from each other. The only variable to completely separate out the three distance bands was average garden size measured at 1 km, for which there was a significant difference between each pair of categories (Table 1).

Landscape heterogeneity: Canonical correspondence analysis
There was a significant global effect of landscape composition and landscape configuration variables for each spatial scale (Table 4). Landscape variables measured at 500 m performed best in explaining variation in beetle community composition. At 250 m only semi-natural habitats (the size of woodland and open habitat patches) and water were identified as key variables. At 500 m the only semi-natural habitat variable retained was the size of scrub patches, otherwise manmade surfaces, arable and improved/amenity grassland best described the landscape in terms of beetle community response. At the 1000 m radius manmade surfaces were highly collinear with other variables. After MANMADE was removed from the model the extent and average size of gardens described the urban gradient whilst the extent of scrub or scattered trees was also important at this scale.
All of the variables from 500 m were retained in the final model, strongly suggesting that this is the best scale at which to assess the response of this beetle community to landscape heterogeneity. Area of water in the surrounding 250 m was the only other variable retained. On the CCA bi-plot (Fig. 4), the first axis (CCA1, F = 5.269, p < 0.001) explained 10.2% of the variation and described a gradient away from more urban landscapes (as measured by manmade surface coverage) and those with greater extent of improved or amenity grassland. The second axis (CCA2, F = 3.271, p = 0.019) explained 6.3% of the variation and described a gradient from landscapes with large patches of scrub or scattered trees and larger amounts of water nearby, to drier landscapes with greater arable cover.
In the variation partitioning, the contribution to adjusted-R 2 of space (significant PCNM), landscape (final model for landscape heterogeneity) and date / stock presence was 7.2%, 5.8% and 4.5%, respectively (Fig. 5). 5.3% of the variation was shared between landscape and space, and 2.3% shared between all three fractions. Overall variation explained was estimated at 25.8%. All individual fractions were significant after the other two were partitioned out (Space: F = 1.5985, p = 0.005; Landscape: F = 1.553, p = 0.005; Date/Stock: F = 1.674, p = 0.012).

Coprophilous beetle community
For the whole coprophilous beetle community, urban distance band explained a relatively small proportion of the dissimilarity between sites, even when using a dissimilarity measure that deliberately gave a high weight to changes in abundance. This suggests that in urban areas this community represents a subset of the wider species pool found in surrounding agricultural landscapes and semi-natural habitat patches beyond, rather than a novel assemblage of species (cf. Hobbs et al. 2009;Kowarik 2011). It also appears that the coprophilous beetle communities in this study were not consistently homogenized by urbanisation as found e.g. by Magura et al. (2010) for ground beetles. This could be connected to dispersal ability: while some brachypterous ground beetles are reluctant to cross roads (Koivula and Vermeulen 2005) and may be absent from urban habitats (Kegel 1990), all of the coprophilous species sampled in our study would have arrived at the traps on the wing.
Several urban sites were similar in species composition to some of the rural or fringe sites, whilst the overall abundance and richness was not significantly lower in urban sites. In other words, the 'best' semi-wooded urban green spaces can host a rich diversity of coprophilous beetle species, sometimes comparable to similar sites further from urban settlements. Similar patterns have been observed in other insect communities (Angold et al. 2006;Croci et al. 2008b).

Dung beetles
Analysing data from the first and second visits separately revealed that patterns of activity in the beetle community were seasonally dependent (Finn et al. 1998). During the first visit, richness and abundance of all saprophagous beetles increased sequentially from urban sites through to rural. In the second visit there was little difference between the urban bands. This can be explained by the active season for dung beetles. These were significantly more abundant at rural sites but only during the first visit in late May and early June, when spring active dung beetles are still abundant (Finn et al. 1998). All of the larger saprophagous species collected were dung beetles belonging to the tribe Aphodiini or the genus Onthophagus. Dung beetles were recorded at only four urban sites during the study; these were all at sites on the outside edge of the Reading conurbation not completely surrounding by built-up areas. There are two possible reasons for this. Either urban areas present dispersal barriers for larger, slow-flying beetles (Magura et al. 2006;Croci et al. 2008b) or their populations are limited by lower resource availability in urban habitat patches (Gibb and Hochuli 2002). Previous casual observations of Aphodiini and Onthophagus species by the authors in a large urban green space in the study area (not sampled in the current study) possibly suggest that resources rather than dispersal limitations are the key limiting factor for urban dung beetles.

Staphylinidae on dung
Very few studies have examined the community of Staphylinidae on dung. Yamamoto et al. (2014) determined that the Staphylinid fauna on deer dung in Japan was dominated by a saprophagous Anotylus species. Three saprophagous Staphylinidae were found in the current study, of which two, Platystethus arenarius (like Anotylus in the subfamily Oxytelinae) and Megarthrus prosseni, were recorded at more than 10 sites, though none were nearly as abundant as the predatory Staphylinids that dominated the highest beetle counts in this study. Platystethus arenarius was present in all urban distance bands but was particularly common in urban sites during the second visit, associated with large patches of scrub in moderately wet landscapes. Megarthrus prosseni was also more frequent at urban sites, but in association with greater improved grassland cover. M. prosseni occurs in all kinds of decaying vegetation (Cuccodoro and Löbl 1997) and has been found frequently in horse dung in one UK county (Lane et al. 2002), which was used as the bait for traps in this study, though in a study in Poland it showed no preference for horse dung over cow (Mroczyński and Komosiński 2014). These species could be more common in the vicinity of urban areas because there tends to be more suitable habitat available, but following a similar experimental design Magura et al. (2013) found that the abundance and species richness of saprophagous Stapylinidae was higher in rural forests than urban ones. Alternatively, they may be experiencing competitive release in the absence of dung beetles or a release from certain predators (MacArthur and Levins 1964;Atkinson and Shorrocks 1981;Kruess and Tscharntke 1994).

Predatory species
The dominant species across all sites were two predatory Anotylus species. Considering the short trapping window of one hour, they were extremely abundant at some sites, with four counts of over 100 individuals recorded. Interestingly, all of the 100+ counts were at urban (3) or fringe (1) sites. In terms of species richness, samples taken during the first visit would confirm the idea that predators can be limited by urbanisation (Hochuli et al. 2004;, though the effect observed here was small, equating to one fewer active species at urban sites than rural ones. On the second visit there were no clear differences between urban, fringe and rural sites,  Table 2 suggesting that phenology of some coprophilous beetles may vary between rural and urban sites, as observed in butterflies (Dennis et al. 2017). Overall, there is no consistent evidence that urbanisation strongly impacts the abundance or richness of predatory beetles associated with dung and other patchy habitats. There is some evidence for a shift in the type of species that occur, however. Although dissimilarity analysis found no significant difference between distance bands, in the canonical correspondence analysis the smallest predators in each family (Anotylus spp., Cercyon pygmaeus) were mostly associated with higher manmade surface cover surfaces, while the larger species (Philonthus spp., larger Cercyon) tended to be associated with small amounts of manmade surfaces. The larger predators often occurred in association with dung beetles, perhaps because the presence of dung beetles indicates a site rich enough in dung or other patch habitats to support good numbers of all saprophagous insect larvae.

Landscape heterogeneity
The combination of higher numbers of dung beetles, larger predators and higher mean body size in rural areas (at least during the first visit, when the larger Scarabs were more active) suggests a shift along the urban gradient from larger specialist species to smaller generalists. This fits the trend frequently observed in ground beetles (Alaruikka et al. 2002;Gaublomme et al. 2008;Magura et al. 2006), even though, unlike ground beetles, this community of insects is unlikely to be limited by dispersal. The main exceptions to this trend were Bisnius fimetarius and Gabrius piliger, which favoured manmade surfaces, water and large patches of scrub, G. piliger very strongly so. Both are associated with any patch habitat that has a high concentration of insect larvae (Lott and Anderson 2011), but little other information is available regarding their ecology. The peri-urban areas of Reading have a good deal of scrub cover, often bordering large water bodies or rivers (Fig. 1); suburbs offer a heterogeneous environment that may be associated with high species richness (Blair and Launer 1997;Blair and Johnson 2008). This could explain the preference of some species for landscape types that fit this description and the higher predicted richness of predators at Fringe and Urban sites compared to Rural in our study area. Some of the trends observed in the beetle community may also be related to interactions with other coprophilous organisms, such as flies or fungi, which it is not possible to quantify with our data. When landscape heterogeneity was considered at a 250 m radius, large patches of woodland and semi-natural grassland were identified as drivers of beetle community composition. They were not included in the final model as other landscape variables measured at 500 m better described the same gradient of variation in the beetle data. However, the fact that woodland and grassland patch sizes were not significantly different between urban distance bands suggests that larger patches of semi-natural habitat are valuable regardless of the location of the site on the urban gradient (Wolf and Gibbs 2004;); Söderström et al. (2001) also found a positive association between tree cover in the surrounding landscape and some dung beetle species recorded in grasslands. The amount of water within 250 m explained some of the variation on the second axis of the correspondence analysis. This may be linked to the preference of some species for damper soils or be a proxy for habitat quality close to water.
The inclusion of improved grassland in the final model might be expected if this cover type was linked to the presence of grazing livestock. In this case it was positively correlated with manmade surfaces, suggesting that the influential grassland type here was more likely to be amenity grassland, e.g., in parks or golf courses. Dung beetles were in fact associated with landscapes with relatively low improved grassland cover and the dung availability index was not included in any model. In this landscape they appear therefore to be more associated with semi-natural habitats, suggesting that livestock farming in the area does not support a high abundance of dung beetles (Hutton and Giller 2003;Webb et al. 2010).
Most of the key landscape variables identified describe landscape composition. The 500 m radius that best explained community composition is comparable to scales used in other landscape studies (Barbaro et al. 2007;Sjödin et al. 2008); in most of these cases measures of landscape composition or habitat amount were again more important than landscape, depending on the mobility of the species group. Given the long dispersal distances reported for some dung beetles (Roslin 2000), the effects of landscape configuration may operate at a larger scale than any measured here or may simply not be important for this community of beetles, especially as none of the species recorded is regarded as particularly specialist (Roslin 2001;Roslin and Koivunen 2001).
Landscape composition in this case is likely to be a proxy for the availability of suitable microhabitats, since the species community is formed around a patchy resource. The final proportion of variation explained by key landscape variables was not particularly high, especially after the portion attributable to spatial autocorrelation was removed, suggesting that local factors may be more important (Söderström et al. 2001;Philpott et al. 2014;Otoshi et al. 2015). The most important gradient in the beetle community is correlated with an increase in manmade surfaces, which is difficult to disentangle with the wider effects of distance from urban area. It is therefore not possible to be certain whether urbanisation, measured as an increase in impervious surfaces, is directly responsible for the trends observed in the beetle community or whether other features of the landscape that are correlated with urbanisation are more influential, as for habitat fragmentation and forest carrion beetles (Wolf and Gibbs 2004). The canonical correspondence analysis in any case usefully illustrates the character of landscapes that host different assemblages of coprophilous beetles, even if all the relationships described are not directly causal.

Summary
This is probably the first study to consider the composition of the whole beetle community on dung on an urban gradient. Our hypothesis that coprophilous beetles are more abundant at rural sites (i) was true for saprophagous species but not for predators. Similarly,species richness (ii) differed only for dung beetles in the early season visit; estimates of total community species richness for saprophagous species suggest that richness is genuinely lower at urban sites while predators are richest at fringe sites. Body size (average body length) followed the expected pattern (iii) in the early season visit, driven by an absence of dung beetles from urban sites. Community composition differed significantly using a distance measure that gives a high weight to abundance, suggesting that our hypothesis (iv) is true depending on how community composition is defined.
While landscape features did not explain a large amount of variation in the coprophilous beetle community, there was some modification of the community by landscape heterogeneity at a 500 m radius that should be taken into account. Species displayed both positive and negative associations with manmade surfaces, again indicating a modification of the community along the urban gradient, while the abundance of others was correlated with other landscape components indicating an effect of heterogeneity beyond that associated with urbanisation.
The experimental design was deliberately robust to disturbance in urban areas, using traps deployed for short time periods that nonetheless detected a difference between the three urban distance bands. This could provide a useful model for other studies. A potential alternative would be to use private gardens as sampling sites, which has been successful for a number of other insect sampling studies, including in the same study area (Orros et al. 2015;). An interesting possibility in suitably undisturbed sites would be to directly assess ecosystem function by quantifying the decomposition rates of dung or other rotting vegetation, in parallel with observations of the species composition of the whole functional community of saprophagous insects and associated predators.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.